OpenFOAM

• 

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

  38
  38
  Posts
  46179
  Views

  X

  @东岳 谢谢东岳老师
• 一

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

  8
  8
  Posts
  84
  Views

  一

  @马乔 您好，我参照你的方法将transportProperties代码写进去，但是很奇怪的是，报错显示，request for dictionary transportPropertie ,而可获得的字典里是transportProperties,两者总是存在单复数的差别（就像你贴的代码里一样）。本来想着直接用令锑度边界条件绕开这个错误，但是最近发现还是得用alphat壁面函数。真是很奇怪呀。
• 七

  关于重叠网格尺寸的问题
  • 七辆战车  

  1
  1
  Posts
  8
  Views

  七

  请问有人知道重叠网格里嵌套网格和背景网格二者之间的尺寸关系需要满足什么条件吗？另外，如果运动物体有薄壳结构（厚度0.02m），那在划分背景网格时尺寸是不是一定要小于0.02m？恳请回复！
• 

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

  4
  4
  Posts
  23
  Views

  

  @zhe DPMFoam.C 这个的108行在算粒子，下面就开始用算过的信息更新alphac了。 108行的定义在这里 KinematicClud.C 673行。这里用到的solve函数就在这个文件的91行。 然后这solve里一大堆追着函数名找，大概就是particle distribution function的过程了。我猜的
• 液

  关于稳态欧拉-拉格朗日计算速度的探讨
  • 液固两相流  

  4
  4
  Posts
  31
  Views

  

  很有可能，cyclicAMI有可能把粒子算丢了
• 

  phiHbyA的计算
  • 东岳  

  1
  1
  Posts
  16
  Views

  

  ddtCorr这个是个老梗了，一直没时间深究，不过今天发现在某些算例下，区别很大 fvc::flux(HbyA) + fvc::interpolate(rho*rAU)*fvc::ddtCorr(U, Uf) 其中第二行加与不加的效果如下，左侧是不加，右侧是加
• 

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

  15
  15
  Posts
  77
  Views

  

  @bestucan 谢谢，已经说的很好了。我之前想过rho=(alpha1rho1+alpha2rho2)/(alpha1+alpha2).这个我觉得是比较合适的。虽然这样定义，直接就说明了alphac=alpha1+alpha2，而不是从后面的1-theta来的，但是我觉得这不会引起什么问题，因为宏观上来说，alphac就是等于alpha1+alpha2. 这样也满足了alpha1+alpha2+theta=alphac+theta=1. 都是归一后的了。我也试了自带的算例。计算出来的结果与之前定义rho=alpha1rho1+alpha2rho2有差别。因为没有真实实验数据可以对比，所以不敢说哪个就是对的，但是后改的这个，更能理解。单独定义rho后，我也就可以重新考虑alphac*rho这个项。我还是觉得对流项中是要考虑alphac的。希望，有朋友可以指出错误，谢谢！
• 七

  关于算例还在运行，但时间步不动了的问题
  • 七辆战车  

  3
  3
  Posts
  26
  Views

  七

  @bestucan 感谢回复！
• L

  结构化网格
  • leilei  

  8
  8
  Posts
  46
  Views

  L

  @bestucan 非常感谢，全明白了. 也谢谢Zhoujq和Samuel-Tu的回复
• 桎

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

  9
  9
  Posts
  50
  Views

  N

  @桎梏 不客气
• X

  关于网格高度变化的动网格设置
  • xiaolin  

  4
  4
  Posts
  20
  Views

  

  @xiaolin 那如果是预定变化轨迹的网格高度的动网格变化是不是应该参照movingConeTopoFvMesh呢 是的，这样简单的多
• 沐

  关于of-V中interCondensatingEvaporatingFoam处理密度疑问
  • 沐风  

  6
  6
  Posts
  23
  Views

  

  多谢老铁对我资料的认可
• K

  time step of twophaseeulerfoam
  • kimy  

  6
  6
  Posts
  28
  Views

  

  Your simulations just blow up, its difficult to debug the reason from your info.
• 一

  openfoam里双时间步如何设置
  • 一颗栗子  

  3
  3
  Posts
  631
  Views

  好

  hello，不知道你后来有设置这个双时间步吗？可能给看一下你的设置吗？谢谢！
• A

  求助 编译后出现“/usr/bin/ld：无法打开输出文件 /opt/openfoam4 ”如何处理
  • as_me  

  4
  4
  Posts
  16
  Views

  A

  @东岳 @bestucan 谢谢两位老师，我把自己打算编译的东西添加到$FOAM_USER_APPBIN路径下了。下次安装openfoam的时候再注意这个问题吧。
• 小

  带源项的相方程离散疑问
  • 小考拉  

  4
  4
  Posts
  62
  Views

  沐

  同学，你好。很高兴发现你也从事interPhaseChangeFoam求解器相关研究工作。我想问一下，在interPhaseChangeFoam求解器主程序中，在pimple.loop()每次循环中中为什么要声明创建surfaceScalarField rhoPhi场，这样做的目的是什么？
• 

  如何在paraview中消除入口段波动的影响
  • chengan.wang  

  1
  1
  Posts
  11
  Views

  

  大家好，如图所示，我在模拟流体通过水渠，水下有障碍物，因此在水面形成波，但水面波的强度较小，并且入口段形成一定的波动，波形不够明显，在paraview中是否可以进行处理，消除背景噪声使得结果更明显？
• 一

  流固耦合alphat使用alphatJayatilleWallFunction报错跪求大佬指点
  alphat边界条件 湍流热扩散系数 湍流自然对流 流固耦合 • • 一颗鸭蛋  

  2
  2
  Posts
  7
  Views

  一

  使用alphatJayatillWallFunction壁面函数的alphat结果。上面贴错了，不好意思。
• N

  openfoam-v2006
  • nanxuan  

  9
  9
  Posts
  39
  Views

  

  用 sudo fdisk -l 看看 swap 分区多大。如果很小或者没有，可能是因为内存。如果很大可以排除内存问题。swap是虚拟内存，物理内存不够用会把硬盘当内存用，只会很慢不会崩。 也可以在编译过程中开着 system monitor 看着内存和 swap 使用状况。 当初装linux的时候分区选择了lvm方式，还能调整swap分区大小。如果不是lvm而是直接分区。只能重装系统修改分区大小。 然后你的源代码放在OpenFOAM1，1??? 不知道会不会影响。
• 董

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

  10
  10
  Posts
  48
  Views

  

  @RebelYoung 这种错误可能是数据本身有问题，不一定是paraview的问题。 如果想确定，可以去paraview的官网下载最新的paraview。如果连它都崩，那就是数据的问题。
• X

  运行pisoFoam的时候，卡住不动了
  • xjwang  

  5
  5
  Posts
  39
  Views

  

  在求解器的各处随机插入 cout << "run in line xx" 重新编译这个求解器，看看具体在哪一步之后就不输出了。 根据描述，断点续算可以，但是直接突破不过去，而且只是不输出而不是不计算。应当是某个变量积累超过输出语句的判断范围。所以还在算但不输出。所以上面提到的具体哪一步就是修改了这个变量超出了范围。 更精准的定位错误就要看输出残差的函数怎么写的了。 OF 的 debug 模式编译有更强大和方便的调试功能。
• L

  求助，算例报错，好像是液滴跑到网格外？
  • ljj  

  4
  4
  Posts
  10
  Views

  L

  @yfclark
• 好

  新方程的预测步和修正步
  • 好好科研  

  3
  3
  Posts
  24
  Views

  好

  谢谢东岳老师
• 莫

  OpenFOAM获取叶片表面压力系数曲线
  • 莫晓柒  

  3
  3
  Posts
  149
  Views

  R

  请问实现了吗？
• 莫

  OpenFOAM 里计算Cp，Cf等参数
  • 莫晓柒  

  6
  6
  Posts
  189
  Views

  R

  @东岳 请问，这个是后处理实现吗？postProcess执行后说找不到这个函数。。
• 

  求助如何用paraview和gnuplot做出空翼的Cl，Cp，Cd分布图
  • yhdthu  

  2
  2
  Posts
  986
  Views

  R

  请问，后来。。。搞出来了吗？？
• E

  动态接触角的问题
  • ethanboem  

  28
  28
  Posts
  364
  Views

  P

  @Peng 大神能不能也分享一下fluent中动态接触角 的udf实现方法给我，十分感谢。penghui0410@foxmail.com
• S

  OF波浪模拟
  • Samuel-Tu  

  7
  7
  Posts
  56
  Views

  摇

  @Samuel-Tu 我以前的工作用过SWAN，SWAN一般用在大范围的波浪模拟上，适用于模拟一个海域的台风浪。它的方程是求解波浪谱，然后根据谱来求解波高、周期等参数。SWAN并不是基于过程的模型，无法给出一帧一帧的波面过程，算出来的结果仅仅是波高、周期一些参数随时间的变化。 shallowWaterFoam是基于浅水方程的模型，一般可以用来模拟像潮流这种周期很长的长波运动，基于静压假设后简化了的NS方程。我们做潮流的一般不用OF，用其他更成熟的浅水方程模型。
• 

  OpenFOAM小代码
  • 东岳  

  8
  8
  Posts
  248
  Views

  霜

  @Zhong-combustion 你好，我目前也在用LES做加热管流的湍流研究，关于湍流入口的设置看了一些资料，如有兴趣，可以加个qq:1191364394，共同学习探讨。
• 林

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

  8
  8
  Posts
  55
  Views

  林

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

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

  1
  1
  Posts
  13
  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
  518
  Views

  X

  @东岳
• 

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

  2
  2
  Posts
  13
  Views

  

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

  二维自适应网格
  • supersoldier  

  21
  21
  Posts
  9154
  Views

  

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

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

  7
  7
  Posts
  26
  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
  11
  Views

  小

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

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

  5
  5
  Posts
  82
  Views

  

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

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

  3
  3
  Posts
  26
  Views

  

  感谢分享！！
• 液

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

  3
  3
  Posts
  13
  Views

  液

  已成功植入emoij
• A

  interFoam的曲率计算
  • anubis  

  2
  2
  Posts
  14
  Views

  

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

  SHM 网格拉伸变形
  • aufalpen  

  1
  1
  Posts
  9
  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

  花

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

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

  7
  7
  Posts
  361
  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
  16
  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速度的意思呢？ 谢谢~