buoyantPimpleFoam相关问题

[李东岳]

问题1，Up, Uf都是未知

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

问题6，U**不是U${}^{t+\mathrm{\Delta }t}$，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 Hby{A}_f^{t+\mathrm{\Delta }t}\cdot S_f$了？因为您又提到了“可以把HbyA理解为速度中间量”，不然图片中间while语句部分的fvc::div(phiHbyA)和您文章中的方程(26)就对不上了。如果真是这样的话，红框②中fvc::flux(HbyA)代表啥呀？fvc::interpolate(rAU)是插值到面上对吗？ fvc::ddtCorr(U, phi)在笔记中没有搜到。

红框③和④

这俩语句作用是啥？

红框⑤

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

红框⑥

我感觉是在计算速度？

红框⑦

我感觉是在更新速度？

[Amadeus]

从微信推送过来看的哈哈哈，东岳老师真的很耐心，楼主也非常细致! 超级棒

[yjc95]

@Amadeus 是呀

[yjc95]

@李东岳 东子老师，您网站上的虚拟机好像不支持win7系统，有支持win7系统的嘛，of版本不是最新的也行呀。因为我的电脑没法联网，我想下载一个离线版本，然后拷贝到win7的电脑里面用

[李东岳]

win7太老了，没计划支持win7了老铁

[yjc95]

@李东岳 东子老师，这几个代码的疑问您抽空看下哈

[李东岳]

7，更新边界条件
6，你说得对
5，更新通量
3、4、调整边界条件，确保边界条件设置的没问题
1、就是HbyA，constrain用于保证HbyA边界条件正确

[yjc95]

@李东岳 红框②中phiHbyA是为了更方便的计算$\sum HbyA\ast S_f$才这样弄的吗？

[李东岳]

fvc::flux(HbyA)就是phiHbyA

[yjc95]

@李东岳 那fvc::flux(HbyA)后面加上fvc::interpolate(rAU)*fvc::ddtCorr(U, phi)是啥作用呀？

[李东岳]

ddtCorr这个是Weller自己提出的防止小时间步长影响时间结果的方法，我明天看看更新到笔记

[yjc95]

@李东岳 老师，我看完icoFoam之后，花了一个周末又重新研究了下buoyantPimpleFoam，又有一些新的疑问。

问题1

方程(20)、(24)、(26)中$U_p^t$前面的系数应该是${\rho }_p^{\ast }/\mathrm{\Delta }t$吧？

问题2

这一部分应该是速度方程吧？

问题3

前文未出现$U_f^{\ast \ast }$和$p^{\ast }$，我觉得应该改成$U_p^{{}^{\prime }}$和$p_{rgh,f}^{{}^{\prime }}$比较好理解吧？

问题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 Hby{A}_f\cdot S_f$。因为在不可压+瞬态这篇算法里面，您说fvc::div(phiHbyA)表示的是$▽\cdot Hby{A}^{t+\mathrm{\Delta }t}$，不然就感觉对不上了。

问题10

pEqn.H文件中const volScalarField psip0(psi*p)，这句代码里面的psip0 和 psi啥意思呀？p我知道是压力。因为后面出现了很多语句都有 psi，不知道的话影响读代码了。

问题11

笔记中这三个红框没看懂，前俩为啥是≈，而第三个却是=了？我觉得三个都应该是=呀。

[yjc95]

@李东岳 东子老师，我真的不是在审稿，我是在认真问问题呀

[李东岳]

今天忙 事多没顾上 明天回你

[李东岳]

问题1 是可压瞬态算法么？没理解

问题2 已更新

问题4 已更新

问题5 是一个物理量

问题6 是修正压力

问题7 是的

问题8 不更新

问题9 笔记里面的截图我看一下

问题10 psi是可压缩性，psi0是更新前的可压缩性

问题11 这三个应该都是约等号，我更新一下

[yjc95]

@李东岳 问题1是可压瞬态算法里面的
问题9的图如下

[yjc95]

@李东岳 东子老师，下面这两张截图是来自可压+瞬态算法，其中：

psi*correction(fvm::ddt(p_rgh)) 对应红框1

fvc::ddt(rho) 对应红框2

那么fvc::div(phi) 对应方程(39)哪一项呢？您在可压+稳态算法中提到，可压求解器中通量phi定义为$\rho U\cdot S_f$，那代入fvc::div(phi) 这句代码中，不应该就是$\rho U\cdot S_f$吗？但方程(39)中跟它很像的只有红框3，只是把$U$替换成了$HbyA$？这样替换的原因是啥呀？

fvm::laplacian(rhorAUf, p_rgh) 我猜是对应红框4，那红框5是被省略了还是合并了？

我发现，我用of8版本，压力泊松方程和您给的不一样，出现了phiHbyA，这也是我在上一个问题问您的。

[李东岳]

方程的3跟5合起来是fvc::div(phiHbyA)

有关fvc::div(phiHbyA)是U还是HbyA，看我那个笔记，昨天更新了一下

[yjc95]

@李东岳 那fvc::div(phi)和of8中的fvc::div(phiHbyA)是代表一个意思，只不过是of版本不同吗？