首选请明确是否调用了湍流模型！

Larsen证明了两方程模型会导致准势流区域的湍动能会指数增长，波高沿程衰减。

建议调用Larsen改进的stablized版本的额湍流模型。

求大佬帮我看一下，跪求

@李东岳
当边界不做判断直接叠加之后是这样的
fc3c9e32-1d39-4a41-b44a-85c2594da99b-image.png

bottom
{
type codedFixedValue;
value uniform 101325;
name face1P;
code #{
scalarField& field = *this;
const label id = this->patch().boundaryMesh().findPatchID("bottom");//在网格中找到inlet
const volScalarField& pp = this->db().objectRegistry::lookupObject<volScalarField>("p");
const scalarField& ppp=pp.boundaryField()[id].patchInternalField();
forAll(field, faceI)
{

field[faceI]=ppp[faceI]+100; } #}; }

@peng23 https://cfd-china.com/topic/4499/openfoam第三方造波包waves2foam安装 你看看这个。waves2foam要注意版本一致性的问题。

另外，GLIBCXX_3.4.29' not found应该是你的ubuntu系统太高。

感觉跟setExprFieldsDict的语法一样

https://www.cfd-china.com/topic/4627/后处理时候singlegraph可以提取两条线的值吗

在solids4foam中运行验证算例（孤立波下斜坡上的直墙变形和受力），但是当波浪到固体交界面时，就会出现浮点报错问题，修改了网格和固体交界面的问题，还是没有改善这个问题。请问有大佬们知道怎么解决这个问题吗？
报错时间步alpha场：（每次到这个地方就出现报错）
2a67acf0-d296-4f56-82eb-f8a159853cca-image.png

报错部分log文件：

Time = 8.933, iteration: 9 Current fsi under-relaxation factor (Aitken) of interface: 3.24861e-06 Maximal accumulated displacement of interface 0: 2.02709e-06 GAMG: Solving for cellMotionUx, Initial residual = 0.27053, Final residual = 4.88471e-05, No Iterations 2 GAMG: Solving for cellMotionUz, Initial residual = 0.276424, Final residual = 7.60589e-05, No Iterations 2 GAMG: Solving for cellMotionUx, Initial residual = 5.19396e-05, Final residual = 4.62925e-07, No Iterations 4 GAMG: Solving for cellMotionUz, Initial residual = 8.68504e-05, Final residual = 8.18649e-07, No Iterations 3 Evolving fluid model: wave2solids PIMPLE: iteration 1 MULES: Solving for alpha1 alpha1 BC on patch inlet alpha1 BC on patch inlet Liquid phase volume fraction = 0.259489 Min(alpha1) = -2191.92 Max(alpha1) = 2360.17 MULES: Solving for alpha1 alpha1 BC on patch inlet alpha1 BC on patch inlet Liquid phase volume fraction = 0.259489 Min(alpha1) = -158062 Max(alpha1) = 163488 Velocity BC on patch inlet 3D_2D Absorption BC on patch outlet "Correction Levels" 1( 0 ) DICPCG: Solving for pd, Initial residual = 0.10778, Final residual = 0.00369688, No Iterations 8 DICPCG: Solving for pd, Initial residual = 6.32049e-06, Final residual = 2.94799e-07, No Iterations 644 DICPCG: Solving for pd, Initial residual = 4.65247e-05, Final residual = 2.1618e-06, No Iterations 637 DICPCG: Solving for pd, Initial residual = 2.0431e-07, Final residual = 9.74494e-08, No Iterations 110 Velocity BC on patch inlet 3D_2D Absorption BC on patch outlet "Correction Levels" 1( 0 ) DICPCG: Solving for pd, Initial residual = 4.64318e-05, Final residual = 2.21746e-06, No Iterations 638 DICPCG: Solving for pd, Initial residual = 2.16117e-07, Final residual = 8.27667e-08, No Iterations 132 Velocity BC on patch inlet 3D_2D Absorption BC on patch outlet "Correction Levels" 1( 0 ) Setting traction on solid interfaces Mapping face values using directMap Total force on fluid interface 0: (4.3321e+10 -2.82466e-19 -0.673933) Total force on solid interface 0: (-4.3321e+10 2.82466e-19 0.673933) Evolving solid solver Solving the updated Lagrangian form of the momentum equation for DD -------------------------------------------------------------------------- mpirun noticed that process rank 6 with PID 168600 on node node03 exited on signal 8 (Floating point exception).

@李东岳 在 piso和pimple 中说：

数值稳定性

李老师您好，我在用buoyantPimple求解器的PSIO模式计算流场压力、密度和速度出现明显的数值振荡，就是随时间步长周期性振荡发展，但是振幅确实是一直在减小的，可否理解为这就是数值振荡问题，通过调整PIMPLE求解器外部动量预测修正次数转化为PIMPLE增强稳定性？
计算时序图如下：
8055ef71-35b7-40e6-b692-ac6aa27a4f2c-image.png
58fdb9e2-331b-4f6b-9a4a-39cb7b94f1bf-image.png
064a77b7-1f57-4077-9480-e6395c761979-image.png

paraview以列表形式导出流场所有数据时导出太慢（每个时间步的数据比较小时无影响）。如果用自动写每个时间步的数据时，则每个时间步的数据都不全。请问各位有遇到过吗？有无解决办法。

学流体的小明 在 压力基求解器在OpenFOAM中的植入问题 中说：

CaseB中，同样的方法，压力p基本都是10000，或者有是9999.99，这和我限制了writePrecision 6有关系。

所以输出流场的时候，最好以binary二进制的格式输出。
后处理的文件，由于要使用其它软件读取，就用ASCII。

最近也在折腾做旋转对称圆柱型计算域，找到个这：https://www.fangzhenxiu.com/post/1826059?code=031TTpFa1IW6xG0skbHa1OsJy23TTpFt&state=afterWeixin
根据这个操作画出来的网格虽然是楔形块，但是楔面并不参照哪一坐标平面对称，依然不能用在OpenFOAM上。
又找到一个这个：https://www.fangzhenxiu.com/post/1829359
作者好像通过两次旋转实现了对称的楔面，但我没有复现出来。
最后还是打算用blockMesh来画了。

@ShaneHEEE 是的，和出口设置一样，就可以了

@HITSC30 你好，请问你这个问题解决了嘛？我也遇到了和你一样的报错。

您能解释一下为什么会有rhoCpPhi=rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2吗？这个是什么物理意义？

类似的案例，我也遇到过，我的案例是在跑别人的代码时出错的，在出现的“No MRF models present/n No finite volume options present"报错，且后面出现了”Kill“字祥时，其中“Kill”字样很重要。
到后面我发现：我的案例文件放错的位置，原先我放在与OF整体文件的所在的主文件夹内，我的解决是通过将案例放在OF的主文件夹内，而不是在OF主文件夹外，这成功。
首先说我为什么这么做：

出现MRF警告时，我也以为时MRF链接出了问题，但事实是其他案例是正常的，故排除MRF问题； 我看到所引用的案例包内有：platfrom文件且内部为”linux64GccDPint32Opt“，这与OF主文件一致，我猜测事路径出了问题，最终事实如此。

@ann 在 chemkin 中说：

在网上找到的C12H26化学反应机理文件，不管是直接运行还是通过chemkinToFoam转换都显示类似格式错误的问题

因为OpenFOAM里chemisrtyReader和chemkinToFoam只会读transportProperties这个文件里的数据，类似下面这种。用在sutherland公式里。

".*" { transport { As 1.512e-06; Ts 120.; } } "H2" { transport { As 6.362e-07; Ts 72.; } } "CO2" { transport { As 1.572e-06; Ts 240.; } }

@ann 在 chemkin 中说：

我在网上找到的trans文件都长类似的样子

chemkin软件的三个文件之一就是这个六列分子特性数据的trans.dat，这才是正确的。只是OpenFOAM因为用sutherland公式就不需要这个文件了。sutherland公式需要的是上面的As和Ts。
使用OpenFOAM自带的chemkinToFoam命令：

chemkinToFoam chem.inp therm.dat transportProperties reactions.dat thermo.dat

大家好，我用 icoUncoupledKinematicParcelFoam求解器模拟静电除尘器中颗粒荷电运动轨迹，颗粒考虑曳力、重力和电场力，计算电场力需要电荷量q（通过方程计算，是动态变化的值）和场强E(矢量场，给定）两个参数，我在submodels/Kinematic/particleforces路径下添加了一个电场力子模型，怎么能做到q的实时读取和电场力的实时更新呢，求大神指点。
0.png

@wangfei9088 邮箱联系您了！期待您的回复

第一列为能量，第二列电离系数为关于能量的函数，能量由能量方程求解，第二列需要根据场中的能量数值插值计算当地的电离系数。
e6d35e3b-505f-4342-ad41-cc3aa7905308-image.png
22ee645c-9f2d-4540-84ef-84a11e942e6f-image.png 替代文字

@李东岳 东岳老师您好，我使用通过结合 waves2foam 和 overInterDyMFoam 的用overWaveDyMFoam求解器计算波浪与旋转模型相互作用的时候,重叠部分区域的流体不随时间改变，重叠的部分已经设置成overset边界条件了，想请教老师是否知道这个是什么原因导致的。
113.png
背景部分如图
111.png
重叠部分如图
112.png
boundary文件设置如下：
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: v1912 |
| \ / A nd | Website: www.openfoam.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

8
(
oversetPatch
{
type overset;
inGroups 1(overset);
nFaces 56;
startFace 222062;
}
inlet
{
type patch;
nFaces 98;
startFace 222118;
}
bottom
{
type wall;
inGroups 1(wall);
nFaces 1132;
startFace 222216;
}
outlet
{
type patch;
nFaces 98;
startFace 223348;
}
atmosphere
{
type patch;
nFaces 1132;
startFace 223446;
}
frontBack
{
type empty;
inGroups 1(empty);
nFaces 221480;
startFace 224578;
}
sides
{
type overset;
inGroups 1(overset);
nFaces 2000;
startFace 446058;
}
cylinder
{
type wall;
inGroups 1(wall);
nFaces 80;
startFace 448058;
}
)

// ************************************************************************* //