如何在src文件中创建标量物理filed并且调用

chengan.wang

大家好，请教一个问题：
我知道如何在求解器中的createFields.H建立field，

    volScalarField summass
    (
        IOobject
        (
            "summass",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        mesh,
        dimensionedScalar ("zero", dimMass, 0.0)
    );

我想进一步在src文件中创建标量物理filed并且调用，比如ParticleCollector.H，ParticleCollector.C等库文件？

谢谢大家。

李东岳

@chengan-wang 在如何在src文件中创建标量物理filed并且调用中说：

src文件中创建标量物理filed并且调用，比如ParticleCollector.H，ParticleCollector.C等库文件？

为什么要必须在src文件中，openfoam的src文件？

chengan.wang

@李东岳感谢李老师回复。
我先在求解器中的createFields.H建立过这个field:Euler-Lagrange的一些解析，sprayFoam

之后在ParticleCollector.C文件添加了

cellI = p.cell();
scalar& PPC = summass_->primitiveFieldRef() [cellI];
PPC += p.nParticle()*p.mass();

但我编译时候出现错误，找不到summass

所以我猜应该是把这个field建立到ParticleCollector.H，ParticleCollector.C等文件中吧？！瞎猜的。

李老师能指点一下吗？

李东岳

@chengan-wang 你在库文件里面要使用求解器层面的volScalarField是调用不起来的。你可以试一下用这个:

const volScalarField& summass = p.mesh().lookupObject<volScalarField>("summass");

scalar& PPC = summass.primitiveFieldRef() [cellI];

chengan.wang

@李东岳谢谢李老师帮忙！
我我的是2.4.0版本，这样写的：

        const volScalarField& summass = p.mesh().lookupObject<volScalarField>("summass");
        const label cellIp = p.cell(); 
        scalar& PPC = summass.primitiveFieldRef() [cellIp]; 
        PPC += p.nParticle()*p.mass();

运行出现如下错误

lnInclude/ParticleCollector.C: In member function ‘virtual void Foam::ParticleCollector<CloudType>::postMove(Foam::ParticleCollector<CloudType>::parcelType&, Foam::label, Foam::scalar, const point&, bool&)’:
lnInclude/ParticleCollector.C:748:77: error: expected primary-expression before ‘>’ token
         const volScalarField& summass = p.mesh().lookupObject<volScalarField>("summass");
                                                                             ^
lnInclude/ParticleCollector.C:750:31: error: ‘const volScalarField {aka const class Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}’ has no member named ‘primitiveFieldRef’
         scalar& PPC = summass.primitiveFieldRef() [cellIp]; 
                               ^
lnInclude/ParticleCollector.C: In member function ‘virtual void Foam::ParticleCollector<CloudType>::postMove(Foam::ParticleCollector<CloudType>::parcelType&, Foam::label, Foam::scalar, const point&, bool&)’:
lnInclude/ParticleCollector.C:748:77: error: expected primary-expression before ‘>’ token
         const volScalarField& summass = p.mesh().lookupObject<volScalarField>("summass");
                                                                             ^
lnInclude/ParticleCollector.C:750:31: error: ‘const volScalarField {aka const class Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}’ has no member named ‘primitiveFieldRef’
         scalar& PPC = summass.primitiveFieldRef() [cellIp];

不知道从哪里入手，麻烦李老师了

李东岳

@chengan-wang 在如何在src文件中创建标量物理filed并且调用中说：

mesh().lookupObject<volScalarField>

p.mesh().objectRegistry::lookupObject<volScalarField>("summass")
这个可以么

chengan.wang

@李东岳在如何在src文件中创建标量物理filed并且调用中说：

p.mesh().objectRegistry::lookupObject<volScalarField>("summass")

lnInclude/ParticleCollector.C: In member function ‘virtual void Foam::ParticleCollector<CloudType>::postMove(Foam::ParticleCollector<CloudType>::parcelType&, Foam::label, Foam::scalar, const point&, bool&)’:
lnInclude/ParticleCollector.C:750:31: error: ‘const volScalarField {aka const class Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}’ has no member named ‘primitiveFieldRef’
         scalar& PPC = summass.primitiveFieldRef() [cellIp]; 
                               ^
lnInclude/ParticleCollector.C: In member function ‘virtual void Foam::ParticleCollector<CloudType>::postMove(Foam::ParticleCollector<CloudType>::parcelType&, Foam::label, Foam::scalar, const point&, bool&)’:
lnInclude/ParticleCollector.C:750:31: error: ‘const volScalarField {aka const class Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}’ has no member named ‘primitiveFieldRef’
         scalar& PPC = summass.primitiveFieldRef() [cellIp];

好像是没有 member named ‘primitiveFieldRef’，我用的是of2.4.0

李东岳

scalar& PPC = summass.internalField() [cellIp]

chengan.wang

@李东岳麻烦李老师了。

        const volScalarField& summass = p.mesh().objectRegistry::lookupObject<volScalarField>("summass");
        const label cellIp = p.cell(); 
        scalar& PPC = summass.internalField()[cellIp]; 
        PPC += p.nParticle()*p.mass();

还是有错误如下

lnInclude/ParticleCollector.C:750:53: error: binding ‘const double’ to reference of type ‘Foam::scalar& {aka double&}’ discards qualifiers
         scalar& PPC = summass.internalField()[cellIp];

李东岳

scalar& PPC = summass.internalFieldRef() [cellIp]

chengan.wang

@李东岳李老师，又出现了之前的问题：

lnInclude/ParticleCollector.C: In member function ‘virtual void Foam::ParticleCollector<CloudType>::postMove(Foam::ParticleCollector<CloudType>::parcelType&, Foam::label, Foam::scalar, const point&, bool&)’:
lnInclude/ParticleCollector.C:750:31: error: ‘const volScalarField {aka const class Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}’ has no member named ‘internalFieldRef’
         scalar& PPC = summass.internalFieldRef() [cellIp]; 
                               ^
lnInclude/ParticleCollector.C: In member function ‘virtual void Foam::ParticleCollector<CloudType>::postMove(Foam::ParticleCollector<CloudType>::parcelType&, Foam::label, Foam::scalar, const point&, bool&)’:
lnInclude/ParticleCollector.C:750:31: error: ‘const volScalarField {aka const class Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}’ has no member named ‘internalFieldRef’
         scalar& PPC = summass.internalFieldRef() [cellIp];

我用grep查了一下，有些函数中用的是internalField()而查不到internalFieldRef()，但如果采用

scalar& PPC = summass.internalField() [cellIp]

错误如下

lnInclude/ParticleCollector.C:750:53: error: binding ‘const double’ to reference of type ‘Foam::scalar& {aka double&}’ discards qualifiers
         scalar& PPC = summass.internalField()[cellIp];

李东岳

2.4版本太老了，我也忘了怎么用了，

scalar PPC = summass.internalField() [cellIp]

这个呢

chengan.wang

@李东岳编译倒是能通过了，但是运行算例log文件显示错误

[12] 
[12] 
[12] --> FOAM FATAL ERROR: 
[12] 
    request for volScalarField summass from objectRegistry region0 failed
    available objects of type volScalarField are

19
(
thermo:mu
thermo:psi
K
h
kappa
air
rho
k
Cp
dpdt
thermo:psi_0
alphat
dQ
p
T
mut
H2O
epsilon
thermo:alpha
)
[12] 
[12] 
[12]     From function objectRegistry::lookupObject<Type>(const word&) const
[12]     in file /home/chengan/OpenFOAM/OpenFOAM-2.4.0/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.
[12]

李东岳

应该是你summass这个代码还没创立的时候就寻找了。把这个代码提前几行

chengan.wang

@李东岳我把

const volScalarField& summass = p.mesh().objectRegistry::lookupObject<volScalarField>("summass");

挪到程序前面了，还是同样的错误

李东岳

@chengan-wang 在如何在src文件中创建标量物理filed并且调用中说：

volScalarField summass
(
IOobject
(
"summass",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar ("zero", dimMass, 0.0)
);

把这个提前

chengan.wang

@李东岳李老师，按照下列代码

        const label cellIp = p.cell(); 
        scalar PPC = summass.internalField() [cellIp]; 
        PPC += p.nParticle()*p.mass();

编译通过了，也能算了，但是好像没有赋值进去，结果是0

FoamFile
{
    version     2.0;
    format      binary;
    class       volScalarField;
    location    "0.0002";
    object      summass;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    walls
    {
        type            calculated;
        value           uniform 0;
    }
}

Tens

@chengan-wang 在如何在src文件中创建标量物理filed并且调用中说：

    const label cellIp = p.cell(); 
    scalar PPC = summass.internalField() [cellIp]; 
    PPC += p.nParticle()*p.mass();

这样只有PPC这个量在变，summass不跟着变。
我感觉如果是在src库中进行调用和操作，还是应该在库中创建这个量

chengan.wang

@tens 应该是这样的，感觉应该带个指针啥的

李东岳

scalar& PPC = summass.internalField().ref() [cellIp]

CFD中文网

如何在src文件中创建标量物理filed并且调用