Numerical Calculations for a Boundary Layer Flow past a Moving Vertical Porous Plate with Suction/Injection and Thermal Radiation

The work presented herein investigates the velocity, heat transfer, Nusselt number and skin friction profiles involved in boundary layer flow past a moving vertical porous plate. Similarity transformations are employed to convert the governing nonlinear unsteady momentum and energy equations from their partial differential equation forms to boundary value ordinary differential equations. The resulting equations are then solved numerically by the Runge-Kutta fourth order method with the help of a shooting technique. Several features of the flow and heat transfer characteristics for different values of problem parameters are analyzed and discussed. These include the effects of the radiation parameter (R), suction and injection parameter (c), Grashof (Gr) and Prandtl (Pr) numbers on the flow and heat profiles. Numerical results show the impact of blowing and sucking as well as radation on boundary layer flows of this type. Both the skin frictions as well as the heat transfer rate are also significantly related to the radiation parameter. For all these cases; the numerical results are found to be in agreement with the physics of the problem.


Introduction
Fluid flow and heat transfer for both moving and stationary media have been studied in the past because of its relevance to thermal insulation, waste disposal, How to cite this paper: Onyejekwe, O.O., Belete, Y.D. and Worku, N.A (2020) Numerical Calculations for a Boundary Layer Flow past a Moving Vertical Porous Plate with Suction/Injection and Thermal Radia-drying of agricultural products, oil recovery, nuclear reactors, geothermal reservoirs and other practical applications in the areas of engineering and geophysical fields. In addition, cases related to boundary layer flow and heat transfer induced by moving surfaces are vastly significant in extrusion of polymers, glass blowing, metal and fiber spinning. A comprehensive list of relevant work in this area can be found in a monograph by Nield and Bejan [1].
Bejan and Khair [2] studied heat and mass transfer effects for the case of boundary layer flow embedded in porous medium. Related studies involving series expansion, suction and blowing for the case of a vertical plate in a saturated porous medium were carried out by Lai and Kulacki [3] [4] and Raptis [5]. A slightly different study involving not only boundary layer flow but air entrainment over a continuous solid surface moving with a constant velocity in an ambient fluid was carried out by Sakiadis [6]. His work was further extended other investigators such as Ishak et al. [7] [8], Seini and Makinde [9] and Chaturvedi [10].
Magnetohydrodymaics (MHD) of boundary layer flow across a vertical porous plate has a wide area of industrial applications in such areas as geophysics, geothermal energy extractions, astrophysics, plasma studies and crude oil purifications. With reference to these applications, the overriding classical problem in this field was introduced by Blasius [11], who was the first to deal with boundary layer flow over a flat plate. Sakiadis [6] later considered boundary layer flow over a continuously moving solid surface with a constant speed. Erickson et al. [12] later followed this up with different heat and mass transfer considerations including suction and injection.
Soundalgekar and Raman [13] examined alterative aspects of this problem by considering the effects of wall temperature, mass transfer and various other fluid properties. Other considerations by Reddy et al. [14] include MHD boundary layer flow in a channel involving radiation, absorption and heat dissipation and generation.
Other interesting work in this area involves the interaction of MHD boundary layer flow and moving objects. Similarity solutions of MHD flow past a vertical plate with convective boundary condition were carried out by Israel-Cookey et al. [15]. Mostafa and Mahmoud [16] investigated the overall effect of hydromagnetic flow across a moving vertical plate embedded in a porous medium. Ramachandran and Bhaksar [17] extended this study by including other effects such as radiation, and mass transfer. An unsteady MHD boundary layer flow of a rotating fluid past a vertical porous plate was studied by Reddy [18]. He applied a perturbation technique to arrive at the effects of pertinent fluid parameters on fluid flow. Relevant work in this area is reported by Makinde et al. [19] and Pantokratoras [20].
In this paper we consider the effect of suction and injection, thermal radiation and other dimensionless flow variables on boundary layer flow across a moving vertical porous plate. The nonlinear coupled partial differential equation describing the temperature and velocity profiles are transformed into ordinary dif-ferential equations using similarity variables and are solved by the fourth order Runge-Kutta shooting technique. The effects of various flow and heat transfer parameters are presented graphically. This includes skin friction coefficient, the Nusselt number as well as the temperature and the velocity profiles. Followed by his is a thorough quantitative description of each of the profiles. This is essential in order to confirm that the solutions obtained are in consonance with the physics which the problem formulations represent.

Problem Formulation
We consider an unsteady incompressible (constant density) viscous fluid with thermal radiation past a vertical moving porous plate. The x-axis is taken as the direction along the plate in the vertical upward direction with the y-axis perpendicular to it. Let u and v be the components of the velocity along the x and y axis.
The governing equations are represented as The above equations represent the continuity, momentum and energy equations respectively. The parameter µ is the kinematic viscosity, k is the thermal conductivity. , , , , , r q g T t T β ∞ ; are the volumetric expansion, radiative heat flux, gravity, temperature, time, and wall temperature respectively. We apply Roseland approximation, to represent the radiative heat flux as: where , s σ are the Stefan-Boltzmann constant and the mean absorption. Taking the temperature differences to be small, we use the Taylor series expansion, to linearize 4 T about T ∞ to yield: Given Equation (4) and Equation (5) Equation (3) can now be written as: The boundary conditions for the velocity and temperature fields at the plate surface and far into the fluid are given as: where , w T T ∞ are the wall and ambient temperatures. In order to write the governing equations and the boundary conditions in dimensionless forms the fol- Substituting Equation (9) into Equation (2) and Equation (3) yields the following nonlinear coupled ordinary differential equations: The corresponding non-dimensional boundary conditions are given as: For the type of boundary layer flow under consideration, both the Nusselt number and the skin friction coefficient are important flow parameters and are described as follows: where w τ is the shear stress along the plate. The temperature field facilitates the computation of the heat transfer coefficient and is expressed non-dimensionally as the Nusselt number.
The shooting method can be applied to solve the Equation (10) and Equation (11) by replacing each of them by two first order ODEs with known functions (Dirichlet or Neumann). Nonlinear terms do not pose any special problem with the shooting technique; however a good initial guess is helpful. The following equations further elucidate the application of this scheme for this study. For Equation (10) To solve these system of IVPs (initial value problems), we need P(0) and z(0). i.e.

Results and Discussion
The primary physical quantities of interest here are, the velocity and temperature profiles, skin friction coefficient and the Nusselt number. In order to arrive at a physical insight into the problem formulation a parametric study is carried out to illustrate the effects of the different governing parameters namely the   The effect of suction on the temperature profile is shown in Figure 2. It can be seen that the temperature decreases with increase in value of the suction parameters while other parameters kept constant. Figure 3 and Figure 4 show that   suction and radiation parameters have a profound effect on the velocity profiles especially very close to the plate. For example there is a rise in the velocity profile with increase in radiation parameter. The opposite is the case for the suction parameter. Most of these observation happen very close to the plate surface and contribute to the overall shapes of the profiles. It can be observed that maximum velocity occurs a short distance from the plate where the overall effect of the plate movement is the greatest and decreases exponentially farthest from the plate.
On the whole it can be seen that the velocity increases as the mass suction increases. The resultant effect on the momentum boundary layer is to make it thinner On the other hand an increase in mass suction decreases the temperature profile.
Other parameters of interest are the skin friction coefficient and Nusselt numbers. Figure 5 shows that an increase in the radiation parameter (R) and suction brings about a decrease in the wall heat transfer rate with other parameter fixed. The asymptotic shape of the graphs suggests that there is a limiting value of (R) that will not result in any change in the wall heat transfer rate. Figure 6 displays an opposite effect for the skin friction for the same set of conditions. An overall increase in rates of skin friction is displayed for the same parameter combinations. Figures 7-10 show the effect of fluid or mass injection and suction. into the boundary layer flow for all the flow and heat transfer parameters considered above.     and radiative heat transfer result in a decrease in the heat transfer on the plate. Figure 10 shows an opposite effect for the case of fluid injection. The same trend is observed in Figure 9.

Coclussion
Numerical computations have been carried out to ascertain the impact of vari- ous physical parameters on fluid flow involving a heat transfer process and the effects have been found to be quite significant. There is enough evidence to confirm that boundary layer profiles exist for as many physical considerations as possible. More investigations are further needed to verify the effects of those physical parameters that were not given full considerations in this study given the relevance of this type of flow in industrial applications. For example, the fact that both the radiation and suction parameters have a profound influence on the rate of heat transfer as well as the skin friction coefficient is an observation that evokes further research.