Skip to content

OpenFOAM

OpenFOAM交流区

5.3k 主题 31.1k 帖子
  • XiFoam中是否可以增加详细反应机理进行计算

    1
    1 帖子
    909 浏览

    诸位大佬,XiFoam中所使用的燃烧模型,是基于反应变量的湍流燃烧模型。
    b=1:未燃混合气b=0:已燃混合气
    是否可以考虑在其中增加详细的化学反应机理,进行计算。

  • PimpleFoam和PisoFoam对于槽道湍流计算的影响。

    3
    3 帖子
    2k 浏览
    2

    @学流体的小明 我的那个倒也没有发散,还在正常的计算,就是感觉要算的时间比pisoFoam要久。

  • 2 帖子
    1k 浏览
    bestucanB

    一个是锥子喷射,一个是锥子孔喷射。

    可能是:

    锥子喷射能保证喷着时是锥形分布。

    锥子孔空射,是保证出口形状是锥子孔,至于能不能喷射成锥形,得看喷射压力和喷射流量了。

    大口径小流量,是会单边喷射的,就像洗澡的花洒没水了。。。。:mianmo:

  • reactingFoam修改化学反应机理后遇到了问题

    5
    5 帖子
    5k 浏览
    Y

    又又又遇到了新问题,一步都没算下去,请问这种问题一般出在哪里:xinlei:
    Create time

    Create mesh for time = 0

    PIMPLE: No convergence criteria found

    PIMPLE: Operating solver in transient mode with 1 outer corrector
    PIMPLE: Operating solver in PISO mode

    Using LTS
    Reading thermophysical properties

    Selecting thermodynamics package
    {
    type heRhoThermo;
    mixture multiComponentMixture;
    transport sutherland;
    thermo janaf;
    energy sensibleEnthalpy;
    equationOfState perfectGas;
    specie specie;
    }

    Reading field U

    Reading/calculating face flux field phi

    Creating turbulence model.

    Selecting turbulence model type RAS
    Selecting RAS turbulence model kEpsilon
    kEpsilonCoeffs
    {
    Prt 0.85;
    Cmu 0.09;
    C1 1.44;
    C2 1.92;
    C3 0;
    sigmak 1;
    sigmaEps 1.3;
    }

    Creating thermophysical transport model

    Selecting thermophysical transport type RAS
    Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
    Creating reaction model

    Selecting combustion model EDC
    Selecting chemistry solver
    {
    solver ode;
    method standard;
    }

    standardChemistryModel: Number of species = 544 and reactions = 2446
    Selecting ODE solver seulex
    Creating field dpdt

    Creating field kinetic energy K

    No MRF models present

    No fvModels present
    No fvConstraints present

    Starting time loop

    Time scales min/max:
    Flow = 9.728569846e-08, 0.0001
    Temperature = 4.494232837e+307, 4.494232837e+307
    Composition = 4.494232837e+307, 4.494232837e+307
    Overall = 9.728569846e-08, 0.0001
    Time = 1e-07

    diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.0001092756986, No Iterations 1
    DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.006862564884, No Iterations 1
    DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.0001167898495, No Iterations 1
    #0 Foam::error::printStack(Foam::Ostream&) at ??:?
    #1 Foam::sigFpe::sigHandler(int) at ??:?
    #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #3 double Foam::standardChemistryModel<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGasFoam::specie >, Foam::sensibleEnthalpy> > >::solve<Foam::Field<double> >(Foam::Field<double> const&) at ??:?
    #4 Foam::combustionModels::EDC::correct() at ??:?
    #5 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
    #6 ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #8 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
    Floating point exception (core dumped)

  • reactingFoam报错

    5
    5 帖子
    3k 浏览
    vbcwlV

    @李东岳 谢谢东岳老师

  • 14 帖子
    8k 浏览
    X

    @lrl3512 要不您加一下我qq 2559727310

  • 气泡槽道流的计算,求指点

    11
    11 帖子
    6k 浏览
    学流体的小明

    更新一下进度
    暂时不用可压缩求解器计算了,因为发现当前的网格算不到声波那种程度。不过我倒是实现了向温度方程内添加源项这件事,就是不太确定公式以及编程对不对。

  • 在fsifoam中运行算例报错

    8
    8 帖子
    5k 浏览
    洱聿

    目前算例移植到计算机集群上计算,模型是水面上漂浮的平板。大约计算到0.8s时出现这个报错,查找原因是内存分配问题;
    8fda3931-8566-4c1d-96e3-08c9fb3e932f-image.png
    有查看了相关的log文件,文件中给出的报错是如下所示
    d5c39409-328d-471d-be7b-801f479d1faa-image.png
    为什么算例运行中会出现nu不被定义的报错呢?是因为计算过程中,网格出现问题嘛?

  • 7 帖子
    4k 浏览
    L

    @lrl3512 自问自答一下,在MULES::explicitSolve()之前调用了

    73804e01-261c-455f-9ddd-b16e1b828204-image.png

    最后一个参数如果是true,

    e99b54b9-951f-4d9b-b586-7be449bc10e7-image.png

    返回是通量的修正值,后面还是要加一个一阶迎风的通量,其实直接false也行,后面就不需要加了

  • 尝试加入神经网络推理,编译成功运行报错

    2
    2 帖子
    1k 浏览
    R

    @Rachel0096 接口部分内容如下:5252f072-1693-4abd-b204-3e6712b25d77-image.png 674fe43f-b3e0-49ad-a2de-daf068772bc0-image.png

  • Paraview颜色问题

    5
    5 帖子
    5k 浏览

    @bestucan 向老师学习了

  • openFoam如何调用热物理库中的气体

    3
    3 帖子
    2k 浏览
    Y

    @vbcwl 请问你调用成功了嘛,我也在尝试firefoam的这里,我想改成正庚烷的

  • 关于浸没边界法

    2
    2 帖子
    1k 浏览
    T

    没有人知道吗?

  • 安装openfoamv1812编译出现问题

    3
    3 帖子
    2k 浏览
    X

    @李东岳 好的 谢谢李老师的回复,我试了一下得换成Ubuntu18版本才可以,22版本就是不可以

  • paraview输出动画或者图片的过程中途报错

    11
    11 帖子
    7k 浏览
    bestucanB

    @dxl 不用放时刻文件夹。文件名中的001 002 003就被当作顺序了

  • paraview如何放大默认的坐标轴和颜色标尺?

    5
    5 帖子
    4k 浏览

    自问自答了

    放大颜色标尺方法:找到下面这个图标,点进去就可以自由发挥了
    0e9ab1e7-51c1-4809-be40-9a0cecb0a87c-1.JPG

  • 请问ReactingFoam可否用于气液两相燃烧计算

    3
    3 帖子
    2k 浏览
    Y

    @bestucan 好的好的,十分感谢!

  • OpenFOAM中的chemistry solve and EulerSI/Stiff ODE

    1
    1 帖子
    944 浏览
    李东岳

    最近在写chemFoam解析。在这里做个记录。

    第一层:在solver求解器层面,输入时间步长dtChem = chemistry.solve(runTime.deltaT().value())并更新。

    第二层:在StandardChemistryModel.C里面,进入solve(c_, Ti, pi, dt, this->deltaTChem_[celli])

    第三层:在ode.C,进入odeSolver_->solve(0, deltaT, cTp_, subDeltaT)

    第四层:ODESolve.C,进入solve(x, y, step)

    第五层:进入某个具体的ODE求解器,比如EulerSI.C,进入adaptiveSolver::solve(odes_, x, y, dxTry)

    第六层:回到EulerSI.C,进入

    Foam::scalar Foam::EulerSI::solve ( const scalar x0, const scalarField& y0, const scalarField& dydx0, const scalar dx, scalarField& y ) const

    这一层是最终的ODE求解步骤。

  • 2维网格拉伸为3维后LES计算报错

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

    @乌云盖雪
    你的k文件写错了,下面是一个我随便找的示例文件

    /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object k; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #include "include/initialConditions" dimensions [0 2 -2 0 0 0 0]; internalField uniform $turbulentKE; boundaryField { // 你的k文件这里是圆括号…… #includeEtc "caseDicts/setConstraintTypes" //- Define inlet conditions #include "include/fixedInlet" outlet { type inletOutlet; inletValue $internalField; value $internalField; } lowerWall { type kqRWallFunction; value $internalField; } motorBikeGroup { type kqRWallFunction; value $internalField; } #include "include/frontBackUpperPatches" } // ************************************************************************* //
  • OpenFOAM冷凝模型在哪

    2
    2 帖子
    1k 浏览
    bestucanB

    https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1compressible_1_1humidityTemperatureCoupledMixedFvPatchScalarField.html

    https://www.cfd-online.com/Forums/openfoam-solving/147129-condensation-using-openfoam.html

    https://www.researchgate.net/publication/279643607_Modeling_of_Direct_Contact_Condensation_With_OpenFOAM

    https://jhmtr.semnan.ac.ir/article_467.html

    https://github.com/NimaSam/phaseChangeHeatFoam