@wg0632
Gauss LUST grad(U) 当中的 grad(U) 对应一个梯度计算格式，会从 gradSchemes 部分取查找。如果 gradSchemes 部分定义了 grad(U) ，则使用该部分的 grad(U) 对应的格式。如果gradSchemes 部分没有定义grad(U) ，则会使用 gradSchemes 部分的 default 对应的格式。如果 default 也没有则报错。

所以，Gauss LUST grad(U) 和Gauss LUST unlimitedGrad(U) 的区别取决于 gradSchemes 里面如何定义的 grad(U) 和 unlimitedGrad(U)。

参考multiphase\overInterDyMFoam\rigidBodyHull

老哥学明白了吗？

https://cfd-china.com/topic/6644

你看看这个，我觉得可以用symmetry边界

@hachikoi 注意到了，似乎在螺旋桨与船连接处存在了一个disk薄片，

有大佬知道如何在V2206版本编译moody系泊插件吗？

@tyhero123 好的谢谢

如果每个网格在同个时刻的平均速度是一样的，更简单的方式是按照C程序写法直接读取这条曲线的离散数据，插值出当前时刻的平均速度，并加到decayingTurbulenceInflowGenerator中的U[faceI]

@wjxhit 这不是网格的问题，是那个位置的粒子正好被注射在了两个cores的中间应该。不知道你具体用的什么注射方式，但只有移除那个粒子才能接着run。

@Jasper-0 你现在能堆积，是因为把youngsModulus系数增大从而增大了碰撞力。这会引起无法堆积的，你可以看你的结果，应该会一直蹦，因为排斥力过大。这个系数因材料而定，不能随便给。比如沙子大概e+06好像。碰撞力的大小需要在计算时调节需要的参数，不能改变粒子的物理属性。

是向后差分，但是我没有用到ddt，是自己写的函数推进时间步

@李东岳 好的，谢谢老师！

我的chemistryProperties是这样设置的：

chemistryType { chemistrySolver ode; chemistryThermo psi; } chemistry on; initialChemicalTimeStep 1e-8; EulerImplicitCoeffs { cTauChem 1; equilibriumRateLimiter off; } odeCoeffs { solver seulex; absTol 1e-12; relTol 1e-1; } // ************************************************************************* //

twoSymm(fvc::grad(U))正解

http://dyfluid.com/tensors.html

@李东岳 感谢！看样子这方面的相关研究确实比较少

好的，谢谢大佬的讲解，我再看看相应求解器下的能量方程。感谢

液池的侧面是给的固壁，上面是开放边界

@李东岳 在 浮力流的初始化问题 中说：

test22.tar.xz

直接blockMesh，然后buoyantFoam跑就可以，我跑到了20多秒没啥问题，流场看起来也正确。

谢谢，我试试看，我一直是发散的。

paraview可以直接对分块的网格后处理

6ad7cb0e-fb0f-4d5d-9c22-b54b6109e592-1700212552480.png
难道还是发散了吗