2018年第一季度CFD/CAE类自媒体影响力分析

本统计只收录2018年第一季度的消息，统计日期为2018年3月30日。

相关数据有证可查（来源于各个公众号），本帖只做事实统计。

• 原创号：内容大部分为原创；很少转载；
• 营销号：内容大部分为转载；偶尔原创；

%(#ff0000)[平均每篇浏览量]

本指数统计公众号发布的头条平均浏览量。计算方法为：总浏览量/篇数

%(#ff0000)[总浏览量前三名]

本指数统计公众号发布的头条的总浏览量。计算方法为：总浏览量

第四名为工程事（浏览量为20343），第五名为风流知音（浏览量为10266），后略。

CFD界头条备案

CFD之道头条概览

ANSYS学习与应用头条概览

南流坊头条

工程事头条

$\alpha_1=1-\alpha_2$是用来保证alpha有界的，因为：
\begin{equation}
0<\alpha<1
\end{equation}
守恒体现在连续性方程中：
\begin{equation}
\frac{\partial \rho\alpha}{\partial t}+\nabla\cdot(\rho\mathbf{U}\alpha)=0
\end{equation}

of中多相流求解器中的alpha方程是套了循环

具体指的是什么？MULES么？

@小龙 OpenFOAM应该不存在相关现成的算例，你尝试过上面的代码了么？能否运行？:big_mouth:

@李东岳 为何我这一直提示证书错误？已经安装必要的库，界面可以打开，crack也是从你的链接卸载的，奇怪了。

@李东岳 还需要安装：

sudo apt-get install csh

否则会提示错误，Ubuntu16.04

:big_mouth: :big_mouth: :big_mouth:

@lsprxd 在 请教几个关于壁面函数的问题，谢谢 中说：

近 壁 区 高 Re数 采 用 的紊 动 各 向同 性 假 设 不 再 适 用,因而k-e紊 流 模式在 此 不 能 应 用,故 应 选 取 距 固 壁某 一 距 离 的 边 界作 为 计 算 边界。

文献出自哪里？

一般是设置了不符合物理的边界条件所致。贴边界条件看看。

还有在这种计算条件下应该如何考虑湍流模型呢？考虑到壁面附近主要为低速高粘的液体，可以认为是层流？而主要空间是较为高速且低粘的空气，会形成涡流影响破碎后的液滴运动轨迹，又还是应该选择湍流模型比如 SST， SAS 这种过渡的模型？求解释。

很难说，很难说具体用什么湍流模型，CFD中很多模型的选用都是很难提前知道，因此好多都是尝试各种模型进行对比。ke和kOmega对比，SST和SA对比，PIV和CFD对比，你应该庆幸幸好没有公认准确的CFD模型，CFD工作者才能有点事干。所以，你可以看看文献用的是什么模型，如果你模拟的数据要和文献对比，那么采用文献的模型，如果文献纯时间，碰巧你研究两相流的湍流模型，对比一下没准能发个Nature :cheeky:

OpenFOAM和我们的想法不太一样，回想有限体积法一般面上的值和周围的网格有关，因此OpenFOAM主要定义面的owner和neighber就足够进行任意迎风的插值格式了。因此你不会看到cell文件，你会看到owner和neighbour，如果想知道cell的信息，你需要知道包围网格的面，然后通过这些面确定cell的点，

看起来网格质量不错，不应该出现问题，你用的哪个求解器？log贴出来一下看看。

是否是这个图的意思？The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM® and Matlab®，page 666。一般这种用来求解大马赫数问题，或者那种那个所谓的全速流动。

@小龙 不，你只需要定义上表面的速度旋转就可以了，不需要别的参考系。

@温宝霖 恭喜。
你模拟的是什么问题？能否发个网格截图出来看看？听起来挺有意思。

@piteqiu 使用laplacianFoam，添加#include "fvOptions.H"，方程中添加源项（参考scalarTransportFoam），在fvOptions中定义codedSource。

src/fvOptions/sources/general/codedSource/CodedSource.H

你需要依据时间确定你的区域添加源项就可以了。可能需要一些网格的东西：mesh.C()返回网格位置，mesh().time()返回当前时间。

    if(RK2SSP_)
    {
        TRK2SSPSolve();
    }
    else
    {
        TEqnSolve();
    }

CFD计算的时候迭代非常多，如果进行10万次迭代，这种判断语句能增加多少计算时间？1s？10s？

这个图我给100分:cheeky: ，你求解的是固体导热还是对流导热？

@温宝霖
你出错的那个例子也是用层流模拟的:upset: 网格snappyHexMesh做的？上来就79万啊。

是的，就是这个意思。:cheeky:

错位网格：

• 进口速度给定，压力修正从进口速度得出
• 出口速度等于内节点速度，压力修正=0

请参考H K Versteeg and W Malalasekera第九章