Servicios Personalizados
Articulo
Indicadores
Links relacionados
- Citado por Google
- Similares en Google
Compartir
R&D Journal
versión On-line ISSN 2309-8988
versión impresa ISSN 0257-9669
R&D j. (Matieland, Online) vol.10 Stellenbosch, Cape Town 1994
On Taylor-Proudman columns and geostrophic flow in rotating porous media
Peter Vadasz
Professor, Department of Mechanical Engineering, University of Durban-Westville, Private Bag X54mf , Durban, 4000 Republic of South Africa
ABSTRACT
Taylor-Proudman columns are a well-known phenomenon in rotating flows in pure fluids (non-porous domains). A theoretical formulation of the problem of incompressible fluid flow in rotating porous media is presented. The criteria for the relative significance of different terms in the equations are identified leading to a formulation which is based on the traditional Darcy's law but extended to include the Coriolis and centrifugal terms resulting from rotation. Finally a proof is provided showing that Taylor-Proudman columns exist in porous media as well. This occurs in the limit of small values of the porous media Ekman number. The corresponding consequences are that a stream function exists in this otherwise three-dimensional flow and this stream function and the pressure are the same in the limit of high rotation rates. This type of geostrophic flow means that isobars represent stream-lines at the leading order for small values of Ekman number.
Nomenclature
Latin symbols
êz a unit vector in the vertical direction
êw a unit vector in the direction of the imposed angular velocity
êg a unit vector in the gravity direction
ên a unit vector normal to the boundary
Ek the porous media Ekman number defined by eq.(3)
Frg gravity related Froude number, equals
Frw centrifugally related Froude number equals
h the bottom topography of the container
k(xi y, z) the dimensionless permeability function
k0a reference value of permeability
lc a macroscopic characteristic length
Pi dimensionless pressure
Pr the dimensionless reduced pressure generalized to include the centrifugal and gravity terms, equals
p rescaled pressure, equals Ek pr
r a coordinate in the radial direction (in a cylindrical system of coordinates)
Re macroscopic Reynolds number, equals
Professor, Department of Mechanical Engineering, University of Durban-Westville, Private Bag X54001, Durban, 4000 Republic of South Africa
Re∆pore size Reynolds number
qc a characteristic filtration velocity
q the dimensionless filtration velocity, relative to the rotating solid matrix
uhorizontal component of filtration velocity in the x direction
V horizontal component of filtration velocity in the y direction
w vertical component of filtration velocity
X position vector, equals xêx + yêy + zêz
xa horizontal coordinate (in a Cartesian system of coordinates)
y a horizontal coordinate (in a Cartesian system of coordinates)
z the coordinate in the vertical direction (in both Cartesian or cylindrical systems of coordinates)
Greek symbols
Φporosity of the porous domain
wc the angular velocity of rotation
v0 the kinematic viscosity
θa coordinate in the angular direction (in a cylindrical system of coordinates)
ψ a stream function defined by eq.(11)
Subscripts
0 reference values
c characteristic values
* dimensional values
Introduction
Rotating flows and heat transfer in porous media have a wide spectrum of applications in engineering and geophysics. The food process industry, chemical process industry and centrifugal filtration processes are some of the traditional applications. More explicitly packed bed mechanically agitated vessels are used in the food processing and chemical engineering industries in batch processes. As the solid matrix rotates due to the mechanical agitation, a rotating frame of reference becomes necessary. The filtration velocity is thus measured relative to this rotating frame of reference which is connected to the solid matrix. Other, modern applications emerged recently as a result of using the porous media approach to non-traditional disciplines including some domains in which the solid matrix is subjected to rotation. Among these applications, the flow of liquid in human tissues like the brain or heart, the development of porous turbine blades and cooling of electronic equipment subject to rotation (e.g. a rotating radar) may serve as examples. Vadasz [1] presented a more detailed discussion of these applications. Nevertheless, no reported research could be found on isothermal flow in rotating porous media. Probably, the main reason behind the lack of interest for this type of flow is that the isothermal flow in homogeneous porous media following Darcy's law is irrotational. However, for a heterogeneous medium with spatial dependent permeability the flow is not irrotational anymore. An example of flow in a rotating heterogeneous porous medium at high values of Ekman number was presented by Vadasz.[1]
In this paper further results are presented showing theoretically that Taylor-Proudman columns, which are a common phenomenon for rotating flows in pure fluids (non-porous domains), exist in porous media as well in the limit of small values of Ekman number.
Problem formulation
Transport phenomena in porous media are represented by a mathematical model at a macroscopic level. This representation is achieved by averaging over a Representative Elementary Volume (REV) the Navier-Stokes and other transport equations which are valid at the microscopic, pore-size scale. Different approaches for averaging have been proposed by Bear & Bachmat,[2] Whitaker,[3] Barrere, Gipoloux & Whitaker [4] and Du Plessis & Masliyah.[5] Eventually a set of equations is obtained at the macroscopic level which are an extension of Darcy's law to include inertial and other effects. As a consequence of the averaging process new variables are defined, e.g. the filtration velocity q is the average (over the REV) of the real velocity, and the pressure in porous media is the average (over the REV) of the real pressure. New properties are introduced as well through the averaging process, like the porosity Φ which is the ratio of the pore volume over the total volume of the porous domain, and the permeability k0, which has the units of square length, representing in principle at the macroscopic level the effective cross-sectional area of the microscopic flow. A major significance of the averaging approach is that it allows one to obtain theoretically the criteria for neglecting terms in the equations. For example, the porous media Reynolds number Re∆defined as Re∆ = ReDa =, controls the validity of Darcy's law. When Re∆is kept small the inertial effects are insignificant. However, the relationship between the familiar Reynolds number in pure fluids Re (non-porous domains) and the porous media Reynolds number is given by the multiplying factor Da, which represents a Darcy number and is defined as Da = This is the square of the ratio between a pore-size length scale and the macroscopic length scale of the problem, lc. As such, Da is typically very small (10-10 - 10-5), hence extending the validity of Darcy's regime to include a wide range of Re number values. However, the averaging techniques were traditionally applied to the Navier-Stokes equations in a non-rotating frame of reference. When the porous medium rotates a rotating frame of reference becomes necessary in order to keep the averaged equations valid, since the filtration velocity is defined relative to the solid matrix and the later rotates as a solid body. As soon as a rotating frame of reference is introduced two additional inertial effects should be incorporated in the model, i.e. the centrifugal and the Coriolis accelerations. The criteria for their relative significance is not controlled by Re∆alone but by other dimensionless groups, like the Ekman number. If Ekman number is very high (Ek →∞) then the Coriolis effect becomes insignificant.
As a result, it was concluded (Vadasz [6]) that the following dimensionless equations govern the incompressible flow in rotating heterogeneous porous media under isothermal conditions.
(ii) Darcy's law extended to include the Coriolis term
where q is the dimensionless filtration velocity, pr is the dimensionless reduced pressure generalized to include the centrifugal and gravity terms, k(x,y, z) is the dimensionless permeability function, êw is a unit vector in the direction of the imposed angular velocity and Ek is the porous media Ekman number defined in the form
where Φis porosity, wc is the angular velocity of rotation, k0 is a reference value of permeability and v0 is the kinematic viscosity.
Equations (1) and (2) are presented in a dimensionless form where the values of are used to scale the filtration velocity and pressure, respectively, and k0is used to scale the permeability function k*.
The Taylor-Proudman theorem in porous media
Equation (2) can be presented in the following form
Multiplying eq.(4) by and rescaling the pressure in the form p = Ek pr yields
Given typical values of viscosity, porosity and permeability one can evaluate the range of variation of Ekman number in some engineering applications. There, the angular velocity may vary from 10 rpm to 10 000 rpm leading to
Ekman numbers in the range from Ek = 1 to Ek = 10-3. The later value is very small, pertaining to the conditions considered in this paper. Therefore, in the limit of Ek → 0, say Ek - 0, and assuming êw = êz equation (5) takes the simplified form
and the effect of permeability variations disappears. Taking the 'curl' of equation (6) leads to
Evaluating the 'curl' operator on the cross product of the left-hand side of equation (7) gives
Equation (8) is identical to the Taylor-Proudman form for pure fluids (non-porous domains); it thus represents the proof of the Taylor-Proudman theorem in porous media and can be presented in the following simplified form
Results and discussion
An example of a Taylor-Proudrnan column
The consequence of the result presented in the previous section can be demonstrated by considering a particular example. Figure 1 shows a closed cylindrical container filled with a fluid saturated porous rnedium. The topography of the bottom surface of the container is slightly changed by fixing securely a small solid object (see Greenspan [7] for the corresponding example in pure fluids). The container rotates with a, fixed angular velocity we. Any forced horizontal flow in the container is expected to adjust to its bottom topography. However since equation (9) upplies for each component of q it applies in particular to w, i.e . But the irnpermeability conditions at the top and bottom solid boundaries require q . ên = 0
represents the bottom topography. The combination of this boundary condition with the requirement that yields w = 0 anywhere in the container. Hence, a flow over the object as described qualitatively in Figure 2 becomes impossible as it introduces a vertical component of filtration velocity. Therefore the resulting flow may adjust around the object as presented qualitatively in Figure 3. However since this flow pattern is also independent of z, it extends over the whole height of the container resulting in a fluid column above the object which rotates a,s a solid body. This is a demonstration of a Taylor-Proudman column in porous media, &s presented qualitatively in Figure 4.
Geostrophic flow in rotating porous media
A further significant consequence of equation (9) is represented by a geostrophic type of flow. It is observed that despite imposing a permeability which is a function of z as well, the flow at high rotation rates, i.e. Ek → 0, is independent of z. In particular , and the continuity equation (1) presented in Cartesian coordinates becomes
A stream function, ψ, can therefore be introduced for the flow in the x - y plane
Through the definition of the stream function ψ, given by equation (11), the continuity equation (10) is identically satisfied. Substituting u and v with their stream function representation given by equation (11) into eq.(6) yields
The conclusion resulting from equations (12) and (13) is that the stream function and the pressure are the same in the limit of high rotation rates (Ek → 0). This type of geostrophic flow means that isobars represent streamlines at the leading order for Ek → 0.
Conclusions
A theoretical formulation and proof of existence of Taylor-Proudman columns in porous media in the limit of small values of the porous media Ekman number was presented. The corresponding consequence leading to a geostrophic type of flow in porous media was discussed. Experimental confirmation of the theoretical results is recommended despite the practical difficulty of reproduction of experiments from pure fluids to porous domains.
Acknowledgement
The author wishes to thank the Foundation for Research Development for funding this research through the Core Programme Rolling Grant.
References
[1] Vadasz P. Fluid flow through heterogeneous porous media in a rotating square channel. Transport in Porous Media, 1993, 12, 43-54. [ Links ]
[2] Bear J & Bachmat Y. Introduction to Modeling of Transport Phenomena in Porous Media, Kluwer Acad. Publ., Dordrecht 1991.
[3] Whitaker S. Flow in porous media I: A theoretical derivation of Darcy's law. Transport in Porous Media, 1986, 1, 3-25. [ Links ]
[4] Barreré J, Gipouloux O & Whitaker S. On the closure problem for Darcy's law. Transport in Porous Media, 1992, 7, 209-222. [ Links ]
[5] Du Plessis JP & Masliyah JH. Mathematical modeling of flow through consolidated isotropic porous media. Transport in Porous Media, 1988, 3, 145-161. [ Links ]
[6] Vadasz P. Fundamentals of flow and heat transfer in rotating porous media. Proceedings of the 10th International Heat Transfer Conference, Brighton, UK, August 1994.
[7] Greenspan HP. The theory of rotating fluids. Cambridge University Press, Cambridge, 1980.
First received April 1994
Final version August 1994