: S. [2 K- a4 x/ n$ m |
算例说明
/ ~+ @3 u- F( i8 r5 M" g& c- s | 本算例求解一180度弯曲管道中的流场和压力场，管径0.5cm，弯曲管道的半径为2.5cm，流体的粘性系数0.04g/(cm.s) ，流体的密度1.0g/cm3 。; Y5 M' H. Q5 Q, r. V
计算的初始条件：给定弯曲管道中的流速为0，压力为0。* g7 x, L# B3 z3 A; i! Q* v; Z+ \
计算的边值条件：入口流速u=60cm/s ，v=0；管壁固定u=v=0。) E' e! u( }1 a, i* d) y7 ?
出口压力为0。
$ @2 K: ^: p: x$ _4 w5 ^本算例能很好的处理对流占优的流体力学问题。
: _5 y' Q" _# t8 o) e) z5 N |
| . s/ N7 y; d; y8 k- I! N
问题的控制方程：
6 \/ P# X, \7 |8 s9 O8 g' L1 [（1）动量方程： http://www.fegensoft.com/images/hydrod5.gif
: J1 z% w5 L7 R$ Z+ [6 U h（2）连续方程： http://www.fegensoft.com/images/hydrod6.gif' ?; {" T4 i" a& I3 g3 a% z5 B
其中u表示流速矢量， ρ 表示流体密度， µ 表示流体粘性系数， f 表示流体体力，p 表示流体压力。
. z2 Q2 {& g8 R& N' C |
|
# v# D. }/ ?) N* ?4 z: L本算例采用分步算法：* \: V7 M# p5 E3 C$ d6 n
: f1 a, F6 J: m. k+ ^9 _% p（1）先计算不含对流项的纳维斯托可斯方程；, |# [3 P- f- U5 F9 h7 V& t
http://www.fegensoft.com/images/hydrod12.gif ( K% I# K* d1 t8 M
& c$ H% s5 s% Z1 c8 {% m6 q
|
| （2）再计算对流项的方程。
$ T! S" [- [ h0 dhttp://www.fegensoft.com/images/hydrod13.gif
* a j/ F* K5 j- {) j
0 B1 @" [& G' b |
结果图
/ S2 d, R2 j3 j3 }: Z4 L4 { |
[复制链接]
发表于 2006-11-30 11:42:41
