6 q9 v* a: ^/ u+ s* ]- V6 I
|
算例说明2 x3 R3 w2 o6 U9 u6 ~) }; @8 g3 O9 s
| 本算例求解一180度弯曲管道中的流场和压力场,管径0.5cm,弯曲管道的半径为2.5cm,流体的粘性系数0.04g/(cm.s) ,流体的密度1.0g/cm3 。- |9 V4 s; @; `; E8 b1 _0 g4 ~6 \
计算的初始条件:给定弯曲管道中的流速为0,压力为0。) r6 E. R. d- E9 a& p- l. e
计算的边值条件:入口流速u=60cm/s ,v=0;管壁固定u=v=0。. g2 p5 Z4 c7 z4 W7 L
出口压力为0。- D5 f9 N3 h9 ^; }" w: L7 [* j2 s7 a
本算例能很好的处理对流占优的流体力学问题。# Z* I; g# a1 Z) R& w8 S9 N
|
|
! V$ w2 x9 I& F7 V* D9 Q( z问题的控制方程:2 k( v5 i8 U) e/ k$ i
(1)动量方程: http://www.fegensoft.com/images/hydrod5.gif ) N3 K; B# t% ^
(2)连续方程: http://www.fegensoft.com/images/hydrod6.gif) H0 _5 F/ e; g$ \% w
其中u表示流速矢量, ρ 表示流体密度, µ 表示流体粘性系数, f 表示流体体力,p 表示流体压力。
9 l9 I2 B3 E5 ~9 A |
| 2 z% c, E4 E3 z; o
本算例采用分步算法:3 _, H; @$ ?- }: D8 L" g$ r
5 J3 r3 k6 g6 t! t: X6 A) r5 g
(1)先计算不含对流项的纳维斯托可斯方程;; M) P5 q8 l6 o& }' k
http://www.fegensoft.com/images/hydrod12.gif
5 }# B1 e' @7 `, i! {" ?! X& l8 \! j! [. Q V2 t/ p6 @
|
| (2)再计算对流项的方程。2 ]* O4 D [2 U/ H- l
http://www.fegensoft.com/images/hydrod13.gif
4 x' j; P; @4 y- p7 S$ Y: Q9 Q7 n: C/ j4 F2 @6 N
|
结果图
" k, r) d# A8 Q/ W |