主任招博士后不，我给您发个简历

压力没有一个绝对值的话，必须需要给一个参考值。但是你是自己写的代码，不好debug

aurora，可以在word里用latex

请问有谁对 bodyPoint（refAttachmentPt_）这个函数了解吗？
这是函数的具体实现：

inline Foam::point Foam::RBD::restraint::bodyPoint ( const point& p ) const { return (model_.X0(bodyID_).inv() && spatialVector(Zero, p)).l(); // 双内积 }

哦，明白了。

@东岳 李老师，我没推出来

确实和概率论有关，谢谢，记录一下

\begin{equation}
\left( \begin{matrix}
k\\
x
\end{matrix}\right )=\frac{k!}{(k-x)!x!}
\end{equation}

@promisegwj 啊不对 是用拉梅系数···

如果不考虑时间项
\begin{equation}
\int^{f+1}_ {f} \frac{\partial}{\partial x}\left(\alpha \frac{\partial \phi}{\partial x}\right)\mathrm{d}x=\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f}
\end{equation}
把后面两项$\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_{f}$把所有的面连起来，看起来是守恒的，剩下的只有边界条件

@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，找了半天没发现错误，更改之后搞定了，谢谢东岳老师

我没太跟随数值传热学中压力修正方程的推导，目前还不能给出推导: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