Literature DB >> 33267198

Numerical Simulation of Entropy Generation for Power-Law Liquid Flow over a Permeable Exponential Stretched Surface with Variable Heat Source and Heat Flux.

Mohamed Abd El-Aziz^1,2, Salman Saleem¹.

Abstract

This novel work explored the second law analysis and heat transfer in a magneto non-Newtonian power-law fluid model with the presence of an internal non-uniform heat source/sink. In this investigation, the motion of the studied fluid was induced by an exponentially stretching surface. The rheological behavior of the fluid model, including the shear thinning and shear thickening properties, are also considered as special case studies. The physical problem developed meaningfully with the imposed heat flux and the porosity of the stretched surface. Extensive numerical simulations were carried out for the present boundary layer flow, in order to study the influence of each control parameter on the boundary layer flow and heat transfer characteristics via various tabular and graphical illustrations. By employing the Shooting Runge-Kutta-Fehlberg Method (SRKFM), the resulting nonlinear ordinary differential equations were solved accurately. Based on this numerical procedure, the velocity and temperature fields are displayed graphically. By applying the second law of thermodynamics, and characterizing the entropy generation and Bejan number, the present physical problem was examined and discussed thoroughly in different situations. The attained results showed that the entropy generation can be improved significantly by raising the magnetic field strength and the group parameter. From an energetic point of view, it was found that the Reynolds number boosts the entropy generation of the fluidic medium and reduces the Bejan number. Also, it was observed that an amplification of the power-law index diminished the entropy generation near the stretched surface. As main results, it was proven that the heat transfer rate can be reduced with both the internal heat source intensity and the magnetic field strength.

Entities: Chemical Disease Gene

Keywords: entropy generation; non-uniform heat source; numerical scheme; permeable surface; power-law fluid

Year: 2019 PMID： 33267198 PMCID： PMC7514973 DOI： 10.3390/e21050484

Source DB: PubMed Journal: Entropy (Basel) ISSN： 1099-4300 Impact factor: 2.524

1. Introduction

The collective necessity of energy with controlled sources has encouraged researchers and engineers to review the devices used for energy exchange and yield advanced techniques for active consumption of partial energy means. This energy requirement is fulfilled by the second analysis of thermodynamics. Entropy generation relates to the measure of eradication of productive energy. The heat transfer phenomenon initiated by its three main sources is the most feasible reason for such energy loss. Additionally, buoyancy and magnetic fields also contribute to entropy analysis. The minimization of entropy generation (EMG) in industrial procedures is among the novel concepts to augment the production of thermal systems. Several sources, such as heat transfer and viscous dissipation, are engaged for the formation of entropy. Additionally, entropy production establishes an augmentation in the power input in the power cycle for power consumption processes, whereas it contracts for power cycle outputs in a power-making system. Production in various systems can be upgraded with entropy generation and the Bejan number. Primarily, the first law of thermodynamics was employed for thermal efficiency in a system. Later, researchers testified that the second law investigation is as precise and significant compared to the first law. Bejan [1] was the primary investigator who adapted the theory of entropy generation. He perceived that conductive and viscous irreversibilities are the two key sources of entropy generation. Analysis of heat and entropy for squeezing nonlinear fluid in the middle of the corresponding walls was discussed by Kaushik et al. [2]. Sheikholeslami et al. [3] offered an investigational study for entropy generation and exergy loss for nano-refrigerant compression progression. Afridi et al. [4] evaluated the combined impact of the Lorentz force and second law analysis on nanofluid dissipative flow about a curved medium. Nouri et al. [5] and Dormohammadi et al. [6] explored entropy analysis for nanofluid flows inside a channel. Convection and heat transfer of nanofluid with analysis of entropy inside a square cavity was inspected by Shermet et al. [7]. Entropy generation in unsteady magneto flow with a non-natural neural network and particle swarm optimization system was explored by Rashidi et al. [8]. Yongbo et al. [9] reported electro magnetohydro dynamics (EMHD) flow with entropy generation in a curved quadrilateral microchannel. Exploration of magneto slip flow under entropy generation on a revolving porous disk with adjustable properties was presented by Rashidi et al. [10]. Rehman et al. [11] proposed a study for the influence of radiation and thermal slip on rotating nanofluid with the entropy phenomenon. Researchers [12,13,14,15,16,17,18,19] have investigated the entropy generation on fluid flow with various physical effects. Non-Newtonian fluids are generally involved in several organic states and have natural and manufacturing uses. Unlike Newtonian fluids, the non-Newtonian fluids own viscosity that is reliant on the shear rate or strain rate history. This compound feature makes non-Newtonian fluid valuable for various particular applications, e.g., coarse particles with transportation [20], while being relatively difficult to guess. In previous years, countless energy has been utilized to streamline the bond between the viscosity and the flow state of the non-Newtonian fluid, as it is crucial to many mathematical models to define this kind of fluid [21,22,23,24,25]. Several organic and natural fluids, such as body fluid and wet sand at the seaside, exhibit shear-thinning and shear-thickening rheological features. More recently, power-law fluid flow over a bluff frame has received considerable attention for some manufacturing products, such as paper-making, processing of food, and heat exchangers. Power-law fluid flow has a wide series of applications, including in food and chemical processing, pharmaceuticals, oil production, power generation devices, and heat exchangers. Power-law fluid [26,27,28,29], by amending its index, can be set to define the comprehensive activities of shear-dependent fluids. The power-law index indicates Newtonian fluid for n = 1, and otherwise characterizes shear-thinning (n < 1) or shear-thickening fluid (n > 1) correspondingly. Mohebbi et al. [30] deliberated on the power-law fluid flow inside two corresponding plates for various Reynolds numbers. Flows of fluid over a permeable medium are of vital significance in power metallurgy, industrial filtration, petroleum technology, groundwater hydrology, ceramic engineering, etc. In the spirals of the geothermal region, water is an electrically conducting fluid due to peak temperature. The abundant power present in the soil’s layer in the geothermal regions has to be brought up to augment fuel output. This is an example to illustrate the flow past porous medium, taking the earth’s surface as a permeable bed. The temperature in the hot springs can be considered with the help of the energy equation. This temperature is useful to operate turbines through a magnetic field to produce electrical energy. In view of the aforementioned importance, investigators started working on porous layers with numerous geometries [31,32,33]. Several important and relevant kinds of literature on exponential stretching flows and different physical aspects, such as a non-uniform heat source and magenetohydro dyanmics (MHD), are listed in the references [34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54]. To the best of the authors’ knowledge, there seems to be no existing document on entropy analysis and a non-uniform heat source on power-law fluid with exponential stretching with heat flux conditions. Therefore, in this scientific report, our goal was to analyze the involvement of entropy generation with nonlinear power-law fluid above an exponential continuous moving surface. Both the shear-thinning and shear-thickening behaviors are discussed. The impact of a non-uniform heat source and suction/injection are also encouraged. The flow was controlled with the magnetohydrodynamics regime. The Runge–Kutta–Fehlberg method helps us to obtain the numerical explanation of the emerging equations. The model outcomes and comprehensive debates are provided in the form of results and discussions. In the conclusion, some key outcomes are also presented. The results obtained and presented in this article will portray that the optimal design and the efficient performance of a flow system or a thermally designed system can be enhanced by taking the suitable values of the physical parameters. This will allow us to reduce the effects of entropy generated within the system.

2. Definition of the Problem

As described schematically in Figure 1, in this paper we propose to study the boundary layer flow and heat transfer characteristics of an incompressible electrically conducting power-law fluid passing over an exponential moving sheet with the presence of a variable heat source/sink . With respect to assumptions, it is supposed that the stretching surface is non-uniformly heated via an imposed heat flux and moved horizontally in the x-direction with a non-uniform velocity , where x is the coordinate measured along the surface of the sheet. In addition, the occurring laminar flow of the non-Newtonian fluid is affected by the presence of horizontal stretching forces and the velocity of suction/injection at the permeable horizontal boundary . An external applied magnetic field vector is functional perpendicular to the stretched surface. The total magnetic field exerted on the fluid flow can be regarded as the sum of the external magnetic field and the induced magnetic field. Physically, the magnetic Reynolds number is defined as the ratio of the advection term to the magnitude of the magnetic diffusion term. For similar flows, the magnetic Reynolds number is very small. Hence, the induced magnetic field can be neglected for the problem under consideration.

Figure 1

Geometry of the physical problem.

The energy equation is more prominent with the variable heat source/sink term. Moreover, an analysis for entropy generation and the Bejan number was also carried out to focus on the importance of the existing flow regime. Under these suppositions, the governing partial differential equations of the present boundary layer flow are stated as with the following appropriate boundary conditions where represents the consistency coefficient, is a characteristic length, refers to the electrical conductivity, , , and are the reference velocity, magnetic field, and heat flux, respectively, and is the velocity of suction/injection, which is negative for suction and positive for injection . The particular problem can be specified in another straightforward form by using the following dimensionless variables Here, is the generalized Reynolds number, where . By virtue of the above considerations, the velocity components can be written in the form From Equation (6), it is obvious that the continuity equation (i.e., Equation (1)) is routinely satisfied. The non-uniform heat source/sink term used in Equation (3) is defined as where A* and B* are the coefficients of space and the temperature-dependent heat source/sink, respectively. By making use of Equation (6), the boundary layer Equations (2) and (3) together with the boundary conditions in Equation (4) are reduced to Here, is the magnetic parameter, is the generalized Prandtl number, is the thermal diffusivity, and is the suction/injection parameter.

3. Relevant Physical Measures

Keeping in mind the above considerations, the local skin friction coefficient and the heat transfer coefficient for the concerned problem are given by

4. Entropy Generation and Bejan Number

By employing the second law of thermodynamics in the present MHD boundary layer flow problem, the volumetric entropy generation is defined locally by [16,17,18,19] Where refers to the temperature gradient vector, denotes the viscous dissipation function, and represents the current density vector, such that Consequently, the volumetric entropy generation can take the following form in which Therefore, Equation (15) becomes The dimensionless form of Equation (19) is called the entropy generation and symbolized by . This thermodynamic quantity is given by where is a characteristic entropy generation, is the Reynolds number, is the Brikman number, and is the temperature difference parameter. Also, the Bejan number is introduced to check the domination of the heat transfer part over the fluid friction part along with MHD. This dimensionless quantity is defined as where is the irreversibility quotient. From Equation (21), it is observed that the heat transfer is higher when , whereas the other parts can dominate when . Also, if , the influence of both parts is of the same magnitude. As proven above, the Bejan number ranges from 0 to 1 (see reference [55] for more details). In addition, the heat transfer irreversibility can dominate when . Moreover, the frictional and magnetic irreversibilities became more prominent when Furthermore, the heat transfer and the other parts have an equal contribution when . Moreover, the manner of the Bejan number is premeditated for the optimal values of the constraints for which the entropy generation takes its lowest value.

5. Solution Methodology and Validation of Results

The nonlinear differential Equations (8) and (9) along with their corresponding boundary conditions (10) and (11) constitute a two-point boundary value problem. These differential equations are solved numerically by utilizing the Shooting Method (SM). For this purpose, the resulting differential system arising from Equations (8)–(11) is converted into an Initial Value Problem. By applying this method, Equations (8) and (9) can be reduced to a system of first-order ordinary differential equations by setting After introducing Equation (22) into Equations (8) and (9), we obtain the following reduced differential equations In this problem, the SM is used as a powerful technique to guess the missing initial conditions and by means of an iterative process until the boundary conditions are satisfied, in such a way that the reduced Equations (23)–(25) were integrated numerically by utilizing the Runge–Kutta–Fehlberg method and taking as the best step size. Then, the computed values of and at were compared with the given boundary conditions and . The guessed values of and were refined with the help of the Newton–Raphson Method to give a better approximation for the desired solutions. The iterative process was repeated until we acquired results with an accuracy level of about of . In order to verify the precise running of the program, the numeric values of the wall shear stress for the viscous fluid case were related to the findings of Magyari and Keller [56], Elbashbeshy [57], Sahoo and Poncet [58], Mukhopadhyay et al. [59], and Sajid and Hayat [60] in their pioneering studies. As expected, it was found that our results are in good agreement with those in the existing literature [56,57,58,59,60] (see Table 1).

Table 1

Comparison between our findings and those of the existing literature results [56,57,58,59,60], in the case where and .

Existing Results for f″(0)					Present Result for f″(0)
[56]	[57]	[58]	[59]	[60]	SRKFM
−1.28180	−1.28181	−1.281811	−1.281812	−1.28181	−1.2818098

6. Results and Discussion

This segment of the present work highlights the impression of related variables on exponential flow of non-Newtonian power-law fluid. Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20, Figure 21 and Figure 22 are plotted to explain the variations of velocity, temperature, entropy analysis, the Bejan number, and flow and heat fluxes in both cases of shear thinning and shear thickening. Figure 2 is dedicated to observing the effect of the magnetic parameter on the velocity of fluid. It was found that the resistive Lorentz force due to a magnetic field reduces the motion of the fluid. The magnitude of momentum boundary layer was higher for the shear thickening case, i.e.,. Figure 3 discloses the effects of suction/injection for velocity with and . It was observed that the thickness of the boundary layer shrinks extra-rapidly for the case. The impact of internal heat source parameters and with and is sketched in Figure 4, Figure 5, Figure 6 and Figure 7, respectively. The temperature of the considered fluid increased in all cases and the enhancement of the thermal boundary layer was maximum for and , respectively. Figure 8 and Figure 9 are presented to display the nature of the temperature for suction/injection with and , respectively. It is clear that this parameter is an inverse function of temperature and there is less variation in magnitude for the and case. The power-law index causes a decrement in temperature (see Figure 10). Figure 11 elucidates the effect of the heat flux index on the temperature profile. Obviously, the thermal boundary layer thickness decreases with a higher The impact of a magnetic constraint on entropy is shown in Figure 12. It is depicted that the Lorentz force created by the magnetic field augmented the entropy of the system. Evidently, is very sensitive to entropy generation augmentation. This is due to the fact that the magnetized fluid forced the dissipative heat energy to thermal diffusion. Figure 13 and Figure 14 were developed to study the variation of entropy generation for the group parameter and the Reynolds number with , respectively. Both parameters are directly proportional to the entropy analysis and the magnitude is almost similar for both cases.