@Do1975 感谢大佬，我去试试

@李东岳 CFDEM采用的是IB求解器，CFD计算采用了pisoFoam求解器

@bestucan 好的，多谢老师解惑，我再好好找找原因

@李东岳 李老师 动图中iter是指迭代步数。目前设置的是计算速度场计算200步完成后，再继续伴随速度场，这边为了减少gif的大小，我只生成了前50步的结果。
在10-15步的时候，此时计算的Ua是正常的，但是随着迭代步增加，计算就不收敛了。一直没搞明白为什么。

这是其他密度场下计算得到的Ua场：
2da7895d-d452-4997-bef2-ce5c1f02afb9-image.png

现在的问题是，有些密度场下可以计算，有些就异常了。（密度场是通过0时刻读入的，作为一个初始值）

volScalarField xh ( IOobject ( "xh", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), mesh, scalar(voluse), zeroGradientFvPatchScalarField::typeName ); xh.read();

@AppleKiller 这个不大清楚。我版本比较低似乎没有这个。system应该是整体系统吧，下面patch就是你自己设置的边界位置。你可以算算， 你设置的这个时间段颗粒等不等于这个system added，估摸着算算

@李东岳 是液面几乎不变时当单相流算，自由液面当刚盖

如果你壁面是水平的，$\bfn_f$指向上，那壁面切应力就是水平的。没问题啊。

两位前辈我现在在多轴组合运动上存在一些问题，想请教一下你们是怎么设置这种组合运动的吗？如果方便能给我一个类似案例设置文件吗？我的问题是一个球体在一个圆柱流域中公转加自转的组合运动。我目前尝试了设置，但是计算出现了发散，我转速很小的情况下库朗数居然会增加到几万。@卡洛 @冷空气LQW

/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2406 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dynamicFvMesh dynamicMultiMotionSolverFvMesh; // Specification of rotation around centre rotation1 { solidBodyMotionFunction rotatingMotion; rotatingMotionCoeffs { origin (0 0 0); axis (0 1 0); omega 5; } } dynamicMultiMotionSolverFvMeshCoeffs { domain { solver solidBody; cellZone domain; solidBodyCoeffs { $rotation1; } } ball { cellZone ball; solver solidBody; solidBodyCoeffs { solidBodyMotionFunction multiMotion; multiMotionCoeffs { rotation_1 { $rotation1; } rotation_2 { solidBodyMotionFunction rotatingMotion; rotatingMotionCoeffs { origin (0 0 0.018); axis (0 1 0); omega -50; } } } } } } // ************************************************************************* //

是想用cantera算一维火焰做FGM表，用在openfoam计算里ToT

MomentumParcelI.H

template<class ParcelType> inline Foam::scalar Foam::MomentumParcel<ParcelType>::Re ( const scalar rhoc, const vector& U, const vector& Uc, const scalar d, const scalar muc ) { return rhoc*mag(U - Uc)*d/max(muc, rootVSmall); }

请问OpenFOAM有关于波浪与柔性结构相互作用的求解器吗？

请教一下各位大佬，我在做双自由度涡激振动的时候，设置了如下的dynamicMeshDict文件，openFoam动网格程序是如何识别水平方向和竖直方向的弹簧的呢？还是说我这样的配置程序根本无法识别，但我这样的设置程序能正常计算，而且两个方向都有运动，不过模拟结果确实差蛮多。我需要怎么去调整呢？
d55f0aa65c8502fd9d43f11f8aea289.png

好的谢谢老师

不好讲故事是比较困境了，哎..

算例是twoSpheresGlowinskiMPI

似乎解决了这个问题，把文本中的单引号改为双引号后，现在切换版本运行目前没啥问题:135:

@benqing 在 SRFSimpleFoam里边界条件的设置？ 中说：

你好，我最近准备使用SRF求解器，但是应用过程中一直存在问题，请问明白前述的问题了吗？
或者，有没有关于SRF求解器的资料可以分享一下？谢谢！

您好，我最近也困于SRF求解器，请问您现在怎么样了，有找到相关资料可以分享一下吗？谢谢

这个是ESI-OpenCFD版本 git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2312.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2306.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2206.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2106.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2012.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2006.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1912.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1906.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1712.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1706.git

第三方包地址：

git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2212.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2206.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2106.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2012.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2006.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v1906.git

@李东岳 好的，谢谢老师

@crispy 大佬，请问你验证的雷诺数是多少呢？现在有比较好的结果了吗？