关于buoyantPimpleFoam的几个问题

[cccrrryyy]

@bestucan 感谢回复！我理解连续性残差其实就是在算求解器得到的密度（质量守恒方程）和热物性更新得到的密度（状态方程）之间的差别，但这个步骤在每次求解完压力并更新速度之后（即pEqn.H的末尾）都已经进行了，所以不太理解为什么每个时间步的开头要计算一次密度。

[bestucan]

@cccrrryyy 惭愧惭愧，您的理解比我深入多了，

我再猜测一下，可能是求解动量守恒方程的时候（夹在物质守恒和能量守恒中间？）要用密度。我记得求解器里的压力是除过密度的压力。这个时候开头算不算密度应该会有影响？

好了，肚里没货了，再和您讨论要明天充充电了

[李东岳]

http://dyfluid.com/rhoPimpleFoam.html 公式5的描述解释了为什么要rhoEqn.H。至于Pimple多次迭代，是因为肯定会更新$A_\rP$，所以没必要开始rhoEqn.H

我昨天想更新buoyantPimpleFoam，发现这俩差不多就没写
不过我整理下heRhoThermo与hePsiThermo然后回复你

[cccrrryyy]

@bestucan 每天都在充电学习，谈不上理解深入或者不深入啥的。

[cccrrryyy]

@东岳 有点理解了！buoyantPimpleFoam和rhoPimpleFoam确实差不多，我这几天想着拿一个小算例用这两个求解器试一下，感觉很有可能会解出来不一样的结果。目前我的猜测是buoyantPimpleFoam用了rhoThermo，而rhoPimpleFoam（低版本）用了psiThermo，用rhoThermo的话求解压力方程的时候密度是已知的，而用psiThermo的话求解压力方程的时候密度是未知的（或者说 psi*p 这一项 p是未知的），但不知道这样会对求解造成什么影响。

东岳老师太辛苦了，是我辈楷模！

[李东岳]

有点耽误 这周事情有点多 但我没忘这件事 当我有时间 我会更新 我先把一段代码写完

[cccrrryyy]

@东岳 你忙~我最近再看Ferziger和Peric那本书的Compressible Flow一章，了解的也更多了。看代码之余还是应该多温习理论！

[lyc]

@cccrrryyy 在 关于buoyantPimpleFoam的几个问题 中说：

做过将buoyantPimpleFoam修改为一个低马赫数求解器（变密度求解器

前辈您好，请问一下您说的将buoyantPimpleFoam写改为低马赫求解器的工作在哪里可以看到？谢谢

[cccrrryyy]

@lyc https://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2018/WillHay/lowMachBuoyantPimpleFoam_Corrected.pdf

不是论文，是他们学校的课程作业，Chalmers有专门的OpenFOAM课

[lyc]

@cccrrryyy谢谢前辈