1 O. B6 b: Y9 I, }/ a) ~" z9 l
|
算例说明9 D( P9 b5 b% S9 H$ O6 p
| 本算例求解一180度弯曲管道中的流场和压力场,管径0.5cm,弯曲管道的半径为2.5cm,流体的粘性系数0.04g/(cm.s) ,流体的密度1.0g/cm3 。. e- _, V- c1 p$ x/ P& f4 S
计算的初始条件:给定弯曲管道中的流速为0,压力为0。4 ~; m8 a& I: U. D3 y/ P) i- h
计算的边值条件:入口流速u=60cm/s ,v=0;管壁固定u=v=0。) W6 A9 T" z3 m- p
出口压力为0。
$ ^6 {0 H2 b* w& r: T. O+ @# U本算例能很好的处理对流占优的流体力学问题。
7 ]+ x' Z! N6 a6 Q0 t7 d |
|
1 F8 B% q0 C, B! e问题的控制方程:3 ^) |; B" [: `' r
(1)动量方程: http://www.fegensoft.com/images/hydrod5.gif
# Z f* e9 z) t3 N/ k4 Y) M(2)连续方程: http://www.fegensoft.com/images/hydrod6.gif3 c) y- U% t+ m
其中u表示流速矢量, ρ 表示流体密度, µ 表示流体粘性系数, f 表示流体体力,p 表示流体压力。$ M, w' q) L" w
|
| ' f( m9 l& V+ k a" c
本算例采用分步算法:. ?6 V# A6 X& ~- l' l
' L3 ]2 X: O# \. L/ @0 p& C
(1)先计算不含对流项的纳维斯托可斯方程;+ o7 }! A( J5 r7 ?5 t& _
http://www.fegensoft.com/images/hydrod12.gif ' K c; f% Y3 o' m8 @' Q
0 S% i( e; A, n V
|
| (2)再计算对流项的方程。
: R9 {: F2 `) T9 r3 z2 S- Vhttp://www.fegensoft.com/images/hydrod13.gif
, G5 P$ M6 O9 r5 W4 k0 P0 R' Z
) N5 N0 f! G7 R1 c3 P |
结果图
0 b' _# t' \1 E' j) U V |