buoyantPimpleFoam相关问题
-
@李东岳 东子老师，我今天花了一整天的时间仔仔细细的看完了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吧？
-
@李东岳 老师，算法文字部分我看懂了，然后我去看源代码，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循环到最后一步时，计算通量的嘛？
红框⑥
我感觉是在计算速度？
红框⑦
我感觉是在更新速度？
-
@李东岳 老师，我看完icoFoam之后，花了一个周末又重新研究了下buoyantPimpleFoam，又有一些新的疑问。
问题1
方程(20)、(24)、(26)中$U_{p}^{t}$前面的系数应该是$\rho_{p}^{\ast}/\Delta t$吧？
问题2
这一部分应该是速度方程吧？
问题3
前文未出现$U_{f}^{\ast \ast}$和$p^{\ast}$，我觉得应该改成$U_{p}^{'}$和$p_{rgh,f}^{'}$比较好理解吧？
问题4
这里应该是(37)吧?
问题5
红框这俩符号应该是一个物理量吧？
问题6
红框部分应该是$p_{rgh}^{\ast}$吧？
问题7
2->3->4->5->6->7->2应该就是一个时间步内的PISO循环吧？
问题8
步骤5中，方程(10)中只更新$U_{f}^{t}$为$U_{f}^{\ast \ast}$，不更新$\rho_{f}^{t}$为$\rho_{f}^{\ast \ast}$吗？
问题9
您在笔记中提到的phiHbyA=$ \sum U_{f} \cdot S_{f}$，我总觉得这里不对，我觉得应该是phiHbyA=$\sum HbyA_{f} \cdot S_{f}$。因为在不可压+瞬态这篇算法里面，您说fvc::div(phiHbyA)表示的是$\bigtriangledown \cdot HbyA^{t+\Delta t}$，不然就感觉对不上了。
问题10
pEqn.H文件中
const volScalarField psip0(psi*p)，这句代码里面的psip0和psi啥意思呀？p我知道是压力。因为后面出现了很多语句都有psi，不知道的话影响读代码了。问题11
笔记中这三个红框没看懂，前俩为啥是≈，而第三个却是=了？我觉得三个都应该是=呀。
