Insight into the Dynamics of Fractional Maxwell Nano-Fluids Subject to Entropy Generation, Lorentz Force and Heat Source via Finite Difference Scheme.

In recent times, the loss of useful energy and solutions to those energy challenges have a wide scope in different areas of engineering. This work focuses on entropy analysis for unsteady viscoelastic fluids. The momentum boundary layer and thermal boundary layer are described under the effects of a magnetic field in the absence of an induced magnetic field. The study of a fractional model of Maxwell nanofluid by partial differential equation using Caputo time differential operator can well address the memory effect. Using transformations, the fractional ordered partial differential equations (PDEs) are transfigured into dimensionless PDEs. Numerical results for fractional Maxwell nanofluids flow and heat transfer are driven graphically. The Bejan number is obtained following the suggested transformation of dimensionless quantities like entropy generation. A mathematical model of entropy generation, Bejan number, Nusselt number and skin friction are developed for nanofluids. Effects of different physical parameters like Brickman number, Prandtl number, Grashof number and Hartmann number are illustrated graphically by MAPLE. Results depict that the addition of nanoparticles in base-fluid controls the entropy generation that enhances the thermal conductivity and application of magnetic field has strong effects on the heat transfer of fractional Maxwell fluids. An increasing behavior in entropy generation is noticed in the presence of source term and thermal radiation parameter.

1. Introduction

The best way to represent the natural phenomena is by using differential equations (DEs) with suitable boundary conditions. Recently, the fractional (non-integral) order of DEs has gained much interest because of their vast scope in many engineering fields like food engineering, oceanography, chemical reactions and chaos. All fractional derivatives like Riemann Levillie, Caputo, Caputo Fabrizio, Antangna Beleanu are used widely in practice [1,2,3], as a deep study can easily be handled using this approach. Natural convection along a vertical wall and cylinder has been explained using Caputo time-fractional derivatives and fractional derivatives [4,5]. The significant difference between fractional and ordinary fluid flow emerges at different times. Furthermore, with large time values, increasing the fractional parameter’s value increases velocity. Fractional calculus covers the complex structure of viscoelastic fluids in various research areas of glass fiber production, exotic lubricants, colloidal solutions, extraction of polymer solutions and cooling processes [6]. A simplified Phan–Thien–Tanner model for viscoelastic fluids was investigated analytically in [7] with different physical parameters. It was established by Sheikh, N.A. et al. [8] that it is not enough to get experimental data by using a conventional derivative model for Maxwell fluids rather than fractional operators. Moreover, it was discovered that a link exists between Maxwell’s constitutive equation and molecular theory [9]. Few other investigations for analytical results have been done in many research articles. But in the rheological perspective, viscoelastic fractional order models have attracted much interest due to their wide range of applications. Fractional Maxwell models are derived by substituting the conventional derivative in the known Maxwell model of stress-strain expression by fractional-order derivatives. Analytical solutions of fractional Maxwell, fractional generalized Maxwell model, fractional second grade, third grade, Oldroyd-B model etc., via analytical techniques of Laplace transform, Fourier transform, Weber transform and Hankel transform are obtained in [10,11,12,13,14].

Many researchers focus on such physical processes involving entropy generation-EG. It is a well-known fact that all physical problems, especially heat transfer, involve entropy generation—EG. Entropy generation-EG plays a vital role in fluid dynamics. The first law of thermodynamics moves around heat transfer processes, whereas the second law of thermodynamics moves around the entropy generation of the system. Entropy generation—EG tells the feasibility and efficiency of the system. In other words, entropy describes the ways a process can control energy loss.

Since entropy generation—EG is important and happens in almost all thermo-dynamical processes, many researchers have worked in this direction. Among the various researchers who have done splendid work on entropy generation–EG, Adrian Bejan, has published many articles and books [15,16,17,18,19,20,21]. Considerable work is done by Pranab Kumar Mondal [22] for irreversibility analysis of Couette flow while applying weak and relatively strong pressure gradient. The fluid dynamics are studied by variation of volumetric entropy generation number and Bejan number.

To control the entropy generation—EG, addition of nanoparticles in base fluid gave a new direction to heat transfer problems (the formation of nanofluids). Nanofluids are nano-sized particles accomplished with base-fluids, i.e., Common nanoparticles Water, engine oil and ethylene glycol are common base fluids. To enhance the thermal conductivity of base fluids, the idea of nanofluids is given by Choi et al. in [23] for the first time. Later on, Tiwari and Das discussed the effects of different shapes of nano-sized particles on thermal enhancement [24]. Using square shape cavity, different aspects of ball-shaped, cylindrical-shaped, and rod-shaped nanoparticles were studied [25]. Analytical results for the effects of nanofluids were driven by [26]. In [27], analytical results are investigated for temperature profile and Nusselt number under the effects of viscous dissipation and porous media. Using a traditional approach, numerical results were obtained, expressing the radial and tangential momentum across the disk decreases for higher Lorentz forces and slip factor [28].

Like many other key sources, entropy generation—EG includes viscous dissipation, chemical reaction, heat and mass transfer, heat convection and conduction and electrical conduction. Considerable work by B. Mahathish et al. has been done via the spectral quasi-linearization (SQL) method for entropy analysis and solving a Williamson model, which provides a base for fractional models. An analytical approach is applied to investigate the effects of viscous dissipation and limiting effects of Nusselt for temperature profile [29]. Bejan number-BN is the ratio of entropy generation—EG due to heat and total entropy generation—EG of system. But Awed in [30] gave a new definition to Bejan number-BN. Bejan number-BN describes the effects of magnetic field irreversibility and fluid friction irreversibility. An investigation has been conducted on laminar falling liquid film along an inclined heated surface in [31]. In [32,33,34] and many other articles, the exact and numerical solution of entropy generation—EG is published. B. Mahathish numerically investigates the effects of quadratic variation of density-temperature (quadratic convection) and the quadratic Rosseland thermal radiation using the modified Bongiorno Model (MBM) [35]. Similar analytical as well as numerical approaches can be seen in [36,37,38,39] for such different viscoelastic models.

Numerical investigations and many others were used in [40,41,42,43]; however, various research gaps are found in these articles, which are still not addressed. Such as how heat transfer can be enhanced by adding the nanoparticle to the base fluid? What is the effect of Lorentz forces on flow dynamics? What do numerical results predict about Skin friction and Nusselt number? What are the formulation of entropy generation and Bejan number in the presence of thermo-physical properties of nanoparticles? In reviewing these research gaps in literature surveys, the main task of this article is to develop a mathematical model of momentum and heat for fractional Maxwell nano-fluids for detailed insights into Lorentz force, heat source/sink, and Nusselt number and entropy generation—EG. A new definition of Bejan number—BN is introduced by adding the coefficients of thermo-physical properties of nanoparticles. Water is taken as base fluid, whereas are the nanoparticles used for graphical results of velocity and temperature profiles. The problem is first modeled fractionally by applying the definition of Caputo time derivative, then using transformation; a dimensionless analysis is done. The resulted model is solved by a numerical technique of finite difference scheme. Plots are drawn for Bejan number—BN against the Brickman number, the square of Hartmann number and dimensionless temperature difference. Moreover, some graphical results are extended to evaluate the Nusselt number and Skin friction . These results are computed via mathematical software MAPLE.

2. Mathematical Model and Formulations

Considering magnetohydrodynamic (MHD) flow of incompressible and unsteady fluid along the infinite vertical plate. The induced magnetic field and pressure gradient are neglected. Initially, at time fluid has velocity zero and has a constant temperature . With the passage of time temperature of the system rise to . The fluid flow is considered along . The magnetic field is applied in the as illustrated in Figure 1 below.

Figure 1

Geometry of the problem.

Taking into account the Boussinesq approximation, the assumptions of the system are,

Flow is unsteady, incompressible and 1-dimensional.

Pressure gradient is absent.

Body force is significant.

Magnetic field is applied (ignoring the induced magnetic field).

Then the equation of continuity is restricted and takes the following form (can be seen in [44]);

But the Navier–Stoke equation [45,46] takes the form, where are dynamic viscosity of nanofluid, Cauchy stress tensor, current density, total magnetic field, gravitational acceleration, thermal expansion coefficient, the temperature of nanofluid and ambient temperature, respectively.

The stress tensor for Maxwell fluids in [47] as, where

In these expressions represents the extra stress tensor, Rivline–Ericksen tensor, kinematic viscosity and material time derivative. Also is expressed in [48] as;

And

By following the assumption of the problem, the Maxwell equation can be written in [49] as;

Keep in view that is the sum of applied and induced magnetic field (neglected).

Using Equations (3)–(7), the Equation (2) takes the form;

Multiplying both sides of Equation (8) by

But the fractional constitutive equation for Maxwell fluids is given by [50],

This constitutive relation contains , Caputo fractional operator is defined in [51] as;

With as the Gamma function defined in [51] by;

Using Equation (10) into Equation (9)

The first law of thermodynamics [46,52] is;

In this equation are constant of heat capacity, coefficient of thermal conductivity, radiative heat flux, and nanofluid thermal conductivity, respectively.

By using the Rosselands approximation for fluids that are considered optically thick, the radiative heat flux expressed in [53] is given as;

In this expression are the Stefan–Boltzmann constant and mean spectral absorption constant, respectively. Approximating by a Taylor’s series expansion in the neighborhood of and neglecting higher power.

(neglecting higher power), so radiative heat flux is

Then Equation (12) becomes

The second law of thermodynamics is given by [46],

In which are volumetric local entropy generation and electrical conductivity.

The proposed boundary and initial conditions of this physical phenomenon are defined in [54] below;

Using the following transformation and thermophysical properties of nanoparticles (can be seen in) [55], the dimensionless governing equation for velocity and temperature profile are obtained;

The dimensionless velocity and temperature profile of the problem is given, and after omitting notation for the sack of brevity of mathematical modeling.

Moreover, the temperature equation takes the following form

In Equations (22) and (23) are the ratio of thermophysical properties given by

Additionally, are the Grashof number, Hartmann number, Prandtl number, radiation parameter and heat generation parameter, respectively, defined in [45] as;

Additionally, the non-dimensional initial and boundary conditions are

The non-dimensional governing equation for velocity and temperature profile in Equations (22) and (23), with dimensionless initial and boundary conditions in Equations (24) and (25), represents the unsteady, incompressible flow fractional Maxwell nanofluids phenomena under the influence of magnetic fields. Water is taken as base fluid, but are the nanoparticles considered for nanofluid preparation. For the numerical results the following Table 1 containing thermo-physical properties of nanoparticles and base fluid will be utilized.

Table 1

The thermophysical properties of different base fluids and nanoparticles at room 25 °C [56].

Material	H_2O	Cu	Al_2O_3
ρ(kgm^−3)	997	8933	3970
C_p(J Kg^−1k^−1)	4197	385	765
k(Wm^−1k^−1)	0.613	400	40
β×10^−5(k^−1)	21	1.67	0.85
σ(Ωm)^−1	0.05	5.96×10^7	2.6×10^6

3. Numerical Procedure

The finite difference scheme is a very efficient and powerful tool to investigate the numerical solutions of the problem arising and mathematical physics and mechanics. In this context, this section is dedicated to extending the finite-difference scheme to tackle the obtained set of fractional-order fluid models and heat transfer. For this, the discretization of the derivative of fractional-order of and are specified as, and the nonlinear term is approximated by means of the following concept where, for Now, the rectilinear grid is assumed to examine the solution of the governing set of fractional-order fluid problems and heat transfer having grid spacing in the direction of space and time separately, where for from . The inner points in the discussed domain are given as and . The discretization of the governing set of fractional-order fluid problems and heat transfer at is given as, for , with the following initial and boundary conditions,

3.1. Numerical Analysis and Discussion

Test Problem. For the validation of the applied scheme, a test problem is considered as

In the discussed problem, the conditions are given below, and the source term can be selected against the choice of fractional-order derivative.

Since such a physical problem contains as the exact solution. Its accuracy has been checked by a number of simulations for the proposed scheme. The plots in Figure 2a,b are drawn for the maximum absolute error (MAE) and computational order of convergence (COC) for different ranges of N, which is N = 10, 20, 40, 80, 160, 320, 640.

Figure 2

Code validation of proposed scheme and varying time mesh size against (a) computational order of convergence (COC) (b) maximum absolute error (MAE), and varying mesh size for (c) time and (d) space against L∞-norm between consecutive solutions.

The convergence of the applied scheme is observed against the selection of each fractional-order derivative, and its convergence order enhances for . Figure 2c,d contains the L∞-norm between consecutive solutions that is and when 0 ≤ i, j ≤ N, M = 500. Again, it is found that the proposed scheme is very efficient, accurate and reliable for this problem. It also demonstrates that the solution is stable against the selection of fractional order and mesh parameters.

3.2. Entropy Generation

For viscous fluid flow in a magnetic field, the volumetric rate of local entropy generation is defined in [57] as;

That is the sum of entropy generation due to heat transfer, due to fluids friction and due to the magnetic field effect, separately mentioned here (can be seen in [33])

Combining all the results

The dimensionless entropy generation calculated is;

where

where , is the characteristic entropy generation rate, is the Brickman number, is the square of the Hartmann Number and is the dimensionless temperature difference.

Bejan Number

The Bejan number is the irreversibility distribution parameter which is expressed mathematically as (can be seen in [58]);

That is

Then the reduced expression for the Bejan number is given as;

With

3.3. Skin Friction and Nusselt Number

For measuring shear stress and heat transfer effects in an ordinary integer order system, local skin friction and Nusselt number are defined in [59] as; and

The skin friction coefficient and local Nusselt number for (FMF) can be written by using Equation (5), that is, the fractional stress tensor for Maxwell fluid on the plate with fractional time Caputo derivative (detail can be seen in [60]).

The non-dimensional form of Equations (38) and (39) is given as

4. Results and Discussion

This section of the article deliberates the detailed results and discusses the plots driven against different physical parameters like representing the magnetic field parameter (the square of Hartmann Number), Prandtl number, radiation parameter, Grashof Number, volumetric fraction of nanoparticle, fractional order parameters and heat generation parameter, respectively. The behavior of these aforementioned physical parameters on dimensionless velocity profile , temperature profile , Bejan number-BN , Skin friction and Nusselt number are drawn by mathematical software MAPLE. The mathematical fractional model of Maxwell nanofluid is developed by using the Caputo time fractional operator. After transforming the governing equations to a dimensionless governing model, the finite difference method is used for the discretization of the model. FDM is a strong tool for dealing with such kinds of problems. The simulation is performed by developing and executing codes.

Results are obtained by solving Equations (22), (23), (32), (39) and (40) with initial and boundary conditions illustrated in Equations (24) and (25) and physical properties of nanoparticles in Equation (21) and Table 1. Various suitable ranges of physical paramet for heat transfer, velocity analysis, skin friction, Nusselt number, entropy generation and Bejan number—BN are considered, and also particular exertion has been done on the effects of these parameters for heat enhancement.

Figure 3 shows the effect of the magnetic field parameter on the velocity profile of fractional Maxwell fluids. It can be seen that velocity decreased by increasing the value of because of the Lorentz’s force. The rise in caused strengthening in the Lorentz force, which increases the internal resistance to flow particles; consequently, fluid velocity decreased. Whereas an increase in Skin friction occurs, as shown in Figure 3b. Since Bejan number is the ratio of total entropy generation to entropy generation due to heat transfer. Thus, advancement in boosts the Bejan number-BN and can be noticed in Figure 3c. An opposite behavior of total entropy generation relative to Bejan number—BN appears against the increasing value of , as can be seen in Figure 3d.

Figure 3

Impact of the Hartmann number on (a) dimensionless velocity profile, (b) coefficient of skin friction, (c) Bejan number and (d) entropy generation when .

Figure 4 depicts the effects of the Grashof number on velocity profile . The enhancement in value of results in increasing the fluid velocity, which can be physically justified as the increasing value of means lower the viscous forces and hence increasing the velocity of fractional Maxwell fluids. While Skin friction increases, the increasing the value of can be noticed in Figure 4b. A slight effect of on Bejan number—BN and entropy generation can be seen in Figure 4c,d, respectively. Similar results are reported by Sarojamma, G., et al. in [61] for comparison. is the ratio of buoyancy forces to viscous force, increasing value of results in laminar boundary layer and vice versa. Because high value of give rise to the temperature of molecules, consequently, the intermolecular forces become weak. Thus, velocity profile also rises. On the other hand, fluid particles collectively gain momentum as increases, so additional heat is lost nearby, that is why skin friction reduces as shown in Figure 4b. Since Bejan number is ratio of entropy generation due to heat to the total entropy generation of the system, that is why reduces the value of Bejan number—BN deliberated in Figure 4c.

Figure 4

Influence of Grashof number on (a) dimensionless velocity profile, (b) coefficient of skin friction, (c) Bejan number and (d) entropy generation when .

Figure 5 shows the effects of the volume fraction parameter on fluid velocity, and it is examined that advancement in controlled volume fraction parameter lowers the velocity profile due to the effect; that is, the viscosity of fractional Maxwell fluids increases by increasing . On the other hand, the addition of nanoparticles in the base fluid causes an improvement in the heat transfer rate at the boundary layer. As the thermal conductivity of base fluid is enhanced, consequently increasing the fluid’s internal temperature can be seen in Figure 5b. Since Nusselt number is the ratio of convective heat transfer to conductive heat transfer, increasing has a decreasing relation with Nusselt number . This is because the skin friction decreases with the passage of time against the volume fraction parameter . These results for skin friction and Nusselt number are viewed in Figure 5c,d, respectively. An increasing trend is noticed for the Bejan number and can be seen in Figure 5e as the heat transfer rate becoming better by increasing . Whereas in Figure 4f entropy generation enhanced by increasing the value of . Because the viscosity of fluid increased by increasing .

Figure 5

Effect of on (a) dimensionless velocity profile, (b) non-dimensional temperature profile, (c) coefficient of skin friction, (d) local Nusselt number , (e) Bejan number and (f) entropy generation when .

In these following plots in Figure 6, results are drawn for velocity profile , temperature profile , skin friction , Nusselt number and Bejan number-BN against fractional parameter .

Figure 6

Impact of fractional parameter on (a) dimensionless velocity profile, (b) non-dimensional temperature profile, (c) Bejan number and (d) Entropy generation when .

The reason is that the gradual increase in fractional parameter gives rise to the viscosity of the nanofluid. This means intermolecular forces between the nanoparticles and base-fluid particles increase; consequently, Brownian motion of particles reduces, that is, a decrease in velocity profile occurs. The consequences of fractional order parameter for velocity profile are inverse as depicted in Figure 6a. On the other hand, have direct relation for , as shown in Figure 6b. Figure 6c shows that for an increasing value of the entropy generation decreases, whereas a rise in varies directly, deliberated in Figure 6d.

Figure 7 shows the effects of thermal radiation parameter on temperature of the fluid. Applying thermal radiations gives rise to the temperature of particles of nanofluids. Hence, the particles’ kinetic energy increases, and the rate of collision between the particles of the nanofluids becomes high, which is why a rising increase in temperature profile occurs. Therefore it is concluded that increasing the value of causes an increase in fluid temperature. is the ratio of convective to conductive heat transfer across the boundary, but enhancement in conduction occurs with the addition of nanoparticle therefore, a decrease occurs in , depicted in Figure 7a,b, respectively. A decrease in Bejan number—BN can be noticed in Figure 7c, whereas entropy decreases initially and then increases gradually, as can be seen in Figure 7d. Similar results are reported in [61].

Figure 7

Influence of thermal radiation parameter on (a) dimensionless velocity profile, (b) non-dimensional temperature profile, (c) Bejan number and (d) entropy generation when .

Prandtl number is the dimensionless number and is the ratio of momentum to thermal diffusivity. It is a fluid property but does not have any dependence on flow type. Thus, an increase in means heat transfer is favored to occur by momentum, not conduction. This parameter controls the relative thickness and thermal boundary layer in heat transfer problems. Lowering the value of means the heat diffuses spontaneously as compared to momentum, which thickens the thermal boundary layer rather than the momentum boundary layer. Therefore, an increase in decreases the temperature profile of fractional Maxwell fluids as expressed in Figure 8a, which expectedly decreases the Nusselt number . Since Bejan number has an inverse relation with entropy generation , due to heat transfer, that is, decreases with the increase in value of as illustrated in Figure 8.

Figure 8

Effect of on (a) dimensionless velocity profile, (b) coefficient of skin friction, (c) Bejan number and (d) entropy generation, when .

Figure 9 shows the effects of the heat source term on temperature profile , since heat source gives rise to temperature profile as shown in Figure 9a, but with the passage of time, it is noticed that the Nusselt number decrease is dependent inversely on the conduction of heat. Thus, entropy generation expectedly decreases as shown in Figure 9d. This phenomenon gives rise to the Bejan number because entropy generation due to heat transfer has an inverse relation with the Bejan number , as shown in Figure 9c.

Figure 9

Influence of heat source on (a) dimensionless velocity profile (b) non-dimensional temperature profile (c) Bejan number and (d) entropy generation when .

5. Conclusions

The graphical analysis of fractional Maxwell nanofluids is made in this article under the influence of a magnetic field (ignoring induced magnetic field). The pressure gradient is supposed to be absent. Effects of different physical parameters are drawn by using the mathematical software MAPLE. The model is formulated by applying the Caputo time derivative. Using suitable transformations, governing equations are made dimensionless.

MHD fractional Maxwell nanofluids are studied numerically; quantities like entropy generation, Bejan number, Skin friction, and Nusselt number are investigated using the finite difference method.

Hence key findings of this study are given below;

It is noted that for increasing the value of fractional order parameter the velocity profile decreases, whereas the temperature profile increases.

The addition of nanoparticles to base fluid enhances the thermal conductivity of fractional Maxwell nanofluids, increasing the value of volume fraction of nanoparticles and decreasing entropy generation .

The magnetic field effect influences the temperature and velocity profile with inverse and direct behavior, respectively.

Nusselt number increases with the variation in and a decrease occurs in with the increase in thermal radiation parameter.

The temperature profile varies directly with the thermal radiation parameter and increasing the value of decreases the Nusselt number. Whereas entropy generation increases, and the Bejan number decrease with a rising value of .

The solution obtained via the finite difference method is excellent in agreement with the test problem and existing results which shows that the finite difference method is a strong and reliable technique to deal with such kind of complex models and it gives a key direction for further study.