Bingham体的连续流和阵性流研究中无量纲参数讨论(刘大有,段新平,章书成,余斌)

热度57票 浏览73次 【共0条评论】【我要评论 时间:2001年1月01日 19:17

Bingham体的连续流和阵性流研究无量纲参数讨论水利论文 u`,peA/f

刘大有1段新平1章书成2余斌2水利论文;p8h(KX:@"JT6RPLT
(1.中国科学院 力学研究所,2. 中国科学院 成都山地灾害与环境研究所

"L G!H jbM0

摘要:本文首先对轴对称管流、二维管流和明渠连续流研究中无量纲参数的选择进行评述,通过对这些有分析解的流动的讨论,说明采用无量纲表达式的优点、无量纲方法的多样性,以及适当选择基本量的重要性。最后讨论了Bingham体阵性流研究中基本量的合理选择。

7jv&k'jS0

关键词:无量纲参数;Bingham体;连续流;阵性流

H c(P,S1q Txz0

基金项目国家自然科学基金(49771003)和中国科学院“特别支持领域 山地灾害 泥石流与滑坡”的资助。水利论文1T;mY.iEH!x3y.P\

作者简介:刘大有(1939-),男,中科院力学所研究员。水利论文3L;p5w D7rn eH

1 引言水利论文*LTH A+rI#`%jcL+__

  高含沙水流、泥石流中的浆体、工业中的水煤浆和其它许多化工原料的流变性质,都接近于Bingham体,它们都有显著的屈服应力τ0和较大的粘度μ0(又称刚度系数η)。

I/sZ"Oz"vo'nQ i vo0

  近年来,人们针对一些典型流型(如二维明渠流动、二维管流和轴对称管流等)对Bingham体的阻力已有很多研究,还引入了各种无量纲参数进行讨论,例如,Reynolds数Re、Hedstrom数He和阻力系数f等。根据实际的重要性,人们主要研究层流运动,但也有一些论文讨论了Bingham体的紊流阻力。水利论文#O\T1X|$i

  本文在一些典型连续流阻力研究的基础上,对无量纲参数的选择进行了讨论,并在此基础上研究了Bingham体阵性流的无量纲表示和独立无量纲量的个数。水利论文dp8Tn8F$Lo

基于以下几点考虑,本文不讨论Bingham体的紊流运动。水利论文4tu:D%v M

  1. Bingham体与牛顿体的差异仅在于是否存在屈服应力。如果流场中平均应力水平较高,远大于屈服应力,那么屈服应力的存在与否对流动的影响就比较小,就可以用牛顿体代替。事实上,许多介质也并非理想的Bingham体,当平均应力水平较高时,Bingham体近似未必比牛顿体近似更优越。当流动进入紊流区时,应力水平往往已很高,失去了当作Bingham体研究的必要性;水利论文1V3P9s2|ne'~7@

  2. 当流动进入紊流区时,总有一定强度的紊流脉动,而Bingham体运动的重要特征是有流核存在。当有紊流脉动时是否还有流核?流核中是否有紊流脉动?流核如何影响紊流脉动?都是远没有研究清楚的问题;

+M+X_ ir){$o0

  3. 高粘度介质在较高速度下运动时,有很大的机械能损失,它需要很大的功率输入,同时产生大量的热,这些热会显著改变介质的流变性质。水利论文}8X ZrR*N0y0wt3Ji

2 几种定常连续流的流速与阻力水利论文MSZ0Bu3MTh

2.1 定常、轴对称管流

b`/y'Z,t,a6X`0

  下面先用量纲理论对定常(/t=0)、轴对称(/φ=0,分速度w=0)、充分发展的管流(/x=0)进行分析。根据定性分析可知:流量和平均速度U依赖于负压梯度(-p/x,以下简写为P)、管径D、介质密度ρ、屈服应力τ0和粘性系数μ0对长度为Δx的一段管子所作的力平衡分析可得到壁面摩阻τw的表达式水利论文[w9u:C&]

(1)
πDτw=πD2P/4, τw=DP/4(2)

如选D、ρ和τ0作为基本量,它们的值分别作为长度、密度和应力的单位(基本量度单位),则时间、速度、加速度、粘性系数和压强梯度的单位分别为(Dρ1/2τ-1/20)、(ρ-1/2τ1/20)、(D-1ρ-1τ0)、(Dρ1/2τ1/20)和(D-1τ0),无量纲的流量[Q]、平均速度[U]、负压梯度[P]、壁面摩阻[τw]和粘性系数[μ0]分别定义为这时,式(1)可表示为[1]

"])b6CH"a%L H0

水利论文Ag)M MB[7~:ew

(3)

8B6{Fax$Z(?js0
(4)

量纲分析可以将依赖于五个参数的函数ξ简化为依赖于两个无量纲量的函数ζ,但是,它给不出该函数的具体形式。量纲分析方法是通用的,它能有效地减少未知函数中独立参数的个数,而未知函数具体形式的寻找,则只能具体问题具体分析了。可以利用该问题的相关理论进行分析求解,也可采用经验的方法。例如,对于牛顿流体,如果人们已知,流量与粘性系数成反比、与负压梯度成正比,就可显著简化该函数:ζ([P],[μ0])=C[P]/[μ0],其中C是一待定常数。对于Bingham体,虽然流量与负压梯度不再成正比,但仍与粘性系数成反比,所以函数ζ仍可得到适度地简化:ζ=ψ(P)/[μ0],其中ψ是一未知的单变量函数,它仍要用实验方法或理论方法来确定。但是它比确定双变量函数ζ简单得多。下面采用流体力学方法求解函数ψ(或者说ζ)。水利论文a Ax _a+~VG

  对于定常、轴对称管流,轴向(x方向)动量方程可简化为(τrx为粘性应力)

,{#hq2B:TS0

r-1(rτrx)/r+P=0水利论文;cI[|#W8P-]b

(5)

x o/g1]&CJ@zh0

解方程得到水利论文J8dk [c0udD+^

τ00u/r≥-τrx(r)=rP/2水利论文#_:e2A W:j

(6)水利论文7I4z7m'wyV0f

u(r)=

8l$sx1\P/}`)}&K0
{P(D-2r)(D-4rcr+2r)/(16μ0)0.5D≥r>rcr

(7)

eU?s-?3d,{0K0
P(D-2rcr)2/(16μ0)

rcr≥r≥0

%D-}x iM m0

rcr水利论文)`n2`T0}1s

{0.5Dτ0/τw=2τ0/P P≥4τ0/D

(8)

pOm+}}W,p0
0.5DP<4τ0/D

:w:A aVK3M0

(9)

k_i0pz-qn0

这就是著名的Buckingham方程[2],其中rcr是流核半径,它们的无量纲形式分别为水利论文`:Y(O2H;m!f

.iA8l~n I*{Z3k+_#~0
{τ0/τw=[τw]-1=4/[P][P]³4

(10)水利论文0f j4p,IQoU7M

1[P]<4

o_rvz0

(11)

X[%U]t^3|d3H+R9S(G0
以上采用流体力学方法给出了函数ψ([P])和ζ([P],[μ0])的具体形式,即式(11)的右边。图1(a)给出了轴对称管流的([μ0[U])~[P]曲线。

Sz xY%bR$FhV'~'P0  无量纲方法不是唯一的,例如,人们也可选μ0、ρ和τ0作为基本量,它们的值分别作为粘性系数、密度和应力的单位,则时间、长度、速度、加速度和压强梯度的单位分别为(μ0τ-10)、(μ0ρ-1/2τ-1/20)、(ρ-1/2τ1/20)、(μ-10ρ-1/2τ3/20)和(μ-10ρ1/2τ3/20),无量纲的流Q*量、平均速度U*、压强梯度P*、壁面摩阻τ*w、管径D*、流核半径r*cr和Reynolds数Re分别定义为水利论文'z-j@ d o*C

图1 轴对称管流中无量纲平均流速 压降曲线
"UfZ(`Hc"h0Fig.1 Relationship between the averaged dimensionless velocity and pressure drop in an axisymmetric pipe flow水利论文6y ["L6jM
水利论文;Qs2F0}$bmMY7D$M

水利论文4p,_#P;D(Fl&Fk

(12)

-W x1E.Q ZhjF0

水利论文z] |"D9^

水利论文@\ q*x[%z'HhE0k

(13)

HCEE$S \'t m0

Eo8Y(Y/?"o!rU0

(14)水利论文+ld)lv;p%@)WU(F M

l Y,JU([e0

(15)

3_Q,Y _5sv!H%@ J0

这时Buckingham方程(9)的无量纲形式为(参见图1(b))水利论文3y%x*v yt

水利论文 G4F [b'c ~"f

    D*P*4

(16)水利论文2a&k[bi'c9B5_;g+c,C J

    D*P*<4

许多文献中引入如下定义的Hedstzom数He和阻力系数f

QgU"I b|0

He=D2τ0r02,f=2τw/rU2水利论文"[N8oNZn

(17)

p/| v`4C7T0

并将Buckingham方程(9)写成如下的无量纲形式[2,3]

UIf/X d9? j0

mt%TO*A~ F:eS0

(18)水利论文z)M?)EO.x%P}

采用这个公式,由负压梯度P(=4τw/D)求平均流速U时很不方便,因为在f和Re中都有未知量U。为便于求解,有些作者提出忽略式(18)中的高次项,其实完全没必要[24]。本来公式(9)并不复杂,写成(18)那样的无量纲公式是把简单的问题搞复杂了。

9UXEq[T0

  引入式(11)和(16)那样的无量纲公式是有意义的,它把公式(1)中涉及很多参量的U与P之间的函数ξ,简化为只依赖单参量的[U]与[P](或者U*与P*)之间的函数ψ。这样简化以后就很容易制成图或表,以备查阅。水利论文"E j$u x(I

  在引入无量纲量时应注意,不宜选未知量作为基本量,因为它会使许多无量纲量都成为未知量,不利于求解。水利论文4~Q3Zy1o9J t

2.2 定常、二维管流

$i v8A+w6DDeZ/a0

对于管高为2h的定常、二维、充分发展的管流,用类似的方法可得到如下各式

2x`N)P"@;h7iYd0

Q=2hU=2xa(P,h,r00)

)RXh,n;rZ0
(1)a

2τw=2hP,τw=hP水利论文9Ut`5N4g {AjPD5Z

(2)a

'C,u6bx0n G7T(X0

τzx/z+P=0

F!KfiPWC p0
(5)a

τ00uu/|z|≥1-τzx(z)|=|z|P水利论文m1KSum| s0ip`1S

(6)a

#D!}/Cpt ^Y7VL1ylp0

u(z)=水利论文SI:D4AMl1p(RZX

{水利论文?0h*xw}

P(h-|z|)(h-2zcr+|z|)(2μ0)h≥|z|>zcr

(7)a

*Ze{cst&?0
P(h-zcr)2/(2μ0)zcr≥|z|≥0

zcr水利论文 f;o XlwW&cK,?~ ?Z

{0/τw=τ0/PP≥τ0/h

(8)a水利论文8}f L8F{4[

hP<τ0/h

二维的Buckingham方程为[2]水利论文 Vl"i/t2K%Q

水利论文U;Z8f"K0[,L2\{`

(9)a水利论文0K1W&m@_5S lJ9X0O*h

以2h、ρ和τ0为基本量,各无量纲量的定义和无量纲的Buckingham方程分别为水利论文I_4S il6R$f Gf

,V1K:m)t Hy&CQ0

(3)a水利论文K q)R;M$Bmh"J

水利论文!QV&}\-g0@~*@7G7c+i

(4)a水利论文$C ~.A?uN}!B

;pos)V3x@0E0
[P]≥2

(10)a水利论文Q*W4g@;Ya

[P]<2

N+q\'I-A0
[P]≥2

(11)a

:l"xj0AYw5U0
[P]<2

如以μ0、ρ和τ0作为基本量,各无量纲量的定义和无量纲的Buckingham方程分别为

3oWC!J`h4a+v(N0

y!Z,oL+x'r P0

(12)a

-cEUhvl)_:x0

.a'~ p?FH0

(13)a水利论文1y|~;y1Wrp6n


a],lB6ffoRj0图2 二维管流(或明渠流)中无量纲流速 压强(或水深)曲线
#_W vZ@0Fig.2 Relationship between the dimensionless velocity and pressure drop (or depth) in a 2-D pipe flow(or channel flow)
水利论文VQ)c L^;H

]7ke Dp_2b,P#z)o0
h*P*≥1

(14)a

U?j&Zd.hl azg0
h*P*<1

水利论文7r.NxP_6q'Q

(15)a

w(R(j,`H0

水利论文T P)}G RU6R6g

h*P*≥1

(16)a

bt9J8r`z$B1Vm0g0
h*P*<1

用式(11)a和(16)a计算的结果分别示于图2。

Vh|3[BY0

  如引入下面(17)a定义的He和f,则无量纲的Buckingham方程如式(18)a所示水利论文Q9F N6toE:V

Xt@!l;Hi ^^%U0

(17)a

[#su`j-gq-V#w*WT0

L%aY&[\)Y?}V0

(18)a水利论文,]Y1w[EQ

2.3 定常二维充分发展的明渠流水利论文B4H"Yi4?KJ

  定常、二维、充分发展的明渠流的流速分布与定常、二维管流的下半部流动完全一样,但是,这时的驱动力P是重力的流向分力,P=ρgsin。依照通常习惯,在二维明渠流研究中,z=0平面位于床面,因此,只要在式(1)a、(6)a、(7)a和(9)a中用(h-z)代替|z|,用τzx代替|-τzx|,用2Q代替Q,就能得到定常二维充分发展的明渠流相应的公式,式(2)a、(5)a和(8)a完全不用改变。水利论文1Nez#az V8xA)fU _V

水利论文h7xQ8X`'B~

(1)b

0@{vbgz0

p R'L0Ul+d0

(6)b

O E)j{:l!D8R0

水利论文Q;y(K{ mW\8l'r

(7)b

p4h| Vg#g0{jR0

水利论文9xo6M;X0{U8k'`

(9)a水利论文w4XT!EZ0g

ri~3~k/Bd)n0G p0若以h、ρ和τ0为基本量引入无量纲量[U],[P],[μ0],[rcr]等,则只要用2[P]代替[P]、用0.5[μ0]代替[μ0]、用[zcr]代替[2zcr],二维定常管流的无量纲式(10)a和(11)a就转换为明渠流各式了,包括图2(a)

Q`(iRG5\0

水利论文;{kr F"Q/Pz+\b

(10)a水利论文!w0Uo[5W3y)b%\6J-Z

O?6u3GU&kVC&k0| @0
[P]≥1

(11)a水利论文h;C:et roa9R

[P]<1

8|efB4Pg0但在明渠流问题中,水深h常是未知的,用它作基本量得到的结果不便于使用。如以μ0、ρ和τ0作为基本量引入无量纲量U*,P*,h*,r*cr等,则各无量纲量的定义和无量纲的Buckingham方程(13)a~(16)a仍基本上适用,包括图2(b),只是Reynolds数的定义和Q*~U*关系需改为

UD'o ])I1T0

水利论文+iei J*^

(15)b

  在明渠流问题中,采用P/ρ(=gsin)、ρ和τ0作基本量可能更方便,它们分别是加速度、密度和应力的单位,则时间、长度、速度和粘度的单位分别为[ρ-1/2τ1/20(gsin)-1]、[ρ-1τ0(gsin)-1]、[ρ-1/2τ1/20]和[ρ-1/2τ3/20(gsin)-1],无量纲的平均流速〈U〉、水深〈h〉、粘度〈μ0〉分别定义为水利论文R4A't)P\+ax0H{C

水利论文y NaMz,f+D

(19)水利论文*W8Gq'q:k

_](p&@!e C1T%U0
{<h>=hrsinq/τ0≥1

(20)水利论文{)vz&c} R:j!te

<h><h>=hrsinq/τ0<1
这时Buckingham方程可表示为(见图3)
    
(21)

  如引入如下定义的Hedstzom数He和阻力系数f,则相应的Buckingham方程如式(18)b所示

"b&V6{:o RNp0

 

YCP(j q:]d0
(17)a

 

/A:Dq(UMd r7C0
(18)a


Jz'S*W5yD5VRN0水利论文,l E8NRc
图3 明渠流中无量纲流速-水深曲线
H.Mv/H:^h7hR0Fig.3 Relationship between the dimensionless velocity and depth in a open channel flow水利论文!W5|KWe3E(u)k?2m

不论由流速求水深,还是由水深求流速,这个方程都不方便,因为水深h和流速U在He,Re和f三个变量中,都至少出现两处。水利论文 DH+aBy

3 二维、充分发展的阵性流水利论文V ?4C(D!\

  粘性泥石流运动多呈阵性,这可能与Bingham体连续流的运动不稳定性有关[5]。阵性泥石流在干床上运动时,由于它在运动后都会留下一定厚度的残留层(能长时间保持不动),因而运动的介质会越来越少,这就是所谓的“铺床过程”。对于一定的介质,在一定坡度(tanθ)的二维渠道内,残留层厚度不可能超过临界值zcr(=τ0/ρgsinθ)。我们的数值计算表明(另文发表),泥石流通过后留下的残留层一般也不小于临界值,基本上等于临界厚度zcr。所以,一定量的泥石流体在已铺过床的渠道内运动时(假设原残留层内的介质与运动的介质相同),流体总量基本上不变。这样,如果渠道的几何条件沿程不变,那么泥石流体的运动速度和形态最终将趋于稳定,这就是充分发展的阵性流。如果有一个跟随阵性流一起运动的坐标系,并在该坐标系内考察运动,那么充分发展的阵性流的运动是定常的。

2l3Km(u(VT H1D0

  由定性分析可知,充分发展的阵性流的运动速度U依赖于介质密度ρ、屈服应力τ0、粘性系数μ0、运动介质的(单宽)总量Q、重力加速度g和倾角,即

M/d#A1qd6HRWV/B;|0

U=F(ρ,τ00,Q,g,θ)

%``!sI'|7x0

(22)水利论文-C\/v1x!F{"_!O

如选ρ、τ0和g作为基本量,它们的值分别作为密度、应力和加速度的单位,则时间、长度、速度和粘性系数的单位分别为(ρ-1/2τ1/20g-1)、(ρ-1τ0g-1)、(ρ-1/2τ1/20)和(ρ-1/2τ3/20g-1),无量纲的阵性流速度〈U〉、运动介质总量〈Q〉、粘性系数〈μ0〉和残留层临界厚度〈zcr〉分别定义为

Zy.w9U"[9f \f0

-b/l SpE5Gu"{ Z0

(23)水利论文 b D2])a6f,B4|GJ

对于阵性流,可采用如下定义的Reynolds数水利论文U"XV7c,H w

水利论文gB^;v g;|}

(24)

C/|Rj| f{OM0

这时式(22)可简化为[1]

Z9B\mucy0

〈U〉=F(1,1,μ*0,Q*,1,θ)=Φ(μ*0,Q*)

.pT;F8O#q&Q0

(25)

C g&V+a)B h-o:Po0

函数Φ的具体形式需通过求解流体力学方程得到。

%s ?;qF~&d{Y IU0

  在实验室坐标系中,非定常、二维、不可压缩流体(无量纲化的)运动方程为(为书写方便,在本段落中略去了表示无量纲的符号〈〉)

|!o\hmbc0

_#p&}p(`c |0
(26)

水利论文] N'M N2z

(27)

c,C*{ F F CCz0
(28)

Bingham体的应力本构关系为[6]水利论文pYb } QX&Y2I(h V

水利论文o'Z(j/Z:Q!xw2h

(29)水利论文Q&epxu~,db

其中,应变率eij=0.5(ui/xj+uj/xi),q=(2eijeji)1/2。由于边界条件都是齐次的,不涉及任何特征量,在此从略;在初始条件中,对充分发展运动状态有影响的仅仅是运动介质的总量〈Q〉,故也不在此详细列出。

8C,[6Do4K pic ~a2E0

  与前面研究过的几种连续流不同,阵性流中某些区域(例如龙头区)有显著的垂向(z向)运动,因而垂向动量方程(28)不能略去。由于垂向动量方程中有因子cosθ,所以,即使在<Q>和<μ0>的定义式(23)2和(23)3中采用gsinθ代替g,方程组(26)~(29)中仍然包含无量纲量,不能减少独立无量纲参数的数目。也就是说,在无量纲函数Φ中,除在变量<μ0>和<Q>中含有因子sinθ外,θ还可能以其它方式影响函数Φ。所以,我们宁愿采用式(23)定义的<μ0>和<Q>,Φ一般说来是<μ0>、<Q>和三个无量纲变量的函数。虽然还可以采用其它的无量纲量,如Reynolds数Re,Hedstzom数He和阻力系数f等,但是,我们推荐采用<μ0>、<Q>和θ,因为它们分别代表了介质物性,阵性流总量和渠道的几何,都有简明的物理或几何意义。例如,在<μ0>的定义中,除不变的常量g以外,都是与介质物理性质有关的量,即密度ρ、屈服应力τ0和粘度μ0。此外,在这些无量纲量的定义中涉及的都是已知量,不包含像运动速度U和龙头高度等未知量,使用比较方便。

*k ohf*r@P0

4 结论

$G)@\qih;X0

  (1)不论是轴对称管流、二维管流、还是明渠流,都可用Buckingham方程(轴对称管流的(9),二维管流的(9)a和明渠流的(9)b)计算平均流速U。有多种方法可将Buckingham方程变成无量纲形式,例如轴对称流的式(11)、(16)和(18),或二维流动的式(11)a、(16)a和(18)a,或明渠流的式(11)b、(16)a、(21)和(18)b,使原来有六个变量的关系变成只涉及三个变量之间的关系,便于制成图表。但是,这些无量纲关系使用的方便程度有很大差别,例如,式(18)、(18)a和(18)b中有两个变量涉及未知量U,(18)b中有两个变量涉及水深h,因而不便于使用。

Y0S q*F3D,m6}0

  无量纲化时,一般不选未知量作基本量。在管流研究中,由于管子的特征长度(管径D或管高2h)通常已知,所以,可选D(或2h)、ρ和τ0作为基本量。水利论文0@z.cN!`U0E}-mn4]

  (2)二维明渠中的连续流与二维管流的下半部流动基本上一致。由于明渠流中的水深h常常是未知的,所以,无量纲化时,可选gsin、ρ和τ0作为基本量。水利论文L \]gu#c$y0l8~

(3)对于二维明渠的阵性流,可选g、ρ和τ0为基本量,这时,无量纲速度〈U〉依赖于<μ0>、<Q>θ和三个无量纲变量,其中,<μ0>代表了运动介质物性,<Q>代表了阵性流总量,而θ则反映渠道的几何特性,都有简明的物理或几何意义。

4X5l6ibD0

参考文献

2A}Ydss0

[1]谢多夫(Седов, Л。И。).力学中的相似方法与量纲理论。科学出版社,1982(沈青等译)。

0n%z [Dkp{c0

[2] 钱宁.高含沙水流运动.清华大学出版社,1989,第112-113页。水利论文/],\-r_ f)Z}^Tx)jO

[3] Chen,C.L.,Generalized viscoplastic modeling of debris flow, J.Hydraulic Engineering, ASCE. 114,pp.237-258.水利论文8YPn3sv R9Z lk"zY

[4] Shamber, D.R.& MacArthur, R.C., One-dimensional model for mudflows, Hydraulics and Hydrology in the Small Computer Age,ASCE.Vol 2, pp.1334-1339.

8fY"G?SA l]-B0

[5]王兆印,林秉南,张新玉。非牛顿体的不稳定流研究.力学学报,22(3),1991.水利论文@4yoh e)F\6E

[6] 陈文芳.非牛顿流体力学.科学出版社,1984.水利论文*`L c$_HF)g$~

TAG: 余斌 段新平 章书成 Bingham 无量纲
顶:4 踩:11
【已经有42人表态】
12票
极差
5票
很差
2票
较差
7票
稍差
7票
稍好
5票
较好
3票
很好
1票
极好
下一篇:长江中下游鹅头型汊道演变规律的分析(马有国,高幼华)
上一篇:中国河口最大浑浊带刍议(沈焕庭,贺松林,茅志昌,李九发)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放