『求助』关于欧拉方程处理亚音速发散的问题？

zhangyinbo · 发表于 2003-10-31 09:52:37

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？注册

x

我做了一个程序只能计算马赫数1.2以上的，小于1.2的不能计算，不知为什么，是不是在处理亚音速时候要关掉人工粘性或者加入焓阻尼措施，望各位大侠指教！！

boling2000 · 发表于 2003-11-2 03:17:49

你用的是什么计算格式？
说的详细一点。

zhangyinbo · 发表于 2003-11-2 21:39:22

我做的是二维非结构网格欧拉方程计算，四步runge-kutta时间推进，中心差分有限体积发的格心格式。加速措施：当地时间步长，焓阻尼结果还是发散

周华 · 发表于 2003-11-2 23:56:25

是不是边界条件的问题？超音速的边界条件设置和亚音速有区别。

boling2000 · 发表于 2003-11-3 08:58:39

如果远场足够远，远场边界条件条件设置成自由来流都应该可以计算的。

zhangyinbo · 发表于 2003-11-3 09:54:02

我计算的是内流，壁面采用有滑移条件，出口给定背压，判断ma是否大于1，若大于1，保持上一步迭代的值不变，小于1，利用背压改变能量方程。进口判断初场ma是否大于1，如果ma>1,进口给定初值不变，否则外推。

周华 · 发表于 2003-11-3 13:11:48

如果不保密的话，帖上来几个图，大家看了才能提建议，现在这样盲人摸象实在猜不透你的问题。

zhangyinbo · 发表于 2003-11-3 15:18:04

[这个贴子最后由zhangyinbo在 2003/11/03 03:24pm 第 1 次编辑]

网格图片，左边进口，右边出口，上下为固体壁面

zhangyinbo · 发表于 2003-11-3 15:18:55

ma＝1.8压力分布图

zhangyinbo · 发表于 2003-11-3 15:22:00

[这个贴子最后由zhangyinbo在 2003/11/03 03:33pm 第 2 次编辑]

马赫数分布图，进口条件总温，总压，马赫数，出口背压，进气角，这个计算程序是我自己动手编的，所以模型取的简单，就是一个斜流板，但是问题很多，程序只能计算ma》1的初场，其他的都不能计算，改变时间步长也不行，结果总是发散。马上就要毕业了，所以现在很着急。谢谢大家的帮助！！

		自动登录	找回密码
密码			注册

『求助』关于欧拉方程处理亚音速发散的问题？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。