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

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

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

最后一个参数如果是true，

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

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

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

@bestucan 向老师学习了

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

没有人知道吗？

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

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

自问自答了

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

@bestucan 好的好的，十分感谢！

最近在写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求解步骤。

@乌云盖雪
你的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" } // ************************************************************************* //

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

可能是网格太大，内存不够用。换个文件尺寸小点的网格文件试试

时间步太大了，库郎数（Courant Number）都飙到三百多了。按说应该小于1

减小时间步，或者降低流速

就是在最大迭代步之前没有达到设定的收敛判据。但是因为达到限制，所以停止迭代。

瞬态计算中这种设置挺常见的。多少个时间步，每个时间步最多迭代多少次。

理想的当然是在达到最大迭代之前，在当前时间步收敛。
但是如果一直不收敛，也不能卡在这个时间步无限迭代，所以要设置最大迭代。

但是最大迭代步数内是不是收敛受多个因素影响：时间步大小（越小越容易收敛）、收敛判据（越大越容易判定收敛）

所以，只是个“警告”，不是个“错误”。因为“达到最大还没收敛”这件事情意味着什么，要结合其他的东西来看待。

Linux里用

paraview --mesa

就是用 CPU 代替 GPU。

win 上也可以在cmd里给paraview.exe 加参数运行。
或者创建快捷方式，在快捷方式属性里加参数

看很多文献中都是将计算域的长、宽、高和二分之一的高度$\delta$做比值等于多少倍的$\pi$，这其中的原理是什么？

@香柏树 请问你用ReactingFoam ras改les成功了吗，想和您请教一下

@myheart 好的好的，十分感谢，我再去研究一下