关于buoyantPimpleFoam的几个问题
-
http://dyfluid.com/rhoSimpleFoam.html http://dyfluid.com/simplefoam.html
我重写了,尤其是rhoSimpleFoam解析,尤其是说明了rho=thermo.rho(),写的比较仓促,你先对一下看看能不能自我回答你的问题,我要出去一趟,过几天把瞬态的更新下,同时更新下thermo,.psi之类,需要很多工作 -
http://dyfluid.com/rhopimplefoam.html 我更新了下rhoPimpleFoam解析,这个跟buoyantPimpleFoam非常接近了,我下周处理buoyantPimpleFoam
-
@东岳 我回来了!不好意思隔了快一周了。看了rhoSimpleFoam和rhoPimpleFoam,目前还有以下几个问题希望和大家一起探讨。
-
高版本的rhoPimpleFoam既可以调用psiThermo也可以调用rhoThermo,而低版本调用的是psiThermo。给人的感觉有点混乱,这两个热物性的区别到底是什么。换句话说,psiThermo和rhoThermo是否可以随意使用,还是某些求解器只适合调用某一种。如果我想自己组装一个新的求解器怎么判断选择psiThermo还是rhoThermo呢?
-
还是关于密度方程的("rhoEqn.H")。我的理解是求解密度方程的原因是为了获得连续性误差(pEqn.H的末尾#include "rhoEqn.H"和#include "compressibleContinuityErrs.H")。如果仅仅是这个原因,为什么每一个时间步的开头要求解密度呢。(对于rhoPimpleFoam,如果PISO修正步数小于等于1则调用"rhoEqn.H";对于buoyantPimpleFoam则没有判定,一定要调用。)这一点其实还好,因为如果每一个时间步严格收敛了,那密度也不会再变化了,求解器密度和热物性密度应当完全一样。
-
综合我看到的关于buoyantPimpleFoam的解析,它适用于求解浮升力引起的温度变化而导致的密度变化。这个是不是类似于所谓的低马赫数假设(密度不随压力变化,而随温度、组分等变化)?如果是(或者不是),原因是什么?
-
-
@cccrrryyy 第二个我知道,求解密度方程是为了求解质量守恒方程,连续性残差是密度残差的反映,具体什么关系没研究过,当做相等也没什么,顶多有个常系数关系。
质量守恒方程所运输的物理量就是密度,所以一直在围着密度转,所以时间步开头就算密度,然后通过迭代减小残差。
-
@cccrrryyy 惭愧惭愧,您的理解比我深入多了,
我再猜测一下,可能是求解动量守恒方程的时候(夹在物质守恒和能量守恒中间?)要用密度。我记得求解器里的压力是除过密度的压力。这个时候开头算不算密度应该会有影响?
好了,肚里没货了,再和您讨论要明天充充电了
-
http://dyfluid.com/rhoPimpleFoam.html 公式5的描述解释了为什么要rhoEqn.H。至于Pimple多次迭代,是因为肯定会更新
,所以没必要开始rhoEqn.H我昨天想更新buoyantPimpleFoam,发现这俩差不多就没写
不过我整理下heRhoThermo与hePsiThermo然后回复你 -
@东岳 有点理解了!buoyantPimpleFoam和rhoPimpleFoam确实差不多,我这几天想着拿一个小算例用这两个求解器试一下,感觉很有可能会解出来不一样的结果。目前我的猜测是buoyantPimpleFoam用了rhoThermo,而rhoPimpleFoam(低版本)用了psiThermo,用rhoThermo的话求解压力方程的时候密度是已知的,而用psiThermo的话求解压力方程的时候密度是未知的(或者说 psi*p 这一项 p是未知的),但不知道这样会对求解造成什么影响。
东岳老师太辛苦了,是我辈楷模!
-
@cccrrryyy 在 关于buoyantPimpleFoam的几个问题 中说:
做过将buoyantPimpleFoam修改为一个低马赫数求解器(变密度求解器
前辈您好,请问一下您说的将buoyantPimpleFoam写改为低马赫求解器的工作在哪里可以看到?谢谢
-
@lyc https://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2018/WillHay/lowMachBuoyantPimpleFoam_Corrected.pdf
不是论文,是他们学校的课程作业,Chalmers有专门的OpenFOAM课
13/19