Large-eddy simulation of highly underexpanded transient gas jets
这个文献是OpenFOAM做的ILES

@李东岳 刚看了您分享的那个帖子 “在runTime++之后，如果T没有进行TEqn.solve()之类的更新，fvc:ddt(T) = 0”
那写这一项就没用了 搞不懂

@李东岳 没有乘号 就是个初始化

在看一个lusgs算法的代码
其中有一行
dRho( -dt* ( fvc::ddt(rho) + fvc::div(dbnsFlux.rhoFlux()) ) );
我理解的fvc是利用当前时间步显式计算 然后返回一个场
那这里的ddt(rho)是怎么计算的呢
求教~

@evensun 有些算例效果还是不太好 这里传不了文件了 你要感兴趣可以留个邮箱 我发给你

@李东岳 确实~ 现在改了改代码 能跑完算例 但是结果还是有点问题 等我完全搞定了 把代码传上来

@starDust boundaryFieldRef()是非const的可以赋值

@Steven_Yuan 同问

@bestucan 找原作者拿到了 ` 谢谢回复

尝试写了一个AUSM+up通量格式的求解器 编译的时候没有报错 但是运行forwardStep算例算两步就发散了
请大神帮我看一下代码 谢谢AUSMFoam.zip