最近在写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 好的好的，十分感谢，我再去研究一下

请问在snappyhex中如何调整参数准确解析结构的边缘，比如下图圆柱端的直角边，以及翼型的尾缘。
5be2513d-fced-459a-aca8-577580c3a9f3-image.png
bbe3744d-7c04-487c-9df3-391365dfbc41-image.png
结构抓取的就非常的难看。我看其他网站也有类似的介绍。
024d5fa1-65db-49b4-b29e-16a71a702bb7-image.png
但是好像改了参数效果还是不好。各位有什么建议？目测应该不是继续加密的问题。

想用wave2foam造波浪+流，但是一直造不出来

@heike256 我对这个东西不太了解，目前也是在慢慢看，谢谢大佬的推荐，我研究研究。

@知行

./Allrun > /dev/null 2>&1

@李东岳 忘了说了，网格都是结构化的，长方体网格:duang:

@李东岳 再请教您一个问题，OpenFOAM中的DPM粒子追踪功能，在运动网格(ALE)上的支持并行计算吗？有没有什么限制？

OpenFOAM中的DPM粒子追踪功能，在运动网格(ALE)上的支持并行计算吗？有没有什么限制？有大佬测试过吗

个人感觉不一定是网格划分的问题，有可能是显示的问题

@学流体的小明 非常感谢，这个问题后面我自己解决了，解决的方法跟您的几乎一样，确实可以通过这种方式得到云图。

@Tens 好的！谢谢您！