buoyantPimpleFoam相关问题

[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版本不同吗？

[李东岳]

是的

[yjc95]

@李东岳 东子老师，根据公式(36)，红色箭头部分应该缺了$1/V_p$吧。

[李东岳]

你说的对，已加上体积，已更新

[yjc95]

@李东岳 老师，这俩没加上体积