Impulsively Started Horizontal Riga Plate Embedded in Unsteady Casson Fluid Flow with Rotation ()
1. Introduction
Casson fluid flows along the Riga plate create a new research area that plays a fundamental role in various industrial and engineering processes, such as MHD generators, thermal nuclear reactors, flow meters and the design of nuclear reactors. Such flows have tremendous applications in civil engineering, mechanical engineering, chemical engineering, food processing and biomechanics. Riga plate is generated by the combination of electrodes and permanent magnets that create a plane surface instead of polarity and magnetization. The flow of Casson fluids with heat transfer is widely used in the processing of chocolate, foams, syrups, toffee and food processing industries. The Riga plate is first induced by Gailitis and Leilausis [1] to generate a wall paralleled Lorentz force to control the fluid flow. Anjum et al. [2] explained the thermally stratified viscous fluid with stagnation point flow dominated by a variable thicked non-linear Riga plate. Ahmed et al. [3] studied the united convection boundary layer flows of a nanofluid along a Riga plate stated vertically with strong suction. The study of a nanofluid flow past a vertical convective heated Riga plate is analyzed by Wahidunnisa et al. [4]. The characteristics of nanofluid boundary layer flow occupied with a Riga plate is concerned by Hayat et al. [5]. Pantokratoras and Magyari [6] proposed an electro-magneto hydrodynamic free convection flow of a weakly conducting fluid from an electromagnetic actuator or Riga plate. Thereafter Pantokratoras [7] investigated the Riga-plate moves with constant velocity or the Riga-plate is situated in a constant free stream. Iqbal et al. [8] investigated the boundary layer Casson fluid flow for a stagnation point along a Riga plate. The incompressible, steady and electrically conducting non-Newtonian Casson fluid flow among two cylinders, where a magnetic field is applied radially is studied by Eldabe et al. [9]. Hydrodynamic impulsively lid-driven flow with heat transfer of a Casson fluid has been described by Attia and Sayed [10]. Hayat et al. [11] analyzed the MHD Casson fluid flow of the Soret and Dufour effects upon a stretched surface. Ramesh and Devakar [12] have considered three fundamental flows Couette, Poiseuille and generalized Couette flow of an incompressible Casson fluid between parallel plates with slip boundary conditions. Mahantesh et al. [13] have considered a two-dimensional flow of non-Newtonian MHD flow of Casson fluid. Raju et al. [14] analyzed the heat source and chemical reaction effect on Casson fluid flow past an exponentially enterable stretching surface with thermal radiation and viscous dissipation. Rao and Sreenadh [15] investigated the steady two-dimensional MHD convective boundary layer Casson fluid flow over an exponentially inclined permeable stretching surface with thermal radiation and chemical reaction. Saidulu and Lakshmi [16] presented the heat and mass transfer boundary layer fluid flow of a non-Newtonian Casson fluid towards a porous stretching sheet with thermal radiation and suction. The Casson fluid flow with heat transfer and thermal radiation along a porous stretching surface is presented by Pramanik [17]. Ahmad [18] described the effect of the Powell-Eyring and Reiner-Phillipoff fluid flow on the Riga plate. Yucel [19] explained the associated free convection flow in a micropolar fluid through a Riga plate. Ghulam Rasool et al. [20] investigated the effect of a chemical reaction of nanofluids in the presence of Lorentz force and thermal radiation.
From the above-mentioned discussion of the authors, we have come to a completely new opinion that has not discussed yet before. The royal object is to investigate the behavior of the fluid particles which are adjacent with an impulsively started horizontal Riga plate in Casson fluid flow where appearance a rotation. Numerical investigations are gained by the finite difference method and its solutions are computed by using the code of Compaq visual FORTRAN 6.6a and also MATLAB R2015a.
2. Problem Formulation
Consider the horizontal Riga plate in laminar, incompressible Casson fluid flows, which is parallel to the x-axis, and the y-axis is normal to it (Figure 1). Also, consider the system rotated about a vertical axis with angular velocity
. Initially, both the plate and fluid are stagnant. At
the plate is set in motion at a constant velocity
in the direction of x-axis. This generates a two dimensional parallel flows near the plate. Instantaneously at time
the temperature of the plate and concentration are raised from
to
and
to
respectively, which are thereafter sustained. Here
,
are temperature and concentration at the wall and
,
designate the temperature and concentration outside the boundary layer. Due to the Riga plate, the Lorentz force is defined as
, which generates magnetic force. According to the Grinberg investigation, this magnetic force is defined as follows:
where,
is the applied current density,
is the magnetization of the permanent magnets and l is the width of magnets and electrodes.
The fundamental equation of Casson fluid can be written as
where
is the shear stress,
is known as plastic viscosity of the non-Newtonian
fluid,
is known as yield stress of the fluid, mathematically expressed as
, where
is the parameter of the Casson fluid,
, (i.e.
is the product of the component of deformation rate with itself, here
is the
component of deformation rate),
is the critical value based on the non-Newtonian model. Dynamic viscosity of Casson fluid is defined as
so that kinematic viscosity may define as
where
.
Within the framework of the above assumptions, the equations related to the unsteady two-dimensional problems governed by the following system of coupled non-linear partial differential equations under the boundary layer approximations, are given as follows:
(1)
(2)
(3)
(4)
(5)
And the corresponding boundary conditions for the problem are given as follows:
(6)
Here, x and y are the Cartesian coordinates systems, u and w are the velocity in the x and y direction respectively,
is the specific heat at the constant pressure, Q is the variable heat source parameter,
is the density of the fluid, k is thermal conductivity,
is the kinematic viscosity,
is the mass diffusivity,
is the coefficient of the thermal diffusivity and
is the Casson fluid parameter.
Similarity Analysis
Now introducing the non-dimensional variables and the necessary transformations are as follows:
Using these quantities into the Equations (1)-(6), removing asterisk sign it is found the dimensionless form of the equations are as follows:
With the boundary conditions:
where,
is the Modified Hartmann number;
is the Rotational parameter;
is the Prandtl number ;
is the Radiative parameter;
is the Eckert number;
is the Heat source parameter;
is the Schmidt number ;
is the Soret number.
3. Method of Solution
The governing non-linear coupled dimensionless form of partial differential Equations (6) to (10) with the associated initial and boundary conditions (11) are solved. In fact for this set of equations, it is not possible to find its exact solution and hence it has been calculated these equations by using the explicit finite difference method. The finite difference schemes with respect to t, x and y are as follows:
Here, the subscripts i and j refer to x and y and the superscript k refers to time t. Same way it has been followed for the other variables. The graphical solutions of the program are illustrated by using both the Compaq visual FORTRAN 6.6a and MATLAB R2015a tools.
Skin Friction, Nusselt Number and Sherwood Number
The effects of various parameters on local and average shear stress from the velocity profile have been investigated. The velocity gradient at the plate is defined as the shear stress; the non-dimensional form of the local shear stress and average shear stress in x-direction is given by the relations
and
respectively. From the temperature profile, the effects of various parameters on local and average Nusselt number have been calculated. The rate of heat transfer at the plate is defined as the Nusselt number; the local Nusselt number and the average Nusselt number are given by
and
respectively. Similarly, the effects of various parameters on local and average Sherwood number from the concentration field have been analyzed. The rate of mass transfer at the plate is defined as the Sherwood number; the local Sherwood number and the average Sherwood number is defined by
and
respectively.
4. Results and Discussion
The behavior of the primary velocity (u), the secondary velocity (w), temperature (
) and concentration distributions (
) have been analyzed for the different values of Rotational parameter (R), modified Hartmann number (
), Prandtl number (
), Radiative parameter (
), Eckert number (
), Heat source parameter (
), Casson fluid parameter (
), Schmidt number (
) and the Soret number (
). Consider the values of
,
,
,
,
,
,
,
and
are fixed as common for all cases except the varied values of respective parameter. The flow characteristics are shown graphically from Figures 2-19.
4.1. Mesh and Time Sensitivity Test
To find the steady-state solution of our problem, the computations have been carried out for different mesh pairs (Figure 2(a)). It has been shown that there are negligible changes for the mesh pairs
,
and
. Also Figure 2(b) depicts that the velocity distributions for the different time
and 8 with the time increment
with mesh pair
. It is observed that after
, there is very minimal change in the velocity distribution. Hence from both figures that there is a suitable choice of maximum time
and mesh pair
with time increment
for steady-state solution.
4.2. Code Sensitivity Test
To verify the results, the computations data are collected by using Compaq visual FORTRAN 6.6a and MATLAB R2015a tools on the velocity distributions with the increasing values of the Casson parameter (
) and its graphical representations are shown in Figure 3(a) and Figure 3(b) respectively. Here it is found that the velocity decreases with the increases of
. The same results have found in both codes.
Figure 2. (a) Mesh sensitivity on primary velocity u; (b) Time sensitivity on primary velocity u.
Figure 3. (a) Primary velocity distribution for u for different values of
[Fortran code]; (b) Primary velocity distribution u for different values of
[Matlab code].
4.3. Primary Velocity Distribution
It is seen from Figure 4 that the nature of fluid velocity for the variation of the Hartmann number, there are very negligible changes in the primary velocity distributions with the increase of
, which is shown in enlarge figure in Figure 4. Figure 5 represents the influence of the rotational parameter (R) on the primary velocity is decreased with the increase of R.
4.4. Secondary Velocity Distribution
It is noticed from Figure 6, with the increasing values of the Casson fluid parameter (
), the secondary velocity w has minor increasing flow near the plate and thereafter it has decreasing flow. It is observed from Figure 7 that w has been very minor increasing effects with the increasing values of Ha. Figure 8 depicts that w is increased with the increase of R.
Figure 4. Primary velocity distribution u for different values of Ha.
Figure 5. Primary velocity distribution u for different values of R.
Figure 6. Secondary velocity distribution w for different values of
.
Figure 7. Secondary velocity distribution w for different values of Ha.
Figure 8. Secondary velocity distribution w for different values of R.
4.5. Temperature Distribution
Figure 9 is displayed that the rising values of the Eckert number (
) lead to raise the fluid temperature. For the influence of
the temperature
is decreased, it is depicted in Figure 10. It is evident from Figure 11 that the effect of
, the temperature
is decreased.
4.6. Concentration Distribution
It is found from Figure 12 that near the plate the concentration distributions has an increasing effect with the increase of Qs, thereafter it has very minor reverse effects. But Figure 13 demonstrates the effect of Sorret number on the concentration profile is increased.
4.7. Local and Average Skin Friction
It is mentioned that the figures (a) and (b) of Figures 14-16 refer to the local
Figure 9. Temperature distribution
for different values of Pr.
Figure 10. Temperature distribution
for different values of Qs.
Figure 11. Temperature distribution
for different values of Qr.
Figure 12. Concentration distribution
for different values of Qs.
Figure 13. Concentration distribution
for different values of S0.
Figure 14. (a) Local shear stress
against R; (b) Average shear stress
against R.
Figure 15. (a) Local shear stress
against R; (b) Average shear stress
against R.
Figure 16. (a) Local shear stress
against
; (b) Average shear stress
against
.
shear stress [
(or
)] and average shear stress [
(or
)] respectively. It is observed from Figure 14(a) & Figure 14(b) that the local (and average) primary shear stresses both are decreases with the increase of R. But for secondary velocity distribution, local and average shear stresses both are increases with R, which is shown in Figure 15(a) & Figure 15(b). The effects of Casson parameter on the local and average primary shear stresses have shown in Figure 16(a) & Figure 16(b). Both figures have shown a very minor decreasing effect with the increases of
.
4.8. Local and Average Nusselt Number and Sherwood Number
Here Figure 17(a) & Figure 17(b) refer to the influence of the heat source parameter on the Local (and average) Nusselt number. It has sheared that
and
both increase with the increase of
. But the figures (a) and (b) of Figure 18 & Figure 19 depict the local (and average) Sherwood number
Figure 17. (a) Local Nusselt number
against
; (b) Average Nusselt number
against
.
Figure 18. (a) Local Sherwood number
against
; (b) Average Sherwood number
against
.
Figure 19. (a) Local Sherwood
number against
; (b) Average Sherwood number
against
.
(and
) respectively. It is noticed in all figures, (
and
) both decrease with the increase of
and
respectively.
5. Conclusions
The impulsively started horizontal Riga plate with rotation of unsteady Casson fluid flow has been discussed. The influence of various values of the Rotational parameter, modified Hartmann number, Prandtl number, Radiative parameter, Eckert number, Heat source parameter, Schmidt number and the Soret number has been analyzed. The non-linear coupled governing equations are solved numerically and the main findings can be summarized as follows:
· The primary velocity u increases with the increase of
, while it decreases with the increase of
and R.
· The secondary velocity w rises with the increase of R and
, while it decreases with the increase of
.
· The fluid temperature
is increased with the increase of
but it decreases with the increase of
and
.
· The concentration
rises for the increasing values of
and
.
·
(or
) is decreased with increase of R and
, while
(or
) is increased with increase of R and
. On the other hand
(or
) is increased with increase of
, while local (or average) Sherwood number is decreased with the influence of
and
.