可不可以结合codestream实现旋转呢？

@Samuel-Tu 请问你解决了这个问题了吗？我是在controlDict中添加code进行计算得到的一些自定义的变量，但是结果在写入到文件中时，除了时间步数之外，其他的参数就一直是6位数或更少，小数点后边的就更少了。也尝试改了李老师说的这个，但是好像没用。

@zhoushan 简单的几何可以用blockMesh直接生成网格，复杂的就要绘制stl文件，用blockMesh搭配snappyHexMesh生成

fluid@zhoushan:~/$WM_PROJECT_USER_DIR/FluidSolidInteraction/src$ ./Allwmake

wmake libso fluidSolidInteraction
./Allwmake: 4: wmake: not found wmake solvers/fsiFoam
./Allwmake: 6: wmake: not found wmake solvers/ampFsiFoam
./Allwmake: 7: wmake: not found wmake solvers/weakFsiFoam
./Allwmake: 8: wmake: not found wmake solvers/fluidFoam
./Allwmake: 9: wmake: not found wmake solvers/solidFoam
./Allwmake: 10: wmake: not found wmake solvers/thermalSolidFoam
./Allwmake: 11: wmake: not found wmake libso utilities/functionObjects/pointHistory
./Allwmake: 13: wmake: not found wmake libso utilities/functionObjects/patchAvgTractionHistory
./Allwmake: 14: wmake: not found wmake libso utilities/functionObjects/centrifugalBodyForce
./Allwmake: 15: wmake: not found
请问上述的not found怎么解决呀 在网上参考了各种方法 都无法成功

最近在做多相流的仿真验证，但是一直找不到LS方法的多相流案例，想请教以下各位大佬openfoam中有没有关于LS方法的多相流案例

如上所述

这个问题我也遇到过，比如对于同一个流场，同样的空化数，来流速度10m/s就不会产生负压，来流速度20m/s就有负压产生。增加SS空化模型的系数n的值或许能解决问题。

今天偶尔看到一个文章。好像意思是通过Optuna + OpenFOAM，可以实现最小的计算时间。感觉类似通过Optuna跑好几次，然后寻找出来最合适的松弛因子。不知道有没有哪个大佬感兴趣玩一玩。

https://ss1.xrea.com/penguinitis.g1.xrea.com/study/OpenFOAM/relax_factor_optimization.html

@wangfei9088 在 paraview 不显示坐标轴数值 中说：

@liujm 假设你用的版本是OpenFOAM-10，其他版本也一样。
路径：OpenFOAM-10/etc/config.sh/
打开paraview文件，将69行取消注释，注释掉72行。
./makeParaView重新编译就是5.0.1版本了。

#export ParaView_VERSION=5.0.1 //69行 #export ParaView_VERSION=5.4.0 #export ParaView_VERSION=5.5.0 export ParaView_VERSION=5.6.3 //72行

感谢！

@LGY 请问一下你这边能够在wave2Foam中生成非定常风了吗

@gooseEast 在 rhoCentralFoam发散，出现Maximum number of iterations exceeded 中说：

在thermoI.H文件第46行标量f是什么

OpenFOAM代码中，这个f是Cp，定压比热容，单位是[J/(kg K)]。这个T函数是通过定压比热容的值迭代求出温度。limit，F，dFdt分别对应的是hConstThermoI.H文件里limit，Cp，dCpdt三个函数。

@李东岳 老师，如果不使用refineMesh的话，结果会好很多，而且和论文中的Cd, Cl指数可以基本对上，但是在overset的边缘处还是会有一点问题。
648bb924-7f1b-4e8c-a39d-48cc7c40906e-image.png

我的oversetInterpolation的method使用的是inverseDistance.

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

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