@星星星星晴 在 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中直接通过”是啥意思呢？麻烦了

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

@chengan-wang 这也简单啊，同理，你不需要判断是否碰到面，而是做一个判断是否在你的区域内即可，比如你需要判断在两个面之间，z-dir

if(p.position(2)<1 & p.position(2)>0)
{
输出
}

但是这种方法会出现什么问题，比如一个parcel在这个区域待了3个$\Delta t_L$, 那就输出3次，就是你的输出文件几何倍数的增长。。

您好，我按照您的提示，思路如下：
先是定义了三维空间的范围，以及坐标质量的数组。

        scalar xmin = -0.01;
        scalar xmax = 0.01;
        scalar ymin = 0;
        scalar ymax = 0.1;
        scalar zmin = -0.01;
        scalar zmax = 0.01;
        scalar dxyz = 5e-3;
        const int nx = floor((xmax - xmin)/dxyz);
        const int ny = floor((ymax - ymin)/dxyz);
        const int nz = floor((zmax - zmin)/dxyz);
        int nxc;
        int nyc;
        int nzc;
        scalar xc[nx];
        scalar yc[ny];
        scalar zc[nz];
        scalar mass3D[nx][ny][nz] = {0};

我只想最后时刻输出结果，所以先用if条件判断是否要输出数据

    if ((this->owner().db().time().value()) == 0.02)
    {

        for (int i=0; i<nx; i++)
        {
            xc[i] = i*dxyz+xmin;
        }

        for (int j=0; j<ny; j++)
        {
            yc[j] = j*dxyz+ymin;
        }

        for (int k=0; k<nz; k++)
        {
            zc[k] = k*dxyz+zmin;
        }

        nxc = floor((p.position().component(0) - xmin)/dxyz);
        nyc = floor((p.position().component(1) - ymin)/dxyz);
        nzc = floor((p.position().component(2) - zmin)/dxyz);

        mass3D[nxc][nyc][nzc] += p.nParticle()*p.mass();


        std::ofstream outfile;
        outfile.setf(ios_base::fixed, ios_base::floatfield);
        outfile.precision(8);
        outfile.open("postProcessingWCA/Mass3D.txt", ios_base::app);

        for (int i=0; i<nx; i++)
        {
            for (int j=0; j<ny; j++)
            {
                for (int k=0; k<nz; k++)
                {
                    outfile << xc[i] << tab << yc[j] << tab << zc[k] << mass3D[i][j][k] << tab << tab  << nl;
                }
            }
        }
    }

之后我的想法就是根据粒子的位置计算索引位置 nxc nyc nzc，带入到mass3D[nxc][nyc][nzc] 相应的位置。
但是存在两个问题：

1. 每一个粒子统计之后都输出到了Mass3D.txt，数据量非常大，而我是想把所有的粒子质量统计到相应的空间点上累积；
2. 另外就是定义方式scalar mass3D[nx][ny][nz] = {0}; 以及累积方法mass3D[nxc][nyc][nzc] += p.nParticle()*p.mass();感觉有问题，从结果上看感觉没有把相同位置的粒子质量加到一起。

想麻烦您帮忙看看，谢谢

@李东岳 我找到如何输出endTime方法了，用this->owner().db().time().endTime().value()就可以。但是放在判断语句中出现如下错误：

note:   ‘Foam::word’ is not derived from ‘const std::fpos<_StateT>’
     if ((this->owner().db().time().timeName()) == (this->owner().db().time().endTime().value()))

/usr/include/c++/5/bits/stl_iterator.h:349:5: note: candidate: template<class _IteratorL, class _IteratorR> bool std::operator==(const std::reverse_iterator<_Iterator>&, const std::reverse_iterator<_IteratorR>&)
     operator==(const reverse_iterator<_IteratorL>& __x,
     ^
/usr/include/c++/5/bits/stl_iterator.h:349:5: note:   template argument deduction/substitution failed:

李老师能指点一下么？

@李东岳 谢谢回复。
我要输出一个物理场，网格数量挺大，200x100x100，要是每时每刻都输出的话，数据量很大，所以我需要加一个判断，只输出我感兴趣的时刻就行

大家好，
我在程序中添加了一段代码如下：

    /*std::ofstream outfile;
    outfile.setf(ios_base::fixed, ios_base::floatfield);
    outfile.precision(8);
    outfile.open("qq.dat", ios_base::app);
    outfile << this->owner().db().time().timeName() << tab
            << "test"                               << tab
            <<  nl;*/

这样的话每次程序执行到这段代码时候都会输出到qq.dat文件，如果运行时间很长，文件会非常大，我又不需要输出所有时间点的数据。如果我想输出特定时间点的数据，比如最后一步，用if条件语句判断，条件式子没想好该怎么写，怎么才能方便地提取controlDict中的endTime呢？

startTime       0;

stopAt          endTime;

endTime         2e-04;

比如写成if( this->owner().db().time().timeName()==endTime)

大家好、各位大佬好，咨询一个编译函数调用问题
我在别的论坛找到了一个在求解器中定义函数编译调用的例子如下：

//fun.H
#ifndef fun_H
#define fun_H

#include "scalar.H"
#include "vector.H"

namespace Foam
{
    vector TestMe(const vector& myvec);
}
#endif

以及

//fun.C
#include "fun.H"

Foam::vector Foam::TestMe(const vector& myvec)
{
    Info<< "Hello from TestMe" << endl;
    Info<< "vector: " << myvec << endl;
}

然后在主函数里面调用 "fun.H"就行了。

#include "fun.H"

using namespace Foam;

int main(int argc, char *argv[])
{
    vector A(1,2,3);
    Info<< "Hello from the main()" << endl;
    return 0;
}

我的思路就是直接在fun.H中写我想调用的代码。比如在“KinematicCloudI.H”文件中定义的

template<class CloudType>
inline const Foam::tmp<Foam::volVectorField>
Foam::KinematicCloud<CloudType>::UPar() const
{
    tmp<volVectorField> tUPar
    (
        volVectorField::New
        (
            this->name() + ":UPar",
            mesh_,
            dimensionedVector(dimVelocity, vector(0,0,0)),
            extrapolatedCalculatedFvPatchVectorField::typeName
        )
    );


    volVectorField& UPar = tUPar.ref();
    forAllConstIter(typename KinematicCloud<CloudType>, *this, iter)
    {
        const parcelType& p = iter();
        const label celli = p.cell();
        UPar[celli] += (p.U())/(p.nParticle());
    }
    UPar.primitiveFieldRef() /= 1;//mesh_.V();
    UPar.correctBoundaryConditions();

    return tUPar;
}

我直接移植过来，编译出现了如下错误：

fun.H:16:41: error: no ‘const Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::KinematicCloud<CloudType>::UPar() const’ member function declared in class ‘Foam::KinematicCloud<CloudType>’
 Foam::KinematicCloud<CloudType>::UPar() const

我尝试把“KinematicCloudI.H”里面的头文件粘上去也通过不了，应该怎么修改以及加一些什么内容呢？
谢谢

@星星星星晴 非常感谢，我再仔细研究一下。我在google以及源代码中没有查到IVT8_StochasticCollisionModel_2022，PPC，它们代表的意思您方便告诉我么？感觉这个PPC应该是粒子的相关信息，又感觉直接用这个名字不对

@星星星星晴 不是要您的code，就是想知道应该改哪个文件而已，您误会了

@星星星星晴 我仔细想了一下，输出文件的方法可能数据量很大，感觉还是转换的方法更合适，
但还是有点没有搞清楚。

您之前告诉我可以放到cloudfunction中，我想知道您具体是放在哪个文件中了，能提供个思路么？我想看看安装包的原始文件，然后再比较一下您的代码，这样可能就清楚该怎么改和调用了。

另外，IVT8_StochasticCollisionModel_2022，PPC这是随便取得名字么？我在使用的时候是不是得改动呢？

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

@chengan-wang 这也简单啊，同理，你不需要判断是否碰到面，而是做一个判断是否在你的区域内即可，比如你需要判断在两个面之间，z-dir

if(p.position(2)<1 & p.position(2)>0)
{
输出
}

但是这种方法会出现什么问题，比如一个parcel在这个区域待了3个$\Delta t_L$, 那就输出3次，就是你的输出文件几何倍数的增长。。

您好，如果加判断条件的话应该是在particleCollector代码中加，然后编译，对吧？这个圆柱面或者矩形面的定义就无所谓了。也就是说单独命名一个输出文件。
还有，我的网格数量挺大，200x100x100，我得定义好多个长方体区域，采用循环的办法，判断粒子与每一个区域的相对位置，如果符合条件就保存数据，跳出循环。但如果我的计算时间很长，每个时间点都保存数据，数据量肯定很大，能否只保存最后一个时间步的数据？也就是说我需要传递controlDict的参数，比如

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         5;

deltaT          0.01;