CFD中文网

@Caijinjin 多孔介质系列的模型本来就是基于表观速度的，如果需要考虑空隙率的的阻碍作用，可以考虑volumeFractionSource。流动也未必是湍流，没有湍流之后的效果可以把该区域强制设为层流区（把k用fixedValueConstaint强制取0）算一遍

使用常用的surfaceFieldValue监控工具无法实现，因为速度场不是面场，有没有办法实现速度场插值到指定的faceZone，或者能够直接监控（不包括修改求解器代码，比如求解器内新建一个面速度场然后写入数据）

    max
    {
        type            surfaceFieldValue;
        libs            ("libfieldFunctionObjects.so");
        regionType      faceZone;
        name           internal1;
        operation       maxMag;
        surfaceFormat   none;  
        fields
        (
            U
        );
    }  

激发态直接给对应的热力学数据值就行，要加E是因为保证电荷守恒。而且激发态的恒压热容和熵应该和原始态有区别吧

红圈的地方分开来，C之前是自定义的信息注释，O 2，代表元素，两个氧; C 1代表一个碳

vscode有时候会有这个问题，比如我明明开了自适应时间步长，但是OF还是按固定步长算的，实际的文件改动根本没同步到wsl。重开下vscode就行

破案了，是VSCode和wsl之间的bug，vscode的自动保存不生效，实际的fvScheme还是老版本，关闭保存再开也没用。重启vscode就没事了

发现一个奇怪的问题，我在算一个组分传递的Case，离散格式刚开始设置如下：

divSchemes
{
    default         none;
    div(phi,Yi_h)         Gauss upwind;
    div(phi,U)      Gauss  upwind;
    div(phi,Yi)   Gauss upwind;
    div(phi,h)      Gauss limitedLinear 1;
    div(phi,K)      Gauss upwind;
    div(phid,p)     Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(phi,k)      Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U)))))   Gauss linear;
}

算了几步基本就开始发散，有时返回矩阵奇异，但是我把div(phi,Yi_h)和div(phi,Yi)位置换了下就好了。有没有大佬解释下