@搬运工不好当 说：

在某文献里看到用OpenFoam划分高速列车流体网格，六面体

目前OpenFOAM处理这种非常简单，snappyHexMesh是一种自动网格划分的生成器。特别擅长画复杂几何。不过从个人经验来看：

你不需要过分追求细节； OpenFOAM不是纯六面体网格；

因此要是发文章我更倾向与使用ICEM

能不能在水槽内部形成一个固体边界？不紧挨着边界处

If you dig deeper, 你会发现ThirdParty里面不是OpenFOAM基金会开发的。例如：paraFoam是paraview的插件，paraview属于Kitware，scotch属于法国的一个...看这里：链接文本。所以称作第三方包。

你是指的取断面为边界入口吗？这样的话可以采用timeVaryingMappedFixedValue；

@搬运工不好当 说：

得到一个特定时间步结果

不得不说如果调整了时间格式为何还要期望获得一样的数据结果？不能期望每一个格式获得同样的数据。

Hi,

这里我并没有发现错误？

Hi,

你的twoPhaseProperties没有定义，如果你是在interFoam中调用mu，在interFoam的createFields.H中已经定义了：

Info<< "Reading transportProperties\n" << endl; incompressibleTwoPhaseMixture twoPhaseProperties(U, phi);

所以可以这样写：

volScalarField& muc(twoPhaseProperties.mu());

试一下。我没有测试，有问题再联系。

@mark Hi mark,

错误提示，缺少twoPhaseMixtureThermo.H文件。可以看看这个文件有没有在你的Make/options中包含

Hi,

我曾经调用过Fortran库来计算矩阵特征值，不过后来我采用Alberto的库而不用Fortran了。在OpenFOAM里面调用Fortran混编需要安装gfortran，可以在终端下通过：

sudo apt-get install gfortran

来安装。安装后下载你所需要的Fortran源代码（目前我还不会写fortran，我用的是现成的包），里面通过makeFile来制定规则，在终端通过make编译就可以。编译出来的是.so文件。OpenFOAM可以直接调用。

也可以在./Allwmake里面把make写为脚本。

使用入口边界条件，turbulentInlet。
或者延长入口使用，maped，类似于前面加一个自循环的pipe效果。
或者使用，timeVaryingMapped，提供入口出的数据。

很早之前用过ssh client之类的东西，可以利用图形界面进行登录，文件拷贝之类的东西，不过运行程序还是需要命令行。

不过还是使用ssh命令行靠谱啊，时间长了效率会比图形界面高的。

这个倒是一个问题，因为最近做到了流化床的燃烧，那么里面肯定是可压缩的，因为燃烧放热，密度变化等，气体为可压缩的。CFDEM在2.9.0好像是可以支持压缩的，但是感觉：流速比较低，处于低马赫数，气体压缩效应并不明显，所以我看绝大多数的冷态CFD-DEM计算都用不可压缩的。我们这里有做高马赫数（激波）等，这时候应该要考虑可压缩性了。

Hi，

你可以试试scotch方法，从他的介绍来说是自动的处理网格疏密问题。可以这样测试一下：分解之后用paraview看一下各个块怎么样。我一直用scotch分解，不过还没测试过。这几天有空我测试下。

@mark Hi,

等我有空把文档写完，那个暂时一句俩句说不清楚。

因为流体力学中常用的是前五个量纲，因此可以默认忽略最后两个量纲，详见openFoamUserManual_PFM第19页。

Hi,

目前我知道的有在windows下面有使用GPU加速计算的，例如Symscape，还有这个：RapidCFD。

前几天看到在CUDA上求解矩阵的算法，感觉很有意思：CUDA。我的理解GPU计算就是在CUDA上用CUDA支持的语言编程来求解矩阵？

暂时知道这么多

Hi，

采用MULES方法求解的相场不应该存在越界问题，不过守恒倒是有一些问题。能否详述一下？

OpenFOAM里面相场越界这个问题已经是一个很重要的问题了。从OpenFOAM-1.4开始OpenCFD公司（目前已注销）就着手采用MULES来处理标量场的守恒问题，并逐步完善比如添加LTS功能等。

@搬运工不好当

Hi,

PMV目前在海尔科学研究院有一套标准，国际上也有一套ISO9000x的标准。在这个标准中，除了公式的值都需要自己设定。然后在PMV框架下迭代求解。可以google一下PMV ISO标准。

Hi,

请注意epsilon的不能赋值为0，因为在kEpsilon方程中的湍流粘度定义为：

\begin{equation}
\nu_t=C \frac{k^2}{\varepsilon}
\end{equation}

明显的如果epsilon为0，湍流粘度不可以计算。

Hi，

我们用方程举例：加入你要求解矩阵：
\begin{equation}
A \phi = b
\end{equation}
直接解法是不能用了大家都知道。一般采用迭代解法。在迭代解法中，如果对矩阵操作一下：
\begin{equation}
A = M-N
\end{equation}
替换到公式(1)中，就是
\begin{equation}
\phi^n = M^{-1} N \phi^{n-1} + M^{-1}b
\end{equation}
这种操作就是预条件，依据不同的M，有不同的预条件方法。OpenFOAM里面采用不完全LU分解。

smoother更倾向是多重网格矩阵求解器的内容，如图：
0_1458467715003_捕获.JPG
多重网格这块暂时还没涉及，以后有涉及随时补充。

推荐看这本书：Numerical Computation of Internal and External Flows - Volume I - Fundamentals of Numerical Discretization