buoyantPimpleFoam相关问题

[李东岳]

问题1这是要做线性化。

问题3，热物理模型库就是状态方程。$\psi=1/R/T$

问题4，一个是离散形式，一个是连续形式。离散形式的是修正量。

问题5，一个是离散形式，一个是连续形式。

[yjc95]

@李东岳
问题3，状态方程不是方程(3)$p=\rho RT$吗？$\psi=1/(RT)$那么不应该$\rho =p/(RT)=p\psi$吗？为什么是$\rho =p/\psi$呢？

[李东岳]

dyfluid上面应该写错了，明天我核实一下告诉你

[李东岳]

rho=p \psi 这个是对的，已更正

[yjc95]

@李东岳 老师，好像没改过来

[李东岳]

你刷新一下看看

[yjc95]

@李东岳 刷新了，还是这样

[李东岳]

我想起了 ，我忘记上传了... 现在上传了

[yjc95]

@李东岳 老师，前面算法理论部分我看懂了，但我在看关键代码这一部分的时候，发现和前面公式有点对应不上。比如高速流transonic中的压力方程通过下面的代码组建：

fvc、fvm是啥意思啊？我查了下无痛苦N-S方程笔记

为啥这里要除以$\bigtriangleup V\bigtriangleup t $ 呢？
另外，$fvc::div(phiHbyA)$ 是对应方程中的哪一项呢？代码中每一项代表方程中的哪一项呢？这个我不是很理解。

[李东岳]

fvc是显性离散，fvm是隐性离散。

除掉时间和体积就是因为代码本身就是这样做的

divPhiHbyA是方程的第一项

[yjc95]

@李东岳 代码等号右边那项是啥呀？

[李东岳]

可以忽略

[yjc95]

@李东岳 东子老师，我看了可压+稳态的算法之后，紧接着看了可压+瞬态的算法，有几个问题想请教您。

问题1

方程(8)等号左边第一项离散后是$ \int (\int \frac{\partial \rho }{\partial t}dt)dV=\bigtriangleup V(\rho ^{t+\bigtriangleup t}-\rho ^{t})$
第二项离散后是
$\int (\int \bigtriangledown \cdot \rho \textbf{U}dV)dt=$ $\bigtriangleup t$ $\sum_{f} \rho _{f}^{t+\bigtriangleup t}$ $\textbf{U} _{f}^{t+\bigtriangleup t}$ $\cdot S$
这俩公式相加化简后应该是
为啥方程(9)左边第一项没有$\bigtriangleup V$呢？而方程(11)又有$V{p}$呢？

问题2

方程(19)左边第一项应该是吧？
另外，今晚回复帖子打字好卡

[李东岳]

问题1你说的是对的，对流项忘记除掉体积，已更新

问题2跟你写的一样吧？

[yjc95]

@李东岳 左括号应该在求和号左边

另外，东子老师，我看代码有点吃力，我需要恶补C++编程吗？之前只接触过python和C#

[李东岳]

是要补充C++编程，但是感觉也没有太多需要补的。编程语言基本都一样，坎习惯了就好

问题2已更新，你说得对

[yjc95]

@李东岳 老师，应该这样吧

[李东岳]

对。。已更新

[yjc95]

@李东岳 东子老师，您在dyfluid上展示的代码是of哪个版本的呀？我打开of8下面的buoyantPimpleFoam.C文件，发现它的顺序和这个求解顺序不一样呀，然后就看不懂了

[李东岳]

怎么不一样