Skip to content

OpenFOAM

OpenFOAM交流区

5.3k 主题 31.4k 帖子
  • 免费帮助debug算例

    4
    4 帖子
    1k 浏览
    liujmL

    @李东岳 李老师好,这里我请教一个问题->链接

    为了简化算例,我把网格三维转化为二维,简化了结构形式,以便排错。

    【算例云图】
    6101495a-adbe-455d-ab41-45e2b1b7b415-image.png
    bb59f1c1-9dde-4169-83d9-cf00c0daf433-image.png

    【简化版算例下载】
    functionalTest_NoPZRefine2_2D.zip
    生成网格createMesh,
    跑算例runCaseOnly,
    全流程runCase

    感谢~

  • 10 帖子
    4k 浏览
    H

    @hy1112006 您好,如果方便的话可以有您的邮箱或者联系方式交流一下使用pisoFOAMopenfast吗?

  • 1 帖子
    281 浏览

    想尝试与一篇文献中实验结果作对比,但根据文中所提供的条件算出的入口速度很小(0.18m/s),以此作为入口的边界条件;同时入口边界的alpha.water也固定为1,计算一段时间后如图所示,发现入口处几乎没有水相的存在了。这是为什么?

    ae0a7d49-651d-4636-a27b-c8a0ff4b9334-image.png

    入口处的边界条件如下:

    U { type fixedValue; value uniform (0.18 0 0); } p_rgh { type zeroGradient; } alpha.water { type fixedValue; value uniform 1; } alpha.air { type zeroGradient; }
  • openfoam中DPM使用LES框架下的颗粒亚格子尺度

    2
    2 帖子
    349 浏览
    李东岳

    我觉得这个是一个自圆其说的东西。如果看sci,肯定有人对比LES和DNS的结果,证明他们的LES或者DNS的结果更好。

    但我认为,本身欧拉拉格朗日里面存在了很多模化,在粒子受力的角度,存在更多的模化,尤其是湍流分散力那面,基本都是参数。本身DPMFoam和DNS就是冲突的东西。所以如果做这方面研究,自圆其说是最重要的。

  • 4 帖子
    631 浏览
    李东岳

    https://cfd-china.com/topic/8018

    能做的话可以上传到这里我给你看看

  • 关于OpenFOAM新field-based Lagrangian系统的询问

    7
    7 帖子
    768 浏览

    @余正东 这个我关注了,未来可期,目前只是初步移植了一些

  • 9 帖子
    2k 浏览
    学流体的小明

    @先先贝 我的方法都在之前写的帖子里面分享出来了。算几个时间步很快的。而且串行计算三个时间步就差不多了,你设置的自适应加密等级也基本到位了。这个时候再运行一次setFields,给加密后的较密的网格赋值,再去分块。
    分块的话也得想想加密的部分会不会在你流场里面移动,你的不同核之间的负载要基本平衡,可能还需要你手动去划分区域。

  • 1 帖子
    189 浏览
    L

    OpenFOAM中使用multiDimAMR编译和运行正常,使用reconstructParMesh也正常,但是之后再使用reconstructParMesh会报错。
    第一个时刻还是正常的,第二个时刻发生报错,还有一个算例是第5个时刻开始报错,报错内容一样。使用paraview中的decomposePar也能查看结果,为什么会出现这个问题呢~
    Reconstructing fields for mesh region0
    Time = 0.01
    Reconstructing FV fields
    Reconstructing volScalarFields
    alpha.water
    alpha.water_0
    p
    p_rgh
    rAU
    Reconstructing volVectorFields
    U
    Reconstructing surfaceScalarFields
    alphaPhi0.water
    phi
    Reconstructing surfaceVectorFields
    Uf
    Reconstructing point fields
    No point fields
    No lagrangian fields
    No FA fields

    Time = 0.02

    --> FOAM Warning :
    From virtual Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
    in file meshes/polyMesh/polyMeshIO.C at line 210
    Number of patches has changed. This may have unexpected consequences. Proceed with care.
    --> FOAM Warning :
    From virtual Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
    in file meshes/polyMesh/polyMeshIO.C at line 210
    Number of patches has changed. This may have unexpected consequences. Proceed with care.

    **--> FOAM FATAL ERROR: (openfoam-2012)
    Processor 2 has a different polyMesh at time 0.02 compared to any previous processors.
    Please check time 0.02 directories on all processors for consistent mesh files.

    From Foam::polyMesh::readUpdateState Foam::processorMeshes::readUpdate() in file processorMeshes.C at line 192.**

    FOAM exiting

  • 3 帖子
    322 浏览
    A

    刚好在用kinematicParcel ,顶大佬一波!!

  • OpenFOAM膨胀喷管边界条件设置

    16
    16 帖子
    3k 浏览
    X

    大佬们,这个算例用rhoCentralFoam能跑通。因为后面想加入颗粒计算,现在在同样的边界条件下我用rhoPimpleFoam计算时,却一直卡在这里,入口的压力也不对,总压给定的是3e6。这是什么问题呢?
    image.png

    image.png

    image.png

  • 7 帖子
    982 浏览
    学流体的小明

    感谢大佬回复@coolhhh

    这两个图分别对应 方法1 和 方法2 的计算结果吗,还是不同方向的结果?

    不同方向的结果。$\xi$是$x$方向的间距,$\zeta$是$z$方向的间距。上面我画的都是相关系数,也就是相关函数进行了归一化,所以在间距$\xi$或$\zeta$为0时,相关系数都为1。

    看起来对于槽道流这种展向和顺流向是均匀+时间平稳的流动是可以用时间平均之后的。问AI,AI也是这么说的。

    我的数据在时间上和空间上也比较多,$\phi(x,z,y=y_0,t)$每个计算网格的值都输出了,时间上也以很小的时间间隔输出了一段时间。无论是时间上的相关系数还是空间上的,计算起来都没什么难度。

  • 16 帖子
    4k 浏览
    A

    东岳流体虚拟机4月更新, 配套ISO镜像同样上传至百度云. 安装方法不再赘述.

  • OpenFOAM自适应网格和并行计算冲突的问题

    14
    14 帖子
    9k 浏览

    @五好青年 大佬,后来者大胆问下这个有解决的好办法或者好思路了么

  • AUSM格式求解发散

    4
    4 帖子
    618 浏览

    @李东岳 是的。我又从头到尾排了一边,好像是求解马赫数的方式有点问题,现在还在验证,希望是这个原因

  • MPPICDyMFoam 中添加颗粒后计算一直卡住

    3
    3 帖子
    350 浏览
    余正东

    @李东岳 感谢李老师,我准备转向org版本进行研究,采用全新求解器重新实操一下

  • OpenFOAM计算复杂外形报错

    1
    1 帖子
    242 浏览
    S

    求助,目前在使用hy2Foam求解器计算类似乘波体外形,已经计算过钝头算例可以得到收敛解。但目前计算外形相对复杂,二维计算,计算域为1*0.55m,网格数量42000。计算过程中经常报错,求问有没有什么调试建议。下面是报错信息。

    Mean and max Courant Numbers = 5.29417e-05 0.0491796 deltaT = 3.99074e-11 Time = 2.017990284e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N2, Initial residual = 0.000161695, Final residual = 3.0494e-13, No Iterations 1 DILUPBiCG: Solving for O2, Initial residual = 0.000148597, Final residual = 1.8029e-13, No Iterations 1 DILUPBiCG: Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 1.07385e-06, Final residual = 1.49309e-06, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 8.83461e-07, Final residual = 1.19425e-06, No Iterations 1000 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 2.94288e-05, Final residual = 3.15326e+222, No Iterations 400 Attempt to use rho2ReactionThermo out of temperature range 11152 times during this iteration. -1.64159e+236 < Tlow: 200 Phase no 1.0 ExecutionTime = 3564.14 s ClockTime = 3564 s Iteration no 4641 (0.79 s) Mean and max Courant Numbers = 5.41311e-05 0.0489379 deltaT = 4.07806e-11 Time = 2.018031064e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N2, Initial residual = 0.000162296, Final residual = 1.67236e-16, No Iterations 2 DILUPBiCG: Solving for O2, Initial residual = 0.000148383, Final residual = 6.70388e-18, No Iterations 2 DILUPBiCG: Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:? #8 Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at ??:? #9 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #10 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #11 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #13 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"
  • 2 帖子
    230 浏览
    李东岳

    这两个不能叠加使用。

  • 如何修改并编译kOmegaSSTBase?

    11
    11 帖子
    4k 浏览
    H

    @杜浚平 十分感谢!

  • mirrorMesh做多区域网格的问题

    1
    1 帖子
    207 浏览
    N

    大家好!
    目前我做一个几何对称的共轭传热模拟,所以我只建模一半后面采用mirrorMesh做镜像,但是遇到问题
    bolckMesh→mirrorMesh→splitMeshRegions有问题,划分不了对称后的区域;
    bolckMesh→splitMeshRegions→mirrorMesh也有问题,区域对称过来了,但是没有划分不同区域;
    具体应该有说明操作,有人有相关经验吗?感谢!

  • Pimple中对速度场的操作

    10
    10 帖子
    3k 浏览
    R

    感谢老师们!这里确实是多变量的计算,类似于把单个组分场拆分为很多不同的场来计算,但是输入相同的源项,所以会有两种不同的Y在同一个方程里。程序就是在方程这一步出现了错误,我试试各位老师的方法,多输出一些内容看看哪里出问题了。