Magnetohydrodynamics (MHD) is concerned with the study of the interaction between magnetic fields and fluid conductors of electricity. The body force acting on the fluid is the Lorentz force that arises when electric current flows at an angle to the direction of an impressed magnetic field. Due to increasing applications of MHD, researchers have been attracted to study the various types of flow to analyse the velocity, temperature and concentration profiles for different dimensionless parameters. Most of the research papers on MHD flow are basically focussed on the analytic solution; however, the accuracy of the analytic solution is very weak. Even though a number of computational tools are available to solve non-linear coupled partial differential equations, accuracy is the problem with the computational tools. To find a solution with the desired degree of accuracy, mathematical modelling has become an important tool currently. The present problem has also been modelled in COMSOL Multiphysics, and the solutions of the partial differential equations have been obtained without sacrificing the relevant physical phenomena.
The effects of slip condition on an unsteady electrically conducting MHD flow past a periodically accelerated horizontal porous plate under the influence of transverse magnetic field and the Hall current have recently been investigated [1]. The effects of optical properties of radiation and the Lorentz force on an unsteady laminar boundary MHD fluid flow and heat transfer along a semi-infinite stretching plate have been numerically investigated [2]. Unsteady squeezing MHD nanofluid flow and heat transfer in the presence of thermal radiation impact and considering the Cattaneo–Christov heat flux model, instead of the conventional Fourier’s law of heat conduction, have also been examined [3]. The effects of temperature-dependent viscosity on the MHD boundary layer flow past a vertical porous plate have been studied in the presence of a magnetic field, thermal radiation, thermophoresis through heat and mass transfer [4]. Heat transfer analysis on MHD nanofluids in a porous channel with stretching walls has been studied [5]. The laminar boundary layer flow driving past a semi-infinite flat plate has been studied taking into account the effects of a magnetic field [6].
The unsteady MHD flow of a Brinkman type nanofluid over a vertical plate embedded in a porous medium with variable surface velocity, temperature and concentration has been investigated [7]. The chemical reaction effect on a boundary layer flow of an electrically conducting micropolar fluid subject to a transverse magnetic field along a vertical plate with variable wall temperature and concentration has been reported [8]. MHD mixed convection and mass transfer flow over a semi-infinite porous inclined plate in the presence of thermal radiation with variable suction and thermophoresis has been analysed numerically [9].
The effect of chemical reactions on an unsteady MHD flow past a semi-infinite vertical porous plate with viscous dissipation has been previously analysed [10]. Analytical solutions for electrically conducting and heat generating/absorbing fluid on a continuously vertical permeable surface in the presence of radiation, a first-order homogeneous chemical reaction and mass flux have been reported [11]. The motion of an electrically conducting fluid in the presence of an external magnetic field has been investigated [12]. The heat generation/absorption and thermo-diffusion of an unsteady free convective MHD flow of a radiating and chemically reactive fluid near an infinite vertical plate through a porous medium with Hall current have been studied [13].
An electrically conducting fluid flow in an array of square ducts, separated by conducting walls of arbitrary thickness, subject to an applied magnetic field has been investigated [14]. An investigation of unsteady hydromagnetic natural convection heat and mass transfer flow past an impulsively moving infinite vertical plate embedded in a uniform porous medium when the temperature of the plate has a temporarily ramped profile has been carried out [15]. The effects of a magnetic field on the three-dimensional flow of a nanofluid having a suspension of ferrous nanoparticles within the framework of a non-uniformly thickened sheet in a slip flow regime has been analysed [16]. An analytical model of a stationary hypersonic MHD shock with an externally applied magnetic field has also been proposed [17]. Decaying homogeneous and isotropic MHD turbulence has been investigated numerically at large Reynolds numbers [18]. A model of a laminar viscous conducting flow, near a dielectric disc in a uniform magnetic field and in the presence of external rotation, has been previously considered [19]. Nanofluid flow due to a rotating disc with variable thickness has been studied with heat and mass transfer subject to volume fraction of nanoparticle and homogeneous– heterogeneous reactions [20]. The convection flow of a viscous fluid over a curved stretching sheet with heat and mass transfer has been reported [21]. Nanofluid flow on an exponentially radiating stretching sheet with heat and mass transfer has also been investigated [22]. Analysis of the boundary layer flow of a nanofluid over a stretching surface has been investigated in the presence of an applied magnetic field and a chemical reaction [23].
In the present work, two-dimensional unsteady laminar flow and MHD free convection with heat and mass transfer is investigated. The viscosity of the magnetic fluid depends on the magnetic field and the fluid is incompressible. It passes through a semi-infinite vertical porous moving plate in a porous medium. The magnetic field is subjected to a transverse direction in the presence of thermal buoyancy, concentration buoyancy and radiation effects in the fluid. The geometry of the flow is given in Fig.1.
The governing equations for the problem are as follows:
In Eq. (2) (Momentum equation), the expressions
and
From Eq. (1), it is clear that the suction velocity is normal to the plate and it is only a function of time. It can be taken in exponential form as follows:
The above equations can be made dimensionless with the help of the following dimensionless quantities:
Using the above dimensionless relations, Eq. (2)–(4) can be written as follows:
The boundary conditions for Eqs. (11)–(13) are as follows:
Eqs. (11)–(13) are modelled in COMSOL Multiphysics. Here, a 1D time-dependent model is selected under the option model wizard and then the geometry is drawn. Then, all the parameters are inserted under the option “global definition”. This modelling is done in the dimensionless mode. The above equations are modelled in COMSOL with the help of the following standard form of differential equation:
where
In Eq. (16)
The present experiment was conducted to know the physical significance of the different parameters of the flow. Figs. 23 and 4 are the unsteady results of the velocity, temperature and concentration profiles. However, Figs. 5, 6, 7, 8, 9, 10, 11, 12, 13 are the steady state results. Previously, researchers [11, 12] have used the analytical methods to solve the Eqs of the type (11), (12), (13) including both steady and unsteady nature problems. However, real physical significance of the flow should be studied through unsteady state problems of MHD flow. The numerical results are obtained here through mathematical modelling in COMSOL without sacrificing the relevant physical phenomena of MHD flow.
Fig. 2 represents the velocity profile for different values of time. In the problem, the plate is moving and the magnetic field is applied in the transverse direction. The Reynolds number is taken to be very small to keep the flow laminar. It is clear from Fig. 2 that on increasing the values of time, the velocity decreases. At
Fig. 3 indicates the temperature profiles for different values of time. Due to MHD free convection in the fluid, heat is transferred through fluid flow. For increasing values of time, the temperatures in the flow also decrease. Comparing Figs. 2 and 3, it is observed that from
Fig. 4 depicts the concentration profiles. Initially, for increasing values of time, the concentration remains approximately the same. However, it is clear from the results that the concentration increases between
Fig. 5 shows the velocity profiles for different values of the magnetic parameter. Increasing the magnetic parameter, the velocity decreases here, but the magnetic parameter does not have a very large impact since the electric field is absent. Fig. 6 represents the temperature distribution for different values of the magnetic parameter. Increasing the magnetic parameter decreases the heat transfer in the fluid. The effect of the magnetic parameter on the concentration profile can also be observed clearly in Fig. 7. Decreasing heat transfer also decreases the concentration.
Figs. 8 and 9 represent the velocity and temperature profiles, respectively. Increasing the values of the parameter
Fig. 11 indicates the effects of the chemical reaction parameter on the concentration profiles. The chemical reaction parameter γ has a negligible effect on the velocity profile and heat transfer. Therefore, it is not plotted here. The velocity and temperature distributions remain the same, as shown in Figs. 8 and 9, for different values of γ . However, increasing the values of γ , the concentration in the flow decreases.
Fig. 12 represents the temperature distributions for different values of the heat generation/absorption parameter (ϕ ) . Increasing the values of ϕ , there is no change in the velocity distribution. Therefore, it is not presented in the results. In the present problem, the parameter ϕ works as the heat absorption parameter. As the values of ϕ increases, the heat transfer decreases in the magnetic fluid. In Fig. 13, the effect of the absorption parameter on the concentration profile can be seen. If heat is absorbed during the flow, the concentration also decreases. For large values of ϕ , the concentration decreases at a linear rate.
A numerical solution is obtained through mathematical modelling. These results are very significant and reliable as compared to previously published results. The unsteady behaviour of MHD flow is also an important aspect of this study, which has been previously avoided by researchers due to its complexity. On detailed analysis, the physical parameters involved in the results might be applied in real-life applications on a priority basis. This kind of numerical approach through mathematical modelling of the problem has not been published yet.