Entropy Analysis for Hydromagnetic Darcy–Forchheimer Flow Subject to Soret and Dufour Effects

: Here, our main aim is to examine the impacts of Dufour and Soret in a radiative Darcy– Forchheimer ﬂow. Ohmic heating and the dissipative features are outlined. The characteristics of the thermo-diffusion and diffusion-thermo effects are addressed. A binary chemical reaction is deliberated. To examine the thermodynamical system performance, we discuss entropy generation. A non-linear differential system is computed by the ﬁnite difference technique. Variations in the velocity, concentration, thermal ﬁeld and entropy rate for the emerging parameters are scrutinized. A decay in velocity is observed for the Forchheimer number. Higher estimation of the magnetic number has the opposite inﬂuence for the velocity and temperature. The velocity, concentration and thermal ﬁeld have a similar effect on the suction variable. The temperature against the Dufour number is augmented. A decay in the concentration is found against the Soret number. A similar trend holds for the entropy rate through the radiation and diffusion variables. An augmentation in the entropy rate is observed for the diffusion variable.


Introduction
Henry Darcy established the framework of homogeneous liquid flow through a permeable medium throughout his work on the progression of water over saturated sand [1]. At higher flow rates, when inertial and boundary impacts arise, then the Darcy law cannot work appropriately. To overcome such an issue, Forchheimer gave the concept of the non-Darcy model through the insertion of a quadratic velocity term in a momentum expression [2]. Later on, the Forchheimer term was so named by Muskat [3]. The moment of fluid flowing through a permeable surface is of keen interest due to its significance in technical, biological and scientific fields such as artificial dialysis, gas turbines, atherosclerosis, catalytic converters, geo energy production and many others. Hayat et al. [4] explored the 2-D Darcy-Forchheimer flow of non-Newtonian liquid with variable properties. Pal and Mondal [5] addressed the convective flow of Darcy-Forchheimer liquid subject to a variable heat sink or source and viscosity. Non-uniform heat conductivity analysis in the reactive flow of Darcy-Forchheimer Carreau nanomaterial subject to a magnetic dipole was presented by Mallawi and Ullah [6]. Alshomrani and Ullah [7] studied the convective flow of Darcy-Forchheimer hybrid nanomaterial subject to a cubic autocatalysis chemical reaction. Seth and Mandal [8] discussed the hydromagnetic effect in rotating the flow of Darcy-Forchheimer Casson liquid toward a permeable space. Little analysis concerning a porous medium is discussed in [9][10][11][12][13][14][15][16][17].
Thermal and solutal transportation in a permeable surface has been a significant consideration of researchers during the last two decades. This is due to its usefulness in geothermal systems, catalytic reactors, nuclear waste repositories, areas of geosciences, chemical engineering, energy storage units, drying technology, heat insulation, heat exchangers for packed beds and many others. Initially, the Dufour effect in liquid was explored by Rastogi and Madan [18]. After that, the diffusion-thermo impact in a homogeneous mixture was investigated in [19,20]. Moorthy and Senthilvalivu [21] studied the Dufour and Soret outcomes in the convective flow of liquid of a non-uniform viscosity subject to a permeable medium. Non-uniform temperature in a convective fluid flow subject to thermal-diffusion and diffusion-thermo effects was illustrated by El-Arabawy [22]. Few reviews with reference to relevant titles have been discussed through certain studies [23][24][25][26][27]. The important reason for entropy production is the conversion of thermal energy in the occurrence of numerous examples of processes, such as fluid friction, kinetic energy, rotational moments, molecular resistance, the Joule Thomson effect, mass transport rate, and molecular vibration. The concept of entropy optimization in liquid flow was initially given by Bejan [28,29]. Entropy analysis in a water-based hybrid nanoliquid subject to mixed convection was discussed by Buonomo [30]. Irreversibility exploration in the dissipative flow of nanomaterial with melting and radiation over a stretching sheet was addressed by Khan et al. [31]. Some investigations of the entropy rate are mentioned in [32][33][34][35][36][37][38][39].
To our knowledge, no study has reported about entropy-optimized radiative Darcy-Forchheimer flows with Soret and Dufour features yet. A porous medium through a Darcy-Forchheimer relation is discussed. Dissipation, Ohmic heating and radiation are scrutinized in energy equations. The physical characteristics for the Soret and Dufour impacts are addressed. A first-order chemical reaction is deliberated. To examine the thermodynamical system performance, we discuss entropy optimization. Nonlinear differential systems are obtained through appropriate transformations. Non-dimensional differential systems are solved through the finite difference method. The performance of appropriate variables concerning the velocity, entropy generation, concentration and thermal field have been scrutinized.

Formulation
An unsteady radiative hydromagnetic Darcy-Forchheimer flow saturating a porous medium is discussed. Joule heating, viscous dissipation and thermal radiation in energy expression have been scrutinized. The Soret and Dufour effects are inspected. The impact of the entropy rate is addressed. Additionally, the flow is subject to a chemical reaction of the first order. A constant magnetic field with a strength (B 0 ) is applied. Consider u = u w = ax as the stretching velocity, with a > 0. The chosen magnetic Reynolds number is small. Figure 1 shows a flow sketch [31]. By taking into account the infinite plate, the term ∂u ∂x becomes zero. Here, the continuity equation becomes ∂v ∂y = 0.
Under the above discussion, the related expression for constant suction becomes with Consider the following formula: ∂θ ∂τ − S ∂θ ∂η with Here the non-dimensional variables are and Br(= Pr Ec).

Nusselt Number
Here, we have with the heat flux q w given by We finally have

Sherwood Number
This is given as in which the mass flux j w is defined as We can write

Entropy
The important reason for entropy production is the conversion of thermal energy in the occurrence of numerous processes, such as fluid friction, kinetic energy, rotational moment, molecular resistance, the Joule-Thomson effect, mass transport rate and molecular vibration. We have the following [30][31][32][33]: The dimensionless expression is In the above expression, the dimensionless parameters are

Solution Methodology
Using the finite difference method, we can solve the nonlinear differential system [40][41][42][43] by writing By employing Equation (23) The entropy generation expression yields

Graphical Results and Review
The physical interpretation of the parameters of the concentration, entropy rate, velocity and temperature have been investigated. The present observations are compared with previous published results in Table 1, and excellent agreement is noticed.

Velocity
The influence of suction (S) upon the velocity (U(τ, η)) is sketched in Figure 2. Obviously, a higher estimation of the suction parameter (S) decays the velocity. As expected, this is in accordance to the physical facts. Figure 3 displays the velocity against a magnetic field. Actually, reduction occurs in the velocity for (M). A physically higher (M) value corresponds to amplifying the Lorentz force which the flow opposes. Hence, velocity decay is guaranteed. Figure 4 was developed in order to recognize the velocity (U(τ, η)) design with variation in the Forchheimer number (Fr). A larger estimation for the Forchheimer number decays the velocity (U(τ, η)). The influence of (λ) on (U(τ, η)) is illustrated in Figure 5. Clearly, U(τ, η) decays against higher λ values.      (θ(τ, η)) through the suction and magnetic variables. Figure 8 portrays the performance of the radiation against the temperature (θ(τ, η)). Larger radiation values lead to the temperature (θ(τ, η)) increasing. Figure 9 displays the performance of the thermal field against the Prandtl number. A larger approximation of (Pr) corresponds to the decay of the thermal diffusivity, and consequently, the temperature (θ(τ, η)) decreases. Figure 10 displays the impact of the thermal field θ(τ, η) on the Eckert number (Ec). A higher estimation for (Ec) corresponds to the temperature being higher.      Figure 11 exhibits the concentration performance against the suction variable (S). Clearly, the concentration (φ(τ, η)) was reduced against larger (S) values. Figure 12 shows the performance of the concentration (φ(τ, η)) versus (Sc). An increment (Sc) decayed the mass diffusivity, and thus the concentration (φ(τ, η)) diminished. An amplification of the Soret number (Sr) led to a decaying value for φ(τ, η) (see Figure 13). Figure 14 comprises the impact of φ(τ, η) on γ. Here, φ(τ, η) decreased against γ.

Entropy Generation Rate
The influence of the entropy rate (N G (τ, η)) via the radiation variable is disclosed in Figure 15. Clearly, a greater Rd value improved the radiation emission, which boosted the collision between the fluid particles, and so N G (τ, η) was enhanced. Figure 16 discloses the impact of L on N G (τ, η)). As predicted, the entropy generation (S G (η)) was greater via the higher approximation of L. A larger approximation of the Brinkman (Br) number enhanced the entropy generation (S G (η)) (see Figure 17). This is because of augmentation through a higher Br value causing the viscous features to improve. As a result, the entropy rate rose. Figure 18 demonstrates the entropy rate for the magnetic parameter. A larger approximation of the magnetic variable led to an increase in the entropy rate.

•
The theermal field and velocity for the magnetic field had opposing trends.