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

没太看明白

@Qing 感谢老师！您说的特别清楚，这对我特别有帮助。再次谢谢您！

大佬们，还在等一个回复呢

也不知算不算高粘度，1-几百mpas都有，感觉速度对的还行，液体表面泡沫层是真的没有

我勒个去，好老的帖子，一晃8年过去了

拉格朗日一般要从粒子碰撞的角度来考虑破碎聚丙啥的

属实没看明白。有一些教材，sci，写的就很晦涩。

很有帮助。

@WeakForm 这是我微信liuzongku

@李东岳 李老师，请问如果求解非守恒型的动量方程，那么对于其中一相为0的情况，是不是就直接求解就好了，不用再做别的处理了，如让相为0的相的速度等于按照模型去求出的速度。

我自己在用TDMA算法，用Matlab

伴随优化我一直想渗入 就是一直没时间详细研究

感觉这个可以发一个mature！

@thegame 感谢老铁，将在下一版本更新

多谢两位老师

@李东岳 感谢东岳老师的回复！

有个方向，叫混沌之类的。详细可以了解下蝴蝶效应和非线性的一些内容。

请问各位前辈，有没有界面蒸发源项的模型啊，找了好多都是壁面成核的，因为要做的工况是在气液界面处产生相变，而不是壁面加热导致相变的，看了一些文章好像找到没有one-group下界面蒸发的模型:136:

就是下面这个
$$
\frac{\p T_{ij}}{\p x_i\p x_j}
$$
其实$T_{ij}=u_iu_j$，所以演变为
$$
\frac{\p u_iu_j}{\p x_i\p x_j}
$$
这个公式是可以写成$\nabla\cdot(\nabla\cdot(\bfU\bfU))$这种形式，$\frac{\p T_{ij}}{\p x_i\p x_j}$好像写不出来 :134:

大佬们看看？

@pengdi 前辈还有个小问题想问问您
您当时求解RP方程后（计算出当前气泡的半径）
6c4261f0-8249-4c72-a1b7-f96df479c043-image.png
是不是需要重新求解平衡状态下的R0和P0？还是直接使用第一次计算得到R0和P0（求解方程如下）
4c110999-f16d-4fa2-ba11-53d27f1e33a9-image.png

因为方程(1.12)中的PB需要用到R0，所以我比较疑惑，想问问您
e2f9fc3e-4af3-40d0-8a8b-0746f3bd41d8-image.png

fluent的udf我不熟悉..

对于URANS，多了个瞬态项，那每个时间步计算的时均速度（时均动量）都是从第一个时间步开始累加后平均后获得么？

对于$\int_t^{t+\Delta t}$的$\Delta t$，从理论上URANS这个$\Delta t$是任意的，如果不考虑附加的湍流封闭模型如ke之类，$\Delta t$足够的小，就变成了LES，$\Delta t$足够的大，就变成了RANS。然而$\Delta t$并不需要显性的给定。所以对URANS取决定性的还是$\mu_t$

@ZhangMengYi 没有，好像只见过Deen那面搞过，不知道他们能不能给你

均匀各向同性湍流，其雷诺应力矩阵非对角项等于0 竖向非均匀各向异性湍流，比如槽道流（水平向均匀），只有$ \left \langle u_{1} u_{2} \right \rangle $有值，其他非对角项均为0。Pope的书中也给出了解释。 现在问题是，均匀各向异性湍流，其雷诺应力矩阵非对角项是否为0？如果参照Pope书中类似推导，在空间3个方向都是均匀的，理论上$ \left \langle u_{1} u_{2} \right \rangle $，$ \left \langle u_{1} u_{3} \right \rangle $，$ \left \langle u_{2} u_{3} \right \rangle $都为0吗？

《Pope, S.B., Turbulent Flows. Vol. 12. 2000: Cambridge University Press.》89页
e0bada21-7ac0-4ae4-acc3-3319a012dfc7-image.png

理论上$\alpha$要么是0，要么是1，因此这个要么就是液相速度，要么就是气象速度，不过写成平均速度也一样，因为$\alpha$要么是0，要么是1

@李东岳 :146: :146: :146: ，谢谢李老师，醍醐灌顶

@WeakForm 是的。不管是多大的相分数，是多大的速度。结果都一样

我没太看过相关的讨论。目前这种只有周期性边界才不会发生这种问题。你说的这个问题再往深了我就没关注过了，遗憾，老铁整的挺深啊:146:

@李东岳 哇，感谢李老师:140:

@hurricane007 我去，忘了。。。这个我当时打算自建一个免费的计算中心。后来发现一直需要给客户装机，电力都放在那面了，这个自己用的一直没搞起来。

我距离百亿市场只是因为电路容量不够 :136:

@李东岳 角动量与C成正比，它在时间上的导数就是受到的力矩，这里面写了粘性力矩; 压力应该没有力矩，因为是势能。另外套上了对流项，建立的方程

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

@星星星星晴 OpenFOAM中的DPM粒子追踪功能，在运动网格(ALE)上的支持并行计算吗？有没有什么限制？

最近看到一种说法，level set方程是简化版的相场方程，结合这三种方法在方程形式上非常相似，不知道从计算精度、计算量等方面来讲，这三个方法的优缺点有没有什么对比性结论？

@李东岳 谢谢李老师，我再多看看文献

@疏影横斜水清浅 自问自答吧，是组合公式
1f3e206f-d460-45ee-9db8-ffeb3cf340df-image.png
C（k，i）

@aiweimo 也行，有链接不

@李东岳 那理解了，谢谢李老师

我对碰撞那面的了解比较浅。据我所知，FVM也可以用来算固体力学，起码小变形可以。OpenFOAM这面算小变形很简单，原理就是1）用FVM离散固体位移方程，2）网格变形，3）下一步推进。

我不太理解的是，按照这个思路，大变形的情况下，网格变的非常非常差，像是左边那种可乐罐子的网格，按照FVM的角度来看，基本是没法算的。

很好奇lsdyna，或者类似的软件是怎么实现的？

替代文字

截断误差首项阶次一样

我推了一下，CD的更小。不一样啊。

捕获.PNG

@李东岳 收到，谢谢李老师。

是的。

openfoam这面PU的求解框架大体一样。

@李东岳 明白了，谢谢李老师耐心细致的解答。

@李东岳 在 DES模型 中说：

这个找找他们做湍流的原始文献应该会有详细的介绍，就是这些DES模型为什么做这些更改。

好的 我再深入学习下，谢谢李老师

“对流”-“扩散”型、偏导数型（对时间t的偏导数加对xyz的偏导数）

没太理解这个描述。

此外，湍流输运方程（如k方程、omega方程）中的dissipation term中也不包含偏导数

在源项里面

@李东岳 感谢李老师，解决了我的心腹大患:xiabanle: :140:

@闹市闲民 再给推荐一本粘性流体力学教材，个人亲测食用最佳:chitang: 7d7a70c8-f0cf-4676-b460-7866a9178ba8-image.png

很多湍流模型的假定就是各向同性，这样模化起来、封闭起来比较简单。各向异性那种，就得上雷诺应力模型了，明显要复杂

@cccrrryyy 谢谢老师