CFD中文网

@李东岳 就是想把离散点数据

    -17.00000000    -59.58333206      0.00000000    302.77685547
    -17.00000000    -62.50000000      0.00000000    302.77685547
    -17.00000000    -65.41666412      0.00000000    302.77685547
    -17.00000000    -68.33333588      0.00000000    302.77685547
    -17.00000000    -71.25000000      0.00000000    302.77685547
    -17.00000000    -74.16666412      0.00000000    302.77682495
    -17.00000000    -77.08333588      0.00000000    302.77679443
    -17.00000000    -80.00000000      0.00000000    302.77664185

整理成这样：

VARIABLES = "X""Y""Z"
 Nodes=39888, Elements=78526, ZONETYPE=FETriangle
    -17.00000000    -59.58333206      0.00000000    302.77685547
    -17.00000000    -62.50000000      0.00000000    302.77685547
    -17.00000000    -65.41666412      0.00000000    302.77685547
    -17.00000000    -68.33333588      0.00000000    302.77685547
    -17.00000000    -71.25000000      0.00000000    302.77685547
    -17.00000000    -74.16666412      0.00000000    302.77682495
    -17.00000000    -77.08333588      0.00000000    302.77679443
    -17.00000000    -80.00000000      0.00000000    302.77664185
       1       2       3
       3       4       1
       5       6       2
       2       1       5
       7       8       6
       6       5       7

@Prometheus10 谢谢哈，我主要目的是建立节点数据之间的关联，三角形，缩放只能是个视觉上的平面。

用paraview或tecplot那能不能把曲面上的数据投影到一个规则的矩形平面上呢？

大家好，
我目前用starccm做了一个仿真，数据能用tecplot导出ZONETYPE=FEPolygon形式，如何进一步转换成ZONETYPE=FETriangle？
谢谢。

大家好，目前我可以输出温度场的三维数据，比如

我想垂直于z轴做100个等间距垂直切面，在每个切面上对温度求积分并取平均值。
\begin{equation}
\int_{A} T dx dy /A
\end{equation}
然后沿着z轴画出一条温度平均值曲线。

谢谢大家！

@李东岳 非常感谢李老师，终于搞定，能输出数据了，但还需要进一步验证是否正确。差点放弃想升级新版本of。

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

希望对以后需要的朋友有点借鉴作用。特别感谢 @星星星星晴 不厌其烦的解答

@李东岳 李老师，这句命令

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

又报错了

lnInclude/ParticleCollector.C:745:71: error: ‘class Foam::objectRegistry’ has no member named ‘lookupObjectRef’
     volScalarField& summass = p.mesh().objectRegistry::lookupObjectRef<volScalarField>("summass");

@李东岳 李老师，还是不行

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

编译出错

lnInclude/ParticleCollector.C:749:38: error: assignment of read-only location ‘(&(& summass)->Foam::GeometricField<Type, PatchField, GeoMesh>::internalField<double, Foam::fvPatchField, Foam::volMesh>())->Foam::Field<double>::<anonymous>.Foam::List<double>::<anonymous>.Foam::UList<T>::operator[]<double>(((Foam::label)cellIp))’
     summass.internalField() [cellIp] += p.nParticle()*p.mass();
                                      ^
lnInclude/ParticleCollector.C: In instantiation of ‘void Foam::ParticleCollector<CloudType>::postMove(Foam::ParticleCollector<CloudType>::parcelType&, Foam::label, Foam::scalar, const point&, bool&) [with CloudType = Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > >; Foam::ParticleCollector<CloudType>::parcelType = Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> >; Foam::label = int; Foam::scalar = double; Foam::point = Foam::Vector<double>]’:
parcels/derived/basicKinematicCollidingParcel/makeBasicKinematicCollidingParcelSubmodels.C:53:1:   required from here
lnInclude/ParticleCollector.C:749:38: error: assignment of read-only location ‘(&(& summass)->Foam::GeometricField<Type, PatchField, GeoMesh>::internalField<double, Foam::fvPatchField, Foam::volMesh>())->Foam::Field<double>::<anonymous>.Foam::List<double>::<anonymous>.Foam::UList<T>::operator[]<double>(((Foam::label)cellIp))’
     summass.internalField() [cellIp] += 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:746:43: error: ‘const InternalField {aka const class Foam::Field<double>}’ has no member named ‘ref’
     scalar& PPC = summass.internalField().ref() [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:746:43: error: ‘const InternalField {aka const class Foam::Field<double>}’ has no member named ‘ref’
     scalar& PPC = summass.internalField().ref() [cellIp]; 

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

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

        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;
    }
}

@李东岳 我把

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

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

@李东岳 编译倒是能通过了，但是运行算例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] 

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

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]; 

@李东岳 麻烦李老师了。

        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]; 

@李东岳 在 如何在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

@李东岳 谢谢李老师帮忙！
我我的是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]; 

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

@李东岳 感谢李老师回复。
我先在求解器中的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等文件中吧？！瞎猜的。

李老师能指点一下吗？

大家好，请教一个问题：
我知道如何在求解器中的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等库文件？

谢谢大家。

@星星星星晴 在 Euler-Lagrange的一些解析，sprayFoam 中说：

@chengan-wang

1. 数据量极大就是这个方法最大的问题，没办法。
2. 你的mass3D没有输出，在结束这一步的迭代以后，在下一步迭代会被清空的。
   你这一步已经很靠近方法1了，你不要用什么mass3D这个变量，你在solver中建立一个field， 比如叫summass，然后在cloudfunction的postmove 这个function中直接通过，就可以对一下field进行赋值了。

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

我的技能也是有限，所以没有办法给你解决所有问题，你可以自己试试，找找资料，如果发现更好的办法，也希望你分享出来。

对于方法2，我个人有个想法，你获得你所要区域的cell number， 然后仅仅对这些cell进行输出也可以的。我自己没做过，大家都是摸着石头过河，可能每个人的解决方法并不通用，很case- sensitive，这也是开源的一个弊端。。

祝你好运～

谢谢。我在createFields.H建立了一个field:

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

我没有理解“在cloudfunction的postmove 这个function中直接通过”是啥意思呢？麻烦了