想问一下，针对复杂几何，大家都用什么软件画然后导入啊？

还没关注过这个问题 :xinlei:

对于使用KomegaSST模型作低雷诺数的湍流管流模拟，当y+值小于1时，是不是可以认为网格足够精细，无需再细化？
另外k omega nut三个量大概在哪个范围内合理？我的结果大概是以下：

k 10的负3次方

omega 最大50w+，最小500左右
nut 10的负10到负5次方

这里整理了一些 找对象 的诀窍:xinxin:

应该是Sb.boundaryField() 返回const reference
Sb.boundaryField().ref() 返回 reference
根据我的经验是这样,如果错误的话把编译错误贴上来
认真看看代码很简单的,编译器的错误已经告诉你所有信息了!

《A Thorough Description Of How Wall Functions Are Implemented In OpenFOAM》一文2.6节描述了nutWallFunction，nutAtmRoughWallFunction中nutw公式与之相同，其他量有改变，我再看看。

最近计算一个方腔流动遇到一些问题，这个流动沿流向（x方向）是周期性边界，yz方向是壁面，流动固定x方向的平均速度，自动求解压力梯度，在采用rans计算时是正常的，可以计算出一个合理的结果

pressure gradient = 0.179097

但是当采用PANS（自己编译的湍流模型）计算的时候，压力梯度会变得非常小，导致计算结果很差，

pressure gradient = 0.0903267

自己编译的模型可以在channel395中正常计算，但是在该算例中出现压力的降低，不知道是什么原因引起的

[2] [1] ##00 [38] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) at ??:? [38] #1 Foam::sigFpe::sigHandler(int) at at ??:???:? [1] [2] ##11 Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHandler(int) at ??:? [38] #2 ? at ??:? [2] #2 ? at ??:? [1] #2 ? in "/lib64/libpthread.so.0" [38] #3 ? in "/lib64/libpthread.so.0" [2] #3 ? in "/lib64/libpthread.so.0" [1] #3 ? in "/lib64/libm.so.6" [38] #4 ? in "/lib64/libm.so.6" [2] #4 ? in "/lib64/libm.so.6" [1] #4 ? in "/lib64/libm.so.6" [38] #5 ? in "/lib64/libm.so.6" [2] #5 ? in "/lib64/libm.so.6" [1] #5 ? in "/lib64/libm.so.6" in "/lib64/libm.so.6" in "/lib64/libm.so.6" _pmiu_daemon(SIGCHLD): [NID 00126] [c0-0c1s15n2] [Tue Jul 9 10:43:27 2019] PE RANK 38 exit signal Floating point exception _pmiu_daemon(SIGCHLD): [NID 00082] [c0-0c1s4n2] [Tue Jul 9 10:43:27 2019] PE RANK 1 exit signal Floating point exception [NID 00126] 2019-07-09 10:43:27 Apid 8697814: initiated application termination

I have tried it again, the same problem occurs.
Last time i have used shear gradient boundary, and this time i just use zeroGraident. However, ....

I guess. the problem is the signal floating exception???

Tks

@队长别开枪

谢谢您的建议，感觉这个方案可行。 我尽早试试，给您答复。。。

谢谢

把你的新库通过controlDict挂载一下，在controlDict里面添加

libs ( "libnew1.so" "libnew2.so" );

有没有老哥用过SPARTA这个软件，有没有什么好的入门读物推荐呀？谢谢啦（不知道帖子该放哪就发到这里了）

@夏雨天 你把7.55-57三个方程加起来，就会发现他们用的是C-N格式，三者加和是7.59。在书上179页最底下也说啦 :duang:

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

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

@东岳 我这边主要想做岩土工程灾害的模拟，比如泥石流、滑坡这种，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随便取