Next Article in Journal
Breakthrough Investigation of Advective and Diffusive Transport in a Porous Matrix with a Crack
Next Article in Special Issue
Parameter Variability in Viscous Convection
Previous Article in Journal
On a Casson Fluid Motion: Nonuniform Width Symmetric Channel and Peristaltic Flows
Previous Article in Special Issue
Linear Stability of a Steady Convective Flow between Permeable Cylinders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical and Analytical Studies of Soret-Driven Convection Flow Inside an Annular Horizontal Porous Cavity

by
Abdelkader Mojtabi
1,*,
Khairi Sioud
1,2,
Alain Bergeon
1 and
Marie Catherine Charrier-Mojtabi
1
1
Institut de Mécanique des Fluides de Toulouse (IMFT), Université de Toulouse, CNRS, 31400 Toulouse, France
2
National Institute of Applied Science and Technology, Université de Carthage, La Marsa 2078, Tunisia
*
Author to whom correspondence should be addressed.
Fluids 2021, 6(10), 357; https://doi.org/10.3390/fluids6100357
Submission received: 18 August 2021 / Revised: 30 September 2021 / Accepted: 1 October 2021 / Published: 11 October 2021
(This article belongs to the Special Issue Convection in Fluid and Porous Media)

Abstract

:
This paper studies the species separation of a binary fluid in a porous cavity between two horizontal concentric cylinders, submitted to a temperature gradient. The thickness of the cavity is e = R o R i , where R i and R o are the internal and external radius, respectively. The numerous previous experiments performed in thermogravitational vertical columns (TGCs) showed that in order to obtain a significant separation, the thickness of the cell must be very small, compared with its height. Therefore, in our configuration, we considered e R i . The solution is assumed to be axisymmetric. Under the assumptions of parallel flow and forgotten effect, an analytical solution is obtained using Maple software, and the results are compared with those found numerically using Comsol Multiphysics. In natural convection, our results are in very good agreement with those evaluated with a regular perturbation method in powers of the dimensionless gap width ε = e R i   of order 15, and with the Galerkin method. The species separation calculated for our configuration is very close to the one obtained in a TGC column of height: H = π R i . One of the main interests of the analytical solution presented here is that it can be used as a basic solution for a stability study analysis.

1. Introduction

Natural convection in a horizontal annular cylinder filled with a monoconstituent fluid or with a porous medium saturated by a pure fluid has been widely studied [1,2,3,4]. Subsequently, the authors focused on the linear stability of the convective flow which arises regardless of the temperature difference between the inner and the outer cylinder. The energy stability of this flow was analyzed by Mojtabi and Caltagirone [5]. The approximate analytical solution of the convective flow, for low-temperature differences or low Rayleigh numbers, was obtained from a two- or three-order expansion of the Rayleigh number [1]. In 1987, Rao et al. [6] used a Galerkin method to obtain the convective flow solution. The approximation orders used by these authors did not exceed 20 in the radial and azimuthal directions. In 1988, Himasekhar and Bau [7] extended the regular perturbation method in powers of the filtration Rayleigh number up to R a 60 . In 1991, Charrier-Mojtabi et al. [8] numerically compared Fourier–Galerkin and collocation–Chebyshev methods and showed that the latter gave a more accurate description of the boundary layers near the inner and outer cylinders. In 1992, Mojtabi et al. [9] calculated the convective flow using a regular perturbation method in powers of the dimensionless gap ε.
ε = R o R i R i = R 1 ,   where   R = R o R i .
A symbolic algebra code allowed the authors to obtain the expanded solution at any order of ε theoretically. In their article, the calculations were carried out up to order 7 because the computer used to run the symbolic software had limited available memory. The results calculated, at O(7), for several radius ratios were compared with those obtained by the Galerkin method. From this comparison, it was found that the agreement between the results obtained by the two methods decreased when the ratio of the radii, R, or the Rayleigh number increased. We have recently extended the development to order O(15). This will be the subject of the first part of this work.
Although the species separation has been extensively studied in vertical annular cylindrical columns [10] or inclined rectangular porous cells [11,12,13,14,15,16], the same cannot be said for the porous horizontal annular column, and, to our knowledge, the only reference is the paper of Abahri et al. [17], in which the authors obtained an approximate solution of the problem using a second-order expansion in Ra and a numerical solution using the Comsol software. In vertical thermogravitational columns, the extent of species separation is proportional to the height of the cavity. The distance between the two opposite isothermal, vertical or inclined surfaces leading to the optimum separation is less than 1 mm [11,12]. It follows that the annular cylindrical column should be of small thickness and large interior radius in order for the streamlines to be long enough to allow optimal species separation. For the study of species separation, it is, therefore, more relevant to write the solution of the problem in the form of an expansion in ε = ( R o R i ) / R i 1   rather than an expansion in the Rayleigh, number, Ra≫1. The results in [9] show that the flow in the horizontal annular cylindrical column is almost conductive. In the continuation of this work, we seek to obtain a solution for the thermogravitational problem when the thermal field is purely conductive.

2. Mathematical Formulation

Thermodiffusion induces a mass fraction flux in binary fluid mixtures subjected to a temperature gradient. In addition to the usual expression for the mass flux density given by Fick’s law, part of the expression due to the temperature gradient (namely, the Soret effect) is used so that [11,12]
J m = ρ   D * C ρ   C ( 1 C ) D T * T ,
where J m is the vector density of mass flux. Darcy’s law is assumed valid when binary fluid and solid phases are in local thermal equilibrium. It is also assumed that viscous dissipation, compressibility, and Dufour effects are negligible. The Boussinesq approximation is used as follows:
ρ = ρ 0 [ 1 β T ( T T 0 ) β C ( C C 0 ) ]
We also assume that the variation of C ( 1 C ) can be neglected, and this combination is replaced by   C 0 ( 1 C 0 ) , where C 0   is the initial value of the mass fraction. Under these assumptions, the equations of continuity, momentum, energy, and species conservation are written in their dimensional form as follows:
{ . V = 0         V = K μ   ( P ρ 0 [ 1 β T ( T T 0 ) β C ( C C 0 ) ] g )         ( ρ c ) * T t + ( ρ c ) f V . T = k *   2 T                     ϵ * C t   +   V . C = D *   2 C + D T * C 0 ( 1 C 0 ) 2 T   ,
The boundary conditions for the horizontal annular cylindrical cell (Figure 1) are
{ r = R i ,         T = T h ,       V . n = 0 ,         J m . n = 0 r = R o ,         T = T c ,       V . n = 0 ,         J m . n = 0 θ = 0 ,   π ,     V θ = T θ = C θ = 0                              
The system of equations is written in a dimensionless form by scaling length, time, velocity, temperature, and pressure, respectively, by:   R i , ( ρ c ) *   R i 2 k * , k * R i ( ρ c ) f   , ( T h T c ) , k * μ K ( ρ c ) f and replacing C C 0   b y   Δ C = Δ T C 0 ( 1 C 0 ) D T * / D * , where D * and D T *   are, respectively, the mass diffusion and thermodiffusion coefficient in a porous medium.
Finally, we obtain the following dimensionless system of equations:
{ . V = 0   V + P =   R a (   T + Ψ   C )   z                           T t + V . T = 2 T                                                               ϵ C t + V . C = 1 L e ( 2 C 2 T )
and the associated dimensionless boundary conditions:
{ r = 1 ,     T = 1 ,   V r . e r = 0 ,     ( C T ) . e r = 0 r = R ,     T = 0 ,     V r . e r = 0 ,     ( C T ) . e r = 0   θ = 0 ,     π ,     V θ = T θ = C θ = 0                                                        
where z = −cos ( θ ) e r + sin ( θ ) e θ .
The filtration Rayleigh number, Ra, and the separation factor, Ψ , are defined by
R a = g β T K R i (   T h T c ) ( ρ c ) f / ν k * and Ψ = β C D T * C 0 ( 1 C 0 ) / β T D * . Here, ν = μ / ρ is the kinematic viscosity of the mixture, L e = k * ( ρ c ) f D *   is the Lewis number, and ϵ = ϵ * ( ρ c ) f   ( ρ c ) *   is the modified porosity.

3. Some Results in Natural Convection

In the present work, the development of the convective solution in the porous horizontal annular space was obtained at order 15 and extended the results previously obtained to order 7 in [9]. The thermal Rayleigh number, R a t h is based on the thickness of the annulus,   ( R o R i ) , and not on the radius, R i . R a t h is deduced from Ra by the following relation: R a t h = Ra ( R 1 ) = ε R a . For this study, we made the following variable change x = ( r 1 ) / ε , so x verifies x [ 0 , 1 ] .
Table 1 presents the average Nusselt number, Nug, obtained from an expansion function of ε at orders 3, 7, and 15 and compares it with those calculated by a Galerkin method for various values of R a t h and for two values of the radius ratio (21/8 and 21/4). For 21/16, the Galerkin method and approximations at orders 3, 7, 15 are in very good agreement for Rayleigh numbers ranging from 50 to 104.
We observed that the global Nusselt number is a function of R a t h 2 , which is in agreement with the results of Himasekhar and Bau [7]. It is also interesting to note that the order 4 approximation obtained by Mojtabi [9]
N u g = 1 + 17 R a t h 2 ( ε 2 ε 3 ) / 40 , 320 + O ( ε 4 ) ,
already provides a good approximation of the global Nusselt number. This approximation can be applied for a wide range of variations of the Rayleigh number. Table 1 shows that the approximation in ε = R o R i R i = R 1   at order 3 provides a very good approximation for low values of ε and moderate values of Rayleigh number,   R a t h , (for R = 2 1 / 32 ,   2 1 / 16 ,   2 1 / 6 ,   2 1 / 4 ,   ε 0.0219 ,   0.0443 ,   0.090 ,   0.189 ) . The Nusselt number, Nug remains close to 1.

4. Soret-Driven Convection: Species Separation

From Table 1, we deduce that the temperature field depends only on the radial coordinate r as, for moderate Rayleigh numbers and for low values of R, the Nusselt number, Nug remains close to 1, which corresponds to a pure conduction regime in the annular space. This result is obtained when the geometrical characteristics of the cavity correspond to the operating conditions, ensuring a good species separation: e << R i . We assume that we are in the presence of a parallel convective flow. This means that apart from the neighborhood of θ = 0 , and θ = π , the streamlines are almost arcs of circles parallel to the cavity boundaries. Under these conditions, we show that | V r | << | V θ | outside the neighborhood of θ = 0 , and θ = π .
When V r  << 1 and T θ = 0 , the energy equation reduces to
2 T = 0
By taking the boundary conditions (6) into account, we obtain the following analytical expression of the temperature field:
T = 1 ln ( r ) / ln ( R )
By introducing the stream function, φ defined by
V r = −∂ φ /r∂θ and V θ = φ /∂r, the continuity equation is automatically verified. The new formulation of the problem in dimensionless form, with the hypothesis of the forgotten effect, previously used by Furry Jones and Onsager (FJO, theory) [18] is given by
2 φ = R a   sin ( θ ) T r
The conservation equation of mass fraction is given by
2 C L e = V r C r + V θ C r θ
Setting φ = R a   sin ( θ   ) f ( r ) Equation (8) leads to
( r f r ) r f r r T r = 0
The circles of radius R i   and R o connected by the vertical segments θ = 0 and θ = π forms a streamline associated with φ = 0, and it implies that f ( r = 1 ) = 0 and   f ( r = R ) = 0 .
By taking these two boundary conditions into account in the differential Equation (10), we have the following analytical expression of the solution
f ( r ) = R 2 ( r 2 1 ) ln ( R ) r 2 ( R 2 1 ) ln ( r ) 2 r ( R 2 1 ) ln ( R )
From the expression of f ( r ) , we deduce the analytical expression of the two components of the velocity
{ V r a = R a   f ( r ) cos ( θ ) / r V θ a = M 2 r 2 ( R 2 1 ) ln ( R ) ,
where the numerator M is given by
M = R a   sin ( θ ) ( R 2 ( r 2 + 1 ) ln ( R ) r 2 ( R 2 1 ) ( ln ( r ) + 1 ) )
Figure 2, shows that the variations of the stream function φ(r,   θ )   with r for θ = π / 2 , Ra = 100 and R =21/16 obtained analytically and from direct numerical simulations, using Comsol Multiphysics software, are in very good agreement. The solid brown line is taken from the analytical results, and the dots in black represent the results of the direct numerical simulation.
The analytical results obtained for low values of the radius ratio, R, and moderate values of the thermal Rayleigh number, Ra, show that V r V θ .
In Figure 3, V θ (r,   θ = π / 2 )   is presented as a function of r for Ra = 100 and R = 21/16.
In Figure 4 ,   V θ (r,   θ = π 2 )   and the norm of the gradient of the stream function,   | φ | are presented as a function of r for Ra = 100 and R = 21/16. From this figure, we deduce that   V r V θ , since: | φ |     V θ ( 1 + ( V r V θ ) 2 ) 1 / 2   V θ (r,   θ = π 2 ) is drawn in brown and | φ | , in black. Figure 4 shows indirectly that   | φ |     V θ . The dots in black are obtained from the direct numerical simulation.
The conservation equation of mass fraction (9) reduces to a linear equation of the second order as follows:
( r C r ) r + 2 C r 2 θ = L e ( r V r C r + V θ C θ ) ,    
with associated boundary conditions:
{ r = 1 ,         C r + 1 r l n ( R ) = 0 r = R ,         C r + 1 r l n ( R ) = 0   θ = 0 ,       π ,     C θ = 0  
Although V r << V θ in the central part of the horizontal annular column, this inequality is not valid for θ in the vicinity of 0 and   π since V θ = 0 for θ = 0 and for θ = π .
Therefore, it was not possible to obtain an explicit analytical relation for the mass fraction C solution of Equation (14). On the other hand, we replaced V r and V θ by their expressions V r a and V θ a given by the system of Equation (12) to obtain an equation where the only unknown is C, which we solved numerically. Then, comparing the values of the mass fraction C obtained from direct numerical simulation (resolution of the system with four conservation Equation (5), without simplifying assumption) and of the numerical resolution of Equation (16) alone, with only C as unknown, associated with boundary conditions (15), we have
( r C r ) r + 2 C r 2 θ = L e ( r V r a C r + V θ a C θ )
The results of direct numerical simulation and of the solution obtained from Equation (16) alone are compared in Figure 5.
The evolution of C as a function of θ for r = R shows that the two approaches are perfectly concordant and that, in the central part of the horizontal column, the evolution of C as a function of θ is linear for any value of r [ 1 , R ] .
Figure 6 illustrates the variation in the mass fraction field obtained numerically for Ra = 50, Le = 100, Ψ = 0.2 and R = 21/4. The colored scale represents the intensity of the mass fraction of the heaviest component. The lines represent the associated streamlines.
Under the effect of thermodiffusion, the convective regime within the annular column leads to a low value of C (blue) at the top of the column and a higher value of C (red) at its lower part. The difference between these two extreme values of C indicates the degree of species separation of the mixture.

5. Conclusions

Recent developments in the field of symbolic computation have made it possible to compute very accurate analytical solutions to non-linear problems. The analytical solution obtained previously to order 7, using a regular asymptotic development in powers of in ε = R o R i R i 1 , was extended here to order 15. The solution obtained for the problem of natural convection in a horizontal porous annular cylinder for a wide range of Rayleigh numbers could be considered as a reference solution for the validation of future numerical 2D codes.
However, for the study of the species separation in a binary fluid, this solution sought in the development of ε is not convenient for carrying out a stability analysis or obtaining analytically, in this case, the optimum of species separation. In the case of the parallel flow approximation and forgotten effect, we showed good agreement between the direct numerical simulation and the analytical expression of the convective flow. We established that the radial velocity component is negligible, compared with the tangential component, except in the vicinity of θ = 0 , π . From the velocity and temperature fields, we were able to calculate the mass fraction, C( r ,   θ ) , using a numerical resolution of a single scalar equation with partial derivatives. This solution is in good agreement with the results of direct numerical simulations.
We obtained for the first time an analytical solution of Soret-driven convection flow in an annular, thin, horizontal porous layer saturated with a binary fluid.
Until now, the species separation was carried out in parallelipipedic or annular vertical columns. In this study, we showed that the species separation of a binary mixture can also be obtained in a horizontal porous annular column.

Author Contributions

Conceptualization, A.M., methodology, software, validation, writing, A.M., K.S., A.B., M.C.C.-M., All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Exclude this statement.

Acknowledgments

This work was supported by CNES, the French National Space Agency.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

CMass fraction
C 0 Initial mass fraction
D Mass   diffusion   coefficient ,   [ m 2 s 1 ]
D * Mass   diffusion   coefficient   in   a   porous   medium   [ m 2 s 1 ]
D T * Thermodiffusion   coefficient   in   porous   medium   [ m 2 s 1 K 1 ] .
e Thickness   of   the   two   cells   [ m ]
g Gravitational   acceleration   [ ms 2 ]
H Height   of   the   cavity   [ m ]
K Permeability   of   porous   medium   [ m 2 ]
k*Effective thermal conductivity [W/m K]
LeLewis number
P′Pressure [Pa]
R Radius   ratio ,   R = R o R i
Ra Filtration   Rayleigh   number ,   R a = g β T K R i (   T h T c ) ( ρ c ) f / ν k *
R a t h R a t h = R a   ( R 1 ) = ε R a  
R i Inner   sphere   radius   [ m ]
R o Outer   sphere   radius   [ m ]
t′Time, [s]
T′Temperature [K]
T c Cold temperature [K]
T h Hot temperature [K]
T 0 Initial temperature [K]
V r Radial velocity
V θ Tangential velocity
Greek symbols
β T Thermal   expansion   coefficient ,   [ K 1 ]
β C Solutal expansion coefficient
εDimensionless gap, ε = R 1
ε*Porosity of the porous medium
ϵ Modified porosity
μDynamic viscosity of the mixture, [m2s−1]
ρDensity of the binary fluid [kg/m3]

References

  1. Caltagirone, J.P. Thermoconvective instabilities in a porous medium bounded by two concentric horizontal cylinders. J. Fluid Mech. 1976, 76, 337–362. [Google Scholar] [CrossRef]
  2. Kuehn, T.H.; Goldstein, R.J. An experimental and theoretical study of natural convection in the annulus between horizontal concentric cylinders. J. Fluid Mech. 1976, 74, 695–719. [Google Scholar] [CrossRef]
  3. Kuehn, T.H.; Goldstein, R.J. An experimental study of natural convection heat transfer in concentric and eccentric horizontal cylindrical annuli. J. Heat Transf. 1978, 100, 635–640. [Google Scholar] [CrossRef]
  4. Burns, P.J.; Yien, C.L. Natural convection in porous medium bounded by concentric spheres and horizontal cylinders. Int. J. Heat Mass Transf. 1979, 22, 929–939. [Google Scholar] [CrossRef]
  5. Mojtabi, A.; Caltagirone, J.P. Energy stability of natural convective flow in horizontal annular space. Phys. Fluid 1979, 22, 1208–1209. [Google Scholar] [CrossRef] [Green Version]
  6. Rao, Y.F.; Fukuda, K.; Hasegawa, S. Study and transient analysis of Natural Convection in Porous Annular with Galerkin method. J. Heat Transf. ASME 1987, 109, 919–927. [Google Scholar] [CrossRef]
  7. Himasekhar, K.; Bau, H.H. Two-dimensional bifurcation phenomena in thermal convection in horizontal concentric annuli containing saturated porous media. J. Fluid Mech. 1988, 187, 267–300. [Google Scholar] [CrossRef]
  8. Charrier-Mojtabi, M.C.; Mojtabi, A.; Azaiez, M.; Labrosse, G. Numerical and experimental study of multicellular free convection flows in an annular porous layer. Int. J. Heat Mass Transf. 1991, 34, 3067–3074. [Google Scholar] [CrossRef]
  9. Mojtabi, A.; Charrier-Mojtabi, M.C. Analytical solution of steady natural convection in an annular porous medium evaluated with a symbolic algebra code. J. Heat Transf. 1992, 114, 1065–1067. [Google Scholar] [CrossRef]
  10. Marcoux, M.C.; Charrier-Mojtabi, M.C.; Azaiez, M. Double diffusive convection in annular vertical porous layer. Int. J. Heat Mass Transf. 1999, 42, 2313–2325. [Google Scholar] [CrossRef]
  11. Seta, B.; Lapeira, E.; Dubert, D.; Gavalda, F.; Bou-Ali, M.M.; Ruiz, X. Separation under thermogravitational effects in binary mixtures. Eur. Phys. J. E 2019, 42, 58. [Google Scholar] [CrossRef] [PubMed]
  12. Platten, J.K.; Bou-Ali, M.M.; Dutieux, J.F. Enhanced Molecular Separation in Inclined Thermogravitational Columns. J. Phys. Chem. 2003, 107, 11763–11767. [Google Scholar] [CrossRef]
  13. Elhajjar, B.; Mojtabi, A.; Costesèque, P.; Charrier-Mojtabi, M.C. Separation in an inclined porous thermogravitational cell. Int. J. Heat Mass Transf. 2010, 53, 4844–4851. [Google Scholar] [CrossRef] [Green Version]
  14. Charrier-Mojtabi, M.C.; Elhajjar, B.; Mojtabi, A. Analytical and numerical stability analysis of Soret-driven convection in a horizontal porous layer. Phys. Fluids 2007, 19, 124104. [Google Scholar] [CrossRef] [Green Version]
  15. Elhajjar, B.; Mojtabi, A.; Charrier-Mojtabi, M.C. Separation of a binary fluid mixture in a porous horizontal cavity. Phys. Rev. E 2008, 77, 2. [Google Scholar] [CrossRef] [PubMed]
  16. Yacine, L.; Mojtabi, A.; Bennacer, R.; Khouzam, A. Soret-driven convection and separation of binary mixtures in a horizontal porous cavity submitted to cross heat fluxes. Int. J. Therm. Sci. 2016, 104, 29–38. [Google Scholar] [CrossRef] [Green Version]
  17. Abahri, O.; Sadaoui, D.; Mansouri, K.; Mojtabi, A.; Charrier-Mojtabi, M.C. Thermogravitational separation in horizontal annular porous cell. Mech. Ind. 2017, 18, 106. [Google Scholar] [CrossRef]
  18. Furry, W.H.; Jones, R.C. On the theory of isotope separation by thermal diffusion. Phys. Rev. 1939, 55, 1083–1095. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of horizontal cylindrical TGC.
Figure 1. Schematic diagram of horizontal cylindrical TGC.
Fluids 06 00357 g001
Figure 2. Stream function φ(r,   θ = π / 2 )   for Ra = 100 and R = 21/16 obtained analytically (continuous line) and using numerical direct simulation (black dots).
Figure 2. Stream function φ(r,   θ = π / 2 )   for Ra = 100 and R = 21/16 obtained analytically (continuous line) and using numerical direct simulation (black dots).
Fluids 06 00357 g002
Figure 3. V θ versus r for Ra = 100 and R = 21/16 obtained analytically (continuous line) and using numerical direct simulations (black dots).
Figure 3. V θ versus r for Ra = 100 and R = 21/16 obtained analytically (continuous line) and using numerical direct simulations (black dots).
Fluids 06 00357 g003
Figure 4. V θ and | φ | versus r for θ = π 2 ,   Ra = 100 and R = 21/16 obtained analytically (continuous line) and using numerical simulations (black dots).
Figure 4. V θ and | φ | versus r for θ = π 2 ,   Ra = 100 and R = 21/16 obtained analytically (continuous line) and using numerical simulations (black dots).
Fluids 06 00357 g004
Figure 5. C ( r = R ,   θ ) versus   θ for Ra = 12.6, Le = 100, Ψ = 0.21 and R = 21/16 obtained analytically (continuous line) and using numerical direct simulation (black dots).
Figure 5. C ( r = R ,   θ ) versus   θ for Ra = 12.6, Le = 100, Ψ = 0.21 and R = 21/16 obtained analytically (continuous line) and using numerical direct simulation (black dots).
Fluids 06 00357 g005
Figure 6. Variation in the mass fraction for R a = 50 ,   L e = 100 , Ψ = 0.2 and R = 2 1 / 4   blue (red) corresponds to low (high) mass fraction values.
Figure 6. Variation in the mass fraction for R a = 50 ,   L e = 100 , Ψ = 0.2 and R = 2 1 / 4   blue (red) corresponds to low (high) mass fraction values.
Fluids 06 00357 g006
Table 1. Average Nusselt number (Nug − 1) for different values of R a t h and R.
Table 1. Average Nusselt number (Nug − 1) for different values of R a t h and R.
R a t h GalerkinO(15)O(7)O(3)
R =21/8
103252 × 10−4252 × 10−4252 × 10−4283 × 10−4
2 × 103946 × 10−4944 × 10−4951 × 10−41132 × 10−4
4 × 103309 × 10−3352 × 10−3373 × 10−3453 × 10−3
R =21/4
200177 × 10−4177 × 10−4177 × 10−4216 × 10−4
500102 × 10−3102 × 10−3102 × 10−3135 × 10−3
750208 × 10−3215 × 10−3232 × 10−3304 × 10−3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mojtabi, A.; Sioud, K.; Bergeon, A.; Charrier-Mojtabi, M.C. Numerical and Analytical Studies of Soret-Driven Convection Flow Inside an Annular Horizontal Porous Cavity. Fluids 2021, 6, 357. https://doi.org/10.3390/fluids6100357

AMA Style

Mojtabi A, Sioud K, Bergeon A, Charrier-Mojtabi MC. Numerical and Analytical Studies of Soret-Driven Convection Flow Inside an Annular Horizontal Porous Cavity. Fluids. 2021; 6(10):357. https://doi.org/10.3390/fluids6100357

Chicago/Turabian Style

Mojtabi, Abdelkader, Khairi Sioud, Alain Bergeon, and Marie Catherine Charrier-Mojtabi. 2021. "Numerical and Analytical Studies of Soret-Driven Convection Flow Inside an Annular Horizontal Porous Cavity" Fluids 6, no. 10: 357. https://doi.org/10.3390/fluids6100357

APA Style

Mojtabi, A., Sioud, K., Bergeon, A., & Charrier-Mojtabi, M. C. (2021). Numerical and Analytical Studies of Soret-Driven Convection Flow Inside an Annular Horizontal Porous Cavity. Fluids, 6(10), 357. https://doi.org/10.3390/fluids6100357

Article Metrics

Back to TopTop