OpenFOAM小代码

李东岳

meanDiameter
    {
        type            coded;
        libs            ("libutilityFunctionObjects.so");
        name            error;
    
        codeExecute
        #{
            const volScalarField& d =
                mesh().lookupObject<volScalarField>("d.alpha.oil");

            scalar d32 = d.weightedAverage(mesh().V()).value();
    
            if (Pstream::master())
            {
                std::ofstream file;
                file.open ("d32", std::ofstream::out | std::ofstream::app);
                file << mesh().time().timeName() << " " << d32 << "\n";
                file.close();
            }
        #};
    }

我最近也凑巧要用这个，我用的上面这种方式。

五好青年

@李东岳
请问东岳老师

(1) 东岳老师，您上面写的：codedFixedValue
(2) yfclark朋友，上面写的：CodeStream

在自定义边界条件方面，这两个有啥区别啊，我研究了几天没想明白

五好青年

@yfclark
朋友，您好：

codedFixedValue和CodeStream

你知道在自定义边界条件方面，这两个有啥区别啊，我研究了几天没想明白

李东岳

codedFixedValue更简单，直接用

code
        #{
            operator==(min(10, 0.1*this->db().time().value()));
        #};

就可以赋值。CodeStream更普适性，我见过一些代码用CodeStream对内部场进行初始化。比如建立U的内部场与空间的函数。

hongjiewang

@李东岳目前我需要将温度边界和压力边界定义为随相含量变化的一个场。我只可以做到变化之后的也是同一类边界条件，但是像下图中的气相区压力是第一类边界条件，液相区温度是第二类边界条件。请问这种变边界类型的边界应该怎么定义~ 1632279639(1).png

hongjiewang

@hongjiewang 在 OpenFOAM小代码中说：

@李东岳目前我需要将温度边界和压力边界定义为随相含量变化的一个场。我只可以做到变化之后的也是同一类边界条件，但是像下图中的气相区压力是第一类边界条件，液相区温度是第二类边界条件。请问这种变边界类型的边界应该怎么定义~

已解决压力的~但是温度的第二类和第三类的转换还没有解决~

疏影横斜水清浅

你好，请问const polyPatch& cPatch = mesh.boundaryMesh()[patchID]；与const fvPatch& cPatch = mesh.boundary()[patchID];有什么不同？

Dingcy

@李东岳在 OpenFOAM小代码中说：

//const fvBoundaryMesh& boundaryMesh = boundaryPatch.boundaryMesh();
//const fvMesh& mesh = boundaryMesh.mesh();

我想用codedFixedValue,写一个壁面的边界条件，壁面都是wall，而不是patch，所以这里想用这两个命令，但是编译过程中总是提示我是不是应该声明点啥呢？我也查了网上的，发现全是在inlet使用这个功能，然后入口是patch

李东岳

@dingcy

const polyPatch& cPatch = mesh.boundaryMesh()[patchID]

类似这个？

Dingcy

@李东岳好像不是这个又显示了mesh没有声明....

vbcwl

@李东岳能请教下这个具体怎么使用么？我没能输出utau

李东岳

@vbcwl 上面那些代码么。上面那些代码都需要自己写到代码里面去。还不是简单的使用。需要有编程基础，这些代码起到一定的提示作用。

2019201300

李老师您好，请问输出努塞尔数代码是直接放进去吗？还是里边的部分数据需要进行修改呢微信截图_20221112113853.png ？

李东岳

@2019201300 这个我太久之前写的。已经忘了。但是这个代码只是提供案例。具体数值大概率需要自己改一下。直接放在controlDict下面。

2019201300

@李东岳好的，谢谢李老师，我再认真研究研究。

疏影横斜水清浅

请问这个代码怎么修改，才可以使用网格的体积,

volScalarField muc
(
    IOobject
    (
        IOobject::groupName("mu", continuousPhaseName),
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    rhoc*laminarTransport.nu()
);

volScalarField Rep("Rep", dp/laminarTransport.nu()*mag(Up - U));
volScalarField CDRe("CDRe", 24/Rep+3.6*pow(Rep,-0.313));
scalar deltaTs = runTime.deltaTValue();

Up = Up + deltaTs*(3/4*muc/rhop/pow(dp,2.0)/mesh.V()*CDRe*(U - Up) + g);

其中的mesh.V()怎么修改才可以。

2019201300

@疏影横斜水清浅参考下下边这个代码吧。

coded
 {
     libs ( "libutilityFunctionObjects.so" );
     enabled         true;
     type            coded;
     name           printMinU;
     executeControl  timeStep;
     executeInterval 200;
     writeControl   timeStep;
     writeInterval  200;

     codeOptions
     #{
         -I$(LIB_SRC)/meshTools/lnInclude
     #};

     codeExecute
     #{
         const volScalarField& T
         (
             mesh().lookupObject<volScalarField>("T")
         );

         const fvPatchList& patches = mesh().boundary();

             std::ofstream file;
             file.open ("nus", std::ofstream::out | std::ofstream::app);
         forAll(patches, patchi)
         {
             const fvPatch& currPatch = patches[patchi];
                 if (currPatch.name() == "bottomWall")
                 {
                     fvPatchScalarField nus = T.boundaryField()[patchi];

                     scalar L = 0.1;
                     scalar Timp = 343;
                     const volScalarField& T = mesh().lookupObject<volScalarField>("T");
                     scalarField Tinternal = T.internalField();
                     scalar sumT = gSum(Tinternal * mesh().V());
                     scalar V = gSum(mesh().V());
                     scalar avgT = sumT / V;
                     nus = T.boundaryField()[patchi].snGrad()*L/(Timp - avgT);
                     //nus2file << "Time: " << mesh().time().timeName() << " " << nus << nl;
         
                     forAll(T.boundaryField()[patchi], facei)
                     {
                         //Pout << mesh().C().boundaryField()[patchi][facei].x()/0.1
                           //   << " " << nus[facei] << nl;
                         file << "Time: " << mesh().time().timeName() << " " << mesh().C().boundaryField()[patchi][facei].x() << " " << mesh().C().boundaryField()[patchi][facei].z() << " " << nus[facei] << std::endl;
                     }
                     
                 }
         }
         file.close(); 
    #};
 }

疏影横斜水清浅

@2019201300 请问，这个代码是能够将每个网格的体积存储在Vcell中吧？

volScalarField Vcell
(
    IOobject
    (
        "Vcell",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedScalar("", dimVolume, 0.0)
);

forAll(Vcell ,i)
{
    Vcell[i] = mesh.V()[i];
}

李东岳

本帖之前写的代码，挑主要的都更新到这里的

http://dyfluid.com/code.html

WangXin

@2019201300 我想问一下，如果我要计算整个壁面的平均努塞尔数，该怎么处理呢

CFD中文网

OpenFOAM小代码

本帖之前写的代码，挑主要的都更新到这里的