[american institute of aeronautics and astronautics 26th thermophysics conference -...

11
v Y AIAA 91-1311 A Compatibility-Satisfying Method for Laminar Channel Flows with a Heated Fence S. M. Liang National Cheng Kung University Tainan, Taiwan, R. 0. C. W. B. Tsai Far East Engineering College Tainan, Taiwan, R. 0. C. AIAA 26th Thermophysics Conference June 24-26, 1991 / Honolulu, Hawaii For permission to copy or republish, contact the American Instie of Aeronaufics and AstrOnOutiCS 370 L'Enfant Romenade, S.W., Washington, D.C. 20024

Upload: wen-bin

Post on 15-Dec-2016

213 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

v

Y

AIAA 91-1311 A Compatibility-Satisfying Method for Laminar Channel Flows with a Heated Fence S. M. Liang National Cheng Kung University Tainan, Taiwan, R. 0. C. W. B. Tsai Far East Engineering College Tainan, Taiwan, R. 0. C.

AIAA 26th Thermophysics Conference June 24-26, 1991 / Honolulu, Hawaii

For permission to copy or republish, contact the American Inst ie of Aeronaufics and AstrOnOutiCS 370 L'Enfant Romenade, S.W., Washington, D.C. 20024

Page 2: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

AIM-91-1311-CP A COMPATIBILITY-SATISFYING METHOD FOR LAMINAR CHANNEL

FLOWS WITH A HEATED FENCE

W

Shen-Min Liang*

National Cheng Kung IJniversity, Tainan, Taiwan, R. 0. C

and

Wen-Bin ’hi**

Far East Engineering Collegc, Tninnn, Taiwan, R.O.C.

A b s t r a c t

The laminar forccd convection associated with its hcat

transfer of a two-dimensional channel with a heated fence

is Ynvestigated numerically. The Reynolds number, IZen,

based on thc channel height, I € , ranges from 100 to 800

and thc blockage ratio, h j H , from 0.25 to 0.4. TIic re-

sults were obtained with air as the working fluid. The

incompressible Navicr-Stokes equations governing the flow

are solvcd on non-staggered grids by using the primitive-

variable approach. The Neumann boundary condithn de-

rived from the normal momentum equation i s prescribed ror

solving thc pressure Poisson equation in which the source

term is corrccted at each iteration The multiple soltions of

pressure resulted from the Neumann condition i s uniquely

determined by fixing the pressure to be unity at the first

point on the upper wall. In order to verify the accuracy

of the prcscnt method, thc computed results arc compared

with the cxperiniental data of Tropea and Gackstatter for

those Ikyno lds numbers at which the flow maintains its

two-dimensionality in experiment. The comparison is sat-

isfactory. The heat transfer ratc of the channel flow is es-

tablished by correlating the overall average Nusselt number,

Nut, thc Reynolds number, Ren, and the blockagr ratio,

h j N , of the form: K t 7 0 . 6 7 4 ( R e ~ ) 0 ~ 3 ’ G ( ( A / H ) ~ ~ ~ ~ ~ ~ z .

‘I Associate professor, Instit,utc of Aeronautics and

‘U

-

.. ~~~ ~ _ ~ _ _

Astronautics. Mcrnbcr A I A A . * * Lecturer.

Copyright 01991 by authors. Published by the American

Institute of Aeronautics and Astronautics, Inc. with per-

mission. -i

Nommclat~llrQ

=::fence height and channel height, rrspcctively

= channel length

L reattachment length

7 local and average Nusselt numbers, respec-

tively

:: pressure of fluid

Prandtl number

= Reynolds number based on channel height

and fence height, respectively

= time and artificial time, respectively

= dimensional temperature

= dimensional temperature on heated elerncnt

:= dimensional temperature of inflow

= maximum velocity at upstreani

= dimensionless velocity components in Cnrte-

sian coordinates

:= Cartesian coordinates

: fence width

=: thermal diffusivity of fluid

= viscosity of fluid

: density of fluid

-~ dimensionless temperature, (T T m ) / ( T w

To3 1

In t roduc t ion

Separated flow near obstacles or steps has brrn a re-

search subject for several decadcs. This problem has rele-

Page 3: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

vancc in a multitude of engineering applications, arid simple

geometry also provides a good meams for testing thoorct-

ical or numerical methods [l-91. In recent ycars, the dc-

sign of high-performance heat exchangers and the cooling

of electronic eqnipments require accurate quantitative in-

formation on hcat transfer rates in separated flow regions.

Numerical simulation can economically meet the need it correct mathematical model could be established.

The present work concerns with a fully-developed, lam-

inar airflow past through a heated fence in a two-dimensional

channel, as shown in Fig. 1. The heat source is maintaincd

a t a constant temperature which is higher than thr up-

stream flow temperature. The fence is considered to be

thin, hence its width is fixed at one-tenth ol the channel

height. The geometry and flow parameters arc the fence

blockage ratio, h / H , and the Reynolds number, Re", based

on the the channel height. The objective of this work is to

develop an accurate numerical method for establishing a

heat transfer correlation of the flow. r . I h c flow over a rectangular obstaclc can, in gcneral,

induce two or three separated flow regions. Uergeles and

Athanmsiadis 11) have observed that a recirculation region

developed ahead of and behind the fence and a third rc-

circulation region with consatnt length is on Lhr top sur-

face of the obstacle, if the obstacle width-to-height ratio,

w l h , is greater than 1. They also showed that the length

of the downstream recirculation region behind the obsta-

cle increases almost linearly with the obstaclc width for

w / h < 4. Later, Tropea and Gackstatter [ZI further in-

vestigated this flow in a fully developed chanrlel Row alld

measured the locations and sizes of recirculation regions for

Reynolds number ranging from 150 to 4500 (bascd on the

channel half-width). They found that an abrupt decreac

in the length of the downstream recirculation region in the

case of high blockago ratio (2 0.5). The abrupt dccreasp

was attributed to an appearance of a recirculation region

on the upper wall. Their results for those Reynolds num-

bers at which the flow maintains its two-dimcnsionality in

experiment are nscd to verify the present numerical results.

Thc hcat transfcr rate in a laminar flow past a surface-

mounted rib was studied numerically by Hsieh and Ilnang

using thc wcll-known SIMPI,E/IID codc /IO]. I n tlrcir stzldy

the temperature on the lower wall and the rib mountcd on

it was niaitained at a constant temperature. A hcat thilris-

fer correlation was established as the form: Nui O . I I ~ : ~ I - W

for Reh 5 2 5 4 and w l h e: I ,

relates t,he overall average Nusselt number, tlici Itrynolds

number and the obstacle aspect ratio. However, thcir corn-

puted data are largely deviated from the corrclution f<jr

large obstacle aspect ratio, ( w l h 2 2 ) . Thus, an a c r ~ d r

method for calculating the incompressible now assoc i ;~ i , ,~ l

with its hcat transfer is desirable.

Recently, Gresho and Sani [ l l ] have analyacd thr wcII-

posed problem of the Navier-Stokes cquations associat,od

with appropriate initial and boundary conditions. mrcy

showed that the Dirichlct and Neumann boundary condi-

tions for the pressure Poisson equation derivcd rrorn t izF

momentum equations give thc same solution. Abclallah

arid Dreyer 112) confirmed the r e d l of Gresho aud Smi

on the inviscid stagnation point flow problem. J l o w e ~ ~ ,

the prescription of a Dirichlet boundary condition hr t , h

pressure Poisson equation seems to be difficult for a ~ l ~ ~ r ~ , ~ , ~ l

Row with ai obstacle. On the other hand, the prrscripi,ion

of the Neumann boundary condition can praducc r n u l t i p l v

solutions for pressure.

W

In this study, the Neumann boundary condition dc-

rived from the normal momentum equation is prescribed.

Th? multiple solutions of pressure resultcd from thc Ntrtr-

mann condition are uniquely determined by assuming thc

upstream pressure to be unity at the first point on the up

per wall. As Abdallah I131 pointed ou t that t h r compati-

bility condition which relates the source term i n tlic prcs-

sure Poisson equation and the Neumann condition rimy not

he automatically satisfied on a non-staggered grid without

using consistent schemes. In order lo satisfy the compati-

bility condition, a uniform correction for tlic source term,

introduced by Briley [ I I ] , is implemented during iter;rt.iou

process.

LIy solving the energy equation, tbc authors have fuoltrlil

that from past numerical experiments, upwind cliffcrcn<-

ing for convective terms was more stable, bnt prrrd~lrnrl

much numerical damping for the thermal plume b c h i u r l .J

Page 4: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

the heated obstacle, compared to central differencing. On

the other hand, ccntral differencing could produce a large

dispersive effect in truncation errors, and hcncc gcncrates

oscillating solutions, leading to negative values of tempera-

ture. Since truncation errors are proportional to mesh sizes,

a uniformly distributed grid with small mesh sizes was uscd

to reduce to truncation errors and hcncc the associated dis-

persive effect.

'W'

Based on the above arguments, central differencing is

used for approximating all the spatial derivatiws, convcc-

tivc terms included, in thc momentum, energy and pressure

Poisson equations. The pressure Poisson equation, modi-

fied by addition of a false-transient term, is integrated in an

artificial time. These discretized equations were solved hy

the two-step Alternating Direction lniplicit (ADI) mathud

1151 as done by Ghia et al. 1161. Steady-state solutions are

obtained through an itcrativc procedure.

M a t h e m a t i c a l Formula t ion

A. Governing Equations

Assuming there are no viscous dissipation and no cx-

ternal forces present, the equations governing the two- di-

mensional, incompressible, laminar channel flow with a

heated fence are the full Navier-Stokes equations:

v

4

Continuit.y Equation D - : V . V = O (1)

Momcntum Equation Vt + V . V V = -Vp+V2? /Rei, (2)

Energy Equation Bt t -?. V? = V 2 6 / ( R e 1 , P 7 ) (3)

-t 4 -->

Thr ahovc equations have been properly nondirnensional-

ked. The spaces variables x , t ~ have been normalizcd by I I ,

the velocity components u , v by l J , the pressure p hy p l J 2 .

Thc energy equation (3) is decoupled from other equations

and can be solved once the velocity distribution is obtained.

Instcad of directly solving Eqs. (1) and (2), an auxil-

iary equation was generally derived to replace the cont,inu-

ity equation by taking the divergence of Eq. (2):

v 2 p = s (4)

with

a S =: - [ - ( u u z + vuy) -1 ( u , ~ ~ -i- u, , ) /Re1~1

+ - [ - (uvz vvy) + ( v m ~t rea] Ilt

az a v aY

The pnrpvsc of retaining the x , y derivatives of u , I, and the

local dilatation D in the source term S is to damp out, the

niimcrical insatbilities during iteration process 1171, becallsc

D may not he zero at the incipient stage. If D is zero, the

continuity equation is satisfied, and Eq. (4) reduces t,o

V2P = 2(u,v, - .,%)

Thus, Rqs.

appropriate boundary conditions.

B. Boundary Conditions

(2)-(4) are the equations to be solvcd with

Since the upper and lower wall surfaces of the chan-

nel are assumed to be adiabat.ic, afJ/aylwa~t = 0. The di-

mensionless temperature on the heated fence is mity. The

boundary condition on the body surface is the no-slip con-

dition, u = v = 0. At the upstream boundary, a fully devel-

oped velocity profile is imposed and the inflow temperature

is zero. At the downstream boundary, zero gradients of

both the velocity and the ternperaturc are prescribed. The

boundary condition used for the pressure Poisson equation

is the normal momentum equation, except the first point

of the upper wall at which the pressure is imposed to he

unity.

The boundary conditions prescribed are summarized

below: (refer to Fig. 1)

1) u = v = 0,

P,Ps and P,P5;

2) u = 4y( l - y),

3) u = v = 0,

8 , = 0, p , = u , , / R e H , on sides P2P3,

v = 0, 0 = I , p = 1, on side P IP2 ;

8 = 1, p , = t ~ , , / R e r ~ , on sidcs PSPB and

p7ps;

4) VI - 1) = 0, 0 = 1,

5) u , = vl: .:- 0 , = 0,

p , - ayy /Re l i , on side P,P,;

p , := ( - u u , + u,,,,)/RerI, on side

p3p.1.

Numer ica l Method

A. Finite-Difference Approximation

Central differences are used to approximate all the spa-

tial derivatives, convective terms included, i n thr momen-

tum, pressure Poisson and energy equations. Thc pressure

Poisson equation is modified by addition of a false-trasient

term

As suggested by Harlow and Welch, 1171 the local dilatation

Page 5: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

D"+' in S is forced to be zero at each iteration. Hence,

the approximation of a D / & reduces to ~ - D " / A t . The D"

ran eliminate numerical instabilities during iteration pro-

cess. The two-step alternating direction implicit method

of Peaceman and Rachford 1151 is used to solve the dis-

cretizewd momentum, pressure Poisson and energy equa-

tions. One-sided difference formaulas of first-order accu-

racy are used for first derivatives. Fictitious points outside

the walls were introduced when central differences were

used to approximate second derivatives. The pressure at

the boundary except at one point is linearly extrapolated

from the known quantities at interior points through the

normal momentum equation

El. Compatibility Condition

Since the pressure Poisson equation is solved for the

pressure, the divergence theorem implies that the compat.

ihility condition:

( 5 )

which relates the source term S and thc Nrumann rondi-

tion. Abdalla has showed that the compatibility conditionn

( 5 ) may not be automatically satisfied on a non-staggered

grid without using a consistent scheme; consequently, a

non-convergence process can occur. Ilence, a correction,

A S , introduced by Briley, is computed at each itcration.

A S = (// Sdxdy - / $ds)/ // dzdy

A C A

which is used to uniformly correct the source term S . When

the iteration process converges, the correction A S

approaches zero and the compatibility condilion is satis-

fied automatically.

C. Solution Procedure

The steps involved in the solution procedure are listed

below.

1) Make initial guesses for u, w , p .

2) Solve Eq. (2) for u, u at n + 1 time step.

3) Evaluate aplan on the boundary from the normal mo-

mentum equation and compute A S .

4) Correct S using A S to obtain new pressure from Eq.

(4).

5 ) Calculate the pressure at the boundary from known

quantities at interior points.

6) Check for convergence criterion

for all i , j , where r is a preassigned value. The proce-

dure is repeated from step 2 to step 6 until c o n d i t i o n

(6 ) is satisfied.

7) Solve Eq. (3) for the temperature from the convergcd

solution (u, u) using the same convergence criterion ( 6 ) .

Results and Discussion

Figure 2 shows the computed streamlines at / I F I l -

400 and 800 for h / H = 0.25 and at R ~ H = 200 and 400

for h / H = 0.4. It is seen that the sizes of rccirculation

regions ahead of and behind the fence increase with t h

Reynolds number. The recirculating intcrrsity of the recir-

culation flow behind the fence also grows with the Rcynolds

number, because more energy is entrained into the recircu-

lating flow from the main stream for highcr Rcynolds num-

ber. By comparing Figs. 2b and 2d, it is indicated t.liat

at the same Reynolds numbcr, the higher fence produces

larger recirculation regions because of stronger bluckagr ef-

fect. The length of the downstream recirculation rcgion,

called the reattachment length, L,, at, various Itcynolds

numbers are plotted in Fig. 3. Notc that thc reatt,acli-

ment Icngths have been normalized by the fcnce Iicight,

h. The comparison of the computed realtachincnt lengths

for h / H L 0.25 with the experimental data of Tropea and

Gackstatter is quite satisfactory. No comparison was rnadr

for the case of h / H = 0.4, since no experimcntal data were

avaiable. The almost linear variation of the reattachrncnt

length with Reynolds number for lower blockagc ratio is

well predicted.

Figurc 4 shows pressure distributions at dilli.rcnt

Reynolds numbers for h / H L 0.25. Note that thc value

of the pressure contour not indicated lies between bwo a&

jacent indicated-value contours. It is clearly seen that thew

are a high pressure region ahead of the fence and a low pres-

sure region developed behind the fencc, resultcd fro111 the

blockage effect. Since the flow is fastly accelcrated nrar

W

Page 6: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

the upstream corner, the flow acceleration leads to a large

pressure gradient in the vicinity of the corner. In the case

of Re" = 200, constant pressures at, downstream stations

wcrc wcll predicted, as shown in Fig. 4a, which incans that

a constant pressurc gradient is formed and the flow is fully

redeveloped at downstream. It is expected that similar re-

sults could be obtained for other Reynolds numbers if thc

computational domain is taken large enough. Moreover, a

lower pressure region behind the fence is developed, and

extends to the upper wall with a positive pressure gradicnt

in the 3 o'clock direction.

W

Figurc 5 shows the surface pressure distribution along

the upper wall. It is seen that a locally low surface pressure

region occurs on the upper wall behind the fencc. Because

of the formation of the low pressure region, the pressure is

decreascd with increasing Rcynolds number and blockagc

ratio. The upper-wall boundary layer at the low pressure

region may results in a flow separation which has bcen ob-

served in experiments 121. The tendancy to flow separa-

tion is also indicatcd in the skin friction coefficient distri-

butions. Figure 6 shows thc skin friction coefficient distri-

bution along the upper wall. The skin friction coefficient

first obtains a maximum and sharply decreases to a mini-

mum which is closer t o zero for larger Reynolds numbers. If

the skin friction coefficient changes its sign from a positive

value to a negative value, a recirculation flow occurs. Fig-

ure 7 shows the lower-surface pressure distributions. The

lower-surface pressure distributions are qualitatively simi-

lar to those (Runs 15 and 17) obtained by Leone and Gresho

1181, using finite element methods.

-

Figure 8 shows the temperature contours at various

Rcynolds numbers. The thermal plume is developed bchind

the heated fence, which is clearly seen in Figs. 8b-8d, in

particular, for larger Reynolds numbers. The zero-isotherm

stretches farther downstream as the Reynolds number in-

creases. The isotherm patterns are very similar for larger

Reynolds numbers.

Define the local heat transfer rate by the Nusselt num-

ber: as an

Nu = --

where n denotes the direction normal to Ibe fcnce surface,

and its average by the average Nusselt number:

- 1 Nu = - l'? Nuds

S

where s denotes the length of the fence surface. Figure 9

shows the distributions of the local Nusselt number along

the front, top and rear faces of the fencc. On the front

face, the local Nusselt number is found to increase with the

Reynolds number, and obtains its maximum at the top cor-

ner, The formation of the maximum local Nusselt number

is resulted from the flow acceleration near the corner where

convection dominates the heat transfer. On thc top face,

the local Nussclt number decreases along the surface. On

the rear face, the local Nusselt number changes very little,

sincc the heat transfer is primarily contributed by ronduc-

tion, not by convection. A similar result of local Nusselt

number distributions was obtained for h / H = 0.4. How-

ever, for h / W = 0.4, lhe peak value of the local Nusselt

number is larger than that for h / H == 0.25 at the samc

Reynolds number. The reason is that the flow accelera-

tion near the corner for h/H I 0.4 is fastcr than that for

h / H = 0.25. The variations of the average Nusselt number

on the front, top, rear faces and the overall average Nus-

selt number with the Reynolds number for h / H = 0.25 are

shown in Fig. 10. In general, the average Nussclt number

on each face increases with thc Reynolds number, except for

the rear face case. In this case, the average Nusselt number

obtains its maximum at Re 600. A result similar 1.0 that

for h l H := 0.25 was obtained for h / N == 0.4. Figure 1 1

shows the variations of the overall average Nusselt number,

Nut, versus the Reynolds number for two different values

of h / H . It is seen that the overall average Nusselt num-

ber for h / H I 0.25 is higher than that for h / H .I.. 0.4 at

a fixed Reynolds number. Namely, the lower the blockage

ratio, the higher the overall average Nusselt numbcr. By

taking the fcnce blockage ratio into account, a correlation

is obtained:

-

%, :: c ( R e , ) a ( h / H ) b for w = O.1H

where a = 0.376,b = - 0 . 2 9 2 , ~ = 0.674. Note that these

parameters may he functions of the fence width. The cor-

relation indicates that the overall average Nusselt uumhcr

Page 7: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

is proportional to the Reynolds number, but inversely to

the blockage ratio.

Conclusion

A compatibility-satisfying method of the primitive-vari

able approach has been developed for solving two- dimen-

siona.1 incompressible channel flows with an isolated, heated

obstacle. The pressure Poisson equation, used to replace

the continuity equation, is solved with the prescription of

the Neumann boundary condition which is derivcd from

the normal momentum equation. The multiple solutions of

pressure resulted from the Neumann condition is uniquely

determined by fixing the pressure to be unity at the first

point of the upper wall. To ensure the satisfication of the

compatibility condition which relates thc source term in

the pressue Poisson equation and the Ncumann condition,

a correction process of Briley is implemented to uniformly

correct the source term on a non-staggcred grid. The ac-

curacy of the present method is verified hy comparing the

computed reattachment lengths with thc cxpcriniental data

of Tropes and Gackstatter. The heat transfcr rate of the

channel flow for w ~ 0.lH is established by the correlation:

%, :~ 0.674(Re~)0~37G(h/H)~~0~292.

Acknowledgement

This work was supported by National Sciencc Coun-

cil of R.O.C. under the contract of NSC 79-4001-3-006-48.

The authors are grateful to Dr. C. Gau for his helpful com-

ments.

References

1. G. Rcrgeles and N. Athanassiadis, “The Flow Past a

Surface-Mounted Obstaclc,” Trans. of the ASME, J .

Fluids Eng., Vol. 105, pp. 461-463, 1983.

2. C. D. Tropca and R. Gackstattcr, “Thr Flow over

Two-Dimensional Surface-Mounted Obstacles at Low

Reynolds Nnmbcrs,” Trans. of the /ISM?:. J . Fluids

Eng., Vol. 107, pp. 489-494, 1985.

3. B. F. Armaly, F. Durst, J. C. F. Pereira arid B.

Schonung, “Experimental and Theoretical Invcstiga-

tion of Backward-Facing Step Flow,’’ J . Fluid Mech.,

Vol. 127, pp. 473-496, 1983.

4. F. Durst and C. Tropea, “Flows over Two-Dirncll.iionaI

Backward-Facing Steps,” Structure of Comp1e:c Tirbu-

lent Shear Flow, IUTAM Symposium Muraeille, 11. I)n-

mas and L. Fulachier (Eds). Berlin, Springer V<!rlag,

1982.

i /

5 . W. Aung, “An Experimental Study of Laminar Ileal

Transfer Downstream of Backsteps,” J . Heat l i n n r l e r ,

Vol. 105, pp. 823-829, 1983.

6. W. Aung, and A. Baron, F. K. Tsou, ‘‘ Wall Indcpcm

dency and Effect of Initial Shear-Layer Thickness i n

Separated Flow and Heat Transfer,” Int. J . I leal and

Mass Transfer, Vol. 28, No. 9, pp. 1757-1771, 1985.

7. C. Berner, F. Durst, and D. M. McEligat, ”Flow Around

Barnes, J . Heat Transfer,” Vol. 106, pp. 743-7119, 1984.

8. I,. P. Hackman, G. D. Raithby, and A. 1%. St,rong,

“Numerical Prediction of Flows Over Backward-Facing

Steps,” Int. J . Numer. Method i n Fluids, 1’01. 4, tip.

711-724, 1984.

9. S. S. Hsieh and D. Y. Hnang, ”Flow Charactcristirs of

Laminar Separation on Surface-Mounted Ribs,“ A / A A

J. , Vol. 2 5 , No. 6, pp. 819-823, 1987.

10. S. S. FIsieh and D. Y . Huang, “Numerical Computnlion

of Laminar Separated Forced Convection on Surfim-

Mounted Ribs,” Nurner. Heat T~ansjer-, Vol. 12, p p

335-348, 1987.

11. P. M. Gresho and R. I,. Sani, “On Pressure I iuunt lary

Conditions for the Incompressible Navier-Stokos ISqua-

t,ions,” Int. . I . Numer. Methods in Fluids, Vol. 7. p p

1111-1145, 1987.

12. S. Abdallah and .J . Dreyer, “Dirichlet and Ncumann

Boundary Conditions for the Pressure Poisson 1Sqiia-

tion of Incompressible Flow,” Int. J . Numer. Methods

in Fluids, Vol. 8 , pp. 1029-1036, 1988.

v

13. S. Abdallah, “Numerical Solutions for the Inrompross-

ible Navier-Stokes Equations in Primitivc Variables Us-

ing a Non-Staggered Grid, 11,” J . Cornput. I’hys., Voi.

70, pp. 193-202, 1987.

14. W . R. Briley, “Numerical Method for Predicting T l i r m -

Dimensional Steady Viscous Flow in Ducts,” .I. Corn-

put. Phys., Vol. 14 , pp. 8-28, 1974. i/

15. D. W. Peaceman and 11. FI. Rachford, “Thc N w n c r i c a l

Page 8: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

A

I

I I 1 '7 % 77-,-r, , , , , , , , , , , , , , , ,- i-T-A ?z?

Page 9: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

20.0

15.0

10.0

5.0

0.0

_t_ h/H = 0.40 . h,,, = 025 } Present results ___ h/H = 025. Experimental data [2]

--?--.--, 7 F--' 200 400 600 800

Reynolds Number, ReH I

Fig. 3 Comparison of computed reattachment lengths with experimental data.

Fig. 4 Pressure contours for diKerent Reynolds numbers, h / H = 0.25: (a) Re, = 200,

(b) Re, = 400, (c) Rei, = 600, (d) ReH = 800.

Page 10: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

'd

w

1.2 -

~ Rei, -200 .- ......... Rei, = 400

Reir = 600 ...... Rei, = 800

, , , , , , 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 0.0 9.0 10.0 11.0 12.0

S, Distance along upper surface

Fig. 5 Pressure distributions along upper surface, h / H = 0.25

0

_ _ ~ _ . Reit = Z O O .... Reit = 400

Rei< = 600 ...... Reti = 000

...............

0.2 , I l I I / I

0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 12.0 Distance along lower surface

Fig. 7 Pressure distributions along the bottom wall, h l H = 0.25.

1.0

Page 11: [American Institute of Aeronautics and Astronautics 26th Thermophysics Conference - Honolulu,HI,U.S.A. (24 June 1991 - 26 June 1991)] 26th Thermophysics Conference - A compatibility-satisfying

/ , , I > ,

7 0.1 0.2 0 3 0.4 0 5 0.6

S, distance along fence surface

(b) R e , = 400, ( c ) Re , = 600, (d) Re,, = 800. r i ~ 9 Variation of local Nusselt number along

fence surface, h j H = 0.25

- --- Overall N u -- Front. Surface t_ Upper Surface A Rear Surface

P

/'.' ___I/

0 no0 500 700 c Reynolds Nirmher. Rei1

20.0

17.5

~ 15.0 a,

5 12.5

a

z a, 10.0

z 7.5

Y

m u,

- 5

a, w g 5.0

5 2.5

0.0

Fig. 10 Variations of average Nusselt numbers on front, top and rear faces and ovcrall average Nusselt number with

Reynolds number, h l H = 0.25.

Fig. 11 A heat transfer correlation for w / l I : 0 .1