" A$ w- y* |4 N4 V1 `/ ?: p
|
算例说明# [! v9 D4 j2 Y8 }* L! H
| 本算例求解一180度弯曲管道中的流场和压力场,管径0.5cm,弯曲管道的半径为2.5cm,流体的粘性系数0.04g/(cm.s) ,流体的密度1.0g/cm3 。( v) {1 H2 v& {. q. `5 @, h
计算的初始条件:给定弯曲管道中的流速为0,压力为0。# T* i- X$ ^# h* z4 k, o
计算的边值条件:入口流速u=60cm/s ,v=0;管壁固定u=v=0。
, Z8 q& p, q$ P( K2 U8 S% v出口压力为0。
" M; x9 ?! X; s* ^0 z! v& G- H: k本算例能很好的处理对流占优的流体力学问题。, o# F/ O% u! o0 m( V
|
|
8 a, h- E* g8 L问题的控制方程:
* M2 V8 w( P9 ^: D(1)动量方程: http://www.fegensoft.com/images/hydrod5.gif
1 ]7 ?+ V; f M6 `(2)连续方程: http://www.fegensoft.com/images/hydrod6.gif5 g0 B$ \- Z" v( [5 a( {' ~
其中u表示流速矢量, ρ 表示流体密度, µ 表示流体粘性系数, f 表示流体体力,p 表示流体压力。' K" n5 k# f# I8 G
|
| # q! I- C Z% o5 `/ m
本算例采用分步算法:6 _' g, U, o" G6 G0 k" N
( \# |) Y6 g a- r# W! w
(1)先计算不含对流项的纳维斯托可斯方程;
# _' \) [1 T& ?& khttp://www.fegensoft.com/images/hydrod12.gif
3 J& v3 \. m+ m1 C G
9 {1 u' Y' z& C9 K$ K4 S |
| (2)再计算对流项的方程。# x6 l! g# \5 @' ^6 e* @! c7 Z
http://www.fegensoft.com/images/hydrod13.gif
+ Y5 G! k B" _7 V. q1 r. H% M" V5 D5 v* M( h% E- f0 m3 \
|
结果图
! ^% U( B/ u( k ]2 J! N- A$ w/ J |