@东岳 不敢班门弄斧，我一直是你的粉丝:ok: :ok: :ok:

谁能帮忙填充一下没考虑到的内容么？

捕获.PNG

捕获.JPG

重写了一个，以后用这个更新

用的Adobe AI

@东岳 谢谢版主的解答。。这个问题暂时搁置了，以后有机会把UDF分享出来。

@FluidGao 嗯嗯，感谢。

前几天OKS的时候有个同学问欧拉方程是不是速度U和压力p按照声速传播。当时我给的答案是不是。在这里详细说一下这个问题。对于双曲系统
\begin{equation}
q_t+Aq_x=0
\end{equation}
其中$q$是求解的矢量，比如$[p,u]$；$_t$是对变量求时间导数，$_x$是对空间求导数，$A$是一个矩阵，但不是对角阵。双曲系统中的$A$是可以对角化的：
\begin{equation}
A=R \Lambda R^{-1}
\end{equation}
其中$R$是特征向量组建的矩阵，$\Lambda$是特征值组件的对角阵，如果定义$w=R^{-1}q$，最原始的方程就可以变成
\begin{equation}
w_t+\Lambda w_x=0
\end{equation}
这个方程式一系列解耦的双曲方程，可以所有的$w$都按照特征值的方向进行传播。因此对于欧拉方程，p和u并不是按照声速传播，按照声速传播的是$w$，$w$可通过$w=R^{-1}q$求出。

$w$也就是特征变量哦。如果感觉困惑，或许可以看看这个：http://homepages.see.leeds.ac.uk/~amt6xw/Distance Learning/CFD5020COMP/node17.html

如果能实现顶盖驱动MAC算法，simple也没有问题啊~可以先试试MAC，速度边界条件都是第一类边界条件，其余都为0.只有顶盖有速度。

@程迪 请问linux如何并行计算，在ubuntu下没有问题，但是在linux下并行，显示这样的错误，希望您指教

OpenFOAM-6(23) mpirun -np 2 interFoam parallel -------------------------------------------------------------------------- The value of the MCA parameter "plm_rsh_agent" was set to a path that could not be found: plm_rsh_agent: ssh : rsh Please either unset the parameter, or check that the path is correct

以前简单看过，好像是研究表面活性剂的方法。大概是在三维曲面上用类似FVM的方法离散控制方程并求解。

Examinination of bubble collisions and coalescence in bubbly flows
D. Broder and M. Sommerfeld
还有这一篇，这好像都是小的那个做的，大的具体尺寸我也不知道，不过什么时候能在做实验不好说，实验仪器快被转移了，大的那个。。

Direct numerical simulation of reactor two-phase flows enabled by high performance computing

捕获.JPG

Simulations of gas-liquid dispersions by four-way coupled E-QBMM with aggregation and breakage processes Abstract Results Main References Job: Data processing Test cases debugging and running Writing paper in English if you are willing to Requirements: Benefits: You can get around 30000 yuan bonus from your host (if you write the manuscript). You can use this paper to get promoted faster (if you write the manuscript). Another 10000 yuan from my side If I write the paper, I can support you 40000 yuan as well. However, I dont have much time to write. So I prefer you writing the paper and to be the first author. Project duration

@东岳 谢谢！

@ibelief 在 颗粒表面化学反应的文献 中说：

基于简单碰撞理论煤粉燃烧动力学模型的研究

非常感谢

moment method对于nucleation和condensation的模拟简单吗？

即为简单。但是溶解/蒸发数值上要复杂些。2017年我在HZDR的时候同事用fixed pivot method模拟增长，我用moment method，相比来看moment method要更简单，并且不需要处理相空间的drift：https://www.sciencedirect.com/science/article/pii/S1004954118301174/pdfft?md5=84216a8842ce3fd6d1a99d4f327c1ee1&pid=1-s2.0-S1004954118301174-main.pdf

带摩擦的运输方程是 u_t = lambda * u_x = eps * u_xx, 通过变化, 我们得到了, 0_1545168405877_3bd01f52-6748-4ece-bac9-e55eeef96e15-image.png ,其中D是导数矩阵, D2 是二阶导数矩阵
然后, 我们对M进行求特征值(EW), 得到了这种图片, 我想知道这是为什么呢?
0_1545168538892_dba5eac1-7474-4ce6-b089-77d3d8903c03-image.png

压差是分析很多工程设备的关键参数，在之前，主要都是针对多孔介质的Ergun模式来处理这一块。
也跟试验对比很好。但是，最近新开发的一个模型，虽然与试验的结果对比比较好，但感觉总是有问题。
于是不知道CFD中静压值到底是怎么得出来的？
对于一块多孔介质好像比较好处理，但是如果这个时候多孔介质表面沉积了灰，这个时候这块的整个压降会升高，然后灰越堆越多整个表面的压降都会升高。然而，这个局部压降和表面灰的分布十分相关，造成整个表面的这个灰分的压降相差一千多Pa。而最后一般的检测压降点，都是压差变松器的两个接口，一般都是设置在容器外壁的。但是如果我模拟检测点和试验一样的时候，这个时候的压差比我的表面计算压差平均值要小得多。（压差是1200多Pa，但是我的多孔介质有800多Pa，加上平均粉尘层的压降有640多Pa，合起来应该有1400多Pa，而不是1200多Pa.）
这就产生了一个问题，这个通过N-S方程求解的压力项到底是怎么处理的?
当他通过一层多孔介质的时候，多孔介质内部有压力源项，而在外部没有，那么当他气体从多孔介质外向内流通的时候，他的这层压力和速度是怎么计算的？如果这个时候我再添加一个源项在多孔介质层，而这个源项在多孔介质外是没有的，那么这层额外的源项是不是也会通过增加流体的静压来表达呢？（通过伯努利方程似乎很好理解，但是N-S是怎么实现的，他的这个\Delta P 是怎么获得的？）

DEM加速比很高
CFD好像一般

二维的vorticity只有一个方向 实际上从涡的形成来看 二维似乎并不适合一些涡量的模拟
但是对于一些简单的流场 实际上2D具备非常大的优势 他的处理简单 很多算法很好自编程实现
并且网格可以完全mapped 最大的优势 还是在于后处理的一目了然
这就是为什么大多数理论研究还是集中2D 但是2D可能不是太适用Fluent Fluent主要接轨工程实际问题

各位专家好：
在多孔介质模拟的过程中，我发现了一个问题，虽然在发表文章中省略了。但是确实是一个直观的问题，我的模型，RANS方法中，不管采用哪一种，在多孔介质的近壁面处的拟合的速度与实测值总是具备一定的差距。这种差距我分析认为是在较高阻力降的多孔介质环境下，气流在流经多孔介质表面时，一部分流体会因为高阻降而产生流向变化，形成比较复杂的近壁面流场。而在这个近壁面流场中，因为Fluent默认是采用的多孔介质域或者是多孔跳跃面，应该是对这层近壁面没有处理的，依然按照充分发展层处理，造成这层近壁面的流速（甚至是流向）与实测值都有较大区别。
虽然能够较好的拟合实时压降（Endo based equations），但是如果在近壁面流场不符合实际的话，那么很显然对于最后的表面覆尘的预算难以符合实际情况。不知道大家对于这种多孔介质的近壁面流是否有过研究，或者能够有比较好的解决思路？欢迎大家提出宝贵的想法和建议，感谢，再次感谢！
下图为提出新模型对于滤管外壁的模拟和实测结果对比：
0_1544421734755_compare.png