OpenFOAM

• 

  合集：OpenFOAM安装准备工作贴
  • 李东岳  

  1
  1
  Posts
  2185
  Views

  

  更换源 默认的源下载什么都非常慢，请更改国内的阿里云源，安装openfoam依赖包5秒完事。 在国内下载OpenFOAM 国外的源也是特别慢，直接在此处下载国内源 算力测试 装好OpenFOAM之后，可以在这里测试下你的算力 切换gcc版本安装OpenFOAM Ubuntu 20.04 安装 OpenFOAM-v1712 杂七杂八 新版ubuntu安装老版本OpenFOAM waves2foam 安装 OpenFOAM常用小代码
• 

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

  39
  39
  Posts
  55506
  Views

  

  可能还是得有需求才能推进学习吧，一直想搞懂OpenFOAM机制上的优势，但一直也还没到需要自己改代码的程度，混，混就完了
• M

  fvSolution问题
  • mark  

  5
  5
  Posts
  1510
  Views

  

  @wwzhao 好的我明白了，您说的外循环指的是pimple算法的循环，内循环指的是线性方程组的迭代求解，pimple循环中每次都要求解一次压力方程，一共求解2-4次左右，是吗？谢谢老师！
• 

  老师们，为什么稳态问题也有显隐式的区分呢？
  • Micro  

  1
  1
  Posts
  5
  Views

  

  众所周知，如上是稳态求解器simpleFoam中的控制方程的描述。可以看到，作为稳态求解器的simpleFoam并没有出现时间项。那么问题来了，显隐式的区别不正是针对时间项的吗？对于未知量如果代入上一时间步的已知值，就是显式，如果视为本时间步的未知量，就是隐式。那么既然如此，simpleFoam中连时间项都没有，为什么也会出现fvm，fvc这样的显隐式的用法呢？谢谢老师@大佬
• B

  使用RapidCFD运行OpenFOAM算例时是否需要对输入文件做相应修改？
  • baojd42  

  11
  11
  Posts
  2629
  Views

  

  @xpqiu 好的谢谢，我按您说的着重看看这两个文件夹
• L

  请教个简单问题
  • leilei  

  6
  6
  Posts
  52
  Views

  L

  @xpqiu : 谢谢老师您的精彩点评，师之教诲，学生记下来 谢谢老师 @veen : 也谢谢veen
• C

  关于postProcess得到张量的排序
  • chszkc  

  2
  2
  Posts
  23
  Views

  X

  @chszkc 在 关于postProcess得到张量的排序 中说： 想请问一下各位OFer，通过命令 postProcess -func "grad(U)" 得到的速度梯度张量共有九个元素，打开数据文件后发现这九个元素是以一维数列的形式保存的，那这九个元素的排列顺序是什么呢？Google了大半天都没找出来..... 之前自己默认顺序是dudx,dudy,dudz,dvdx,dvdy,dvdz,dwdx,dwdy,dwdz，但算出来的结果有点怪，所以怀疑自己是不是在排列顺序上理解错了 应该不是这个顺序，而是转置过来这个顺序
• Y

  turbulence->divdDevReff(U)
  • yuhx  

  6
  6
  Posts
  945
  Views

  

  @xpqiu gzt200361@163.com 非常感谢
• 

  pisoFoam中的UEqn.relax();也进行方程松弛？
  • aiweimo  

  2
  2
  Posts
  12
  Views

  X

  @aiweimo UEqn.relax() ，即使松弛因子是1，其实也会对 UEqn 的 系数矩阵进行一些调整，以使它更满足对角占优。fvSolution 里面如果没有定义松弛因子，则默认值为1。
• 中

  DES边界条件
  • 中国轰炸机飞机斯曼  

  7
  7
  Posts
  3670
  Views

  X

  @aiweimo k-w方程是低雷诺数方程，本身不需要壁面函数 严格来说，应该是 $k-\omega$ 方程可以一直算到 $y^+<1$，但是，$k-\omega$ 方程也可以结合壁函数来算$y^+>30$ 的网格。相比之下，标准$k-\varepsilon$模型的$\varepsilon$ 方程则只在log-law region之外才严格有效，不可直接用来算很小 $y^+$ 的网格。 OF 里面的那些湍流壁面边界条件，有一点容易引起误解的是他所有的名字都带 wallFunction，而大部分人看到 wall function，想到的都是标准壁函数，也就是 log-law region的 $U^+=\frac{1}{\kappa}\ln(Ey^+)$ 以及层流底层的 $U^+=y^+$ 。 实际情况下原比这个复杂，因为很难保证整个region 都满足 $y^+ > 30$ 或者 $y^+<1$，所以，对于实际问题，湍流模型的壁面处理需要考虑如何处理$1<y^+<30$ 的区域。文献中有很多相关的研究，openfoam 的 omegaWallFunction 其实就是一种，里面构造了一个blending function 把适合 $y^+ > 30$和 $y^+<1$的 omega 公式blend一下，来让omega 在任意 $y^+$ 下都能算出来一个相对合理的值。 所以，OF 里面那些 wallFunction 的边界，可以理解为定义多种不同的 wallTreatment 方式，具体怎么用，需要结合着实际情况来。2 楼提到的就是一种常见的用法。文献中壁面处理的方式太多太多了。。。 至于 $\omega_{wall}$ 指定为零梯度或者固定值，这种做法不太合理。因为 $\omega$ 在壁面附近需要满足某种渐近特性（asymptotic），零梯度或者固定值都不满足这个条件，虽然可以可以算出来结果，但是不太符合物理。
• D

  OpenFOAM中MRF的代码
  • dzw05  

  3
  3
  Posts
  1833
  Views

  

  我最近看了下， 上面这个方程应该是稳态算法的方程。里面MRF.DDt(U)用于处理$\Omega\times\bfU_I$。 但是上面这个方程不是OpenFOAM植入的方程。在瞬态情况下应该是下面这个： $$ \frac{\p \bfU_I}{\p t}+\nabla\cdot(\bfU_I\bfU_R)+\Omega\times\bfU_I=... $$ 下面的瞬态代码分别对应： tmp<fvVectorMatrix> tUEqn ( fvm::ddt(U) //第1项 + fvm::div(phi, U) //第2项 + MRF.DDt(U) //第3项 + turbulence->divDevReff(U) == fvOptions(U) );
• 

  当我们增加自定义库的时候，为什么不需要重新编译上层应用？
  • Micro  

  1
  1
  Posts
  10
  Views

  

  众所周知，当我们编译好一个OF中新的自定义库libxxx.so的时候，我们只需要做的就是在controlDict文件夹中增加libs("libxxx.so")的一句话，此时我们就可以使用该库了。我想问的是，此处可以这样做的深层机制是什么？为什么我们可以通过增加libs("libxxx.so"）的方式来避免重新编译上层应用呢？谁能解释下OF在此处用到的"黑魔法"呢？谢谢诸位！
• 

  请问高阶格式的矩阵还能用LDU格式存储吗？
  • Micro  

  1
  1
  Posts
  8
  Views

  

  众所周知，OF中对于矩阵存储采用的是LDU格式，可是这种格式有一种限制，就是只能存储对称（元素值不一定相等，但是非零元素位置对称）矩阵。那么这就产生了一个问题，对于某些采用高阶迎风格式离散后的矩阵,必然不再是对称矩阵，OF能不能存储呢？或者说OF是以什么样的格式来存储的呢？还是说OF根本不支持高阶迎风呢？谢谢各位大佬!
• S

  local index to global index
  • Samuel-Tu  

  4
  4
  Posts
  124
  Views

  

  @Samuel-Tu 在 local index to global index 中说： 但是写入solver的时候调用不了这个函数。说没有定义此函数，如何调用这个函数呢。 大佬您好，请问您解决这个问题了吗？我也想问一下怎么调用这个函数，感恩！
• 

  OpenFOAM小代码
  • 李东岳  

  20
  20
  Posts
  2126
  Views

  H

  @李东岳 FoamFile { version 2.0; format ascii; class scalarField; location "constant/../postProcessing"; object totalLiquid; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 1(0.2) 麻烦老师再帮我看一下，为什么每次都只会输出一个时刻的值，不会随时间增加输出的
• S

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

  3
  3
  Posts
  58
  Views

  

  大佬您好，为啥 那请问局部索引应该怎么显示呢？谢谢！
• 李

  如何在paraview中显示颗粒浓度
  • 李梦瑶  

  1
  1
  Posts
  7
  Views

  李

  请问各位老师，如何在paraview中显示颗粒浓度kinematicCloud.theta 做颗粒后处理时，只有alpha.air,如图所示，
• R

  并行效率疑问
  • random_ran  

  62
  62
  Posts
  21446
  Views

  

  @浪迹天大 好的，我自己看下，有进展再向您汇报
• S

  functionObject里面可以加自定义的程序吗
  • Samuel-Tu  

  9
  9
  Posts
  255
  Views

  H

  @Samuel-Tu 您好，您例子中输出的是一个场量，如果我只是想把一个数值输出到postProcessing中，应该定义一个什么样的类型（类似IOField）。 const volScalarField& b = mesh().lookupObject<volScalarField>("alpha.liquid"); scalar c= b.weightedAverage(mesh().V()).value();
• 

  关于icoFoam解析的一个小问题
  • aiweimo  

  3
  3
  Posts
  27
  Views

  

  多谢！已更新
• 

  window下的paraview如何读取拉格朗日场？
  • 孟豪龙  

  3
  3
  Posts
  19
  Views

  

  @李东岳 谢谢李老师的回复
• D

  关于icoFOAM中通量phi的更新
  icofoam • • dabiu  

  3
  3
  Posts
  78
  Views

  D

  @李东岳 感谢李老师的回复，谢谢
• Z

  IBM浸没边界法生成拉格朗日标记点
  • zp气固两相流  

  8
  8
  Posts
  1315
  Views

  O

  @zp气固两相流 在 IBM浸没边界法生成拉格朗日标记点 中说： A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries 好的谢谢，我去看看这两篇文章
• H

  动网格细化植入interFoam
  • Hunahum  

  6
  6
  Posts
  198
  Views

  F

  @李东岳 好的，谢谢东岳老师！
• V

  运行时报错
  • Vortex  

  2
  2
  Posts
  54
  Views

  

  这个就是发散了。具体原因需要更深层次的debug
• U

  怎么读取openFoam中系数矩阵
  • upc_ngh  

  3
  3
  Posts
  54
  Views

  U

  @李东岳 岳 东岳老师好，上面公式S其实表达的意思不一样的，我这里用错了，现在存在的困难是怎么求解这个压力方程值p'，感觉好难写出来只关于p'的方程，再进行求解，希望能得到东岳老师和各位前辈的指点，谢谢！
• A

  interFoam + 非正交网格计算
  • anubis  

  2
  2
  Posts
  64
  Views

  

  你这个最好用wedge网格来计算就不会出现这个问题了
• V

  可以直接提取energy dissipation rate吗？
  • Voynich  

  2
  2
  Posts
  30
  Views

  

  要自己算，还挺麻烦的 https://www.cfd-china.com/topic/4320/求助-openfoam中的les-dns如何输出湍流耗散率/2?_=1618104636253
• H

  多组分模型multiComponentPhaseModel，组分的分数分层？
  • Henry Zhou  

  1
  1
  Posts
  15
  Views

  H

  大家好，我用multiphaseEulerFoam做一个案例，对某一相使用多组分模型multiComponentPhaseModel，但是相中组分分数分层。因为组分会扩散，正常情况是不应该分层的，调了一段时间，只是换了分层的位置。不知道大家遇到过这种情况不？谢谢各位老师的指点！
• 

  Ubuntu 20.04 安装 OpenFOAM-v1712
  • 李东岳  

  3
  3
  Posts
  829
  Views

  

  李老师您好， 我按照您的提示： sudo update-alternatives --install /usr/bin/gcc gcc /usr/bin/gcc-5 100 目前默认gcc-5 gcc (Ubuntu 5.3.1-14ubuntu2) 5.3.1 20160413 Copyright (C) 2015 Free Software Foundation, Inc. This is free software; see the source for copying conditions. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. 我想切换回g9，用命令 sudo update-alternatives --config gcc 得到： There is only one alternative in link group gcc (providing /usr/bin/gcc): /usr/bin/gcc-5 Nothing to configure. 没有切换提示，咋办？ 我用命令 ls /usr/bin/gcc* /usr/bin/gcc /usr/bin/gcc-ar /usr/bin/gcc-nm /usr/bin/gcc-ranlib /usr/bin/gcc-5 /usr/bin/gcc-ar-5 /usr/bin/gcc-nm-5 /usr/bin/gcc-ranlib-5 /usr/bin/gcc-9 /usr/bin/gcc-ar-9 /usr/bin/gcc-nm-9 /usr/bin/gcc-ranlib-9 应该是g9也在的，不知道为啥不能切换，跟百度不一样
• 冠

  在编译求解器时引用弹性项
  • 冠竹  

  5
  5
  Posts
  77
  Views

  冠

  @sungda 十分感谢您
• 

  矩阵求解器
  • 肖恩曼  

  7
  7
  Posts
  1036
  Views

  

  @王慧博 在 矩阵求解器 中说： LUsolve 您好。请问LUsolve是OF自带的函数吗？ 还有一个问题，为啥LUsolve可以直接调用而不是通过对象访问成员函数的方式调用呢？难道LUsolve是一个全局函数吗？谢谢大神！！！
• 搬

  OpenFOAM中的residual？
  • 搬运工不好当  

  4
  4
  Posts
  1501
  Views

  

  主要输OF中的残差是线性方程组求解的残差，不是时间步迭代的残差，您说对吗？李老师。另外，是不是只有稳态问题我们才会关注线性方程组求解的残差残差，瞬态问题的residual 有意义吗？
• F

  OpenFoam湍流模型修改增加条件语句，运行报错
  • fangyuanaza  

  3
  3
  Posts
  46
  Views

  F

  已解决。。。先给出ADD的表达式，然后对其进行限制。修改代码为：
• W

  TopoSet和RefineMesh局部加密网格
  • WYing  

  1
  1
  Posts
  31
  Views

  W

  使用topoSet和refineMesh进行了两次网格加密，但可视化之后发现在topoSet没有指定的顶部区域也进行了加密，检查blockMesh，topoSet和refineMesh的过程都没有问题，网格区域选取也正确。唯一能想到的原因是top和bottom的边界条件是cyclic，而第一次加密的区域包含了最底层网格，所以根据cyclic条件，第一次加密也自动包括了最顶层网格？
• E

  大涡模拟壁湍流雷诺剪切应力
  • ewz  

  6
  6
  Posts
  906
  Views

  

  @李东岳 WALE不需要vanDriest吧？
• 搬

  如何显示涡量图
  • 搬运工不好当  

  12
  12
  Posts
  3370
  Views

  如

  @huntersong 哇，好详细，我去试试， 十分感谢！
• 明

  关于reactingFoam中reactionsGRI的问题
  • 明小米  

  7
  7
  Posts
  55
  Views

  明

  @wangfei9088 哇，太谢谢了。我也考虑过是单位的问题，但还未在代码中找到数据验证猜想，非常感谢大佬！
• Z

  在曳力模型中植入网格中心到边界距离的问题
  • Zhy2022  

  3
  3
  Posts
  38
  Views

  Z

  在dargmodel.C这个衍生类中没有mesh这个成员对象，不知道该如何引入
• 

  paraFoam是否可以生成旋转叶片的表面流线和周围流线？
  • tidedrinker  

  12
  12
  Posts
  142
  Views

  C

  @xpqiu 好的，谢谢您，我去查看测试一下
• J

  关于理解setFieldsDict
  • Joann  

  5
  5
  Posts
  105
  Views

  J

  @xpqiu 感谢您的答疑 原来如此，是两个box相当于框定了一个水域区范围间的意思。 zoneID场的值设置明白了，还想问一下zoneID这个场是不是因为重叠网格所以要设置的，具体作用是什么呢？
• 

  H2/CO燃烧
  • liujunCFD  

  3
  3
  Posts
  55
  Views

  

  我大部分情况下模拟结果都是这种 上面呢种是因为激波前移叠加了？ 不太懂啊
• 

  MPPIC中，怎么实现的alpha-alphaPacked?
  • zhe  

  1
  1
  Posts
  30
  Views

  

  如题，最近一直想不明白，在MPPIC中，怎么可以进行alpha和alphaPacked做差的。一个是fieldScalar，一个是scalar.不知道是不是哪里漏掉了，希望有人帮忙解惑，谢谢。
• W

  Doxygen中各模块关系图看不懂？
  • warnerchang  

  7
  7
  Posts
  78
  Views

  W

  @bestucan 好的，谢谢老师的耐心解答！
• Y

  fsiFoam安装
  • Yan  

  8
  8
  Posts
  350
  Views

  L

  @sr0922msc 我没有尝试过运行4.0之前的版本，即使是4.0版本，我在运行tutorials案例的时候，也遇到了一些问题，有些可以运行，有些不可以
• 好

  计算添加的方程想用更小的时间步，为什么不能在C文件下改动
  • 好好科研  

  2
  2
  Posts
  59
  Views

  

  因为这不是C 这是C++。 你在这里修改，但是其他过程在其他地方修改了。而且这个赋值八成也是不成功的。给类成员赋值一般都由成员函数完成。 改起来还是挺麻烦的。
• 

  openfoam中cyclic周期性边界的问题
  • 余正东  

  22
  22
  Posts
  8172
  Views

  T

  @dzw05 最近在研究processorFvPatch，来挖个坟补充下~~ 对于并行耦合边界条件，在方程组迭代求解中的“耦合”过程在矩阵-向量乘法（Amul）中的实现如下： 1、initMatrixInterfaces：只负责初始化边界数据，将耦合边界内侧（本地）的数据传送到外侧（相邻进程）。 2、本地数据与向量相乘，参考程迪的描述。 3、updateMatrixInterfaces：使用第一步初始化得到的数据来耦合外侧变量（$C\cdot x_n $）. 并行边界条件中，对于一个本地进程来说耦合过程是单向的，只将外侧数据耦合到本地，但是对于一个耦合边界面来说耦合过程是双向的，即面两侧的进程都需要耦合外侧数据（$C=C_L+C_U$ 一侧是$C_L$另一侧是$C_U$）。 耦合边界类使用的接口是initInterfaceMatrixUpdate，updateInterfaceMatrix，这两个接口的实现在对应的子类边界条件中，如何想研究何以找找。
• J

  怎么设置使snappyHexMesh的网格颜色和blockMesh的网格颜色不一样
  • Joann  

  2
  2
  Posts
  44
  Views

  

  显示变量那里，选vtkblockcolors或者vtkcompositeindex都行，再设定颜色
• 星

  编译遇到的问题
  • 星星星星晴  

  2
  2
  Posts
  48
  Views

  

  掩耳盗铃法：mkdir那些不存在的文件 偷梁换柱法：把of4里这里的代码复制过来 根治：修改代码中调用of的接口，不要调用这些已经不存在或者改名的模块
• 

  paraview输出尾涡
  • chengan.wang  

  3
  3
  Posts
  931
  Views

  如

  @李东岳 老师好，我想问一下，为什么我用tecplot的Q准则算出来后，都是一个颜色。