CFD中文网

    CFD中文网

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

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

      1
      1
      帖子
      9148
      浏览

      李东岳

      更换源

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

      在国内下载OpenFOAM

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

      算力测试

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

      切换gcc版本安装OpenFOAM

      Ubuntu 20.04 安装 OpenFOAM-v1712

      杂七杂八

      waves2foam 安装

      OpenFOAM常用小代码

      CFDEM+OpenFOAM+LIGGGHTS安装过程

    • 李东岳

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

      46
      46
      帖子
      78506
      浏览

      黑山雁

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

    • 杨

      OpenFOAMoam7虚拟机求解器编译报错
      • 杨英狄

      2
      2
      帖子
      5
      浏览

      李东岳

      这个求解器应该是调用了原来位置的相对路径下面的createAlphaFluxes.H文件。你看一下Make/options,换成绝对路径之类的

    • T

      OpenFoam实时监测流场某一截面上组合变量问题
      • Try harder

      2
      2
      帖子
      185
      浏览

      H

      @Try-harder 请问这个问题解决了吗

    • H

      openfoam监测某个面上的涡量数据
      • hongjiewang

      1
      1
      帖子
      2
      浏览

      H

      请教大家:
      openfoam2106版本中如何监测某个面上的涡量值
      我的思路是在controdict中加入surfaceSampling和vorticity两个函数,发现确实时间文件里面有了涡量的结果,而且也出现了postprocessing文件夹,但是该文件夹的数据里面只有x y z坐标,而没有涡量的值,请问大家有什么办法吗,具体的Function如下

      surfaceSampling
      {
      type surfaces;
      libs ( "libsampling.so" );
      enabled true;
      timeStart 0;
      writeControl runTime;
      writeInterval 0.06;
      interpolationScheme cellPointFace;
      surfaceFormat raw;
      fields (vorticityField);
      surfaces
      (
      plane
      {
      type plane;
      planeType pointAndNormal;
      pointAndNormalDict
      {
      normal (0 0 1);
      point (0 0 0.1) ;
      }
      bounds (0.1 0.1 0.08) (0.6 0.4 0.12);
      }
      );
      }

      vorticity1
      {
      type vorticity;
      libs (fieldFunctionObjects);
      field U;
      result vorticityField;
      region region0;
      enabled true;
      log true;
      timeStart 0.25;
      timeEnd 1000;
      executeControl timeStep;
      executeInterval 5;
      writeControl runTime;
      writeInterval 0.06;
      }

    • 悬

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

      15
      15
      帖子
      4107
      浏览

      李东岳

      你得壁面函数用的哪个壁面函数?标准壁面函数还是Spalding那个

    • H

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

      7
      7
      帖子
      35
      浏览

      H

      @李东岳 谢谢李老师。

      这应该是TJunctionFan这个tutorial吧,这个确实能得到局部的速度增加,但是如果管道足够长的话,管道的入口和出口处仍然是相同的速度和流量。

      回复您之前的问题,是inlet,outlet这两个地方。我想实现的结果就是在outlet的速度比inlet速度高,这确实是一个不符合物理的问题。。流量不守恒。没有外力的情况下是不可能实现的。硬要做的话,好像是可以通过添加momentum source实现,但是确实是一个非物理的问题,打扰了。。

    • W

      求助:interFoam模拟2微米液滴碰撞平面,液滴内部产生气泡
      • wr_z

      3
      3
      帖子
      19
      浏览

      W

      @李东岳 谢谢老师的回复,我再排查一下错误!

    • L

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

      3
      3
      帖子
      43
      浏览

      学流体的小明

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

    • 沐

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

      2
      2
      帖子
      149
      浏览

      李东岳

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

    • W

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

      9
      9
      帖子
      67
      浏览

      李东岳

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

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

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

    • 四

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

      2
      2
      帖子
      29
      浏览

      李东岳

      Info改成Pout试试?

    • 李东岳

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

      5
      5
      帖子
      174
      浏览

      李东岳

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

    • 正

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

      2
      2
      帖子
      33
      浏览

      李东岳

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

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

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

      7
      7
      帖子
      80
      浏览

      李东岳

      @fangyuanaza :146: :146: :146:

    • M

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

      6
      6
      帖子
      48
      浏览

      李东岳

      @miying :146: :146: :146:

    • H

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

      60
      60
      帖子
      2415
      浏览

      C

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

    • 尚

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

      3
      3
      帖子
      24
      浏览

      尚

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

    • J

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

      10
      10
      帖子
      82
      浏览

      O

      @jasper-0 我觉得不需要

    • W

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

      3
      3
      帖子
      43
      浏览

      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
      帖子
      44
      浏览

      李东岳

      只能祝收敛..

    • 一

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

      14
      14
      帖子
      6714
      浏览

      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
      帖子
      187
      浏览

      沉

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

    • Y

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

      8
      8
      帖子
      2697
      浏览

      J

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

    • F

      python进行OpenFOAM流场后处理
      • fangyuanaza

      30
      30
      帖子
      1426
      浏览

      田畔的风

      @yuhxfoamer

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

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

      4
      4
      帖子
      106
      浏览

      李东岳

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

    • S

      自定义速度边界
      • soulx7

      2
      2
      帖子
      60
      浏览

      C

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

    • Z

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

      15
      15
      帖子
      1717
      浏览

      D

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

    • L

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

      22
      22
      帖子
      268
      浏览

      L

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

      谢谢李老师

    • 缄

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

      3
      3
      帖子
      42
      浏览

      缄

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

    • C

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

      12
      12
      帖子
      124
      浏览

      C

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

    • 流

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

      3
      3
      帖子
      80
      浏览

      流

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

    • 氕

      openfoam监测数据
      • 氕氘-川

      2
      2
      帖子
      56
      浏览

      李东岳

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

    • 星

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

      62
      62
      帖子
      3667
      浏览

      流

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

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

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

    • L

      blastXiFoam
      • Leeseungnan

      5
      5
      帖子
      288
      浏览

      李东岳

      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
      帖子
      81
      浏览

      S

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

    • 杨

      OpenFOAM上实现VOF-DPM
      • 杨英狄

      2
      2
      帖子
      80
      浏览

      杨

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

    • Y

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

      8
      8
      帖子
      1383
      浏览

      杨

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

    • chengan.wang

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

      1
      1
      帖子
      46
      浏览

      chengan.wang

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

      谢谢大家!

    • C

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

      3
      3
      帖子
      95
      浏览

      C

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

    • W

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

      8
      8
      帖子
      117
      浏览

      李东岳

      你先备份一下。

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

    • W

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

      3
      3
      帖子
      74
      浏览

      W

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

    • C

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

      4
      4
      帖子
      122
      浏览

      C

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

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

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

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

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

    • 学流体的小明

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

      4
      4
      帖子
      113
      浏览

      学流体的小明

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

    • D

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

      31
      31
      帖子
      5439
      浏览

      W

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

    • C

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

      10
      10
      帖子
      174
      浏览

      C

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

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

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

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