OpenFOAM

• 

  OpenFOAM初级入门建议
  • 东岳  

  34
  34
  Posts
  45933
  Views

  Z

  @东岳 老师写的太好了
• L

  openfoam计算时，显示chemkin文件报错，有人出现过类似错误吗？真诚请教！
  • lts  

  3
  3
  Posts
  17
  Views

  L

  谢谢！！已好
• C

  两相流中设置一面只有固相能感受到的墙
  • cxzhmg  

  2
  2
  Posts
  7
  Views

  

  欧拉欧拉？欧拉拉格朗日？欧拉拉格朗日比较好处理。欧拉欧拉的话，目前能想到的就是多孔介质。因为欧拉欧拉压力方程是公用的。速度可以单独设置无穿透。但是压力不能单独设置。 如果有大佬看到相关的sci，感谢分享
• 雾

  求教！foam-extend算例中提取力时报错
  • 雾雨魔理沙  

  2
  2
  Posts
  8
  Views

  

  你把nu放在括号外试试
• 1

  波浪与圆柱作用的动力响应
  • 19490  

  2
  2
  Posts
  28
  Views

  H

  这个估计很难了。你这个涉及流固耦合的问题。CFD这边做不鸟。但是你可以这么做，流场获得你需要的压强时程曲线，然后你在固体场中把压强加上去。这样是可以的。
• S

  OF波浪模拟
  • Samuel-Tu  

  6
  6
  Posts
  38
  Views

  H

  波浪得求解器很多得。一般三个，waves2Foam，olaFlow, ihFoam。我现在可能开始用olaFlow
• S

  并行计算能否直接引用到case的地址
  • Samuel-Tu  

  3
  3
  Posts
  10
  Views

  S

  去源代码里面找了一下caseConstant()对于串行和并行，都可以返回case里面的constant的地址，我试了一下，达到我想要的效果了
• 

  MPPICInterFoam速度公式與interFoam用的一樣嗎?結果很不同.
  • zhe  

  7
  7
  Posts
  19
  Views

  

  @东岳 感謝东岳老師的回覆.首先說個題外話,最近在英國這邊一個公衆號上看到您發的一片勸解博後選擇的文章.因爲是轉發的,起初不知道是您發的,後來越看越像您的口吻,最後在下面看到"轉自李东岳"的時候,心裏很是激動.謝謝东岳老師對後輩們的勸解! 我又做了一個用MPPICInterFoam在沒有粒子的情況下,繞柱(光滑壁面)的模擬.這次用了LES湍流模型.效果很好.並且並行運算時,LES比K-omega SST運行更快.大概是因爲這種情況更符合實際情況吧.跟實驗數據比較,甚至比之前interFoam做的還要好.圖中我做了MPPICInterFoam和interFoam的對比.上面的是MPPICInterFoam的.在柱後水流形成的湍流旋渦非常明顯.在不同高度的對比,效果也很好.至少說明,在水動力這塊兒,這個求解器配合LES,可以得到不錯的效果.當然,邊界條件也很重要. 有問題的大家一起討論啊.
• 

  文件夹大小 文件大小 查看
  • 东岳  

  2
  2
  Posts
  69
  Views

  

  删除算例目录下除了m开头的所有文件，注意备份！！注意会把constant文件夹内的也会删除 find . ! -name 'm*' -exec rm -f {} \;
• 一

  自己写的后处理的functionObjects怎么直接调用？
  • 一二  

  6
  6
  Posts
  17
  Views

  

  in file db/functionObjects/functionObjectList/functionObjectList.C at line 278 看起来你改了functionObjectList.C的代码？
• M

  请教扩散现象不随时间变化的原因？
  • merlinizlf  

  5
  5
  Posts
  17
  Views

  M

  @Samuel-Tu 非常感谢您的回复。讲的非常明白，我之前没有想到这些。其实，我想实现的就是一个热源点慢慢加热最终提高出口温度，或许我应该试试改成热源强度再试一试看。非常感谢您。
• S

  OF如何输出header
  • Samuel-Tu  

  9
  9
  Posts
  22
  Views

  

  感谢分享！！！收录在此 https://www.cfd-china.com/topic/3499/openfoam小代码
• S

  自动中止命令
  • Samuel-Tu  

  3
  3
  Posts
  15
  Views

  S

  @xpqiu 感谢，我去查一下!不过目前我用break也能顺利终止程序，并行和串行都可以
• Z

  请问怎么由3D湍流模型的速度场，得到在特定位置上的湍谱(turbulence spectra)？
  • zhangxc0223  

  4
  4
  Posts
  825
  Views

  X

  没有实际做过，但是这个链接提供的方法看起来很靠谱。https://www.physixfan.com/如何计算湍流能谱/
• 王

  可压缩浸没边界法的实现
  • 王慧博  

  5
  5
  Posts
  49
  Views

  王

  @东岳 非常感谢东岳老师！按照您的指导，我设置了新的边界条件immersedboundaryEnthalpyFvPatchField，其他部分的设置目前都正常，但是想要使用浸没温度边界上的函数来更新焓值时会报错，我的程序部分如下所示： const label patchi = patch().index(); wordList type = thermo.T().boundaryField().types(); if(isA<immersedBoundaryFvPatchScalarField>(type[patchi])) { immersedBoundaryFvPatchScalarField& TIb = thermo.T().boundaryField()[patchi]; scalarField & Tref = TIb.refValue(); } 之后程序就会报错， error: invalid initialization of reference of type ‘Foam::immersedBoundaryFvPatchScalarField& {aka Foam::immersedBoundaryFvPatchField<double>&}’ from expression of type ‘const Foam::fvPatchField<double>’ thermo.T().boundaryField()[patchi]; 应该是无法将基类转换为派生类，想和老师请教一下遇到这个问题要怎么解决呢，麻烦老师了！！！
• C

  fvm::ddt(rho, T)与rho*fvm::ddt(T)的区别
  • ChimneyFoam  

  3
  3
  Posts
  18
  Views

  

  fvm::ddt(rho,T) = $\frac{\p \rho T}{\p t}$, rho*fvm::ddt(T) = $\rho\frac{\p T}{\p t}$
• S

  function里的scalarTransport中resetOnStartUp
  • Samuel-Tu  

  1
  1
  Posts
  6
  Views

  S

  问题背景是在simpleFoam计算中添加被动标量运输，因此在controlDict里添加了function的scalarTransport。但是在设置参数的时候发现resetOnStartUp。 在源代码中的discription是Reset scalar to zero on start-up。 现在比较疑惑的是，如果打开resetOnStartUp: １．是在每个时间增量计算开始时把标量reset为０ ２．还是一整个计算结束后，再次重新计算时把标量设为０ 大家有经验吗
• S

  并行计算错误
  • Samuel-Tu  

  3
  3
  Posts
  18
  Views

  S

  经过一上午的排查，找到原因了。在我使用并行的时候，往scalarField里面加数据的时候超过了定义的大小
• 西

  如何在dynamicMeshDict里设置restraint，实现浮体运动状态的控制？
  • 西极若木  

  2
  2
  Posts
  13
  Views

  C

  有现成的，可以看看6DOF模块的代码，里面又控制方程，也有写好的阻力
• C

  如何将6DOF和动边界一起使用？
  • CFDngu  

  5
  5
  Posts
  26
  Views

  C

  @Sloan 感谢回复，overset是另一个路子，我这里没有涉及
• Y

  checkMesh fail with openfoam2.3
  • yfclark  

  7
  7
  Posts
  282
  Views

  S

  @东岳 这位是我后来一步一步检查，用checkMesh得到得具体的问题 如果感觉没有那么复杂的不规则边界，那么基本就是几何拓扑关系出了问题吧
• Q

  fvVectorMatrix Ueqn
  • qingdong_wang  

  5
  5
  Posts
  30
  Views

  

  pEqn里面用到了rAUc
• S

  分享: from faceId of patch to global index
  • Samuel-Tu  

  2
  2
  Posts
  10
  Views

  

  感谢分享！
• S

  local index to global index
  • Samuel-Tu  

  1
  1
  Posts
  6
  Views

  S

  用并行计算的时候，想通过某个face的local index得到其global index，查了一下似乎是这个函数faceProcAddressing，但是写入solver的时候调用不了这个函数。说没有定义此函数，如何调用这个函数呢。。ｗiki上说globalIndex=mag（faceProcAddressing[localIndex]）-1
• D

  如何设置一个入口边界条件：入口边界值为出口平均值的一半
  • D.Benjamin  

  22
  22
  Posts
  102
  Views

  S

  @东岳 东岳老师，decomposed后的patch里面的face是否还是按分解前的顺序排列呢。例如分解前的patch中有6个face编号为0，1，2，3，4，5。分解的时候刚好从中间分开了，分解成两个processor里面，则processor1和processor2里的face编号都变成了 0，1，2。假如processor2里面的face是未分解时候的4，5，6，那么现在face 0，1，2一定对应的是4，5，6吗，也就是说这个顺序会不会改变，会不会0，1，2对应成原来的5，6，4了
• S

  调用某个patch上的数据
  • Samuel-Tu  

  3
  3
  Posts
  17
  Views

  S

  @东岳 感谢!！!！!！
• P

  openFOAM7安装时配置环境变量一直失败的问题
  • pcxu  

  6
  6
  Posts
  18
  Views

  

  @东岳 踩坑踩多了而已
• H

  颗粒在切向的运动方向中的受力，要同时考虑曳力升力么
  • hotairballon1997  

  17
  17
  Posts
  29
  Views

  Z

  @东岳 好的，谢谢老师～～
• 史

  OpenFOAM多个并行算例同时运行导致计算效率下降
  • 史浩  

  25
  25
  Posts
  1631
  Views

  S

  @史浩 嗯嗯，好的，很详细了，谢谢您的回复！我看看mpi里的参数设置，琢磨下看看
• 炎

  俯仰翼型网格总是出现负体积求助
  • 炎冰灬灬  

  3
  3
  Posts
  17
  Views

  炎

  @东岳 这个是出现问题前的最后一步，网格变形太厉害了，我在考虑是不是可以用圆尾的好一些
• 

  请教如何划分才能得到这种网格？
  • DY大世界  

  3
  3
  Posts
  23
  Views

  

  @东岳 感谢李老师回复！圆柱周围的网格细化我明白，但是在整个水槽的自由液面处网格细化也是用snappyHexMesh做的吗，而且这里的自由液面区域和圆柱周围区域的衔接也很好，看上去像是一个整体加密区域。但是我总感觉自由液面处的加密像是通过blockMesh实现的，我在这里还是有点困惑，李老师您能详细的指点我一下吗。
• 玉

  codedFixedValue边界条件实现入口颗粒通量等于出口颗粒通量
  • 玉玺  

  6
  6
  Posts
  12
  Views

  玉

  @东岳 而且我也不要把其他场映射过去，只需要颗粒速度场
• Y

  变梯度边界条件
  • Yuchuan Guo  

  6
  6
  Posts
  14
  Views

  Y

  @东岳 @Hunahum ，我在cfd-online上看到这样一段代码，想尝试着写出来，不过总是报错，请问你知道问题在哪儿吗
• H

  动网格细化植入interFoam
  • Hunahum  

  3
  3
  Posts
  12
  Views

  H

  @东岳 谢谢东岳老师的关注，第一个报错应该是关于内存释放的，有时候还会报错free() invalid next size (normal).尝试了把EEqn.H能量方程中的relaxfactor项删除，就可以避免这个报错，但是会对计算结果有影响…… 还尝试过用gdb运行程序，就不会碰到这个报错，能量方程可以顺利求解（至少第一个时间步可以），但还是会提示下面那个网格数量不匹配的错误
• S

  请教如何只提取某处的液面高度时程图，就像实验放置在水槽中的wave height probe那样
  • Sloan  

  6
  6
  Posts
  39
  Views

  S

  @东岳 @xiezhuoyu 谢谢！
• S

  想请教下setFields不工作，很奇怪
  • Sloan  

  4
  4
  Posts
  35
  Views

  S

  @东岳 谢谢！
• P

  paraview如何衡量长度
  • pcxu  

  3
  3
  Posts
  34
  Views

  P

  @东岳 ，好的，万分感谢
• 西

  关于OpenFOAM里自带的waves算例的疑惑
  • 西极若木  

  6
  6
  Posts
  22
  Views

  

  @西极若木 邮箱：mechanics_xzy@163.com 微信：xzy11235813 多交流多交流 共同进步
• P

  关于interFoam dambreak算例的一点疑问
  • pcxu  

  2
  2
  Posts
  37
  Views

  S

  setFields里alpha.water=1的区域里就是water，性质是给你说的1项的。OpenFOAM基于FVM，一个cell里相分数alpha=1表示都是水，alpha=0表示空气。供参考
• S

  foam-extend的solidMechanics
  • Samuel-Tu  

  3
  3
  Posts
  22
  Views

  S

  @Samuel-Tu !就像这样，我显示的是位移，但是总有一层初始的膜（初始时刻的位移图）一样把我的结果给包住了，看不见我想要的结果。。( 图片地址)
• Z

  LES圆管湍流的速度曲线和DNS不符，请问问题可能出在哪？
  • zhangxc0223  

  8
  8
  Posts
  1175
  Views

  Z

  @霜染丹枫 入口用的是turbulentInlet，出口是inletOutlet
• Y

  多相流射流求解
  • yfclark  

  4
  4
  Posts
  69
  Views

  L

  此处的压力无反射边界条件是哪一个？
• S

  多版本OF共用paraview
  • Samuel-Tu  

  3
  3
  Posts
  26
  Views

  S

  @东岳 目前我是这么做的感觉不是很方便。讨论群里说应该添加paraview到PATCH这个环境变量的最前面。。没查到相应的命令
• 

  动网格和overset物体速度快速震荡导致崩溃的问题
  • tower  

  12
  12
  Posts
  70
  Views

  

  @CFDngu 谢谢建议！！
• 

  DMD动态模态分解
  • 红豆沙  

  1
  1
  Posts
  19
  Views

  

  有懂的可以交流下么
• 炎

  动网格的旋转运动边界方法求助
  • 炎冰灬灬  

  1
  1
  Posts
  14
  Views

  炎

  各路大神好，想请教个问题。使用dynamicMotionSolverFvMesh网格运动求解器和displacementLaplacian动网格求解器时，设定什么样的时均运动边界条件可以指定为旋转运动，让旋转的角速度随时间均匀变化并且可以中途改变旋转方向。比如：从0度开始顺时针匀速旋转到10度，随后又逆时针匀速旋转到0度，这样往复旋转。
• Y

  devReff()通过devRhoReff()实现其功能 （rho =1）？（OF4.x）
  • Yongbo  

  3
  3
  Posts
  15
  Views

  S

  我记得构造哪个class的时候就是默认了1
• 炎

  pimpleDyMFoam求解器发散求助
  • 炎冰灬灬  

  5
  5
  Posts
  19
  Views

  炎

  @羽之下 动网格用的自带算例wingMotion为模板的。目前问题找打了，是因为网格变形过大。
• 一

  从objectRegistry中无法获得输运字典问题
  • 一颗鸭蛋  

  7
  7
  Posts
  73
  Views

  一

  @马乔 我这个报错已经解决了，是流体区零文件下的alphat文件的壁面边界条件不能使用alphatJayatilleWallFunction函数。但是具体具体为什么不能使用我并不知道其中缘由。也许你可以调试一下边界条件。 该壁面函数我在单域方腔绝缘算例中是能够使用的，但是在流固耦合求解器计算的流固耦合算例中却会报类似上述错误。希望了解该函数的朋友能够讲解一番，不胜感激。
• 

  求指点---颗粒模拟
  • 红豆沙  

  1
  1
  Posts
  14
  Views

  

  研究流体携带的砂粒在一定的工艺参数条件（比如注入速率、携砂液粘度、砂子粒径、砂子密度、裂缝宽度、裂缝角度等等）下的运移和分布规律，除了使用DEM外，有没有什么较为便宜的方法，可以得到相类似的结果呢，DEM跑起来太。。。。。