CFD中文网

    CFD中文网

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

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

      2
      2
      帖子
      63
      浏览

      李东岳

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

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

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

    • C

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

      13
      13
      帖子
      1618
      浏览

      李东岳

      后来这个解决了吗 :-)

    • WeakForm

      OF7中的多相流求解器
      • WeakForm

      2
      2
      帖子
      75
      浏览

      李东岳

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

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

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

      用reacting类求解器

    • 李东岳

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

      33
      33
      帖子
      9107
      浏览

      李东岳

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

    • 知

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

      尚

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

    • 悬

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

      16
      16
      帖子
      4291
      浏览

      学流体的小明

      nut 在壁面上使用壁面函数 nutUSpaldingWallFunction,会出现上面那样的速度分布图
      但另外一个${Re_{\tau}}$也等于1000的算例,得到的速度分布也差不多
      我等等重开一个帖子说明下问题

    • 是

      关于snappyHexMesh网格划分问题
      • 是嘿嘿牛啊

      4
      4
      帖子
      81
      浏览

      李东岳

      @是嘿嘿牛啊 haha 理解理解

    • vbcwl

      openfoam如何将结果文件写到其他文件夹?
      • vbcwl

      4
      4
      帖子
      78
      浏览

      李东岳

      你在机械硬盘跑,数据少的话,不影响跑的速度

    • N

      interFoam中怎样标识不同入口的流体
      • nbyjn

      1
      1
      帖子
      40
      浏览

      N

      000.jpg

      如图,是一个立方体空腔,腔体内是空气,立方体顶层有四个入口,每个入口注入的都是同种液体(比如水)。

      请问怎样能够比较方便的标识出入口信息呢?就是把inlet1注入的液体标识为1,inlet2注入的液体标识为2,以此类推。(不需要考虑液体混合的问题,只模拟一小段时间,到混合前就停止了)

      烦请李老师和各位大佬解答 @李东岳

    • 糊

      DPMFoam计算发散的问题
      • 糊里糊涂

      8
      8
      帖子
      743
      浏览

      D

      @李东岳 有两个问题,第一,计算稳态流场的情况下buoyantBoussinesqSimpleFoam可以采用这种方式吗?第二,这种方法是不是没计算离散相方程。

    • 一

      多重网格算法在代码实践问题
      • 一颗鸭蛋

      1
      1
      帖子
      43
      浏览

      一

      用多重网格算法求解线性方程组Ax=b,算法求解过程中按理只需要给定系数粗细网格的系数矩阵就可以了,但在看算法理论时,都涉及到粗细网格的位置关系的判定来确定插值矩阵、限定矩阵。这个就需要在算法上导入网格具体的物理位置了吧。这样好像不对(按理是只要给定系数矩阵A,b,就能求解Ax=b)对于这点一直想不明白:mihu:

    • 小

      超算上如何运行老版本openfoam
      • 小考拉

      7
      7
      帖子
      894
      浏览

      H

      @hy1112006 是的李老师,在cfd-online上找到试了一个方法,虽然有系统警告,但暂且管用:
      cafe1fe5-612f-45bc-adac-914ea277d477-image.png

      警告如下,目前不影响运行:
      7997c7d3-2fcc-47c1-a44b-c92d5a4c78f1-image.png

    • W

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

      11
      11
      帖子
      135
      浏览

      李东岳

      捕获.PNG

      看起来属于这部分内容,我之前手算过类似的东西(附录A.3 扩散系数的影响)。是,那是要多计算域求解。但是你这个还涉及到一个动边界的问题。其他sci里面怎么处理的?

    • 李东岳

      CFD性能怪兽之双机并联!
      • 李东岳

      2
      2
      帖子
      3190
      浏览

      R

      @李东岳 李老师真是无私啊,正好公司想把两台工作站并联,你这直接有现成的教程!

    • 杨

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

      2
      2
      帖子
      28
      浏览

      李东岳

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

    • T

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

      2
      2
      帖子
      218
      浏览

      H

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

    • H

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

      1
      1
      帖子
      29
      浏览

      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;
      }

    • H

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

      7
      7
      帖子
      79
      浏览

      H

      @李东岳 谢谢李老师。

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

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

    • W

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

      3
      3
      帖子
      55
      浏览

      W

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

    • L

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

      3
      3
      帖子
      79
      浏览

      学流体的小明

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

    • 沐

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

      2
      2
      帖子
      176
      浏览

      李东岳

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

    • 正

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

      2
      2
      帖子
      46
      浏览

      李东岳

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

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

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

      7
      7
      帖子
      109
      浏览

      李东岳

      @fangyuanaza :146: :146: :146:

    • M

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

      6
      6
      帖子
      82
      浏览

      李东岳

      @miying :146: :146: :146:

    • 尚

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

      3
      3
      帖子
      38
      浏览

      尚

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

    • J

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

      10
      10
      帖子
      111
      浏览

      O

      @jasper-0 我觉得不需要

    • W

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

      3
      3
      帖子
      61
      浏览

      W

      @李东岳 收到 多谢李老师

    • J

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

      2
      2
      帖子
      48
      浏览

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

      李东岳

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

    • 缄

      残差曲线
      • 缄默

      2
      2
      帖子
      68
      浏览

      李东岳

      只能祝收敛..

    • 一

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

      14
      14
      帖子
      6805
      浏览

      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}}$

    • Y

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

      8
      8
      帖子
      2743
      浏览

      J

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

    • F

      python进行OpenFOAM流场后处理
      • fangyuanaza

      30
      30
      帖子
      1695
      浏览

      田畔的风

      @yuhxfoamer

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

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

      4
      4
      帖子
      139
      浏览

      李东岳

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

    • S

      自定义速度边界
      • soulx7

      2
      2
      帖子
      85
      浏览

      C

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

    • Z

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

      15
      15
      帖子
      1805
      浏览

      D

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

    • L

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

      22
      22
      帖子
      349
      浏览

      L

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

      谢谢李老师

    • 缄

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

      3
      3
      帖子
      75
      浏览

      缄

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

    • C

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

      12
      12
      帖子
      156
      浏览

      C

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

    • 流

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

      3
      3
      帖子
      112
      浏览

      流

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

    • 氕

      openfoam监测数据
      • 氕氘-川

      2
      2
      帖子
      80
      浏览

      李东岳

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

    • 星

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

      62
      62
      帖子
      4111
      浏览

      流

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

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

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

    • L

      blastXiFoam
      • Leeseungnan

      5
      5
      帖子
      311
      浏览

      李东岳

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

    • 氕

      后处理问题
      • 氕氘-川

      3
      3
      帖子
      114
      浏览

      氕

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

    • D

      DPMFoam咋后处理看见颗粒
      • Dingcy

      2
      2
      帖子
      65
      浏览

      D

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

    • S

      codedfixedvalue使用问题
      • soulx7

      8
      8
      帖子
      110
      浏览

      S

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

    • 杨

      OpenFOAM上实现VOF-DPM
      • 杨英狄

      2
      2
      帖子
      113
      浏览

      杨

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

    • Y

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

      8
      8
      帖子
      1429
      浏览

      杨

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