Next Article in Journal
Classical Chaos Described by a Density Matrix
Next Article in Special Issue
Instability of Vertical Throughflows in Bidisperse Porous Media
Previous Article in Journal
Instability of Traveling Pulses in Nonlinear Diffusion-Type Problems and Method to Obtain Bottom-Part Spectrum of Schrödinger Equation with Complicated Potential
Previous Article in Special Issue
Spatially Developing Modes: The Darcy–Bénard Problem Revisited
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vadasz Number Effects on Convection in a Vertical Rotating Porous Layer, Placed Far from Axis of Rotation, and Subjected to Internal Heat Generation and Centrifugal Jitter

1
Eskom Holdings SOC, Generation, Asset Management, Johannesburg 2000, South Africa
2
Department of Engineering, School of Mechanical Engineering, Howard College Campus, University of Kwa Zulu, Durban 4001, South Africa
Submission received: 30 June 2021 / Revised: 20 August 2021 / Accepted: 23 August 2021 / Published: 1 September 2021

Abstract

:
The flow and heat transfer in a rotating vertical porous layer, placed far from the axis of rotation, and subjected to internal heat generation and centrifugal jitter, is considered. The linear stability theory is used to determine the convection threshold, in terms of the critical Rayleigh number. Typical liquids used in engineering applications and heavy liquid metals are used to demonstrate conditions at which the Vadasz number is sufficiently small to warrant the retention of the time derivative in the momentum equation. When considering low amplitude and high frequency approximation, the results show that vibration has a stabilizing effect on the onset of convection. The impact of increasing the Vadasz number is to stabilize the convection, in addition to reducing the transition point from synchronous to subharmonic solutions. In summary, when the Vadasz number is large, centrifugal jitter has no impact on the convection stability criteria. In contrast, when the Vadasz number is small, centrifugal jitter impacts the convection stability criteria.

1. Introduction

The classical problem of gravity driven free convection in fluid saturated porous media has been widely researched for various configurations. The engineering applications include, inter alia, nuclear reactor technology, food and chemical processing industry, thermal storage systems, etc. Pioneering past studies for a constant basic temperature gradient have been investigated by past researchers. Previous pioneering study [1,2,3,4] for non-rotating porous media set the foundation for researchers in heat and mass transfer in porous media. Later, rotation was analyzed in detail, and these studies are presented in [5,6,7,8,9,10,11]. The Vadasz number was first proposed in [12], specifically for the body of study presented in [9]. In almost all of the papers dealing with flow and heat transfer in porous media, the Vadasz number is neglected on the basis that it is understood to be large for the typical fluids used in engineering. The analysis in [9] considers scenarios where the Vadasz number could assume small values, thereby warranting its retention together with the time derivative in the momentum equation. The results presented in [9] discover that only when the Vadasz number is small to moderate, is the oscillatory mode of convection possible. Later study [13,14,15,16] introduced the effects of g-jitter on the stability of convection in porous media. The study in [13,14], the shows that the transition from synchronous to subharmonic solutions is characterized by a spike in the curve of the synchronous mode, just prior to the transition. Furthermore [16] shows that (gravity) g-jitter effects only affect the convection stability criteria for small to moderate Vadasz numbers. For large Vadasz numbers, the g-jitter has no impact on the convection stability criteria, in the absence of internal heat generation.
The focus of the current study is to consider the effects of centrifugal jitter for a vertical porous layer with internal heat generation and placed far from the axis of rotation. The readers are referred to pioneering studies for flow and heat transfer with internal heat generation with constant gravity [17,18]. Early study on g-jitter on porous media with internal heat generation are presented in [19]. Research involving specific configurations and internal heat generation in porous media are presented in [20,21,22]. The readers are also referred to an important book in porous media and applications [23] for further review/reading, together with an excellent reference list for further reading.
In typical engineering applications, one needs to cater for the impact of centrifugal jitter in a design or perhaps design the system to use the effects of centrifugal jitter to control the heat transfer. In the scenario, where the influence of centrifugal jitter on the stability of convection needs to be dialed out, one may opt to choose a fluid that allows for large Vadasz numbers. For large Vadasz numbers [16], vibration (g-jitter) was shown to have no influence on the stability of convection. The current study investigates this and determine if the same result is applicable for centrifugal jitter in the presence of internal heat generation. There could be scenarios where the stabilizing effect of vibration (g-jitter or centrifugal jitter) may be required in an engineering application. The fluid properties, and perhaps even the medium, may need to be designed in such a way that the resulting Vadasz number is small enough, to allow the effects of vibration (g-jitter or centrifugal jitter) to influence the stability criteria for convection.
The present research recovers the basic conduction flow and temperature profiles and then use the linear stability theory to solve for the characteristic Rayleigh number, for a rotating vertical porous layer, placed far from the axis of rotation, and subjected to internal heat generation. The equation for the centrifugal-jitter amplitude is cast into the canonical form of the Mathieu equation to ensure that the critical Rayleigh number (in terms of the Vadasz number) may be determined (for the low amplitude, high frequency vibration case). The results for the rotating vertical porous layer, placed far from the axis of rotation and subjected to centrifugal jitter is then be compared to the case of a horizontal porous layer subjected to g-jitter, with internal heat generation [24]. In principle, the current study extends the study of [24] to a vertical rotating porous layer that is placed far away from the axis of rotation, thereby allowing for a constant external centrifugal body force (that may be likened to gravitational body force in [24]). The study developed in [11] discovered that gravity plays a passive role in the stability of convection for the configuration studied in the current study. Therefore, although the gravity force is be shown in the formulation, its role is be passive. The stability criteria developed in the current study is then compared to the previous study [24] and a few important are drawn regarding the impact of the Vadasz number on the stability criteria. Based on the review of literature performed, the author is unable to locate study similar to that presented in this paper.

2. Problem Formulation

In the current study, the porous layer is placed far from the axis of rotation, in order to recover the case when the centrifugal effect is constant due to the offset distance X 0 * . When X 0 * > > L * , where L * is the porous medium length, the variation of the centrifugal force within the porous layer is small compared to the centrifugal force due to X 0 * . From hereon, variables with the superscript “*” refer to dimensional quantities, whilst those without the “*” refer to dimensionless variables. For the special case of the porous layer placed far away from the axis of rotation, the constant centrifugal force due to offset distance X 0 * is analogous to the case of convection in a differentially heated porous layer subjected to g-jitter. In the following, a comparison is made between the current case and the case involving gravity only to illustrate this point. Figure 1 shows a vertical fluid saturated porous layer, placed far away from the axis of rotation, and subjected to vibration and internal heat generation. In addition to internal heat generation, the porous layer is also differentially heated, with a precise definition of the boundary conditions later in this section. Internal heat generation refers to heat generation (as a source) from within the porous layer per unit volume. A typical example of internal heat generation in an engineering application could be a nuclear reactor. The porous layer is sandwiched between two impermeable vertical boundaries (a distance L * apart), subjected to centrifugal effects Ω * 2 X 0 * , and vibration b * ω * 2 sin ( ω * t * ) . At x * = 0 , T * = T C * , and at x * = L * , T * = T H * , and the Boussinesq approximation is applied to model the effects of the density variations. Here, T H * and T C * are the hot wall and cold wall temperatures, Ω * is the angular velocity, ω * represents the vibration frequency, and t * is the time variable. This results in the following system of dimensional equations for continuity, momentum, and energy:
* V * = 0 ,
( ρ o * ϕ * t * + 1 ) V * = χ * μ * { * p * + ρ * [ Ω * 2 X 0 * + b * ω * 2 sin ( ω * t * ) ] e ^ x ρ * g * e ^ z } ,
T * t * + V * * T * = κ * ( * 2 T * + q * k m * ) .
In the system (1)–(3), V * is the fluid velocity, p * is the pressure, the density is denoted by ρ * , the temperature is denoted by T * , the permeability of the porous medium is denoted by χ * , μ * is the dynamic viscosity, k m * is the fluid thermal conductivity, and e ^ x and e ^ z are the unit vectors in the x- and z- directions.
The governing equations may be non-dimensionalized using the scaling variables κ * / L * , μ * κ * / k 0 * , and Δ T * = T H * T C * for the filtration velocity components ( u * , v * , w * ) , reduced pressure ( p * ) , and the temperature variations T * T C * . For this scaling, κ * is the thermal diffusivity, and k 0 * is a characteristic permeability associated with the porous medium. The height of the porous medium L * is used to scale all the spatial lengths according to the relations x = x * / L * , y = y * / H * , and z = z * / L * . Applying the scaling factors to Equations (1)–(3) yields the following governing equations for the specific case of constant permeability (i.e., χ = 1 ):
V = 0 ,  
( 1 V a t + 1 ) V = { p R a ω 0 [ 1 + δ sin ( ω t ) ] } T e ^ x + R a g T e ^ z ,
T t + V T = 2 T + 1 .
The key non-dimensional parameters that emanate from the rescaling of the Equations (1)–(3) are the gravitational Rayleigh number, R a g = [ ρ 0 * β * g * H * k 0 * / ( μ * κ * ) ] ( q * H * 2 / k m * ) ; the centrifugal Rayleigh number, R a ω 0 = [ ρ 0 * β * Ω * 2 X 0 * H * k 0 * / ( μ * κ * ) ] ( q * H * 2 / k m * ) ; and the vibration amplitude δ = η F r ω 2 , where η = b * / H * and the Froude number F r = κ * 2 / ( Ω * 2 X 0 * H * 3 ) . Here, H * is the height of the horizontal porous layer. In Equation (5), Va is the Vadasz number, defined as V a = ϕ * Pr / D a (where Pr is the Prandtl number and Da is the Darcy number, defined as D a = k 0 * / L * 2 ); and the symbols V, T, and p r represent the dimensionless filtration velocity vector, temperature, and reduced pressure, respectively. In a previous study [9], the author has put forward a motivation for specific cases when the Vadasz number is small and is retained. In the current study, typical fluids used in engineering applications and selected heavy liquid metals were considered to determine the Vadasz number, Va, as a function of the porosity, as shown in Figure 2 and Figure 3. Here, we use the Kozeny equation of the form k * = L * 2 ϕ * 3 / [ 172.8 ( 1 ϕ * ) 2 ] to estimate the permeability when calculating the Darcy number, and the Vadasz number.
Figure 2 shows that when the porosity is low, the Vadasz number is very large, and this leads to the neglect of the time derivative in Equation (5). In the instance when the Vadasz number is very large, the presence of centrifugal jitter has no effect on the stability of convection. Figure 2 also shows that when the porosity is large, the Vadasz number could attain values that allow for the retention of the time derivative in Equation (5). In the case of high porosity media, the Brinkman model may need to be adopted in the formulation of the momentum equation. The introduction of the Brinkman term in Equation (2) is the subject of current study, and the study presented in [25] for a rotating vertical porous layer, excluding gravity and Coriolis effects (for both the Darcy and Brinkman models), is consulted in that analysis.
In contrast, for heavy liquid metals, it is clearly demonstrated in Figure 3 that that Vadasz number is orders of magnitude smaller across the entire range of porosity values. Figure 3 shows that across the range of porosity values indicated, the Vadasz number is of an order of magnitude that allows for the retention of the time derivative in Equation (5).
Therefore, for systems using typical liquids with high porosity media, liquid metals allow for the retention of the time derivative in Equation (5), which, in turn, allows the centrifugal jitter to impact the convection stability criteria. The solutions for the basic temperature and flow field are given as T B = 1 / 2 ( 3 x x 2 ) and w B = 1 / 2 R a g ( x 2 3 x + 7 / 6 ) . In providing a solution to the system of equations, it is convenient to apply the curl operator ( × ) twice on Equation (5) and, using Equation (4), consider only the horizontal x-component of the result to yield,
( 1 V a t + 1 ) 2 u R a ω 0 [ 1 + δ sin ( ω t ) ] V 2 T + R a g 2 T x z = 0 ,  
where V 2 = 2 / y 2 + 2 / z 2 .

3. Linear Stability Analysis

Assuming small perturbations around the basic solution of the form T = T B + T and w = w B + w , the linearizing Equations (4)–(7) yield the following equations:
( 1 V a t + 1 ) 2 u R a ω 0 [ 1 + δ sin ( ω t ) ] V 2 T + R a g 2 T x z = 0 ,  
d T B d z u = ( 2 T T t w B T z ) .
In Equations (8) and (9), the primes denote the perturbed quantities. Here, in any event, the cases of large and small to moderate Vadasz numbers are considered separately to demonstrate the effect of the Vadasz number on the stability of convection.

3.1. Large Vadasz Numbers

When using typical engineering working fluids in systems with low porosity, large Vadasz numbers result. For large Vadasz numbers, Equations (8) and (9) may be written as follows:
2 u R a ω 0 [ 1 + δ sin ( ω t ) ] V 2 T + R a g 2 T x z = 0 ,
d T B d z u = ( 2 T T t w B T z ) .
Substituting Equation (11) in Equation (10) and simplifying yields the following equation for T :
2 ( 2 T T t w B T z d T B d z ) R a ω 0 [ 1 + δ sin ( ω t ) ] V 2 T + R a g 2 T x z = 0 .  
Assuming an expansion into normal modes in the y- and z-directions of the form yields the following:
T = exp ( i s y y + i s z z ) k = 1 N a k ( t ) sin ( k π x ) .
which satisfies the boundary conditions T = 0 and w = 0 at x = 0 and x = 1 . In Equation (10). s 2 = s y 2 + s z 2 , and k = 1, 2, 3, … N. It was shown in [11] that Equation (12) only yields real solutions when s z = 0 , thus indicating that gravity plays a passive role in the convection threshold criteria. Taking note of the findings in [11], substituting Equation (13) into Equation (12), and then multiplying the result by orthogonal functions and integrating over x [ 0 , 1 ] yields the following system of equations:
k = 1 N { ( k 2 π 2 + s 2 + t ) ( 1 2 + ( k 2 π 2 + s 2 ) ( D 2 ( D 1 ) + 1 6 1 π 2 ( k + l ) 2 ) ) R a ω 0 s 2 ( 1 + δ sin ( ω t ) ) ( 1 2 ( D 1 ) ( D 2 D 2 + 1 2 ) + 1 8 + 3 2 ( 1 2 D ) π 2 ( l + k ) 2 ) } ξ l k a k + k = 1 N { ( k 2 π 2 + s 2 ) ( k 2 π 2 + s 2 + t ) 2 k l ( 2 D 1 ) π 2 ( k 2 l 2 ) 2 R a ω 0 s 2 ( 1 + δ sin ( ω t ) ) ( 6 k l ( D 1 ) 2 π 2 3 k l ( 1 2 D ) π 2 6 π 4 ( 1 ( l k ) 4 1 ( l + k ) 4 ) ) } ξ l + k , 2 p 1 a k = 0 .
where a is the vibration amplitude, D = 3 / 2 , ξ is the Kronecker delta function defined as ξ l k = 1 when l = k, and ξ l + k , 2 p 1 is 1 when ( l + k ) is odd, and zero otherwise. As can be observed, the system shown in Equation (14) becomes rather complicated when solved to higher ranks of N. A numerical solver may be used to solve Equation (14) to higher ranks, but that is outside the scope of the research. Useful information may be drawn from rank N = 1. For the current study, we consider rank N = 1, and Equation (14) may be presented as follows:
d a 1 d t + F 0 s 2 [ ( R a ω 0 R 0 ) R a ω 0 δ sin ( ω t ) ] a 1 = 0 ,
where
F 0 = 5 8 3 4 π 2 [ ( s 2 + π 2 ) ( 13 24 1 4 π 2 ) 1 ] ,  
and
R 0 = ( s 2 + π 2 ) [ ( s 2 + π 2 ) ( 13 24 1 4 π 2 ) 1 ] s 2 ( 5 8 3 4 π 2 ) ,
where Equation (17), incidentally, represents the Rayleigh number corresponding to internal heat generation with no vibration. Solving the first order differential equation in Equation (15) yields the following:
a 1 = b 0 exp { F 0 s 2 [ σ t R a ω 0 δ ω cos ( ω t ) ] } ,
where b 0 is an integration constant and σ = R a ω 0 R 0 . For stationary convection, σ = 0 , therefore,
R a ω 0 = R 0 ,
and
a 1 = b 0 exp [ F 0 s 2 R a 0 κ F r ω cos ( ω t ) ] .
In Equation (19), R 0 represents the Rayleigh number when the vibration frequency is zero. From Equation (19) it is clear that for large Vadasz numbers, the stability criteria for convection are the same as that for the case of no vibration. This simply means that vibration has no effect on the stability criteria for a porous layer with internal heat generation when the Vadasz number, is large. In engineering practice, if the centrifugal jitter effect is required to stabilize the convection, then for a given fluid (as per Figure 2), one would need to increase the porosity of the matrix to reduce the Vadasz number. Conversely, if the effects of centrifugal jitter are present and one does not want it to affect the convection threshold, then one may choose a working fluid and porosity combination to ensure that a large Vadasz number, results to ensure that the time derivative is neglected, thus resulting in the criteria presented by Equation (19).

3.2. Small to Moderate Vadasz Numbers

When using non-traditional working fluids, such as liquid metals (Figure 3), small to moderate Vadasz numbers result across the range of porosity values. For small to moderate Vadasz numbers, substituting Equation (9) in Equation (8) and simplifying yields the following equation for T :
( 1 V a t + 1 ) 2 ( 2 T T t w B T z d T B d z ) R a ω 0 [ 1 + δ sin ( ω t ) ] V 2 T + R a g 2 T x z = 0 .
We assume an expansion into normal modes in the y- and z- directions as per Equation (13); substituting in Equation (21) and then multiplying the result by orthogonal functions and integrating over x [ 0 , 1 ] yields the following system of equations:
k = 1 N { ( 1 V a t + 1 ) ( k 2 π 2 + s 2 + t ) ( 1 2 + ( k 2 π 2 + s 2 ) ( D 2 ( D 1 ) + 1 6 1 π 2 ( k + l ) 2 ) ) R a ω 0 s 2 ( 1 + δ sin ( ω t ) ) ( 1 2 ( D 1 ) ( D 2 D 2 + 1 2 ) + 1 8 + 3 2 ( 1 2 D ) π 2 ( l + k ) 2 ) } ξ l k a k + k = 1 N { ( k 2 π 2 + s 2 ) ( 1 V a t + 1 ) ( k 2 π 2 + s 2 + t ) 2 k l ( 2 D 1 ) π 2 ( k 2 l 2 ) 2 R a ω 0 s 2 ( 1 + δ sin ( ω t ) ) ( 6 k l ( D 1 ) 2 π 2 3 k l ( 1 2 D ) π 2 6 π 4 ( 1 ( l k ) 4 1 ( l + k ) 4 ) ) } ξ l + k , 2 p 1 a k = 0 .
As before, useful information may be drawn from rank N = 1, which may be presented as follows:
d 2 a 1 d t 2 + ( V a + s 2 + π 2 ) d a 1 d t + V a F 0 s 2 [ ( R a ω 0 R 0 ) R a ω 0 δ sin ( ω t ) ] a 1 = 0 ,
where F 0 and R a 0 are as defined before. Equation (23) may be transformed into the Mathieu equation by taking a 1 = e σ τ X 1 ( τ ) to ensure that the resulting equation may be presented as follows:
X 1 + [ A + 2 Q cos ( 2 τ ) ] X 1 = 0 .
In Equation (24), the coefficients A and 2Q are defined for stationary convection as follows:
A = 4 V a s 2 ( 5 8 3 4 π 2 ) ω 2 [ ( s 2 + π 2 ) ( 13 24 1 4 π 2 ) 1 ] ( R a ω 0 R 0 ) ,
2 Q = 4 V a s 2 ( 5 8 3 4 π 2 ) R a ω 0 δ ω 2 [ ( s 2 + π 2 ) ( 13 24 1 4 π 2 ) 1 ] .
In previous studies [14,20], authors have related A and Q in Equations (25) and (26) via an indirect numerical method. Whilst that method recovers both the synchronous and subharmonic modes, it is quite cumbersome. In this paper, we propose the following asymptotic expansions for low amplitude and high frequency vibrations in an attempt to link A and Q:
X 1 = X 0 + Q X 1 + Q 2 X 2 + ,
A = A 0 + Q A 1 + Q 2 A 2 + ,
where Q, as defined in Equation (26), assumes small values for high vibration frequencies. The methodology used in the derivation of Equations (29) and (30) involves the use of simple asymptotic expansions, i.e., Equations (27) and (28), that are then applied on the Mathieu equation, Equation (24). The resulting equations are solved to the various orders in Q, noting that one needs to apply the solvability condition to the resulting differential equations. The result shown in Equation (30), therefore, is the solvability condition. The relation between A and Q are then used in conjunction with their definitions in Equations (25) and (26) to yield the characteristic equation for the Rayleigh number. In the derivation of Equations (29) and (30), the coefficients of the exponential growth terms are set to zero and the constant, c 0 > 0 , is an integration constant (and a real number). The resulting solutions to rank N = 1 and the solvability condition may be represented as follows:
a 1 = X 1 = c 0 [ 1 Q 2 cos ( 2 τ ) ] ,
A = 1 2 Q 2 .
Using Equations (25) and (26) and substituting in Equation (30) yields the following characteristic equation for the Rayleigh number:
V a s 2 ( 5 16 3 8 π 2 ) ω 2 [ ( s 2 + π 2 ) ( 13 24 1 4 π 2 ) 1 ] ( κ F r ω ) 2 R a ω 0 2 R a ω 0 + R a 0 = 0 .
From Equation (31), one may observe that for the case when the vibration frequency approaches zero, i.e., ω 0 , the Rayleigh number approaches the definition, R a ω 0 R a 0 . Equation (31) is solved for both zero and non-zero vibration frequencies. The Vadasz number and the Rayleigh numbers are henceforth scaled as γ = V a / π 2 , R = R a ω 0 / π 2 , and R 0 = R a 0 / π 2 .

4. Results and Discussion

When the vibration is zero, the critical Rayleigh number for internal heat generation to rank N = 1 for the horizontal porous layer subjected to gravity (and no rotation) is R 0 , c r . = 3.819 π 2 , whilst the critical wavenumber is s c r . = 0.49 π . In contrast, for the vertical rotating porous layer, placed far from the axis of rotation and subjected to internal heat generation, the critical Rayleigh number is R 0 , c r . = 3.391 π 2 and the critical wavenumber is s c r . = 0.81 π . These critical Rayleigh and wavenumbers are then compared to the critical Rayleigh number and wavenumber for the classical Benard convection stated as R 0 , c r . = 4 π 2 and s c r . = π . Figure 4 shows the three mentioned cases for comparison purposes, and it is clear from the results that the critical Rayleigh number for the rotating porous layer presents the most unstable case in terms of the convection threshold.
Previously, the authors [24] showed that for the horizontal porous layer, subjected to internal heat generation and g-jitter, the critical Rayleigh number is significantly lower than the modified internal Rayleigh number of 470 ( 47.621 π 2 ) predicted in [17]. When comparing the critical Rayleigh number for the rotating porous layer, placed far from the axis of rotation and subjected to internal heat generation, one observes that although comparable to [24], the magnitude is still significantly smaller than [17].
Figure 5 shows the critical Rayleigh number versus the vibration frequency for selected values of the scaled Vadasz number, for the vertical porous layer placed far away from the axis of rotation and subjected to internal heat generation (black lines). Superimposed on Figure 5 are the data generated in [24] for a horizontal porous layer subjected to gravity (no rotation) and internal heat generation (gray lines). It is observed in Figure 5 that each of the curves end abruptly at specific frequency asymptotes, beyond which there are no more real solutions. These points also denote the end of the synchronous frequencies and the start of the subharmonic frequencies. As stated initially, the current solution approach does not recover solutions corresponding to subharmonic frequencies. The solution approach and methodology presented in [13,14] recovers the transition point from synchronous to subharmonic solutions. The impact of increasing the Vadasz number causes the onset of the subharmonic solutions at lower frequency values, as indicated in Figure 5. At very low-scaled Vadasz numbers, circa γ = 0.01 , the onset of the subharmonic frequencies occurs at very large vibration frequencies.
The results show that the onset of convection for the vertical rotating porous layer placed far from the axis of rotation and subjected to centrifugal jitter occurs sooner across the bandwidth of vibration frequencies for a particular Vadasz number. As observed in Figure 5, increasing the Vadasz number stabilizes the convection.

5. Conclusions

The results for internal heat generation within a rotating porous layer placed far from the axis of rotation and subjected to centrifugal jitter is presented. Both the typical fluids used on engineering applications and heavy liquid metals were analyzed in the current study. The results of this analysis showed that for the typical fluids used in engineering applications, when the porosity is medium to low, the Vadasz number is large, thus neglecting the effects of the time derivative in the momentum equation. In principle, for large Vadasz numbers, the presence of centrifugal jitter has no impact on the stability criteria for convection. However, even with typical engineering working fluids, it was established that a small Vadasz number could result when the porosity is high. In the instance of high porosity media, the Brinkman model may need to be considered when determining the stability criteria. For the instance when heavy liquid metals are considered, it was established that the Vadasz number is many orders of magnitude smaller across the entire range of porosity values when compared to the case of typical engineering working fluids. For large Vadasz numbers, it was shown that the convection stability criteria are independent of the centrifugal jitter.
For small to moderate Vadasz numbers, the results show that in comparison to the case of the horizontal (non-rotating porous layer) subjected to internal heat generation [24], the convection occurs sooner (i.e., at a lower Rayleigh number). However, although these two cases are comparable, the magnitude of the Rayleigh numbers are still significantly smaller than those predicted in [17].

Funding

This research received no external funding.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

The author dedicates this work to his second daughter Vinayia-Chrisinthi Govender, for her interest and questions on the topics that I research.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Horton, C.W.; Rogers, F.T. Convection currents in a porous medium. J. Appl. Phys. 1945, 16, 367–370. [Google Scholar] [CrossRef]
  2. Lapwood, E.R. Convection of a fluid in a porous medium. Math. Proc. Camb. Philos. Soc. 1948, 44, 508–521. [Google Scholar] [CrossRef]
  3. Wooding, R.A. Rayleigh instability of a thermal boundary layer in a flow through a porous medium. J. Fluid Mech. 1960, 9, 183–192. [Google Scholar] [CrossRef]
  4. Palm, E.; Weber, J.E.; Kvernvold, O. On steady convection in a porous medium. J. Fluid Mech. 1972, 54, 153–161. [Google Scholar] [CrossRef]
  5. Vadasz, P. Natural convection in porous media induced by the centrifugal body force—The solution for small aspect ratio. J. Energy Resour. Technol. 1992, 114, 250–254. [Google Scholar] [CrossRef]
  6. Vadasz, P. Centrifugally generated free convection in a rotating porous box. Int. J. Heat Mass Transf. 1994, 37, 2399–2404. [Google Scholar] [CrossRef]
  7. Vadasz, P. Stability of free convection in a narrow porous layer subject to rotation. Int. Commun. Heat Mass Transf. 1994, 21, 881–890. [Google Scholar] [CrossRef]
  8. Vadasz, P. Stability of free convection in a rotating porous layer distant from the axis of rotation. Transp. Porous Media 1996, 23, 153–173. [Google Scholar] [CrossRef]
  9. Vadasz, P. Coriolis effect on gravity driven convection on a rotating porous layer heated from below. J. Fluid Mech. 1998, 376, 351–375. [Google Scholar] [CrossRef]
  10. Govender, S. Oscillatory convection induced by gravity and centrifugal forces in a rotating porous layer distant from the axis of rotation. Int. J. Eng. Sci. 2003, 41, 539–545. [Google Scholar] [CrossRef]
  11. Govender, S. Coriolis effect on the linear stability of convection in a porous layer placed far away from the axis of rotation. Transp. Porous Media 2003, 51, 315–326. [Google Scholar] [CrossRef]
  12. Straughan, B. A sharp nonlinear stability threshold in rotating porous convection. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2001, 457, 87–93. [Google Scholar] [CrossRef]
  13. Govender, S. Stability of convection in a gravity modulated porous layer heated from below. Transp. Porous Media 2004, 57, 103–112. [Google Scholar] [CrossRef]
  14. Govender, S. Linear stability of gravity modulated convection in a porous layer heated from below: Transition from synchronous to subharmonic solutions. Transp. Porous Media 2005, 59, 215–225. [Google Scholar] [CrossRef]
  15. Govender, S. Stability of gravity driven convection in a cylindrical porous layer heated from below. Transp. Porous Media 2006, 63, 489–502. [Google Scholar] [CrossRef]
  16. Govender, S. Vibration effects on convection in a rotating fluid saturated porous layer distant from the axis of rotation. Int. J. Heat Mass Transf. 2019, 141, 112–115. [Google Scholar] [CrossRef]
  17. Gasser, R.D.; Kazimi, M.S. Onset of convection in a porous medium with internal heat generation. J. Heat Transf. 1976, 20, 49–54. [Google Scholar] [CrossRef]
  18. Tveitereid, M. Thermal convection in a horizontal porous layer with internal heat sources. Int. J. Heat Mass Transf. 1977, 20, 1045–1050. [Google Scholar] [CrossRef]
  19. Bhadauria, B.S.; Hashim, I.; Siddheshwar, P.G. Study of heat transport in a porous medium under g-gitter and internal heating effects. Transp. Porous Media 2013, 96, 21–37. [Google Scholar] [CrossRef] [Green Version]
  20. Yadav, D.; Wang, J.; Lee, J. Onset of Darcy-Brinkman convection in a rotating porous layer induced by purely internal heating. J. Porous Media 2017, 20, 691–706. [Google Scholar] [CrossRef]
  21. Mahajan, A.; Nandal, R. Stability of an anisotropic porous layer with internal heat source and Brinkman effects. Spec. Top. Rev. Porous Media 2019, 10, 65–87. [Google Scholar] [CrossRef]
  22. Storesletten, L.; Rees, D.A.S. Onset of convection in an inclined anisotropic porous layer with internal heat generation. Fluids 2019, 4, 75. [Google Scholar] [CrossRef] [Green Version]
  23. Nield, D.; Bejan, A. Convection in Porous Media, 4th ed.; Springer: New York, NY, USA, 2017; pp. 1–983. [Google Scholar]
  24. Govender, S. Vadasz number effects on convection in a horizontal porous layer subjected to internal heat generation and g-jitter. Fluids 2020, 5, 124. [Google Scholar] [CrossRef]
  25. Saravanan, S.; Brindha, D. Global stability of centrifugal filtration convection. J. Math. Anal. Appl. 2010, 367, 116–128. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Vertical rotating porous layer, placed far from axis of rotation, and subjected to centrifugal jitter and internal heat generation. T H * and T C * are the dimensional hot wall temperature and cold wall temperatures, t * is the dimensional time variable, and Ω * and ω * are the rotational speed and vibration frequency. The length the porous layer is L * and the offset distance is X 0 * whilst the gravity is defined as g * .
Figure 1. Vertical rotating porous layer, placed far from axis of rotation, and subjected to centrifugal jitter and internal heat generation. T H * and T C * are the dimensional hot wall temperature and cold wall temperatures, t * is the dimensional time variable, and Ω * and ω * are the rotational speed and vibration frequency. The length the porous layer is L * and the offset distance is X 0 * whilst the gravity is defined as g * .
Physics 03 00044 g001
Figure 2. Vadasz number illustrated for typical liquids, used in engineering applications. Here, ϕ is the porosity.
Figure 2. Vadasz number illustrated for typical liquids, used in engineering applications. Here, ϕ is the porosity.
Physics 03 00044 g002
Figure 3. Vadasz number, illustrated for heavy liquid metals. Here, ϕ is the porosity.
Figure 3. Vadasz number, illustrated for heavy liquid metals. Here, ϕ is the porosity.
Physics 03 00044 g003
Figure 4. Comparison of characteristic Rayleigh number curves for the horizontal non-rotating porous layer [24], and vertical rotating porous layer, with internal heat generation, to the classical Benard convection characteristic curves. α denotes the wavenumber.
Figure 4. Comparison of characteristic Rayleigh number curves for the horizontal non-rotating porous layer [24], and vertical rotating porous layer, with internal heat generation, to the classical Benard convection characteristic curves. α denotes the wavenumber.
Physics 03 00044 g004
Figure 5. Critical Rayleigh number versus vibration frequency for selected values of Vadasz number for a vertical porous layer placed far away from the axis of rotation (black lines). The results for the non-rotating horizontal porous layer with internal heat generation [24] (gray lines) superimposed.
Figure 5. Critical Rayleigh number versus vibration frequency for selected values of Vadasz number for a vertical porous layer placed far away from the axis of rotation (black lines). The results for the non-rotating horizontal porous layer with internal heat generation [24] (gray lines) superimposed.
Physics 03 00044 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Govender, S. Vadasz Number Effects on Convection in a Vertical Rotating Porous Layer, Placed Far from Axis of Rotation, and Subjected to Internal Heat Generation and Centrifugal Jitter. Physics 2021, 3, 728-738. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3030044

AMA Style

Govender S. Vadasz Number Effects on Convection in a Vertical Rotating Porous Layer, Placed Far from Axis of Rotation, and Subjected to Internal Heat Generation and Centrifugal Jitter. Physics. 2021; 3(3):728-738. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3030044

Chicago/Turabian Style

Govender, Saneshan. 2021. "Vadasz Number Effects on Convection in a Vertical Rotating Porous Layer, Placed Far from Axis of Rotation, and Subjected to Internal Heat Generation and Centrifugal Jitter" Physics 3, no. 3: 728-738. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3030044

Article Metrics

Back to TopTop