原实验论文里面把这个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的定义来求。你需要写个后处理程序

@我是河滩 不需要显式设置密度，自己心里知道就行了。

@mohui 您也用interfoam算过水渠流动问题吗

谢谢老师，已成功，继续学习案例

两个比较好的可用来测试得实验。气泡溃灭。把他包含在OKS
几何简单，数学模型复杂，验证算法利器

Numerical simulation of a collapsing bubble subject to gravity

Bubble collapse near a solid boundary: a numerical study of the influence of viscosity

Numerical investigation of an underwater explosion bubble based on FVM and VOF

QQ图片20190623213234.png

采用interFoam修改的求解器计算算例，遇到网格太小带来的壁面处压力问题（图中为第一个输出时间步（约0.0001 s）的压力分布，出现了不正常最大和最小压力）。
算例的壁面（BOTTOM）为锯齿状壁面，图中所示部分为倾斜向上的一小部分，壁面附近的网格精度约为1微米，压力场初始为0，压力和速度壁面处边界条件如下：
0001.png
p_rgh

BOTTOM { type fixedFluxPressure; value uniform 0; }

U

BOTTOM { type fixedValue; value uniform (0 0 0); }

大家有没有遇到过这样的问题，有什么解决思路吗？

谢谢，格式为 fluentMeshToFoam xxx.msh -writeZones

@guohuiqun 我的手机号，18971612958，微信：IAM315

李老师，我试了试，mergemeshes没弄出来: 上传文件，网格文件50m，网站限制发不了，我发您邮箱了。。。到弄了好久还没弄出来:mihu: 您看看用什么方法能合并下

发现错误了，是文件名字的问题，文件名本来是大写字母，我写成小写了。警告对结果没影响，或者使用fin函数进行输入，就没有警告了。

@Hungryandfool 您说的我也想过，但我主要是应付不可预测的爆发性计算需求，超算中心这种流程繁琐，而且灵活性不足。

当时买的时候，搜索了这里。没有搜索CFD-ONLINE，是我的错误。

可以看CFD-ONLINE，相关讨论。 可以看chip-hell，相关讨论。 可以看知乎，相关讨论。

结论：AMD的EPYC，真yes。应该点击一下。

用滤镜calculator可以计算一些，不知道能否满足你的需求

说的对，这篇文章中我没写，你可以看看这个 https://www.jianguoyun.com/p/DWnVKFIQ9s3ZBhjYsrkB

等我回去有空更新下中文的解析 多谢反馈

@东岳 在 讨论 Fluent计算效率优于OpenFOAM 中说：

你需要深挖一下

给你举个例子，下面是运算的一个log文件，时间步长非常小，2e-5，算2天之后，结果完全可以是正确的，但我觉得有问题

捕获.JPG

在自适应调节时间步长的情况下，时间步长会自动跳转满足CFL标准。其实对于我这个，如果在代码里加上几行，输出Ur的话，会发现下图中全场计算域中，就这么一个地方Ur的值特别大，导致局部库郎数过大。

捕获2.JPG

我建议你看看你的算例是不是也是这样，我怀疑跟我差不多

想请教一下遇到这种情况应该怎么办呢？特别算vof的时候，界面某些点出经常出现很大的速度，导致时间步长小的无法接受。

@yishanbb2018 3个太少了，我这面见过文献是一两百个

@史浩 在 multiRegion计算中，interface附近场量后处理显示异常 中说：

将固定路径的数据通过OpenFOAM中后处理工具以txt的格式输出来，然后通过其他软件作图即可。毕竟paraview中的云图还能凑合看，写文章编辑也发现不了什么

我也遇到过这个问题，我一般用sample工具然后用gnuplot作图，目前我只知道这么做了:yes:

我找到原因了，就是不同的网格，对应参数的求解方法，以及求解方法的smoother也要做相应的修改，我这个算例的速度求解没问题，压力求解方法的smoother 出了问题，导致计算终止。谢谢各位老师关注

这个文章已经建立合作，研究将于7月末开始做，计划10月份投出去。

@东岳 是啊我也看到了，改成压力基说是在low mach时更好。感觉各种拼凑。。

两个计算域从5000网格映射到300000网格，网格不相容，重叠的只是一小部分，选的mapNearest方法，运行特别慢，十多分钟才跑完。
感觉不应该这么慢吧，哪里出问题了？

@刘雄国 对，我是说没有功能上的作用

@kimy 在 In solid-water flow 中说：

interphase slip agebrism

Never heard of that

That's internal field, a good approximation is good for convergence. From your depiction I feel you need to learn more about the general CFD knowledge. You can definitely feel that the lacking of CFD theory is holding you back now

@hangsz 非结构有限体积也有高精度格式

@wwzhao 原来是这样，icem软件的六面体导出网格前都要转成非结构，应该就是为了适应fluent中的非结构计算。

@刘雄国 我知道这个网站，之前就看了一些，现在是OpenFOAM的版本差的有点多，之前的cloudFunctions是单独编译的，这样就很好弄，但是现在的是和liblagrangianIntermediatemy.so放在一起编译的，包含很多程序，我不知道在/src/lagrangian/intermediate/Make/file 里面哪个地方是指定的对于cloudFunctions 的编译。

在Openfoam中的时变边界条件中的多项式设置如图
a244f322-80ca-4bf8-99ab-670cd38326c9-image.png
如何实现在多项式上添加时间段的限制
例如：0＜t＜1，velocity = 0.1 + 1.3x^2 + 2.7x^3
1≤t＜2，velocity = 2x^2 + 22x^3
这种应该如何添加时间段限制

拉格朗日taylor-green涡中运动，好酷炫

image.jpg

@东岳 试了一下，等待结果中~

多谢，已修正方程6、8、11、12、13 :xiexie:

回头试了下，用ICEM画就可以转换成功。

是不是用户的问题，切换其他用户试试？

不太清楚，可能需要去看一下代码了。

@东岳 wo我把除了inlet和outlet的patch都改为stick，然而一样的结果，一样的报错信息。