Loading in 2 Seconds...

Dr. R. Nagarajan Professor Dept of Chemical Engineering IIT Madras

Loading in 2 Seconds...

- 138 Views
- Uploaded on

Download Presentation
## PowerPoint Slideshow about 'Dr. R. Nagarajan Professor Dept of Chemical Engineering IIT Madras' - pelham

**An Image/Link below is provided (as is) to download presentation**

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript

Module 4 - Lecture 16

Momentum Transport: Flow in Porous Media &

Packed Beds

Dr. R. Nagarajan

Professor

Dept of Chemical Engineering

IIT Madras

FLOW IN PORUS MEDIA AND PACKED BEDS

- Volumes of interest may contain a solids fraction, f, made up of:
- Granular particles (sand, pebbles)
- Wool (steel wool, fiberglass, etc.)
- Gauzes, screens (woven metals)
- Porous pellets (absorbent, catalyst support)
- Void fraction, e = 1 – f
- Usually, packing geometry is random

FLOW IN PORUS MEDIA AND PACKED BEDS

- e> 0.8:
- Large void fraction
- Flow about each object may be considered as “perturbed” external flow
- Each wetted object contributes drag
- e < 0.5:
- View as internal flow through tortuous ducts between particles
- Used in following module

FLOW IN PORUS MEDIA AND PACKED BEDS

- In case of flow through a straight duct of noncircular cross-section, effective duct diameter
- For flow through a packed bed, effective average interstitial (duct) diameter

where

FLOW IN PORUS MEDIA AND PACKED BEDS

Effective diameter of each particle in bed

Therefore:

Appropriate Reynolds number for internal flow

where interstitial mass velocity

FLOW IN PORUS MEDIA AND PACKED BEDS

- fbed dimensionless momentum transport coefficient

Function of Re

- For a single straight duct of short length, dimensionless momentum-transfer coefficient

where

FLOW IN PORUS MEDIA AND PACKED BEDS

Hence:

and

Correlates well with experimental data (next slide)

Ergun’s approximation:

FLOW IN PORUS MEDIA AND PACKED BEDS

Experimentally determined dependence of fixed-bed friction factor fbed on

the bed Reynolds’ number(adapted from Ergun (1952))

FLOW IN PORUS MEDIA AND PACKED BEDS

- Laminar flow region:
- Rebed < 10
- fbed ≈ 150/Rebed
- Darcy’s law:
- Linear relationship between G0and (-dP/dz)
- Effective local permeability = G0n/ (-dP/dz)
- >> intrinsic permeability of each particle
- Fully-turbulent asymptote: (Burke-Plummer)
- Rebed > 1000
- fbed ≈ 1.75

FLOW IN PORUS MEDIA AND PACKED BEDS

- Above equations are basis for most practical pressure-drop calculations in quasi-1D packed ducts
- Can be used to estimate incipient fluidization velocity (by equating –Dpto bed weight per unit area)
- Can be generalized to handle multidimensional flows through isotropic fixed beds
- Can be used to estimate inter-phase forces between a dense cloud of droplets and host (carrier) fluid

For a combustion turbine materials-test program it is desired to expose specimens to a sonic but atmospheric pressure jet of combustion-heated air with a post- combustion (stagnation) chamber temperature of 2000 K.

- What should the pressure in the combustion chamber (upstream of the nozzle) be?
- What should the shape of the nozzle be?
- How much air flow (g/s) must be handled if the exit jet must be 2.5 cm. in diameter?

d. What will the exit (jet) velocity be?

e. What will be the flow rate of axial momentum at the nozzle exit?

f. By what factor will the gas density change in going from the nozzle inlet to the nozzle outlet?

- Upstream (Chamber) Pressure

Isentropic gas flow from T0 = 2000 K to sonic speed @ 1 atm. Assume

Therefore,

and

Consider a dilute acetylene-air mixture at atmospheric pressure and 300 K with a C2H2 mass fraction of 4.5%. If the heat of reaction of acetylene is 11.52 kcal/gm C2H2, estimate:

- The speed with which a Chapman- Jouguetdectonation would propagate through this mixture (km/s);
- The stagnation pressure (atm) immediately behind such a detonation wave;
- The corresponding Chapman –Jouguet deflagration speed (km/s)

SOLUTION 2

Therefore

And, since Ma 2 =1(Chapman-Jouguet Condition)

conclusion: Detonation Speed=1.98 km/s ((Ma)1 =5.7). Stagnation Pressure:

PROBLEM 3

Consider the steady axisymmetric flow of hot air in a straight circular tube of radius aw and cross sectional area A

Conditions (at exit):

p=1 atm (uniform)

T= 1500 K (uniform)

aw= 5 mm

PROBLEM 3

Suppose it has been observed that the axial-velocity profile is, in this case, well described by the simple equation:

where U=103 cm/s. Using this observation, the conditions above, answer the following questions:

a. If the molecular mean-free-path in air is approximately given by the equation:

PROBLEM 3

Estimate the prevailing mean free path l and the ratio of l to the duct diameter- i.e., relevant Knudsen number for the gas flow:

What conclusions can you now draw concerning the validity of the continuum approach in this case?

b. Calculate the convective mass flow rate (expressed in g/s) through the entire exit section. For this purpose assume the approximate validity of the “perfect “ gas law, viz,:

PROBLEM 3

Here p is the pressure (expressed in atm), M is the molecular weight (g/g-mole)(28.97 for air), R=82.06 (univ gas const), and T is the absolute temperature (expressed in kelvins). Also note that for this axisymmetric flow a convenient area element is the annular ring sketched below

( where is the unit vector in the z- direction).

PROBLEM 3

c. Also calculate the average gas velocity at the exit section and the corresponding Reynolds’ number;

PROBLEM 3

d. Calculate the convective axial momentum flow rate (expressed in g. cm/s2) through the exit section. Is your result equivalent to Why or why not?

e. Calculate the convective kinetic-energy flow rate (expressed in g. cm2/s3). Is your result equivalent to Why or why not?

f. If, in a addition to the axial component of the velocity vz, the air in the duct also has a swirl component how would this influence your previous estimates ( of mass flow rate, momentum flow rate, kinetic energy flow rate)? Briefly discuss.

g. If the local shear stress is given by the following degenerative form of Newton’s law:

at what radial location does maximize? Calculate the maximum value of and express your result in dyne cm-2 and Newton m-2. Calculate the skin-friction coefficient, cf (dimensionless), at the duct exit. At what radius does take on its minimum value? Can- be regarded as the radial diffusion flux of axial momentum? Why or why not? Does the rate at which work is done by

the stress maximize at either of the two locations found above? Why or why not?

h. Characterize this flow in terms of flow descriptors and defend your choices.

SOLUTION 3

therefore

h. Descriptors:

Continuum Laminar

“Incompressible” Quasi-one-dimensional

Newtonian (viscous) Internal

Steady Single-Phase

- Estimate the drag force (Newtons) per meter of length for each of the following long objects of transverse dimension 5 cm if placed in a heated air stream with the following properties.
- a. A circular cylinder.
- b. A thin “plate” perpendicular to the stream (i.e., at 90o incidence).
- c. A thin plate aligned with the stream (i.e., at 0o incidence).

- In each case qualitatively discuss how the drag is apportioned between “form” (pressure difference) drag and “ friction” drag
- For Part ( c), is the application of laminar boundary-layer theory likely to be valid? (Briefly discuss your reasoning) If so, what would be the estimated BL thickness, , at the trailing edge of the plate, i.e., at x=L? Suppose two such adjacent plates were separated by a distance much greater than --- would they strongly “interact” with respect to momentum transfer?

d. Justify the use of an incompressible Newtonian fluid CD(Re, shape) curve to solve Part (a) ( involving the gas air) by showing that is small enough under these conditions to neglect .

SOLUTION 4

Momentum Transfer to (Drag on) Immersed Objects

Drag/meter of axial length=? for objects of transverse dimension 5 cm. in U =10 m/s, air @ 1 atm., 1200 K.

a. Cylinder in Cross flow

SOLUTION 4

Most of this drag is due to the p(q) distribution- that is, “ form” drag.

b. Plate Normal to flow: check literature

c. Plate Aligned with flow:

SOLUTION 4

In this case

and

Since ReL<106 (approx.) we expect flow in the momentum defect Boundary Layer to be laminar. Then

where

SOLUTION 4

This drag is entirely due to - i.e., it is “friction drag”

- Reconsider the sonic jet test facility specified in problem1 from the viewpoint of turbulent jet momentum exchange (mixing) with the surrounding atmosphere, and the entrainment of that atmosphere
- a. Calculate the Reynolds’ number at the nozzle exit and compare it to the value above which such jets almost certainly lead to turbulent mixing with the surrounding atmosphere.

b. Estimate the appropriate value of by assuming that the relevant density is about the arithmetic mean between . How much larger is the effective turbulent momentum diffusivity, , than the intrinsic momentum diffusivity of the jet fluid ?

c. Estimate the downstream distance at which the time-averaged velocity (axial momentum per unit mass) along the jet centerline will be reduced to 10% of the initial jet velocity (axial momentum per unit mass) as a result of momentum diffusion. Compare this to the result that would have been obtained had the jet

remained laminar (with kinematic viscosity ).

d. At this location what is the approximate ratio between the entrained (laboratory air) mass flow and the “primary” (combustion-heated air) jet?

SOLUTION 5

Momentum Transfer: Turbulent Round Jet

SOLUTION 5

We find: z=88.6 cm =0.89 m

(This would have been 0.46 km if there had been no turbulent enhancement in momentum diffusion.)

where

SOLUTION 5

Exercises:

1. Calculate the time-averaged profile at this location.

2. Can you estimate the centerline, time averaged temperature and CO2(g) concentration at this point? (Itemize and discuss the underlying assumptions.)

Download Presentation

Connecting to Server..