像大家请教一个问题，算低速流体（10m/s以下）且求解域很小（1cm1cm1cm量级）当用可压求解器求解好像CFL不管设多小，算足够长时间都会发散，除开用不可压求解器以外还有什么其他解决方式么？

@李东岳 想明白了，多谢东岳大神！

是中国人，还没改国籍。

互联网的某个角落 or 提问。

@搬运工不好当 每个边界都要赋值，格式大概是
nonuniform List<scalar>
xxx
(
aaa
bbb
ccc
.....
);
xxx 要等于边界面元的数目，括号里的数字的个数要等于xxx。

另一个用法是vtk可以存储部分数据，而不是整个的场。大型算丽非常有效。

另外有些paravirw的后处理需要vtk做源source.

@dyj19901127 hi.

我觉得一个好的习惯就是禁止复制构造函数。除非你每个类都写，但一般也用不到。henry应该是采取不写的做法吧。

这个在c++ primer plus 和effictive c加加里面提到过。书并不在手边你可以看看那个类里面存在指针的时候，复制构造函数怎么泄露的。我更觉得这是个和编程有关的问题和cfd关系不大？

@搬运工不好当 你先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,

这里我并没有发现错误？