vor.png

@haaDoo 物性参数都可以变，温度不一样物性参数就不一样，或者换不同的流体介质

@东岳 对于$\nabla{\mathbf{U}}$的矩阵形式应该是什么这个问题，从张量的角度来看我站第一种写法，更规范。:huahua:
但是从CFD角度来看，公式这东西本就是规则，制定了规则就是交流用的， 可读性和易理解性还是更为重要的。看个人理解，我觉得并无好坏。:mianmo:

@流动魅力 我比较好奇你的相间作用力的直径怎么给的

想开个贴讨论下多相流界面追踪类方法目前的前沿水平，例如能算到的最高参数，各类方法都有什么限制等等。目前多相流算法里，计算界面追踪的多相动力学，比如droplet collision, droplet splashing, droplet impact on solid wall，这几个典型案例，目前有用VOF，level-set，phase field，LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态，实验都做出来了，但感觉计算还没跟上。

调研过不下百篇的文献，基本上关键词相关的都翻了一遍，在文献中看到droplet collision大概最好的也就是Re=6500,We=1500，密度比DR=1000左右的，大致是这个级别，像droplet impact有篇08年的用TransAT算到了RE=11000，而最近16年左右的用levelset这些也算到Re=8000，We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟，而且网格也都比较密，计算资源消耗也都比较大，估计都得上超算，当然和实验也都对的上，基本都是到了prompt splattering这样的形态。

想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪？ 因为多相流在算到高RE和We的时候会数值发散，稳定性应该是主要限制，还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况，需要单独使用人工识别薄液层再稳定的方式，对不规则碰撞并不实用。

@东岳 我最近用reactingMultiphaseEulerFoam做算例，也是设置的时候某一相为0，湍流模型一直出现divide by zero，找了半天没发现错误，更改之后搞定了，谢谢东岳老师

[N-S方程前传：普罗米修斯] 剧情简介

21世纪末，人类的科技水平已高度发达，计算器求解N-S方程早已实现，不再是梦想。与此同时，许多科学家仍孜孜不倦追索着N-S方程的秘密与真相。通过对许多古老文明的考察与对比，科学家伊丽莎白·肖和查理·赫洛维发现，N-S方程可能是来自一个遥远星系的更底层的外星文明创造的。在dyfluid公司资助下，他们乘坐宇宙飞船普罗米修斯号前往那颗未知的星球。

经过对当地的考察，科学家的设想成功得到印证。可是他们贸然探查N-S起源的神秘行为，也将自己引入万劫不复的算法黑洞……

我没太跟随数值传热学中压力修正方程的推导，目前还不能给出推导:tishizi:

我自己的推导 http://dyfluid.com/docs/RhieChow.html

谢谢李老师

新人发帖，还不懂如何编辑帖子，请各位老师见谅。
捕获.JPG
该公式是安德森的空气动力学基础一书中对等熵条件下准一维流动推导出的速度面积关系式。根据该关系式，当管道流通面积A不变时，速度u是不变的，则根据质量守恒关系可以知道密度和压力都是不变的。但是我在观察真实条件下的抽气管道流动中，在管道流通面积不变的情况下，上游管道静压大于下游管道静压（下游安装有抽气机组）。请问这是否说明等熵条件下准一维流动速度面积关系式不适用于真实条件下的均匀管道流动？
流体力学小白，还希望有老师能指点一下。

$\frac{\mathrm{d} A}{A}=\left ( M^{2}-1 \right )\frac{\mathrm{d} u}{u}$

@东岳 原来是连续方程的引用啊，谢谢李老师啦

感谢分享！

用于求解不可压的pseudo-compressible method, 被认为有更好的收敛特性，见A. J. Chorin, JCP, 1967, pp. 12-26

http://www.cfd-china.com/topic/2744

我重新推出来一个带偏导数的压力方程..

ef44221d-062f-4315-b9c1-50e96688d88e-image.png

一个公式真是折腾死了。

我看了好几个人讲关于von Neumann稳定性分析。

12步入门Python这个女老师讲的。 MIT qiqi wang讲的PDE课程

还是有些地方很模糊。

@G8S7

究竟能不能用欧拉-欧拉模型来说明这个体系下固体颗粒的运动情况，

如果有实验图片的话更好判断。取决于固体颗粒粒径，相分数，粘度等。如果用你所说的VOF+欧拉-欧拉才能描述这一个过程，也是可以的。另一种方式是三相欧拉，在皮层熔体中对界面进行尖锐挤压处理。

另外，用什么模型取决你你关注的重点，如果关注内部的运动，皮层溶体不是重点用欧拉无妨，VOF也是有耗散。

At last, 几何很简单，是验证算法的利器，不用OpenFOAM可惜了。

@余正东 我都没听过 。。。

好的，谢谢

@肖艺 只是看过相关的影片.. :zoule:

我没说清楚。P23和P25的页码写在了下面。而前面所有的奇数页都在上面。没统一。