buoyantPimpleFoam相关问题

[李东岳]

buoyantPimpleFoam解析

你看的是哪个链接？网站已经不存在这个链接了

[yjc95]

@李东岳 http://dyfluid.com/buoyantPimpleFoam.html
在这个链接里

[李东岳]

那个是旧的链接，写的不完善。目前我得网页里面不提供公开访问权限，你找到的应该是搜到的，我先删掉了。你看看新版本的。

http://dyfluid.com/compsimple.html

http://dyfluid.com/compiso.html

[yjc95]

@李东岳 谢谢东子老师

[yjc95]

@李东岳 东子老师，我把可压+稳态算法仔仔细细的看完了，大体上看懂了，但有一些疑问向您请教。您上次说我排版不利于阅读，这次我专门看了您提供的帮助里面的教程。

问题1

您在文章里说：上标n表示为当前迭代步（已知），上标表示下一个迭代步（待求），n+1也表示下一个迭代步。在稳态条件下， 应该等于 吧。方程(18)-(14)等号左边不应该是0吗？为啥是 呢？上标、**、n与n+1是啥关系呀？我猜测*、**是不是n到n+1之间的迭代步？

问题2

方程(23)没有 呀？另外，为啥是 而不是 或者 呢？

问题3

由方程(3)可得 ，又 ，所以
但由方程(27)得到的却是图中红框部分，如何得出方程(28)中蓝框部分的呢？方程(36)是我算的结果，方程(42)又是红框的结果了。

问题4

方程(19)由(18)-(14)得来， ，而图中这个推导过程：
是 ，是不是矛盾了呢？

同样的，方程(43)中， ；
方程(25)中，
我被这些上标n、'、*搞懵了

问题5

方程(48)中，红框部分是不是少了一个 ?

[李东岳]

我昨天凑巧更新了一下，http://dyfluid.com/compsimple.html 你再看看

打公式的时候，用美元符号最好 $\alpha$

[yjc95]

@李东岳 谢谢东子老师

[yjc95]

@李东岳 东子老师，我看完了您更新的文章，还有几个小问题。另外，有的时候我用 $\alpha$打不出来我想要的公式，可能前后公式会不统一，请见谅。

问题1

方程(4)为啥把速度处理成 ${\mathbf{\text{U}}}^n{\mathbf{\text{U}}}^{\ast }$呢？为啥不是${\mathbf{\text{U}}}^n{\mathbf{\text{U}}}^n$或者${\mathbf{\text{U}}}^{\ast }{\mathbf{\text{U}}}^{\ast }$呢？

问题2

方程(25)密度为啥写成${\rho }_f^n+{\rho }_f^{{}^{\prime }}$,为啥要修正呢？速度为啥用${\mathbf{\text{U}}}_f^{\ast \ast }$呢？用${\mathbf{\text{U}}}_f^{\ast }$不行吗？

问题3

方程(15)处，您说更新热物理模型库中的密度为${\rho }^{\ast }=p^n/{\psi }^n$，这和方程(3)不一致呀？方程(43)处您又说是根据状态方程更新密度中间量。俩方程长的一样，一个是热物理模型库，一个是状态方程？

问题4

方程(42)中，您是用${\rho }_f^{{}^{\prime }}$来替代方程(41)中的$▽\rho$吗？

问题5

方程(48)写成(49)，等号右边这一项为啥可以用
来代替$▽{\rho }_{rhg}^{\ast }$呢？

[李东岳]

问题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方程笔记

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

[李东岳]

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

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

divPhiHbyA是方程的第一项

[yjc95]

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

[李东岳]

可以忽略