【讨论】高Ra 数，高Pr数自然对流模拟遇到的问题

yzgdy

自然对流模拟中，Ra=gβ△TH^3/(αν)，如Ra增大，温差参数设置不变，H不变，Pr数不变的情况，
（1）如热膨胀系数也不变，则一般调小粘性系数或热扩散系数，但这样松弛时间变小，为得到一定稳定状态或者某一无量纲时间点，计算量增大。如松弛时间过小，不稳定问题也可能出现
（2）如粘性系数或热扩散系数不变，即松弛时间系数不变，则增大热扩散系数，计算中则发现计算出感觉“非物理”的现象
（3）如其它参数不变，提高H,则带来计算量的飞速增大的问题
这些问题，请大家给些建议，解决方法或计算方法出现的错误或问题

[ 本帖最后由 yzgdy 于 2010-3-8 06:42 编辑 ]

yzgdy

另一个相关问题：高Pr 数的问题，对于气体Pr =0.71左右，但对于液体系统，则可能很高，
如Pr 数=50，Pr=ν / α, 如用一般的双分布函数模型，粘性系数及热扩散系数都相关到两松弛时间参数，必然一个大，一个小，大带来精度问题，而小带来计算量问题及稳定性问题
敬请大家给些建议，

ywang

我比较懒，没看完，，，汗，，，

就题目来回答吧，高Ra数的话，有些方法可以搞，比如用插值、fdlbm、fvlbm，还有一些忘记了，你去查一查，有一段时间，出来了不少paper搞高Re数和高Ra数的。

原帖由 yzgdy 于 2010-3-8 14:24 发表

                               
登录/注册后可看大图

自然对流模拟中，Ra=gβ△TH^3/(αν)，如Ra增大，温差参数设置不变，H不变，Pr数不变的情况，
（1）如热膨胀系数也不变，则一般调小粘性系数或热扩散系数，但这样松弛时间变小，为得到一定稳定状态或者某一无量纲时 ...

ywang

你这儿调Pr用的双分布函数法，，，我们算变物性也是这个法子

此外，调Pr还可以用MRT，这个思路我们也曾打算做，08年左右郭老师他们有篇文章里解决了这个问题，你查Du rui + Guo

原帖由 yzgdy 于 2010-3-8 14:41 发表

                               
登录/注册后可看大图

另一个相关问题：高Pr 数的问题，对于气体Pr =0.71左右，但对于液体系统，则可能很高，
如Pr 数=50，Pr=ν / α, 如用一般的双分布函数模型，粘性系数及热扩散系数都相关到两松弛时间参数，必然一个大，一个小，大带 ...

yzgdy

谢谢！指点迷津!