@搬运工不好当 你先wclean一下试试:）

@搬运工不好当 shm生成的网格有的时候非常的tricky，不错如果你采用shm成功的生成了mesh okay的网格，openfoam大部分求解器都能算收敛。of里面单相求解器已经非常稳健了。

如果对于现有几何，一定要生成一个非常好的网格，我建议你

简化几何 换专业的网格生成程序 慢慢的，充分的调节shm里面的网格控制参数，主要是网格密度，迭代次数，甚至删除边界层，这个步骤需要很长时间。

话说回来，shm对于复杂几何生成的是混合网格。虽然单相求解器对混合网格非常稳健，但某些多相求解器还算不了。

@hongfu2233 不，湍流粘度并不是一个直接求解的量。并没有湍流粘度的方程，因此不需要初始值。

@liliang2009223 我都那个丸不是很了解。这俩个互溶么？液体

@mark dgdt分析已经上传

对于第一类边界条件，OF里面壁面法向导数只有一介精度吧？

Foam::tmp<Foam::Field<Type> > Foam::fvPatchField<Type>::snGrad() const { return (*this - patchInternalField())*patch_.deltaCoeffs(); }

请问有没有二阶精度的sngrad?

请教一下，hooke/stiffness model怎么使用，有成功的吗？能否分享下。

@搬运工不好当 说：

在某文献里看到用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分解，不过还没测试过。这几天有空我测试下。