@东岳 谢谢东岳老师的关心，后来还是有点进展，不过进度不是很大，现在所走的路就是修改求解器，文章还没有出来。等出来之后我会详细说下这个问题。

看起来只不过是模拟准确不准确的问题哦？可能下图模拟出来就是这样啊

@东岳 我这边主要想做岩土工程灾害的模拟，比如泥石流、滑坡这种，DPMFoam可能不合适。

原实验论文里面把这个case称作‘subturbulent’，雷诺数是6020，T0 = 297.15 K,
P0 = 0.0927 MPa, D = 0.0274 m，V0=3.73 m/s，文中好像有讲到流动有趋于层流的趋势。

在做多物理场耦合时用到弹性力学边界条件，对位移向量，需要对不同方向的分量分别施加不同的边界条件，如下图这种：在一个方向上为fixedValue，另外两个方向为fixedGradient
TIM截图20190705111025.png
想请教一下大家：在OF里有没有类似的边界条件，如果自己实现大概是什么思路

在.H文件中找到了：
Reference:
D.M. Hargreaves and N.G. Wright, "On the use of the k-epsilon model
in commercial CFD software to model the neutral atmospheric boundary
layer", Journal of Wind Engineering and Industrial Aerodynamics
95(2007), pp 355-369.

最近在模拟有传热情况下额气泡上升问题，求解器为compressibleInterFoam。当连续计算时，出现如下错误，有没有大神指导一下？不胜感激！！
Lookup solver for TFinal
smoothSolver: Solving for T, Initial residual = 0.0076371, Final residual = 8.40008e-09, No Iterations 5
[60] #0 Foam::error::printStack(Foam::Ostream&)[36] #0 Foam::error::printStack(Foam::Ostream&)[15] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
at ??:?
[15] #1 Foam::sigFpe::sigHandler(int)[36] #1 Foam::sigFpe::sigHandler(int) at ??:?
[60] #1 Foam::sigFpe::sigHandler(int) at ??:?
[15] #2 ? at ??:?
[36] #2 ? at ??:?
[60] #2 ?Lookup gradScheme for nHat
Cache: Calculating nHat, originating from alpha.water event No. 624
in /lib64/libc.so.6
[15] #3 in /lib64/libc.so.6
[36] #3 in /lib64/libc.so.6
[60] #3 ???Lookup interpolationScheme for interpolate(nHat)
at ??:?
[15] #4 Foam::H2O::rho(double, double) const at ??:?
[36] #4 Foam::H2O::rho(double, double) const at ??:?
[60] #4 Foam::H2O::rho(double, double) const at ??:?
[15] #5 Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::Es(double, double) const at ??:?
[36] #5 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy> > >::he(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
[60] #5 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy> > >::he(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
[15] #6 Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::T(double, double, double, double (Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::)(double, double) const, double (Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::)(double, double) const, double (Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::*)(double) const) const at ??:?
[36] #6 Foam::twoPhaseMixtureThermo::correctThermo() at ??:?
[60] #6 Foam::twoPhaseMixtureThermo::correctThermo() at ??:?
[36] #7 at ??:?
[60] #7 at ??:?
[15] #7 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy> > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool)?? at ??:?
[15] #8 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy> > >::correct() at ??:?
[36] #8 __libc_start_main at ??:?
[60] #8 __libc_start_main at ??:?
[15] #9 Foam::twoPhaseMixtureThermo::correctThermo() in /lib64/libc.so.6
[36] #9 in /lib64/libc.so.6
[60] #9 at ??:?
[15] #10 ??? at /home/abuild/rpmbuild/BUILD/glibc-2.22/csu/../sysdeps/x86_64/start.S:121
at /home/abuild/rpmbuild/BUILD/glibc-2.22/csu/../sysdeps/x86_64/start.S:121
at ??:?
[15] #11 __libc_start_main in /lib64/libc.so.6
[15] #12 ? at /home/abuild/rpmbuild/BUILD/glibc-2.22/csu/../sysdeps/x86_64/start.S:121

@东岳 cfd online的回答说有时候用随机值。。。。我正在测试算例，随便用的0.01:zoule:
尤其让我无法理解的就是无滑移边界uwall=0，uTheta随便取

@Irisch 谢谢回复！我查一下研究下。

还真是，大意了@东岳

谢谢。

icoFoam twoPhaseEulerFoam

@东岳 好的谢谢！
我想问下圆周沿线要怎么画？:chouchou: ，另外流向的平均就是流向沿线每个点的u+相加起来然后除以点数吗？
那y方向上的位置要选多少个？谢谢！

@东岳 谢谢 岳老师

据小弟了解：湍流数值模拟方法包括DNS，RANS，DES和LES等，其基本控制方程都是N-S方程耦合自身封闭方程。 请问：
1.对于数值模拟中SPH方法，和上述方法的归类情况如何？
2.还有那些常见的数值方法？

谢谢各位！

你运行完decomposePar后的显示是什么？手动分块容易造成有的核分了很多网格，有的分了很少，分的多的跑的很慢，所以其他核还是要等它算完了才能继续。我一般用scotch，不用输入任何参数，openfoam自动会把分区做成1、每个分区cell数目接近；2、分区之间的processor面尽量的小，从而减小数据传输的需求。

batch系统调度感觉你需要再细看一下？node是结点数，你24个分块，如果比如超算上一个结点就有24个核（4个6核cpu），那你node数量只需要1吧？这个超算是什么样的，每个结点有几个cpu，每个cpu几核呀？

微信图片_20190627084835.jpg
outlet该用什么边界条件好一些？
zerogradient的话，油层会上升，然后排出计算域，freestream的话油层会下降。正常的物理现象是油层会波动
然后outletphaseinletoutletvolicity 的话，开始的时候计算的挺好，到后面回流特别严重，

@东岳 哈哈，好吧:xiexie:

等图 :laile:

就是如何从LES计算epsilon吧？目前据我所知只能通过epsilon的定义来求。你需要写个后处理程序