CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新
    1. 主页
    2. OpenFOAM
    登录以发表
    • 从新到旧
    • 从旧到新
    • 回复最多
    • 最多赞同
    • Most Views
    • 王

      求UnnamedMoose / BasicOpenFOAMProgrammingTutorials 编程交流
      • 王杰

      1
      1
      帖子
      40
      浏览

      王

      题目就是想要问的,有没有大佬也在做这,求交流啊 。
      我C++可以说基础等于无,囫囵吞枣看了一本陈志泊的《面向对象的程序设计语言》,明白基本原理,高级的技术都还不会呢。所以做这个代码就等于抄写。要是有也在做这个的,且恰巧没有人交流的,交流交流哈。

    • 五好青年

      请教关于浸入边界法的问题,IBM方法的问题
      • 五好青年

      35
      35
      帖子
      2000
      浏览

      E

      @李东岳 delta函数

    • 风大仙

      湍流模型编译通过但是在计算的时候显示不能加载动态库
      • 风大仙

      3
      3
      帖子
      32
      浏览

      风大仙

      @李东岳 谢谢李老师,计算之前controlDict已经事先挂载过了

      libs ("test1_libIncompressibleTurbulenceModels.so");
    • Z

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

      11
      11
      帖子
      2539
      浏览

      Z

      @eason 如果是multi-body 问题,然后其中一个body需要IBM,就直接用Gmsh做其他body 的conforming mesh就可以。只有一个body且需要IBM,就可以做一个enclosed的正方体或者圆柱体,画网格就行。

    • T

      拖曳力、升力
      • ThomasShi

      11
      11
      帖子
      1213
      浏览

      J

      @xpqiu 老师您好,我有一个关于监测多孔介质区域体积力的问题想请教下您,我在波浪场中前后设置两块多孔介质区域porosity1和porosity2,然后用function函数分别输出porosity1和porosity2的受到的体积力,但是两个输出的力是一模一样的,经过计算我发现输出的力是porosity1和porosity2的合力,也就是说openfoam没有办法分别输出porosity1和porosity2各自的受力?谢谢老师!

    • C

      关于ubuntu18.04下paraFoam黑屏的问题
      • cyberk

      4
      4
      帖子
      38
      浏览

      李东岳

      感谢分享。

      在没有独显的服务器上,windows的paraview也是黑屏。

    • 四

      求助:如何在边界a上索引边界b的值,且可以并行(包含多种边界值获取方式)
      • 四季之中的隐星

      3
      3
      帖子
      55
      浏览

      四

      @李东岳 谢谢李老师的回复,这个问题我最后发现是一个很简单的错误导致的。因为decompsePar没有把边界场进行拆分:136: ,导致各个处理器文件夹下的场定义缺失所以发生了错误

    • M

      OF模拟结果涡脱频率是正常结果的四分之一
      • Mania

      3
      3
      帖子
      61
      浏览

      drink water

      @Mania 怀疑fft变换时参数选错了,OF不大可能错,因为你说之前的都没问题。

    • R

      OpenFoam边界条件设置
      • Rachel0096

      5
      5
      帖子
      160
      浏览

      李东岳

      @Chenwei 这个确实是个问题。主要就是这个边界条件的组合会影响计算的结果。你有尝试过全部totolPressure么。或者进口速度给一个稍微小的值,如1e-3,然后流量吧这个减掉?或者干脆不剪掉,差得多么。

      原始算例应该考虑的是稳定性的问题。

    • T

      Wilson面
      • Turbulence

      1
      1
      帖子
      20
      浏览

      T

      请教各位老师,Wilson面的概念和物理意义以及表达形式,搜索到下面的论文可是下载不下来,谢谢:xinxin:
      Bohn D , Kerpicci H , Ren J , et al. Homogeneous and heterogeneous nucleation and wet-steam Flow in a nozzle and a LP-turbine cascade. 2002.

    • L

      OPENFOAM中有些后处理工具是否只能在计算时运行,计算完毕后无法运行?
      • lllwonderliquid

      7
      7
      帖子
      2258
      浏览

      S

      @李东岳 我这试了下system下创建fieldAverage文件并复制上述代码,删掉functions和fieldAverage1可以计算完成后处理,controDict不需要挂载,完整的文件内容如下

      /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 9 \\/ M anipulation | ------------------------------------------------------------------------------- Description Calculates and writes the time averages of given list of fields. \*---------------------------------------------------------------------------*/ type fieldAverage; libs ( "libfieldFunctionObjects.so" ); writeControl writeTime; timeStart 0; timeEnd 100; resetOnOutput false; fields ( U { mean on; prime2Mean on; base time; } );
    • ZY-GONG

      液态水的EOS
      • ZY-GONG

      3
      3
      帖子
      38
      浏览

      李东岳

      :135: :135:

    • 洱聿

      在fe40上安装waves2foam环境变量出错
      • 洱聿

      2
      2
      帖子
      32
      浏览

      洱聿

      简单记录一下,两个问题
      1.ocean-wave3d这个安装包没有下载好
      2.fe40需要gcc5版本,虚拟机里没有统一成一个版本。

    • 学流体的小明

      气泡槽道流的计算,求指点
      • 学流体的小明

      2
      2
      帖子
      63
      浏览

      李东岳

      气泡群加入到槽道中,有周期性边界,让它们一直输运,我想要考察壁面上的脉动压力。

      只有特别微观的槽道,10几个气泡的能算。这个我之前算过,算半个月也没算出来个统计平均值。后来就不整了。目前也没遇到过什么好办法。VOF算这种千万网格打底。除非换欧拉方法,特别快。

      看看其他大佬有啥意见

      Tryggvason那面算了非常多的这种东西,我估计你肯定知道。

    • L

      网格加密后的奇怪现象
      • leilei

      3
      3
      帖子
      103
      浏览

      L

      @bestucan
      先谢谢老师的回复
      老师说的对。。。。

      谢谢老师

    • R

      rhoCentralFoam发散,出现Maximum number of iterations exceeded
      • Rachel0096

      6
      6
      帖子
      68
      浏览

      R

      应该是我的边界条件设置的不合理,我修改了物性和边界条件以后,目前是能继续算下去的。
      贴一下我参考的网址:
      https://www.researchgate.net/post/Can-anyone-see-this-problem-in-openfoamcan-anyone-solve-it/1

    • R

      浸没边界法+超音速自定义求解器出现“浮点数例外,核心已转储”
      • Rachel0096

      12
      12
      帖子
      194
      浏览

      李东岳

      好的好的

      我发现这个欧拉场到拉格朗日场,拉格朗日场到欧拉场的映射,有各种不同的方法..

    • strive_DUT

      请问LES模拟后平均值与真实值不对应是什么原因导致?
      • strive_DUT

      1
      1
      帖子
      35
      浏览

      strive_DUT

      如题,请教各位老师同学,我在LES模拟后,对数据进行后处理时,发现一个问题:定压比热的场值经过平均后被缩小了,但是其他场的平均值没有这个问题,请问各位这是怎么回事呢?平均后的峰值(8000)明显低于真实值(120000),如图所示:
      比热.jpg
      比热云图.png

    • vbcwl

      槽道流添加源项时不考虑流场内部的能量耗散?
      • vbcwl

      5
      5
      帖子
      31
      浏览

      vbcwl

      @李东岳 only 外力,就是被这个误导了才没转过来弯:136:

    • V

      rhoCentralFoam求解器的问题
      • vivian

      10
      10
      帖子
      7546
      浏览

      R

      @vivian 我也遇到了一样的问题,请问您的问题解决了吗,可以指点一下吗?

    • L

      HPC 上编译 waves2foam 的问题
      • leilei

      3
      3
      帖子
      29
      浏览

      L

      大家好,
      还有如下问题
      /lib64/libstdc++.so.6: version `GLIBCXX_3.4.29' not found
      怎么解决?

      谢谢

    • J

      openfoam 如何在双cpu工作站并行计算
      • Jasper 0

      3
      3
      帖子
      69
      浏览

      J

      @李东岳 谢谢李老师答复。我现在就在用-oversubscribe方法跑,但是我不知道这个方法具体是怎么调用的核数,比如是不是我的两个cpu都调用了,还是说只调用了一个cpu里的38个线程(一个cpu总共48线程)

    • 同学博

      blockMesh生成复杂模型
      • 同学博

      12
      12
      帖子
      1758
      浏览

      李东岳

      三年了。继续膜拜blockMesh大佬

      :haqi:

    • 搬

      有人在mac上装of么?
      • 搬运工不好当

      30
      30
      帖子
      4715
      浏览

      M

      @田畔的风 在 有人在mac上装of么? 中说:

      @一颗鸭蛋 默认情况下,你编译求解器生成的可执行文件会保存到$FOAM_APPBIN,比如我的默认在<上级文件夹>/OpenFOAM-v2206/platforms/darwinARM64ClangDPInt32Opt/bin,但是这个App的加载磁盘是只读的,所以会报错。

      解决方案有两个:

      修改在Make/files中你的求解器保存路径 EXE = $(FOAM_APPBIN)/xxxFoam -> EXE = <新路径>/xxxFoam

      然后在环境变量中为$FOAM_APPBIN附加这个路径,以便让OpenFOAM在执行时能搜索到这个文件

      export $FOAM_APPBIN=<新路径>:$FOAM_APPBIN 创建一个可读写的分区,直接编译安装。参考 https://github.com/BrushXue/OpenFOAM-AppleM1。

      您的意思是把openfoam安装到创建的可读写分区吗?我用的也是M1的Mac,现在想在Mac上编译自己修改的求解器,大佬能不能给点详细的指导?谢谢🙏

    • 李东岳

      中性大气环境湍流动能的自保持 | 附有算例下载
      • 李东岳

      27
      27
      帖子
      834
      浏览

      C

      @李东岳 或许是我这边fluent有某个设置错了,杨老师论文中顶面用的slip结果不会出现k衰减。我再试了fluent中顶面设置no-slip,结果是比较符合预期,接近顶面壁面平均速度减小,湍动能增大
      Ux剖面.png TKE剖面.png

    • 卡

      waves2Foam编译问题
      • 卡洛

      9
      9
      帖子
      922
      浏览

      卡

      @GenhouSUN 两台电脑,一个是9.5.0,一个是9.4.0。而且我别的版本安装不上,只在OpenFOAM-v1912版本上用,waves2Foam是2137版本

    • 李东岳

      添加动量方程源项导致速度持续增加?
      • 李东岳

      6
      6
      帖子
      126
      浏览

      李东岳

      @Vortex 大佬太厉害了。用语言描述就是,地表速度低,科氏力弱,压力梯度导致风向偏移;高度上升后,气流速度增加,科氏力变强,将偏移拉了回来,并最终与压力梯度产生的压力平衡。

    • 就

      OpenFOAM中怎么才能在LES模拟时获取湍动能分布?
      • 就是个胖子

      7
      7
      帖子
      2107
      浏览

      疏影横斜水清浅

      @李东岳 好的,明白了,谢谢李老师

    • A

      fvOptions中添加的源项需要考虑网格的体积吗?
      • allanZhong

      2
      2
      帖子
      282
      浏览

      李东岳

      @allanZhong 2年过去了。。。是需要乘以体积。例如我刚写了个科氏力源项:

      nergySource { type coded; selectionMode all; field U; codeInclude #{ #}; codeAddSup #{ Pout<< "**codeAddSup**" << endl; const volVectorField& U = eqn.psi(); vectorField& USource = eqn.source(); vector omega(0, 0, 5e-5); const DimensionedField<scalar, volMesh>& V = mesh().V(); USource += (2.0*omega)^U*V; #}; codeAddRhoSup #{ Pout<< "**codeAddRhoSup**" << endl; #}; codeAddAlphaRhoSup #{ Pout<< "**codeAddAlphaRhoSup**" << endl; #}; }
    • Y

      paraview如何绘制二维流线
      • yfclark

      7
      7
      帖子
      1573
      浏览

      知

      这个方法很好,不过请注意,这个链接中的SurfaceVectors 应该是作为Seed Source 而maskPoints作为Input才对,作者笔误了.亲测可用.原文链接如下 https://zhuanlan.zhihu.com/p/592092340
      还有一个带有箭头的流线图链接如下,不过我并没有试过哈,也贴出来.
      https://blog.csdn.net/weixin_44571596/article/details/122187429

    • C

      气液搅拌釜顶部出口速度异常,导致epsilon等结果发散
      • Caijinjin

      5
      5
      帖子
      50
      浏览

      C

      @李东岳 好的,我先尝试着修改下再去跑跑看

    • 1

      求教,CENTOS无网络安装OPENFOAM出现问题
      • 18126359560

      9
      9
      帖子
      50
      浏览

      1

      @李东岳 可以安装了

    • C

      OpenFOAM中壁面函数问题
      • coolhhh

      5
      5
      帖子
      113
      浏览

      C

      @李东岳 李老师,这是我的邮箱 coolclw67@163.com。之前没研究RANS,但挺感兴趣的,最近我也研究下,谢谢李老师

    • vbcwl

      某个面开始是封闭的,当压力变大到一定值后变为流入流出边界
      • vbcwl

      3
      3
      帖子
      35
      浏览

      vbcwl

      @李东岳 谢谢东岳老师,看了一下完美满足需求:146:

    • 流

      sprayFoam中的液滴蒸发速率
      • 流体中二少年

      1
      1
      帖子
      39
      浏览

      流

      本人用sprayFaom进行一些模拟,请问如何获得液滴的蒸发速率?
      质量方程:1676896511168.jpg
      代码:

      solve ( fvm::ddt(rho) + fvc::div(phi) - parcels.Srho(rho) );

      论文中说Smass是蒸发速率,那我需要把代码中的parcels.Srho(rho)给输出出来,如何实现呢?

    • 香柏树

      sprayFoam计算等压环境下点火延迟时间
      • 香柏树

      8
      8
      帖子
      287
      浏览

      流

      前辈,你好!请问sprayFoam中的液滴蒸发率怎么获得?

    • J

      传热撞击流模拟
      • jinlinna

      4
      4
      帖子
      735
      浏览

      李东岳

      @wpx http://dyfluid.com/buoyantSimpleFoam.html

      @jinlinna 我明天看看

    • C

      计算y+,Δx+,Δz+的疑问
      • coolhhh

      13
      13
      帖子
      392
      浏览

      李东岳

      @coolhhh 在 计算y+,Δx+,Δz+的疑问 中说:

      但最后还是回归到网格除了法向外,另外两个方向怎么定义x和z方向

      是的。一般sci里面怎么弄?我个人怀疑sci就是按照流向方向来了。

    • H

      求助:openfoam三维算例中如何沿着某一个方向做空间上的平均
      • hongjiewang

      3
      3
      帖子
      54
      浏览

      H

      @学流体的小明 在 求助:openfoam三维算例中如何沿着某一个方向做空间上的平均 中说:

      我算槽道流,是有一个专门的postChannel后处理程序的,可能你得自己编程做了。
      可以参考一下这个后处理程序。
      https://develop.openfoam.com/Development/openfoam/-/tree/OpenFOAM-v2012/applications/utilities/postProcessing/miscellaneous/postChannel

      感谢,我发现了一个叫fluidfoam的后处理python工具,已经植入到openfoam中,需要py3.8以上版本,of的版本估计也得是2106以上,具体可以参照他的官方说明,该工具功能极其强大且便携

    • Y

      LES计算时初边值条件中的B和s是什么
      • Yzhang

      4
      4
      帖子
      55
      浏览

      Y

      @coolhhh 这里的$S_k$指的是k输运方程的源项,应该不是0文件夹中的s

    • L

      openfoam里面的煤燃烧intrinsicRate model方程植入单位不统一
      • Leeseungnan

      2
      2
      帖子
      28
      浏览

      李东岳

      @Leeseungnan 在 openfoam里面的煤燃烧intrinsicRate model方程植入单位不统一 中说:

      我发现两个变量的单位不一致,但是在方程植入的时候,程序里面却将两项做的差值处理,这种操作可以吗?

      理论上是不可以的。但是还见过有类似的操作。我也很懵。举个例子,下面这个方程,ln()里面不能带单位,但是里面确是有单位的。类似下面这个公式在很多sci里面都有,我只是举例。搞得我很懵。

      捕获.PNG

      类似你那种两个变量相减,a-b这种,就更需要单位一致了。这个问题我认为你可以问一下openfoam.com那面:1)要么是他们单位写错了,2)要么是有问题

    • 卡

      滑移网格多轴旋转
      • 卡洛

      5
      5
      帖子
      206
      浏览

      冷

      @卡洛 前辈方便的话能加下您的微信吗?想向您具体请教下,谢谢。我的微信是13652081916

    • tidedrinker

      interFoam短周期波浪衰减过于严重
      • tidedrinker

      10
      10
      帖子
      199
      浏览

      L

      @tidedrinker
      我没验,你验证下,看到底用哪个理论。。。 深h对T影响有限,我更倾向深水用stokes1

      还有,请教下:
      你熟悉不规则波浪反射的估算吗?

    • 杨

      可凝气体的冷凝数值模拟
      • 杨英狄

      2
      2
      帖子
      63
      浏览

      李东岳

      @杨英狄 在 可凝气体的冷凝数值模拟 中说:

      高温水蒸气在液态水中的冷凝及流动

      不太好理解。水蒸气在水里冷凝然后流动?

    • C

      关于多相流的电势方程(Laplace方程)的求解问题
      • cfd_lilili

      13
      13
      帖子
      1615
      浏览

      李东岳

      后来这个解决了吗 :-)

    • WeakForm

      OF7中的多相流求解器
      • WeakForm

      2
      2
      帖子
      73
      浏览

      李东岳

      主要是因为OpenFOAM那面多相流的工作量比较大,改进也比较多,新老版本都存在,所以就同时有了这么多东西。

      twoPhaseEulerFoam: 最老版本的双流体求解器,稳健性一般 multiphaseEulerFoam: 最老版本的多流体求解器,稳健性不好 reactingTwoPhaseEulerFoam: 末期版本的双流体求解器,稳健性比较好 reactingMultiPhaseEulerFoam: 末期版本的多流体求解器,稳健性比较好 twoLiquidMixingFoam 不属于上面几类

      想用IATE和类方法计算不含相变和传热的问题

      用reacting类求解器

    • 李东岳

      PK! OpenFOAM VS Fluent计算效率
      • 李东岳

      33
      33
      帖子
      9105
      浏览

      李东岳

      偶尔误入这个帖子,还挺有意思。

    • 知

      twoLiquidMixingFoam interMixingFoam 互溶流体动网格模拟
      • 知行

      3
      3
      帖子
      35
      浏览

      知

      @李东岳 十分感谢李老师,后面课题组有培训经费一定要来上您的编程课:huahua:

    • 吴

      如何在库文件中修正中心差分的表面格式
      • 吴琼

      1
      1
      帖子
      29
      浏览

      吴

      你好。想问一下在openfoam中离散过程中对于面插值的中心差分格式的源文件代码openfoam/src/finiteVolume/interpolation/surfaceInterpolation/surfaceInterpolationScheme/surfaceInterpolationScheme.C
      其中对于点乘面差值格式的代码形式主要由Foam::surfaceInterpolationScheme<Type>::dotInterpolate 函数决定的,我现在的主要想法是想通过修正插值格式在代码sfi[fi] = Sfi[fi] & (lambda[fi]*(vfi[P[fi]] - vfi[N[fi]]) + vfi[N[fi]] ) 后面加体心梯度项点乘偏移量向量,如图所示:
      b2b2b5cc-a725-48e9-ad35-ba8bc36d7208-1675996388057.png

      我该如何修改代码,或者说大致的思路,请各位指点,拜托了。
    • 尚

      湍流分辨率质量估计--OpenFOAM v2212的一个新功能
      • 尚善若水

      3
      3
      帖子
      66
      浏览

      尚

      @李东岳 是的李老师,这个功能挺有意思的。我前不久写过Pope和Celik Index的类似后处理代码,看到之后新功能就分享给大家了。当然我实现的比较简单,适用性没有那么好,但是足够用了。您说的写到controlDict的思路也可以,这个其实也就是个简单验证,也不需要多优雅~