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
. 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
ε.
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
rather than an expansion in the Rayleigh, number, R
a≫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]
where
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:
We also assume that the variation of
can be neglected, and this combination is replaced by
, where
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:
The boundary conditions for the horizontal annular cylindrical cell (
Figure 1) are
The system of equations is written in a dimensionless form by scaling length, time, velocity, temperature, and pressure, respectively, by:, , , , and replacing , where and are, respectively, the mass diffusion and thermodiffusion coefficient in a porous medium.
Finally, we obtain the following dimensionless system of equations:
and the associated dimensionless boundary conditions:
where
−cos (
.
The filtration Rayleigh number, Ra, and the separation factor, , are defined by
(/ and . Here, is the kinematic viscosity of the mixture, is the Lewis number, and 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,
is based on the thickness of the annulus,
and not on the radius,
.
is deduced from
Ra by the following relation:
= Ra. For this study, we made the following variable change
, so
verifies
.
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
and for two values of the radius ratio (2
1/8 and 2
1/4). For 2
1/16, the Galerkin method and approximations at orders 3, 7, 15 are in very good agreement for Rayleigh numbers ranging from 50 to 10
4.
We observed that the global Nusselt number is a function of
, 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]
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
at order 3 provides a very good approximation for low values of ε and moderate values of Rayleigh number,
, (for
. 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 << . We assume that we are in the presence of a parallel convective flow. This means that apart from the neighborhood of
, and
, the streamlines are almost arcs of circles parallel to the cavity boundaries. Under these conditions, we show that
| | << |
| outside the neighborhood of
, and
.
When
<< 1 and
, the energy equation reduces to
By taking the boundary conditions (6) into account, we obtain the following analytical expression of the temperature field:
By introducing the stream function, defined by
−∂/r∂θ and
∂/∂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
The conservation equation of mass fraction is given by
Setting
Equation (8) leads to
The circles of radius and connected by the vertical segments = 0 and = forms a streamline associated with φ = 0, and it implies that and .
By taking these two boundary conditions into account in the differential Equation (10), we have the following analytical expression of the solution
From the expression of
we deduce the analytical expression of the two components of the velocity
where the numerator M is given by
Figure 2, shows that the variations of the stream function
φ(r,
with
r for
,
Ra = 100 and
R =2
1/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 .
In
Figure 3,
(r,
is presented as a function of
r for
Ra = 100 and
R = 2
1/16.
In
Figure 4(r,
and the norm of the gradient of the stream function,
are presented as a function of
r for
Ra = 100 and
R = 2
1/16. From this figure, we deduce that
, since:
(r,
is drawn in brown and
, in black.
Figure 4 shows indirectly that
. 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:
with associated boundary conditions:
Although << in the central part of the horizontal annular column, this inequality is not valid for in the vicinity of 0 and since = 0 for 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
and
by their expressions
and
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
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 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 .
Figure 6 illustrates the variation in the mass fraction field obtained numerically for
Ra = 50,
Le = 100,
and
R = 2
1/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 , 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 . From the velocity and temperature fields, we were able to calculate the mass fraction, C(, 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.