@vbcwl @hhh 我觉得$\tau_w$是负数是正常的，负号代表方向

是这样的，非常感谢您

公式来一个看看，老铁

@李东岳 谢谢老师:140:

如何在OpenFOAM里面实现对颗粒各种力的添加，具体的步骤应该怎么做

可以，没问题，自带的算例就有一个类似的。wingMotion2D_pimpleFoam

没玩过liggghts，知识盲区了

是的，数量2，内核16，可以分32
不过虚拟机性能损失挺大，核数多也没啥大用

@李东岳 这个算例我设置的45s，单核跑到4.2s左右报错

@李东岳 我的模型两侧还有静止区域，只是这个简单模型还没加上，所以SRF是不适用的。

@qingdong_wang 还有这个pimple算法这，按我理解不应该是速度预测，压力修正，速度修正，迭代不足回到压力修正这种？为什么您这个是速度直接求出来了，然后直接就压力值了，迭代不成功回压力修正

@sxz0823 似乎没有自带的，只能自己编写，这难度不小

大佬做过粒径或者浓度不同高度注射吗

@李东岳 在 不使用湍流模型，而是直接求解器中实现湍流计算 中说：

有一些代码不知道有啥作用，那就留着，留着总比删了更保险。

好像是的，对结果不影响的都会跟着作者留下来了。:papa:

@李东岳 在 不使用湍流模型，而是直接求解器中实现湍流计算 中说：

这个在extend那面很常见。

extend我只是听过，但没有用过。现在我常使用的版本一个是Org 6.0的，另外一个是ESI v2312的。

@chenboyao 在 寻求帮忙看看采用重叠网格对平台做自由衰减后自由衰减曲线不符合衰减曲线的原因 中说：

我采用重叠网格对浮式平台做自由衰减，结果自由衰减的数据曲线到最后不符合自由衰减的趋势，而是还向下偏移了，不知道是什么情况造成的，希望有人帮我看看一下，为什么会出现这种情况。
Screenshot from 2024-06-17 09-20-41.png

您好，请问您的自由衰减实验做成功了嘛？我现在也在做浮式平台的自由衰减实验，也遇到过类似的问题，想很您交流一下:chouchou:

这个是当我把水的物性从常数改为多项式后遇到的问题
Courant Number mean: 0.134868655 max: 0.237479859
deltaT = 3.4820232e-06
Time = 0.00378807

PIMPLE: iteration 1
MULES: Solving for alpha.gas
MULES: Solving for alpha.liquid
gas fraction, min, max = 0.00155200072 1e-15 0.0672465558
liquid fraction, min, max = 0.998447999 0.932753444 1
Phase-sum volume fraction, min, max = 1 0.999994574 1.00026614
MULES: Solving for alpha.gas
MULES: Solving for alpha.liquid
gas fraction, min, max = 0.00155951826 1e-15 0.0674835707
liquid fraction, min, max = 0.998440482 0.932516429 1
Phase-sum volume fraction, min, max = 1 0.999994794 1.00025528
MULES: Solving for alpha.gas
MULES: Solving for alpha.liquid
gas fraction, min, max = 0.00156703579 1e-15 0.0677203557
liquid fraction, min, max = 0.998432964 0.932279644 1
Phase-sum volume fraction, min, max = 1 0.999995014 1.00025493
Constructing momentum equations
Tf.gasAndLiquid: min = 350.131001, mean = 456.035361, max = 549.778465
iDmdt.gasAndLiquid: min = -0.959612341, mean = 11.8287784, max = 447.419943, integral = 0.000182163188
wDmdt.gasAndLiquid: min = 0, mean = 138.094307, max = 5727.35791, integral = 0.00212665234
smoothSolver: Solving for e.gas, Initial residual = 0.104953394, Final residual = 1.40390231e-09, No Iterations 2
smoothSolver: Solving for e.liquid, Initial residual = 0.282448508, Final residual = 6.30884149e-18, No Iterations 2
[stack trace]

#1 Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#4 void Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#6 Foam::heatTransferModels::RanzMarshall::K(double) const in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#7 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::BlendedInterfacialModelFoam::heatTransferModel::evaluate<double, Foam::fvPatchField, Foam::volMesh, double>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > (Foam::heatTransferModel::*)(double) const, Foam::word const&, Foam::dimensionSet const&, bool, double) const in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#8 Foam::BlendedInterfacialModelFoam::heatTransferModel::K(double) const in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#9 Foam::ThermalPhaseChangePhaseSystem<Foam::PhaseTransferPhaseSystem<Foam::TwoResistanceHeatTransferPhaseSystem<Foam::MomentumTransferPhaseSystemFoam::multiphaseSystem > > >::correctInterfaceThermo() in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#10 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/reactingMultiphaseEulerFoam
#11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#12 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/reactingMultiphaseEulerFoam

@coolhhh 存在的，每个保存的时刻里面都有Qdot

李老师，今天问题解决了，原算例的初始场是用mapFields从另一个静态网格的中间时刻得到的，如果不用mapFields设置初始场似乎就一切正常了，不过我确实想不明白为什么mapFields会产生这样的影响……

@李东岳 在 关于进行瞬态计算矩阵迭代过程的数据操作 中说：

是要改很多私有成员变量啥的，来回嵌套。我想想有没有什么好办法。

谢谢李老师关注，确实我本来也以为改一改矩阵求解器就行，结果发现OF里嵌套继承的类太多了，越改越麻，所以才来提问，希望有更好的方法。

kOmega模型是高雷诺数湍流模型。如果做了壁面解析，需要保证在1左右，如果使用壁面函数，可以在30以上。