https://github.com/OpenFOAM/OpenFOAM-6/blob/master/src/TurbulenceModels/turbulenceModels/linearViscousStress/linearViscousStress.C 这里面计算的散度

template<class BasicTurbulenceModel> Foam::tmp<Foam::fvVectorMatrix> Foam::linearViscousStress<BasicTurbulenceModel>::divDevRhoReff ( volVectorField& U ) const { return ( - fvc::div((this->alpha_*this->rho_*this->nuEff())*dev2(T(fvc::grad(U)))) - fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U) ); }

各位老师好，我最近研究了一下reactingmultiphaseEulerFoam，对其中的表面张力的计算有一些困惑：
在pEqn.H中定义了：

phigFs.set ( phasei, ( alpharAUfs[phasei] *( ghSnGradRho - (fvc::interpolate(phase.rho() - rho))*(g & mesh.Sf()) - fluid.surfaceTension(phase)*mesh.magSf() ) ).ptr() );

最后通过该式更新各项的速度：

mSfGradp = pEqnIncomp.flux()/rAUf; forAll(phases, phasei) { phaseModel& phase = phases[phasei]; phase.U() = HbyAs[phasei] + fvc::reconstruct ( alpharAUfs[phasei]*mSfGradp - phigFs[phasei] ); phase.U().correctBoundaryConditions(); fvOptions.correct(phase.U()); }

我个人推导的公式为：
78b21ff7-7858-4542-a777-3f2e993f3957-image.png
代码里面实现的公式为：
efa31825-beb7-4ee1-8483-980d46400369-image.png
表面张力项多了一个alpha，请大佬看看，不知道是我推导错了还是代码有bug
@东岳
reactingtwophaseEulerFoam中的式（20)和式（21）有个地方符号有问题：
d07eb439-d134-41ab-862d-944dc3a3efd9-image.png
最后压力项前应该为$-\alpha_{d}$

多次尝试，终于搞定了，打扰了

80a6ed67-6d71-4c4f-9c0b-de7dd78ac863-image.png
如图左侧进口压力为0梯度，来流速度为抛物线分布，上下两边为0，中间最大。中点处速度从零开始平缓增大至最大值，其随时间的变化规律如下式，请问这种情况如何设置0文件夹下的边界条件呢？谢谢
4956dd8a-7642-4853-9570-74be9dab8073-image.png

自问自答一波:
我现在感觉好像是检查constant/porosityProperties 文件里面的字典. 只要里面存在字典 cellZone(键) 对应的 name_(值)的block就被当作是多孔介质, 然后通过

cellZoneID = mesh.cellZone().indices(name_)

这个函数来确认该区域的ID
上述内容是在
porosityModel.H porosityModel.C中感觉到的.但是没有找到实质的证据.

@东岳 明白了，谢谢东岳老师

@王金成 我用的1812，可以编译SWAK，不过需要自己做些配置，不然会报错，过段时间我写个教程。

@浪迹天大 谢谢您！ 问题解决了

@东岳 把基类重写了一遍，实现了我想要的功能，虽然麻烦，但是相当于在Cloud里面添加功能，对原有的代码结构基本没有影响，还算不错的思路。

@qjh888 您好，不知道您如何更新的solver？方便讲一下吗？最近我也在做相关的事情。

各位CFD朋友，目前我想学习一下firefoam，但是baidu搜到的资料实在有限，如果方便的话能否提供一些算例资料给我，谢谢大家。

在求解可压缩流动时，出现了如下情况：导致最终结果发散，可能是哪个参数的问题呢？

rhoEqn max/min : 1.21419 0

@东岳 请问东岳老师，可以用速度基求解器比如rhoCentralFoam求解这个算例么？还得加个温度去解压力？
没算过DNS，望指点！谢谢！

@lllwonderliquid 已更新simpleFoam解析 http://dyfluid.com/simplefoam.html

已经解决了。谢谢您。可能是自己安装的经验不够吧。

滤镜-time annodate
即可

:tishizi:

谢谢，谢谢。

@程迪 程老师您好，请问您知道如何把Prt定义成某个表达式吗？

东岳老师好，大家好：
我是正在学习OpenFOAM的新手
最近在学习模拟湍流时运行pisoFoam下自带算例porousBlockage时出现错误：

Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes No MRF models present Creating finite volume options from "constant/fvOptions" Selecting finite volume options model type explicitPorositySource Source: porosity1 - selecting cells using cellZone porousBlockage --> FOAM FATAL ERROR: Cannot find cellZone porousBlockage Valid cellZones are 0() From function void Foam::fv::cellSetOption::setCellSet() in file cellSetOption/cellSetOption.C at line 142. FOAM exiting

请问是什么原因呢？

此外，还想问一下，如果使用标准k-epsilon模型（RAS算例下）模拟湍流，进出口使用什么条件比较合适呢？