buoyantPimpleFoam相关问题

[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"
        }

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

[李东岳]

我看了你的回帖，我认为我写的跟代码都能对得上。我觉得你需要从icofoam看起，再来看这个。现在理解这个不好理解。

[yjc95]

@李东岳 好的，谢谢东子老师

[yjc95]

@李东岳 东子老师，我今天花了一整天的时间仔仔细细的看完了icoFoam（CFD: 不可压 + 瞬态）的算法，在看源代码之前有几个疑问。

问题1

红框部分应该是$U_{f}^{*}$吧？

问题2

无痛苦N-S笔记方程(5.138)我有点疑问，从红框部分接着计算，

分子部分

所以最后的结果应该为
$\frac{\phi _{WW}+\phi _{EE}-2\phi _P}{2\Delta x}\Delta y\Delta z$
所以分母系数应该是2吧？

问题3

红框部分是咋来的呀？

问题4

红框部分是否应该为$U_{p}^{t}$？

问题5

应该是方程(21)相对于方程(20)存在一定的滞后吧？

问题6

方程(22)中，
$U_{p}^{'}= U_{p}^{t+\Delta t}- U_{p}^{\ast}$，
方程(29)处
$U_{p}^{\ast\ast}= U_{p}^{'}+ U_{p}^{*}$，您是想表达迭代两步之后就收敛了，所以$ U_{p}^{t+\Delta t}= U_{p}^{**}$？

问题7

这里是(33)吧？

问题8

红框应该是PIMPLE吧？

[李东岳]

问题1，Up, Uf都是未知

问题3，就是梯度的常规计算方法

问题6，U**不是U$^{t+\dt}$，29不是最终跳出时间步的方程，还要继续迭代

问题8，是SIMPLE。PIMPLE=SIMPLE+PISO

问题2，你说的对，已更正
问题4，你说的对，已更正
问题5，你说的对，已更正
问题7，你说得对，已更正

[yjc95]

@李东岳 老师，算法文字部分我看懂了，然后我去看源代码，icoFoam文件夹下面主要有icoFoam.C和creatFields.H俩文件，我只看了icoFoam.C文件，有几个疑问，都在PISO循环部分。

红框①

您展示的关键代码没有这行代码，感觉像是定义HbyA，constrain啥意思呀？括号里U，P有啥作用？

红框②

我看您在无痛N-S方程里面提到过phiHbyA=$\sum U_{f}\cdot S_{f}$，明显phiHbyA不等于HbyA，但您在笔记中又提到了“同时需要注意，在组建phiHbyA的时候，可以通过fvc::flux(HbyA)函数调用不同的离散格式”，意思是fvc::div(phiHbyA)=$\sum HbyA_{f}^{t+\Delta t}\cdot S_{f}$了？因为您又提到了“可以把HbyA理解为速度中间量”，不然图片中间while语句部分的fvc::div(phiHbyA)和您文章中的方程(26)就对不上了。如果真是这样的话，红框②中fvc::flux(HbyA)代表啥呀？fvc::interpolate(rAU)是插值到面上对吗？ fvc::ddtCorr(U, phi)在笔记中没有搜到。

红框③和④

这俩语句作用是啥？

红框⑤

我感觉这两行代码是在PISO循环到最后一步时，计算通量的嘛？

红框⑥

我感觉是在计算速度？

红框⑦

我感觉是在更新速度？