@浪迹天大
需要一个动网格求解器，然后添加动网格字典，就是动网格~

@vivian
错误显示迭代步数超出了。

我采用了debug模式在rhoCentralFoam.C的215行加入断点，想输出e,结果屏幕输出了大量信息，我应该如何输出才能得到想要的场信息呢？

e是个场，你要输出场的话，就会输出所有的场的值，不清楚你要看什么，也可以单独输出场的最大值最小值看看越界否。如果一定要查看整个场的信息，输出到log再看，这样log文件会变得很大。

rhoCentralFoam &> log

如果我想梳理出程序运行的主线，搞清楚从场的离散点到通过离散方法，获得最终求解的代数方程矩阵，我认为最好的办法是从建立场开始，关注场中存储的离散点的信息变化。那么请问如何输出存储点压力，面通量信息的矩阵呢？

从离散的矩阵系数如fvScalarMatrix组建我们常规的稀疏矩阵不是一俩行代码可以搞定（估计不到100行？）因为openfoam里面存储的是稀疏矩阵，并用指针+指针。所以，可能得用点时间。

更简单的方法，是回到问题本身，很少是因为矩阵系统解不出来的原因，而是设置有问题导致发散。

So what is your question?

@MX
可以通过壁面函数nutkroughwallfunction来实现

对于一个面，要计算体积流量，另一种方法是通过phi乘以面积来获得。

对于面速度矢量，压力和每相流量是怎么给出的?

面速度=fvc::interpolate(U)是插值过来的。不需要后处理，再求解器里面已经处理好了。所以显然的

• phi=fvc::interpolate(U) & S，
• Q=fvc::interpolate(U) & S * mag(S)，
• phia=alpha1f * fvc::interpolate(U) & S，
• Qa=phia * mag(S)

多谢反馈。手边目前没有linux系统。你用的什么浏览器？没有这个标志？

@cfd-green

maxCo取决于时间格式、流动状态、网格质量、求解模型。

显性格式maxCo可以最大为1；
隐形格式maxCo可以适当扩大如2；
规则几何和流型库朗数可以适当扩大；
多相流库朗数可以最大为1，levelset可以为0.5或者更低；
动网格要确保点移动求解器每个时间步不能过大引起负体积；等等

所以是个很复杂的话题，并没有统一标准。可以逐步加大库朗数，知道发散为止。然后再降低。:sunglasses:

泊松方程如果给定Neumann边界条件，需要给定一个附加条件。在NS方程中，这个附加条件就是就是要保证速度边界条件是守恒的。

例如：我们有1D压力方程：
\begin{equation}
\frac{\partial }{\partial x} \frac{\partial p}{\partial x}=k
\end{equation}
其解为：
\begin{equation}
p=\frac{1}{2}kx^2+C_1 x+C_2
\end{equation}
如果给定Neumann条件，$p'(0)=\alpha$，$p'(1)=\beta$
带入后我们有：
\begin{equation}
C_1=\alpha, C_1=\beta-k
\end{equation}
这个就是附加条件，在此情况下，$C_2$无解，因此在压力均为Neumann边界条件的时候，需要设定参考压力，并且满足公式3. 公式3在NS方程中即为adjustphi()的目的。

从代码看这个函数只修正了边界上的流量，没有碰内部场，这样有什么意义吗？就是为了强制边界上总的通量是守恒的

修正边界条件确保压力泊松方程有解。

它不检查存在给定压力边界（p.needReference()为false）的算例，是默认这种算例能最终达到边界通量的守恒？

如果存在Dirichlet边界条件，则p.needReference()为false，因此不需要满足附加条件。

@CFDci

我同意你的观点，但是les模拟和网格关联性很大。也可以在大尺度下做les模拟，比如，同一个工况，不同的网格尺度，执行les模拟，在这种情况下，过滤的涡是不同的:cheeky:

这个是不是更简单点

好像目前还是最后一点用的比较多。模拟都要和实验对比，才能说明你算得准。
模拟和模拟对比的文章我没怎么遇到。

模拟终究逃脱不了实验啊！

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的点，