1. Introduction
Flow in porous media is essential in the metallurgical industry, particularly in the foundry industry, in which molten metals are poured into sand molds. The gases generated during the pouring of the molten metal into the sand molds must be able to escape to the outside so that the solidified pieces become free of pores and blowholes [1] [2] . The property of molding sands that measures their ability to allow gases to pass through them is called permeability. This property depends both on the physical characteristics of the sand (size and shape of the grains, grain size distribution, degree of compaction, porosity, tortuosity, humidity, additives, and so on), and on the physical properties of the gases [3] [4] .
For laminar flow of gases in sands, velocity and pressure gradients have a linear relationship expressed by Darcy’s Law [5] . However, for turbulent flows Darcy’s Law is no longer valid. In the case of turbulent flow, several modifications to Darcy’s Law have been proposed that introduce nonlinearity between gas velocity and pressure gradients. Nonlinear expressions between velocity and pressure gradients previously published in the literature are used for turbulent flows [6] [7] . Reported numerical results show that for Reynolds numbers greater than 1 the deviations from Darcy’s Law begin to be significant and the linearity between the gas velocity and the pressure gradients is no longer valid. In this work, the effect of the Reynolds number on the validity of Darcy’s law in molding sands is numerically studied using the Ergun equation. An expression that explicitly quantifies the dependence between the Reynolds number and the non-Darcy flow fraction of the air flow in molding sands is obtained.
2. Darcy Flow
Parameters such as humidity, additives, porosity, compactness, and tortuosity, determine the permeability of molding sands. Sometimes the Kozeny-Carman equation is proposed to determine the permeability of a porous media based on the parameters mentioned above [8] [9] [10] [11] . The Kozeny-Carman equation is a semi-empirical model that is obtained, among other considerations, from the exact solution for a laminar flow around a round tube [10] . In [12] the following expression is proposed to determine tortuosity τ as a function of porosity ϕ:
(1)
where A = 1 and n = 2 for molding sands. In [13] it is proposed that
(2)
For low flow rates through porous media, i.e. laminar flows, Darcy’s Law establishes a proportional relationship between flow and pressure gradient. Darcy’s Law is a semi-empirical expression which was originally represented as follows [1] [14] :
(3)
where Q is the volumetric flow of fluid through the porous media, k is the specific permeability, A is the cross-sectional area, μ fluid viscosity, ΔP is the pressure drop, and L is the thickness of the porous layer. In vector form, Darcy’s law is expressed as follows:
(4)
(5)
where
is the velocity vector, and
is the pressure gradient.
In [15] it is stated the Darcy’s Law is valid only for flow rates for which the Reynolds number is less than 1, i.e.:
(6)
where u is the mean microscopic velocity, and d is the mean pore diameter. Estimating pore diameter from the expression
(7)
then, in accordance to [15] , the criterion for the validity of Darcy’s Law becomes
(8)
3. Non-Darcy Flow
For turbulent gas flow in molding sands, Darcy’s Law losses applicability since the relationship between velocity and pressure gradients is no longer linear. For porous materials such as molding sands, a macroscopic Reynolds number is defined as [1]
(9)
where u = Q/A is the macroscopic velocity and dp is the particle diameter. From Equation (9),
(10)
To handle turbulent flows, some corrections to Darcy’s Law have been proposed. An example is the Forchheimer equation [16] [17] for spherical particles:
(11)
where x is the flow direction and
(12)
(13)
Ergun equation can be applied for non-uniform size of sand particles and large pressure drops [1] :
(14)
where λ is the shape factor. It is assumed that λ = 1 for spherical particles. The above equation can be rewritten as follows:
(15)
Pressure gradient for turbulent flows in Equation (15) can be calculated by considering the sum of the linear contribution of the velocity given by Darcy’s Law, DF, and the non-linear contribution of the velocity, NDF:
(16)
where
(17)
(18)
The linear and non-linear fraction contributions of velocity to the pressure gradient are defined here in this way:
(19)
(20)
Of course, it is verified that
+
= 1.
Substitution of Equations (10), (17) and (18) into Equation (20) yields
(21)
where
(22)
For laminar flows in which Reynolds number is low,
tends to zero, while for turbulent flows in which Reynolds number is high the value of
tends to 1, therefore
(23)
(24)
Table 1 shows some values of Non-Darcy flow fraction XNDF for several values of the Reynolds number obtained from Equations (21) and (22) for spherical particles and air using λ = 1.0, ϕ = 0.1, ρ = 1.223 kg·m−3, μ = 1.9 × 10−5 kg·m−1·s−1, dp = 1.0 × 10−4 m. It is observed that as the Reynolds number increases, the value of the Non-Darcy contribution of velocity to the flow is increased too. According to Table 1, a Reynolds number equal to 1 gives a contribution of 0.0115, i.e. 1.15%, of the non-Darcy flow to the total air flow through the sand grains. This result seems to be compatible with the one reported in [15] which establishes that Darcy’s law is valid for Reynolds number less than 1.
Table 1. Non-Darcy flow fraction as function of Reynolds number for spherical grains, using λ = 1.0, ϕ = 0.1, ρ = 1.223 kg · m−3, μ = 1.9 × 10−5 kg·m−1·s−1, dp = 1.0 × 10− 4 m .
Figure 1 depicts results of Table 1. Non-Darcy flow fraction of air in molding sands increases as the Reynolds number is increased. Figure 1 also shows that the non-Darcy flow fraction depends strongly on the shape factor of the sand grains. As the shape of the grains moves away from sphericity, that is, as the value of the shape factor increases, the contribution of the non-Darcy flow decreases. Figure 1 suggests that for sand grains with a shape factor greater than 5, the validity of Darcy’s law is still maintained even for Reynolds number much greater than 1.
Figure 2 shows the nonlinear contribution of velocity, i.e. non-Darcy flow, to the air flow for a Reynolds number equal to 1. The nonlinear contribution becomes significantly smaller for shape factors with values greater than 6. The following expression was obtained through a statistical fit of the data shown in Figure 2:
(25)
Figure 1. Non-Darcy flow fraction of air flow through sand particles as function of Reynolds number for several values of the shape factor.
Figure 2. Non-darcy flow fraction of air as a function of the shape factor of the sand grains for a Reynolds number equal to 1.
The red line in Figure 2 depicts an exponential decay of the non-linear flow with the shape factor. For example, for flake graphite with a factor form of 7.96 Equation (25) predicts a value of 0.0022 for non-linear flow fraction, which is extremely low. According to [1] , the shape factor for screened sand grains varies from 1.15 to 2.54, which means that for values of the Reynolds number less than 1, Darcy’s law adequately represents the flow of gases through molding sands.
4. Conclusions
The effect of the Reynolds number on the validity of Darcy’s law in molding sands was numerically studied in this work using the Ergun equation. From numerical results, the following conclusions arise:
1) An expression that explicitly quantifies the dependence between the Reynolds number and the non-Darcy flow fraction of the air flow in molding sands was obtained and reported.
2) For values of the Reynolds number less than 1 the contribution of non-linearity between air velocity and pressure gradient in molding sands with spherical particles is only 1.15%.
3) This percentage of non-linearity decreases even more with the increase in the form factor of the sand grains, reaching 0.27% for a form factor equal to 7.