Next Article in Journal
Data-Targeted Prior Distribution for Variational AutoEncoder
Next Article in Special Issue
Numerical and Analytical Studies of Soret-Driven Convection Flow Inside an Annular Horizontal Porous Cavity
Previous Article in Journal
Semi-Implicit Finite Volume Procedure for Compositional Subsurface Flow Simulation in Highly Anisotropic Porous Media
Previous Article in Special Issue
Analysis of Temperature Anomalies during Thermal Monitoring of Frozen Wall Formation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Linear Stability of a Steady Convective Flow between Permeable Cylinders

Institute of Applied Mathematics, Riga Technical University, Zunda Embankment 10, LV-1048 Riga, Latvia
*
Author to whom correspondence should be addressed.
Fluids 2021, 6(10), 342; https://doi.org/10.3390/fluids6100342
Submission received: 24 August 2021 / Revised: 23 September 2021 / Accepted: 24 September 2021 / Published: 28 September 2021
(This article belongs to the Special Issue Convection in Fluid and Porous Media)

Abstract

:
Linear stability analysis of a steady convective flow in a tall vertical annulus caused by nonlinear heat sources is conducted in the paper. Heat sources are generated as a result of a chemical reaction. The effect of radial cross-flow through permeable porous walls of the annulus is analyzed. The problem is relevant to biomass thermal conversion. The base flow solution is obtained by solving nonlinear boundary value problem. Linear stability analysis is performed, using collocation method. The calculations show that radial inward or outward flow has a stabilizing effect on the flow, while the increase in the Frank–Kamenetskii parameter (proportional to the intensity of the chemical reaction) destabilizes the flow. The increase in the Reynolds number based on the radial velocity leads to the appearance of the second minimum on the marginal stability curves. The rate of increase in the critical Grashof number with respect to the Reynolds number is different for inward and outward radial flows.

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 F * (proportional to the intensity of the chemical reaction) such that there are two solutions in the domain 0 < F < F * , one solution for F = F * and no solutions in the region of F > F * .
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 D = { R 1 < r ˜ < R 2 , 0 φ < 2 π , < z ˜ < + } between two infinitely long concentric cylinders, where ( r ˜ , φ , z ˜ ) 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 θ 0 . Heat is generated inside the annulus, due to the exothermal chemical reaction described by Arrhenius’ law:
Q = Q 0 e E / ( R T ˜ ) ,
where E is the activation energy, R is the universal gas constant, Q 0 is a constant and T ˜ 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 r ˜ = R 1 and r ˜ = R 2 . 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:
v t + G r ( v · ) v = p + Δ v + T e k ,
T t + G r v · T = 1 P r Δ T + F P r e T ,
· v = 0 ,
where v is the velocity of the fluid, T is the temperature, p is the pressure and e k = ( 0 , 0 , 1 ) . 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, h = ( R 2 R 1 ) / 2 , time, h 2 / ν , velocity, g β h 2 R θ 0 2 / ( ν E ) , temperature, R θ 0 2 / E , and pressure, ρ g β h R θ 0 2 / E . 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 η = R 1 / R 2 , r 1 = R 1 / h and r 2 = R 2 / h .
Equations (2)–(4) admit a steady solution of the following form:
v 0 = ( U 0 ( r ) , 0 , W 0 ( r ) ) , T 0 = T 0 ( r ) , p 0 = p 01 ( r ) + p 02 ( z ) .
The radial inflow or outflow through porous cylinders r = r 1 and r = r 2 is described by the function U 0 ( r ) . It follows from (4) that U 0 ( r ) satisfies the following equation:
d U 0 d r + U 0 r = 0 .
Hence,
U 0 = D r ,
where D is an arbitrary constant. The boundary conditions are as follows:
U 0 | r = r 1 = R e G r , U 0 | r = r 2 = R e η G r .
The following four dimensionless parameters are used to describe the problem: the Grashof number, G r = g β R θ 0 2 h 3 / ( ν 2 E ) , the Prandtl number, P r = ν / κ , the Frank–Kamenetskii parameter, F = [ ( Q 0 k 0 E h 2 ) / ( κ R θ 0 2 ) ] exp [ E / ( R θ 0 ) ] , and the radial Reynolds number R e = U h / ν , where the radial component of the base flow in dimensionless form is given by the following:
U 0 ( r ) = R e r G r .
Here, U is the velocity of the fluid at r ˜ = R 1 . 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 ( R e > 0 ) corresponds to the case where fluid enters the domain through the inner boundary r = r 1 and leaves it at r = r 2 . Radial inflow ( R e < 0 ) 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:
W 0 + W 0 r ( 1 R e ) + T 0 = C ,
T 0 + T 0 r ( 1 P r R e ) + F e T 0 = 0 ,
where C = d p 02 / d z . The function p 01 ( r ) is not determined since it will not be used in sequel.
The boundary conditions are the following:
W 0 ( r 1 ) = 0 , W 0 ( r 2 ) = 0 , T 0 ( r 1 ) = 0 , T 0 ( r 2 ) = 0 .
It is assumed that the annulus is closed so that the total fluid flux through the cross-section of the annulus is zero:
r 1 r 2 r W 0 ( r ) d r = 0 .
The nonlinear boundary value problem of (8)–(11) is solved, using Matlab routine bvp4c for different values of the parameters, F, P r , and R e . 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 η = 0.6 and F = 0.5 , respectively.
Instability is expected here since the velocity profiles contain inflection points. For R e = 0 , 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:
v = v + W 0 e k + U 0 e r , T = T + T 0 , p = p + p 0 ,
where v , T , and p are small unsteady perturbations and e r 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:
v ( r , z , t ) = u ( r ) exp ( i k z λ t ) ,
T ( r , z , t ) = θ ( r ) exp ( i k z λ t ) ,
where k is the wave number, λ = λ r + i λ i is a complex eigenvalue and u = ( u ( r ) , 0 , w ( r ) ) . The flow (8)–(11) is linearly stable if all λ r > 0 and is unstable if at least one λ r < 0 . The flow (8)–(11) is marginally stable if one eigenvalue has λ r = 0 , 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:
u ( 4 ) + 2 r U 0 G r u 3 r 2 + 2 k 2 + U 0 G r U 0 G r r i k W 0 G r u + 3 r 3 2 k 2 r U 0 G r r + 2 U 0 G r r 2 i k G r W 0 r + k 2 G r U 0 u + ( 2 k 2 r 2 3 r 4 + k 4 + U 0 G r r 2 2 G r U 0 r 3 i k G r W 0 r + i k G r W 0 r 2 + k 2 G r U 0 + i k 3 G r W 0 + i k G r W 0 ) u i k θ = λ u + u r u r 2 k 2 u ,
θ + θ r k 2 θ + F exp ( T 0 ) θ P r G r ( u T 0 + U 0 θ + i k W 0 θ ) = λ P r θ .
The boundary conditions are as follows:
u ( r 1 ) = u ( r 2 ) = 0 , u ( r 1 ) = u ( r 2 ) = 0 , θ ( r 1 ) = θ ( r 2 ) = 0 .
The eigenvalue problem of (15)–(17) is solved numerically, using the collocation method based on Chebyshev polynomials. In particular, the interval [ r 1 , r 2 ] is transformed to the interval [ 1 , 1 ] by means of the following transformation:
r = r 2 r 1 2 ξ + r 2 + r 1 2 ,
where ξ [ 1 , 1 ] . The functions u and θ (in terms of the variable ξ ) are approximated as follows:
u ( ξ ) = m = 0 N a m ( 1 ξ 2 ) 2 T m ( ξ ) , θ ( ξ ) = m = 0 N b m ( 1 ξ 2 ) T m ( ξ ) ,
where T m ( ξ ) = cos m arccos ( ξ ) is the Chebyshev polynomial of the first kind of order m. The collocation points are the following:
ξ j = cos π j N , j = 0 , 1 , , N
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, k = 1 , P r = 2 , R = 0.6 , R e = 2 , F = 0.7 , 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 G r , 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 N = 60 for all cases considered in the paper.
All stability characteristics are calculated below for the case η = 0.6 and η = 0.7 . 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 0 < η < 0.3 . This is the reason we restrict ourselves with axisymmetric perturbations.
Marginal stability curves for η = 0.6 and different values of F and R e 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 G r and k, denoted by G r c and k c , respectively. Thus, the flow is linearly stable for all k if G r < G r c .
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 R e = 0 (no radial cross-flow), the marginal stability curve has one minimum. As the Reynolds number increases (see the range 2 R e 4 ) , the second minimum appears on the curves. The increase in the Reynolds number (up to R e = 4 ) 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 0 < R e < R e * , where the value of R e * depends on F, R, and P r . Calculations show that for the case η = 0.6 , P r = 2 , and F = 0.5 , we have R e * = 3.625 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 R e . In order to check what happens for larger values of R e , we perform calculations for F = 0.5 in the range 5 R e 11 . The results are shown in Figure 4.
It is seen from Figure 5 that for large R e , the marginal stability curves have the same shape as for R e = 0 (with one minimum). In addition, the increase in R e stabilizes the flow (the critical value G r c increases as R e grows).
The effect of both outward and inward radial flows is investigated further. Figure 6 plots the marginal stability curves for the case F = 0.5 , η = 0.7 and both negative and positive Reynolds numbers. It is seen from the graph that the outward radial flow (positive R e ) is less stable than the inward radial flow (negative R e ).
The critical Grashof numbers versus R e are plotted in Figure 7. The Reynolds number has a stabilizing effect on the convective flow since G r c is increasing as R e grows. However, the rate of growth is not the same. The critical Grashof numbers increase faster in the range 0 < R e < R e * = 3.625 where instability is associated with small wave number perturbations. Then, there is a relatively small growth in G r c for R e * < R e < 6 , and then G r c increases faster.
Critical wave numbers versus R e are shown in Figure 8. The finite jump at R e = R e * 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 R e in the range 8 R e 8 for four different values of the Frank–Kamenetskii parameter F, namely, F = 0.1 , 0.3 , 0.5 and 0.7. The values of R and P r are fixed at η = 0.6 and P r = 2 , 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 R e can be identified. The first interval (approximately from R e = 2 to R e = 3.6 ) is associated with relatively strong stabilization of the base flow. The second interval ( R e > 3.6 ) 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 R e increases further, stabilization becomes stronger (the rate of increase in G r c with respect to R e increases). A similar situation takes place around the value R e = 2 . Here, again, transition to a larger wave number takes place. As a result, the rate of growth of G r c with respect to R e decreases, but then increases again as the Reynolds number becomes more and more negative. Third, the rate of increase in G r c with respect to R e 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 F = 0.1 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 0 < F < F * (see [10] for the case R e = 0 ), where F * depends on R and R e , and there is no steady solution for F > F * .

4. Discussion

Linear stability of a steady convective flow caused by nonlinear heat sources in a tall vertical annulus is analyzed in the paper. It is assumed that there is a radial inward or outward flow through the permeable walls of the annulus. The base flow is obtained numerically as a solution of the nonlinear boundary value problem for the system of ordinary differential equations. The linear stability problem is solved by the collocation method for different values of the parameters of the problem. The Prandtl number is fixed at P r = 2 . Analysis of the base flow velocity profiles suggests that the Reynolds number (proportional to radially inward or outward velocity) will stabilize the flow. This fact is confirmed by calculation of the marginal stability curves. It is shown that the second minimum appears in the region of smaller wave numbers as the Reynolds number increases. This minimum may or may not be the global minimum. However, for large Reynolds numbers, the second minimum disappears. The other parameter that affects stability characteristics is the Frank–Kamenetskii parameter (proportional to the intensity of the rate of a chemical reaction). The increase in the Frank–Kamenetskii parameter destabilizes the flow. In addition, there exists the value F * such that steady flow exists only in the range 0 < F < F * . There is no steady flow in the domain F > F * (this fact is denoted as the thermal explosion in the literature). In the region where steady flow exists, the stability characteristics are determined by the concurrence of the two factors: (a) stabilizing effect of the radial flow (both inward and outward) and (b) destabilizing effect of the intensity of the chemical reaction (the Frank–Kamenetskii parameter).
Wide gaps between the walls of the annulus are not considered in the present study (the values of η that are used in the paper are η = 0.6 and η = 0.7 ). As was shown in [1,20], asymmetric perturbations (depending on the angular coordinate φ ) are the most unstable for very large gaps (small values of η ), which is why only axisymmetric perturbations are considered in the present paper. In the future, we plan to investigate the role of the radius ratio on the stability boundary as well as to consider different Prandtl numbers by analyzing both axisymmetric and asymmetric perturbations. In addition, weakly nonlinear theory can be used to analyze the development of instability beyond the threshold. The authors are currently working on these topics.

Author Contributions

Conceptualization, A.K.; methodology, M.Z., A.K. and I.I.; software, M.Z., I.I. and A.K.; validation, M.Z. and I.I.; formal analysis, M.Z. and A.K.; writing—original draft preparation, M.Z. and A.K.; writing—review and editing, M.Z., A.K. and I.I.; supervision, I.I. and A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Latvian Council of Science project No. lzp-2020/1-0076.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors thank the editor and anonymous referees for their constructive comments, which helped us to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kolyshkin, A.A.; Vaillancourt, R. Stability of internally generated thermal convection in a tall vertical annulus. Can. J. Phys. 1991, 69, 743–748. [Google Scholar] [CrossRef]
  2. Kolyshkin, A.A.; Vaillancourt, R. On the stability of convective motion caused by inhomogeneous heat sources. Arab. J. Sci. Eng. 1992, 17, 655–662. [Google Scholar]
  3. Saravan, S.; Kandaswamy, P. Thermal instability of a nonunoformly heat generating annual fluid layer. Int. J. Heat Mass Transf. 2006, 49, 269–280. [Google Scholar] [CrossRef]
  4. Gershuni, G.Z.; Zhukhovitskii, E.M.; Iakimov, A.A. Two kinds of instability of stationary convective motion induced by internal heat sources. J. Appl. Math. Mech. 1973, 37, 544–548. [Google Scholar] [CrossRef]
  5. Rogers, B.B.; Yao, L.S. The importance of Prandtl number for mixed-convection instability. Trans. ASME C 1993, 115, 482–486. [Google Scholar] [CrossRef]
  6. Lewandowski, W.M.; Ryms, M.; Kosakowski, W. Thermal biomass conversionb: A review. Processes 2020, 8, 516. [Google Scholar] [CrossRef]
  7. Yamakawa, C.K.; Quin, F.; Mussatto, S.I. Advances and opportunities in biomass conversion technologies and biorefineries for the development of a bio-based economy. Biomass Bioenergy 2018, 119, 54–60. [Google Scholar] [CrossRef]
  8. Drazin, P.G.; Reid, W.H. Hydrodynamic Stability, 2nd ed.; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  9. Bebernes, J.; Eberly, D. Mathmatical Problems from Combustion Theory; Springer: New York, NY, USA, 1989. [Google Scholar]
  10. Gritsans, A.; Koliškins, A.; Ogorelova, D.; Sadirbajevs, F.; Samuilika, I.; Yermachenko, I. Solutions of Nonlinear Boundary Value Problem with Applications to Biomass Thermal Conversion. In Proceedings of the Engineering for Rural Development, Proceedings, Latvia, Jelgava, 26–28 May 2021. accepted. [Google Scholar]
  11. Deka, R.K.; Takhar, H.S. Hydrodynamic stability of viscous flow between curved porous channel with radial flow. Int. J. Eng. Sci. 2004, 42, 953–966. [Google Scholar] [CrossRef]
  12. Wronsky, S.; Molga, E.; Rudniak, L. Dyniamic filtration in biotechnology. Bioprocess Eng. 1988, 4, 99–104. [Google Scholar] [CrossRef]
  13. Kroner, K.H.; Nissinen, V.; Ziegler, H. Improving dynamic filtration of microbial suspensions. Biotechnology 1987, 5, 921–926. [Google Scholar]
  14. Min, K.; Lueptow, R.M. Hydrodynamic stability of viscous flow between rotating porous cylinders with radial flow. Phys. Fluids 1994, 6, 144–151. [Google Scholar] [CrossRef]
  15. Kolyshkin, A.A.; Vaillancourt, R. Convective instability boundary of Couette flow between rotating porous cylinders with axial and radial flows. Phys. Fluids 1997, 9, 910–918. [Google Scholar] [CrossRef]
  16. Ilin, K.; Morgulis, A. On the stability of the Couette-Taylor flow between rotating porous cylinders with radial flow. Eur. J. Mech.-B/Fluids 2019, 80, 174–186. [Google Scholar] [CrossRef] [Green Version]
  17. Ghidaoui, M.S.; Kolyshkin, A.A.; Vaillancourt, R. Convective stability of mixed flow between permeable cylinders with internal heat generation. In Proceedings of the Eleventh International Heat Transfer Conference, Kyongju, Korea, 23–28 August 1998; Lee, J.S., Cotta, R.M., Hakhoe, T.K., Eds.; Begel House: Kyongju, Korea, 1998; Volume 3, pp. 263–268. [Google Scholar]
  18. Barletta, A.; Celli, M.; Rees, D.A.S. On the stability of parallel flow in a vertical porous layer with annual cross-section. Transp. Porous Media 2020, 134, 491–501. [Google Scholar] [CrossRef]
  19. Barletta, A.; Celli, M.; Rees, D.A.S. Buoyant flow and instability in a vertical porous slab with permeable bloundaries. Int. J. Heat Mass Transf. 2020, 157, 119956. [Google Scholar] [CrossRef]
  20. Kolyshkin, A.; Koliskina, V.; Volodko, I.; Iltins, I. Convective instability of a steady flow in an annulus caused by internal heat generation. Proc. Latv. Acad. Sci. Sect. B 2020, 74, 293–298. [Google Scholar] [CrossRef]
  21. Iltins, I.; Iltina, M.; Kolyshkin, A.; Koliskina, V. Linear stability of a convective flow in an annulus with a nonlinear heat source. JP J. Heat Mass Transf. 2019, 18, 315–329. [Google Scholar] [CrossRef]
  22. Barmina, I.; Purmalis, M.; Valdmanis, R.; Zake, M. Electrodynamic control of the combustion characteristics and energy production. Combust. Sci. Technol. 2016, 188, 190–206. [Google Scholar] [CrossRef]
  23. Kalis, H.; Marinaki, M.; Strautins, U.; Zake, M. On numerical simulation of electromagnetic field effects in the combustion process. Math. Model. Anal. 2018, 23, 327–343. [Google Scholar] [CrossRef] [Green Version]
  24. Frank-Kamenetskii, D.A. Diffusion and Heat Exchange in Chemical Kinetics; Princeton: Princeton, NJ, USA, 1955. [Google Scholar]
  25. White, D., Jr.; Johns, L.E. The Frank-Kamenetskii transformation. Chem. Eng. Sci. 1987, 42, 1849–1851. [Google Scholar] [CrossRef]
Figure 1. Base flow velocity distribution for F = 0.5 and η = 0.6 .
Figure 1. Base flow velocity distribution for F = 0.5 and η = 0.6 .
Fluids 06 00342 g001
Figure 2. Base flow temperarture distribution for F = 0.5 and η = 0.6 .
Figure 2. Base flow temperarture distribution for F = 0.5 and η = 0.6 .
Fluids 06 00342 g002
Figure 3. Marginal stability curves for different values of F and R e .
Figure 3. Marginal stability curves for different values of F and R e .
Fluids 06 00342 g003
Figure 4. Marginal stability curves for F = 0.5 , η = 0.6 and R e = 3.625 .
Figure 4. Marginal stability curves for F = 0.5 , η = 0.6 and R e = 3.625 .
Fluids 06 00342 g004
Figure 5. Marginal stability curves for F = 0.5 and η = 0.6 (the Reynolds numbers are in the range 5 R e 11 ).
Figure 5. Marginal stability curves for F = 0.5 and η = 0.6 (the Reynolds numbers are in the range 5 R e 11 ).
Fluids 06 00342 g005
Figure 6. Marginal stability curves for F = 0.5 and η = 0.7 (the Reynolds number is in the range 3 R e 3 ).
Figure 6. Marginal stability curves for F = 0.5 and η = 0.7 (the Reynolds number is in the range 3 R e 3 ).
Fluids 06 00342 g006
Figure 7. Critical Grashof numbers versus R e for F = 0.5 .
Figure 7. Critical Grashof numbers versus R e for F = 0.5 .
Fluids 06 00342 g007
Figure 8. Critical Grashof numbers versus R e for F = 0.5 .
Figure 8. Critical Grashof numbers versus R e for F = 0.5 .
Fluids 06 00342 g008
Figure 9. Critical Grashof numbers versus R e in the range 8 R e 8 for different values of F.
Figure 9. Critical Grashof numbers versus R e in the range 8 R e 8 for different values of F.
Fluids 06 00342 g009
Table 1. The values of the Grashof number G r for different number of collocation points N.
Table 1. The values of the Grashof number G r for different number of collocation points N.
N Gr
30977.416134
40977.387524
50977.392292
60977.392292
70977.392292
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zigunovs, M.; Kolyshkin, A.; Iltins, I. Linear Stability of a Steady Convective Flow between Permeable Cylinders. Fluids 2021, 6, 342. https://doi.org/10.3390/fluids6100342

AMA Style

Zigunovs M, Kolyshkin A, Iltins I. Linear Stability of a Steady Convective Flow between Permeable Cylinders. Fluids. 2021; 6(10):342. https://doi.org/10.3390/fluids6100342

Chicago/Turabian Style

Zigunovs, Maksims, Andrei Kolyshkin, and Ilmars Iltins. 2021. "Linear Stability of a Steady Convective Flow between Permeable Cylinders" Fluids 6, no. 10: 342. https://doi.org/10.3390/fluids6100342

APA Style

Zigunovs, M., Kolyshkin, A., & Iltins, I. (2021). Linear Stability of a Steady Convective Flow between Permeable Cylinders. Fluids, 6(10), 342. https://doi.org/10.3390/fluids6100342

Article Metrics

Back to TopTop