|
|

楼主 |
发表于 2012-6-7 16:51:00
|
显示全部楼层
原帖由 hilllike 于 2012-3-22 21:03 发表
李老师,请教柱坐标的奇异点怎么消除的?
我用有限体积写了个程序,压力在奇异点附近的控制体积是个三角形,解方程却发现不对.
在奇性轴周围,利用物理量的旋转特性,给出奇点另一侧的值,从而在奇点周围也可以采用高精度差分法计算。
例如在柱坐标系下,r=0是奇点。 在r=0处不布置网格点(如rj=(j-1/2)*h , j=1,2....N). 对应坐标线(seta=const), 利用柱坐标的性质,给出奇点另一侧的值,从而可以在j=1 点也可以采用高精度差分法计算。
当然,这种方法稳定性不是太好,奇性轴周围的流场容易变形。
最好的方法是在奇点采用有限体积法计算。 |
评分
-
1
查看全部评分
-
|