@lyc 我也是菜鸡，这玩意我师兄去年也想过，一直没成功，今年我做课题时候要添加自己颗粒属性方程可能接触一点。但是如何实现我也在研究中，你可以看看我之前和其他大佬的一些疑问，他们指出过在哪些文件里面进行属性添加修改，不过都不全，导致我目前现在想要修改内容也还差一点

如果我10月份有时间的话，或许会写个解析，更新到 http://dyfluid.com/cfd.html
就是我这个时间比较紧

@李东岳 在 关于LES模拟中采用大时间步长的问题 中说：

你可以放大库朗数，或者你可以用URANS更快。这些都没问题。真正分析的时候库朗数放小就行。并且库朗数一般是局部库朗数的最大值。如果1000个网格里面有1个网格特别小。会出现999个网格的库朗数都在1以下，但这1个网格在10左右。都不是问题。

好的，谢谢李老师！我明白了

@李东岳 感谢回复。我是想先写个简单的特定算例学习一下算法，格式边界条件啥的都直接写死，然后再慢慢做的通用一点，能搞出来肯定开源。

我今天看一篇论文，一个简单的方形腔，通过添加金属颗粒来改变腔内相变材料的热性质。文中一个研究点是改变金属颗粒的体积分数，按理说相同的纳米颗粒体积分数下混合流体的热性质都是一样的，根据Gr数计算公式，Gr数不应该变化呀，但是相同的颗粒体积分数下文中又设置了10000与100000两个量级的Gr数，不知道怎么弄的，求大佬解惑。
对了，我找遍那篇文章都没发现方形腔的物理尺寸，不知道这个与Gr的控制有没有关系

居然还发公众号了

三维的类似情况，我最后一个打错字了。

@李东岳 明白了，谢谢李老师。排版排了好久，清楚一点方便讨论:papa:

As an application of this geometric formulation we prove that the solution map of the incompressible Euler equation, associating intial data in Hs(Rn) to the corresponding solution at time t > 0, is nowhere locally uniformly continuous and nowhere differentiable. by Hasan Inci (2013)

teoria.png
如图所示，这是非常经典的Le Méhauté (1976)波浪理论适用图，但我有点迷糊，主要有两方面疑问：

图中红色部分看起来是Stokes三阶波/四阶波与五阶波的交接处，那么我应该选择用三阶波/四阶波还是五阶波呢？同样地，图中浅蓝色和五阶波也有交集，那么在既是浅蓝色又是五阶波的区域应该用cnoidal还是五阶波呢？ 我还注意到图中五阶波下面还有一个流函数理论，但我看olaFlow的用户手册，用流函数的话需要输入一大堆参数（如下图所示），这些参数应该怎么选取呢？本菜鸡在此请教一下各位:135:
2024-07-08 21-55-00 的屏幕截图.png
2024-07-08 22-01-59 的屏幕截图.png

推荐一本CFD-DEM基础书籍，感觉其中的unresolved CFD-DEM的理论公式说明比较清晰。
Norouzi 等 - 2016 - Coupled CFD‐DEM Modeling: Formulation, Implementation and Application to Multiphase Flows

@李东岳 好的，谢谢李老师 ：-）

@youhaoyu 1719302824186.jpg 看了本比较通用的书，里面把我第一张图里（1.15）直接写在方程里面的（我感觉这样写比较合适）。那么这个R0和Pg0我还是不理解什么意思，有没有大佬解释解释

opm_set_num_threads应该是设置openmp并行线程数量的，具体错误原因还需要看详细代码

@AppleKiller 感谢，我再看看，以前也接触过这些资料，感觉用欧拉法求解器为基础修改的拉格朗日求解器不大好，没怎么深入去学习。我再看看这些思路，感谢

@李东岳 东岳老师您好，我最近发现了一个现象可能造成了我上述提出的问题哈哈哈哈哈，我最开始使用的是DPMFoam默认的的Ergunwenyuforce曳力，我后来换成了sphereforce曳力就能稳定了，而且对于颗粒运动和流场有挺大影响，看来这和我的设置也有一些关系。（最开始图方便就直接使用的默认的模板了，看来还是要根据实际情况来改变）

哪位老师知道目前哪一款CFD软件能够精确模拟这种现象？感激不尽！

@RQLI_CFD 感谢，我研究研究

一般泰勒格林涡的文章都会指定计算域长度与涡的大小啊。比如这个：https://www.jianguoyun.com/p/DblYnzoQ9s3ZBhil274FIAA 方程6.2

我感觉有必要擦除流体与固体的试算结果，再进入下一个迭代步。固体好说，只有一个D，我令D=D.oldTime();就能把D的试算结果消去，但流体就很麻烦，因为有很多个变量，我试了半天也不知道该怎么把所有变量的值都擦除掉，恢复迭代前的值

没太看明白