A New Exact Solution for the Flow of a Fluid through Porous Media for a Variety of Boundary Conditions

: The viscous ﬂuid ﬂow past a semi-inﬁnite porous solid, which is proportionally sheared at one boundary with the possibility of the ﬂuid slipping according to Navier’s slip or second order slip, is considered here. Such an assumption takes into consideration several of the boundary conditions used in the literature, and is a generalization of them. Upon introducing a similarity transformation, the governing equations for the problem under consideration reduces to a system of nonlinear partial di ﬀ erential equations. Interestingly, we were able to obtain an exact analytical solution for the velocity, though the equation is nonlinear. The ﬂow through the porous solid is assumed to obey the Brinkman equation, and is considered relevant to several applications.


Introduction
The flow of a fluid through a porous medium has numerous applications in industries dealing with polymer extrusion process, glass blowing, metallurgical processes, and geophysical and allied areas (see [1]). A variety of equations have been used to describe the flow of a fluid through a porous medium as it is one of the important key factors in maintaining the temperature in the medium. These equations due to [2][3][4][5] and others, are merely approximations to the appropriate balance laws. A variety of ideas have been suggested to model the flow of mixtures, and one such approach is that which follows from the seminal works of Darcy and Brinkman and has been given a formal structure by [6,7]. Several specific problems have been solved using such an approach (see [8][9][10][11][12][13][14][15][16][17][18][19]). Here, we study the flow of a fluid through a porous media that is governed by the Brinkman equation (see [20][21][22][23][24][25][26][27][28][29][30]) for a discussion of the status of the Brinkman equation within the context of mixture theory). The fact that we are able to obtain an analytical solution to the problem makes the study all the more interesting. Despite the fact that advanced computing facilities are available to obtain the numerical solution, investigators around the world are much more interested in providing the analytical solution due to their accuracy, relevance, and convenient to analyze physical process, in comparison to numerical solutions. The analytical solution can provide a better assessment of consistency and parameter estimates. Many authors have investigated the fluid flow through porous media and provided analytical solution (see [31][32][33][34][35][36][37]). The novelty of this study is our use of a variety of boundary conditions that subsumes those that have been considered earlier, in addition to new conditions concerning slip and proportional shearing at the boundary.
In this study, we consider the flow of a fluid through a semi-infinite porous media with one boundary subject to the slipping or adherence of the fluid, the solid being proportionally sheared, and the fluid being injected at the boundary (see Figure 1). We are able to obtain an analytical solution by introducing a similarity variable that greatly simplifies the governing equation. The effects of the boundary conditions on the flow through the porous media are determined. conditions that subsumes those that have been considered earlier, in addition to new conditions concerning slip and proportional shearing at the boundary. In this study, we consider the flow of a fluid through a semi-infinite porous media with one boundary subject to the slipping or adherence of the fluid, the solid being proportionally sheared, and the fluid being injected at the boundary (see Figure 1). We are able to obtain an analytical solution by introducing a similarity variable that greatly simplifies the governing equation. The effects of the boundary conditions on the flow through the porous media are determined.

Theoretical Model
Two dimensional laminar, steady, incompressible fluid flow through a porous media is considered. The x-axis is taken along the stretching of the sheet in the direction of the motion, and y-axis is perpendicular to the slit. In order to confine the fluid flow in the region y>0, two forces of equal strength are applied along the x-axis. u and v denote the axial as well as transverse velocities in the flow field. Figure 1 depicts the physical flow problem subjected to the boundary conditions.
We considered the flow of the classical incompressible Navier-Stokes [38] fluid through a porous half-space. We assumed that the equations governing the flow are those given by the Brinkman equation for flow through porous media which assumes that the fluid is incompressible and, hence, the conservation of mass reduces to and the conservation of linear momentum that takes the form where eff μ represents the effective viscosity of the fluid(see [39,40]) provides the definition of the other parameters. The Brinkman equation can be shown to be obtained as a systematic approximation using mixture theory by assuming special structures for the interaction forces between the porous solid and fluid, and assuming the porous solid is rigid (see [41] for a detailed derivation).The transformed governing equations for the conservation of mass and the balance of linear momentum are given as

Theoretical Model
Two dimensional laminar, steady, incompressible fluid flow through a porous media is considered. The x-axis is taken along the stretching of the sheet in the direction of the motion, and y-axis is perpendicular to the slit. In order to confine the fluid flow in the region y > 0, two forces of equal strength are applied along the x-axis. u and v denote the axial as well as transverse velocities in the flow field. Figure 1 depicts the physical flow problem subjected to the boundary conditions.
We considered the flow of the classical incompressible Navier-Stokes [38] fluid through a porous half-space. We assumed that the equations governing the flow are those given by the Brinkman equation for flow through porous media which assumes that the fluid is incompressible and, hence, the conservation of mass reduces to and the conservation of linear momentum that takes the form where µ e f f represents the effective viscosity of the fluid(see [39,40]) provides the definition of the other parameters. The Brinkman equation can be shown to be obtained as a systematic approximation using mixture theory by assuming special structures for the interaction forces between the porous solid and fluid, and assuming the porous solid is rigid (see [41] for a detailed derivation).The transformed governing equations for the conservation of mass and the balance of linear momentum are given as and ν e f f = µ e f f ρ f . The Forchheimer term in the interaction is neglected as it produces little impact on the fluid flow in a porous medium governed by the Brinkman equation (see [42]). Also, the pressure gradient is neglected, and the time factor is zero for the steady case.
The governing boundary conditions are (see [43,44]) Here, d is the parameter of proportional shearing at the boundary, with d 0 and d = 0 corresponding to the boundary at y = 0 and being either proportionally sheared or being fixed. The constants A and B represent the first-and second-order slip coefficients, respectively. Also, the mass transpiration parameter, v c , represents suction or injection depending on v c > 0 or v c < 0, respectively. In order to carry out the analysis, the physical stream functions in terms of similarity variables f and η are introduced as follows: where In terms of physical stream function ψ, the axial and transverse velocities can be rewritten as follows: The Equation (8) satisfies the continuity equation. Upon substitution of Equations (6) and (7) into Equation (4), we obtain Here, the second term indicates the Jacobian and subject. The appropriate boundary conditions are (see [43,44]) as follows: Here, Λ = µ e f f µ is the Brinkman number or viscosity ratio. Using Equations (9) and (10) with Equation (6), the following transformed equation with constant coefficient is derived: The governing boundary conditions for Equation (11) are given as where Γ 1 = A α ν > 0 and Γ 2 = B α ν < 0 are the first-and second-order slip parameters, and K 1 = νφ 2 αK is the reciprocal of Darcy number Da = l 2 K , with l = φ ν α and V c = v c √ α ρ as the mass suction/injection parameters. The subscript denotes the derivative with respect to η.

The Analytical Solution
The flow problem considered is the generalization of the classical works of [43][44][45][46][47][48]. In our problem, the viscous flow with first-and second-order velocity slips over a porous half-space that is stretched or contracted at the boundary (see Figure 1), and the flow being governed by the Darcy-Brinkman model is considered. One obtains nonlinear partial differentiation from Equations (3) and (4), which are mapped into systems of the nonlinear ordinary differential Equation (11), with a constant coefficient by means of similarity transformation subjected to the imposed boundary (12)- (13).The analytical solution for the velocity distribution is determined.
The exact analytical solution of Equation (11) subjected to the governing boundary conditions Equations (12) and (13) is derived. The condition in Equation (13) suggests choosing the equation of the form where β > 0 is to be determined later. Also, A 1 and B 1 are constants that are to be determined by using Equation (12): It follows from Equation (11), (14), and (15) that Here, β > 0 is one of the real roots (see [49,50]). By using the transformation variable ξ = β + a 3 4 , Equation (16) transforms into where p = a 2 − 3 8 a 2 3 , q = A 1 − 1 2 A 2 A 3 + 1 8 A 3 3 , and r = a 0 − 1 4 a 1 a 3 + 1 16 a 2 a 2 3 − 3 256 a 4 3 , and The four corresponding roots of the algebraic Equation (17) are where and Equation (18) gives the complete solution of Equation (17). However, it should be noted that there is only one feasible solution for the Equation (17) when Γ 2 < 0, and based on the flow field Equation (16) has feasible solutions for β > 0.

Results and Discussion
In this paper, we were able to establish the impact of various physical parameters on the velocity distribution. The solutions obtained are in good agreement with that of the classical works when suitably restricted to specific conditions. The main emphasis of this study is the effect of boundary conditions on the flow through porous media. There are several possibilities at the boundaries, namely, the fluid meeting the no-slip adherence condition, the Navier slip condition, the second-order slip condition, as well as the possibility of blowing of the fluid.
The effects of physical parameters such as the mass transpiration parameter (V C ), first-order Navier slip (Γ 1 ), second-order slip (Γ 2 ),Brinkman ratio (Λ), and proportional shearing parameter (K 1 )are discussed graphically. As the velocity distribution is an exponential function with a negative argument, it decreases with the increase in η. Since β is a function of mass transpiration parameter (V C ), first-order Navier slip (Γ 1 ), second-order slip (Γ 2 ), Brinkman ratio (Λ), and proportional shearing parameter (K 1 )both axial as well as transverse velocities are forced to decrease exponentially.
The solution domain of β in Equation (16)has only complex roots when D 1 < 0, only real roots when D 1 > 0,and real repeated roots when D 1 = 0.              Furthermore, under the given slip parameter and mass suction parameter, one can see the increasing velocity boundary thickness with the increase in proportional shearing parameter. Also, the reduction in mass suction leads to the decrease in velocity boundary for other physical parameters fixed. However, for the case of mass injection, the velocity boundary increases with increasing values of slip parameter. Thus, the flow geometry and the rate of change of velocity boundary layer thickness are significantly influenced by the slip parameter.             Figure 8a,b depicts the effect of various physical parameters on the shear stress profile. In all of these plots, there are crossover points for the shear stress profiles, and the combined effects on the porous solid can be observed. The increase in the values of Γ 1 , Γ 2 , and K 1 results in increasing shear at the wall boundary. In the case of mass injection, the shear wall boundary decreases faster for a smaller value of second-order slip parameter. Interestingly, however, the increasing value of Brinkman ratio leads to decreasing shear wall boundary as seen in Figure 8a.

Concluding Remarks
In conclusion, the viscous fluid flow past a porous solid wherein the flow is governed by the Brinkman equation with first-and second-order slip in the presence of mass transpiration was solved, and the exact analytical solution for the governing nonlinear partial differential equation was obtained. The solution was analyzed for the effect of slip parameters, the mass transpiration parameter, the Brinkman ratio, and the extent of shearing or contraction. In the case of the boundary contracting, the solution branches (Figure 5a,b), whereas in the case of the boundary stretching there is only one branch of the solution, and depending on the mass transpiration parameter, the solution is branched (Figure 5a,b).