^{1}

^{*}

^{2}

^{1}

The aim of this work is the study of the magnetohydrodynamic (MHD) unsteady free convective flow of water near 4 ℃ past an infinitely vertical plate moving with constant velocity. The influence of constant uniform suction was also considered. The partial differential equations (PDEs) and their initial and boundary conditions, describing the problem under consideration, are dimensionalized and the numerical solution is obtained by using the finite volume discretization methodology which is suitable for Fluid Mechanics applications. The numerical results for the velocity and temperature fields are shown in figures for different dimensionless parameters entering in the problem under consideration, such as the magnetic parameter, M and the Grashof number, Gr . This study predicts the effects of a constant magnetic field and uniform suction on the free convective flow of water near 4 ℃ , when the water is electrically conductive. Analysis of the results showed that the velocity and temperature profiles are noticeably influenced by these parameters.

Free convection flow past an infinite vertical plate is an important application from a technological point of view. It becomes a more attractive problem when the fluid is water near 4˚C, electrically conductive, and the flow is subjected to a transverse and constant magnetic field.

It is known that for a fluid like air or water at ordinary temperature and atmospheric pressure the variation of the density with the variation of the temperature is given by

where at 20˚C. However, for temperature variations of magnitude ±4˚C away from 4˚C, the variations in density are very closely given by

with equal to. From the above it is apparent that for small temperature variations, free convection in water near 4˚C would be different from that at 20˚C [

Many researchers have studied the steady free convective flow of water near 4˚C past a vertical plate. Govindarajulu [

The combined convection flow of water near 4˚C through a porous medium bounded by a vertical plate was studied by Raptis and Pop [

The MHD free-convection effects on the oscillatory flow of water near 4˚C past an infinite porous plate was studied by Georgantopoulos et al. [

In this work we consider the unsteady free convective flow of water near 4˚C in the laminar boundary layer over a vertically moving permeable plate, under the influence of a constant transverse magnetic field. The presented results were obtained after dimensionalization of the PDEs using a numerical approach. This approach is based on the finite volume (FV) discretization scheme which is suitable for Fluid Mechanics applications. The discretization was performed with the help of a specialized symbolic package created by the authors in Mathematica.

We consider the unsteady free convective and MHD flow of water near 4˚C past an infinite vertical moving plate with uniform suction. The -axis is taken along the plate in the vertical upward direction and the -axis normal to the plate,

Continuity equation

Equation of motion

Energy equation

where is the velocity component at the -direction, is the normal velocity at the plate,: the time;: the acceleration due to gravity;: the coefficient of thermal expansion;: the kinematic viscosity;: the electrical conductivity;: the magnetic induction;: the fluid density;: the fluid temperature;: the

fluid temperature at infinity;: the specific heat at constant pressure and: the thermal conductivity.

The initial and the boundary conditions are

where u_{0} is the velocity of the plate, is the constant normal suction and is the temperature of the plate. From (3) it is obtained that, for every.

We introduce the dimensionless quantities

Using the dimensionless quantities in (8), Equations (4) and (5) become respectively:

The initial and boundary conditions in dimensionless form are as follows:

Finally, the problem under consideration is described by the system of Equations (9) and (10), subjected to the initial and boundary conditions (11) and (12).

Contrary to the technique used in previous works of the authors for numerically solving the problem at hand [

The analysis begins by obtaining the discretized form of the system of Equations (9)-(10) by using a symbolic package developed by the authors for that purpose. To discretize the coupled set of PDEs the finite volume method on a collocated grid is used (all variables are discretized at the center of the control volume,

where and are the unknown quantities at the center of the control volume as shown in

Having obtained the discretized systems (13)-(14), we construct the system of algebraic equations that constitute the grid for each time step. By grid and time in-

dependence studies with different grid sizes, it is established that the results are time and space independent for equal to 0.05 and equal to.

Then, for each time step, the system is solved algebraically by Mathematica’s function Solve in respect to the grid values of the functions, after first substituting the values of the previous time steps to the system. The procedure is optimized for speed by keeping all the grid values for all the time steps in memory.

The system of equations under consideration is a highly coupled system of two nonlinear equations, whose solution can be tedious. However, with the advancement of numerical techniques such as the finite volume (FV) methodology, which is suitable for Fluid Mechanics applications, problems like this can be solved. A thorough analysis of the problem under consideration includes the study of the velocity and temperature fields under the influence of the dimensionless parameters entering in the problem, such as the magnetic parameter, and the Grashof number,. The Prandtl number was the same for all cases and equal to 11.4.

Initially, the suction velocity was considered equal to zero, and the influence of the dimensionless parameters entering in the problem was studied. Next, the effect of a constant suction velocity, , on the velocity and temperature profiles under the influence of a constant transverse magnetic field was investigated.

The velocity and temperature distributions for different dimensionless times are shown in

The effect of the magnetic parameter, M, on the

velocity field for different times, (and) are shown in Figures 4 and 5. For the magnetic parameter, M, decreases the velocity,

Figures 6 and 7 show the effect of Grashof number, , on the velocity profiles for and. The velocity increased noticeably near the plate as the Grashof number increased, due to the increase of the buoyancy forces compared to the viscous forces acting on the fluid (water) near 4˚C.

The effect of a constant suction velocity on the velocity and temperature fields under the influence of a constant transverse magnetic field was also studied. The effects of suction velocity upon the velocity and temperature fields for dimensionless time t = 0.25 are shown in Figures 8

and 9. It is evident from the figures that as the absolute value of suction velocity increases, the profiles of the velocity and temperature decrease. These findings are in accordance with previously published results [

The unsteady free convective flow of water near 4˚C, past a vertical plate, moving with constant velocity, and subjected to a transverse magnetic field and constant suction velocity was numerically studied.

In the present paper an analytic/symbolic approach was chosen for solving the system of Equations (9) and (10), subjected to the initial and boundary conditions (11)

and (12). A symbolic package developed in Mathematica is used for obtaining the discretization of the problem and for constructing the algebraic system to be solved for each grid point. Then, for each time step the system was solved analytically using the numerical data of the previous time step, hence avoiding errors in accuracy due to the use of a numerical solver such as a Newton-like method. As a consequence, the obtained solution is more accurate for greater values of t.

Moreover, further development of the symbolic package for the fast and accurate discretization of dynamical systems using the Finite Volume approach (with the option of different grid choices available to the user) will be an invaluable tool for any researcher using the Finite Volumes discretization scheme, especially for systems in multi-dimensions where the calculations by hand are not only time consuming but also error prone.

The results of this study showed that the velocity and temperature fields increase away from the flat plate, as time increases. For, the velocity decreased throughout the boundary layer when the magnetic parameter was increased while it increased noticeably near the plate when the Grashof number was increased. Finally, as the value of the suction velocity increased, the profiles of velocity and temperature decreased.