divShemes离散格式的顺序问题

[han1112]

发现一个奇怪的问题，我在算一个组分传递的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)位置换了下就好了。有没有大佬解释下

[shirouchen]

我记得李老师之前的介绍里是说，后会覆盖前。但是你这div(phi,Yi_h)和div(phi,Yi)的离散格式不都是迎风吗，按理不会有区别吧。

[李东岳]

这个就很奇怪。你要是能准备一个小算例，我可以看看，如果是bug，可以提bug。但是这个连bug都不像。像是你还是改了什么，然后忘记了

[han1112]

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

[han1112]

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

[李东岳]