Investigation of thermal performance of Maxwell hybrid nanofluid boundary value problem in vertical porous surface via finite element approach.

The study of thermo-physical characteristics is essential to observe the impact of several influential parameters on temperature and velocity fields. The transportation of heat in fluid flows and thermal instability/stability is a charming area of research due to their wider applications and physical significance because of their utilization in different engineering systems. This report is prepared to study thermal transportation in Maxwell hybrid nanofluid past over an infinite stretchable vertical porous sheet. An inclusion of hybrid nanofluid is performed to monitor the aspects of thermal transportation. Keeping in mind the advantages of thermal failure, non-Fourier theory for heat flux model is utilized. Aspects of external heat source are also considered. The mathematical formulation for the considered model with certain important physical aspects results in the form of coupled nonlinear PDEs system. The obtained system is reduced by engaging boundary layer approximation. Afterwards, transformations have been utilized to convert the modeled PDEs system into ODEs system. The converted nonlinear ODEs system is then handled via finite element method coded in symbolic computational package MAPLE 18.0. Grid independent survey is presented for the validation of used approach and the comparative analysis has been done to confirm the reliability of obtained solution. The obtained solution is discussed and physical aspects have been explored and recorded against numerous involved influential variables. Motion into hybrid nanoparticles and nanoparticles becomes slow down versus higher values of Forchheimer and Darcy's porous numbers. Thermal growth is enhanced for the case of hybrid nano-structures rather than for case of nanofluid. Thickness regarding momentum layer is dominated for hybrid nanoparticles rather than case of nanoparticles.

Introduction

Modeling of complex phenomenon to understand the mathematical physics behind them is the growing and useful tool for researchers, mathematicians and engineers. Several important empirical relations have conducted and proposed to model the nature of complex phenomenon arising in mathematical physics.

Several important contributions have been performed so far on Maxwell model. For instance, Mushtaq et al.[1] examined the buoyancy force effects by considering generalized heat flux theory under variable and constant wall temperature past over an elastic vertical plate. They examined the impact on opposing and assisting flow for different emerging parameters. They recorded the increase in temperature field against buoyancy parameter and decrease in velocity field against fluid relaxation time. Zhang et al.[2] computed exact solution of unsteady fractional Maxwell model and thermal transport mechanism is highlighted via Laplace transform procedure. They have plotted several graphs and analyzed the contribution of numerous emerging parameters in the presence/absence of slip factor. They have shown the decline in temperature profile against fractional parameters and upsurges is monitored in velocity against the said parameters. Numerical approach coupled with Lie group investigation has been applied by Shafiq et al.[3] to examine the flow of Maxwell model past over a stretched penetrable sheet. They noticed the opposite behavior in thermal profile for heat absorption/generation parameter. Moreover, they noticed the increase in fluid velocity for Deborah number. Raza and Ullah[4] used two different fractional derivatives approaches to study thermal transportation in unsteady Maxwell model. They computed the solution via Laplace procedure. They examined the increase in fluid velocity against Maxwell fluid parameter and Grashof number. Kaushik et al.[5] examined the flow of Maxwell model in rotating channel. They obtained streamlines for velocity solution and impact of numerous parameters on solution is explored and explained in detail. Thermal transportation in MHD fractional Maxwell model with dissipation was studied by Bai et al.[6]. They used numerical and analytical approach to discuss the solution behavior against numerous involved parameters. They monitored the decline in velocity for magnetic parameter and similar bearing of Prandtl number in thermal field. Iqbal et al.[7] worked on convective Maxwell model obeying peristaltic transport phenomenon in a channel. Transportation of heat and mass is also investigated by them. They recorded the trapping phenomenon. They examined the depreciation in concentration against Schmidt and Soret numbers. Furthermore, opposite contribution of Biot number is monitored for concentration and temperature profiles.

Extensive utilization of nanoparticles[8-10] in different mechanisms including medicines, cancer treatment, machineries and different other household instruments make a remarkable attention of this topic. Growing researches and development in technology accept the involvement and significance of nanoparticles. Several researchers pay their attention and contributed on the topic. For instance, Abdelsalam and Sohail[11] studied bi-directional bio-convective stretched flow with variable thermophysical attributes. They obtained the resulting reduced PDEs after utilizing boundary layer theory. They solved the converted ODEs from PDEs with the addition of suitable transformation by OHAM tool. They established that higher order approximations reduce the error. They recorded the depreciation in velocity field against magnetic parameter and upsurges in temperature profile. Furthermore, decline in concentration profile is observed against Schmidt number. Mixing of graphene nanoparticles in Eyring Powell MHD unsteady model with momentum slip past over a stretched sheet was examined by Khan et al.[12]. Ghadikolaei et al.[13] worked on magneto nanofluid Casson model in the presence of inclined non-uniform magnetic field. They considered radiation, viscous dissipation and Joule heating in temperature expression. Shooting approach has been utilized for the solution of nonlinear coupled ODEs. They recorded the augmentation in temperature field against Casson parameter and decline in fluid velocity. Ramzan et al.[14] addressed the mixed convection phenomenon in stretched nanofluid model of viscoelastic material with Soret and Dufour effects. They monitored the resistance in thermal profile against mounting values of Prandtl number. Rotating MHD CNTs mixed Casson fluid flow under radiation and heat generation effects was studied by Muhammad et al.[15]. Hady et al.[16] studied saturated yield exhibiting flow past over a porous surface. They treated nonlinear modeled equations numerically. They recorded the depreciation in Nusselt and Sherwood numbers against yield stress parameter. Nayak et al.[17] studied stretched MHD radiated viscous fluid over an exponential stretching sheet numerically via shooting method. They have examined the involvement of different nanoparticles and monitored the heat transportation rate. Aman et al.[18] examined the inclusion of different nanoparticles mixture in Maxwell model. Additionally, the involvement is SWCNTs and MWCNTs have been presented to notice thermal performance. They handled the resulting nonlinear equations analytically. Few important latest contributions containing several important aspects are covered in Refs.[19-21].

The manufacturing of hybrid nanoparticles has been developed by scattering nanofluids while hybrid nanofluid can boost ability of thermal properties in nanofluids. Main role of hybrid nanofluid has been improved thermal conductivity, low cost and stability. Such role of hybrid nanofluid makes it favorite in field of nanotechnology. Therefore, hybrid nanofluid has been produced great interest for researchers. Some study related to hybrid nanofluid is mentioned here. For example, Hanif et al.[22] considered heated cone to investigate thermal features under magnetic field imposing hybrid nanoparticles whereas they have estimated entropy generation. Khan et al.[23] studied thermal aspects using silicon dioxide and molybdenum disulfide over a 3D heated surface. Hanif et al.[24] studied performance of hybrid nanofluid under thermal radiation in heated cone in the presence of magnetic field numerical solved by Crank-Nicolson approach. Hanif et al.[25] scrutinized a novel investigation of viscosity model for unsteady flow using magnetic field imposing hybrid nanofluid. Saqib et al.[26] discussed heat energy characterizations in the presence of hybrid nanoparticles in form of fractional DEs (differential equations) using integral transforms. Saqib et al.[27] analyzed energy transfer under MHD flow in rheology of Maxwell liquid in view of Cattaneo–Friedrich model. Hanif et al.[28] investigated magneto-hydrodynamic characterizations in view of nanoparticles towards a vertical porous heated cone. Jamil et al.[29] investigated features of Maxwell martial considering effect of chemical reaction in form fractional derivatives using thermal radiation. Anwar et al.[30] studied microscopic investigation under magnetic field of inertial motion considering porous geometry. Dinarvand et al.[31] discussed improvement in micro-circulatory process inserting role of hybrid nano-structures over a porous surface. Mousavi et al.[32] investigated performance of Casson liquid and dual solutions of developed model using nanoparticles over an expanding surface. Dinarvand et al.[33] analyzed swirling flow in mass-based model using hybrid nanoparticles via von Kármán’s theory. Aghamajidi et al.[34] discussed MHD motion and connective heat transfer using Tiwari-Das model related to nanoparticles over a cone. Dinarvand et al.[35] developed boundary layer model inserting hybrid nano-structures. Dinarvand et al.[36] analyzed an effect of magnetic field using MgO–Ag nanoparticles over moving heated needle.

Two-dimensional model related to heat transfer for unsteady flow is developed past a vertical surface. A phenomenon for hybrid nanoparticles is imposed in the presence of non-Fourier’s theory for unsteady flow along with heat generation term. A rheology of Maxwell martial is studied with Darcy’s Forchheimer theory using various thermal properties. A finite element method is implemented to simulate numerical results. It is noticed that such developing model is not investigated yet. The conducted research is organized as: literature survey is included in “Introduction” section, modeling is covered in “Physical statement” section with important physical quantities, an effective and convergent as well as stable numerical approach is presented in “Procedure of finite element method” section, results with physical interpretation are mentioned in “Results and discussion” section and consequences related considered problem have been reported in “Consequences of current problem” section.

Physical statement

Characteristics of thermal energy in Maxwell liquid inserting dispersion of nanomaterials and hybrid nanomaterials are considered past melting vertical surface as shown in Fig. 1. Phenomena of heat energy characteristics are occurred due to heat generation and magnetic field. Mechanism of non-Fourier’s law is modeled in energy equation. Theory of Darcy’s Forchheimer is observed. The fluid runs over melting surface because of movement of wall. The behavior of fluid rheology is captured by Fig. 1. It is noticed that movement of wall is occurred along y-direction and constant magnetic field is assumed along normal direction. Amount of nanoparticles and hybrid nano-structures is considered above x-direction. A vertical surface is considered to analyze feature of considered effects. Governing physical system of equations contains the utilization of following:

Figure 1

Geometry of current model.

two dimensional unsteady stretching porous sheet;

mixed convective flow;

single phase model;

hybrid nanoparticles;

viscoelastic material (Maxwell model);

heat generation/absorption;

non-Fourier heat flux;

no slip theory;

boundary layer analysis.

The PDEs are modeled trough concept of boundary layer for the conservation laws and modeled PDEs[38] arewhere velocity is is the relaxation time, is the inertia coefficient, permeability is denoted by is the temperature, thermal relaxation time is heat generation is temperature in ambient fluid is is the thermal conductivity, fluid density is represented by specific heat is and stands hybrid nanofluid, gravitational acceleration is and is thermal expansion coefficient.

Boundary conditions regarding problem are

Here, and are constant number whereas and have unit. The relationships among thermo-physical properties of pure fluid, mono, hybrid nanoparticles and mixture fluid are mentioned[21,38] below and their values are mentioned in Table 1.

Table 1

Thermal properties[37,38] for silver, ethylene glycol and molybdenum disulfide.

MoS_2	Ag	EG
ρ_MoS_2=5060	ρ_Ag=10480	ρ_EG=1113.5
(C_p)_MoS_2=397.21	(C_p)__Ag=235	(C_p)__EG=2430
k_MoS_2=904.4	k__Ag=429	k__EG=0.253
σ_MoS_2=2.09×10^-5	σ__Ag=6.30×10^7	σ__EG=4.3×10^-5

Following change of variables is suggested for non-dimensional PDEs

The porosity, Deborah, Forchheimer, Prandtl, heat generation, thermal relaxation and unsteadiness numbers are captured below

Coefficient associated with drag force is

Rate of heat energy is

Reynolds parameter is illustrated by

Procedure of finite element method

Finite element method is implemented and computer code is developed to handle the complex coupled (ordinary differential equations) ODEs which results after modeling the Maxwell fluid model with heat transport. The procedure consists of following key steps.

Residuals of current ODEs are integrated and weighted and residuals of considered problem are

Here, , and are weight functions and unknown variables (, and ) are formulated in the presence of shape function () as

The weak procedures are made via weighted residuals.

Stiffness matrices are calculated by using Galerkin approximations (into weak procedures) whereas Stiffness matrices are also used in assembly approach. The stiffness elements are constructed as

System of algebraic equations (nonlinear) is obtained.

The algebraic system is linearized and solved iteratively. Convergence and mesh analysis are done and flow fields are recorded under the variation of physical parameters raised as a result of dimensional analysis. The error analysis of developed model is

Mesh free study is carried out and convergence is ensured.

Mesh free analysis is displayed in Table 2 considering is 8. It is observed that investigation related to mesh-free for 300 elements is simulated. Elements of desired domain are increased 30 elements to 300 elements whereas solution is converged for 300 elements and outcomes are not changed (affected) after 300 elements. It is concluded that solution is converged and solution becomes repeatable at 300 elements. Therefore, computation investigation is simulated within 300 elements of considered problem.

Table 2

Simulations for mesh-free in temperature and velocity observing 300 elements[38] when

Number of elements	f^′η_max2	θη_max2
30	0.0075101	0.53284
60	0.0070855	0.51706
90	0.0070129	0.50973
120	0.0069580	0.50590
150	0.0066293	0.50517
180	0.0069777	0.50629
210	0.0065570	0.50628
240	0.0069647	0.50196
270	0.0071784	0.50850
300	0.0071764	0.50858

Stopping condition and tolerance

Numerical simulations related to temperature gradient and skin friction coefficient are captured against various indicated numerical values of parameters. An indigenous computer is running and simulate considered problem in view of iterative manner. It is noticed that exact solution of considered problem is not available. Therefore, stooping condition is illustrated by an error . Here, is very small which is and is almost equal to . Numerical values of Nusselt number and skin friction coefficient are noted when above condition is satisfied. Errors related various parameters are also simulated by Table 4.

Table 4

Numerical values of − and − versus and [21] when

		Nanoparticles					Hybrid nanoparticles
		(Re)^1/2C_f	Error	(Re)^1/2Nu	Error	(Re)^1/2C_f	Error	(Re)^1/2Nu	Error
γ	0.0	0.68347	6×10^-8	0.20378	7×10^-8	5.1854	5×10^-8	3.50901	7×10^-8
	0.4	0.78369	6×10^-8	0.20132	6.6×10^-8	5.2147	4.7×10^-8	3.30984	6.04×10^-8
	0.8	0.79257	5.7×10^-8	0.10463	5.2×10^-8	5.5147	4.1×10^-8	3.21069	5.31×10^-8
Fr	0.0	0.74058	5×10^-8	0.14246	5.1×10^-8	5.3441	3.7×10^-8	3.14617	4.3×10^-8
	0.7	0.77769	4.6×10^-8	0.14438	3.9×10^-8	5.5152	3.1×10^-8	3.14914	3×10^-8
	1.3	0.79600	4×10^-8	0.14680	2.3×10^-8	5.8441	2.5×10^-8	3.15157	2.73×10^-8
Hs	− 1.3	0.70793	3.7×10^-8	0.55061	2×10^-8	2.3238	2.2×10^-8	3.42032	2×10^-8
	0.3	0.39701	3.3×10^-8	0.33040	1.8×10^-8	2.2411	1.7×10^-8	3.20134	1.74×10^-8
	0.7	0.19601	3×10^-9	0.14039	1.2×10^-8	2.1953	1×10^-8	3.11143	1.1×10^-8
δ	0.0	0.08757	2.7×10^-9	0.17348	1×10^-8	2.4216	0.8×10^-9	3.17887	0.93×10^-9
	0.2	0.08212	1×10^-9	0.20652	0.91×10^-9	2.2149	0.6×10^-9	3.01252	0.5×10^-9
	0.5	0.08110	0.9×10^-9	0.14790	0.09×10^-8	2.3841	0.2×10^-9	3.00199	0.23×10^-9

Results and discussion

The constructing model including various characteristics is observed past a heated vertical porous plate. A hybrid approach along with non-Fourier’s theory inserting heat generation number provides complex model. Such complex model is solved by finite element method. It is mentioned that mixture of and is known as hybrid nanofluid and is investigated as a nanofluid. Graphs related velocity and thermal energy are plotted among comparison nanofluid and hybrid nanoparticles versus physical parameters. The graphical impacts are given below.

Visualization of flow characteristics via change in physical parameters

Measurement of velocity distribution is visualized versus change in Deborah number, Forchheimer number and porosity number. Related impacts of flow behaviors are observed by Figs. 2, 3 and 4. Distribution of fluid particles against the variation in Forchheimer number and porosity number is captured by Figs. 2 and 3. The concept of parameters related to and is formulated in momentum equation due to Forchheimer porous effect on flow and resistive force (Darcy’s porous). It is noticed that velocity has linear relationship against Darcy’s while Forchheimer has nonlinear relation against flow. Phenomena of fluid particles experience resistive force during the flow for both cases of Forchheimer and Darcy’s porous. Hence, retardation force is investigated during flow generating by Forchheimer and Darcy’s porous. Moreover, layers related to MB (momentum boundary) are decreasing function versus variation in and Flow is generated by hybrid nanoparticles are higher than for the case of nanostructures. Figure 4 reflects the role Deborah number on the velocity curves including hybrid nanomaterials and nano-structures. The motion of nanoparticles and hybrid nanomaterials is reduced versus the change in . This decrement behavior of flow is generated because of elastic nature (Maxwell fluid). Due to this nature of Maxwell liquid restores more deformation in view of flow analysis. MBLs (momentum boundary layers) are decreased due applying higher values of Deborah number. It is also investigated that approach related to hybrid nanoparticles are observed more efficient in view of maximum flow phenomena rather than approach related nanoparticles. Decreasing nature of flow characteristics is occurred due to elastic nature. It is also demonstrated that Maxwell liquid is more heated up as compared for the case of viscous fluid. Approach related to hybrid nanoparticles is observed more efficient to obtain maximum heat energy as compared for the case of approach related to nano-structures.

Figure 2

Velocity curves against when

Figure 3

Velocity curves against when

Figure 4

Velocity curves against when

Visualization of heat energy characteristics via change in physical parameters

The parameters called heat generation (), Deborah number (), unsteadiness () and time relaxation () are tested on temperature profile including approach of hybrid nanoparticles. The comparative simulations among nanoparticles and hybrid nano-structures are visualized by Figs. 5, 6, 7 and 8. The variation of heat energy against Deborah number is taken out by Fig. 5. Heat energy increases when time relaxation number is enhanced. The term related to relaxation time is arisen in dimensionless energy equation using concept of non-Fourier’s law. From Fig. 5, it is noticed that heat energy due to Fourier’s law of heat conduction is less than heat energy for the case of concept of non-Fourier’s law. An increment in heat energy is noticed when time relaxation number is increased. Figure 6 is plotted for measurement of heat energy with respect to large values of heat generation number. In this figure, positive values indicate heat generation while negative values reveal demonstrates heat absorption mechanism. Moreover, heating and cooling phenomena are indicated by positive and negative value, respectively. Maximum heat generation is made the reason to obtain maximum heat energy. External heat source is occurred due heat generation number. Thickness related to (MB) is inclined using higher values of heat generation number. A role of on temperature curves is estimated by Fig. 7. The concept of is formulate due to temperature difference in momentum equations. Temperature curves are enhanced using higher values of . Temperature curves for the case nano-structures are reduced temperature curves for hybrid nanoparticles. Hence, is investigated as a suitable parameter to generate more heat energy. Figure 8 illustrates the impact of unsteadiness parameter on heat energy profile including approach of hybrid nanoparticles. The unsteady flow is occurred due to unsteadiness parameter. Mathematically, it is the ratio of constant parameter (due to time) and parameter (due to stretching of wall). Flow becomes slow down by increasing the values of unsteadiness parameter because of inverse relation. Fluid is heated up applying large values of unsteadiness parameter. In comparative point of view, hybrid nanoparticles are observed more efficient to obtain higher amount for heat energy.

Figure 5

Temperature curves against when

Figure 6

Temperature curves against when

Figure 7

Temperature curves against when

Figure 8

Temperature curves against when

Behavior of temperature gradient and surface force

Table 3 is prepared to study the impact of Prandtl number on heat transfer rate. From the table, it is clear that rate of heat transfer enhances against Prandtl number. Furthermore, it is recorded that the obtained data is an excellent settlement with the studies reported in Ref.[39]. Distribution in temperature gradient (Nusselt number) and surface force (skin friction coefficient) is observed versus change in and . These impacts are simulated by Table 4. It is observed that reduction is found in temperature gradient using large values of and . Moreover, rate of heat transfer for the case of hybrid nanoparticles is higher than rate of heat transfers for nanofluid. Surface force increases versus the large values of and Further, approach related to hybrid nanoparticles are observed as significant in view of surface force as compared for the case of nanofluid.

Table 3

Validation of results regarding temperature gradient when

	Bilal et al.[39]	Present results
Pr
0.20	0.1693	0.1688371
0.70	0.4540	0.4541030
2.00	0.9114	0.9113893

Consequences of current problem

Characteristics of heat energy in Maxwell fluid suspending by hybrid nanoparticles pas melting vertical surface are modeled. Theories related to non-Fourier’s law and Darcy’s Forchheimer law is visualized in heat transfer phenomena along with heat generation. Finite element method is used to know results of current study. The main findings are listed below:

Thermal growth is enhanced for the case of hybrid nano-structures rather than for case of nanofluid.

Maximum thermal energy is achieved when heat generation number is increased.

Convergence analysis of considered problem is confirmed up to 300 elements.

Heat energy and TBL (thermal boundary layer) are augmented against large values of unsteadiness number and relation time parameter.

Motion into hybrid nanoparticles and nanoparticles becomes slow down versus higher values of Forchheimer and Darcy’s porous numbers.

Thickness of MBL (momentum boundary layer) is dominated for the case of hybrid nanoparticles rather than case of nanoparticles.

Minimum temperature gradient is measured versus variation in time relaxation and fluid number while skin friction coefficient is increased against variation heat generation and Forchheimer numbers.