CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新
    1. 主页
    2. 最新
    登录以发表
    • 全部板块
    • 全部主题
    • 新建主题
    • 关注主题
    • 未回复主题
    • 洱聿

      并行计算浮点报错
      OpenFOAM • • 洱聿

      2
      0
      赞同
      2
      帖子
      15
      浏览

      李东岳

      @洱聿 在 并行计算浮点报错 中说:

      单核可以运算。并行运算到某个时间步库朗数突然增大

      类似这种问题。有的时候会发生。很难debug。这个一般来自于矩阵求解器,你可以试试换一下,比如GAMG换成PCG。但不保证解决

    • G

      使用A100和V100GPU运行RapidCFD
      OpenFOAM • • gtian

      13
      0
      赞同
      13
      帖子
      291
      浏览

      李东岳

      @gtian 是否可以测试下200万网格的效果?https://www.cfd-china.com/topic/3988/200万网格并行算力测试-openfoam版本

    • 李东岳

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

      40
      0
      赞同
      40
      帖子
      1450
      浏览

      李东岳

      @疏影横斜水清浅 是,HW07本身计算雨很大,可以放大计算域,Yang2009这个我记得我尝试也不行,现在有点记不清了。

    • D

      鼓泡塔液位模拟
      Fluent • • dongsc

      16
      0
      赞同
      16
      帖子
      6298
      浏览

      阿

      @dongsc 感觉是速度过大的原因,可以设置小点试试看。另外,你的这个折算是按照Re相同折算的表观气速?

    • WeakForm

      低含气率下的泡状流
      Algorithm • • WeakForm

      4
      0
      赞同
      4
      帖子
      25
      浏览

      李东岳

      @WeakForm 意义不大。

    • M

      求助~Fluent PBM UDF-Luo聚并模型
      Fluent • • Maud

      7
      0
      赞同
      7
      帖子
      1011
      浏览

      C

      请问luo破碎核的UDF怎么编写呢?

    • Y

      进行cfd计算时对进出口进行延长是必要的吗
      Fluent • • Yan

      1
      0
      赞同
      1
      帖子
      18
      浏览

      Y

      纯理论问题的探讨。因为看文献时,同一个问题,有人延长有人没有延长。延长进出口是为了增加计算的稳定性还是单纯让进出口边界的流动更充分?或者是为了避免回流?如何判断是否需要对进出口进行延长?如果延长了进出口,那么进行后处理分析时需要去除延长部分吗?

    • T

      使用paraview时出现segmentation fault错误
      OpenFOAM • • TStones

      4
      0
      赞同
      4
      帖子
      32
      浏览

      李东岳

      @TStones 在 使用paraview时出现segmentation fault错误 中说:

      这个等于是子系统自身的原因造成的问题吗?

      是的。类似linux远程+x也会有这个问题。

    • L

      怎样输出与一个patch发生碰撞的颗粒的速度,位置信息
      OpenFOAM • • lys

      1
      0
      赞同
      1
      帖子
      18
      浏览

      L

      怎样输出与一个patch发生碰撞的颗粒的速度,位置信息

    • B

      8000w网格算力测试结果
      Fluent • • BINFeng

      9
      0
      赞同
      9
      帖子
      89
      浏览

      李东岳

      大佬, 你那个GPU占用率好像一些能到92%,一些能到71%,这个是在变动的么

      然后这个内存,好像GPU跑的时候都是占一半?狂吃内存啊。

    • 四

      multiphaseEulerFoam对相变的计算
      OpenFOAM • • 四季之中的隐星

      1
      0
      赞同
      1
      帖子
      24
      浏览

      四

      各位老师和学者好,我最近在尝试使用multiphaseEulerFoam来计算水的闪蒸过程。但是遇到了一些问题。对于纯水的两相流动(没有空气),应该是调用ThermalPhaseChangePhaseSystem

      type thermalPhaseChangeMultiphaseSystem; phases (steam water); phaseChange on; steam { type purePhaseModel; diameterModel isothermal; isothermalCoeffs { d0 3e-3; p0 1e5; } residualAlpha 1e-3; } water { type purePhaseModel; diameterModel constant; constantCoeffs { d 1e-4; } residualAlpha 1e-3; }

      在此基础上,后续的换热模型选用了RanzMarshall,但是我发现根本不会发生闪蒸(我在能量方程里对water相直接补充了热量,保证了温度上升)

      in EEqns.H

      forAll(fluid.anisothermalPhases(), anisothermalPhasei) { phaseModel& phase = fluid.anisothermalPhases()[anisothermalPhasei]; const volScalarField& alpha = phase; tmp<volScalarField> tRho = phase.rho(); const volScalarField& rho = tRho(); tmp<volVectorField> tU = phase.U(); const volVectorField& U = tU(); volScalarField q = alpha*rho*(U&g)-alpha*rho*(U&g); if(phase.name() == "water") { q[93] = 5.0E8; q[94] = 5.0E8; q[95] = 5.0E8; q[96] = 5.0E8; q[97] = 5.0E8; } ···········

      不会发生相变的原因确定于代码中,对于质量源项的计算有

      ThermalPhaseChangePhaseSystem.C

      // Interfacial mass transfer update { volScalarField& dmdtf(*this->dmdtfs_[pair]); volScalarField& Tf(*this->Tfs_[pair]); const volScalarField Tsat(saturationModelIter()->Tsat(thermo1.p())); const volScalarField L ( volatile_ != "none" ? this->Li(pair, volatile_, dmdtf, Tsat, latentHeatScheme::upwind) : this->L(pair, dmdtf, Tsat, latentHeatScheme::upwind) ); volScalarField H1(this->heatTransferModels_[pair].first()->K(0)); volScalarField H2(this->heatTransferModels_[pair].second()->K(0)); volScalarField dmdtfNew((H1*(Tsat - T1) + H2*(Tsat - T2))/L);

      其中,换热系数H1与2的计算
      RanzMarshall.C

      Foam::tmp<Foam::volScalarField> Foam::heatTransferModels::RanzMarshall::K(const scalar residualAlpha) const { volScalarField Nu(2 + 0.6*sqrt(pair_.Re())*cbrt(pair_.Pr())); return 6 *max(pair_.dispersed(), residualAlpha) *pair_.continuous().thermo().kappa() *Nu /sqr(pair_.dispersed().d()); }

      可以发现,当离散相为水蒸气(初始条件下为0),连续相为水时,输入的残差为0时,

      volScalarField H1(this->heatTransferModels_[pair].first()->K(0)); volScalarField H2(this->heatTransferModels_[pair].second()->K(0));

      这个H1和H2必然为0,因此计算所得到的相变量dmdtfNew也必然为0,无论液体的温度高过饱和温度多少。

      请问有老师碰到过这种情况吗,任何点拨都感激不尽:xinlei:

    • T

      浸没边界法边界条件报错
      OpenFOAM • • ThomasShi

      3
      0
      赞同
      3
      帖子
      47
      浏览

      T

      @李东岳 好的,谢谢李老师,我去看一下

    • 小

      多个网格合并问题?
      Meshy • • 小龙

      13
      0
      赞同
      13
      帖子
      9393
      浏览

      M

      @XIN 请问您解决了吗,我也遇到这个问题了

    • 数

      简版NS方程用于CFD有完备性问题?这篇论文有看法
      Algorithm • ns方程 本构方程 角动量 不守恒 • • 数星少年

      5
      0
      赞同
      5
      帖子
      86
      浏览

      数

      @李东岳 角动量与C成正比,它在时间上的导数就是受到的力矩,这里面写了粘性力矩; 压力应该没有力矩,因为是势能。另外套上了对流项,建立的方程

    • R

      openfoam中添加湍流激励
      OpenFOAM • • Rachel0096

      8
      0
      赞同
      8
      帖子
      103
      浏览

      李东岳

      @杨英狄

      pturb=0.39*(密度*湍动能)

      类似这种,可以通过这个实现,不需要复杂的操作

      const volScalarField& k = mesh_.lookupObject<volScalarField>("k");
    • B

      大家怎么看待十沣科技封装软件疑似套壳OpenFOAM
      C斯达克 • • bleach2012

      4
      0
      赞同
      4
      帖子
      1470
      浏览

      C

      @bleach2012 这两个链接现在都,失效了:chouchou:

    • Y

      centos平台运行并行setRefinedField及blastFoam报错
      OpenFOAM • • yyif

      7
      0
      赞同
      7
      帖子
      68
      浏览

      Y

      @李东岳 好的!十分感谢您

    • 杨

      如何在求解器中声明湍流相关信息,如湍动能k
      OpenFOAM • • 杨英狄

      1
      0
      赞同
      1
      帖子
      36
      浏览

      杨

      大家好,想请教一下各位,如何在求解器中,或者说是在createFields.H中声明湍流相关的场,如湍动能、湍流耗散率等。

      因为目前需要考虑湍流对汽液传质的影响,需要在压力方程中用到湍动能,所以想知道怎样在求解器中声明。

      感谢各位!

    • R

      出现Negative initial temperature,但全场并未出现负温
      OpenFOAM • • Rachel0096

      5
      0
      赞同
      5
      帖子
      40
      浏览

      R

      好的,李老师,我在回去多调试一下:zoule:

    • 李东岳

      200万网格并行算力测试(ANSYS Fluent版本)
      Fluent • • 李东岳

      29
      0
      赞同
      29
      帖子
      5811
      浏览

      百

      @李东岳 最近在跑作业,后面停歇时试试

    • 阿

      基于Fluent 中的UDF做双向流固耦合。
      Fluent • • 阿波罗

      5
      0
      赞同
      5
      帖子
      1329
      浏览

      C

      @吴优 您好,不知道可否得到您的联系方式,想向您请教一下。谢谢

    • U

      openFoam与Fluent模拟简单的Laval蒸汽喷管
      OpenFOAM • • upc_ngh

      20
      0
      赞同
      20
      帖子
      693
      浏览

      U

      @李东岳,感谢李老师挂念,其实当时算的不是喷管,而是引射器,喷管只是其中一个部分,后面因为引射器的整体结构我们无法设计,所以很难判断是模拟的问题还是设计结构的问题,后来就搁浅了。看看后面业主能不能提供现场数据,再算算看。

    • T

      Foam-extend 并行计算报错
      OpenFOAM • • ThomasShi

      7
      0
      赞同
      7
      帖子
      57
      浏览

      T

      @李东岳 好的吧,谢谢李老师!

    • vbcwl

      DPMFoam给粒子添加温度等
      OpenFOAM • • vbcwl

      4
      0
      赞同
      4
      帖子
      90
      浏览

      vbcwl

      @upc_ngh 感谢,用的点粒子模型,这样就够了

    • 李东岳

      bubble colume算例一个
      OpenFOAM • • 李东岳

      1
      0
      赞同
      1
      帖子
      55
      浏览

      李东岳

      bubblePlume.tar.xz

      相关描述请参考 http://dyfluid.com/openfoam.html

    • M

      Fluent S2S辐射模型一步收敛且残差为0
      Fluent • • miying

      1
      0
      赞同
      1
      帖子
      15
      浏览

      M

      在使用Fluent计算传热时,对于空腔(填充了空气介质)辐射,使用S2S辐射模型,计算不会报错,但是log文件中有个问题,辐射计算总是会一步收敛,且残差为0,觉得非常奇怪,请各位老师见仁见智!

      “Radiosity converged after 1 iterations
      Final radiosity residual = 0.000000e+00”

    • S

      OpenFOAM跨节点并行计算问题
      OpenFOAM • • Stan同学

      3
      0
      赞同
      3
      帖子
      32
      浏览

      S

      @李东岳
      感谢东岳老师的分享,我先照着该教程试一试

    • R

      运行算例出现segmentation fault
      OpenFOAM • • Rachel0096

      4
      0
      赞同
      4
      帖子
      35
      浏览

      李东岳

      Cp_[cellI]输出一下

    • K

      仿真计算得到的出口处的压力是指流体在出口处达到的压力还是说该处的压力恒定为工作压力?
      Fluent • • Kirknwpu

      7
      0
      赞同
      7
      帖子
      73
      浏览

      K

      @李东岳 刚才自己也计算了两个例子,确实如此,谢谢老师的指导,看来我一直对压力出口的理解有错误。

    • L

      HPC 节点反馈
      OpenFOAM • • leilei

      8
      0
      赞同
      8
      帖子
      73
      浏览

      L

      @李东岳
      谢谢李老师
      您的意思是:
      如“大网格”,那OF会自动超线性加速? 无需在提交HPC任务时开超线程?

    • 氕

      openfoam怎么获得最大温度位置处的速度?
      OpenFOAM • • 氕氘-川

      3
      0
      赞同
      3
      帖子
      56
      浏览

      氕

      @李东岳老师,这个是求某个标量场内最值的吧? 我是想先求出最大温度处的位置,然后根据这个位置在求出该位置处的速度

    • 学流体的小明

      贴个对volSymmTensorField进行操作的代码
      OpenFOAM • • 学流体的小明

      4
      0
      赞同
      4
      帖子
      60
      浏览

      李东岳

      感谢分享!!:146: :146: :146:

    • Z

      有大神精通Basilisk的吗?
      C斯达克 • • zhaocx

      10
      0
      赞同
      10
      帖子
      855
      浏览

      H

      @DoubleTrump goodbountyhunter

    • G

      【求教】OpenFoam模拟波浪的壁面函数设置
      OpenFOAM • 壁面函数 波浪 • • GenhouSUN

      3
      0
      赞同
      3
      帖子
      83
      浏览

      O

      个人觉得没必要考虑吧

    • 李东岳

      2D黎曼问题OpenFOAM算例
      OpenFOAM • • 李东岳

      1
      0
      赞同
      1
      帖子
      72
      浏览

      李东岳

      shockTube.tar

      riemann2d.png

    • M

      twoPhaseMixture类中添加muf()函数报错求助
      OpenFOAM • • Miss Zhou_DUT

      4
      0
      赞同
      4
      帖子
      671
      浏览

      R

      @Miss-Zhou_DUT 您好,请问您解决这个问题了嘛

    • 学流体的小明

      如何在controlDict中创建一个场,不需要修改求解器
      OpenFOAM • • 学流体的小明

      1
      0
      赞同
      1
      帖子
      52
      浏览

      学流体的小明

      话说这个latex代码怎么转义不了
      使用LES计算气泡槽道流,因为是多相流,某一点可能会有气泡经过,解析的雷诺应力$u\prime_i u\prime_j$会受到很大影响,所以我想对 $u\prime_i u\prime_j$ 进行加权处理。
      具体方法是:生成张量场,计算每一个cell每一个瞬时的$\alpha u\prime_i u\prime_j$,$(1-\alpha) u\prime_i u\prime_j$,$\rho u\prime_i u\prime_j$,然后使用fieldAverage进行系综平均。
      我想在controlDict中生成场$\rho u\prime_i u\prime_j$,然后每个时刻计算一下这个场,怎样通过functions{ }操作呢?因为是要每一个时间步都计算以获得正确的值,所以不太想按一定时间间隔输出数据之后再进行后处理。
      当然也可以直接在solver的createFields.H中创建,我就是想找一个更通用的不用修改求解器的方法,直接很简便地在controlDict - functions{ }中实现。
      毕竟functions{ }可以用writeObjects输出

      /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Version: v2012 \\ / A nd | Website: www.openfoam.com \\/ M anipulation | ------------------------------------------------------------------------------- Description Writes out specified objects, e.g. fields, stored on the case database. \*---------------------------------------------------------------------------*/ type writeObjects; libs (utilityFunctionObjects); objects (<object names>); writeControl writeTime; // ************************************************************************* //
    • 洱聿

      mpirun并行显示运算中,但是实际log文件无输出
      OpenFOAM • • 洱聿

      6
      0
      赞同
      6
      帖子
      108
      浏览

      洱聿

      @洱聿 目前主动停止算例运行,会出error文件(但是之前以为是主动停止报的错误,就一直没有注意过),显示错误如下:

      [manage01:228000] [[60655,0],0] tcp_peer_recv_connect_ack: received different version from [[60655,1],1]: 1.8.8 instead of 1.10.2

      目前集群使用的mpi版本是1.10.2,是因为mpi版本太高,所以foam extend4.0不适用嘛?有fe40专门对应的openmpi版本嘛?

    • R

      mvConvection->fvmDiv
      OpenFOAM • • Rachel0096

      2
      0
      赞同
      2
      帖子
      32
      浏览

      李东岳

      组分方程一般会有几十个。这样写的话,可以通过div(phi,Yi_h)定义相同的格式。

    • F

      如何获得固体边界上水相部分某一场的数值
      OpenFOAM • • Foamer

      2
      0
      赞同
      2
      帖子
      30
      浏览

      李东岳

      我觉得你可以写个判断语句。对整个边界patch做循环,比如边界上alpha>0.5的就是水箱,然后提取这个压力。

    • T

      cyclicAMI边界
      OpenFOAM • • ThomasShi

      2
      0
      赞同
      2
      帖子
      67
      浏览

      李东岳

      把你的boundary文件贴一下看看

    • 书生kao

      时间步长对计算结果的影响
      OpenFOAM • • 书生kao

      2
      0
      赞同
      2
      帖子
      107
      浏览

      李东岳

      带化学反应的时候存在刚性方程组。所以还不能像是普通流场那样放大时间步。

    • L

      关于icoReactingMultiphaseInterFoam中的组分方程
      OpenFOAM • • lrl3512

      3
      0
      赞同
      3
      帖子
      65
      浏览

      李东岳

      这在以前的openfoam里面是出现过这种问题的。我详细的跟踪过twoPhaseEulerFoam的算法演变。你这个求解器应该是ESI那面的吧?可以问问他们什么时候能把alpha、rho搞进去

    • I

      希望下版本openfoam能重写multiphaseEulerFoam的传质扩散模型
      OpenFOAM • • IO-Turbulent

      2
      0
      赞同
      2
      帖子
      71
      浏览

      李东岳

      这个肯定是没法提说不出口的。咱们要知道,官方开发openfoam40年了。他们对自己的代码/算法非常自信。除非明显的植入错误。他们不会做类似处理。我比较熟悉基金会那面。类似这种需求,做一个统一化的架构整合,属于技术支持,提出来他们会需要你提供fund,然后他们会处理。

      ESI那面我联系不多。不过这种大概率来看。自己整合一个然后提交给他们会比较好。

    • 1

      fluent动画设置如何添加一个时间标识?
      Fluent • • 13673619804

      4
      0
      赞同
      4
      帖子
      99
      浏览

      I

      @13673619804 如果是图片的话,可以用imagJ 时间戳功能

    • 建

      ICEM-CFD出现这种情况咋调整过来呀
      Fluent • • 建筑外环境通风模拟小学生

      3
      0
      赞同
      3
      帖子
      111
      浏览

      I

      @建筑外环境通风模拟小学生 多截几张图

    • L

      关于在compressibleInterFoam中添加组分方程
      OpenFOAM • • lrl3512

      4
      0
      赞同
      4
      帖子
      140
      浏览

      R

      您好,请问一下您实现在compressibleInterFoam中添加组分方程了嘛:chitang:

    • S

      大佬们这句该怎么礼貌而又得体的喷死他
      C斯达克 • • shrine

      4
      0
      赞同
      4
      帖子
      86
      浏览

      李东岳

      @shrine 在 大佬们这句该怎么礼貌而又得体的喷死他 中说:

      本文应指出本文工作内容与Openform软件结合的优越性和必要性,包括二次开发等介绍

      类似的,肯定会有审稿人说:

      本文应指出本文工作内容与CFD结合的优越性和必要性

    • 李东岳

      2000万网格OpenFOAM算例
      OpenFOAM • • 李东岳

      1
      0
      赞同
      1
      帖子
      167
      浏览

      李东岳

      2000w.tar.xz

      加特林,2000步,906秒

      20000w.tar.xz

    • F

      后处理得到momentum thickness
      OpenFOAM • • fangyuanaza

      7
      0
      赞同
      7
      帖子
      442
      浏览

      F

      @luofq-sysu 第一个问题:是的,我取的是流向速度时均值。
      第二个问题,是的,是不同流向位置。我当时做的是一个平板流动,平板上流向的网格113个。150个不是一个必须的值,看你算例的情况跟需求。我的是二维流动,因此没有展向的问题。