OpenFOAM

4.5k 主题 26.2k 帖子
  • forAll boundary

    3 帖子
    1k 浏览

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

  • 关于nutkAtmRoughWallFunction中nutw公式的疑问

    2 帖子
    1k 浏览

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

  • 压力驱动的方腔流动

    1 帖子
    697 浏览

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

    pressure gradient = 0.179097

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

    pressure gradient = 0.0903267

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

  • bubbly flow printStack error

    6 帖子
    2k 浏览
    [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

  • 9 帖子
    3k 浏览

    @队长别开枪

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

    谢谢

  • 怎样在waves2Foam中增加新的自定义波浪类型?

    3 帖子
    1k 浏览

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

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

    1 帖子
    595 浏览

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

  • 关于分步法和PISO算法

    5 帖子
    3k 浏览

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

  • 数值模拟时间与实验时间不符合的问题

    5 帖子
    3k 浏览

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

  • 2 帖子
    801 浏览

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

  • 关于CFDEM学习资料

    3 帖子
    1k 浏览

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

  • DNS中间结果疑惑点

    5 帖子
    2k 浏览

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

  • 1 帖子
    728 浏览

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

  • 2 帖子
    931 浏览

    在.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 error

    1 帖子
    729 浏览

    最近在模拟有传热情况下额气泡上升问题,求解器为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

  • 液滴动态接触角如何设置

    3 帖子
    2k 浏览

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

  • interFoam中设置两相的初速度

    7 帖子
    2k 浏览

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

  • 模版函数编译报错

    3 帖子
    920 浏览

    还真是,大意了@东岳

  • 湍流模型

    3 帖子
    1k 浏览

    谢谢。

  • SCI中的OpenFOAM求解器

    1 帖子
    772 浏览
    icoFoam twoPhaseEulerFoam
  • 看完东岳老师写的DNS案例的疑问。

    7 帖子
    2k 浏览

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

  • openfoam里怎么写求Nu面平均的方程

    3 帖子
    1k 浏览

    @东岳 谢谢 岳老师

  • CFD数值模拟方法分类

    1 帖子
    811 浏览

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

    谢谢各位!

  • 有关Openfoam在超算并行运算的问题

    18 帖子
    5k 浏览

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

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

  • multiphaseInterfoam 边界条件设置

    1 帖子
    782 浏览

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

  • 3 帖子
    1k 浏览

    @东岳 哈哈,好吧:xiexie:

  • 并行运算出现错误

    2 帖子
    784 浏览

    等图 :laile:

  • 如何在FiledAverage里设置得到Taylor length scales

    2 帖子
    787 浏览

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

  • 不可压计算雷诺数的计算问题

    4 帖子
    2k 浏览

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

  • interfoam+大涡模拟

    3 帖子
    1k 浏览

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

  • OF新手询问----顶盖驱动流的网格查看问题

    5 帖子
    2k 浏览

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

  • compressibleInterFoam测试算例

    1 帖子
    784 浏览

    两个比较好的可用来测试得实验。气泡溃灭。把他包含在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

  • 网格小引起的壁面处压力问题

    1 帖子
    705 浏览

    采用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); }

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

  • fluent mesh转化为Openfoam网格的问题

    3 帖子
    2k 浏览

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

  • openfoam两相流—波浪与结构物相互作用!

    8 帖子
    2k 浏览

    @guohuiqun 我的手机号,18971612958,微信:IAM315

  • openfoam网格导入问题

    3 帖子
    1k 浏览

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

  • 2 帖子
    969 浏览

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

  • 5 帖子
    2k 浏览

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

  • 1 帖子
    841 浏览

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

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

    结论:AMD的EPYC,真yes。应该点击一下

  • 2 帖子
    907 浏览

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

  • reactingTwoPhaseEulerFoam 空泡份额在哪求解的

    2 帖子
    926 浏览

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

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

  • 讨论 Fluent计算效率优于OpenFOAM

    30 帖子
    18k 浏览

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

    你需要深挖一下

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

    捕获.JPG

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

    捕获2.JPG

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

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

  • 2 帖子
    1k 浏览

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

  • 2 帖子
    947 浏览

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

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

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

  • simpleFoam 求解出错,

    2 帖子
    983 浏览

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

  • 合作一篇模拟,写一篇小文

    9 帖子
    4k 浏览

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

  • rhoCentralFoam 为什么做成半隐式?

    6 帖子
    2k 浏览

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

  • mapFields运行很慢吗

    1 帖子
    519 浏览

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

  • 请教一个代码IOdictionary_H

    4 帖子
    2k 浏览

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

  • In solid-water flow

    6 帖子
    1k 浏览

    @kimyIn 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