buoyantPimpleFoam相关问题

[yjc95]

@李东岳
问题3，状态方程不是方程(3)$p=\rho RT$吗？$\psi=1/(RT)$那么不应该$\rho =p/(RT)=p\psi$吗？为什么是$\rho =p/\psi$呢？

[李东岳]

dyfluid上面应该写错了，明天我核实一下告诉你

[李东岳]

rho=p \psi 这个是对的，已更正

[yjc95]

@李东岳 老师，好像没改过来

[李东岳]

你刷新一下看看

[yjc95]

@李东岳 刷新了，还是这样

[李东岳]

我想起了 ，我忘记上传了... 现在上传了

[yjc95]

@李东岳 老师，前面算法理论部分我看懂了，但我在看关键代码这一部分的时候，发现和前面公式有点对应不上。比如高速流transonic中的压力方程通过下面的代码组建：

fvc、fvm是啥意思啊？我查了下无痛苦N-S方程笔记

为啥这里要除以$\bigtriangleup V\bigtriangleup t $ 呢？
另外，$fvc::div(phiHbyA)$ 是对应方程中的哪一项呢？代码中每一项代表方程中的哪一项呢？这个我不是很理解。

[李东岳]

fvc是显性离散，fvm是隐性离散。

除掉时间和体积就是因为代码本身就是这样做的

divPhiHbyA是方程的第一项

[yjc95]

@李东岳 代码等号右边那项是啥呀？

[李东岳]

可以忽略

[yjc95]

@李东岳 东子老师，我看了可压+稳态的算法之后，紧接着看了可压+瞬态的算法，有几个问题想请教您。

问题1

方程(8)等号左边第一项离散后是$ \int (\int \frac{\partial \rho }{\partial t}dt)dV=\bigtriangleup V(\rho ^{t+\bigtriangleup t}-\rho ^{t})$
第二项离散后是
$\int (\int \bigtriangledown \cdot \rho \textbf{U}dV)dt=$ $\bigtriangleup t$ $\sum_{f} \rho _{f}^{t+\bigtriangleup t}$ $\textbf{U} _{f}^{t+\bigtriangleup t}$ $\cdot S$
这俩公式相加化简后应该是
为啥方程(9)左边第一项没有$\bigtriangleup V$呢？而方程(11)又有$V{p}$呢？

问题2

方程(19)左边第一项应该是吧？
另外，今晚回复帖子打字好卡

[李东岳]

问题1你说的是对的，对流项忘记除掉体积，已更新

问题2跟你写的一样吧？

[yjc95]

@李东岳 左括号应该在求和号左边

另外，东子老师，我看代码有点吃力，我需要恶补C++编程吗？之前只接触过python和C#

[李东岳]

是要补充C++编程，但是感觉也没有太多需要补的。编程语言基本都一样，坎习惯了就好

问题2已更新，你说得对

[yjc95]

@李东岳 老师，应该这样吧

[李东岳]

对。。已更新

[yjc95]

@李东岳 东子老师，您在dyfluid上展示的代码是of哪个版本的呀？我打开of8下面的buoyantPimpleFoam.C文件，发现它的顺序和这个求解顺序不一样呀，然后就看不懂了

[李东岳]

怎么不一样

[yjc95]

@李东岳 我按步骤顺序来说我看到的有啥不一样哈，为了方便比较，我就用截图的方式了。

步骤1

#include "rhoEqn.H"来求解$\rho^{*}$
of8中buoyantPimpleFoam.C头文件并没有和它类似的

步骤2

组建速度方程在buoyantPimpleFoam.C文件中也没有找到，反而是在该文件夹下的UEqn.H文件里面找打了

步骤3

用thermo.correct()更新$\psi^{*}$,用rho = thermo.rho()更新$\rho^{\ast}$
前者我是在EEqn.H中找到的，这是个啥文件呀？后者我是在pEqn.H里面找到的。

步骤4

用

fvScalarMatrix p_rghDDtEqn
(  
    psi*correction(fvm::ddt(p_rgh))
  + fvc::ddt(rho) 
  + fvc::div(phi)
  - fvm::laplacian(rhorAUf, p_rgh)
);

组建压力柏松方程，我在pEqn.H文件中找到了一个类似的，但长得不一样

步骤5

用

phi += p_rghEqn.flux()

来更新速度，我在pEqn.H文件中找到了一个类似的，但长得不一样

步骤6

我在pEqn.H文件中找到了密度更新

步骤7

我在pEqn.H文件中找到了判断连续性误差

为啥会分散到这么多文件中呀？除了您展示的关键代码，其他部分例如

if (correctPhi)
        {
            divrhoU = new volScalarField
            (
                "divrhoU",
                fvc::div(fvc::absolute(phi, rho, U))
            );
        }

        if (LTS)
        {
            #include "setRDeltaT.H"
        }
        else
        {
            #include "compressibleCourantNo.H"
            #include "setDeltaT.H"
        }

确实看的一脸懵圈，心态处于将崩不崩的边缘