OpenFOAM

• 

  OpenFOAM初级入门建议（2020年更新版）
  • 东岳  

  38
  38
  Posts
  48142
  Views

  X

  @东岳 谢谢东岳老师
• H

  请问在RANS模拟中，可以使用自定的速度曲线作为入口条件吗？
  • hotairballon1997  

  4
  4
  Posts
  21
  Views

  S

  codeFixedValue
• J

  CFDEM安装
  • Jacob Z  

  8
  8
  Posts
  157
  Views

  F

  @Jacob-Z 在 CFDEM安装 中说： @fubianhanshu CFDEM®coupling 3.8.0=LIGGGHTS® 3.8.0+OF 5.x 好的，谢谢！我也尝试下，以后一起交流
• C

  在质量守恒方程中添加源项对湍动能的影响
  • CFDngu  

  11
  11
  Posts
  69
  Views

  C

  @CFDngu 我不做波浪所以不清楚一般大家都怎么做，所以我觉得看一下湍流粘度比值会有帮助。RANS计算的话，对一般的流动而言这个比值应该较大，如果比值很小那说明其实没有太多湍流的作用。也许真的像你们说的，其实用层流就没问题。 网格我是想问你有没有做过无关性分析之类的，另外这种场里面突然有一个小点突然算出来很大的量以至于发散，有可能是网格正好碰到了什么奇怪的数值点。当然究其根本肯定是程序哪里有什么问题，但有时候换一套网格有可能（极小概率的哈）能避开这个问题。
• 

  200万网格管道流qDNS并行算力测试
  • 东岳  

  25
  25
  Posts
  1527
  Views

  K

  @东岳 这是我44核心的log文件的最后几个时间步：之前的迭代次数也差不多： Time = 9.8 Courant Number mean: 0.733368 max: 1.27353 GAMG: Solving for p, Initial residual = 0.0998056, Final residual = 0.00082233, No Iterations 6 time step continuity errors : sum local = 4.80691e-07, global = 1.99762e-20, cumulative = 1.59798e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000123775 GAMG: Solving for p, Initial residual = 0.0809999, Final residual = 0.000584259, No Iterations 6 time step continuity errors : sum local = 3.45394e-07, global = 5.63939e-20, cumulative = 1.60362e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120435 GAMG: Solving for p, Initial residual = 0.0189678, Final residual = 9.4483e-07, No Iterations 39 time step continuity errors : sum local = 5.51633e-10, global = 3.45318e-20, cumulative = 1.60707e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119638 ExecutionTime = 170.51 s ClockTime = 171 s Time = 9.9 Courant Number mean: 0.73339 max: 1.27886 GAMG: Solving for p, Initial residual = 0.100098, Final residual = 0.000824461, No Iterations 6 time step continuity errors : sum local = 4.82558e-07, global = 8.66441e-20, cumulative = 1.61574e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000123823 GAMG: Solving for p, Initial residual = 0.0809731, Final residual = 0.000577791, No Iterations 6 time step continuity errors : sum local = 3.42103e-07, global = 4.07868e-20, cumulative = 1.61981e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120471 GAMG: Solving for p, Initial residual = 0.0189619, Final residual = 9.31189e-07, No Iterations 39 time step continuity errors : sum local = 5.44568e-10, global = 7.2382e-20, cumulative = 1.62705e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119674 ExecutionTime = 172.22 s ClockTime = 172 s Time = 10 Courant Number mean: 0.733413 max: 1.27933 GAMG: Solving for p, Initial residual = 0.0999699, Final residual = 0.000818608, No Iterations 6 time step continuity errors : sum local = 4.79761e-07, global = 8.13808e-20, cumulative = 1.63519e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00012387 GAMG: Solving for p, Initial residual = 0.081007, Final residual = 0.000576888, No Iterations 6 time step continuity errors : sum local = 3.41793e-07, global = 5.45899e-20, cumulative = 1.64065e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120508 GAMG: Solving for p, Initial residual = 0.0189675, Final residual = 9.02577e-07, No Iterations 39 time step continuity errors : sum local = 5.28191e-10, global = 1.77208e-20, cumulative = 1.64242e-17 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00011971 ExecutionTime = 173.99 s ClockTime = 174 s End 这是我32核心的： Time = 9.8 Courant Number mean: 0.733369 max: 1.27363 GAMG: Solving for p, Initial residual = 0.0998053, Final residual = 0.000728134, No Iterations 6 time step continuity errors : sum local = 4.26347e-07, global = 1.16708e-20, cumulative = 5.86486e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000123776 GAMG: Solving for p, Initial residual = 0.0810209, Final residual = 0.00054928, No Iterations 6 time step continuity errors : sum local = 3.24835e-07, global = 4.55333e-20, cumulative = 5.9104e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120434 GAMG: Solving for p, Initial residual = 0.0189497, Final residual = 8.88034e-07, No Iterations 38 time step continuity errors : sum local = 5.18253e-10, global = 2.24608e-20, cumulative = 5.93286e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119638 ExecutionTime = 167.35 s ClockTime = 167 s Time = 9.9 Courant Number mean: 0.73339 max: 1.27894 GAMG: Solving for p, Initial residual = 0.100093, Final residual = 0.00073128, No Iterations 6 time step continuity errors : sum local = 4.28756e-07, global = 1.40838e-20, cumulative = 5.94694e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000123823 GAMG: Solving for p, Initial residual = 0.0809929, Final residual = 0.000542798, No Iterations 6 time step continuity errors : sum local = 3.21522e-07, global = -2.66223e-20, cumulative = 5.92032e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120471 GAMG: Solving for p, Initial residual = 0.0189464, Final residual = 8.9418e-07, No Iterations 38 time step continuity errors : sum local = 5.22767e-10, global = 2.26167e-20, cumulative = 5.94294e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119674 ExecutionTime = 169.06 s ClockTime = 169 s Time = 10 Courant Number mean: 0.733413 max: 1.27934 GAMG: Solving for p, Initial residual = 0.0999826, Final residual = 0.000732193, No Iterations 6 time step continuity errors : sum local = 4.29802e-07, global = 1.57268e-20, cumulative = 5.95866e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00012387 GAMG: Solving for p, Initial residual = 0.0810305, Final residual = 0.00054485, No Iterations 6 time step continuity errors : sum local = 3.22957e-07, global = 2.11478e-20, cumulative = 5.97981e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000120508 GAMG: Solving for p, Initial residual = 0.0189583, Final residual = 8.97181e-07, No Iterations 38 time step continuity errors : sum local = 5.24979e-10, global = 2.66036e-20, cumulative = 6.00642e-18 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00011971 ExecutionTime = 170.82 s ClockTime = 171 s End
• 

  如何导入*.dat文件的网格
  • wuyukai  

  1
  1
  Posts
  15
  Views

  

  请问一下，有哪位大佬用过datToFoam这个命令吗，*.dat文件需要写成什么格式呢？
• 

  OF5版本以后重心坐标与绝对坐标的相互转换的问题。
  • 上级  

  7
  7
  Posts
  128
  Views

  

  上面那个公式引用是公式（8）
• 

  新编译的湍流模型求解算例的时候epsilon方程发散
  • strive_DUT  

  3
  3
  Posts
  37
  Views

  

  @东岳 我重新编译试过了，还是这个问题呢
• 一

  方腔中自然对流中的K、Omega等边界条件的初始值如何确定？
  • 一颗鸭蛋  

  2
  2
  Posts
  39
  Views

  

  就是从定义来求解，k还好说，epsilon和omega牵涉到一个特征长度的问题，这个取起来就不容易了，只能多调试。
• 

  网格整个运动的边界条件问题
  • wuyukai  

  1
  1
  Posts
  27
  Views

  

  运用pimpleFoam，求解一个平板的运动问题，想让整个网格一起运动，程序顺利运行了，但是结果明显是错误。猜测是边界条件设置的问题。想请教一下各位大佬，这种情况下，边界条件如何设置？或者是否有其他方面的问题。
• S

  of2006安装完找不到paraview和paraFoam
  • Samuel-Tu  

  5
  5
  Posts
  28
  Views

  S

  奇怪了，虽然用起来没什么问题，但是打开的时候还是　会有warning Invalid $PV_PLUGIN_PATH and -plugin-path= not defined No supplementary ParaView/OpenFOAM reader modules Using builtin reader: paraFoam -vtk Created temporary 'platePimple.foam'
• S

  paraview后处理的清晰度不够
  • Samuel-Tu  

  10
  10
  Posts
  103
  Views

  S

  @xpqiu 还真是，我用paraview5.6.3可以完美显示出来了
• S

  of2006的setExprBoundaryFields似乎很好用
  • Samuel-Tu  

  2
  2
  Posts
  32
  Views

  S

  链接文本，研究了一下用法，大家可以看看
• 

  《OpenFOAM用户指南》勘误
  • 东岳  

  27
  27
  Posts
  5943
  Views

  B

  李老师太客气了。免费看了这么久，作为新手其他的也没办法出力，能找几个小bug，也算是贡献微薄之力。 159页5.4.4中的“参见图5.12”应为”图5.11“。 我上贴最后一条写错了。不过新版UserGuide里5.4.5小节后增加了一小节网格细化功能的介绍，所以李老师更新时这些小bug应该会注意到。 161页5.4.6，“图5.14有一个例子”，应为“图5.13”。5.4.7，“参见图5.15”应为“参见图5.14”。 163页5.4.7，“边界层划分通过snappyHexMesh”应为“snappyHexMeshDict”。表中“minMedianAxisAngle 中轴点角度 130”应为“90”。 164页5.4.8， 表遗漏了minTetQuality关键字。（是因为V5版本后引入新功能取代该关键字的缘故吗？）
• S

  of 2006安装
  • Samuel-Tu  

  3
  3
  Posts
  64
  Views

  S

  @bestucan 看了官方的，主要就是怕需要一些依赖的包。。安装一下试试，看看新功能爽一下
• 

  pbm模型
  • zhanghan  

  2
  2
  Posts
  48
  Views

  

  何止有啊，简直一大堆，搜搜 OpenQBMM。
• 

  编译新的湍流模型遇到no match for ‘operator<’
  • strive_DUT  

  4
  4
  Posts
  53
  Views

  

  @东岳 是的李老师，p2是一个体标量场，0是整型，那这个地方应该怎么处理呢，求指点
• 胖

  求教对于一个四周固壁的开敞水面边界条件设定问题，谢谢~（仅仅只有水，不涉及空气）
  • 胖小帅  

  2
  2
  Posts
  40
  Views

  F

  @胖小帅 速度是不是要设置成slip?你试试
• C

  请问mesh.magSf() 怎么在内部网格使用呀
  • CarmeloSun  

  3
  3
  Posts
  47
  Views

  C

  @东岳 李老师。我是要在壁面附近一层网格区域加源项，每个网格需要加的质量计算的公式，需要乘以壁面这个网格的面积，再除以网格的体积，请问我需要怎么写呀
• 小

  Ubuntu18.04安装openfoam3.0.0报错
  • 小考拉  

  3
  3
  Posts
  52
  Views

  小

  @bestucan 感谢前辈的回复，目前问题已经解决 在编译OpenFOAM之前，需要先编译ThirdParty，但OpenFOAM WiKi中的Ubuntu18.04版安装教程里面好像并没有“编译ThirdParty这一步”，希望以后大家注意吧
• 

  time step continuity errors
  • zhanghan  

  8
  8
  Posts
  2695
  Views

  

  OpenFOAM中有pbm群体平衡模型的求解器吗？？？求指点
• L

  小白对于openfoam k-epsilon湍流模型程序编写的困惑
  • lllwonderliquid  

  11
  11
  Posts
  2589
  Views

  

  @东岳 谢李老师！
• 液

  关于颗粒属性的输出
  • 液固两相流  

  6
  6
  Posts
  97
  Views

  F

  @液固两相流 我是在计算域内初始化一定量的颗粒，输出计算过程中每个颗粒的位置和速度信息，这个可以实现吗？
• 

  请问OF里有没有融化和凝固相关的求解器
  • 麦迪文  

  7
  7
  Posts
  109
  Views

  

  @dzw05 谢谢
• 

  求问一个公式在OpenFOAM中的表示方法？
  • strive_DUT  

  9
  9
  Posts
  441
  Views

  

  @东岳 谢谢李老师！
• B

  关于CreateTurbulenceFieldsLES疑问
  • benqing  

  2
  2
  Posts
  42
  Views

  

  箭头是对指针的成员的引用，找到这个类的定义：turbulenceModel.H epsilon()函数声明在205行，是个虚函数，所以在 turbulenceModel.C 里没有定义，而是在使用其作为模板创建新类的时候再各自定义 在上面那个窗口里搜 epsilon() 函数： 看你的情况是用的哪种定义，其实里面 epsilon() 函数定义大差不差， 你用的 LES，咋找上层流的了，你找的那个确实是层流的对于 epsilon() 的一个实现，在我发的那个图片里也有 脉动速度都用上了，那应该不是层流了，速度分成脉动速度和平均速度这就是用了湍流模型了
• L

  湍流模型问题请教：
  • leilei  

  6
  6
  Posts
  139
  Views

  闻

  @leilei 我用RSM做过旋流场气液流模拟，收敛性还可以，只是刚开始微调一点就行，但是计算速度比很多其他RANS慢很多。
• 

  关于传热边界条件的问题
  • win  

  1
  1
  Posts
  36
  Views

  

  我使用的buoyantBoussinesqSimpleFoam，边界条件设置如下 出现的问题如下
• S

  压力泊松方程中A和HbyA怎么说比较好
  • Samuel-Tu  

  4
  4
  Posts
  100
  Views

  S

  最近看了以下jasak的论文，找到了相关描述
• Y

  闲来无事，做一个Axial Mesh的教程吧
  • yuan_neu  

  6
  6
  Posts
  2912
  Views

  Y

  @random_ran 这位大佬，我也尝试了您说的这个方法，是不是这样产生的网格都是大小一致的，没办法局部加密网格？
• F

  关于MPPICInterFoam的周期性边界的疑问？
  • fubianhanshu  

  5
  5
  Posts
  54
  Views

  F

  文献有点大，上传不了，麻烦东岳老师下载吧（1）http://dx.doi.org/10.1016/j.ijsrc.2017.09.006（2） 10.1002/2013JF002911文献里都是用DEM方法模拟的，绝大多数文献里MPPIC都是用于气固或气液，想着试试用它做液固，跟别人的对比看看如何，增加些模型什么的，目前刚开始接触欧拉拉格朗日法，没什么思路
• H

  在集群上foam-extend-4.0以及extend-bazaar中流固耦合包的安装问题
  • huningdong  

  9
  9
  Posts
  107
  Views

  S

  @雾雨魔理沙 不好意思哈，回复的比较晚...我的qq是891931269
• 

  DNS出口边界条件
  • 东岳  

  17
  17
  Posts
  519
  Views

  

  @Calf-Z-DNS 你好有关于对流传热的问题想要咨询下你，就是对流出口边界条件的设置问题，是只对U使用的，还是全部都用呢？
• R

  并行效率疑问
  • random_ran  

  49
  49
  Posts
  20138
  Views

  刘

  @刘雄国 在 并行效率疑问 中说： 大牛说 下一代cfd 需要 大规模计算机网格 + Algorithm + Coding 整体考虑 我觉得很有道理 我说错了, 下一代cfd 需要 大规模计算机网格 + CFD Algorithm + CFD Coding + physical网格 整体考虑
• S

  提取边界为stl
  • Samuel-Tu  

  7
  7
  Posts
  82
  Views

  

  感谢分享
• 

  并行计算与单核计算的结果不能完全吻合吗？
  • DY大世界  

  5
  5
  Posts
  83
  Views

  

  @东岳 哈哈，明白，谢谢李老师~
• 

  请教各位前辈一个问题啊，OpenFoam可以同时计算机翼外流场、机翼结构和机翼中燃油吗？
  • Jerry  

  5
  5
  Posts
  61
  Views

  

  谢谢李老师啦，我先看一下这几个求解器，了解一下思路
• Y

  fvMatrix 的系数正负求解稳定性
  • yfclark  

  7
  7
  Posts
  64
  Views

  Y

  @东岳 ok，明白了，谢谢东岳老师
• S

  blockMESH文件里可以写函数吗
  • Samuel-Tu  

  2
  2
  Posts
  59
  Views

  R

  @Samuel-Tu 記得好像是用巨集 #calc
• S

  关于流固耦合求解器fsiFoam的问题
  • shepherd  

  5
  5
  Posts
  61
  Views

  S

  @benqing 会不会是因为固体域没有时间步，所以网格不发生变化？
• 

  CFDEM安装测试
  • 东岳  

  7
  7
  Posts
  3327
  Views

  

  今天在编译liggghts和cfdem的时候也出现如下报错： 结果证明需要如下依赖包： sudo apt-get install libvtk6-dev 另外好像OF8没法完成和liggghts的耦合，我用的是5.x的。
• 

  关于paraview后处理颗粒的问题
  • 上级  

  8
  8
  Posts
  155
  Views

  

  @闻久STU 谢谢哦，您这个方法我也记在小本本上了，以后肯定能用上
• 

  OpenFOAM中的Lax-Friedrichs方法
  • 东岳  

  1
  1
  Posts
  70
  Views

  

  http://dyfluid.com/rhocentralfoam.html 参考这个 把Lax-Friedrichs方法植入进去了 耗散的很 不过是个OKS好材料
• S

  shm的stl可以不封闭吗
  • Samuel-Tu  

  2
  2
  Posts
  44
  Views

  S

  简单验证了一下，可以用shm，没问题
• 

  如何CFD自学成才
  • 东岳  

  5
  5
  Posts
  224
  Views

  闻

  @walten 在 如何CFD自学成才 中说： @闻久STU 乐观点， 六年 这样才学得深入 不能提 长期抗战是肯定的了
• Z

  在颗粒流中的动量转移大概有多少？
  • zhangxc0223  

  1
  1
  Posts
  39
  Views

  Z

  抱歉问的问题可能太泛了，具体举个例子来说。比如说在直径D=0.01m的空气管流中，速度Uf=15m/s，雷诺数Re=10000。颗粒的斯托克斯数St=10，直径是大概dp=50×10^(-6)m，密度为ρp=1000kg/m^3，颗粒-空气的质量比是Φ=0.1，双向耦合，只考虑颗粒曳力。 在这种情况下，我查看OpenFOAM计算得到的kinamaticCloude:UTrans文件中，动量转移都是在10^(-14)量级（单位是kg×m/s），请问这合理吗，为什么只有这么小。。？ 据我所知，这个UTrans文件表征的就是颗粒→流体相的动量转移，也就是颗粒相通过源项Sp，施加在流体相的曳力作用。不知道这个理解对不对？
• S

  修改文件内容的shell
  • Samuel-Tu  

  8
  8
  Posts
  73
  Views

  S

  最近多次使用foamDictinary感觉真是太用了，自己在文件里面加的keys也能识别出来并修改
• W

  拉格朗日颗粒喷射
  • wangtao121  

  1
  1
  Posts
  40
  Views

  W

  各位大佬问一下有关拉格朗日颗粒喷射的问题，我如何往喷射模型中（比如manualInjection模型）动态的传递喷射点坐标，坐标每一个时间步都会发生变化，怎么实现内部传递呢。
• 

  请问，MPPIC这个方法的particle distribution function的运作
  • zhe  

  11
  11
  Posts
  169
  Views

  

  @zhe 在 请问，MPPIC这个方法的particle distribution function的运作 中说： 固相是由particle distribution function表示的，并且固相的相分数也是通过它进行求解的。那么这个方程在OpenFOAM里有体现吗？ 没有，你说的这个方程是一个非常基本的控制方程，其实就是玻尔兹曼方程，很多模型都是从这个模型上展开对他进行封闭的。比如 双流体模型，就是粒子分布函数传输方程取2次矩的模型。 矩方法就是取高阶矩的模型。 MPPIC在这个方程上对碰撞进行BGK封闭。可以更简单的处理碰撞。看这个文献的方程2 A model for collisional exchangeing as /liquid/solidfluidizedbeds 《无痛苦N-S方程笔记》4,8节也有讨论。这个方程没人去解。自由度太多了。 可以彼此融会贯通 左右逢源一下
• S

  of输出continuity error残差
  • Samuel-Tu  

  3
  3
  Posts
  39
  Views

  S

  @马乔 有没有什么办法输出这个误差呢，一般大家是怎么判断pisoFoam或者pimpleFoam是否收敛呢? 我发现误差降到１e-5以下就很困难了。。。所以我是设置的监测点，发现关注的变量变化不大了就认为收敛。。