100 , E c = 1 , S c = S r = 0.2 , D f = P r = 0.5 , γ = 0.1 .

Table 6. Variation in Nusselt number and Sherwood number with varying α when R e = 25 , H a = 100 , E c = 1 , S c = S r = 0.2 , D f = P r = 0.5 , γ = 0.1 .

Table 7. Variation in Nusselt number and Sherwood number with varying H a when R e = 25 , E c = 1 , S c = S r = 0.2 , D f = P r = 0.5 , γ = 0.1 .

Table 8. Variation in Nusselt number and Sherwood number with varying S c when R e = 25 , E c = 1 , S r = 0.2 , D f = P r = 0.5 , γ = 0.1 .

Table 9. Variation in Nusselt number and Sherwood number with varying S r when R e = 25 , E c = 1 , S c = 0.2 , D f = P r = 0.5 , γ = 0.1 .

Table 10. Variation in Nusselt number and Sherwood number with varying P r when R e = 25 , E c = 1 , S r = S c = 0.2 , D f = 0.5 , γ = 0.1 .

Table 11. Variation in Nusselt number and Sherwood number with varying E c when R e = 25 , S r = S c = 0.2 , D f = P r = 0.5 , γ = 0.1 .

Table 12. Variation in Nusselt number and Sherwood number with varying D f when R e = 25, E c = 1, S r = S c = 0.2, P r = 0.5, γ = 0.1 .

Table 13. Variation in Nusselt number and Sherwood number with varying γ when R e = 25, E c = 1, S r = S c = 0.2, D f = P r = 0.5 .

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 ( α > 0 ).

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 f ( η ) for divergent channel causes more effect at the middle channel as well as it represents as a maximum position at the central line (when η = 0 ). It also has least effect in part near the walls (when η = 1 , 1 ). 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 f ( η ) , β ( η ) and ϕ ( η ) are observed in Figure 3. The velocity field f ( η ) is decreasing with increasing Reynold number and clearly shows that the highest level reaches the central part (when η = 0 ). 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 S r .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 D f 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 ( α < 0 ).

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 R e is quite opposite to the behavior of f ( η ) , β ( η ) 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. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = D f = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , when the angle α is varied.

Figure 3. f ( η ) , β ( η ) , ϕ ( η ) for the value H a = 100 , P r = D f = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Reynolds number R e is varied.

Figure 4. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , P r = D f = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Harmann number H a is varied.

Figure 5. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = D f = 0.5 , E c = 1 , S c = 0.2 , γ = 0.1 , α = 3 when the Soret number S r is varied.

Figure 6. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = D f = 0.5 , E c = 1 , S r = 0.2 , γ = 0.1 , α = 3 when the Schmidt number S c is varied.

Figure 7. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , D f = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Prandtl number P r is varied.

Figure 8. β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = D f = 0.5 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Eckert number E c is varied.

Figure 9. β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Dufour number D f and chemical reaction parameter γ are varied.

Figure 10. f ( η ) , β ( η ) , ϕ ( η ) for the value H a = 100 , P r = D f = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Reynolds number R e is varied.

Figure 11. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = D f = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , when the angle α is varied.

Figure 12. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 30 , P r = D f = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Harmann number H a is varied.

Figure 13. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = D f = 0.5 , E c = 1 , S c = 0.2 , γ = 0.1 , α = 3 when the Soret number S r is varied.

Figure 14. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = D f = 0.5 , E c = 1 , S r = 0.2 , γ = 0.1 , α = 3 when the schmidt number S c is varied.

Figure 15. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , D f = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Prandtl number P r is varied.

Figure 16. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = D f = 0.5 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Eckert number E c is varied.

Figure 17. f ( η ) , β ( η ) , ϕ ( η ) for the value R e = 50 , H a = 100 , P r = 0.5 , E c = 1 , S c = S r = 0.2 , γ = 0.1 , α = 3 when the Dufour number D f 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 α , R e , E c , P r and D f for diverging channel.

・ Hartmann number H a 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 R e .

・ 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 P r , E c , S r , S c , D f 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.

Cite this paper
Al-Saif, A. and Jasim, A. (2019) 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. Applied Mathematics, 10, 268-300. doi: 10.4236/am.2019.104020.
References

[1]   Eckert, E.R.G. and Drake, R.M. (1972) Analysis of Heat and Mass Transfer. McGraw-Hill Book Company, New York.

[2]   Olajuwon, B.I. (2011) Convection Heat and Mass Transfer in a Hydromangetic Flow of a Second Grade Fluid in the Presence of Thermal Radiation and Thermal Diffusion. International Communication in Heat and Mass Transfer, 38, 377-382.
https://doi.org/10.1016/j.icheatmasstransfer.2010.11.006

[3]   Kumar, A.G.V., Goud, Y.R. and Varma, S.V.K. (2012) Thermal Diffusion and Radiation Effects on Unsteady MHD Flow through Porous Medium with Variable Temperature and Mass Diffusion in the Presence of Heat Source Sink. Advance in Applied Science Research, 3, 1494-1506.

[4]   Davidson, P.A. (2006) An Introduction to Magnetohydrodynamics. Cambridge University Press, Cambridge.

[5]   Gailitis, A. and Lielausis, O. (2006) On a Possibility to Reduce the Hydrodynamical Resistance of a Plate in an Electrolyte. Applied Magnetohydrodynamics, 12, 143-146.

[6]   Henoch, C. and Stace, J. (1995) Experimental Investigation of a Salt-Water Turbulent Boundary-Layer Modified by an Applied Stream Wise Magneto Hydrodynamic Body Force. Physics of Fluids, 7, 1371-1383.
https://doi.org/10.1063/1.868525

[7]   Ellahi, R., Raza, M. and Vafai, K. (2012) Series Solutions of Non-Newtonian Nanofluids with Reynolds’ Model and Vogel’s Model by Means of the Homotopy Analysis Method. Mathematical and. Computer Modelling, 55, 1876-1891.
https://doi.org/10.1016/j.mcm.2011.11.043

[8]   Khan, U., Ahmed, N. and Tauseef Mohyud-Din, S. (2016) Thermo-Diffusion, Diffusion-Thermo and Chemical Reaction Effects on MHD Flow of Viscous Fluid in Divergent and Convergent Channels. Chemical Engineering Science, 141, 17-27.
https://doi.org/10.1016/j.ces.2015.10.032

[9]   Hatami, M., Sheikholeslami, M., Hosseini, M. and Ganji, D.D. (2014) Analytical Investigation of MHD Nanofluid Flow in Non-Parallel Walls. Journal of Molecular Liquids, 194, 251-259.
https://doi.org/10.1016/j.molliq.2014.03.002

[10]   Hatami, M., Sheikholeslami, M. and Ganji, D.D. (2014) Nanofluid Flow and Heat Transfer in an Asymmetric Porous Channel with Expanding or Contracting Wall. Journal of Molecular Liquids, 195, 230-239.
https://doi.org/10.1016/j.molliq.2014.02.024
https://doi.org/10.1016/j.molliq.2016.10.073

[11]   Turkyilmazoglu, M. (2014) Extending the Traditional Jeffery-Hamel Flow to Stretchable Convergent/Divergent Channels. Computers and Fluids, 100, 196-203.
https://doi.org/10.1016/j.compfluid.2014.05.016

[12]   Asadullaha, A.M., Khan, U., Ahmed, N. and Tauseef Mohyud-Dinb, S. (2016) Analytical and Numerical Investigation of Thermal Radiation Effects on Flow of Viscous Incompressible Fluid with Stretchable Convergent/Divergent Channels. Journal of Molecular Liquids, 224, 768-775.
https://doi.org/10.1016/j.molliq.2016.10.073

[13]   Asadullah, M., Khan, U., Manzoor, R., Ahmed, N. and Mohyud-din, S.T. (2013) MHD Flow of a Jeffery Fluid in Converging and Diverging Channels. International Journal of Modern Mathematical Sciences, 6, 92-106.

[14]   Hayat, T., Nawaz, M. and Sajid, M. (2012) Effect of Heat Transfer on the Flow of a Second-Grade Fluid in Divergent/Convergent Channel. International Journal for Numerical Methods in Fluids, 64, 761-776.

 
 
Top