Abstract

Little is known in the literature about the concept of nonuniform heat source/sink and higher-order chemical reaction for the dynamics of Oldroyd-B nanoparticles. Therefore, the present article addresses the nonuniform heat source/sink and higher-order chemical reaction features in nonlinear mixed convection bidirectional MHD dynamics of Oldroyd-B nanoparticles with thermal radiation aspects through porous space. Stratification effects for both the temperature and concentration setups are also used in the mathematical model with the significance of random movement and thermodiffusion of nanoparticles. Shape-preserving transformations have been employed to convert the transport equations into solvable forms. An innovative analytical tactic, namely, homotopy analysis method, has been adopted to find the solution of the modeled problem. Behaviors of pertinent parameters on thermal and concentration profiles have been discussed through various graphs. Inspection of heat/mass transport against appropriate varieties of pertinent parameters has been made and explained physically. Thermal profile is augmented with the higher estimations of space and temperature-dependent heat source/sink links. Concentration profile is diminished with the augmentation of higher-order chemical reaction parameter. Sherwood number is improved with the estimation of and is reduced with the growth of . Nusselt number is declined with the upgraded amounts of and .

1. Introduction

Rheology of non-Newtonian fluids is based on their multiphase nature, and they have got great importance on account of their wide spread involvement in industrial and technological applications. Shear-thinning, shear thickening, and viscoelasticity are the main characteristics of non-Newtonian liquids. These characteristics are described by nonlinear relationship between shear force and shearing deformation. Several models of non-Newtonian materials have been suggested, like Maxwell, Jeffery, and Oldroyd-B models. An Oldroyd-B model is a rate type non-Newtonian material having relaxation and retardation time effects.

Thermal energy has a significant role in almost every field of science, engineering, industries, biomedicine, plant processing, transports, power houses, generating stations, and many others. The effectiveness of all the above cited applications is greatly dependent upon thermal conductivity of the fluids involved in the heat interchanging processes. Many fluids like water, ethylene, and oil bear the restrictions in applications due to low thermal conductivity. Hence, nanoparticles of average size of about 1–100 nm are inserted into the base fluid to augment the thermophysical and concentration circumstances of the fluid, and the obtained colloidal suspension is regarded as nanofluid. The novel attributes of nanofluids make them positively functional in the numerous applications of heat transference along with microelectronics, energy cells, and domiciliary refrigeration and heating equipment, in disintegrating and cracking units, machining, hybrid powered engines, medications, and heat exchange units. Choi and Eastman [1] firstly proposed the term nanofluid through a permeable prolonging sheet along with consumption/injection. Then, Xuan and Li [2] extended the concept of nanofluids by proposing that fluids with solidified particles of metalloid and polymeric origin in a base fluid can also be termed as nanofluids. Buongiorno [3] formulated a mathematical tool for the analysis of convective transference in nanofluids with the consideration of Brownian movement and thermal diffusion properties. Copper, titanium, silicon, aluminum, zinc, magnesium, and graphene oxide are the best well-known particulates used for the development of nanofluids, while water, oil, glycol, and ethylene are the most regular working liquids.

Boundary-layer flows through stretching surfaces have been considered in studies because of their countless applications from household practices to aerospace such as in industries, irrigation waterways, environmental production, aerodynamics, sports, enhancement of heat transference, enhancement of mixing, transport of species, formation of rubber and plastic sheets, and fiber glass manufacturing. Boundary-layer flow produced by continuous elongation of obstacle was primarily studied by Sakiadis [4]. Crane [5] introduced the development of boundary-layer stream through stretching of surface and found exact solutions for the flow field. The bidirectional flow caused as a result of linearly plane bidirectional stretching of sheet was studied by Wang [6]. Ariel [7] deduced approximate analytical and numerical solutions of bidirectional steady flow over an elongating sheet. Sajid and Hayat [8] premeditated the aspects of thermal deposition and heat transportation in the boundary-layer dynamics of a viscid fluid influenced by an exponential elongation of obstacle. An analysis of the motion of a nanofluid produced by stretching of surface using connective model was developed by Khan and Pop [9].

Bhattacharyya and Layek [10] investigated the effects of chemically receptive solute convergence in magnetohydrodynamic boundary-layer stream. Ahmad et al. [11] continued the work accomplished by Liu and Anderson [12] by estimating the Darcy resistance impact and applied magnetic field. Stratification is a predominant phenomenon which has gained remarkable consideration because of its inclusion in the geophysical flows such as waterways, within lakes, oceans, ground-water supplies, warm energy stockpiling frameworks, and so on. A lot of researchers experimented with stratification with different effects. Ibrahim and Makinde [13] considered the doubly stratified boundary-regime flow of nanomaterials through an upright plate. Loganathan et al. [14] inspected the effect of second-order slip in a convective dynamic of an Oldroyd-B material with the aspects of chemical reaction and thermal radiation. Sandeep and Reddy [15] studied magnetohydrodynamic (MHD) flow having double stratified and cross diffusion aspects for an Oldroyd-B fluid. Waqas et al. [16] described the features of mixed convection progression of an Oldroyd-B nanomaterial via accumulating the impact of heat creation as well as heat/mass stratification.

The notion of heat production/consumption is convenient in many practices including heat disposal of atomic fuel wreckage, underground removal of radiative garbage substances, food storages, and disconnecting materials in packed bed reactors. In this way, nanocomposites have a capacity to assemble incident radiation. Turkyilmazoglu and Pop [17] explored the concurrent effects of heat and mass transport with thermal radiation through an unsteady nanomaterial progression. Rashidi et al. [18] analyzed the effect of second law of thermodynamics in magnetohydrodynamic dynamics of nanoliquid over a revolving disk. Moradi et al. [19] performed an analysis to check the effects of nanoparticles in Jeffery–Hamel flow. Makinde et al. [20] analyzed the buoyancy-driven dynamics of nanoliquid near a stagnated domain with connective boundary conditions. Mahanthesh et al. [21] worked toward approximate solution for the dynamics of an Oldroyd-B material with the significance of thermal deposition and heat consumption/generation impacting upon a nonlinearly expanding obstacle.

The subject of fluid motion with the involvement of porous medium has developed significantly in many fields of science and engineering like flow through water rocks, skin pronouncement, liquid purification processes, chemical trash, crude oil production, porous insulation, grain storage, mastic transport modeling, and underground removal of atomic waste. The interdependence of the voids in a porous medium permits certain liquids to flow through material. Ghosh and Sana [22] studied the time dependent dynamics of an Oldroyd-B material influenced by reformed sine pulses. Mukhopadhyay [23] considered and studied the time dependent heat transference and mixed connective dynamics produced by a penetrable stretching obstacle. The magnetohydrodynamic dynamics of tiny-sized nanoparticles passed over a permeable wedge is scrutinized by Kandasamy et al. [24]. Tripathy et al. [25] analyzed the impact of chemical reaction on magnetohydrodynamic free connective sheet over an upright moving plate through permeable material. Ali et al. [26] presented the analytical solution for flowing fluid through exponential stretching of permeable surface including heat flux in permeable medium utilizing homotopy analysis method (HAM). Goud [27] introduced the magnetohydrodynamic flow through an upright appendicular plate including radiation and chemical reaction effects through a permeable medium. Some more impactful scientific contributions related to the bidirectional dynamics of nanofluids with various mechanical and thermal aspects have been made by some researchers [28–36].

The above literature survey reveals that the study of non-Newtonian nanofluid performs a vital role in the industrial and engineering advancements. Moreover, a huge gap regarding the combined study of nonlinear aspects of mixed convection, chemical reaction, and heat consumption/generation is found during literature review. Hence, the novelty of present exploration is to discuss the steady dynamics of 3D radiative non-Newtonian (Oldroyd-B) nanofluid with the nonlinear aspects of mixed convection, chemical reaction, and heat source/sink. Additionally, significance of porous medium, magnetohydrodynamics, and double stratification has been addressed to achieve the thermal engineering relevancy of the mathematical model. To the best of our knowledge, no such attempt has been made previously. Mathematical modeling of the physical model has been completed with the help of boundary-layer approximations and basic laws of fluid dynamics. With the functionality of similarity transformations, transport equations have been converted into one parameter family of solvable equations based on the physical domain. Analytical analysis of the mathematical model has been presented by following the procedure of homotopy analysis method [37–44]. Some more solution techniques (Keller-Box method, Lie group analysis, etc.) have also been offered by some researchers/scholars [45–50], but here preference is given to homotopy analysis method because of its compatibility with the modeled problem. Finally, the obtained outcomes have been discussed physically with the support of various charts and tables against the fluctuating choices of involved parameters.

2. Model Development

An incompressible, doubly stratified, and three-dimensional steady dynamic of radiative Oldroyd-B magneto-nanomaterial impacting upon an expanding obstacle through porous space is considered with the nonlinear aspects of mixed convection, chemical reaction, and heat consumption/generation. Effects of random motion and thermodiffusion of nanoparticles have also been provided via Buongiorno nanofluid model. The Lorentz force is perpendicular to the -plane, whereas gravitational force is parallel to the -direction. Influence of chemical reaction has been included in the nanoparticle mass transport equation, whereas radiation effects have been involved in the heat equation of the nanofluid. The expansion velocity along -direction is explained as , whereas expansion velocity along the -direction is elaborated as . The relations for wall and free stream temperatures are settled according to the rule of thermal stratification, whereas the relations for wall and free stream nanoparticle mass concentration are established according to the regulation of solutal stratification. The graphical abstract of the present investigation is presented via Figure 1. The governing equations for the present analysis supported by [16, 22, 27, 31, 32, 44] are composed as

Here, .

The relevant boundary conditions supported by [15, 16, 44] are

Using the relation of Rosseland approximation, radiative heat-flux is determined as (see [48]) with .

Here, , and describe the space coordinates; , , and express parts of velocity in the direction of , , and , respectively; and are used to represent relaxation time and retardation time, respectively; denotes kinematic viscosity with as dynamic viscosity and as density; represents the gravitational acceleration; exhibits the involvement of linear thermal expansion; represent involvement of nonlinear thermal expansion; and depict, respectively, linear and nonlinear concentration coefficients; denotes Stefan–Boltzmann constant; stands for magnetic field strength with electrical conductivity ; is used to express the porosity constant; the permeability constant; is the temperature, is the concentration of nanoparticles; denotes the thermal diffusivity with thermal conductivity ; and and are the symbols of random and thermodiffusion coefficients of nanocomposites, respectively. Efficient heat capability of nanoparticles and efficient heat capability of the fluid are symbolized by and , respectively; and represent the first- and second-order chemical reaction parameters, respectively; is thermal deposition; is nonuniform heat consumption/absorption; is temperature at the surface; denotes the ambient temperature; and and serve as nanoparticle concentrations at the obstacle and far from the obstacle. Here, the expansion velocities, temperature, and nanoparticle mass concentration at the surface are described as (see [15, 16, 44])

Here, are dimensional constants.

We utilize the accompanying similarity transformations supported by [37, 38, 44] as follows:

Via using the above-mentioned local similarity set, (1) is straightforwardly fulfilled and simplified forms of (2-5) are mentioned below:

The boundary conditions are confined as

In the above-mentioned equations, and are the dimensionless velocities along - and -direction, respectively; and are the dimensionless temperature and concentration, respectively; and are the Deborah numbers corresponding to relaxation and retardation times, respectively; and indicate the nonlinear convection constraints for temperature and concentration, respectively; is the magnetic constraint; is used to express the porosity constraint; stands for buoyancy ratio parameter; and denote the mixed convective parameters for thermal and concentration setups, respectively; is the Lewis number; signifies the Prandtl number; denotes the random movement involvement; corresponds to thermodiffusion parameter; and denote, respectively, the quantities of space-dependent and temperature-dependent heat consumption/generation; is the rate of chemical reaction; is the dimensionless chemical reaction parameter; represents radiation parameter; represents the Rosseland mean absorption constant; is the expansion ratio parameter; and are the representatives of the heat and mass stratified parameters, respectively; and prime represents differentiation taking as an independent similarity variable.

The amount of mass transference is presented by local Sherwood number while the amount of heat transference is described by local Nusselt number , and these numbers can be expressed in the form of Reynolds number as below:

3. Homotopic Series Solution

This section deals with the series solutions of the mathematical modeled dimensionless equations (9)–(12) based on boundary restrictions via (8) and (9) with the procedure of homotopy analysis method (HAM). For this purpose, primary guesses are approximated in accordance with boundary conditions (8) and (9) and are shown as

The linear operators are denoted by , , and as follows:

In the above relations, , , , and are auxiliary constraints, whereas , , , and serve as nonlinear operators, and is the embedding parameter. The BCs (boundary conditions) in the coding language can be expressed as

The general solutions , , , and of (23)–(26) with regard to particular solutions are expressed by

The general constants in the above equations are calculated by applying boundary conditions and are stated as

4. Convergence Analysis

The admissible range for unknown boundary restrictions is the first task to obtain the solution via homotopy analysis method. This range can be traced by identifying those portions of curves lying exactly parallel to the horizontal axis as sketched in Figure 2. The solid-green line curve expresses with admissible range of . The admissible range for the curve of is , and it is depicted by the black dashed line. Similarly, the admissible ranges of and are estimated as and , and they are represented by red and blue curves, respectively. For convenience, throughout the analysis we have chosen . Hence, the convergence table for appropriate amounts of emerging parameters is generated, and the corresponding results are shown through Table 1. Table 1 illustrates that at least sixty-four orders of approximation are needed to obtain the missing boundary conditions for and , whereas at least seventy orders of approximation are required to attain the momentum boundary conditions. To maintain the convergence and stability of the solution, the numbers of approximations are increased up to eight.

5. Results and Discussion

This section describes the graphical illustrations of temperature fluctuation, concentration fluctuation, local Nusselt number, and local Sherwood number under the consequences of important involved parameters. Figure 3 explains the impact of thermal stratification parameter on dimensionless temperature fluctuation when the similarity variable is increased from to with the constant choices of other involved parameters. Maximum temperature fluctuation occurs at , minimum temperature fluctuation is achieved at , and thermal stratification phenomenon is generated for . Temperature profile is reduced with the improvement in the choice of . Physically, higher choice of develops the ambient temperature and reduces the temperature at the surface. As a result, the temperature profile is reduced. The thickness of thermal layer is improved with the improvement in the value of the similarity variable. Figure 4 illustrates the trend of solutal stratification parameter on dimensionless concentration fluctuation when the similarity variable is increased from to with the fixed amounts of other influential parameters. Maximum concentration is achieved at , minimum concentration is obtained at , and solutal stratification phenomenon is traced for . Concentration setup is diminished with the higher estimation of . The physical reason behind this outcome is the downfall of nanoparticle concentration at the wall and improvement of nanoparticle concentration away from the stretching device. Concentration layer thickness is also increased with the improvement of the similarity variable.

Figure 5 indicates the role of space-dependent heat source/sink parameter in dimensionless temperature fluctuation with the increase of the similarity variable from to . Minimum temperature is attained for heat consumption case, ; maximum temperature is achieved for heat generation case, ; and moderate temperature is reached in the absence of heat generation/consumption, . Mathematically, is involved as a multiple of velocity profile in (7); therefore, the velocity profile also increases with the higher estimation of , and hence temperature profile is improved. Therefore, it is called space-dependent heat source/sink influence on thermal profile. The thickness of thermal layer is also improved with the growth of the similarity variable. Figure 6 illuminates the potential of temperature-dependent heat source/sink parameter on dimensionless temperature fluctuation with the growth of the similarity variable from to . Minimal temperature is accomplished for heat ingestion case, ; maximal temperature is realized for heat production case, ; and modest temperature is touched in the nonappearance of heat production/consumption, . Precisely, is involved as a multiple of in (7); therefore, the thermal profile is also augmented with the greater estimation of , and hence it is called temperature-dependent heat source/sink influence on thermal profile. The depth of thermal layer is also enriched with the progress of similarity variable. Figure 7 establishes the connection between dimensionless concentration fluctuation and chemical reaction parameter with the growth of the similarity variable up to . Decreasing trend in concentration fluctuation is observed with the growth in the estimation of , maximal concentration is achieved for smaller values of , and minimal concentration is generated for greater estimations of . Physically, higher estimation of increases the second-order chemical reaction and reduces the first-order chemical reaction at the same time. As an outcome, negative trend in the concentration profile is attained. The thickness of concentration layer is also developed with the improvement in the similarity variable. Figure 8 explores the effect of thermal radiation parameter on thermal fluctuation with the evolution of the similarity variable from to . A huge improvement in the thermal fluctuation is attained with the intensification in . Maximal temperature is measured for higher , and minimal temperature is obtained for basic/lower . The thickness of thermal layer is also improved and it is found to be higher for the middle values of . Physically, radiation is basically the transference of energy in the form of electromagnetic waves. More electromagnetic waves are produced with the improvement of , and hence dominant temperature fluctuation is observed when reaches maximal estimation. Figure 9(a) deals with the relationship between relaxation time Deborah number and temperature distribution for the similarity variable . Here, thermal environment is upgraded with the higher estimation of . Basically, is generated in the momentum equation due to the inclusion of Oldroyd-B nanoparticles and defined as the product of expansion rate and the relaxation time . Hence, higher estimation of enhances the value of as well as the expansion rate. As an outcome, improvement in thermal setup is noticed with the upgraded width of thermal layer. Figure 9(b) shows the effect of on concentration fluctuation. Concentration is tremendously promoted with the advanced values of . Physically, upgraded value of enhances the expansion rate that increases the movement of Oldroyd-B nanoparticles, and so concentration setup is boosted. Figure 10(a) establishes the link between retardation time Deborah number and thermal fluctuation for the similarity variable . At this time, thermal environment is degraded with the refined approximation of . Fundamentally, involves in the momentum equation due to the presence of Oldroyd-B nanoparticles and is defined as the mathematical product of the stretching rate and the retardation time . Therefore, estimation of augments the value of as well as the expansion rate. As a conclusion, reduction in thermal setup is perceived with the progressive width of thermal stream. Figure 10(b) illustrates the influence of on concentration flux. Concentration is enormously decreased with the advancement of . Actually, improved value of develops the retardation rate, and so concentration setup is reduced.

Nusselt number measures the rate of heat transference whereas Sherwood number measures the rate of mass transference from hot region to cold region. Figure 11 describes the rate of heat transport with the combined inspiration of random motion and thermodiffusion parameters. Rate of heat transport is declined with the upgraded amounts of these parameters related to the nanofluid. Maximal rate of heat transport is measured for smaller amounts of these nanoparticle parameters, whereas minimal rate of heat transport is attained for greater amounts of and . Maximum value of Nusselt number is approximated as , whereas minimum value of Nusselt number is estimated as for the adopted ranges of Brownian and thermodiffusion parameters. Physically, a higher value of enhances the Brownian diffusion coefficient, and a higher value of enhances the thermodiffusion coefficient, whereas kinematic viscosity of Oldroyd-B nanofluid is reduced with the escalation of these parameters. As a whole, these physical changes are the main causes for the reduction of heat transport. Figure 12 shows the rate of mass transport with the collective stimulation of random motion and thermodiffusion parameters. Rate of mass transport is degenerated with the progress in the amount of , whereas it is improved with the advancement in the amount of . Highest rate of mass transport is stated for smaller choice of and higher choice of . Maximum value of Sherwood number is formulated as , whereas minimum value of Sherwood number is predicted as for the implemented ranges of Brownian and thermodiffusion factors. Physically, a higher value of enhances the heat capacity induced by nanoparticles, and a higher value of develops the temperature difference, whereas density of Oldroyd-B nanofluid is upgraded with the intensification of these forces. Overall, these physical variations are the main reasons for the improvement of Sherwood number with the positive tendency in and reduction of Sherwood number with the positive approach of . Figure 13(a) indicates the rate of heat transport with the communal effect of thermal convection parameter and concentration convection parameter . The maximum value of heat transport is reported as for smaller communal effect of convection parameters, whereas the minimum value of heat transport is calculated as for larger communal effect of convection parameters. Nusselt number is reduced with the higher estimation of , and it is improved with the growth of . On the other hand, Figure 13(b) explains the rate of mass transport with the combined effect of convection parameters. Sherwood number is improved with the estimation of , and it is reduced with the growth of . Maximum mass transport is observed at the right top of the figure and minimal mass transport is noticed at the left bottom of the figure. Physically, thermal convection parameter depends on the temperature difference, and thermal environment is reduced with its escalation; hence, rate of heat transport is improved. Moreover, concentration convection parameter relies on the concentration difference, and thermal environment is boosted with its growth; hence, rate of heat transport is reduced for this situation.

Table 2 is arranged to compare the present outcomes with the previously published work in the absence of porosity, magnetic environment, dual stratifications, radiation, chemical reaction, and space-dependent heat source/sink factors. The reduced results for thermal and concentration gradients are best matched with the previously published report by Khan et al. [47] via homotopy analysis method.

6. Conclusions

The significance of higher-order chemical reaction and nonuniform heat source/sink influence on the dynamics of Oldroyd-B nanoparticles have been analytically investigated in this contribution. For physical relevancy, thermal radiation, mixed convection, and double stratification have also been incorporated in the bidirectional flow of the nanoliquid with the help of Buongiorno nanofluid model. Chief findings/outcomes of the ongoing exploration are listed below:(i)Level of heat exposure is improved with the intensification in thermal stratification parameter , whereas level of mass exposure is reduced with the improvement in solutal stratification parameter .(ii)Level of mass exposure is developed with the higher selections of Brownian motion and chemical reaction parameters, whereas level of heat exposure is condensed with the escalation in the choices of heat source/sink parameters and radiation factor .(iii)Nusselt and Sherwood numbers are decreased with greater amounts of concentration convection parameter, whereas opposite trend is noticed for thermal convection parameter.(iv)Higher estimations of Brownian and thermodiffusion parameters reduce the rate of heat transport.(v)Temperature and concentration profiles are reduced with the positive tendency of retardation time Deborah number , and opposite trends are detected for relaxation time Deborah number .

Finally, the outcomes obtained through this investigation are helpful in the improvement of thermal instruments. This mathematical model has some applications in biomedicine and chemical and mechanical engineering like coating a sheet with nanomaterial, paper production process, and treatment of cancer.

Data Availability

The raw data supporting the conclusion of this paper will be made available by the corresponding author without undue reservation.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors contributed equally to this work and read and approved the final version of the paper.