CFD中文网

    CFD中文网

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

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

      1
      1
      帖子
      9135
      浏览

      李东岳

      更换源

      默认的源下载什么都非常慢,请更改国内的源,安装openfoam依赖包5秒完事。

      在国内下载OpenFOAM

      国外的OpenFOAM安装包下载也是特别慢,直接在此处下载国内OpenFOAM安装包

      算力测试

      装好OpenFOAM之后,可以在这里测试下你的算力

      切换gcc版本安装OpenFOAM

      Ubuntu 20.04 安装 OpenFOAM-v1712

      杂七杂八

      waves2foam 安装

      OpenFOAM常用小代码

      CFDEM+OpenFOAM+LIGGGHTS安装过程

    • 李东岳

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

      46
      46
      帖子
      78465
      浏览

      黑山雁

      @李东岳 李老师好啊,这个页面上几个 dyfluid.com 的链接是不是需要更新下啊,好像路径已经发生变化了,打不开了:136:

    • H

      请问一个流速随距离增加的管流如何实现?
      • hotairballon1997

      5
      5
      帖子
      24
      浏览

      李东岳

      你这个fan那面的外圈 留一点空隙 我估计会好一些 要不然质量不守恒

    • L

      求助:两相VOF求解器中当一相是可压的,另一相不可压时,压力泊松方程如何构建
      • lrl3512

      3
      3
      帖子
      35
      浏览

      学流体的小明

      @lrl3512 compressibleInterFoam可以参考这个
      https://zhuanlan.zhihu.com/p/66444729
      我也是看这个学的

    • 悬

      q-DNS计算槽道流遇到了一些问题,求大神们指点
      • 悬铃神木

      14
      14
      帖子
      4099
      浏览

      学流体的小明

      @XieXiaoyang 您好,我想请教一下,像这种 DNS / LES 槽道流的结果,我想画速度剖面,进行正则化时候的摩擦速度${u_{\tau}}$和粘性长度${\delta_{\nu}=\nu/u_{\tau}}$是怎样计算的?
      我当前的做法是:使用postChannel(就像channel395那个算例一样)将空间时间都做平均之后,可以得到${u(y)}$这样一个分布,然后计算壁面上的速度曲线的斜率,计算壁面摩擦应力,计算摩擦速度:
      \begin{equation}
      {u_\tau } = \sqrt {\frac{{{\tau _w}}}{\rho }} {\rm{ = }}\sqrt {\nu \frac{{\partial u}}{{\partial y}}}
      \end{equation}
      但算出来的摩擦速度和预设的相差很大。所以最后的曲线是这样的
      8544ca30-1b5b-4552-9e71-567bb4fe6134-channel1000DirectCalUtao.jpg
      算了很久,对数率那里也无法和理论解对上。计算结果的整个曲线也与理论解的形状不符合。
      感觉我使用的方法不太对,想请问一下您的做法。

    • 沐

      Tobias书中关于pimple的final松弛因子的设定
      • 沐风

      2
      2
      帖子
      142
      浏览

      李东岳

      pFianl控制连续方程的最终残差来保证速度一进一出是守恒的。在瞬态的情况下要保证时间步准确,pFianl不能添加任何低松弛。PIMPLE通过松弛确实可以实现大时间步的跨越,如果外迭代有足够多的次数,可以保证大时间步后的收敛解。这样可以解释的通。但保守起见,我更倾向把fianl的松弛去掉。

    • W

      OF两相流中是否能够实现只在气相求解某一个方程
      • WangQi 0

      9
      9
      帖子
      66
      浏览

      李东岳

      @WangQi-0 openfoam里面存在一些多计算域模拟,比如chtMultiRegionFoam。

      就是想把界面的以及界面以上的液相电势全部设置为0

      你具体这个是算什么?求解液相电势方程么。上面是空气还是其他的?

    • 四

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

      2
      2
      帖子
      27
      浏览

      李东岳

      Info改成Pout试试?

    • 李东岳

      中性大气环境湍流动能的自保持
      • 李东岳

      5
      5
      帖子
      165
      浏览

      李东岳

      @Hope 多谢多谢!我也研究研究。:duang:

    • 正

      OpenFOAM的密度和通量量纲不对应问题
      • 正应力

      2
      2
      帖子
      30
      浏览

      李东岳

      你的phi需要乘以密度才行。目前来看你的phi就是网格乘以速度

      phi *= fvc::interpolate(rho);
    • F

      考虑动量源项后残差变差
      • fangyuanaza

      7
      7
      帖子
      77
      浏览

      李东岳

      @fangyuanaza :146: :146: :146:

    • M

      关于OpenFOAM-9中limitTemperature的用法
      • miying

      6
      6
      帖子
      47
      浏览

      李东岳

      @miying :146: :146: :146:

    • H

      钝体建筑扰流的大涡模拟
      • HITSC30

      60
      60
      帖子
      2392
      浏览

      C

      @李东岳 感谢李老师分享,看了这篇文章主要适用于RANS模型,目前对这块还没有研究。我主要研究的是LES大涡模拟入口湍流生成方法,由于入口湍流生成方法有很多种,生成的风场性质也都不同,湍流发展过程可能也有不一样的特点。针对不同入口风场,源项、或壁面函数、或粗糙元布置方式应该也要对应修改,使得每个瞬时的近地面速度发生预期的调整,最后要求计算得到的平均风和湍流强度符合预期值,目前感觉实现一种通用的设置方式还是比较困难。

    • 尚

      定义了一个新的输出量,但是编译错误error: no viable overloaded '='
      • 尚善若水

      3
      3
      帖子
      24
      浏览

      尚

      @李东岳 是的,编译通过了。谢谢李老师!

    • J

      植入拉格朗日粒子的Ueqn中需不需要加phicForces?
      • Jasper 0

      10
      10
      帖子
      79
      浏览

      O

      @jasper-0 我觉得不需要

    • W

      求助:两相流求解拉普拉斯方程时,气液粘性差距太大如何提高精度
      • WangQi 0

      3
      3
      帖子
      37
      浏览

      W

      @李东岳 收到 多谢李老师

    • J

      如何让injectionModels在每一个时间步都读取kinematicLookupTableInjection
      • Jasper 0

      2
      2
      帖子
      34
      浏览

      J

      @jasper-0 更新一下目前进度,我尝试在KinematicLookupTableInjection.C里将injectors_设置为MUST_READ_IF_MODIFIED,

      template<class CloudType> Foam::KinematicLookupTableInjection<CloudType>::KinematicLookupTableInjection ( const dictionary& dict, CloudType& owner, const word& modelName ) : InjectionModel<CloudType>(dict, owner, modelName, typeName), inputFileName_(this->coeffDict().lookup("inputFile")), duration_(this->coeffDict().getScalar("duration")), parcelsPerSecond_(this->coeffDict().getScalar("parcelsPerSecond")), randomise_(this->coeffDict().getBool("randomise")), injectors_ ( IOobject ( inputFileName_, owner.db().time().constant(), owner.db(), IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE ) ),

      但是在运行solver的时候提示不支持自动重新读取,有没有大佬知道这种情况咋解决呢?

      GlobalIOList<kinematicParcelInjectionData> kinematicLookupTableInjection constructed with IOobject::MUST_READ_IF_MODIFIED but GlobalIOList<kinematicParcelInjectionData> does not support automatic rereading.
    • D

      颗粒在壁面没有反弹?
      • Dingcy

      2
      2
      帖子
      24
      浏览

      李东岳

      应该是drag太大了,颗粒没弹起来

    • 缄

      残差曲线
      • 缄默

      2
      2
      帖子
      43
      浏览

      李东岳

      只能祝收敛..

    • 一

      LES介绍的文章的一个公式
      les • • 一二

      14
      14
      帖子
      6705
      浏览

      G

      @一二 在 LES介绍的文章的一个公式 中说:

      嗯,下面是我推的
      \begin{equation}
      \overline{\mathbf{S}}:\mathbf{\tau}+C_{e}\frac{k_{sgs}^{1.5}}{\Delta}=0 \
      left = \overline{\mathbf{S}}:[\frac{2}{3}k_{sgs}\mathbf{I}-2\mu_{sgs}\mathbf{dev}(\overline{\mathbf{S}})]+k_{sgs}^{1.5}\frac{C_{\epsilon}}{\Delta}
      \end{equation}
      因为
      \begin{equation}
      \mu_{sgs}=C_{k}\Delta\sqrt{k_{sgs}} \
      = \overline{\mathbf{S}}:[\frac{2}{3}k_{sgs}\mathbf{I}-2C_{k}\Delta\sqrt{k_{sgs}}]+k_{sgs}^{1.5}\frac{C_{\epsilon}}{\Delta} \
      = \sqrt{k_{sgs}}(\frac{C_{\epsilon}}{\Delta}k_{sgs}+\frac{2}{3}\mathbf{tr}(\overline{\mathbf{S}})\sqrt{k_{sgs}}-2C_{k}(\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}})) \
      = ak_{sgs}+b\sqrt{k_{sgs}}-c\
      = right = 0
      \end{equation}
      其中
      \begin{equation}
      a=\frac{C_{\epsilon}}{\Delta} \
      b=\frac{2}{3}tr(\overline{\mathbf{S}}) \
      c=2 \Delta C_{k}(\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}} \
      \sqrt{k_{sgs}}=\frac{-b+\sqrt{b^2+4ac}}{2a}
      \end{equation}
      当为不可压缩流体时$tr{\overline{\mathbf{S}}}=0$,那么$b=0$、$c=2 \Delta C_{k} (\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}})$,就可以得到$k_{sgs}=\frac{c}{a}=\frac{2C_{k}||\overline{\mathbf{S}}||^{2}{}\Delta^2}{C_{\epsilon}}$

    • C

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

      2
      2
      帖子
      185
      浏览

      沉

      类似分离涡这种大尺度的涡结构,WMLES即可,涉及转拟湍流,湍流结构等,三方向的尺度需满足计算要求。

    • Y

      怎么编写法向二阶梯度为0的边界条件
      • yang山青

      8
      8
      帖子
      2697
      浏览

      J

      @yang山青 您好,请问您这个问题最后解决了嘛?

    • F

      python进行OpenFOAM流场后处理
      • fangyuanaza

      30
      30
      帖子
      1414
      浏览

      田畔的风

      @yuhxfoamer

      首先pathlib这个模块在最终的代码中并没有被使用到,所以你可以将这个模块和对应的代码移除,这并不会影响代码运行结果。 然后谈谈你这个报错。readData.py这个文件是用pvpython执行的,这个python解释器是paraview自带的,对应的第三方库也是在安装paraview时配置好的,所以错误来源于你的paraview附带的pvpython解释器(可能是安装问题,可能是系统环境问题,也可能是软件版本问题,请自行排查),和系统python解释器无关。你使用pip和conda只是在安装你系统python解释器的第三方库。 尽管最终版本这个模块被弃用,不过这里再补充一下:pathlib在本代码中的用途是创建空文件夹(pathlib.Path),在上面的讨论中也提到了,因此可以使用os.mkdir替代。
    • R

      请教关于浸入边界法(IBM方法)的问题
      • Rachel0096

      4
      4
      帖子
      99
      浏览

      李东岳

      你这个看起来就是压力波的传递导致的。近音速可能会这样。openfoam有无反射边界条件,不过用起来效果一般。你可以试试调整出口的边界条件,也可以试试把计算域延长一些

    • S

      自定义速度边界
      • soulx7

      2
      2
      帖子
      59
      浏览

      C

      @soulx7 可以参考这篇文章实现:OpenFOAM|20 自定义边界条件

    • Z

      能不能设置两组颗粒(一组one-way,一组four-way)同时模拟?
      • zhangxc0223

      15
      15
      帖子
      1713
      浏览

      D

      @dingcy 关掉耦合应该是oneway,也可以这么用,要是DPMFoam满足你使用要求的话就可以。

    • L

      HPC上运行OF
      hpc上运行of • • leilei

      22
      22
      帖子
      262
      浏览

      L

      @李东岳 :
      太感谢李老师了,明白了,谢谢。也谢谢曾帮过我的朋友。
      一个玩笑:
      把CFD大佬比为大城市美女,那咱乡下穷傻小子的心情就如那句歌词:
      。。。would it be a sin, when i cannot help falling in love with u ........(有可能记忆有出入)

      谢谢李老师

    • 缄

      OpenFOAM并行计算中暂停后如何继续计算
      • 缄默

      3
      3
      帖子
      41
      浏览

      缄

      @李东岳 谢谢李老师回复,麻烦您了。我知道我哪里错了是我太粗心了。

    • C

      多个作业分配到同一个节点上会导致计算效率大幅降低
      • cccrrryyy

      12
      12
      帖子
      122
      浏览

      C

      @李东岳 可以可以,我找时间试试看。

    • 流

      移植sprayFoam液相功能
      • 流体中二少年

      3
      3
      帖子
      79
      浏览

      流

      @tens 不好意思,隔这么久回复你,求解器已经改好了,问题在于矩阵的形式不同,接下来准备验证求解器的可靠性:chigua:

    • 氕

      openfoam监测数据
      • 氕氘-川

      2
      2
      帖子
      55
      浏览

      李东岳

      还真没这么玩过。看看其他大佬怎么说

    • 星

      Euler-Lagrange的一些解析,sprayFoam
      • 星星星星晴

      62
      62
      帖子
      3639
      浏览

      流

      刚接触sprayfoam,请教三个个问题:

      初始场指定区域怎么布置液滴?setFields? 液滴的喷射位置和喷射速度怎么设置为跟随inlet边界? 我模拟的是预气化的液滴,比如C7H16既有液滴也有气体,我在算例里面该怎么区分他们?

      谢谢前辈,希望解答或者提供一些解决问题的方向:140:

    • L

      blastXiFoam
      • Leeseungnan

      5
      5
      帖子
      287
      浏览

      李东岳

      blastFoam这个求解器里面的算法还挺有意思的。求解器的作者也很厉害。我一直在关注。

    • 氕

      后处理问题
      • 氕氘-川

      3
      3
      帖子
      83
      浏览

      氕

      @李东岳 效果非常好,非常感谢李老师:140:

    • D

      DPMFoam咋后处理看见颗粒
      • Dingcy

      2
      2
      帖子
      47
      浏览

      D

      自问自答,我粗略的整出来了,如图1,要注意画红圈地方的选择:chitang:
      746340cd-5bc6-47e8-a92f-b507dd49d9a8-image.png
      图1

    • S

      codedfixedvalue使用问题
      • soulx7

      8
      8
      帖子
      79
      浏览

      S

      @wangfei9088 你好,哈哈昨天晚上搞得太晚神智有点不清了,刚改了编译通过了,多谢

    • 杨

      OpenFOAM上实现VOF-DPM
      • 杨英狄

      2
      2
      帖子
      78
      浏览

      杨

      @杨英狄 就是说想以这样的方式注入拉格朗日粒子:每个时间步注入流场的粒子的数量、位置都与当前时间步的流场有关,这样的注入方式怎么实现呢?

    • Y

      拉格朗日粒子包(parcel)入注是什么时候添加到粒子云(cloud)的
      • yfclark

      8
      8
      帖子
      1378
      浏览

      杨

      @yfclark 你好,你这个功能实现了吗,可以请教一下如何实现的吗?

    • chengan.wang

      paraview做多个切面,然后在每个面上对物理量积分
      • chengan.wang

      1
      1
      帖子
      43
      浏览

      chengan.wang

      大家好,目前我可以输出温度场的三维数据,比如
      QQ截图20221222151315.png
      我想垂直于z轴做100个等间距垂直切面,在每个切面上对温度求积分并取平均值。
      \begin{equation}
      \int_{A} T dx dy /A
      \end{equation}
      然后沿着z轴画出一条温度平均值曲线。

      谢谢大家!

    • C

      openfoam允许在一段程序建立两个或多个fvMatrix吗?
      • chaoscfd

      3
      3
      帖子
      94
      浏览

      C

      @李东岳 李老师您好,D()是一个volVectorField。主要是发现建立两个fvMatrix时,前面一个fvMatrix会对后面的fvMatrix产生影响,主要体现在源项上。按理说这是两个独立的fvMatrix对象,不应该互相影响啊

    • W

      在OpenFOAM7中无法编译Benjamin的OF7的makeAxialMesh
      • WangGH

      8
      8
      帖子
      116
      浏览

      李东岳

      你先备份一下。

      然后你wclean一下,看一下这些东西会不会消失。消失的话再编译一下

    • W

      OF能否将计算的每时刻的速度势进行输出?
      • WangGH

      3
      3
      帖子
      73
      浏览

      W

      @李东岳 好滴,感谢李老师,我编一下

    • C

      VOF方法中如何准确计算自由表面曲率?
      • CFDngu

      4
      4
      帖子
      115
      浏览

      C

      @李东岳 在 VOF方法中如何准确计算自由表面曲率? 中说:

      这个方法在函数calculateK()里面,可直接调用。计算结果在3232网格上能得到光滑的结果,但不对。在6464网格上出现了数值振荡。

      VOF结果在网格比较细的时候结果不对?这个结论很有挑战性

      李老师,我用interIsoFoam算了一下,还是不对。在32$\times$32的网格上可以得到光滑的结果,但是和解析解不一样;在64$\times$64的网格上结果出现了振荡。interfoam也跑了一下,结果和interisofoam不一样,也不对。。。(这个算例的设置参照和解析解在主贴的参考文献[1]中)

      算例在这里,大家有兴趣可以跑一下。要是算例设置有问题也可以指出来。另外,我是在v2006环境下跑的(这个链接可能需要注册坚果云,我没找到其他的上传方法:xinlei: )
      capillaryWave

    • 学流体的小明

      paraview查看自适应加密网格出错
      • 学流体的小明

      4
      4
      帖子
      111
      浏览

      学流体的小明

      @李东岳 并行文件夹里面的2e-08文件夹是有pointProcAddressing文件的,但是后面算出来的2e-02,4e-02就没有这个文件了。所以reconstructPar说找不到这个文件,报错了。

    • D

      请问.msh的2维旋转轴对称网格如何转换到OpenFOAM
      • D.Benjamin

      31
      31
      帖子
      5437
      浏览

      W

      @exthan 您好,我刚刚也下载了OF7版本的makeAxialmesh,但好像没找到编译文件,没有Allmake文件或者wmake文件,请问您是怎么编译的呢d27bf862-45db-4120-8ba1-050e18d46128-image.png

    • C

      每个时间步文件夹下没有数据写入,求大神!
      • Chen_hao

      10
      10
      帖子
      171
      浏览

      C

      @bestucan 谢谢老师,我试了下,还是不行~

      他是由于none量 他把它删了,看来我的问题就是那个omega,只要用

      const volScalarField& omega_ = ( this->mesh_.objectRegistry::template lookupObject<volScalarField>("omega") );

      就会找不到,估计是结构问题,难道只能重写了。

    • R

      分区用simpleFoam求解器出现的问题,求问各位这样的报错应该如何解决?
      • Rachel0096

      3
      3
      帖子
      80
      浏览

      R

      @李东岳 李老师您好,感谢您提出的意见,现在只是算的有些慢但是似乎问题已经解决了,谢谢。33.png

    • J

      Not enough slots available
      • Jacob Z

      2
      2
      帖子
      463
      浏览

      S

      @jacob-z

      mpirun --oversubscribe -np 20 xxFoam -parallel

      应该是进程超限锁死,需要添加申请允许超限,就可以了

    • 李

      OpenFoam后处理读取颗粒错误
      • 李梦瑶

      8
      8
      帖子
      876
      浏览

      T

      还有种方法,通过一个工具将OF中颗粒坐标转换成旧版本的格式,然后就能被paraview和tecplot等读取了
      https://github.com/blueCFD/lagrangianExtraFunctionObjects

    • W

      关于源项是放速度方程中还是压力泊松方程中处理的问题
      • WangQi 0

      3
      3
      帖子
      76
      浏览

      W

      @李东岳 多谢李老师,我的物理问题中速度是一定会求解的。所以放在泊松方程中能够更加稳定,避免震荡。
      明白了,多谢您