Range of Applying the Boundary Condition at Fluid/Porous Interface and Evaluation of Beavers and Joseph’s Slip Coefficient Using Finite Element Method

: In this work, Finite Element Method (FEM) is applied to obtain the condition at the boundary of the interface between a channel and a porous medium. The boundary conditions that should be applied to the inhomogeneous interface zone between the two homogeneous regions of free fluid and porous medium are derived. The comparison has been performed for porous material characterizations to provide the velocity at the inhomogeneous interface zone with variable permeability between the two homogeneous regions of free fluid and porous medium. Also, the dependence of the slip coefficient on the thickness of the transition zone is established and the values of the thickness are so justified that the numerical results and the numerical results of our proposed technique are found to be in good agreement with experimental results in the literature. 𝑥 (cid:3036) ( 𝑖 = 1,2,3) . Numerical evaluations show that and M.A.F.; software, M.A.F.; validation, A.H.S. and M.A.F.; formal analysis, A.H.S.; investigation, A.H.S. and M.A.F.; resources, A.H.S. and M.A.F.; data curation, A.H.S. and M.A.F.; writing—original draft preparation, A.H.S.; writing—review and editing, M.A.F.; visualization, A.H.S. and M.A.F.; supervision,


Introduction
Several authors have discussed the boundary condition at the interface between a free fluid and a porous medium [1], boundary conditions at a naturally permeable wall [2], boundary condition at the interface of a porous medium [3], and numerical simulations of pressure jump interface law for Stokes-Darcy coupling [4]. Also, the justification of pore level Navier stokes equations is noticed and discussed in Reference [5] and Reference [6]. Czochra and Mikelic [7] studied the effective pressure interface law for transport phenomena between an unconfined fluid and a porous medium using homogenization. Jäger and Mikelic [8] discussed the interface boundary conditions by analyzing beavers. Jäger et al. [9] studied the asymptotic analysis of the laminar viscous flow over a porous bed. Upscale Navier-Stokes equations are rigorously derived in Reference [10]. The single phase Poiseuille flow over a permeable block is studied in Reference [11] and the boundary conditions that must be applied to the inhomogeneous interface zone between the free fluid and porous medium are derived using the matched asymptotic expansions method, without specifying the porosity-dependent function and permeability-dependent function at the interface zone. The outputs for nano-refrigerant (R600a/oil/CuO) boiling heat transfer within flattened channels utilizing an experimental method are presented in Reference [12]. The authors discussed the influence of a flattened percentage, flow rate, and vapor quality as well as the mass fraction of CuO on boiling heat transfer (h). The Finite element method has been employed in Reference [13] to show that variation of the energy storage efficiency of Copper oxide nanoparticles and V shaped fins is involved in a storage unit to expedite the solidification. In Reference [14], the finite element method is recalled to obtain the outputs, which are the roles of radiation parameter (Rd), Darcy number (Da), nanofluid volume fraction (Φ), Rayleigh number (Ra), and supplied voltage (Δφ).
In this work the authors show the extent to which this condition can be applied. In other words, the range of application of this condition is discussed. The materials suggested by Beavers and Joseph are considered and the thickness of the interface that make the condition applicable is evaluated. The condition deduced here using the technique of Finite Element Method (FEM) is compared with the analytical solution and experimental results obtained by Beavers and Joseph [2], which emphasize the range of validity of this condition.
In this work, the boundary condition at the interface between a free fluid and a porous medium is derived using FEM, which is where ₃ is vertical axis perpendicular to the flow direction, is the velocity in the channel, ₀ is Darcy's velocity, and (0) is the velocity at the upper boundary of the interface. Also ₁ and ₂ (obtained using the Mathematica) depend on the permeability ₀ and the thickness of the interface . Equation (1) is compared with the Beavers and Joseph [2] condition, namely where and are the velocity in the channel and at = 0, respectively, and is the volume flow rate per the unit cross sectional area. Also, is a dimensional constant that distinguishes the geometry of the interface region between the free incompressible fluid and the porous medium (see Reference [11]). The problem which is described above starting from the up-scaled Navier-Stokes equation of the porous medium. Then, the required condition (1) is obtained using the Galerkin finite element scheme [15] and the velocity in the transition zone is deduced. Also, is evaluated and the numerical solution is compared to the experimental results of Beavers and Joseph [2].

Description of the Problem
In order to model the inhomogeneous interface zone between a channel and a porous medium, we shall consider the Cartesian coordinates system , , (as depicted in Figure 1) and unit vectors ₁, ₂, ₃ with ₃ vertically up, where the flow is in the direction ₂. The slab is taken to occupy the region ₃ ≤ 0 ( ₃ → −∞) , the upper channel boundary is taken to be = ℎ which is assumed to be impermeable, and a transition zone is set at (− ≤ ₃ ≤ 0), where is the thickness of the transition zone.
where is a constant pressure gradient in the flow direction, ( ₃) is the permeability, and is the fluid viscosity. Specifically, the flow is considered to be

Channel Flow
In the channel ( ) = 0 and then Equation (3) reduces to Thus where the non-slip condition has been applied at the boundary = ℎ and the slip condition is = at ₃ = 0.

Flow in the Strict Interior of the Porous Medium
The strict interior of the slab is homogeneous, then where ₀ is the slab permeability in the direction flow e₂ and is the thickness of the interface.
Thus, in the strict interior for a slab of infinite thickness the (bounded) solution is is Darcy's velocity.

Flow in the Interface
According to Reference [10], Equation (3) in the interface becomes where ( ) is assumed to be varying as ⁄ , where is the thickness of the interfacial region. As expected, this only implies that the function is rapidly varying in the interfacial region, as expected. The arbitrary constants arise from solving this second order linear differential equation together with the expressions of and in Equation (9) may be evaluated by invoking the global smoothness of and ′.
According to Reference [10], the linear differential Equation (11) has been solved together with the expressions of and in Equation (9), considering the following relations with

Finite Element Formulation
The unknown variable in Equation (14) is approximated over an element as where is the velocity shape function and , are the appropriate boundary conditions of Equation (14). Applying the Galerkin weighted residual [12] criterion results in subject to the continuity of the velocity gradient at the upper and the lower boundaries of the interface. According to the standard Galerkin method and governing equations (Appendix), the global equations are Here ₁ , ₂ and ₃ are the system nodal values at = 0, , − . Solving for and applying ₁ = , the velocity at the upper boundary of the interface yields The expression of the boundary condition (1) is obtained by dividing the interfacial region into at certain values of for each material considered by Beavers and Joseph. The boundary condition (1) can be considered to evaluate in Equation (6): Then, This expression is identical to that calculated by Beavers and Joseph [2] after identifying their parameter with . When the interfacial region is divided into 4 − elements bounded by 5 − values ( = 1 − 5), then the boundary condition is where and are obtained in the same way as and . Accordingly, using the last identification, the boundary condition (1) becomes The boundary condition (21) has the same form of Equation (2) for = .

The Velocity in the Interface
Matching the velocities in the interface with those associated with the strict interior of the porous medium, the constant in Equation (9) is evaluated.
Recalling (15), the velocity in the interface becomes where ₁ and ₂ are evaluated using MATHEMATICA.

Comparisons with Experimental Data and Analytical Solution
According to Beavers and Joseph [2], the fractional increase in mass flow rate due to the presence of the porous wall is: where and are the mass flow rate within the free fluid region for porous medium and impermeable wall, respectively. The obtained result of is The comparison with the experimental results is performed using the characteristics of the porous materials used by Beavers and Joseph in their experiments and recalled in Table 1. These results show that the increase in mass flow rate is sensitive to the value of the .  Table 1 shows that the range of ϵ to apply the boundary condition is obtained by the finite element method and is used to evaluate the fractional of increase in mass flow rate φ. Here, a two-element solution is formulated by taking equally spaced nodes at x = 0, ϵ 2, ϵ ⁄ . Note that the arrangement of  (Figures 6 and 7). With respect to the experiments, it is not possible to find values of ϵ for these materials that allow the curves to pass through all the experimental points simultaneously. However, the appeared values of ϵ give the best fit with the experimental points. The highest values for the thickness ϵ are reached by the Aloxite B. As depicted in Figure 8, the experiment and the numerical results agree for Fometal C at σ < 20 while for Aloxite A is partially agree for 20 < σ < 80. The fourth column in the table shows the range of α corresponding the thickness and compared with values of α proposed by B&J. The comparison between Aloxite A and Fometal C ( Figure 8) shows that the experimental results of the Aloxite A spread over a range of σ much greater than the Fometal C which provides an explicit dependence of φ on the thickness of the interface.  Table 2 shows that the range of ϵ to apply the boundary condition that is obtained by the finite element method. Four element solution is formulated by taking equally spaced nodes at x = 0 , ϵ 4, ϵ 2 ⁄ ⁄ , 3ϵ 4, ϵ ⁄ . A comparison was made between the two-element and four elements in FEM. It is clear from the Figure 9, the agreement between them and the analytical results is at value ϵ = 2.0 × 10 and ϵ = 0.1 × 10 in the first and second method, respectively, for k = 1.1 × 10 . For the permeability k = 1.5 × 10 the curves in both methods apply with the analytical method when the value of epsilon is equal to ϵ = 2.0 × 10 and ϵ = 2.0 × 10 ( Figure 10). Finally, the agreement between them and the analytical results at the value 5.7 × 10 and 5.0 × 10 . in the first and second method respectively at k = 2 .4 8 × 10 ( Figure 7) . Figure 11 shows the relationship between ϵ and α for different values of k. The reverse relationship between ϵ and α values is clear in Figure 11.

Discussion
In this study, we used the materials listed in the table so we could compare the numerical results obtained by the analytical and experimental results conducted by Beavers and Joseph. Table 1 shows the numerical results using the FEM fit with the analytical solution and experimental data of Beavers and Joseph [2] for the Fometal and partially for the Aloxite. It is noticed that the range in case of Aloxite is greater than that for Fometal, which could be due to the irregular nature of Aloxite, as mentioned above. Relation (19) is true only within the range of that appeared in the table. The numerical data indicated that values of Γ₂ range from 21.830 to 35.552 for Aloxite, from 165.663 to 201.374 for Fometal with low permeability, and 295.077 to 1203.67 for Fometal with a high permeability.
As mentioned above, = ₂ is true within the values of the indicated and is compared with the values considered by Beavers and Joseph [2]. Figure 8, where α is plotted against , shows that α decreases with the increase of , which the outcomes of values of from the numerical and experimental results agree with. The comparison with Beavers and Joseph [2] for Fometal , , and Aloxite shows that the values of α, drop within a numerically evaluated range. For example, for Fometal B, the value of the thickness that makes For Tables 1 and 2, alpha has upper and a lower bounds, whereas there is a specific value in the results derived by Beaver and Joseph. Beaver and Joseph indicated that the alpha values differ according to the fluid used in the experiment. Then it is expected that if different types of fluids are used in the same porous material, then alpha values will fall in the specified period obtained by the FEM. Figure 5 compares the curve that represents for the thickness = 2 × 10 when dividing the interfacial region into the two elements in Figure 5a and the four elements in Figure 5b. It illustrates the convergence of numerical computations in both cases to the analytical solution developed by Beaver and Joseph. Figure 7 shows the affinity of the numerical results with the analytical and experimental cases in the case of dividing the interface into two and four elements, but for a different thickness of the interfacial area.
Another discussion concerns the velocity in the interfacial zone. According to relation (25), the velocity profile in the interfacial zone (− ≤ ₃ ≤ 0) is a straight line in ₃ for the range of that appeared in Table 1.

Conclusions
In this work, flow in a channel over a porous medium under a constant pressure gradient is studied. The boundary conditions that should be applied to the inhomogeneous interface zone with variable permeability between the two homogeneous regions of free fluid and porous medium are derived. The up scaled Navier-Stokes equation that is used to describe the flow in the two homogeneous regions is assumed to hold in the inhomogeneous interface zone with variable permeability. The finite element method has been used to solve the problem in the whole domain, where the velocity and the stress are continuous across the interface. The computational results show that the stress condition is related to the thickness of the inhomogeneous interface zone with variable permeability.
Good agreement is obtained between the computational results and the Beavers and Joseph [2] experiment by adjusting the values of the thickness of the transition zone. In particular, it was able to correctly predict the slope of the fractional increase in mass flow rate φ versus the non-dimensional height of the free fluid region for the Aloxite experiments. Despite this, the values of the thickness had to be adjusted to obtain good agreement with the available data.
This study shows that the obtained boundary condition using the Finite Element Method is true at a specified thickness values of the transition zone for both Fometal and Aloxite materials. Furthermore, this methodology is applicable to a wide range of transport phenomena in engineering problems [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32] and of transport phenomena in biological problems [33][34][35]. Finally, an expression of α is obtained and the numerical results on both the permeability and the thickness of the transition zone. Besides, a comparison with the suggested α values of Reference [2] is carried out.