Effect of Rayleigh Accelerations Applied to an Initially Moving Fluid

  • PDF / 495,074 Bytes
  • 13 Pages / 414.72 x 648 pts Page_size
  • 45 Downloads / 203 Views

DOWNLOAD

REPORT


203

EFFECT OF RAYLEIGH ACCELERATIONS APPLIED TO AN INITIALLY MOVING FLUID ROBERT F. DRESSLER NASA Headquarters, Washington, D. C. S. J. ROBERTSON AND L. W. SPRADLEY Lockheed Missiles & Space Company, Huntsville, Alabama ABSTRACT Although the instability of an initially motionless fluid subjected to a Rayleigh-type acceleration is well known, this type of excitation applied to an initially moving fluid has not been studied previously. We have therefore performed numerical analyses for two-dimensional unsteady thermal convection in circular cylinders under variable low-g conditions associated with space flight. When an acceleration vector is applied parallel to the thermal gradient for a fluid at rest, no convection results for the stable direction, and an instability leads to Rayleigh convection for the opposite direction. However, when the acceleration has a component which is orthogonal to the gradient, convection always results at any Rayleigh number; this is the usual situation during space flight. We therefore study the resultant effect on convection when both types of acceleration are applied, concurrently or sequentially, and when the resultant vector varies in direction with time. Our results indicate that for space flight conditions, the Rayleigh accelerations impose significant, but not dominating, alterations in the established convection even when the Rayleigh number is less than

critical. NOMENCLATURE Symbol

Description

D g

cylinder diameter = 2 R acceleration due to gravity (sea level)

r

radial distance

R

cylinder radius

Ra

Rayleigh number = egATD

RaH

horizontal Rayleigh number

RaV

vertical Rayleigh number

T

temperature

To0initial AT

3

3 ega0D p

mid-point temperature temperature difference across circular cylinder

204

y

R

Cgold Side

Hot Side

I

I

I

\Ig

I

I)

/

H E

T + AT

T AT

Fig. 1.

Geometry for circular cylinder enclosure (Parallel dotted lines are initial isotherms.)

y4

Fig. 2.

Geometry and computational grid for numerical simulation.

205 Ra

S1.0

= 5000

-•Cold

.o

S id e.



Hot S ide

S0.8 E Hg

0.6 S• S•

2000 1000 0 -IOOO

0.4 >



-5000

0o.2 o

Rav= 1000

0.0

2.0

1.0

0

3.0

Time, t (dimensionless) Fig. 3a.

Spatial maximum velocity for various Ra values of excitation superimposed on established'steady state flow for RaV = 1000.

0.5 -•

U) .0

••RaH

5000

0.4

(2000

5 x

1000 0

0.3 ---

-1000 -5000

S0.2

Cold Side

'5 0.1

Hot Side

RaV = 5000 En 0.0

Vg 0.25

0.5

1.0

1.25

Time, t (dimensionless) Fig. 3b.

Spatial maximum velocity history for various RaH values of excitation superimposed on previously established steady state flow for RaV = 5000.

206 Description

Symbol t t*

tdimensionless

time transient response time to reach 63% (1 steady state 2 time = vt/R

e- 1 ) of

dimensionless response time Vt /R v

velocity

vmax

spatial maximum velocity

v max

vmax for zero horizontal acceleration component

AVmax

Vmax - Vmax

v

dimensionless velocity -

0

128v D2 v

x, y

rectangular coordinates (Fig. 1)

(X

thermal diffusivity

P