Analytical Modeling of Mhd Flow over a Permeable Rotating Disk in the Presence of Soret and Dufour Effects: Entropy Analysis

The main concern of the present article is to study steady magnetohydrodynamics (MHD) flow, heat transfer and entropy generation past a permeable rotating disk using a semi numerical/analytical method named Homotopy Analysis Method (HAM). The results of the present study are compared with numerical quadrature solutions employing a shooting technique with excellent correlation in special cases. The entropy generation equation is derived as a function of velocity, temperature and concentration gradients. Effects of flow physical parameters including magnetic interaction parameter, suction parameter, Prandtl number, Schmidt number, Soret and Dufour number on the fluid velocity, temperature and concentration distributions as well as entropy generation number are analysed and discussed in detail. Results show that increasing the Soret number or decreasing the Dufour number tends to decrease the temperature distribution while the concentration distribution is enhanced. The averaged entropy generation number increases with increasing magnetic interaction parameter, suction parameter, Prandtl number, and Schmidt number.


Introduction
Rotating disk flows have received much attention in several industrial and engineering processes.They have feasible applications in many industries, such as rotating machinery, lubrication, oceanography and computer storage devices.Von Karman [1] was, to the best of our knowledge, the first who studied fluid flow due to an infinite rotating disk.He introduced his famous appropriate transformations, resulting in ordinary differential equations, which are a reduced form of the governing partial differential equations.
In several studies, Dufour and Soret effects were assumed to be negligible on heat and mass transfer according to the effects described by Fourier's and Fick's laws [2].These effects are more important when the density differences exist in the flow regime [3].When heat and mass transfer happen simultaneously in a moving fluid, the energy flux can be generated by temperature gradients as well as composition gradients.The energy flux caused by a composition gradient is named the Dufour or diffusion-thermo effect and also the mass fluxes can be developed by the temperature gradients which are called the Soret or thermal-diffusion effect [4].In this matter, several studies have been carried out.Das et al. [5] displayed the effect of mass transfer on the free convective flow and heat transfer of a viscous incompressible electrically conducting fluid past a vertical porous plate through a porous medium.Rashidi et al. [6] presented an analytical solution for steady magnetohydrodynamics (MHD) convective and slip flow due to a rotating disk in the presence of viscous dissipation and Ohmic heating.Hayat et al. [7] illustrated thermal-diffusion and diffusion-thermo effects on 2D MHD axisymmetric flow of a second grade fluid in the presence of Joule heating and first order chemical reaction.Osalusi et al. [4] studied Soret and Dufour effects on combined heat and mass transfer of steady hydromagnetic convective and slip flow due to a rotating disk in the presence of viscous dissipation and Ohmic heating numerically using a shooting method.Pal and Talukdar [8] investigated effects of thermal radiation and first-order chemical reaction on unsteady MHD convective flow past a semi-infinite vertical flat plate in the presence of transverse magnetic field under oscillatory suction and heat source in a slip-flow regime.In another study, Turkyilmazoglu and Pop [9] studied heat and mass transfer characteristics of unsteady electrically conducting fluid flow past a suddenly started vertical infinite flat considering Soret and heat source effects.
In recent years, efficiency calculation of heat exchanger systems was restricted to the first law of thermodynamics in many studies.In many industrial systems, various mechanisms that account for irreversibility compete with each other.Thermodynamic optimization has become the concern of several researchers in recent years and is also the condition of the most desirable trade-off between two or more competing irreversibilities [10].Entropy generation minimization has been comprehensively covered by Bejan [11], specifically in the fields of refrigeration, heat transfer, storage, solar thermal power conversion, and thermal science education.Entropy generation minimization method is employed to optimize the thermal engineering devices for higher energy efficiency.In order to access the best design of thermal systems, one can employ the second law of thermodynamics by minimizing the irreversibility [12,13].The performance of engineering equipment in the presence of irreversibilities is reduced and the entropy generation function is a measure of the level of available irreversibilities in a process.Since entropy generation is a criterion for measurement of available work destruction of systems, reduction of entropy generation is essential to obtain the optimal design of energy systems [14].Moreover, entropy generation causes the systems to decrease useful power cycle outputs for a power production device or increase the power input to the cycle for power consumption devices.It is important to emphasize that the second law of thermodynamics is more reliable than the first law of thermodynamics analysis, because of the limitation of the first law efficiency in heat transfer engineering systems [15].The evaluation of entropy generation is carried out to improve the system performance.In addition, heat transfer, mass transfer, viscous dissipation, finite temperature gradients, etc. can be introduced as the sources of entropy generation [16].
Many researchers have been motivated to conduct applications of the second law of thermodynamics in the design of thermal engineering systems in recent decades.Rashidi et al. [17] investigated analysis of the second law of thermodynamics applied to an electrically conducting incompressible nanofluid flowing over a porous rotating disk.Jafari and Freidoonimehr [18] studied the second law of thermodynamics over a stretching permeable surface in the presence of a uniform vertical magnetic field in a slip nanofluid regime.In another study, Abolbashari et al. [19] used Homotopy Analysis Method (HAM) to study entropy analysis in an unsteady magnetohydrodynamic nanofluid regime adjacent to an accelerating stretching permeable surface.In addition, Abolbashari et al. [20] investigated heat and mass transfer and entropy generation for steady laminar non-Newtonian nanofluid flow induced by a stretching sheet in the presence of velocity slip and convective surface boundary conditions using Optimal HAM.Further, Ellahi et al. [21] studied natural convection boundary layer flow along an inverted cone.They also considered effect of the shape of nanosize particles on entropy generation with base fluid.
Some strongly nonlinear equations used to describe physical systems in the form of mathematical modeling have no exact solutions.Numerical or analytical methods such as Runge-Kutta Entropy 2016, 18, 131 3 of 19 method [22,23], differential transform method (DTM) [24][25][26] and the present employed analytical method, HAM, can be applied to solve these nonlinear equations.Numerical methods have many disadvantages in comparison with analytical methods.Thus in this article, we apply one of the analytical methods named Homotopy analysis method (HAM) to solve the system of nonlinear differential equations.The HAM was firstly introduced by Liao [27][28][29] to offer a general analytical method for nonlinear problems.Awad [30] presented a simple method for calculating heat transfer from a rotating disk to fluids for a wide range of Prandtl numbers using asymptotic analysis.Recently, HAM was used to solve many nonlinear problems in fluid dynamics and heat transfer [31][32][33][34].
Recently, many researchers have been interested to study MHD and its practical applications [35][36][37][38][39][40][41][42].The object of this paper is to study the second law of thermodynamics of steady MHD flow over a permeable rotating disk in the presence of Soret and Dufour effects analytically via HAM.The effects of various parameters such as magnetic interaction parameter, suction parameter, Prandtl number, Soret number, Dufour number, and Schmidt number on the fluid velocity, temperature and concentration distributions as well as the averaged entropy generation number are analyzed.

Mathematical Formulation
We assume steady, axially symmetric, incompressible flow of an electrically conducting fluid with heat and mass transfer flow past a rotating permeable disk.Consider that the fluid is infinite in extent in the positive z-direction.The fluid is assumed to be Newtonian.The external uniform magnetic field B 0 which is considered unchanged by taking the small magnetic Reynolds number is imposed in the direction normal to the surface of the disk.The induced magnetic field due to motion of the electrically-conducting fluid is negligible.Uniform suction is also applied at the surface of the disk.The flow description and geometrical coordinates are shown in Figure 1.The governing equations, respectively, of continuity, momentum, energy and species diffusion in laminar incompressible flow are given by: Bu Br `u r `Bw Bz " 0, By using cylindrical polar coordinates (r, φ, z), the disk rotates with constant angular velocity Ω and is placed at z = 0, where z is the vertical axis in the cylindrical coordinate system with r and φ as the radial and tangential axes.The components of the flow velocity (u, v, w) are in the directions of increasing (r, φ, z), respectively.P is pressure, ρ is the density of the fluid, and T and C are the fluid temperature and concentration, respectively.ν is the kinematic viscosity of the ambient fluid, σ is the electrical conductivity, k is the thermal conductivity, c p is the specific heat at constant pressure, D is the molecular diffusion coefficient, K T is the thermal diffusion ratio, C s is the concentration susceptibility, and T m is the mean fluid temperature.The appropriate boundary conditions subject to uniform suction w 0 through the disk are: susceptibility, and Tm is the mean fluid temperature.The appropriate boundary conditions subject to uniform suction w0 through the disk are: In order to obtain the non-dimensional form of Equations ( 1)-( 6), the following dimensionless variables are introduced as 1/ 2 ( ), ( ), ( ) ( ), ( ), ( ).

Entropy Generation Analysis
According to [16,[43][44][45], in the presence of axial symmetry and magnetic field and considering the mass transfer effect, the volumetric rate of local entropy generation is defined as where r∇Ts " r∇Cs " Φ " 2 We assume that the electric force per unit charge E is negligible in comparison with V ˆB, in Equations ( 17) and ( 21) and we also consider that the electric current J is much greater than QV.Thus, by applying the above assumption and after reduction, Equation ( 17) can be further simplified as follows .where Ws = w0/(ν•Ω) 1/2 is the suction/injection parameter and Ws < 0 shows uniform suction at the disk surface.

Entropy Generation Analysis
According to [16,[43][44][45], in the presence of axial symmetry and magnetic field and considering the mass transfer effect, the volumetric rate of local entropy generation is defined as ,  .
We assume that the electric force per unit charge E is negligible in comparison with 17) and ( 21) and we also consider that the electric current J is much greater than QV.Thus, by applying the above assumption and after reduction, Equation ( 17) can be further simplified as follows 2

Thermal irreversibility
Fluid friction irreversibility The right-hand side of Equation ( 22) consists of four parts; the first part is local entropy generation due to heat transfer irreversibility, the second part arises due to fluid friction irreversibility, the third part denotes the magnetic effects and entropy generation due to diffusion can be calculated from the fourth part which denotes diffusive irreversibility.It must be noted that Thermal irreversibility where Ws = w0/(ν•Ω) 1/2 is the suction/injection parameter and Ws < 0 shows uniform suction at the disk surface.

Entropy Generation Analysis
According to [16,[43][44][45], in the presence of axial symmetry and magnetic field and considering the mass transfer effect, the volumetric rate of local entropy generation is defined as ,  .
We assume that the electric force per unit charge E is negligible in comparison with 17) and ( 21) and we also consider that the electric current J is much greater than QV.Thus, by applying the above assumption and after reduction, Equation ( 17) can be further simplified as follows 2

Thermal irreversibility
Fluid friction irreversibility The right-hand side of Equation ( 22) consists of four parts; the first part is local entropy generation due to heat transfer irreversibility, the second part arises due to fluid friction irreversibility, the third part denotes the magnetic effects and entropy generation due to diffusion can be calculated from the fourth part which denotes diffusive irreversibility.It must be noted that where Ws = w0/(ν•Ω) is the suction/injection parameter and Ws < 0 shows uniform suction at the disk surface.

Entropy Generation Analysis
According to [16,[43][44][45], in the presence of axial symmetry and magnetic field and considering the mass transfer effect, the volumetric rate of local entropy generation is defined as ,  .
We assume that the electric force per unit charge E is negligible in comparison with 17) and ( 21) and we also consider that the electric current J is much greater than QV.Thus, by applying the above assumption and after reduction, Equation ( 17) can be further simplified as follows 2 Thermal irreversibility Fluid friction irreversibility The right-hand side of Equation ( 22) consists of four parts; the first part is local entropy generation due to heat transfer irreversibility, the second part arises due to fluid friction irreversibility, the third part denotes the magnetic effects and entropy generation due to diffusion can be calculated from the fourth part which denotes diffusive irreversibility.It must be noted that Magnetic effects where Ws = w0/(ν•Ω) 1/2 is the suction/injection parameter and Ws < 0 shows uniform suction at the disk surface.

Entropy Generation Analysis
According to [16,[43][44][45], in the presence of axial symmetry and magnetic field and considering the mass transfer effect, the volumetric rate of local entropy generation is defined as ,  .
We assume that the electric force per unit charge E is negligible in comparison with 17) and ( 21) and we also consider that the electric current J is much greater than QV.Thus, by applying the above assumption and after reduction, Equation ( 17) can be further simplified as follows 2 Thermal irreversibility Fluid friction irreversibility The right-hand side of Equation ( 22) consists of four parts; the first part is local entropy generation due to heat transfer irreversibility, the second part arises due to fluid friction irreversibility, the third part denotes the magnetic effects and entropy generation due to diffusion can be calculated from the fourth part which denotes diffusive irreversibility.It must be noted that Diffusive irreversibility .(22) The right-hand side of Equation ( 22) consists of four parts; the first part is local entropy generation due to heat transfer irreversibility, the second part arises due to fluid friction irreversibility, the third part denotes the magnetic effects and entropy generation due to diffusion can be calculated from the fourth part which denotes diffusive irreversibility.It must be noted that the entropy generation due to diffusion is the sum of a cross term with both thermal and concentration gradients and a pure term which involves concentration gradient only.The non-dimensional form of the entropy generation rate is the entropy generation number and indicates the ratio of the actual entropy generation rate ´. S 3 gen to the characteristic entropy generation rate ´. S 3 0 ¯.By using the Von Karman similarity variables given in Equations ( 9) and (10), the entropy generation number pN G q becomes: where α " ∆T{T w is the dimensionless temperature difference, Br " µ Ω 2 L 2 {k∆ is the rotational Brinkman number, Re " ΩL 2 {ν is the rotational Reynolds number, β " ∆C{C w is the dimensionless concentration difference, λ " R g D∆C{k is the diffusive constant parameter, and gen { pkΩ∆T{ν T w q is the dimensionless entropy generation rate.As the Brinkman number (a measure of the effect of viscous dissipation) increases, entropy generation due to viscous and magnetic effects also increases.Increasing the Reynolds number results in a decrease in the effect of fluid friction irreversibility.
The averaged entropy generation number, which is an important measure of total global entropy generation can be evaluated using: where @ is the volume considered.In order to consider the effects of velocity, thermal and concentration boundary-layers, we calculate the volumetric entropy generation in a large finite domain.Thus, the integration of Equation ( 24) is obtained in the domain 0 ď R ď 1 and 0 ď η ď m; where m is a sufficiently large number.

Optimal Convergence Control Parameters
Convergence of the series in Equations ( 40)-( 44) forcefully depends on the auxiliary parameter, as mentioned by Liao [27].It is essential to select an appropriate value of the auxiliary parameter to control and speed up convergence of the approximation series with the assistance of the so-called -curve.Obviously, the valid regions of correspond to the line segments nearly parallel to the horizontal axis.The -curves of F 1 p0q, G 1 p0q, H 2 p0q, θ 1 p0q and ϕ 1 p0q obtained by the 20th order approximation are shown in Figure 2.
Convergence of the series in Equations ( 40)-( 44) forcefully depends on the auxiliary parameter, as mentioned by Liao [27].It is essential to select an appropriate value of the auxiliary parameter to control and speed up convergence of the approximation series with the assistance of the so-called  -curve.Obviously, the valid regions of  correspond to the line segments nearly parallel to the horizontal axis.The -curves of (0), (0), (0), (0) obtained by the 20th order approximation are shown in Figure 2. The averaged residual errors are defined as: For example, using the above residual error equations (Equations ( 64) and (65)), one can obtain the optimum values of auxiliary parameters using the below equations [46][47][48]: where ∆x " 10{K and K " 20.For the given order of approximation m, the optimal value of is given by the minimum values of the ∆ F, m and ∆ G, m corresponding to nonlinear algebraic equations: In order to check the accuracy of the method and determine the optimal values of , the residual errors displayed in Equations ( 64) and (65) for the 20th order approximation HAM solutions are presented in Figure 3.
In order to check the accuracy of the method and determine the optimal values of  , the residual errors displayed in Equations ( 64) and (65) for the 20th order approximation HAM solutions are presented in Figure 3.  64) and (b) Equation ( 65) obtained by 20th order approximation of the HAM solution

Results and Discussion
The nonlinear ordinary differential Equations ( 11)-( 15) subject to boundary conditions ( 16) have been solved via HAM for some values of the magnetic interaction parameter M, Prandtl number Pr, Schmidt number Sc, Soret number Sr, Dufour number Du, and suction parameter W s .The values of the physical flow parameters are mentioned in each of the graphs and tables.Tables 1 and 2 illustrate a comparison between our results and those reported by Turkyilmazoglu [49] and Kelson and Desseaux [50] for F 1 p0q and G 1 p0q as well as different values of magnetic interaction parameter and suction parameter.Excellent agreement is observed.The diluting chemical species of most common interest have Schmidt number between 0.1 and 10.0.Thus, we chose Schmidt number 0.22, 0.64, 0.78, and 1, which represent the Schmidt number of helium, ammonia, carbon monoxide, and carbon dioxide, respectively.The effect of magnetic interaction parameter on the velocity components in radial, tangential and axial directions, temperature distribution as well as concentration profile is presented in Figure 4. Infliction of the vertical magnetic field to the electrically conducting fluid causes a drag-like force named the Lorentz force.This force has the tendency to slow down the flow around the disk at the expense of increasing its temperature and concentration.Thus, as the magnetic field becomes stronger, the velocity profiles in radial, tangential and axial directions decrease and the thermal boundary layer and concentration field increase.It is important to note that as the vertical magnetic field increases, increased resistance on the fluid particles applies, which causes heat to generate in the fluid.suction parameter.
Figure 6 demonstrates the effect of Prandtl number on the temperature distribution as well as the effect of Schmidt number on the concentration profile.The thermal boundary-layer thickness decreases with increasing Prandtl number.This physically means that the flow with large Prandtl number prevents spreading of heat in the fluid.As the Schmidt number increases, the concentration boundary layer thickness decreases.In other words, molecular diffusion decreases, as Schmidt number increases.
Figure 7 indicates the simultaneous effects of Soret and Dufour number on temperature distribution as well as concentration profile.The Soret effect is a mass flux due to a temperature gradient and the Dufour effect is enthalpy flux due to a concentration gradient and appears in the energy equation.It should be mentioned that Dufour and Soret numbers are arbitrary constants provided that their product remains constant [6,7,51].Moreover, Du = 0 and Sr = 0 correspond to the condition when thermal diffusion and diffusion thermo effects are of smaller order of magnitude than effects described by Fourier's and Fick's laws [7,52].The thermal boundary layer increases by increasing Dufour number or simultaneously decreasing Soret number.As the Dufour number increases or Soret number decreases, the rate of mass transfer (concentration boundary layer thickness) decreases at the disk.(  H ()   The effect of the suction parameter on all velocity components as well as temperature distribution and concentration profile is illustrated in Figure 5.When suction is applied at the disk surface, the radial, tangential and axial velocity profiles decrease.Applying suction leads to draw the amount of fluid particles into the wall and consequently the velocity boundary layers decrease.In addition, the radial velocity component becomes very small for large values of the suction parameter.The usual decay of temperature and concentration profiles occurs for larger values of the suction parameter.
-0.5 0  Figure 6 demonstrates the effect of Prandtl number on the temperature distribution as well as the effect of Schmidt number on the concentration profile.The thermal boundary-layer thickness decreases with increasing Prandtl number.This physically means that the flow with large Prandtl number prevents spreading of heat in the fluid.As the Schmidt number increases, the concentration boundary layer thickness decreases.In other words, molecular diffusion decreases, as Schmidt number increases.Figure 7 indicates the simultaneous effects of Soret and Dufour number on temperature distribution as well as concentration profile.The Soret effect is a mass flux due to a temperature gradient and the Dufour effect is enthalpy flux due to a concentration gradient and appears in the energy equation.It should be mentioned that Dufour and Soret numbers are arbitrary constants provided that their product remains constant [6,7,51].Moreover, Du = 0 and Sr = 0 correspond to the condition when thermal diffusion and diffusion thermo effects are of smaller order of magnitude than effects described by Fourier's and Fick's laws [7,52] Figures 8-11 display the averaged entropy generation number as a function of suction parameter, Prandtl number, Schmidt number, and Soret and Dufour numbers for a wide range of magnetic interaction parameter.It should be mentioned that the values of α, β and λ are chosen as unity.Figure 8 shows that as the suction through the disk surface increases, the averaged entropy generation number also increases.Moreover, increasing the Prandtl number increases N G, av (Figure 9).It is clear that Schmidt number follows the same trend as Prandtl number.From Figure 10 maximum values of the averaged entropy generation number occur when the values of both Soret and Dufour numbers are maximized simultaneously.Finally, all entropy generation related figures reveal that as the magnetic interaction parameter increases, the averaged entropy generation number also increases.
unity. Figure 8 shows that as the suction through the disk surface increases, the averaged entropy generation number also increases.Moreover, increasing the Prandtl number increases NG, av (Figure 9).It is clear that Schmidt number follows the same trend as Prandtl number.From Figure 10 maximum values of the averaged entropy generation number occur when the values of both Soret and Dufour numbers are maximized simultaneously.Finally, all entropy generation related figures reveal that as the magnetic interaction parameter increases, the averaged entropy generation number also increases.unity.Figure 8 shows that as the suction through the disk surface increases, the averaged entropy generation number also increases.Moreover, increasing the Prandtl number increases NG, av (Figure 9).It is clear that Schmidt number follows the same trend as Prandtl number.From Figure 10 maximum values of the averaged entropy generation number occur when the values of both Soret and Dufour numbers are maximized simultaneously.Finally, all entropy generation related figures reveal that as the magnetic interaction parameter increases, the averaged entropy generation number also increases.

Conclusions
In the current study, a mathematical formulation has been derived for the second law of thermodynamics analysis of MHD fluid flow due to a permeable rotating disk.HAM was used to solve the system of ordinary differential equations.The presented semi numerical/analytical simulations agree closely with previous studies for some special cases.HAM has been shown to be a very strong and efficient technique in determining analytical solutions for nonlinear differential equations.Effects of the six key thermo-physical parameters governing the flow, i.e., magnetic interaction parameter, Prandtl number, Schmidt number, Soret number, Dufour number, and suction parameter on the all dimensionless velocity components, temperature distribution and concentration profile as well as the averaged entropy generation number have been depicted graphically and interpreted in detail.The main results of the present analysis are listed below:

Conclusions
In the current study, a mathematical formulation has been derived for the second law of thermodynamics analysis of MHD fluid flow due to a permeable rotating disk.HAM was used to solve the system of ordinary differential equations.The presented semi numerical/analytical simulations agree closely with previous studies for some special cases.HAM has been shown to be a very strong and efficient technique in determining analytical solutions for nonlinear differential equations.Effects of the six key thermo-physical parameters governing the flow, i.e., magnetic interaction parameter, Prandtl number, Schmidt number, Soret number, Dufour number, and suction parameter on the all dimensionless velocity components, temperature distribution and concentration profile as well as the averaged entropy generation number have been depicted graphically and interpreted in detail.The main results of the present analysis are listed below:

Figure 1 .
Figure 1.Configuration of the flow and geometrical coordinates.

Figure 1 .
Figure 1.Configuration of the flow and geometrical coordinates.

Figure 2 .Figure 2 .
Figure 2. The  -curves obtained by 20th order approximation of the HAM solution.h

Figure 3 . 40 Figure 3 .
Figure 3.The residual error of (a) Equation (64) and (b) Equation (65) obtained by 20th order approximation of the HAM solution

Figure 5 .Figure 6 .Figure 7 .Figure 7 .
Figure7indicates the simultaneous effects of Soret and Dufour number on temperature distribution as well as concentration profile.The Soret effect is a mass flux due to a temperature gradient and the Dufour effect is enthalpy flux due to a concentration gradient and appears in the energy equation.It should be mentioned that Dufour and Soret numbers are arbitrary constants provided that their product remains constant[6,7,51].Moreover, Du = 0 and Sr = 0 correspond to the condition when thermal diffusion and diffusion thermo effects are of smaller order of magnitude than effects described by Fourier's and Fick's laws[7,52].The thermal boundary layer increases by increasing Dufour number or simultaneously decreasing Soret number.As the Dufour number increases or Soret number decreases, the rate of mass transfer (concentration boundary layer thickness) decreases at the disk.

DuFigure 11 .
Figure 11.Change of N G,av with respect to magnetic interaction parameter for different values of Soret and Dufour numbers when Pr " 0.71, Sc " 0.78, W s " ´1 and Re " Br " 5.
(a) HAM is shown to demonstrate excellent potential, convergence and accuracy for simulating flow over rotating disk problems.(b) As the magnetic field becomes stronger, the velocity profiles in radial, tangential and axial directions decrease and the thermal boundary layer and concentration field increase.(c) When suction is applied at the disk surface, the radial, tangential and axial velocity profiles decrease.The usual decay of temperature and concentration profiles occurs for larger values of the suction parameter.(d) The thermal boundary-layer thickness decreases with increasing Prandtl number.Furthermore, as the Schmidt number increases, the concentration boundary layer thickness decreases.(e) The thermal boundary layer increases by increasing Dufour number or simultaneously decreasing Soret number.As the Dufour number increases or Soret number decreases, the rate of mass transfer (concentration boundary layer thickness) decreases at the disk.(f) The averaged entropy generation number increases by increasing the magnetic interaction parameter, suction parameter, Prandtl number, and Schmidt number.In addition, the maximum values of averaged entropy generation number occur when the values of both Soret and Dufour numbers are maximized simultaneously.
is the Schmidt number, Sr = D¨(T 8 ´Tw )¨K T /ν¨T m (C 8 ´Cw ) is the Soret number, Du = D¨(C 8 ´Cw )¨K T /C s ¨cp ¨ν¨(T 8 ´Tw ) is the Dufour number, and F, G, H, θ, and ϕ are the dimensionless functions of modified dimensionless vertical coordinate η.

Table 1 .
Numerical values of the radial skin friction coefficient F 1 p0q.

Table 2 .
Numerical values of the tangential skin friction coefficient ´G1 p0q.
Change of N G,av with respect to magnetic interaction parameter for different values of suction parameter when Du " 0.2, Sr " 0.25, Pr " 0.71, Sc " 1, and Re " Br " 5.