New Analytical Study of the Effects Thermo-Diffusion, Diffusion-Thermo and Chemical Reaction of Viscous Fluid on Magneto Hydrodynamics Flow in Divergent and Convergent Channels ()
1. Introduction
The importance of thermal-diffusion and diffusion-thermo effects for various fluid flows has been studied by Eckert and Drake [1] . Olajuwon [2] examined convection heat and mass transfer in a hydromagnetic flow of a second grate fluid past a semi-infinite stretching sheet in the presence of thermal diffusion and thermal radiation. Kumar et al. [3] have investigated thermal diffusion and radiation effects on unsteady magneto hydrodynamics (MHD) flow through porous medium with variable temperature and mass diffusion in the presence of heat source or sink. Magnetohydrodynamics is the study of the interaction between magnetic fields and moving, conducting fluids [4] and the behavior of an electrically conducting fluid in the presence of a magnetic field. In this case, a force is produced inside the fluid which is proportional to fluid velocity and this force always opposes the flow. Another way to produce a force inside a flowing fluid, not known widely, is the application of an externally applied magnetic as well as an externally applied electric field. This force is called Lorentz force and can be generated by a strip wise arrangement of flush mounted electrodes and permanent magnets of alternating polarity and magnetization. The Lorentz force which acts parallel to the plate can either assist or oppose the flow. The idea of using a Lorentz force to stabilize a boundary layer flow over a flat plate belongs probably to Gailitis and Lielausis [5] [6] . It is a known fact that the temperature and concentration gradients present mass and energy fluxes, respectively. Concentration gradients result in Dufuor effect (diffusion-thermo) but Soret effect (thermal-diffusion) is due to temperature gradients. The heat and mass transfer with chemical reaction plays an important role in designing of chemical processing equipment, damage of crops due to frost, formulation and dispersion of fog. The mass transfer can be defined as a phenomenon when there is an escape of vapors into the atmosphere while heat transfer happens when there is heating or cooling of a liquid or fluid. That is, both of these phenomena play an important role in the industry. Because nonlinearity of the equations for these problem exact solutions is known, so many analytical techniques have been studied. Homotopy analysis method [7] [8] and Adomian’s decomposition method [9] [10] [11] [12] are also analytical techniques used to solve the nonlinear equations. In this article the governing equations of the problem contain a system of partial differential equations which are transformed by usual transformation into a non-dimensional system of partial coupled non-linear differential equations. The purpose of this article is to investigate the diffusion-thermo and thermal diffusion effects on converging and diverging channel in the presence of chemical reaction. In present problem we apply a new technique to solve the equation governing the flow of viscous fluid in diverging and converging channels called a new algorithm. This new algorithm that includes the use of several steps, first integration and then we use Taylor expansion in addition to the last step involves extracting the value of derivatives. MHD, Soret, Dufour and chemical reaction effects are taken into account. Influences of physical parameters on temperature and concentration profiles are discussed for both diverging and converging channels with the help of graphs. The structure of this paper is organized as follows: Section 2 definitions of Mathematical formulation. Section 3 explains Description of the new scheme. Section 4, we apply new algorithm method to solve the magneto hydrodynamic (MHD) flow of viscous fluid in a channel with non-parallel plates show its ability and efficiency in finding new approximate solutions. Section 5 evidences that a new algorithm is converged through new theorems with its application. Section 6 discusses the effect of physical parameters on velocity, temperature, concentration profiles by help tables and graphics. Section 7 introduces conclusions of the present work.
2. Mathematical Formulation
Consider the flow of an incompressible fluid due to source or sink that is located at the intersection of two rigid plane walls angled
apart. Radial and symmetric nature of the flow is taken into consideration. Induced magnetic field is ignored and an applied magnetic field is considered that is applied across the flow direction. Under the aforesaid assumptions velocity field takes the form
, where
is a function of both r and
. Soret and Dufour effects are also considered that are incorporated in energy and concentration equations respectively. The fluid is also assumed to be chemically reacting. Also, the temperature and concentration are also the function of both r and
. While the results of angle opening
on concentration profile show that the increase in angle gives a decrease in concentration profile. The governing equations for mass, motion, energy, and mass transfer in polar coordinates under imposed assumptions become [13] [14] (Figure 1).
(1)
(2)
(3)
(4)
(5)
the boundary conditions are,
Figure 1. Schematic diagram of the flow problem.
(6)
where p is the fluid pressure,
is kinematic viscosity, D in order are the
specific heat and coefficient of mass diffusivity.
correspondingly are the thermal conductivity, thermal-diffusion ratio and coefficient of mass diffusivity. Further
represent the concentration susceptibility, temperature at wall, mean fluid temperature, concentration at the wall and the chemical reaction constant respectively. From the continuity Equation (1), we can write
(7)
with the use of dimensionless parameters award [6]
(8)
Eliminating p from Equations (1) and (2) using Equations (7) and (8), we get a system of nonlinear ordinary differential equation for the normalized velocity profile
, temperature profile
and concentration profile
(9)
(10)
where,
,
,
.
,
,
,
, represent Reynolds, Hartmman, Prandtl, Eckert, Dufour, Schmidt and Soret number respectively while
is the first order chemical reaction parameter. The local Nusselt and Sherwood numbers are define by
(11)
(12)
3. Description of the New Algorithm
This section describes how to obtain a new scheme to calculate the coefficients of the power series solution resulting from solving nonlinear ordinary differential equations to find analytical-approximate solution. These coefficients are important basis to construct the solution formula, therefore they can be computed recursively by differentiation ways. To illustrative the computation and operations for these coefficients and derivation the new scheme, we summarized the detail a new outlook in the following steps.
Step (1): Consider the non-linear differential equation as follows:
(13)
integrating Equation (13) with respect to
on
yield
(14)
where,
(15)
Step (2): We take Taylor series expansion of the function
about
as follows
(16)
rewriting the Equation (16)
(17)
Now, we assume that
and substituting Equation (17) in Equation (14), we obtain
(18)
where,
(19)
Step (3): We focus on computing the derivatives of G with respect to
which is the crucial part of the proposed method. Let start calculating
.
(20)
(21)
(22)
(23)
The calculations are more complicated in the second and third derivatives because of the product rules. Consequently, the systematic structure on calculation is extremely important. Fortunately, due to the assumption that the operator G and the solution f are analytic functions, then the mixed derivatives are equivalence.
We note that the derivatives function to f unknown, so we suggest the following hypothesis
(24)
Therefore Equations (20)-(23) are evaluated by
(25)
(26)
(27)
(28)
Step (4): Substituting Equations (25)-(28) in Equation (18) we get the required analytical-approximate solution for the Equation (13).
4. Application of the New Algorithm to the Magneto Hydrodynamic (MHD) Flow of Viscous Fluid in a Channel with Non-Parallel Plates
The new algorithm described in the previous section can be used as a powerful solver to the nonlinear differential Equations (9)-(10) and to find new an analytical-approximate solution. From step (1) we have
(29)
rewrite the Equation (29) as follows
(30)
where,
and
(31)
From the boundary conditions the Equation (30) becomes
(32)
From step (2) suppose that
, yield
(33)
and the analytical-approximate solution are
(34)
From step (3) yields
(35)
(36)
(37)
(38)
We note that the derivatives of f with respect
that are given in (19), can be computing by Equations (35)-(38) as
(39)
(40)
(41)
(42)
The extraction of the first derivatives of G can be represented as:
(43)
(44)
(45)
from Equation (33) by using Equations (39)-(42), gives the following,
(46)
(47)
(48)
(49)
From step (4) substitution Equations (46)-(42) in Equation (34), the analytical-approximate solution can be resulted as follows:
(50)
5. The Analysis of Convergence
Here, the analysis of convergence for the analytical-approximate solution (50) that was resulted from the application of new power series algorithm for solving the problem has been extensively studied.
Definition (1): Suppose that H is Banach space, R is the real numbers and
is a nonlinear operators defined by
. Then the sequence of the solutions generated from a new algorithm can be written as
(51)
Definition (2): 1 Suppose that
satisfies Lipschitz condition such that for
,
, we have
(52)
Now, we assume that
for simplify with
,
yield,
(53)
The sufficient condition for convergent of the series analytical-approximate solutions
,
,
is given in the following theorems.
Theorem (1): 2 The series of the analytical-approximate solution
generated from new algorithm converge if the following condition is satisfied:
(54)
Proof. From the above definition, the next equation can be written as
(55)
since G[F,H,P] satisfies Lipschitz condition. Let
, then
(56)
hence,
(57)
from Equation (57) we get
(58)
By using triangle inequality, we find that as
, we have
, then
is a Cauchy sequence in Banach space H3.
Theorem (2): 3 Let
be a nonlinear operator satisfies Lipschitz condition from H3 to H3. If the series analytical-approximate solution
converges, then it is converged to the solution of the problem (9)-(10).
Proof.
Therefore, from the Banach fixed-point theorem, there is a unique solution of the problem (9)-(10). We will prove that
converges to S.
In practice, the theorems (1) and (2) suggest to compute the value of
, as described in the following definition.
Definition (1): 4 for
(59)
Now, the definition (1) can be applied on the magneto hydrodynamic (MHD) flow of viscous fluid in a channel with non-parallel plates to find convergence, then to obtain for examples as below.
If we choose
,
,
,
,
,
,
,
then obtain:
Also, if we get
,
,
,
,
,
.
Then
,
and
converge to the solutions
,
and
respectively when
,
.
6. Results and Discussions
This section is dedicated to study the influence of various non dimensional physical parameters on velocity field
, temperature field
and concentration field
. Also the influence of different parameters on rate of heat transfer and rate of mass transfer are under observation for diverging and converging channels. In Table 1 and Table 2 proof convergence the values
,
and
of initial solutions. The stability of these values can be clearly distinguished from the fourth approximation. Table 3 and Table 4 present a comparison of the solutions obtained by new algorithm and numerical algorithm. These tables show that an excellent between the solutions. On the other hand Tables 5-13 are explained to analyze the behavior of Nusselt number and Sherwood number with variation parameters. As for can say that impotent to mention that Nusselt number gives a description of heat transfer rate at the wall, while Sherwood number represents the rate of mass transfer at the wall. In Table 5 and Table 6, Nusselt number with very Reynolds number Re and channel opening
are discussed. One can clearly observe that for diverging channel, Nusselt number a rise with increasing in Reynolds number, but for increasing Re lead to the rate of heat transfer a drop for converging channel. Sherwood number appears to be increasing as the Reynolds number grows for converging channel and decreases for diverging channel. The behavior of Nusselt number in channel opening is quite opposite for Reynolds number. Table 7 gives description of variation in heat transfer rate and mass transfer rate at the wall with the rising values of Hartmman number. Increase in heat transfer rate and mass transfer rate are observed for increasing Hartmman number Ha in the case of both converging and diverging channels. Table 8 and Table 13 portray that the rate of heat transfer increasing for both diverging and converging channels an increase in Schmidt, Soret, Prandtl, Eckert, Dufour numbers and chemical reaction parameter respectively. Furthermore these tables for both channels clear that the rate of mass transfer at the walls decreases with an increase in Schmidt, Soret, Prandtl, Eckert, Dufour numbers and chemical reaction parameter.
Table 1.
.
Table 2.
.
Table 3. Comparison between new scheme and
scheme for the analytical solutions
,
,
when
,
,
,
,
,
,
.
Table 4. Comparison between new scheme and
scheme for the analytical solutions
,
,
when
,
,
,
,
,
,
.
Table 5. Variation in Nusselt number and Sherwood number with varying
when
.
Table 6. Variation in Nusselt number and Sherwood number with varying
when
.
Table 7. Variation in Nusselt number and Sherwood number with varying
when
.
Table 8. Variation in Nusselt number and Sherwood number with varying
when
.
Table 9. Variation in Nusselt number and Sherwood number with varying
when
.
Table 10. Variation in Nusselt number and Sherwood number with varying
when
.
Table 11. Variation in Nusselt number and Sherwood number with varying
when
.
Table 12. Variation in Nusselt number and Sherwood number with varying
when
.
Table 13. Variation in Nusselt number and Sherwood number with varying
when
.
In addition to this section highlights the major outcomes of the analytical study presented by new algorithm. The analysis of the variations in temperature and concentration profiles for different parameters is prepared. For that purpose, Figures are plotted for varying several parameters. Moreover we have been divided this section into two subsections follow as.
・ Channel divergent (
).
In Figures 2-9 are plotted to show the behavior curves of velocity, temperature and concentration profiles under the impact of different physical parameters. An increasing the opening angle
gives variations in velocity, temperature and concentration profiles as displayed in Figure 2. The influence of parameter
on the velocity field
for divergent channel causes more effect at the middle channel as well as it represents as a maximum position at the central line (when
). It also has least effect in part near the walls (when
). High temperature in the central region of the channel is conspicuously clear and the maximum temperature lies there too. While the results of angle opening
on concentration profile show that the increase in
gives a decreased concentration profile. The central portion of the channel is more effect, while the portion of the near of walls the concentration is less affected. Effects of the increase Reynolds number Re in
and
are observed in Figure 3. The velocity field
is decreasing with increasing Reynold number and clearly shows that the highest level reaches the central part (when
). A rise in temperature profile
leads to an increase in Re as for the concentration profile
, it is in a state of decrease when there is an increase in Re. The prominent appearance in the central part can be observed in the lowest effect. The behavior of temperature distribution under the influence of Hartmann number Ha can be seen from Figure 4. The velocity distribution remains unchanged with the increase Hartmann. This figure gives a clear picture of how a stronger magnetic field can lead to a change in the temperature of the fluid. But Hartmann number gives a simple rise to the concentration
although this increase is simple, however can be a way to control the concentration of the fluid. In Figure 5 the change in temperature and concentration with an increase in Soret number Sr are plotted. A rise in temperature and a drop in concentration at the central portion of the channel are observed. In fact that stronger viscous forces are responsible for these phenomena. Impact of the Schmidt number Sc on temperature and concentration are demonstrated in Figure 6, with note that the effect Sr is similar to the effect of the Soret number
.Figure 7 is plotted the effect of growing values of Prandtl number on temperature profile. Arise in temperature of the fluid with increasing Pr can be seen. Figure 8 explains that change in temperature with rising Eckert number Ec. Again the increase in temperature was observed to increase Ec. The effect of Eckert number on
for diverging channel, an increase in Ec decreases the concentration profile. This shows, the strong viscous forces are responsible for a rise temperature of the fluid and this rise is a lot of effect at the central portion. In Figure 9 the changes of chemical reaction
and Dufour number
on
and
are presented. Dufour number arises due to the concentration gradient present in energy equation. The main variations in temperature are in central portion of the channel, where the variation in temperature near the walls is almost negligible, which can be said, the stronger concentration results in higher temperature values. Also this figure demonstrates that the increasing of
lead to decrease the concentration profile.
・ Channel convergent (
).
For the converging channel, the variations in velocity, temperature and concentration profile due to the varying parameters are depicted in Figures 10-17, the behavior of velocity and temperature for changing angle opening
and Reynolds number
is quite opposite to the behavior of
,
and
in diverging channel as seen in Figure 10 and Figure 11. On that other hand, Figures 11-17 tell the effect of Harmann, Eckert.
Figure 2.
,
,
for the value
,
,
,
,
,
, when the angle
is varied.
Figure 3.
,
,
for the value
,
,
,
,
,
when the Reynolds number
is varied.
Figure 4.
,
,
for the value
,
,
,
,
,
when the Harmann number
is varied.
Figure 5.
,
,
for the value
,
,
,
,
,
,
when the Soret number
is varied.
Figure 6.
,
,
for the value
,
,
,
,
,
,
when the Schmidt number
is varied.
Figure 7.
,
,
for the value
,
,
,
,
,
,
when the Prandtl number
is varied.
Figure 8.
,
for the value
,
,
,
,
,
when the Eckert number
is varied.
Figure 9.
,
for the value
,
,
,
,
,
,
when the Dufour number
and chemical reaction parameter
are varied.
Figure 10.
,
,
for the value
,
,
,
,
,
when the Reynolds number
is varied.
Figure 11.
,
,
for the value
,
,
,
,
,
, when the angle
is varied.
Figure 12.
,
,
for the value
,
,
,
,
,
when the Harmann number
is varied.
Figure 13.
,
,
for the value
,
,
,
,
,
,
when the Soret number
is varied.
Figure 14.
,
,
for the value
,
,
,
,
,
,
when the schmidt number
is varied.
Figure 15.
,
,
for the value
,
,
,
,
,
,
when the Prandtl number
is varied.
Figure 16.
,
,
for the value
,
,
,
,
,
when the Eckert number
is varied.
Figure 17.
,
,
for the value
,
,
,
,
,
,
when the Dufour number
and chemical reaction parameter
are varied.
Prandtl and Dufour numbers on the temperature profile is similar for the effect in diverging channel. Also these Figures demonstrated the concentration profile possess same effect when there are changing in Hartmann number, Schimdt number, Soret number and chemical reaction parameter. Physical explanations can be provided the temperature profile show that the temperature at the central region increases with increasing angle opening. This can be attributed to that for fixed Reynold number, increasing angle opening leads to increase the cross-sectional flow area. This in turn leads to decrease the flow velocity and this mean the flow will be decelerated. Therefore, the heat dissipation will be reduced which leads to increase the temperature of the fluid. The inertia force of the fluid increases with increasing Reynold number which leads to enhance the parabolic behavior (increasing central temperature) for diverging channel with the opposite view in converging channel. Hartmman number increase in this case Lorentz force is also increasing for diverging and converging channels. This force imports extra drag to the flow. Therefore the temperature profile becomes more flat which means decreasing the temperature within the central region. The thin boundary layers that are near to the wall lead to that the temperature gradient at the highest level. Furthermore to the existence of the thick boundary layer in central region lead to that the temperature gradient at low level. Eckert number increases with increased temperature and thus produces an increase in Kinetic energy. The change of the temperature profile with Prandtl number, and the increase of temperature with Prandtl number result from increasing of the momentum diffusivity. The Dufour number shows to increase shows less effect on temperature, increasing Dufour leads to increase the thermal energy of the fluid thus the temperature increase. The rate of most chemical reactions increases with a decrease the concentration of reactants. As for temperature, it increases if a reaction is heat-emitting and decreases when the reaction absorbs heat.
7. Conclusions
In this paper, the unsteady and two-dimensional magneto hydrodynamic (MHD) flow of viscous fluid in a channel with non-parallel plates is studied analytically using a new algorithm. The solution obtained by new algorithm is an infinite power series for appropriate initial approximation. The construction of this algorithm possessed good convergent series and the convergence of the results is explicitly shown. Graphical results and tables are presented to investigate the influence of physical parameters on velocity, temperature and concentration. Analysis of the converge confirms that the new algorithm is an efficient technique as compared to Range-Kutta algorithm with help of Shooting algorithm. The new algorithm that is widely applied to solve ordinary differential equations lead to the solutions resulting from this algorithm is compatible with numerical solution. Effects of different parameters on temperature and concentration profiles are analyzed and presented graphically. The conclusions can be drawn from the analysis presented:
・ The behavior of temperature and concentration profiles are the same results
,
,
,
and
for diverging channel.
・ Hartmann number
can be used to reduce the temperature of the flow fluid. Also, concentration of the fluid can also be controlled by employing a strong magnetic field.
・ For converging channel, the variations in temperature are opposite for diverging channel with an increase in channel opening
and
.
・ For diverging channel, Nusselt number drops with a rise in angle opening and increases with a rise in Reynolds number and behaves oppositely for convergent channel.
・ Increase in heat transfer rate is observed for increasing
,
,
,
,
and
in both channels.
・ Increase in Reynolds number and Angle opening gives a drop to mass transfer rate for diverging channel and a rise for the converging channel.
・ The rate of mass transfer decreased for both channels with an increase in Schmidt, Soret, Prandtl, Eckert, Dufour numbers and chemical reaction parameter.
・ Results obtained by new algorithm are in excellent agreement with numerical solution obtained.