Characterization of Marangoni Forced Convection in Casson Nanoliquid Flow with Joule Heating and Irreversibility

The Marangoni forced convective inclined magnetohydrodynamic flow is examined. Marangoni forced convection depends on the differences in surface pressure computed by magnetic field, temperature, and concentration gradient. Casson nanoliquid flow by an infinite disk is considered. Viscous dissipation, heat flux, and Joule heating are addressed in energy expressions. Thermophoresis and Brownian motion are also examined. Entropy generation is computed. The physical characteristics of entropy optimization with Arrhenius activation energy are discussed. Nonlinear PDE’s are reduced to highly nonlinear ordinary systems with appropriate transformations. A nonlinear system is numerically computed by the NDSolve technique. The salient characteristics of velocity, temperature, concentration, entropy generation, and Bejan number are explained. The computational results of the heat-transfer rate and concentration gradient are examined through tables. Velocity and temperature have reverse effects for the higher approximation of the Marangoni number. Velocity is a decreasing function of the Casson fluid parameter. Temperature is enhanced for higher radiation during reverse hold for concentration against the Marangoni number. The Bejan number and entropy generation have similar effects for Casson fluid and radiation parameters. For a higher estimation of the Brinkman number, the entropy optimization is augmented.


Introduction
An investigation of Marangoni forced convection is of great interest for the dissipative boundary layer between a two-phase liquid flow like liquid-liquid and gas-liquid interfaces. Marangoni convection depends on the difference in surface pressure determined by the concentration gradient, magnetic field, and temperature gradient. These gradients occur when the liquid boundary layer has different characteristics. Few significant applications of the Marangoni forced convection impact include thin-film stretching, material sciences, applied physics, silicon wafers, nanotechnology, semiconductor processing, soap films, etc. In addition, melting and welding are much more efficient manufacturing applications of the Marangoni convection concept. Heat transfer in the Marangoni boundary layer flow is also comprehensively deliberated. For instance, the magnetohydrodynamic Marangoni convection flow of Casson axisymmetric nanomaterials with Joule heating is presented by Shafiq et al. [1]. Rasool et al. [2] examined magnetohydrodynamic (MHD) Marangoni forced convection flow of second-grade nanomaterials with Brownian diffusion and thermophoresis. Heat flux in Marangoni convective flow of a carbon nanotube was examined by Hayat et al. [3]. Imai et al. [4] developed A nonlinear ordinary differential system is computed by the NDSolve technique. Characteristics of different involved parameters for velocity, temperature, concentration, entropy optimization, and Bejan number are examined. The computational outcome of Nusselt and Sherwood number are numerically deliberated. Reasonable agreement is found between the present and previous results in a limiting sense (see Table 1).

Statement of Problem
Marangoni forced convective inclined magnetohydrodynamic flow of Casson nanoliquid by an infinite disk is studied. Flow generated due to concentration and temperature gradients is also studied. Thermal radiation, Joule heating, and viscous dissipation in energy equation are considered. The thermal and solutal capillaries' behaviors are key components in Marangoni convection of liquid and nanoliquid. Impacts of thermophoresis and Brownian diffusion are discussed. Salient features of entropy optimization and Arrhenius activation energy are explained. The current density (J) is Ohm's laws is J = σ(E + V × B). Here, E represents the strength of the electric field, which is neglected; B is magnetic field strength, and σ is electrical conductivity. Induced magnetic field and Hall effects for low-magnetic Reynolds numbers are neglected. The ionslip effects are ignored. The magnetic field of constant strength (B 0 ) is exerted. The linearized form of thermal radiation is taken. The geometry of the problem is presented in Figure 1.  Table 1).

Statement of Problem
Marangoni forced convective inclined magnetohydrodynamic flow of Casson nanoliquid by an infinite disk is studied. Flow generated due to concentration and temperature gradients is also studied. Thermal radiation, Joule heating, and viscous dissipation in energy equation are considered. The thermal and solutal capillaries' behaviors are key components in Marangoni convection of liquid and nanoliquid. Impacts of thermophoresis and Brownian diffusion are discussed. Salient features of entropy optimization and Arrhenius activation energy are explained. The current density (J) is Ohm's laws is J = σ(E + V × B). Here, E represents the strength of the electric field, which is neglected; B is magnetic field strength, and σ is electrical conductivity. Induced magnetic field and Hall effects for low-magnetic Reynolds numbers are neglected. The ionslip effects are ignored. The magnetic field of constant strength ( 0 B ) is exerted. The linearized form of thermal radiation is taken. The geometry of the problem is presented in Figure 1. Employing continuity, momentum, energy, and concentration expressions [2,3,6]: The boundary conditions are where µ indicates the dynamic viscosity, u and w are the velocity components, ρ is density, r and z are cylindrical coordinates, T w and C w are temperature and concentration of disk, β is Casson fluid parameter, σ ⊗ is electrical conductivity, α * is inclination angle, σ is surface tension, k is thermal conductivity, T is temperature, c p is specific heat, T ∞ is ambient temperature, τ is ratio of heat capacities, k * is mean absorption coefficient, D T is thermophoresis diffusion coefficient, σ * is Stefan-Boltzman constant, D B is Brownian movement coefficient, C is concentration, k r is reaction rate, C ∞ is ambient concentration, E a is activation energy, n is fitted rate constant, K = 8.61 is Arrhenius function, and a is positive dimensional constant. We consider surface tension as a linear function of concentration and temperature [1,2]: with where γ T , σ 0 , and γ C denote the positive constants. Taking we obtain Marangoni number, γ T A Marangoni ratio parameter, Pr = ν α Prandtl number, Brownian diffusion parameter, Br(= PrEc) Brinkman number, r Ω chemical reaction parameter, E = E a KT ∞ activation energy parameter, and α 1 = T 0 T ∞ temperature difference variables.

Entropy Modeling
Entropy generation is defined as [30,31,34]: it should be noted that in the above expression, the entropy generation is because of four effects. The first term on RHS is due to the heat transfer across a finite temperature difference. It is a local volumetric entropy generation known as heat transfer irreversibility (HTI). The second term, known as fluid friction irreversibility, occurs because of viscous dissipation. The third term occurs in view of the magnetic field. The last two terms indicate diffusion irreversibility. These are because of diffusion or mass transfers across finite concentration differences.
In dimensionless variables we have

Bejan Number
The Bejan number is defined as Be = Heat and mass transfer irreversibility Total irreversibility , (15) or where indicates the entropy rate, L = RDC 0 k the diffusion variables, and α 2 = C 0 C ∞ the concentration ratio parameter.

Nusselt Number
The Nuselt number is defined as with heat flux q w given by From Equations (17) and (18) and dimensionless variables are obtained

Sherwood Number
We have with mass flux j w given by now we finally write Here, Re x = ru w ν , shows the local Reynold number.

Solution Methodology
The method NDSolve in "Numerical Differential Equations" allows us to get computational convergent solutions differential systems. The NDSolve technique deals with single differential systems and simultaneous differential systems. It also deals with partial differential systems. In the scheme of ordinary systems, there are unknown functions like Yi. These functions depend on a distinct "independent variable" x, which is the same for individual functions. It is a standard technique built in Mathematica 9.0 [37]. Table 1 is prepared for the comparison of the present results with a limiting sense and past studies. Table 1 focuses on the assessment of the Nusselt number for different estimations of the Prandtl number (Pr) with studies [35,36]. The results are very good in agreement.

Validation of Result
The above tabulated values are in reasonable agreement.

Discussion
Here, we applied the NDSolve method to develop computational outcomes for the given differential equations. Salient features of various parameters on velocity, concentration, Bejan number, entropy optimization, and temperature are studied. Velocity gradient and Nusselt and Sherwood numbers are numerically computed.

Velocity
The influence of the Casson fluid parameter (β) on f (ξ) is seen in Figure 2. Clearly, the velocity is decreased for higher (β). An increment in Casson fluid parameter yields more disturbances to fluid flow. Therefore, velocity is decreased. Physically, Casson fluid parameter (β) leads to resistance for fluid particles. Figure 3 is drawn to show the effect of (M a ) on ( f (ξ)). One can find that a higher (M a ) gives rise to velocity. Physically, a higher M a corresponds to less viscosity. A less viscous force has a tendency to speed up the fluid motion. Figure 4 reveals the outcome of (M) on f (ξ). Here, velocity decreases against higher (M). In fact, the magnetic field applied in transverse direction opposes the transport phenomena. That is why the application of such a magnetic field yields a drag-like force, namely, the Lorentz force. This force acts in the opposite direction, which results in the reduction of velocity. Such an observation qualitatively agrees with the expectations. The influence of (α * ) on f (ξ) is portrayed in Figure 5. The results of α * and M on f are similar in a qualitative sense. applied in transverse direction opposes the transport phenomena. That is why the application of such a magnetic field yields a drag-like force, namely, the Lorentz force. This force acts in the opposite direction, which results in the reduction of velocity. Such an observation qualitatively agrees with the expectations. The influence of ( α * ) on ( ) f ξ ′ is portrayed in Figure 5. The results of * and on ′ are similar in a qualitative sense.             Figure 6 examines the effect of (M a ) on θ(ξ). The temperature is increased for a higher estimation of (M a ). There is no doubt that the velocity grows in view of less viscosity for a higher (M a ). Temperature is the average kinetic energy of a fluid particle. As a result, the temperature enhances for a higher M a . The behavior of thermal radiation parameter (Rd) on (θ(ξ)) is portrayed in Figure 7. For higher radiation parameter, the temperature θ(ξ) increases. Physically, the rise of Rd has the tendency to enhance conduction impact and to increase temperature at each point away from the surface. In other words, there is a higher heat flux for a larger Rd. Figures 8 and 9 displayed the behaviors of (Nt) and (Nb) on θ(ξ). It is clear from Figures 8 and 9 that the temperature augmented for higher (Nt) and (Nb). Clearly, the velocity grows at a higher N t and N b . These parameters lead to an enhancement of thermal-boundary layer thickness. In view of such reasons, the temperature is enhanced. The behavior of (Pr) on temperature is shown in Figure 10. For higher (Pr), thermal diffusivity reduces and consequently θ(ξ) decreases. As expected, this is in accordance with the fact that when θ is at higher P r , the fluid has a thinner thermal boundary layer and it yields more temperature gradient. The velocity reduces and momentum boundary layer thickness. The impact of (M) on temperature is illustrated in Figure 11. Clearly, the fluid temperature accelerates by increasing the strength of the magnetic field. This leads to the argument that an applied magnetic field tends to heat the liquid, and thus heat transfer from the disk enhances. All of this occurs due to Lorentz force in view of the transverse applied magnetic field opposing the transport phenomena and slowing down the fluid motion. No doubt, the magnetic can be employed as a useful agent for controlling the flow and heat transfer characteristics. Figure 12 displays the influence of Brinkman number on temperature (θ(ξ)). One can find that the temperature increased against (Br). An augmentation in Br creates more kinetic energy that enhances the temperature field. Physically, Eckert number in Brinkman number is responsible for the increase in temperature. Such an increase in temperature through higher Eckert number is due to more pronounced viscous heating. Obviously, a high Eckert number gives rise to wall cooling and thus transfer of heat to fluid, so fluid temperature rises.

Temperature
displays the influence of Brinkman number on temperature . One can find that the temperature increased against ( Br ). An augmentation in Br creates more kinetic energy that enhances the temperature field. Physically, Eckert number in Brinkman number is responsible for the increase in temperature. Such an increase in temperature through higher Eckert number is due to more pronounced viscous heating. Obviously, a high Eckert number gives rise to wall cooling and thus transfer of heat to fluid, so fluid temperature rises.

Concentration
Characteristics of Marangoni number (M a ) on φ(ξ) are shown in Figure 13. For higher (M a ), the concentration decays because of surface tension produced by temperature and concentration gradients. The reduction in concentration is accompanied by simultaneous decay in the concentration of boundary layer thickness. Figures 14 and 15 display the effects of (Nb) and (Nt) on φ(ξ). Opposite impacts of (Nb) and (Nt) are seen on concentration. Such behavior totally depends upon the nanoparticles volume fraction. When Brownian motion parameter N b enhances, nanoparticles volume fraction decays and so φ reduces. On the other hand, for a higher N t , the volume fraction of nanoparticles increases. Figure 16 exhibits the effect of (E 1 ) on φ(ξ). In fact, higher activation energy parameter reduces the modified Arrhenius function, which consequently increases the generative chemical reaction. As a result, the concentration is augmented. Figure 17 shows the effect of (γ) on concentration. Clearly, for a higher (γ) the liquid becomes thicker, and therefore concentration decreases. This happens in view of higher viscosity. Figure 18 is provided to study the effect of (Le) on concentration. One can find that concentration decays versus (Le). It is noticed that high Lewis number (L e ) decays the nanoparticles volume fraction. Lewis number is the ratio of thermal diffusivity to mass diffusivity. A higher value of Lewis number enhances the thickness of the thermal boundary layer. The concentration boundary layer thickness decays. nanoparticles volume fraction decays and so φ reduces. On the other hand, for a higher , the volume fraction of nanoparticles increases. Figure 16 exhibits the effect of ( 1 E ) on ( ) φ ξ . In fact, higher activation energy parameter reduces the modified Arrhenius function, which consequently increases the generative chemical reaction. As a result, the concentration is augmented. Figure 17 shows the effect of ( γ ) on concentration. Clearly, for a higher ( γ ) the liquid becomes thicker, and therefore concentration decreases. This happens in view of higher viscosity. Figure 18 is provided to study the effect of ( Le ) on concentration. One can find that concentration decays versus ( Le ). It is noticed that high Lewis number ( ) decays the nanoparticles volume fraction. Lewis number is the ratio of thermal diffusivity to mass diffusivity. A higher value of Lewis number enhances the thickness of the thermal boundary layer. The concentration boundary layer thickness decays.

Entropy Optimization and Bejan Number
Characteristics of Brinkman ( Br ) number on G N and Be are elucidated in Figures 19 and   20. In fact, Eckert number in definition of increases the temperature, and consequently entropy enhances. Bejan number has the opposite scenario for higher . In other words, there is a direct relationship of with heat transfer by molecular diffusion and heat generated by fluid friction. For higher , there is more production of heat in the system. This raises the system disorderliness. Here, more disorderliness of the system is due to a higher diffusivity in the fluid particle. This fact leads to an increase in the diffusive variable, the entropy generation, and the Bejan number. Figures 23 and 24 are displayed to study the variation of ( Rd ) on G N and Be . For a higher estimation of ( Rd ), the G N increases. In fact, for a higher ( Rd ) the radiation enhances. As a result, there is a growth in temperature, and therefore the entropy rate increases. Be also boosts up via higher ( Rd ), because of the fact that heat transfer irreversibility dominates over total entropy optimization. Influence of ( β ) on G N and Be is explained in Figures 25 and 26. Clearly, for a higher approximation of ( β ), both G N and Be are augmented. Actually, for a higher estimation of ( β ), the more resistive force is created in flow region, which improves the thermal energy due to the collision of liquid particles. As a result, the entropy

Entropy Optimization and Bejan Number
Characteristics of Brinkman (Br) number on N G and Be are elucidated in Figures 19 and 20. In fact, Eckert number in definition of Br increases the temperature, and consequently entropy enhances. Bejan number Be has the opposite scenario for higher Br. In other words, there is a direct relationship of Br with heat transfer by molecular diffusion and heat generated by fluid friction. For higher Br, there is more production of heat in the system. This raises the system disorderliness. Figures 21 and 22 elucidate the characteristics of (L) on N G and Be. Clearly, both N G and Be are increased for a higher diffusion parameter (L). Here, more disorderliness of the system is due to a higher diffusivity in the fluid particle. This fact leads to an increase in the diffusive variable, the entropy generation, and the Bejan number. Figures 23 and 24 are displayed to study the variation of (Rd) on N G and Be. For a higher estimation of (Rd), the N G increases. In fact, for a higher (Rd) the radiation enhances. As a result, there is a growth in temperature, and therefore the entropy rate increases. Be also boosts up via higher (Rd), because of the fact that heat transfer irreversibility dominates over total entropy optimization. Influence of (β) on N G and Be is explained in Figures 25 and 26. Clearly, for a higher approximation of (β), both N G and Be are augmented. Actually, for a higher estimation of (β), the more resistive force is created in flow region, which improves the thermal energy due to the collision of liquid particles. As a result, the entropy generation is enhanced.

Engineering Quantities
Salient features of different variables on Nusselt and Sherwood numbers are discussed in Tables  1 and 2. 8.5.1. Nusselt Number

Engineering Quantities
Salient features of different variables on Nusselt and Sherwood numbers are discussed in Tables 1 and 2. 8.5.1. Nusselt Number Table 2 elucidates the numerical approximation of (Nu x ) for (M a ), (Nb), (Nt), and (Pr). Nusselt number increases via (Nb), while the reverse is observed through (Nt). For a higher (M a ), the Nusselt number is decreased. (Nu x ) increases against higher (Pr). There is a thinner thermal boundary layer for higher P r . As a result, the rate of heat diffusion enhances.   Table 3. Here, (Sh x ) boosts up for higher estimations of (Nb) and (Le). Sherwood number (Sh x ) is a decreasing function of (Nt).

Conclusions
The main points of the presented analysis are listed below.
• There is a reduction in velocity for a higher Casson fluid parameter and magnetic field. However, the scenario is different for velocity through the Marangoni number. The present attempt is basic, and modeling and simulation about such problems can be extended for the Williamson fluid subject to Soret and Dufour effects, viscous dissipation with rheological impact, and porous medium employing modified Darcy's law.
Nothing is yet known about such facts in the existing literature.