@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
难道还是发散了吗

在paraview中，采用plotoninterectioncurves方法分析边界面的数据时，
73802fe2-b805-4040-a061-4e64402d5cba-a73bd1166d3518015a17532210b349d.png 、
这两个Mesh框是否勾选，输出的结果也不同。应该如何选择呢？
求各位老师指教

4472302b-0218-4b15-9bf2-4e4a46f06c78-image.png

请问你解决了嘛？我也遇到了同样的问题

放了固定结构物，只能跑0.几秒，没有报错，库朗数小于1，集群提交作业里面显示还在运行，但是dat文件不再更新。
300b9d82cdd8eb2cce0d25e75e1747a.png

你说的没毛病，老铁

@李东岳 在 LES计算中入口 湍流度/脉动速度 如何设置和控制？ 中说：

后续的有关LES入口的可以参考这个帖子 LES定义入口速度的问题（DSRFG方法）

用的动态拉格朗日模型

@aiweimo 大佬你当时为什么用这个湍流模型？有特殊的原因么

李老师，没有特殊原因，当时用的是第三方的SOWFA代码包，刚好用了这个模型，没有对亚格子模型进行研究过