CFD中文网

@Exthan 前辈您好，我最近也使用openfoam做超临界二氧化碳流动换热这块，和您出现了类似的问题，计算时连续性那块总是很差也很难收敛。想请教下对于类似的问题您目前有啥好的解决方法嘛

@王慧博 大佬您好，我最近也使用openfoam做超临界二氧化碳流动换热这块，和您出现了类似的问题，计算时连续性那块总是很差也很难收敛。然后我也尝试过使用定物性，计算则能很快的收敛，所以想向您请教下是不是buoyantPimpleFoam求解器本省对这种变物性问题不能很好的求解呢，或者目前有什么方法可以提高稳定性呢。

哦 我知道了 我用了nparticle

@星星星星晴 我看了，但是orgid编号怎么才到1000+，而且很多颗粒orgid是一样的

@星星星星晴 对的，我发现了这个id， ，我去试试，谢谢你

@星星星星晴 我在一个颗粒的时候就直接Info输出了，因为有的时候需要的变量在不同的函数里面，比如对流换热功率，表面反应的功率，对流换热计算在ThermoParcel.c里面，直接在计算换热的时候info输出一下就行，但是很多颗粒的时候，每个颗粒都调用一下，直接输出是不行的，请问有什么好方法吗

求助，使用了lagrangian库，算例中有很多颗粒，但是我想选择一个代表性的颗粒，跟踪这个一个颗粒，得到它质量，温度等这些属性随时间的变化情况，有什么好办法吗，我尝试过foamToVTK再paraview中去后处理，但是发现颗粒的id好像不是恒定的

@xr 如果你想继续研究的话，我们可以加一下联系方式，一起研究一下

@xr 我昨天又重新看了一下代码，发现Yi并不是质量分数，而是体积分数，所以如果你直接在下面扩散方程里面加扩散系数，应该是有问题的

@xr 你对 https://cfd-china.com/topic/6398/关于icoreactingmultiphaseinterfoam中的组分方程 怎么看

为什么icoReactingMultiphaseInterFoam组分方程求解的是 而不是
质量守恒

不可压缩把密度去掉 难道不是关于aU的方程吗

@xr 你说的好像对 我说算出来怎么扩散不明显

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

最后一个参数如果是true，

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

组分扩散在最下面啊
在constant里面有是否打开扩散的开关

@李东岳 东岳老师 ，我看了您关于MULES的解析，MULES的主要就是用两套通量去构建一个既保持高精度又不越界的通量，那具体要调用MULES::explicitSolve()要怎么输入参数呢，如果我要求解$\frac{\partial \alpha}{\partial t}+\nabla_{\bullet}\left(\alpha \vec{U}\right)+\nabla\cdot(\alpha_1\alpha_2c|\mathbf{U}|\frac{\nabla\alpha}{|\nabla\alpha|})=0$

template<class RhoType, class SpType, class SuType>
void Foam::MULES::explicitSolve
(
    const RhoType& rho,
    volScalarField& psi,
    const surfaceScalarField& phiPsi,
    const SpType& Sp,
    const SuType& Su
)

我在第三个参数phipsi要传入$\nabla_{\bullet}\left(\alpha \vec{U}\right)_{high}+\nabla\cdot(\alpha_1\alpha_2c|\mathbf{U}|\frac{\nabla\alpha}{|\nabla\alpha|})+\nabla_{\bullet}\left(\alpha \vec{U}\right)_{low}$这个总的通量吗

@李东岳 东岳老师，那加这个迎风的通量目的是为了什么，在icoReactingMultiphaseFoam里面前面是一个alpha的通量flux::fvc(alpha,phi)的，后面又加了一个 ，原作者注释掉了upwind<scalar>(mesh_, phi).flux(alpha1)，实际上alpha1Eqn.flux的结果和upwind<scalar>(mesh_, phi).flux(alpha1)是一样的，同样在compressmultiphaseinterfoam里面也加上了这个通量，而在intermixingFoam里面是先减去在加上这个

还是就是Sf * Uf * alpha

请问一下这个代码是什么意思，alpha插值到面上，用迎风格式吗？我看求解相方程里面都有这个，加这一个标量有什么意义吗

@Tens 我参考了icoreactingmultiphasefoam办法，采用mules求解组分方程，您看一下这样可行吗？
https://cfd-china.com/topic/6366/关于在compressibleinterfoam中添加组分方程

请问求解这方程理论上是可行的吗，推导似乎和alpha方程没什么区别