# On the Nondimensionalization Process in Complex Problems: Application to Natural Convection in Anisotropic Porous Media

Authors: I. Alhama; M. Cánovas; F. Alhama
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The nondimensionalization of the equations governing a given problem is a direct, relatively easy, and low-cost way to provide interesting information, the dimensionless groups that rule the problem and define its final patterns of solution. In complex problems, however, this technique frequently does not provide the precise and complete set of monomials we are looking for. The use of discrimination in the process of nondimensionalization, the detailed application of which is explained in this paper, always leads to a minimum set of parameters, which, separately, determine the solution of the problems. In addition, the physical meaning and order of magnitude of these discriminated monomials are also provided by the discrimination. The technique is applied to the coupled problem of natural convection between horizontal plates heated from below, containing an anisotropic porous medium.


---

## Body

## 1. Introduction

The search for the minimum group of parameters that rule the solution to a given problem, especially if it is complex and can only be solved numerically, is perhaps the main goal of experimental, numerical, and theoretical research and great effort has been put into the same. Apart from experimental work, there are basically two techniques that are used: the application of dimensional analysis and the nondimensionalization of the governing equations. The cost of these nonexperimental techniques is very low in light of the advantages they provide.The first procedure [1–3], which is followed in many textbooks for academic purposes, starts from the correct set of relevant quantities whose potential flows balance each other in the conservation laws of the problem, as well as of physical characteristics that related the relevant quantities in the conservation laws. By choosing a suitable dimensional basis, formed by the fundamental quantities of the theory in which the problem is involved, and applying the pi theorem, it is possible to derive the final dimensionless groups. However, as mentioned above, this technique is scarcely used and few works that use it can be found in the scientific literature [4–9]. Two fundamental reasons justify this absence and its progressive abandonment. Firstly, a deep understanding of the theory involved in the physical process is required for correct formulation of the complete and unique list of variables; the researcher frequently does not include in this list hidden quantities that are not mentioned in the statement of the problem but which, however, take part in the balances or conservation laws; this occurs, for example, in 2D scenarios in which one of the lengths extends beyond the limits of where the phenomenon is mainly located. Secondly, even when the list of relevant variables is correctly formulated, the classical application of pi theorem does not always provide the independent dimensionless groups unless discrimination is taken into account in the nondimensionalization procedure. The concept of discrimination as an extension of dimensional analysis is introduced by Madrid and Alhama [10] in the field of convective heat transfer and used by these and other authors for the derivation of the correct dimensionless groups [11–13].The use of the more formal nondimensionalization technique, however, is extensively used both in advanced textbooks and in the scientific literature in general [14–17]. In its classical form, this technique does not distinguish between scalar and vectorial quantities but chooses any length, independently of its spatial direction, to make all the independent spatial variables dimensionless. The same occurs with velocities and other quantities of a vectorial character. This nondiscriminated (classical) procedure frequently leads to results that can be improved. Many of the dimensionless numbers, such as Reynolds, Rayleigh, Nusselt, and Grashof (or their reverse), have very high or very low values compared with unity and lack any sense in terms of quantity balances that counteract the physical process; indeed, this would make an interesting research subject. In contrast, the use of discrimination clearly improves the solution in fundamental aspects: the number of resulting groups is less, their order of magnitude is unity, and (as a consequence) they can be immediately interpreted physically as balances of quantities in the domains where they are applied. Madrid and Alhama [17] and Madrid and Alhama [18] applied the discriminated nondimensionalization to solve the problems of laminar forced convection along a vertical plate and natural convection along an isothermal vertical plate, providing very precise results.In this paper, mathematically orientated from the point of view of equations management but without losing the physical understanding of the problem, discriminated nondimensionalization is applied to a coupled, nonlineal, very complex problem involving many physical parameters, in order to demonstrate not only the formal rules of application, but also some complementary and interesting aspects derived from this technique, such as the potential existence of hidden quantities necessary for the definition of dimensionless dependent or independent variables. Subsequently, when dimensionless groups have been derived, the meaning and even the order of magnitude of these quantities can be elucidated. The problem, studied by Holzbecher and Yusa [19] and Holzbecher [20], refers to the natural heat convection between large horizontal plates, heated from below, containing an anisotropic porous medium. In this kind of problem, the flow is organized in close structures regularly distributed, named cells (or Bénard’s cells), along the horizontal coordinate.After the discriminated dimensionless groups have been obtained, many scenarios in which these groups retain the same value but in which the physical characteristics or geometrical parameters are changed are solved numerically in order to demonstrate the validity of the solutions. By making use of the physical meaning of the resulting groups, the order of magnitude of the principal unknown of the problem can be obtained.


## 2. Governing Equations and Boundary Conditions

The scenario studied is shown in Figure1. It refers to a fluid-flow and heat-transport process—a natural convection problem—occurring between horizontal plates subjected to a thermal gradient and impermeable to fluid. The scenario is 2D rectangular, and the porous medium assumes anisotropy for the hydraulic and thermal conductivities.Figure 1
Natural convection between horizontal plates. Scheme of the problem.


### 2.1. Dimensional Form

In the steady state, the set of governing equations is formed by (i) the expressions of the horizontal and vertical components of the flow velocity,u and v, respectively, derived from Darcy’s law, (ii) the mass conservation equation, and (iii) energy conservation equation. For anisotropic 2D media, using rectangular coordinates, these are
(1)u=-(Kxμ)(∂P∂x),v=-(Kyμ){(∂P∂y)+ρg},(2)(∂u∂x)+(∂v∂y)=0,(3)(ρfcp,f)u(∂T∂x)+(ρfcp,f)v(∂T∂y)=km,x(∂2T∂x2)+km,y(∂2T∂y2).
In these equations K is the hydraulic conductivity (m/s), P the pressure (N/m2), ρ the fluid density (kg/m3), g the gravitational acceleration (m/s2), μ the fluid viscosity (kg m−1 s−1), cp,f the specific heat (J/Kg°C), k the thermal conductivity (W/m°C), and T the temperature (°C); finally, x(m) and y(m) are the independent variables. Subscript m means average value, f and s refer to fluid and solid matrix, respectively, and x and y denote the anisotropic character of the parameter. Average thermal conductivity is defined in terms of the porosity, ϕ (dimensionless), and the fluid and solid matrix conductivities (kf and ks, resp.), km=[ϕkf+(1-ϕ)ks].Eliminating the pressure in the Darcy equation and making use of the Boussinesq approximation (a hypothesis widely used in most processes in this field, Holzbecher [20]),
(4)ρ=ρ0-ρ0βΔT=ρ0(1-βΔT),orΔρ=-ρ0β(ΔT),
where β is the thermal expansion coefficient (°C−1), subscript 0 refers to a reference temperature, and the momentum equations are combined in the form
(5)(KyKx)(∂u∂y)-(∂v∂x)=-(ρKygβμ)(∂T∂x).
This reduces the mathematical models to (2), (3), and (5). As regards the pressure variable, (1) would be used for deriving its changes from the characteristic velocities obtained in this work. To complete the problem, the following boundary conditions are assumed: temperature at the top and bottom boundaries is isothermal (Tb>Tt), with ΔT=Tb-Tt, while both left and right sides are adiabatic; in addition, all sides are impermeable to water flow:
(6)T(x,y=0)=Tb,T(x,y=H)=Tt,(∂T∂x)left=(∂T∂x)right=0,q(x,y=0)=q(x,y=H)q(x=0,y)q(x=L,y)=0,
where q is the fluid flow (m/s) and L and H define the geometry of the rectangular domain. The aspect ratio of the problem L/H is sufficiently large, L≫H.


### 2.2. Discriminated Dimensionless Form

#### 2.2.1. The Reference Quantities

Discrimination forces the choice of different characteristic values, according to their spatial direction, to make the variables of vectorial character, velocity (u and v), and coordinates (x and y) dimensionless. On the one hand, as regards the characteristic lengths, although two discriminated lengths appear in the statement of the problem, one of them (L) is not relevant since its value is large enough to neither influence the problem nor rule the solution. It is easy to deduce that the steady state patterns cannot be dependent on L when L is large enough; that is, patterns may cover the totality of the domain or a part of it, but its aspect remains invariable as L increases. As a consequence, L cannot be chosen as a horizontal reference quantity, and, since no more references exist, we assume a “hidden” length (lx*) for this propose, an unknown whose physical meaning and value will be investigated later, after derivation of the dimensionless groups. In short, the reference quantities to make x and y dimensionless are (lx*) and H, respectively. On the other hand, as regards velocities, the vertical reference can be derived directly from the first addend of the left part of (1) making Δρ=Δρmax⁡, since the set of parameters “(Ky/μ)(Δρmax⁡)g” provides a characteristic value related to the vertical coordinate. Finally, we assume a new hidden quantity, ux*, as the horizontal reference velocity since this cannot be established directly or indirectly by the set of parameters contained in the statement. The length and velocity references are as follows:for thex coordinate: lx*,for they coordinate: ly*=H,for the vertical velocity,v: v*=-(Ky/μ)(Δρ)g=-(Ky/μ)ρ0gβ(ΔT),for the horizontal velocity,u: u*,for the temperature,T: ΔT=Tb-Tt.


#### 2.2.2. The Dimensionless Variables

Based on the above references, the following dimensionless variables can be defined:(7)x′=xlx*,(8)y′=yH,(9)u′=uu*,(10)v′=vv*=vμ[ρ0gβ(ΔT)K],(11)T′=TΔT.
We assume that the dimensionless variables x′ and u′ range in the interval [0,1], even though the references used in their definition are unknowns; this is a hypothesis that will be tested later. Furthermore, it is clear that y′ and v′ also range in the same interval. This property of the dimensionless variables will be used later in the dimensionless governing equation to assign a suitable order of magnitude to the discriminated groups we are looking for and to give a physical meaning and an order of magnitude to the “hidden” references.


#### 2.2.3. The Dimensionless Equations

Firstly, we study the mass conservation equation because of its simplicity. Substituting (7)–(10) into (2) yields
(12)(u*∂u′lx*∂x′)+(v*∂v′H∂y′)=0
or
(13)(u*∂u′lx*∂x′)+([ρ0gβ(ΔT)Ky]∂v′μH∂y′)=0.
Reorganizing this equation by grouping all the parameters in the same addend, we can write
(14)(μHu*[ρ0gβ(ΔT)Ky]lx*)(∂u′∂x′)+(∂v′∂y′)=0.
Note that the expressions within the parentheses are dimensionless.Secondly, we will work with the momentum equation (5). Again, using the dimensionless variables, this equation can be written in the form
(15)(KyKx)(u*∂u′H∂y′)-(v*∂v′lx*∂x′)=-(ρKygβμ)((ΔT)∂T′lx*∂x′)
or, making use of (12) and simplifying in the dimensionless form of the momentum equation,
(16)(Kylx*2KxH2)(∂v′∂y′∂x′∂y′)+(∂v′∂x′)=(∂T′∂x′).
Finally, the nondimensionalization of the energy equation (3), which is carried out in a similar way, leads to
(17)(ρfcp,f)(u*lx*)(u′∂T′∂x′)+(ρfcp,f)(v*H)(v′∂T′∂y′)=km,xlx*2(∂2T′∂x′2)+km,yH2(∂2T′∂y′2).


## 2.1. Dimensional Form

In the steady state, the set of governing equations is formed by (i) the expressions of the horizontal and vertical components of the flow velocity,u and v, respectively, derived from Darcy’s law, (ii) the mass conservation equation, and (iii) energy conservation equation. For anisotropic 2D media, using rectangular coordinates, these are
(1)u=-(Kxμ)(∂P∂x),v=-(Kyμ){(∂P∂y)+ρg},(2)(∂u∂x)+(∂v∂y)=0,(3)(ρfcp,f)u(∂T∂x)+(ρfcp,f)v(∂T∂y)=km,x(∂2T∂x2)+km,y(∂2T∂y2).
In these equations K is the hydraulic conductivity (m/s), P the pressure (N/m2), ρ the fluid density (kg/m3), g the gravitational acceleration (m/s2), μ the fluid viscosity (kg m−1 s−1), cp,f the specific heat (J/Kg°C), k the thermal conductivity (W/m°C), and T the temperature (°C); finally, x(m) and y(m) are the independent variables. Subscript m means average value, f and s refer to fluid and solid matrix, respectively, and x and y denote the anisotropic character of the parameter. Average thermal conductivity is defined in terms of the porosity, ϕ (dimensionless), and the fluid and solid matrix conductivities (kf and ks, resp.), km=[ϕkf+(1-ϕ)ks].Eliminating the pressure in the Darcy equation and making use of the Boussinesq approximation (a hypothesis widely used in most processes in this field, Holzbecher [20]),
(4)ρ=ρ0-ρ0βΔT=ρ0(1-βΔT),orΔρ=-ρ0β(ΔT),
where β is the thermal expansion coefficient (°C−1), subscript 0 refers to a reference temperature, and the momentum equations are combined in the form
(5)(KyKx)(∂u∂y)-(∂v∂x)=-(ρKygβμ)(∂T∂x).
This reduces the mathematical models to (2), (3), and (5). As regards the pressure variable, (1) would be used for deriving its changes from the characteristic velocities obtained in this work. To complete the problem, the following boundary conditions are assumed: temperature at the top and bottom boundaries is isothermal (Tb>Tt), with ΔT=Tb-Tt, while both left and right sides are adiabatic; in addition, all sides are impermeable to water flow:
(6)T(x,y=0)=Tb,T(x,y=H)=Tt,(∂T∂x)left=(∂T∂x)right=0,q(x,y=0)=q(x,y=H)q(x=0,y)q(x=L,y)=0,
where q is the fluid flow (m/s) and L and H define the geometry of the rectangular domain. The aspect ratio of the problem L/H is sufficiently large, L≫H.


## 2.2. Discriminated Dimensionless Form

### 2.2.1. The Reference Quantities

Discrimination forces the choice of different characteristic values, according to their spatial direction, to make the variables of vectorial character, velocity (u and v), and coordinates (x and y) dimensionless. On the one hand, as regards the characteristic lengths, although two discriminated lengths appear in the statement of the problem, one of them (L) is not relevant since its value is large enough to neither influence the problem nor rule the solution. It is easy to deduce that the steady state patterns cannot be dependent on L when L is large enough; that is, patterns may cover the totality of the domain or a part of it, but its aspect remains invariable as L increases. As a consequence, L cannot be chosen as a horizontal reference quantity, and, since no more references exist, we assume a “hidden” length (lx*) for this propose, an unknown whose physical meaning and value will be investigated later, after derivation of the dimensionless groups. In short, the reference quantities to make x and y dimensionless are (lx*) and H, respectively. On the other hand, as regards velocities, the vertical reference can be derived directly from the first addend of the left part of (1) making Δρ=Δρmax⁡, since the set of parameters “(Ky/μ)(Δρmax⁡)g” provides a characteristic value related to the vertical coordinate. Finally, we assume a new hidden quantity, ux*, as the horizontal reference velocity since this cannot be established directly or indirectly by the set of parameters contained in the statement. The length and velocity references are as follows:for thex coordinate: lx*,for they coordinate: ly*=H,for the vertical velocity,v: v*=-(Ky/μ)(Δρ)g=-(Ky/μ)ρ0gβ(ΔT),for the horizontal velocity,u: u*,for the temperature,T: ΔT=Tb-Tt.


### 2.2.2. The Dimensionless Variables

Based on the above references, the following dimensionless variables can be defined:(7)x′=xlx*,(8)y′=yH,(9)u′=uu*,(10)v′=vv*=vμ[ρ0gβ(ΔT)K],(11)T′=TΔT.
We assume that the dimensionless variables x′ and u′ range in the interval [0,1], even though the references used in their definition are unknowns; this is a hypothesis that will be tested later. Furthermore, it is clear that y′ and v′ also range in the same interval. This property of the dimensionless variables will be used later in the dimensionless governing equation to assign a suitable order of magnitude to the discriminated groups we are looking for and to give a physical meaning and an order of magnitude to the “hidden” references.


### 2.2.3. The Dimensionless Equations

Firstly, we study the mass conservation equation because of its simplicity. Substituting (7)–(10) into (2) yields
(12)(u*∂u′lx*∂x′)+(v*∂v′H∂y′)=0
or
(13)(u*∂u′lx*∂x′)+([ρ0gβ(ΔT)Ky]∂v′μH∂y′)=0.
Reorganizing this equation by grouping all the parameters in the same addend, we can write
(14)(μHu*[ρ0gβ(ΔT)Ky]lx*)(∂u′∂x′)+(∂v′∂y′)=0.
Note that the expressions within the parentheses are dimensionless.Secondly, we will work with the momentum equation (5). Again, using the dimensionless variables, this equation can be written in the form
(15)(KyKx)(u*∂u′H∂y′)-(v*∂v′lx*∂x′)=-(ρKygβμ)((ΔT)∂T′lx*∂x′)
or, making use of (12) and simplifying in the dimensionless form of the momentum equation,
(16)(Kylx*2KxH2)(∂v′∂y′∂x′∂y′)+(∂v′∂x′)=(∂T′∂x′).
Finally, the nondimensionalization of the energy equation (3), which is carried out in a similar way, leads to
(17)(ρfcp,f)(u*lx*)(u′∂T′∂x′)+(ρfcp,f)(v*H)(v′∂T′∂y′)=km,xlx*2(∂2T′∂x′2)+km,yH2(∂2T′∂y′2).


## 2.2.1. The Reference Quantities

Discrimination forces the choice of different characteristic values, according to their spatial direction, to make the variables of vectorial character, velocity (u and v), and coordinates (x and y) dimensionless. On the one hand, as regards the characteristic lengths, although two discriminated lengths appear in the statement of the problem, one of them (L) is not relevant since its value is large enough to neither influence the problem nor rule the solution. It is easy to deduce that the steady state patterns cannot be dependent on L when L is large enough; that is, patterns may cover the totality of the domain or a part of it, but its aspect remains invariable as L increases. As a consequence, L cannot be chosen as a horizontal reference quantity, and, since no more references exist, we assume a “hidden” length (lx*) for this propose, an unknown whose physical meaning and value will be investigated later, after derivation of the dimensionless groups. In short, the reference quantities to make x and y dimensionless are (lx*) and H, respectively. On the other hand, as regards velocities, the vertical reference can be derived directly from the first addend of the left part of (1) making Δρ=Δρmax⁡, since the set of parameters “(Ky/μ)(Δρmax⁡)g” provides a characteristic value related to the vertical coordinate. Finally, we assume a new hidden quantity, ux*, as the horizontal reference velocity since this cannot be established directly or indirectly by the set of parameters contained in the statement. The length and velocity references are as follows:for thex coordinate: lx*,for they coordinate: ly*=H,for the vertical velocity,v: v*=-(Ky/μ)(Δρ)g=-(Ky/μ)ρ0gβ(ΔT),for the horizontal velocity,u: u*,for the temperature,T: ΔT=Tb-Tt.


## 2.2.2. The Dimensionless Variables

Based on the above references, the following dimensionless variables can be defined:(7)x′=xlx*,(8)y′=yH,(9)u′=uu*,(10)v′=vv*=vμ[ρ0gβ(ΔT)K],(11)T′=TΔT.
We assume that the dimensionless variables x′ and u′ range in the interval [0,1], even though the references used in their definition are unknowns; this is a hypothesis that will be tested later. Furthermore, it is clear that y′ and v′ also range in the same interval. This property of the dimensionless variables will be used later in the dimensionless governing equation to assign a suitable order of magnitude to the discriminated groups we are looking for and to give a physical meaning and an order of magnitude to the “hidden” references.


## 2.2.3. The Dimensionless Equations

Firstly, we study the mass conservation equation because of its simplicity. Substituting (7)–(10) into (2) yields
(12)(u*∂u′lx*∂x′)+(v*∂v′H∂y′)=0
or
(13)(u*∂u′lx*∂x′)+([ρ0gβ(ΔT)Ky]∂v′μH∂y′)=0.
Reorganizing this equation by grouping all the parameters in the same addend, we can write
(14)(μHu*[ρ0gβ(ΔT)Ky]lx*)(∂u′∂x′)+(∂v′∂y′)=0.
Note that the expressions within the parentheses are dimensionless.Secondly, we will work with the momentum equation (5). Again, using the dimensionless variables, this equation can be written in the form
(15)(KyKx)(u*∂u′H∂y′)-(v*∂v′lx*∂x′)=-(ρKygβμ)((ΔT)∂T′lx*∂x′)
or, making use of (12) and simplifying in the dimensionless form of the momentum equation,
(16)(Kylx*2KxH2)(∂v′∂y′∂x′∂y′)+(∂v′∂x′)=(∂T′∂x′).
Finally, the nondimensionalization of the energy equation (3), which is carried out in a similar way, leads to
(17)(ρfcp,f)(u*lx*)(u′∂T′∂x′)+(ρfcp,f)(v*H)(v′∂T′∂y′)=km,xlx*2(∂2T′∂x′2)+km,yH2(∂2T′∂y′2).


## 3. Resulting Dimensionless Groups and Solutions

The dimensionless groups that rule the steady state solution patterns of the problem must be investigated in light of the final dimensionless equations (14), (15), and (17). Starting with the first, the following monomial containing two unknowns (lx* and u*) emerges:
(18)π1=(μHu*[ρ0gβ(ΔT)Ky]lx*).
The physical meaning of this monomial is clear since it comes from the mass conservation equation; it refers to the mass balance in steady state within the domain defined by the reference quantities lx* and H, but we can investigate even deeper the approximate value of π1. To do that, we estimate the order of magnitude (~) of each derivative term appearing in (14). Making use of the range assigned above to the dimensionless variables, we can confirm that
(19)(∂u′∂x′)~Δu′Δx′~1,(∂v′∂y′)~Δv′Δy′~1,
a result that, when substituted in (14), provides the order of magnitude of the monomial
(20)π1=(μHu*[ρ0gβ(ΔT)Ky]lx*)~1.
This means that, in terms of the reference quantities,
(21)(u*lx*)~(v*H).
Let us continue with the dimensionless momentum of (16). This provides the new monomial
(22)π2=(Kylx*2KxH2)
whose meaning is a kind of hydraulic permeability ratio corrected by the domain geometry parameters lx* and H. It is interesting to note that the pure ratio Kx/Ky is dimensionless from the classical point of view; however, when the discrimination is assumed, this ratio is not a dimensionless group. Besides, based on similar reasoning as regards the order of magnitude of π2, it can immediately be asserted that
(23)(∂v′∂y′∂x′∂y′)~(∂v′∂x′)~(∂T′∂x′)~1
and, consequently,
(24)π2=(Kylx*2KxH2)~1.
Finally, to derive the monomials of energy equation (17), we will divide all the terms by, for example, the last one (km,y/H2), rewriting the equation in the form
(25)(u*lx*)(ρfcp,fH2km,y)(u′∂T′∂x′)+(v*H)(ρfcp,fH2km,y)(v′∂T′∂y′)=(km,xH2km,ylx*2)(∂2T′∂x′2)+(∂2T′∂y′2).
This provides three new dimensionless monomials
(26)π3=(u*lx*)(ρfcp,fH2km,y),π4=(v*H)(ρfcp,fH2km,y),π5=(km,xH2km,ylx*2),
two of which, π3 and π4, based on (21), are identical. Particularly, π3 and π4 are a kind of discriminated Rayleigh number. Besides, from the arguments used for the conservation and momentum equations, we can declare that these monomials must also be of the order of magnitude unity.In short, the problem is ruled by three dimensionless groups,π2, π3 (or π4), and π5. Because of their properties of being dimensionless and of the order of magnitude unity, these monomials can be combined by mathematical manipulation and presented in different forms. Doing this, since we are interested in the unknown lx* appearing in only one group for this variable to be expressed as a function of the other parameters, we will change the monomial π5 for that resulting from the product π5π2, which does not contain lx*. The final renamed monomials are
(27)πI=(Kylx*2KxH2),πII=(Kykm,xKxkm,y),πIII=([ρ0gβ(ΔT)Ky]Hμαm,y),
where α (m2/s) is the thermal diffusivity and αm,y=km,y/ρfcp,f. Now, the solution for lx* is
(28)lx*=HKxKyΨ[(Kykm,xKxkm,y),([ρ0gβ(ΔT)K]Hμαm,y)]
with Ψ an unknown function of the arguments πII and πIII.


## 4. Validation

Several tests to validate the results of the above section are presented. These consist of a selection of scenarios (or cases) for which some of the parameters are suitably changed, retaining the same values for the final dimensionless groups (πII and πIII) that rule the problem or determine its steady state (flow and temperature) patterns. The complexity of the problem requires numerical techniques for the solution and we have chosen the network method as the most suitable tool. This method whose basis can be found in González-Fernández [21] has been applied for the solution of a large number of nonlineal problems in different fields of science and engineering such as heat transfer [22], vibrations [23], electrochemistry [24], elasticity [25], and others [26]. The network model of the problem is run on the software Pspice [27]. For all simulations presented in this section, a grid of 80 (horizontal) × 20 (vertical) volume elements was chosen.For a better illustration of the results, the dependent variable velocity has been substituted by the stream function variable,ψ (m2/s); so, the steady state patterns of this variable directly show the way that fluid particles flow. The connection between velocity and stream function variables is given by the expressions
(29)u=(∂ψ∂y),v=-(∂ψ∂x).
The first simulation extends to twelve cases, for all of which the discriminate monomials πII and πIII, (27), retain the same value. They are organized, in turn, into four groups (three cases per group) in such a way that the parameters H, Kx, and Ky are forced to change, retaining the same value for H(Kx/Ky)1/2 in each group. In this form the characteristic length, given by (28), must have the same value for each group. Table 1 shows the value of the parameters of all cases in terms of their meaningful changes in relation to case 1. The geometry of the scenario is L=360 (m) and H=90 (m), while the temperatures are Tb=1°C,  Tu=0°C; finally, the stream function at all boundaries has the constant value ψb=ψu=ψl=ψr=0, allowing the impermeability condition to be applied.Table 1
Values of the physical parameters, S.I. (cases 1–12).
Case
K
x
K
y
α
x
α
y
μ
Δ
ρ
g
v
*
Group 1
1
K
x
,
1
K
y
,
1
α
x
,
1
α
y
,
1
µ
1
δ
ρ
1
g
1
2.645
e
-
5
5.75
e
-
13
2.3
e
-
12
7.5
e
-
6
3
e
-
5
2
e
-
4
230
10
2
3
K
x
,
1
3
K
y
,
1
2
α
x
,
1
2
α
y
,
1
1.5
µ
1
2
δ
ρ
1
0.5
g
1
5.29
e
-
5
1.725
e
-
12
6.9
e
-
12
1.5
e
-
5
6
e
-
5
3
e
-
4
460
5
3
0.5
K
x
,
1
0.5
K
y
,
1
1.5
α
x
,
1
1.5
α
y
,
1
2
µ
1
2
δ
ρ
1
3
g
1
3.968
e
-
5
2.875
e
-
13
1.15
e
-
12
1.125
e
-
5
4.5
e
-
5
4
e
-
4
460
30
Group 2
4
2
K
x
,
1
K
y
,
1
2
α
x
,
1
α
y
,
1
µ
1
δ
ρ
1
g
1
2.645
e
-
5
1.15
e
-
12
2.3
e
-
12
1.5
e
-
5
3
e
-
5
2
e
-
4
230
10
5
4
K
x
,
1
2
K
y
,
1
4
α
x
,
1
2
α
y
,
1
2
µ
1
δ
ρ
1
2
g
1
5.29
e
-
5
2.3
e
-
12
4.6
e
-
12
3
e
-
5
6
e
-
5
4
e
-
4
230
20
6
K
x
,
1
0.5
K
y
,
1
α
x
,
1
0.5
α
y
,
1
0.5
µ
1
0.5
δ
ρ
1
g
1
1.323
e
-
5
5.75
e
-
13
1.15
e
-
12
7.5
e
-
6
1.5
e
-
5
1
e
-
4
115
10
Group 3
7
5
K
x
,
1
K
y
,
1
5
α
x
,
1
α
y
,
1
µ
1
δ
ρ
1
g
1
2.645
e
-
5
2.875
e
-
12
2.3
e
-
12
3.75
e
-
5
3
e
-
5
2
e
-
4
230
10
8
10
K
x
,
1
2
K
y
,
1
10
α
x
,
1
2
α
y
,
1
0.5
µ
1
0.5
δ
ρ
1
g
1
5.29
e
-
5
5.75
e
-
12
4.6
e
-
12
7.5
e
-
5
6
e
-
5
1
e
-
4
115
10
9
2.5
K
x
,
1
0.5
K
y
,
1
25
α
x
,
1
5
α
y
,
1
0.5
µ
1
δ
ρ
1
5
g
1
1.323
e
-
5
1.438
e
-
12
1.15
e
-
12
1.875
e
-
4
1.5
e
-
4
1
e
-
4
230
50
Group 4
10
K
x
,
1
3
K
y
,
1
α
x
,
1
3
α
y
,
1
µ
1
δ
ρ
1
g
1
7.935
e
-
5
5.75
e
-
13
6.9
e
-
12
7.5
e
-
6
9
e
-
5
2
e
-
4
230
10
11
2
K
x
,
1
6
K
y
,
1
α
x
,
1
3
α
y
,
1
2
µ
1
2
δ
ρ
1
0.5
g
1
7.935
e
-
5
1.15
e
-
12
1.38
e
-
11
7.5
e
-
6
9
e
-
5
4
e
-
4
460
5
12
(
1
/
3
)
K
x
,
1
K
y
,
1
2
α
x
,
1
6
α
y
,
1
0.5
µ
1
δ
ρ
1
3
g
1
1.587
e
-
4
1.92
e
-
13
2.3
e
-
12
1.5
e
-
5
1.8
e
-
4
1
e
-
4
230
30With the values of Table1, where α=km/ρfcp,f is the thermal diffusivity, πII=1 and πIII=79.4 for all cases (the high value of πIII will be justified later). As regards πI, we expect it to retain the same value in all cases and its order of magnitude to be unity.Figures2 and 3 show the temperature and stream function patterns for the cases 1 to 3 (Group 1). As expected, the temperature patterns, which extend to half of the domain, are identical for the three cases; also, the stream function patterns, which extend to the two cells adjoining the left side, are identical in aspect, but their values of the stream function are proportional to the Darcy velocity, a parameter that differs for each case (Table 1).Figure 2
Temperature patterns for cases 1–3 (Group 1).Figure 3
Stream function pattern for cases 1–3 (Group 1).Similar results, reflecting the expected patterns, are obtained by numerical simulation for the three other groups. On the one hand, Figures4, 6, and 8 show the temperature patterns for Groups 2 (cases 4–6), 3 (cases 7–9), and 4 (cases 10–12), respectively; in all the simulations identical patterns are obtained for the three cases of each group. On the other hand, Figures 5, 7, and 9 show the stream function patterns of the respective groups; as occurring in the cases of Group 1, these patterns retain the same aspect, but their stream function values once again depend on the Darcy velocity. Note that cases 10 and 11 (of group 4) have the same stream function pattern due to their identical Darcy velocity.Figure 4
Temperature patterns for cases 4–6 (Group 2).Figure 5
Stream function pattern for cases 4–6 (Group 2).Figure 6
Temperature patterns for cases 7–9 (Group 3).Figure 7
Stream function pattern for cases 7–9 (Group 3).Figure 8
Temperature patterns for cases 10–12 (Group 4).Figure 9
Stream function pattern for cases 10–12 (Group 4).Table2 shows the number of cells and the characteristic length resulting from each simulation. Since the monomials πII and πIII retain the same value (79.35 and 1, resp.) for all cases, the value of the unknown function Ψ[πII,πIII] must also be the same and this, according to (28), provides the same value for lx* in each group. In turn, within each group, the factor H(Kx/Ky)1/2 remains constant and, as a consequence, both the characteristic length and the number of cells of the domain must also retain a constant value. A fundamental result of discrimination is that the unknown function of the arguments formed by the dimensionless groups of a problem, except for the group that contains the unknown variable in which we are interested, has an order of magnitude unity regardless of the particular values of the arguments. The numerical simulations confirm this result: Ψ(πII=79.35,πIII=1)=0.8, 0.71, 0.89, and 0.87 (all of the order of magnitude unity) for Groups 1, 2, 3, and 4, respectively.Table 2
Number of cells and characteristic lengths for cases 1–12.
Group
1
2
3
4
H
(
K
x
/
K
y
)
1
/
2
45
63.64
100.62
25.98
Number of cells
10
8
4
16
l
*
x
36
45
90
22.5
Ψ
[
π
II
,
π
III
]
=
π
I
0.87
0.71
0.89
0.87The second simulation extends to ten cases for whichπII=165.6 and πIII=1. These are organized in five groups (two cases for each group) with constant values for the factor H(Kx/Ky)1/2 within each group. The scenario has the same geometry and boundary conditions as those of the first simulation. Table 3 shows the parameters for all cases and any meaningful change compared with case 13.Table 3
Values of the physical parameters, S.I. (cases 13–22).
Case
K
x
K
y
α
x
α
y
μ
Δ
ρ
g
v
*
Group 5
13
K
x
,
13
K
y
,
13
α
x
,
13
α
y
,
13
µ
13
Δ
ρ
13
g
13
5.520
e
-
6
1.94
e
-
13
4.8
e
-
13
1.21
e
-
6
3
e
-
6
2
e
-
4
230
10
14
K
x
,
13
K
y
,
13
2
α
x
,
13
2
α
y
,
13
0.5
µ
13
0.5
Δ
ρ
13
2
g
13
1.104
e
-
5
1.94
e
-
13
4.8
e
-
13
2.42
e
-
6
6
e
-
6
1
e
-
4
115
20
Group 6
15
2
K
x
,
13
K
y
,
13
2
α
x
,
13
α
y
,
13
µ
13
Δ
ρ
13
g
13
5.520
e
-
6
2.743
e
-
13
4.8
e
-
13
1.711
e
-
6
3
e
-
6
2
e
-
4
230
10
16
2
K
x
,
13
K
y
,
13
2
2
α
x
,
13
2
α
y
,
13
0.5
µ
13
0.5
Δ
ρ
13
2
g
13
1.104
e
-
5
2.744
e
-
13
4.8
e
-
13
3.422
e
-
6
6
e
-
6
1
e
-
4
115
20
Group 7
17
2
K
x
,
13
K
y
,
13
2
α
x
,
13
α
y
,
13
µ
13
Δ
ρ
13
g
13
5.520
e
-
6
3.886
e
-
13
4.8
e
-
13
2.425
e
-
6
3
e
-
6
2
e
-
4
230
10
18
2
K
x
,
13
K
y
,
13
4
α
x
,
13
2
α
y
,
13
0.5
µ
13
0.5
Δ
ρ
13
2
g
13
1.104
e
-
5
3.886
e
-
13
4.8
e
-
13
4.85
e
-
6
6
e
-
6
1
e
-
4
115
20
Group 8
19
4
K
x
,
13
K
y
,
13
4
α
x
,
13
α
y
,
13
µ
13
Δ
ρ
13
g
13
5.520
e
-
6
7.772
e
-
13
4.8
e
-
13
4.85
e
-
6
3
e
-
6
2
e
-
4
230
10
20
4
K
x
,
13
K
y
,
13
8
α
x
,
13
2
α
y
,
13
0.5
µ
13
0.5
Δ
ρ
13
2
g
13
1.104
e
-
5
7.772
e
-
13
4.8
e
-
13
9.7
e
-
6
6
e
-
6
1
e
-
4
115
20
Group 9
21
0.5
K
x
,
13
K
y
,
13
0.5
α
x
,
13
α
y
,
13
µ
13
Δ
ρ
13
g
13
5.520
e
-
6
9.715
e
-
14
4.8
e
-
13
6.0625
e
-
7
3
e
-
6
2
e
-
4
230
10
22
0.5
K
x
,
13
K
y
,
13
α
x
,
13
2
α
y
,
13
0.5
µ
13
0.5
Δ
ρ
13
2
g
13
1.104
e
-
5
9.715
e
-
14
4.8
e
-
13
1.21
e
-
6
6
e
-
6
1
e
-
4
115
20Figures10 and 11 show the temperature and stream function patterns for cases 13 to 14 (Group 5). Again, as expected, the temperature patterns are the same, while the stream function patterns, despite retaining the same aspect, provide a net flow velocity proportional to the Darcy velocity. The patterns of the rest of the groups are shown in Figures 12-13 (Group 6), Figures 14-15 (Group 7), Figures 16-17 (Group 8), and Figures 18-19 (Group 9).Figure 10
Temperature patterns for cases 13-14 (Group 5).Figure 11
Stream function pattern for cases 13-14 (Group 5).Figure 12
Temperature patterns for cases 15-16 (Group 6).Figure 13
Stream function pattern for cases 15-16 (Group 6).Figure 14
Temperature patterns for cases 17-18 (Group 7).Figure 15
Stream function pattern for cases 17-18 (Group 7).Figure 16
Temperature patterns for cases 19-20 (Group 8).Figure 17
Stream function pattern for cases 19-20 (Group 8).Figure 18
Temperature patterns for cases 21-22 (Group 9).Figure 19
Stream function pattern for cases 21-22 (Group 9).Table4 organizes the rest of the information: characteristic length, total number of cells, factor H(Kx/Ky)1/2, and Ψ[πII,πIII]. It is interesting to note that, despite the appreciable change in the monomial π2 from 79.35 to 165.6, the value of the function Ψ[πII,πIII] scarcely changes from approximately 0.80 (for π2=79.35) to 0.55 (for π2=165.6), determining the same order of magnitude for the monomial πI. Of course lx* depends on the factor H(Kx/Ky)1/2 and has the same value for both cases of each group.Table 4
Number of cells and characteristic lengths for cases 13–22.
Group
5
6
7
8
9
H
(
K
x
/
K
y
)
1
/
2
57.26
68.043
80.98
114.52
40.49
Number of cells
12
10
8
6
16
l
*
x
30
36
45
60
22.5
Ψ
[
π
II
,
π
III
]
=
π
I
0.524
0.529
0.556
0.524
0.556A comparison between the results for Groups 2 and 6 confirms the relatively small influence of the monomialπ2 (a kind of discriminated Rayleigh number) in the length of the cell, compared with the appreciable change in π2. While π2 changes from 79.35 to 165.6, lx* changes from 45 to 36 m.To finish, we will return to the question why we assumed large values forπ2 (the discriminated Ra number). Although the scientific literature for isotropic scenarios uses a wide range of values (from ≈40 to more than 1000; indeed, convection in isotropic media takes starts around π2=40), according to the formal protocol used for the derivation of the dimensionless groups, these values should be of the order of magnitude unity. Further investigations are needed to throw light on this point.In complex coupled problems such as the application studied in this paper, besides the global balances occurring in the complete domain, new balances could take place in small regions or subdomains, in which neither the geometrical parameters nor the governing equations are necessarily the same as those that rule the global problem. While the governing equations of a subdomain are a simplification of the governing equations of the complete problem and cannot provide new dimensionless numbers through nondimensionalization, the geometrical parameters (the fractions of the total lengths that define the subdomain) should reduce the order of magnitude of a monomial from its high value (when referring to the complete domain) to unity (when applied to a suitable subdomain), especially if the geometrical parameters appear in the monomial as a potential dependence. Indeed, this is the most probable explanation to justify the high controversial value ofπ2 in this problem.


## 5. Conclusions

Despite the fact that it is currently applied by researchers to investigate the dimensionless groups that determine the solution patterns of a given problem, classical nondimensionalization may not lead to the most precise solution, that is, to the expressions and number of dimensionless groups that independently rule the solution, particularly in complex problems. This is the case, for example, in conjugate processes with anisotropic domains. In these cases, we have to take into account the use of discrimination in the process of nondimensionalization. The formal application of this technique considers the quantities of vectorial character (lengths, velocities, etc.) as if they were independent of each spatial direction and uses suitable references (one for each spatial direction) to make each one of the dependent and independent variables of the problem dimensionless. In addition, the values of the references must allow the variables to range in the interval[0,1]. If no information is given at the statement of the problem regarding one of these references, it must be established as if the references were known.By doing this, the resulting groups, besides ruling the solution of the problem, have two interesting inherent properties: they are of the order of magnitude unity and have a clear physical meaning in terms of balance of the quantities that balance in the process. In addition, since, in terms of the dimensionless groups, the solution can be expressed by putting one group in terms of the rest by means of an unknown function of their arguments, it is always possible to determine the order of magnitude of the above-mentioned unknown reference.In the context of the problem of natural convection between horizontal plates heated from below, in a hydraulic and thermal anisotropic medium, three dimensionless groups emerged that could not be obtained without discrimination. The order of magnitude of the characteristic length of the typical cell is directly derived from these groups.



---

*Source: 796781-2014-06-22.xml*