1. Introduction
The analysis of instability of a steady convective flow generated by internal heat sources is important for many problems in science and engineering. The stability of a steady convective flow in an annulus caused by heat sources with constant density is analyzed in [
1]. The case of heat sources of non-uniform density is considered in [
2,
3]. It is shown in these papers that for annual flow, two types of instability exist: (a) shear instability and (b) buoyant instability, discovered earlier in [
4,
5] for the vertical planar layer and circular pipe, respectively. Recent interest in biomass thermal conversion [
6,
7] stimulated research in the stability of flows driven by nonlinear heat sources. The base flow solution in this case cannot be found analytically as in many classical problems in hydrodynamic stability (see [
8], for example). Bifurcation theory [
9] can be used to determine the number of solutions and the region in the parameter space where solutions exist. In a recent paper [
10], bifurcation theory is used to analyze nonlinear boundary value problem in an annulus. It is found that for each radius ratio, between the radii, there exists the value of the Frank–Kamenetskii parameter
(proportional to the intensity of the chemical reaction) such that there are two solutions in the domain
, one solution for
and no solutions in the region of
.
The stability of viscous flow between two concentric cylinders with a radial flow through permeable porous walls of the annulus is analyzed in [
11]. Such models (in the presence of rotation) are motivated by applications to dynamic filtration devices and vortex flow reactors [
12,
13]. The stability of such flows is investigated in [
14,
15,
16]. In the case of non-isothermal conditions, the presence of radial flow can alter the perturbation dynamics. The stability of a combined flow caused by internal heat sources with constant density and a radial inflow or outflow between the permeable walls of the annulus is analyzed in [
17], where it is shown that radial flow can stabilize or destabilize the convective flow. In a recent study, the linear stability of a flow in porous medium with permeable boundaries in a vertical layer of annular cross-section is analyzed by [
18] for the isothermal case and in [
19] for the case where the walls of the layer are maintained at different temperatures. It is shown that for non-isothermal case the flow becomes more and more unstable as the aspect ratio increases.
Linear stability of a steady convective flow in an annulus generated by non-homogeneous heat sources is studied in [
20]. The effect of nonlinear internal heat sources on the stability boundary of a vertical flow in an annulus is investigated in [
21] for both asymmetric and axisymmetric perturbations for wide range of radius ratios, Prandtl numbers and Frank–Kamenetskii parameters (proportional to the intensity of a chemical reaction). Such models are used to investigate processes in biomass thermal conversion. Different factors (as shown experimentally in [
22]) affect the process: external electric field, the degree of swirl of the flow, and convection in the chamber. Complex physical processes are analyzed in nature and engineering, using three basic approaches: (a) experimental analysis, (b) numerical modeling ([
23]), and (c) stability analysis. In the present paper, we consider linear stability of a steady convective flow in an annulus caused by nonlinear heat sources. In addition, there is a radially inward or outward flow through porous walls of the annulus. In this case, we are looking for flow instability since it results in more intensive mixing and (hopefully) more efficient energy conversion. The stability of a combined flow (vertical flow due to heat sources and radial flow through the walls) is analyzed. Wide gaps are not used in the study; therefore, the linear stability is analyzed with respect to axisymmetic perturbations. The corresponding nonlinear boundary value problem for the base flow is solved numerically. Using the method of normal modes, we reduce the linear stability problem to the solution of a system of ordinary differential equations with variable coefficients. The collocation method based on Chebyshev polynomials is used to discretize the problem. Calculations show that both inward and outward radial flows stabilize the convective flow in the vertical direction. On the other hand, the Frank–Kamenetskii parameter has a destabilizing influence on the flow. It is shown that a second minimum appears on the marginal stability curve as the Reynolds number increases. For larger values of the Reynolds number, the second minimum disappears.
2. Mathematical Formulation of the Problem
Suppose that a viscous incompressible fluid is situated in the region
between two infinitely long concentric cylinders, where
is a system of cylindrical polar coordinates centered at the axes of the cylinders. The cylinders’ walls are maintained at a constant and equal temperature
. Heat is generated inside the annulus, due to the exothermal chemical reaction described by Arrhenius’ law:
where
E is the activation energy,
R is the universal gas constant,
is a constant and
is the absolute temperature. Models with internal heat generation given by (1) are used in practice in order to describe processes during biomass thermal conversion [
20,
21,
23] with the objective to obtain cleaner and more efficient sources of energy. Instability is a desirable phenomenon in this case since it enhances fluid mixing, which, in turn, leads to more efficient energy conversion. One of the factors that affects the stability characteristics is the radial flow in the transverse direction. We assume that there is a radial inward or outward flow through the permeable walls
and
. The following convention is used in the paper: the variables with tildes are dimensional while the variables without tildes are dimensionless.
The problem is described by the system of the Navier–Stokes equations under the Boussinesq approximation. The dimensionless form of the system is as follows:
where
is the velocity of the fluid,
T is the temperature,
p is the pressure and
. The Frank–Kamenetskii transformation is used [
24] to transform the source term in (1). The idea is to expand the exponent in the Taylor series and take into account only the linear terms of the expansion. The accuracy of the Frank–Kamenetskii transformation is analyzed in [
24,
25], where it is shown that for typical values of the parameters in (1), the transformation is quite accurate. The advantage of using it is related to the fact that the source term in (3) is mathematically easier to work with than the term in (1).
The following values are chosen as the measures, respectively, of length, , time, , velocity, , temperature, , and pressure, . Here, is the density of the fluid, g is the acceleration due to gravity, is the coefficient of the thermal expansion and is the viscosity of the fluid. In addition, we introduce the notations , and .
Equations (2)–(4) admit a steady solution of the following form:
The radial inflow or outflow through porous cylinders
and
is described by the function
. It follows from (4) that
satisfies the following equation:
Hence,
where
D is an arbitrary constant. The boundary conditions are as follows:
The following four dimensionless parameters are used to describe the problem: the Grashof number,
, the Prandtl number,
, the Frank–Kamenetskii parameter,
, and the radial Reynolds number
, where the radial component of the base flow in dimensionless form is given by the following:
Here, U is the velocity of the fluid at . It follows from (6) that the velocity at the inner boundary (the first boundary condition in (6)) is fixed. In this case, it is required from continuity (Equation (4)) that the second condition in (6) must be satisfied. Radial outflow ) corresponds to the case where fluid enters the domain through the inner boundary and leaves it at . Radial inflow corresponds to the flow in the opposite direction.
Substituting (6) and (7) into (2)–(4), we obtain the system of equations describing the base flow:
where
. The function
is not determined since it will not be used in sequel.
The boundary conditions are the following:
It is assumed that the annulus is closed so that the total fluid flux through the cross-section of the annulus is zero:
The nonlinear boundary value problem of (8)–(11) is solved, using Matlab routine bvp4c for different values of the parameters,
F,
, and
. The graphs of the base flow velocity and temperature distributions are shown in
Figure 1 and
Figure 2, respectively, for three values of the Reynolds number. The parameters
R and
F are fixed at
and
, respectively.
Instability is expected here since the velocity profiles contain inflection points. For , there is one flow upstream in the middle portion of the channel and two flows downstream near the boundaries. As the Reynolds number increases, the intensity of the downstream flow near the outer boundary decreases. Velocity gradients also become smaller as the Reynolds number grows. Thus, it is expected that radial flow will have a stabilizing influence on the stability boundary. This fact is confirmed later by numerical calculations. The base flow temperature distribution becomes more asymmetric as the Reynolds number grows.
3. Numerical Results
Consider a perturbed motion of the following form:
where
,
, and
are small unsteady perturbations and
is the unit vector in the positive
r-direction. Using a standard linearization procedure, we represent the perturbed quantities in the form of axisymmetric normal modes as follows:
where
k is the wave number,
is a complex eigenvalue and
. The flow (8)–(11) is linearly stable if all
and is unstable if at least one
. The flow (8)–(11) is marginally stable if one eigenvalue has
, while all other eigenvalues have positive real parts. Eliminating the pressure and longitudinal velocity perturbations from the linearized system, we obtain the following system of ordinary differential equations:
The boundary conditions are as follows:
The eigenvalue problem of (15)–(17) is solved numerically, using the collocation method based on Chebyshev polynomials. In particular, the interval
is transformed to the interval
by means of the following transformation:
where
. The functions
u and
(in terms of the variable
) are approximated as follows:
where
is the Chebyshev polynomial of the first kind of order
m. The collocation points are the following:
In order to estimate the number of collocation points needed for accurate determination of the Grashof number, we performed calculations for one set of parameters, namely,
,
,
,
,
, and different number of collocation points
N. The results are shown in
Table 1. It is seen from the table that 50 collocation points are sufficient for the calculation of
, accurate to within 6 decimal places after the decimal point. Similar calculations are performed for other sets of parameters. The calculations show that it is sufficient to use
for all cases considered in the paper.
All stability characteristics are calculated below for the case
and
. It is shown in [
1] that only axisymmetric perturbations are the most unstable for small and moderate gaps, while the first asymmetric mode is the most unstable for the range
. This is the reason we restrict ourselves with axisymmetric perturbations.
Marginal stability curves for
and different values of
F and
are shown in
Figure 3. The dots on the curves correspond to the calculated points, while solid lines are the interpolating curves.The flow is linearly stable in the regions below the curves and unstable in the regions above the curves. On the marginal stability curve, one eigenvalue has zero real part, while all other eigenvalues have positive real parts. The point of an absolute minimum on each marginal stability curve corresponds to the critical values of the parameters
and
k, denoted by
and
, respectively. Thus, the flow is linearly stable for all
k if
.
Several conclusions can be drawn from the graphs in
Figure 3. First, the increase in
F has a destabilizing influence on the flow (the critical Grashof numbers decrease as
F increases). Second, for each fixed
F, there is a continuous transformation of marginal stability curves as the Reynolds number increases. For
(no radial cross-flow), the marginal stability curve has one minimum. As the Reynolds number increases (see the range
, the second minimum appears on the curves. The increase in the Reynolds number (up to
) leads to a shift of the minimum point to the region of smaller
k. Whether the second minimum is the global minimum depends on the value of the Reynolds number.
It is seen from
Figure 3 that the minimum corresponding to smaller
k is the global minimum in the range
, where the value of
depends on
F,
R, and
. Calculations show that for the case
,
, and
, we have
The corresponding marginal stability curve is plotted in
Figure 4, where the presence of two equal minima is clearly seen.
The second minimum, which is seen in
Figure 3, seems to disappear for higher
. In order to check what happens for larger values of
, we perform calculations for
in the range
. The results are shown in
Figure 4.
It is seen from
Figure 5 that for large
, the marginal stability curves have the same shape as for
(with one minimum). In addition, the increase in
stabilizes the flow (the critical value
increases as
grows).
The effect of both outward and inward radial flows is investigated further.
Figure 6 plots the marginal stability curves for the case
,
and both negative and positive Reynolds numbers. It is seen from the graph that the outward radial flow (positive
) is less stable than the inward radial flow (negative
).
The critical Grashof numbers versus
are plotted in
Figure 7. The Reynolds number has a stabilizing effect on the convective flow since
is increasing as
grows. However, the rate of growth is not the same. The critical Grashof numbers increase faster in the range
where instability is associated with small wave number perturbations. Then, there is a relatively small growth in
for
, and then
increases faster.
Critical wave numbers versus
are shown in
Figure 8. The finite jump at
is associated with the transition to perturbations with larger
k as indicated in
Figure 4.
Figure 9 plots the critical Grashof numbers versus positive and negative values of
in the range
for four different values of the Frank–Kamenetskii parameter
F, namely,
and 0.7. The values of
R and
are fixed at
and
, respectively. Several conclusions can be drawn from the graphs in
Figure 9. First, both inward and outward radial flows (negative and positive Reynolds numbers) have a stabilizing influence on the stability boundary. Second, for each
F, three different intervals characterizing the rate of increase in the Grashof number with respect to
can be identified. The first interval (approximately from
to
) is associated with relatively strong stabilization of the base flow. The second interval (
) appears right after the transition to a larger wave number takes place (see
Figure 4 and
Figure 5 for details). The critical Grashof numbers continue to grow, but at a lower rate. However, as
increases further, stabilization becomes stronger (the rate of increase in
with respect to
increases). A similar situation takes place around the value
. Here, again, transition to a larger wave number takes place. As a result, the rate of growth of
with respect to
decreases, but then increases again as the Reynolds number becomes more and more negative. Third, the rate of increase in
with respect to
is not the same for positive and negative Reynolds numbers. Fourth, base flow stabilization also depends on the value of the Frank–Kamenetskii parameter
F. Stabilization is much stronger for small
F (see the graph for
in
Figure 9) and less pronounced for larger
F. Fifth, as
F increases, the critical Grashof number approaches zero. The last conclusion is consistent with the fact that a steady convective flow in the vertical direction generated by internal heat sources exists only in the range
(see [
10] for the case
), where
depends on
R and
, and there is no steady solution for
.