CFD中文网

同问,我也正在做类似的东西,但是算到最后连续性方程都是收敛的,但湍流方程总不收敛

@guanghuilineu 那看起来发散是个很正常的结果，如果相变区速度不为0时，也变速也变密度。

如果只变速，和周围空间交换能量，即动能转换为其他形式的能量，没问题。
如果速度和密度一起变，动能是依赖于质量的。这时候，质量也变了，不守恒了，就崩了。

消失的那部分质量，带有的动能如何传递出去。

压力就是传递能量的桥梁。在传递的时候一头突变了，就断了。可不可以先使相变区域停止流动再变密度。如果这样，不会出现动能突变的状况。或者在变密度操作后立刻更新一下速度；变前计算动能；变后用动能，根据变后的质量，计算新的速度。

或者最粗鲁的方法，用极小的时间步，这样突变也可以稳住。

有些东西就挺好意思。CFD原版书这个玩意，绝对属于精神食粮。但是偏偏就是国内买不到。就比如这本Computational Methods for fluid dynamics，国内就是运不进来。他们不给国内发货啊。

没有精神食粮怎么搞。。我就帮大家搞吧。http://dyfluid.com/profs.html 可以点击一下看看我对这个书的评价。我个人是强推这一本书。还有其他书，大家也可以看看。

这个价格就是这样算：国外amazon的原价，每一本加上10欧的运费（我找人带回）。CFD原版书一般一本在100欧左右（一本书800元左右），贵的150/200欧（1600元一本），便宜的70/80欧（600元一本）。

货期就不一定了。一个月之内没问题。但是我相信大家这个也不着急。一般看书没有着急的。并且还有电子版。一般买CFD原版书，主要是收藏用。因为不一定啥时候有人回国。如果发dhl，那还很贵。如果不差钱的可以告诉我发DHL。

也可以开发票哦，比如服务器租赁费测试费培训费啥的。感兴趣的回帖联系吧。如果是有过合作的客户。直接邮件联系也行。

目前我们可以这样操作。反正就是在这里报名吧。想要的把联系方式留一下我让助理联系你。

纯给大家折腾玩，生命不息，折腾不止。

同问，目前在用PBM处理气液两相流，用discrete方法求解PBM，也想试试qmom方法，上面几个矩代表什么意义？要模拟气液流中的气泡，该怎样设置粒径，数量还是。。。？

@李东岳 在 求流体凝固为固体的算法 中说：

@北斗 最近在研究多孔介质。目前一种非常简单的凝固模型，就是利用多孔介质来处理。固体区域给定一个非常小的渗透率。

血小板聚集形成血栓

这种看起来像是拉格朗日粒子的黏连

请问李博士，您的意思是通过多孔介质来实现凝固吗，如何实现，这方面资料可以发出来吗

在验证质量守恒，计算域内的质量可以求出，想请教一下，出口液相的质量应该怎么计算呢

@东岳 谢谢教授

请教大家一个问题，笛卡尔网格，已知网格8个顶点的坐标，以及网格中心的物理量的值，但网格排序是杂乱无章的，有什么好方法进行后处理呀？

按照李东岳老师的case修改即可完成本项工作。:huahua:

@hurricane007 老铁，发个联系方式

@cccrrryyy 有可能, 但是确定流体域中大部分为湍流，只有细缝附近不是很确定，这种在Fluent里面如果一部分为层流一部分为湍流，感觉就不怎么好处理了。

感觉DDPM比较适合你的例子，就是计算慢了些，或者用双欧拉

本人按在实际出口下方5cm处建立一个出口边界，实际出口则设为内部边界，最后调整参数，得到了与文献相似的结果。感谢各位的帮助！！

该书11.8节（中文版305页，英文版411页）中的例2，推导源项的表达式时，高阶格式减去一阶迎风格式，一阶格式的值取上上游的值φu，为什么这里不取φc呢？先谢谢了。
7a735b6e-01c9-4088-971a-359144da8189-image.png

文章题目为《The inviscid impingement of a jet with arbitrary velocity profile》，主要是分析了法向二维冲击射流流场的解析解，但是个人比较疑惑的是为什么控制方程只包括了压力的泊松方程（如下）而没有包括二维的欧拉方程？
d64e75d4-5884-42a4-b208-20c90325cc29-image.png

全文见附件，希望有朋友可以解答。The inviscid impingement of a jet with arbitrary velocity profile.pdf

自问自答一下
C-AMG和SA-AMG的区别在于粗化方案不同以及插值矩阵的构造是不一样的，C-AMG基于细网格生成粗网格，可以理解为更精细，因此计算收敛性会更好；SA-AMG根据聚类的思想来的，所以在setup阶段会大大缩减时间，但是收敛性不如C-AMG；而两者都比较鲁棒性，因此并不是只针对某种方法有效。
而后一个问题目前调研到的是，AMG在商软中均可以独自使用或者作为预条件使用；而Fluent、COMSOL、STAR-CCM+这三个软件的官网都有对应的介绍，基本上作为预条件都是与CG/BiCGSTAB联用，其他的还有GMRES。
找了好久也算是把资料找到了

@blgy 非常的感谢，你这个是在哪看到的啊？你这个图片也看不出来用什么软件，用什么方法划分的啊，仅仅是欣赏一下:ok:

在破碎模型中，破碎的概率和子气泡分布一般涉及到积分式，比如说经典的Luo模型270626e3-d11f-4d47-a461-cdc3cbde1bb6-image.png 这里边openfoam官方直接用大佬提出来的化积分为多项式公式替代求解：e4c4c71f-d206-4464-bb15-5ee960dcdba8-image.png 有学者就提出来了这个exp内的判据不精准，使用表面能和韦伯数多重盘踞，然后这个积分式就这样了：
6f926600-85cf-4d78-a936-ef166ab439ab-image.png 现在这样的话积分还怎么求解呢，因为min和max函数的存在原来的积分式公式还能使用么。还是需要分段求积。原论文是植入到fluent中大佬没写怎么处理的.

加上来了

@冰冰玩 好的，谢谢！我之后尝试一下。