OpenFOAM中的方程类为fvMatrix,该类为一个基于有限容积的稀疏矩阵类,更确切的说它是一个类模板。对于定义标量方程 fvMatrix<scalar>,也可以写成fvScalarMatrix,两者一样,因为openfoam运用了 typedef进行了类型别名定义
typedef fvMatrix<scalar> fvScalarMatrix ; //标量
typedef fvMatrix<vector> fvVectorMatrix ; //矢量
fvMatrix中有两个基本变量A 和 b初始化为0,用以存储系数矩阵A和右边值b。
Z
Best posts made by zwl
-
RE: fvVectorMatrix与fvScalarMatrix有什么区别?
Latest posts made by zwl
-
paraview数据后处理
请问paraview能够统计网格数吗?比如液相体积分数在0.5-0.6的网格数有多少。或者其它软件或者程序能够计算出来的方法,有了解的大神指点下,谢谢。
-
RE: reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界
是的,求解器自带的是蒸发,我想计算空化,二者的计算原理大不相同,我更改了质量源相,导致质量分数越界。
-
RE: reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界
@东岳 我看MULES代码,感觉他就是为求解相方程定制的,没有扩散相的处理,所以我不知道这里的扩散相放在哪个位置,怎么处理。
-
RE: reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界
@东岳 那么下面带拉普拉斯项的组分输运方程如何采用MULES求解呢?谢谢。
fvm::ddt(alpha, rho, Yi) + fvm::div(alphaRhoPhi, Yi, "div(" + alphaRhoPhi.name() + ",Yi)") - fvm::laplacian ( fvc::interpolate(alpha) *fvc::interpolate(this->muEff()/Sct_), Yi ) == alpha*this->R(Yi) + Source + fvc::ddt(residualAlpha_*rho, Yi) - fvm::ddt(residualAlpha_*rho, Yi)
-
RE: reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界
@yfclark 非常感谢你的建议。源相我确实是直接添加的,没有再做处理,这个Sp会有什么作用(虽然有时我也会模仿着这样处理)?越界的主要原因是什么?