你是指的取断面为边界入口吗？这样的话可以采用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

Hi,

这个问题很简单，举例：

1核计算需要24小时；

2核计算需要12小时；

4核计算需要6小时；

8核计算需要4小时；

16核计算需要3小时；

32核机选需要2小时；

64核计算需要1.5小时；

很明显，核多了会减小计算时间，但是在超过某个限制之后，就不会是线性的。并且，并行计算和CPU的通信关系很大，有些玩家自己组建的通过网卡连接的工作站，一个千兆网卡是必须的。

你好，

这个其实区别不大，最简单的当然是deb安装。但是由于网络的问题国内经常安装不上。有关安装，我写了个非常详细的说明请参考：OpenFOAM安装说明

Yes,

不过需要你自己定义拉格朗日求解器。目前twoPhaseeulerfoam是欧拉方法下的双欧拉模型。如果需要添加一个无质量对流体流动无影响的颗粒，仅仅做一个颗粒跟踪。有俩种方法：

简单的捆绑dpmFoam中的拉格朗日粒子，我之前在pisoFoam中实现过这个，

通过paraFoam后处理来进行，这个在cfd-online有人讨论过，我当时也尝试实现不过没成功，目前还不确定到底paraview能不能实现这个方法。