The fluids that differ from Newtonian fluids in behavior and characteristics in the sense of not obeying Newtonian’s law are termed as non-Newtonian fluids, which include honey, paste, ketchup, and grease lubricant. There are many applications of the non-Newtonian fluid flow in modern industries and technology such as printing technology, biological solution, polymer, braking and damping devices, production of foods, and reduction agents in dragging. These fluids are considered to be the most effective in heat transmission phenomena (Ogunseye, Salawu, Tijani, Riliwan, Sibanda; Salawu and Ogunseye, 2020). Sharma and Shaw (2022) calculated the nanofluid flow over an expanding surface by assuming viscous dissipation and non-linear radiation and have concluded that the drag force has been augmented by an upsurge in the magnetic factor. Kumar and Sahu (2022) inspected the non-Newtonian fluid flow past an elliptical rotary cylinder through a laminar flow stream and have investigated the flow phenomenon numerically. Khalil et al. (2022) inspected the influences of fluctuating fluid properties of the double-diffusive model over the dissipated non-Newtonian liquid flow on a stretched surface. Sneha et al. (2022) appraised the magnetohydrodynamic (MHD) radiated nanoliquid flow toward a stretchy and shrinking sheet subject to the impact of carbon nanotubes and concluded that the velocity of the fluid declined, while the temperature had an upsurge with growing values of the magnetic parameter. Hu et al. (2022) used non-Newtonian fluids in a square channel to discuss the polydispersal, migration, and formation chain of particles. Islam et al. (2020) inspected the impacts of the MHD radiated micropolar fluid flow in a channel with the influence of hybrid nanoparticles. Waini et al. (2022) investigated the thermally radiative flow across an extending sheet by using magnetic field effects.

The branch of science that deals with magnetic characteristics of electrically conducting materials is termed as magnetohydrodynamics (MHD). This field of science provides a basis for many scientific, industrial, and technological applications such as liquid metals, cooling systems for automobiles, cooling of electronic chips, and production of chemicals. Sohail et al. (2020) inspected the MHD Casson fluid flow and entropy production, subject to variable heat conductivities past a non-linear bi-directional stretched surface, and deduced that the upsurge values of the magnetic factor have supported concentration and thermal profiles. Reddy et al. (2022a) used the MHD fluid flow with a porous medium to use the influence of radiation, thermal, and velocity slips and highlighted that the width of the boundary layer weakened with the growth in slip and heat factor parameters. Mishra et al. (2022) numerically explored the Williamson MHD nanofluid flow, subject to variable viscosities over a wedge. Reddy et al. (2016) debated the effect of thermal radiation over the MHD nanoparticle-based liquid flow past an extending surface and compared their results with a fine agreement to those results established in the literature. Bejawada and Yanala (2021) inspected Soret and Dufour impacts upon the time-dependent MHD liquid flow past an inclined surface placed vertically. Reddy et al. (2022b) scrutinized the influence of different slip effects over the MHD liquid flow past a stretchy sheet subject to Soret and Dufour effects. Sandeep et al. (2022) discussed the influence of the non-linearly radiated MHD hybrid nanofluid fluid flow using a heat source and concluded that the fluid flow declined and the thermal flow had an upsurge with a growth in the magnetic factor. Sandeep and Ashwinkumar (2021) studied the impact of different nanoparticles’ shapes upon the MHD fluid flow over a thin movable needle. Ashwinkumar et al. (2021) explained a 2D MHD hybrid nanoparticle flow using two different geometries of a cone and plate and proved that the flow and temperature incrimination are more visible in the case of the plate than that in the case of the cone. Mabood et al. (2022) inspected the influence of the non-linearly radiated 3D time-based MHD hybrid nanofluid flow. Readers can further study about the impact of MHD on mass and heat transmission in Sulochana et al. (2018), Alshehri et al. (2021), Mabood et al. (2021), Bejawada et al. (2022), Kumar et al. (2022), and Nalivela et al. (2022).

The mass and thermal flow problems with the impact of chemical reactions play a pivotal role in numerous fluid flow models. They have captivated more consideration due to its widespread utilization in many engineering and natural phenomena such as refrigeration, aerodynamic extrusions, and human transpiration. Sharma and Mishra (2020) documented the MHD nanoliquid flow using an internal thermal source. Singh et al. (2021) numerically solved the flow of a liquid past an enlarging sheet with the impact of chemical reactions and concluded that an upsurge in the stretching factor declined the diffusivities of heat and mass. Khan et al. (2021) discussed the bioconvection micropolar nanoparticle flow past a thin needle subject to binary chemical reactions and highlighted that mass diffusion declined with an upsurge in the chemical reaction factor and Brownian motion. Kodi et al. (2022) inspected the MHD Casson nanofluid flow past a vertically placed permeable plate subject to the impact of thermal diffusivity and chemical reactions. Kumar and Sharma (2022) discussed the influences of Stefan blowing on a fluid flow past a rotary disc subject to chemical reactions. Raghunath et al. (2022) inspected the time-dependent MHD flow of a liquid over an inclined permeable plate using magnetic impacts and chemical reactions.

Brownian motion and thermophoretic effects are responsible for controlling mass and thermal diffusivities subject to concentration and temperature gradients. Both these effects have numerous applications in different areas of science such as aerosol technology, nuclear safety phenomena, atmospheric pollution, aerospace technology, and hydrodynamics. Irfan (2021) considered the collective influence of Brownian and thermal diffusivity over the nanoparticle flow past a sheet with varying thickness, subject to slip conditions, and concluded that the augmentation of the Brownian number and thermophoresis factor has an upsurge in thermal profiles. Upreti et al. (2022) described the Casson fluid flow past a Riga plate subject to the effects of microorganisms. Saleem et al. (2022) studied the motion of water carrying three different types of nanoparticles subject to thermophoretic effects and the Brownian motion. Kiyani et al. (2022) inspected the MHD micropolar nanoparticle flow past an exponentially radiated surface using thermal radiations, Brownian motion, and thermophoretic effects upon the fluid flow system. Mehta and Kataria (2022) inspected the MHD fluid flow through the shrinking surface along with thermal radiations. Tayyab et al. (2022) have numerically studied the three-dimensional rotary nanoliquid flow subject to bio-convective activation energy.

The boundary-layer flows of nanofluids caused by stretched surfaces have attracted researchers’ attention recently (Andersson et al., 1994; Xu and Liao, 2009; Prasad et al., 2012; Mukhopadhyay, 2013). Their enormous significance in engineering and industrial applications has been the key driver behind this. These uses are particularly common in extrusion operations, paper and glass fiber manufacture, electronic chip manufacturing, paint application, food preparation, and the transfer of biological fluids. There is not a single constitutive connection between stress and the rate of strain that can be used to investigate all non-Newtonian fluids. The diversity of these fluids, their constitutive behavior, and simultaneous viscous and elastic properties make it nearly impossible to distinguish between effects resulting from a fluid’s shear-dependent viscosity and effects resulting from the fluid’s elasticity. A few mathematical models have been explained that closely match the experimental findings (Wu and Thompson, 1996). The Maxwell model is utilized for relaxation time in some highly concentrated polymeric fluids.

In this work, the authors have considered to present a semi-analytical solution of the Maxwell fluid flow over a bi-directional stretching sheet. Additionally, the thermal convective and mass flux conditions are taken into consideration. The mathematical framework of the existing problem is constructed on highly non-linear PDEs. Suitable similarity transformations are used for the conversion of PDEs into ODEs, which is presented in section 2. The flow problem is tackled with a homotopy analysis method, which is capable of solving higher-order non-linear differential equations, presented in section 3. The convergence of the HAM technique is also shown in section 4. Different flow profiles against various flow parameters are discussed physically, as shown in section 5. Finally, the concluding remarks are presented in section 6.

We consider the steady, laminar, and incompressible three-dimensional flow of a Maxwell fluid over a bi-directional linearly extending surface. The surface stretches along x and y directions and with velocity v w x = b y , where both a and b are constants. A magnetic field of strength B 0 is applied normal to the fluid flow. The surface temperature is denoted by T w , T f represents the reference temperature, and T ∞ shows the ambient temperature. In addition, the surface concentration is denoted by C w and C ∞ , showing the ambient concentration. Brownian motion, thermophoresis, and chemical reaction impacts are considered in this analysis. Additionally, the thermal convective and mass flux conditions are taken into consideration, as shown in Figure 1. Under the aforementioned suppositions, the principle equations are as follows (Bilal Ashraf et al., 2016; Dawar et al., 2021): ∂ u ∂ x + ∂ v ∂ y + ∂ w ∂ z = 0 , (1) u ∂ u ∂ x + w ∂ u ∂ z + v ∂ u ∂ y = ν ∂ 2 u ∂ z 2 − σ B 0 2 ρ u + λ w ∂ u ∂ z − λ u 2 ∂ 2 u ∂ x 2 + w 2 ∂ 2 u ∂ z 2 + v 2 ∂ 2 u ∂ y 2 + 2 u v ∂ 2 u ∂ x ∂ y + 2 u w ∂ 2 u ∂ x ∂ z + 2 v w ∂ 2 u ∂ y ∂ z , (2) u ∂ v ∂ x + w ∂ v ∂ z + v ∂ v ∂ y = ν ∂ 2 v ∂ z 2 − σ B 0 2 ρ v + λ w ∂ v ∂ z − λ u 2 ∂ 2 v ∂ x 2 + w 2 ∂ 2 v ∂ z 2 + v 2 ∂ 2 v ∂ y 2 + 2 u v ∂ 2 v ∂ x ∂ y + 2 v w ∂ 2 v ∂ y ∂ z + 2 u w ∂ 2 v ∂ x ∂ z , (3) u ∂ T ∂ x + w ∂ T ∂ z + v ∂ T ∂ y = k ρ C p f ∂ 2 T ∂ z 2 + ρ C p p ρ C p f D B ∂ T ∂ z ∂ C ∂ z + D T T ∞ ∂ C ∂ z 2 , (4) u ∂ C ∂ x + w ∂ C ∂ z + v ∂ C ∂ y + k 1 C − C ∞ = D B ∂ 2 C ∂ z 2 D T T ∞ ∂ 2 T ∂ z 2 , (5) with boundary conditions, given as follows (Bilal Ashraf et al., 2016; Dawar et al., 2021): u = a x , v = b y , w = 0 , D B ∂ C ∂ z + D T T ∞ ∂ T ∂ z = 0 , − k ∂ T ∂ z = − h T f − T , a t   z = 0 , v → 0 , u → 0 , C → C ∞ , T → T ∞   a s   z → ∞ . (6)

In the aforementioned equations, u , v , a n d   w are the velocity components; x , y , a n d   z are the coordinate axes; σ is the electrical conductivity; ρ is the density; λ is the relaxation time; B 0 is the magnetic field strength; C p is the specific heat; k is the thermal conductivity; D B is the Brownian diffusion coefficient; D T is the thermophoretic coefficient; T is the temperature; C is the concentration; k 1 is the chemical reaction coefficient; h is the heat transfer coefficient; and a and b are the velocity constants.

The similarity transformations are defined as follows (Bilal Ashraf et al., 2016; Dawar et al., 2021): v = a y g ′ η , u = a x f ′ η , w = − a ν g η + f η ,   θ η = T − T ∞ T f − T ∞ , ϕ η = C − C ∞ C w − C ∞ = , η = z a ν . (7)

The leading equations are transformed by using similarity transformations defined in (7): f ‴ η + 1 + M β f η f ″ η + g η f ″ η − f ′ η 2 − M f ′ η + β   2 f η f ′ η f ″ η + 2 g η f ′ η f ″ η − f ‴ η f η 2 − f ‴ η g η 2 − 2 f ‴ η g η f η = 0 , (8) g ‴ η + 1 + M β f η g ″ η + g η g ″ η − g ′ η 2 − M g ′ η + β   2 g ″ η f η g ′ η + 2 g ″ η g η g ′ η − g ‴ η f η 2 − g ‴ η g η 2 − 2 g ‴ η f η g η = 0 , (9) 1 Pr θ ″ η + g η θ ′ η + f η θ ′ η + N b θ ′ η ϕ ′ η + N t θ ′ η 2 = 0 , (10) ϕ ″ η + S c f η ϕ ′ η + S c g η ϕ ′ η + N t N b θ ″ η − S c K ϕ η = 0 , (11) with boundary conditions given as follows: f 0 = 0 , f ′ 0 = 1 , g ′ 0 = α , g 0 = 0 , θ ′ 0 = γ θ 0 − 1 , N t θ ′ 0 + N b ϕ ′ 0 = 0 , g ′ ∞ = 0 , f ′ ∞ = 0 , θ ∞ = 0 , ϕ ∞ = 0 . (12)

The embedded parameters are discussed as follows: M = σ B 0 2 ρ a , β = λ a , S c = ν D B , α = b a , Pr ⁡ = μ C p f k , K = k 1 a , γ = h f k ν a , N t = ρ C p p D B T w − T ∞ ρ C p f ν T ∞ , N b = ρ C p p D B C w − C ∞ ρ C p f ν . (13)

Here, M is the magnetic factor, β is the Deborah number, S c is the Schmidt number, α is the stretching constant, Pr is the Prandtl number, K is the chemical reaction factor, γ is the thermal Biot number, N t is the thermophoresis factor, and N b is the Brownian motion factor.

The Nusselt, Sherwood, and density numbers are defined as follows: N u x R e x = − θ ′ η η = 0 , S h x R e x = − ϕ ′ η η = 0 , (14) where R e x = x u w x ν is the local Reynolds number.

For an analytical simulation of the existing model, the HAM technique is considered. The initial guesses are given as follows: f 0 η = 1 − e − η , g η = α 1 − e − η , θ 0 η = γ γ + 1 e − η , ϕ 0 η = − γ 1 + γ N t N b e − η . (15) The linear operators are taken as follows: L f η = f ‴ − f ′ , L g η = g ‴ − g ′ , L θ η = θ ″ − θ , L ϕ η = ϕ ″ − ϕ , (16) with the following properties: L f ℧ 1 + ℧ 2 exp η + ℧ 3 exp − η = 0 , L g ℧ 4 + ℧ 5 exp η + ℧ 6 exp − η = 0 , L θ ℧ 7 exp η + ℧ 8 exp − η = 0 , L ϕ ℧ 9 exp η + ℧ 10 exp − η = 0 , (17) where ℧ i i = 1,2,3 , … , 10 are the constants.

R ∈ 0   1 shows the entrenching factor, and ℏ shows the auxiliary parameter. Then, the zero-order problems are constructed as follows: 1 − A L f f η ; A − f 0 η = A ℏ f N f f η ; A , g η ; A , (18) 1 − A L g g η ; A − g 0 η = A ℏ g N g g η ; A , f η ; A , (19) 1 − A L θ θ η ; A − θ 0 η = A ℏ θ N θ θ η ; A , f η ; A , g η ; A , ϕ η ; A , (20) 1 − A L ϕ ϕ η ; A − ϕ 0 η = A ℏ ϕ N ϕ ϕ η ; A , f η ; A , g η ; A , θ η ; A , (21) f 0 ; A = 0 , g 0 ; A = 0 , f ′ 0 ; A = 1 , g ′ 0 ; A = α , θ ′ 0 ; A = γ 0 ; A θ − 1 , N t θ ′ 0 ; A + N b ϕ ′ 0 ; A = 0 , f ′ ∞ ; A → 0 , θ ∞ ; A → 0 , g ′ ∞ ; A → 0 , ϕ ∞ ; A , (22) N f f η ; A , g η ; A = ∂ 3 f η ; A ∂ η 3 + 1 + M β f η ; A ∂ 2 f η ; A ∂ η 2 + g η ; A ∂ 2 f η ; A ∂ η 2 − M ∂ f η ; A ∂ η − ∂ f η ; A ∂ η 2 + β 2 f η ∂ f η ; A ∂ η ∂ 2 f η ; A ∂ η 2 + 2 g η ; Α ∂ f η ; A ∂ η ∂ 2 f η ; A ∂ η 2 − f η ; A 2 ∂ 3 f η ; A ∂ η 3 − g η ; A 2 ∂ 3 f η ; A ∂ η 3 − 2 g η ; A f η ; A ∂ 3 f η ; A ∂ η 3 , (23) N g g η ; A , f η ; A = ∂ 3 g η ; A ∂ η 3 + 1 + M β f η ; A ∂ 2 g η ; A ∂ η 2 + g η ; A ∂ 2 g η ; A ∂ η 2 − ∂ g η ; A ∂ η 2 − M ∂ g η ; A ∂ η + β 2 f η ; A ∂ g η ; A ∂ η ∂ 2 g η ; A ∂ η 2 + 2 g η ; A ∂ g η ; A ∂ η ∂ 2 g η ; A ∂ η 2 − f η ; A 2 ∂ 3 g η ; A ∂ η 3 − g η ; A 2 ∂ 3 g η ; A ∂ η 3 − 2 f η ; A g η ; A ∂ 3 g η ; A ∂ η 3 , (24) N θ θ η ; A , f η ; A , g η ; A , ϕ η ; A = 1 Pr ∂ 2 θ η ; A ∂ η 2 + g η ; A ∂ θ η ; A ∂ η + f η ; A ∂ θ η ; A ∂ η + N b ∂ θ η ; A ∂ η ∂ ϕ η ; A ∂ η + N t ∂ θ η ; A ∂ η 2 , (25) N ϕ ϕ η ; A , f η ; A , g η ; A , θ η ; A = ∂ 2 ϕ η ; A ∂ η 2 + S c f η ; A ∂ ϕ η ; A ∂ η + S c g η ; A ∂ ϕ η ; A ∂ η + N t N b ∂ 2 θ η ; A ∂ η 2 − S c K ϕ η ; A . (26)

For Α = 0 and Α = 1 , we obtain the following: f η ; 0 = f 0 η , f η ; 1 = f η g η ; 0 = g 0 η , g η ; 1 = g η θ η ; 0 = θ 0 η , θ η ; 1 = θ η ϕ η ; 0 = ϕ 0 η , ϕ η ; 1 = ϕ η . (27)

Using the Taylor series, we obtain the following: f η ; Α = f 0 η + ∑ Χ = 1 ∞ f Χ η R Χ   a s   f Χ η = 1 Χ ! ∂ Χ f η ; A ∂ A Χ Α = 0 , (28) g η ; Α = g 0 η + ∑ Χ = 1 ∞ g Χ η R Χ   a s   g Χ η = 1 Χ ! ∂ Χ g η ; A ∂ A Χ Α = 0 . (29) θ η ; Α = θ 0 η + ∑ Χ = 1 ∞ θ Χ η R Χ   a s   θ Χ η = 1 Χ ! ∂ Χ θ η ; A ∂ A Χ Α = 0 . (30) ϕ η ; Α = ϕ 0 η + ∑ Χ = 1 ∞ ϕ Χ η R Χ   a s   ϕ Χ η = 1 Χ ! ∂ Χ ϕ η ; A ∂ A Χ Α = 0 . (31)

The Χ th-order deformation problems can be written as follows: L f f Χ η − λ Χ f Χ − 1 η = ℏ f R Χ f η , (32) L g g Χ η − λ Χ g Χ − 1 η = ℏ g R Χ g η , (33) L θ θ Χ η − λ Χ θ Χ − 1 η = ℏ θ R Χ θ η , (34) L ϕ ϕ Χ η − λ Χ ϕ Χ − 1 η = ℏ ϕ R Χ ϕ η , (35) f Χ 0 = f Χ ′ 0 = f Χ ′ ∞ = 0 , g Χ 0 = g Χ ′ 0 = g Χ ′ ∞ = 0 , θ Χ ′ 0 − γ θ Χ 0 − 1 = θ Χ ∞ = 0 , N b ϕ Χ ′ 0 + N t θ Χ ′ 0 = ϕ Χ ∞ = 0 , , (36) R Χ f η = f Χ − 1 ‴ η + 1 + M β ∑ n = 0 Χ − 1 f Χ − 1 − j η f Χ − 1 ″ η + ∑ n = 0 Χ − 1 g Χ − 1 − j η f Χ − 1 ″ η − ∑ n = 0 Χ − 1 f Χ − 1 − j ′ η f Χ − 1 ′ η − M f Χ − 1 ′ η + β 2 ∑ n = 0 Χ − 1 f Χ − 1 − j η ∑ l = 0 n f Χ − l ′ η ∑ p = 0 l f l − p ″ η + 2 ∑ n = 0 Χ − 1 g Χ − 1 − j η ∑ l = 0 n f Χ − l ′ η ∑ p = 0 l f l − p ″ η + ∑ n = 0 Χ − 1 f Χ − 1 − j η ∑ l = 0 n f Χ − l η ∑ p = 0 l f l − p ‴ η − ∑ n = 0 Χ − 1 g Χ − 1 − j η ∑ l = 0 n g Χ − l η ∑ p = 0 l f l − p ‴ η − 2 ∑ n = 0 Χ − 1 g Χ − 1 − j η ∑ l = 0 n f Χ − l η ∑ p = 0 l f l − p ‴ η , (37) R Χ g η = g Χ − 1 ‴ η + 1 + M β ∑ n = 0 Χ − 1 f Χ − 1 − j η g Χ − 1 ″ η + ∑ n = 0 Χ − 1 g Χ − 1 − j η g Χ − 1 ″ η − ∑ n = 0 Χ − 1 g Χ − 1 − j ′ η g Χ − 1 ′ η − M g Χ − 1 ′ η + β 2 ∑ n = 0 Χ − 1 f Χ − 1 − j η ∑ l = 0 n g Χ − l ′ η ∑ p = 0 l g l − p ″ η + 2 ∑ n = 0 Χ − 1 g Χ − 1 − j η ∑ l = 0 n g Χ − l ′ η ∑ p = 0 l g l − p ″ η + ∑ n = 0 Χ − 1 f Χ − 1 − j η ∑ l = 0 n f Χ − l η ∑ p = 0 l g l − p ‴ η − ∑ n = 0 Χ − 1 g Χ − 1 − j η ∑ l = 0 n g Χ − l η ∑ p = 0 l g l − p ‴ η − 2 ∑ n = 0 Χ − 1 g Χ − 1 − j η ∑ l = 0 n g Χ − l η ∑ p = 0 l g l − p ‴ η , (38) R Χ θ η = 1 Pr θ Χ − 1 ″ η + ∑ n = 0 Χ − 1 g Χ − 1 − j η θ Χ − 1 ′ η + ∑ n = 0 Χ − 1 f Χ − 1 − j η θ Χ − 1 ′ η + N b ∑ n = 0 Χ − 1 θ Χ − 1 − j ′ η ϕ Χ − 1 ′ η + N t ∑ n = 0 Χ − 1 θ Χ − 1 − j ′ η θ Χ − 1 ′ η , (39) R Χ ϕ η = 1 Pr ϕ Χ − 1 ″ η + S c ∑ n = 0 Χ − 1 f Χ − 1 − j η ϕ Χ − 1 ′ η + S c ∑ n = 0 Χ − 1 g Χ − 1 − j η ϕ Χ − 1 ′ η + N t N b θ Χ − 1 ″ + − S c K ϕ Χ − 1 η , (40) where λ Χ = 0 , Χ ≤ 1 , 1 , Χ > 1 . (41)

The factors ℏ f , ℏ g , ℏ θ , and ℏ ϕ are called auxiliary factors, regulating the homotopic convergence. At the 23rd order of approximations, convergence regions for primary velocity, secondary velocity, temperature, and concentration distributions are shown in Figure 2. The convergence area of f ″ 0 is − 2.1 ≤ ℏ f ≤ 0.0 , g ″ 0 is − 2.25 ≤ ℏ g ≤ 0.2 , θ ′ 0 is − 2.1 ≤ ℏ θ ≤ 0.0 , and ϕ ′ 0 is − 2.25 ≤ ℏ ϕ ≤ 0.2 .

The physical investigation of the flow of the Maxwell fluid with the occurrence of the magnetic effect past an extending surface is explored in this section. With the occurrence of the thermal and mass diffusivity, the role of heat and mass transport is analyzed. The HAM procedure is used for the simulation of the existing model. Impacts of various flow parameters on flow distributions of the nanofluid are computed and discussed. Table 1 shows the validation of the present results with the published results. Here, a close relation between both the results is found, and we can validate our present analysis. The influences of N b , N t , and γ on the Nusselt number R e x − 1 2 N u x are investigated in Table 2. Table 2 shows that by increasing N b , N t , and γ , R e x − 1 2 N u x also increases. In Table 3, the variation in Sherwood number R e x − 1 2 S h x versus flow constraints such as N b , N t , S c , and K is examined. In this analysis, it is observed that greater N b , S c , and K values decrease R e x − 1 2 S h x , while increasing N t augments R e x − 1 2 S h x . Figures 3, 4 display variations in the primary velocity distribution of the nanofluid via the increasing Maxwell fluid factor β and magnetic parameter M , respectively. It is detected that the primary velocity distribution declines with growing values of β . Incidentally, increasing values of β correspond to the higher viscosity of the fluid, which, consequently, reduces the velocity of the fluid flow. Thus, the velocity distribution declines with the increase in β . Additionally, β = 0.0 corresponds to the Newtonian fluid. Thereafter, it is found that the Newtonian fluid is less viscous than the non-Newtonian fluid. Figure 4 is drawn to determine the role of the primary velocity distribution against augmenting values of M . In this observation, a decreased performance in the primary velocity distribution is found for the increasing M value. With the increasing magnetic field, the Lorentz force shows a retarding behavior against the flow behavior. Therefore, the Lorentz force opposes the fluid motion, which, consequently, decreases the boundary layer thickness and velocity distribution of the fluid flow. Furthermore, the increase in the magnetic parameter shows an upsurge of frictional forces between particles of fluids. This explains why velocity distribution is lower for higher magnetic factors. The effects of β and M on secondary velocity distributions are analyzed in Figures 5, 6. Since, the surface stretches linearly along both x − and y − directions. Therefore, similar impacts of the Maxwell fluid parameter β and magnetic parameter M are also found along the secondary velocity distribution. Figures 7–9 are plotted for the assessment of nanoliquid temperatures against increasing values of N b , N t , and γ . The consequence of N b on an energy profile is shown in Figure 7. An increase in the nanofluid temperature is examined for expanding the values of N b . Brownian motion refers to the movement of particles; therefore, an increased production of heat occurs, which raises the energy profile. Figure 8 shows the effect of N t on the temperature distribution. Figure 8 describes that the nanoliquid temperature is enhanced for intensifying values of N t . In the case of thermophoresis, liquid elements are quickly transformed from the hot region to the cold region with a rising the thermophoresis parameter N t which, consequently, shows a surge in the temperature distribution. The outcome of the nanoliquid temperature for increasing values of the thermal Biot number γ is observed in Figure 9. In this figure, the increasing behavior in the temperature distribution due to γ is observed. For greater values of γ , the heat transfer coefficient is enhanced because the heat transfer coefficient is directly related to the thermal Biot number γ . Therefore, the temperature distribution of the nanoliquid increases for the higher thermal Biot number γ . Figures 10–13 are displayed to discuss variations in the nanoliquid concentration distribution with respect to expanding values of N b , N t , S c , and k r . The result of N b on the concentration distribution is shown in Figure 10. It is examined that, a rise in N b reduces the nanofluid concentration distribution. As N b increases, the concentration gradually falls. The explanation behind this is because higher values of the Brownian parameter enhance fluid particle collisions and lower the viscosity of nanofluids. Figure 11 explains the role of N t on the concentration distribution. In this figure, it is perceived that enhancing values of N t increase the nanofluid concentration distribution. This is because the thermodiffusion coefficient and the thermophoresis parameter are closely related. Increased diffusion coefficients are implied by higher values of N t , which intensify the concentration distribution. The impact of S c on the concentration distribution is discussed in Figure 12. It is observed that the nanofluid concentration is lower with the expansion of the Schmidt number S c . The Schmidt number S c is the ratio of momentum diffusivity and mass diffusivity. So, the mass diffusivity of the fluid decreases with the increase of the Schmidt number S c . This is because the Schmidt number and mass diffusivity are inversely related to each other. Thus, a decrease in mass diffusivity decreases the concentration distribution. The consequence of K on the concentration distribution is shown in Figure 13. It should be noted that the growth in the values of K decays the nanofluid concentration. Furthermore, it can be perceived that molecular diffusivity is lower for higher chemical reactions. Therefore, a lower molecular diffusivity decreases the boundary layer thickness and concentration distribution of the nanoliquid. Figures 14, 15 show streamline patterns of Newtonian and Maxwell fluids, respectively. It should be noted that the analysis streamlines for Maxwell and Newtonian models in the current study are indeed distinct from one another.

This article examines the 3D flow of a Maxwell nanofluid across a bi-directional stretching surface with magnetic field applications. In this approach, the effects of Brownian motion, thermophoresis, and chemical reactions are taken into account. The conditions for mass flux and thermal convection are also taken into account. The modeled problem is solved using the HAM technique. The HAM convergence is also demonstrated. Simulations and detailed discussions are carried out to determine the effects of various physical factors on flow profiles and quantities of interest. Key findings of the current problem are as follows:

• The magnetic and Maxwell fluid parameters determine the decreasing functions of primary and secondary velocities.