OpenFOAM

• 

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

  38
  38
  Posts
  46149
  Views

  X

  @东岳 谢谢东岳老师
• 

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

  1
  1
  Posts
  1
  Views

  

  大家好，请问有没有人知道MPPIC这个方法的particle distribution function具体怎么运作的?我查找了很多文献，也在OpenFOAM上看了好久的codes，但是没有说要解这个方程。那么这个方程到底在这个方法里的贡献提现在哪里啊？希望有人可以解答，谢谢！
• K

  time step of twophaseeulerfoam
  • kimy  

  3
  3
  Posts
  11
  Views

  K

  @东岳 Thanks for your reply. The simulation was not crashed but the time step is small, says 1e-5. I found most of researchers used a smaller time step even than mine in this phenomenon. Thus, it is usual to be like this? If this, I need a lot of computational time?
• 董

  paraView5.6不适用于openfoam-dev了？
  • 董力阳  

  8
  8
  Posts
  35
  Views

  R

  @东岳 我发现是我不会新建文件，而建立了文件夹。。导致报错。。。 第一次用linux，有点懵。
• 林

  能否在集群的个人文件夹内安装OF
  • 林之流风  

  8
  8
  Posts
  42
  Views

  林

  @bestucan 非常感谢老师的详细建议。因为集群之前就已经安装过OF，我是在个人账户里安装OF7版本的，所以可能就差这一个包。我已经使用了跑了计算，目前来看没什么问题
• 

  多面体网格转换
  • 土豆王子  

  1
  1
  Posts
  11
  Views

  

  用fluent Meshing划分的多面体网格转换成OpenFOAM格式出现错误 /--------------------------------------------------------------------------- | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: v1912 | | \ / A nd | Website: www.openfoam.com | | \/ M anipulation | | *---------------------------------------------------------------------------*/ Build : _f3950763fe-20191219 OPENFOAM=1912 Arch : "LSB;label=32;scalar=64" Exec : fluent3DMeshToFoam 000.msh Date : Jul 03 2020 Time : 09:25:27 Host : King PID : 514 I/O : uncollated Case : /mnt/d/OpenFOAM/run/test nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Reading header: "ANSYS(R) TGrid(TM) 3D, parallel 18.1.0" Dimension of grid: 3 Number of points: 34807 --> FOAM Warning : Found unknown block of type: "11" on line 16 Number of faces: 39160 Number of cells: 6053 PointGroup: 1 start: 0 end: 15. Reading points...done. PointGroup: 8 start: 16 end: 23. Reading points...done. PointGroup: 216 start: 24 end: 307. Reading points...done. PointGroup: 218 start: 308 end: 3281. Reading points...done. PointGroup: 59 start: 3282 end: 9500. Reading points...done. PointGroup: 61 start: 9501 end: 34806. Reading points...done. --> FOAM Warning : Found unknown block of type: "11" on line 34838 --> FOAM Warning : Found unknown block of type: "11" on line 34852 --> FOAM Warning : Found unknown block of type: "11" on line 34866 --> FOAM Warning : Found unknown block of type: "11" on line 35156 FaceGroup: 2 start: 39148 end: 39159. Reading uniform faces...done. FaceGroup: 7 start: 0 end: 3253. Reading mixed faces...done. FaceGroup: 62 start: 3254 end: 39147. Reading mixed faces...done. CellGroup: 63 start: 0 end: 6052 type: 1 --> FOAM Warning : Found unknown block of type: "316" on line 74613 --> FOAM Warning : Found unknown block of type: "71" on line 74632 Zone: 63 name: zone000:000--- type: fluid. Reading zone data...done. Zone: 62 name: interior-p-62 type: interior. Reading zone data...done. Zone: 7 name: geom-zone000:000---:7 type: wall. Reading zone data...done. Zone: 2 name: geom-zone000:000--- type: geometry. Reading zone data...done. Zone: 64 name: edge-p-64 type: boundary-edge. Reading zone data...done. Zone: 214 name: edge-p-214 type: boundary-edge. Reading zone data...done. Zone: 9 name: _zone000:000---:feature:9 type: boundary-edge. Reading zone data...done. Zone: 3 name: _zone000:000---:feature type: boundary-edge. Reading zone data...done. --> FOAM Warning : Found unknown block of type: "73" on line 78743 FINISHED LEXING Creating patch 0 for zone: 2 name: geom-zone000:000--- type: geometry Adding polyPatch for unknown Fluent type geometry Creating patch 1 for zone: 7 name: geom-zone000:000---:7 type: wall Creating cellZone 0 name: zone000:000--- type: fluid Creating faceZone 0 name: interior-p-62 type: interior faceZone from Fluent indices: 3254 to: 39147 type: interior patch 0 from Fluent indices: 39148 to: 39159 type: geometry patch 1 from Fluent indices: 0 to: 3253 type: wall --> FOAM FATAL ERROR: Illegal cell label -1 in neighbour addressing for face 39148 From function void Foam::polyMesh::initMesh() in file meshes/polyMesh/polyMeshInitMesh.C at line 67. FOAM exiting 
• 

  关于OF并行分区的问题
  • aiweimo  

  11
  11
  Posts
  515
  Views

  X

  @东岳
• 

  openfoam 自适应网格通量修正设置
  • 南风  

  2
  2
  Posts
  11
  Views

  

  我这边后面主要涉及到水汽两相流，空化相关的问题，谢谢
• 

  MPPICInterFoam 的Ueqn中到底該不該再乘以alphac
  • zhe  

  11
  11
  Posts
  42
  Views

  

  @东岳 好的,多謝回覆.我會接着再試試重新定義rho,等您有時間了,希望能多指導.
• 

  二维自适应网格
  • supersoldier  

  21
  21
  Posts
  9148
  Views

  

  你好，自适应网格这块dynamicMeshDict字典关于通量修正的的设置，里面的这些项表示的意义不太明白，什么情况下需要修正那些东西呢，应该如何设置，非常感谢，或者有推荐的资料学习也可以
• B

  请教用openfoam模拟大气边界层流动的问题
  • beidu  

  7
  7
  Posts
  24
  Views

  

  你直接用openfoam-v2006就好，不需要用5.0版本的编译了 参考最底部：https://www.openfoam.com/releases/openfoam-v2006/solver-and-physics.php
• 小

  编译求解器出现 undefined reference
  • 小考拉  

  7
  7
  Posts
  38
  Views

  小

  @bestucan 感谢前辈的热心回复，目前问题已经解决，问题出现在类phaseChangeTwoPhaseMixture的make文件的files文件里面， LIB = $(FOAM_USER_LIBBIN)/libphaseChangeTwoPhaseMixtures //会报错 LIB = (FOAM_USER_LIBBIN)/libtwoPhaseChange //随便叫一个名字均可，均可编译通过，但不能是libphaseChangeTwoPhaseMixtures 自己也是偶然发现的，但也比较困惑为什么这么神奇 不过最终问题还是解决了，很感谢前辈耐心的解答 ，自己后面也应该好好补一点C++知识
• 小

  求解相方程越界问题
  • 小考拉  

  2
  2
  Posts
  8
  Views

  小

  @小考拉 alpha方程如下：
• 明

  关于waves2foam造波衰减问题
  • 明月奴516  

  5
  5
  Posts
  81
  Views

  

  补充一下学长 @臧伟 的内容 2.关于时间步长的问题，考虑到波浪信息在单个时间步长内传播距离不宜太大，一般取T/500左右。以及上文提到的网格尺寸问题，可参考： 查晶晶. 基于OpenFOAM的数值造波与消波模型及其应用[D]. 上海交通大学, 2011. 如果单纯造波消波的话，时间项采用二阶离散格式也有利于抑制波幅衰减的问题。 3.如果是阻尼消波的话，阻尼区域建议为波长1~2倍。
• 液

  OpenFoam中如何设置总压进出口
  • 液固两相流  

  3
  3
  Posts
  23
  Views

  

  感谢分享！！
• 液

  多重坐标系下如何考虑颗粒在旋转区域的科式力？
  • 液固两相流  

  3
  3
  Posts
  13
  Views

  液

  已成功植入emoij
• A

  interFoam的曲率计算
  • anubis  

  2
  2
  Posts
  12
  Views

  

  或许可以把方程写出来对比一下嘛？
• A

  SHM 网格拉伸变形
  • aufalpen  

  1
  1
  Posts
  8
  Views

  A

  请问各位大神，SHM里哪个参数导致生成的网格如图所示的拉伸变形？多谢！
• 夜

  自己编译安装libccmio库
  • 夜阑烟寒  

  7
  7
  Posts
  1421
  Views

  

  https://portal.nersc.gov/cfs/visit/third_party/libccmio-2.6.1.tar.gz 连接挂掉了
• 花

  请教压力测量问题
  • 花之水聿  

  1
  1
  Posts
  11
  Views

  花

  采用动态压力传感器测量实验对象的压力脉动，在与理论结果进行对比时，应该对比理论结果的静压还是绝对压力还是其它类型的压力？
• 桎

  自适应网格加入了dynamicMeshDict但是不起作用
  • 桎梏  

  8
  8
  Posts
  43
  Views

  桎

  @nuaa_zp 你好。我的问题已经解决了 我的网格是400w，maxCells是200w，所以不起作用 我把maxCells调大就可以正常使用了，谢谢这两天的回复
• 人

  使用集群并行计算时出现错误
  • 人来人往  

  7
  7
  Posts
  360
  Views

  U

  @人来人往 您好，请问您是石油大学华东这边的吗，最近也准备在集群上安装OF，完全没有头绪，可以请教下吗？谢谢
• U

  DPM并行计算
  • upc_ngh  

  13
  13
  Posts
  62
  Views

  U

  @东岳 好的，我去了解一下，谢谢各位老师的指点
• 李

  snappyHexMesh定义了多域网格（MultiRegion）失败。
  multiregion • • 李子橙  

  6
  6
  Posts
  31
  Views

  

  感谢分享
• 吴

  如何修改并编译kOmegaSSTBase？
  湍流模型，编译 • • 吴优  

  3
  3
  Posts
  14
  Views

  吴

  @东岳 谢谢老师的回复，我是按照在turbulentTransportModels.C加入了 #include "kOmegaSST2.H" makeRASModel(kOmegaSST2); 然后在kOmegaSST2.H中 #include "kOmegaSST2Base.H" 实现这一编译过程，但是还是会报错，似乎要增加一个新的Base需要修改的地方非常多。 我目前暂时采取的解决方案是直接在修改原base，增加了一个Switch，在算例的constant/turbulenceProperties中设置true来启动新功能，这样也算保留了原模型，打算后面如果找到更好的方法我再改过来
• 

  一个关于 Axisymmetric Subsonic Jet 的问题
  • King_RVM  

  31
  31
  Posts
  79
  Views

  

  @King_RVM 感谢分享！！
• Z

  praView5.4.1版本中的矢量速度显示疑问
  • Zhoujq  

  1
  1
  Posts
  8
  Views

  Z

  在调节是矢量的时候： 主图如下： 其中分为3个部分： 1部分 在Vector中有U，也就是说里面所有带箭头的都是表示速度； 2部分 但是在这里，并没有看见速度的选项，只有一个GlyphVector 3部分 在数标图中显示的为GlyphVector Magnitude 请问下，这里的GlyphVector Magnitude是否就是U速度的意思呢？ 谢谢~
• 

  openfoam自适应网格使用Q准则加密失败
  • 南风  

  4
  4
  Posts
  25
  Views

  

  @东岳 谢谢李老师。我这边加进去了
• Q

  OpenFoam如何施加垂直边界平面的边界条件，及用其他软件生成网格导入至OF中的一些疑问
  • qz  

  8
  8
  Posts
  3199
  Views

  Y

  您好，我最近也在使用这个边界条件，有个问题想要请教您，关于refValue的取值是直接等于具体问题的取值吗？？期待您的回复
• 

  ParaView并行体渲染速度跟工作站硬件配置的关系？
  • aiweimo  

  8
  8
  Posts
  43
  Views

  

  我发现每次都是这个OpenGLProjectedTetrahedraMapper卡进程，然后CPU利用率就飙升。GPU主要在旋转视角的过程中（不松鼠标）才会占用到。有限网格量，开很多核并不能跑满CPU，但是可以跑满主频，所以相同浮点运算能力下，可能提升主频对于体渲染更有用。
• 炎

  将翼型网格替换fsiFoam自带算例网格后运行报错问题
  • 炎冰灬灬  

  6
  6
  Posts
  17
  Views

  炎

  @东岳 好的，谢谢东岳大神
• 

  一个关于使用SnappyHexMesh控制边界层生成的小程序
  • wg0632  

  8
  8
  Posts
  81
  Views

  

  感谢分享！
• 

  rhoSimpleFoam求解器算例
  • 东岳  

  1
  1
  Posts
  37
  Views

  

  3.1 二维叶栅 3.2 三段翼亚音速流动 3.3 亚音速轴对称射流 有关算例可以在这里讨论
• C

  并行计算结果与串行结果差异很大的问题
  • ChimneyFoam  

  7
  7
  Posts
  29
  Views

  

  buoyantPimpleFoam大体上可以参考rhoPimpleFoam 不过既然你们都说要buoyantPimpleFoam 我可以写一写 http://dyfluid.com/rhoPimpleFoam.html 网上也有一些其他的资料应该，不过自从我自己写自己的推导之后，就没看过别人的 这个求解器只不过就是可压缩求解器+传热，没啥大的特殊
• 炎

  fsiFoam自带算例报错问题。
  • 炎冰灬灬  

  7
  7
  Posts
  31
  Views

  炎

  @东岳 好的，谢谢东岳大神
• 我

  如何给定网格内部面上的流动参数？
  • 我是河滩  

  15
  15
  Posts
  73
  Views

  我

  @xiaolin 是的
• 宝

  如何在给定总压入口的同时给定速度的方向
  • 宝丁  

  4
  4
  Posts
  390
  Views

  液

  @宝丁 请问这个边界条件中指定的压力值是总压吗？
• N

  compressibleInterFoam运行中读入热物性参数报错
  • nuaa_zp  

  4
  4
  Posts
  3
  Views

  N

  @nuaa_zp 自己的计算模型及温度边界条件设置： 报错信息：
• C

  关于buoyantPimpleFoam的几个问题
  • cccrrryyy  

  16
  16
  Posts
  87
  Views

  C

  @东岳 你忙~我最近再看Ferziger和Peric那本书的Compressible Flow一章，了解的也更多了。看代码之余还是应该多温习理论！
• C

  作东岳老师OF编程入门case时程序引用自定义库出现函数未定义的引用问题
  • CUGB_Zhou  

  3
  3
  Posts
  20
  Views

  C

  @东岳 问题已解决，感谢东岳老师回答，是我自己代码输入有问题，不好意思。
• H

  请问DPMFoam的source reference是哪篇文章？
  • hotairballon1997  

  6
  6
  Posts
  16
  Views

  L

  @hotairballon1997 在 请问DPMFoam的source reference是哪篇文章？ 中说： @东岳 所以您知道DPMFoam的参考文献是哪篇嘛？ 方程3、4是用在OpenFOAM里DPMFoam的表达式嘛？我看这里有考虑体积分数α，好像普通的颗粒方程没有α这一项吧？ 不考虑颗粒相的体积分数也是可以的，针对稀疏颗粒流就可以不用考虑
• Z

  关于forceCoeffs的一点问题
  • Zhoujq  

  4
  4
  Posts
  25
  Views

  

  问题描述清楚的好处，除了“使人理解”（描述不清无法使人理解），还有“使更多的人理解”。有价值的问题值得转化成文本与世界讨论。所以下一步才是体现“描述清楚问题”价值的一步：问题翻译成英文，去cfd online找个合适的板块发帖
• F

  twoLiquidMixingFoam计算时，设置成laminar，就是DNS模拟吗？
  • fubianhanshu  

  3
  3
  Posts
  21
  Views

  F

  @bestucan 好的，谢谢
• 

  关于OpenFOAM里simple.loop()的咨询
  • win  

  3
  3
  Posts
  31
  Views

  大

  你好，请问这个问题你解决了吗？我也遇到类似问题。如果A程序是稳态问题，B程序是瞬态问题。而且AB还要循环多次求解。请问这个问题该如何解决呢？
• 

  如何利用setField设置非均匀分布的初场？
  • yhdthu  

  3
  3
  Posts
  38
  Views

  好

  你可以考虑swak4Foam这个插件，可以设置不均匀场
• 爱

  Hua-Klein parameter---一个涡量场参数的含义
  • 爱学习的LDH  

  2
  2
  Posts
  27
  Views

  爱

  找到了，是一个二维湍流下漩涡的精准判据。
• 

  OpenFOAM中kOmegaSST和yPlusRAS
  • qjh888  

  7
  7
  Posts
  3529
  Views

  青

  @sibo 现在问题解决了吗
• Z

  关于paraview中VOF的奇怪现象
  • Zhoujq  

  13
  13
  Posts
  66
  Views

  Z

  @东岳 哈哈，老师，不知道啥情况，本来就一台电脑有这种情况，昨天测试后两台电脑都出现了，最后连虚拟机都变成这样了~幸亏提前找到了病根~要不都不知道啥问题~ alpha.water和alpha.water.orig文件里的object只要都是object alpha.water;在paraview就只有一个alpha.water选项，它实际读取的是alpha.water.orig文件里的object alpha.water;所以显示的是单相； 把alpha.water.orig文件里的object改为object alpha_.water;，然后复制文件去掉orig变成alpha.water文件，这时里面的还是object alpha_.water;但是在终端执行setFields后，alpha.water文件里的object alpha_.water;自动变成了object alpha.water;，并没有人工修改~ 这时查看paraview，同时出现了alpha.water和alpha_.water；本来带有.orig的文件不该读取，不知道为什么读取了~这个bug太坑人了；还是要再次感谢 @bestucan ；
• 

  如何用extrudeMesh生成360度的axisymmetric 3D网格
  • yhdthu  

  1
  1
  Posts
  20
  Views

  

  如题，目前的方法是blockMesh先生成方形单层mesh，再利用extrudeMesh生成360的3D网格。 但是出现个问题，里面的面不会自动的成为interior，不太清楚怎么处理。 以下是extrudeMeshDict constructFrom patch; sourceCase "$FOAM_CASE"; sourcePatches (side1); exposedPatchName side2; extrudeModel sector; nLayers 10;//90; sectorCoeffs { point (0 0 0); axis (0 1 0); angle 20;//360; } flipNormals false; mergeFaces false;
• 沐

  采用焓-多孔介质方法模拟相变材料熔化过程交流讨论
  • 沐风  

  5
  5
  Posts
  34
  Views

  

  挺好挺好， buoyantPimpleFoam的解析那个我还没更新 https://www.cfd-china.com/topic/3665/关于buoyantpimplefoam的几个问题/12 rhoPimpleFoam目前更完善。 我这面带了一个学生目前打算做这方面工作，我还没找好具体的方向，不过也是类似的。熔化和凝固。我得先看看这方向是不是可以做的，然后再让他做