In recent times, numerous fascinating energy applications have made use of nanotechnology to provide more effective and eco-friendly products and services. Moreover, one crucial aspect that is necessary to bring down the cost of maintenance or installation is the effective heat transport distribution across a power network. The optimal transmission of thermal qualities, such as density, specific heat, viscosity, and thermal conductivity, is essential for the optimum thermal efficiency. In this context, we consider the graphene nanoparticles to improve the thermal characteristics of the host fluid (CMC−water). In order to increase the rate of heat exchange, the flow patterns of nanofluids across thermal energy technologies are crucial. Several nanofluids have very engrossing applications in several fields, such as electronics cooling, automotive engine cooling, heat exchangers, boiling heat transport, solar collectors, medicine, and nuclear system cooling [1]. For the first time, [2] introduced the idea of nanofluids and empirically supported his concept. Gamachu and Ibrahim [3] investigated the mixed convection hybrid viscoelastic nanofluid. The heat exchange of the viscoelastic nanofluid past the stretched surface with many impacts was studied by [4]. In [5], a corrugated channel was considered to study the heat transport of the graphene nanofluid. [6] examined the influence of viscous dissipation on the viscoelastic nanofluid due to a circular cylinder. Numerical examination of the graphene nanofluid because of the stretchable wall was conducted by [7]. [8] directed the enhanced thermal progress of the graphene hybrid nanofluid past a stretched sheet. The impact of zero nanoparticle flux and hall current regarding the movement of the nanofluid above a stretched wall was discussed by [9]. Several other researchers studied the graphene nanofluid in [10–15]. Recently, [16] examined the MHD nanofluid past a convectively heated Riga plate placed horizontally, embedded in the porous (Darcy−Forchheimer) medium. The EMHD flow of a heterogeneous micropolar mixture with different concentrations of water, ethylene glycol, and copper oxide nanoparticles was discussed by [17]. Due to the distributions of the polymer/CNT structure nanocomposite material, morphological nanolayers have an influence on the movement of hybrid nanofluids which was studied by [18]. [19] analyzed the tetragonal nanoparticles with a variety of density and conductivity characteristics flowing in 3D using non-linear Boussinesq and Rosseland estimations. Recent research investigates how the presence of nanoparticles can significantly improve heat exchange in a variety of physical geometries [20–22]. [23] examined the analytical answers for a thermal conductor force peristaltic flow to temperature-dependent nanofluid viscosity. Utilizing two distinct methods for nanofluid analysis, the silver−water nanofluid electro-osmotic movement that is controlled by peristalsis was investigated by [24]. [25] studied the electro-osmotically assisted peristaltic pressurization of the MoS₂ Rabinowitsch nanofluid and produced randomness. The use of novel nanofluids throughout clinical isolates to fight Staphylococcus aureus was investigated by [26]. Many researchers talked about boundary layer problems in various media [27–29].

The investigation of the convective boundary layer flows across a stretched surface has a broad range of uses in the polymer production sector. For example, it is employed in manufacturing fibers, polymeric extruder, and fiber glass, sheet manufacturing, condensing of aqueous films, food production, growth of crystals, and fabrication of artificial films [4]. The idea of a boundary layer flow above a stretching sheet was promoted by [30]. [31] examined the velocity slip impact on the MHD nanofluid because of the stretched sheet. The study of the magnetite nanofluid provoked by a stretching wall was observed by [32]. [33] introduced the system disorder inside the nanofluid past a stretchable wall. The MHD nanofluid with two different nanoparticles due to a heated stretched wall was analyzed by [34]. [35] revealed the linearly stretched sheet-induced double-diffusive time-dependent magnetohydrodynamic flow of the nanofluid involving convection boundary constraints. Optimizing entropy production with an unstable stagnating Casson nanofluid flow having dual chemical changes across a stretched surface was discussed by [36]. The investigation of the nanofluid with many effects above a stretching sheet was presented in [37–39].

From the potential applications mentioned previously, we analyzed the MHD graphene−CMC−water nanofluid past a stretchable wall with Joule heating and velocity slip impact. The closed-form solutions were derived analytically. To understand the behavior of various physical parameters on the fluid temperature, Nusselt number, velocity, and skin friction, the numeric tables and graphs were portrayed.

In framing the model, the flow of an incompressible and steady MHD viscoelastic nanofluid (Walter’s liquid B type) across a stretchable sheet in two dimensions is investigated. The graphene and carboxymethyl cellulose (CMC) solution−water are taken as nanoparticles and host fluid, respectively. The fluid is made up of the region y > 0, where the x-axis runs alongside the stretched surface in a flow pattern and the y-axis runs perpendicular to the flow. With a velocity of u ̲ = c x + β ∂ u ̲ ∂ y , the surface is pulled in the x dimension. Furthermore, the magnetic field B ₀ is supplied normally toward the flowing fluid, and the consequences of thermal radiation are also utilized. The basic equations regulating the flow are as follows [6]: ∂ u ̲ ∂ x + ∂ v ̲ ∂ y = 0 , (1) u ̲ ∂ u ̲ ∂ x + v ̲ ∂ u ̲ ∂ y = μ gnf ρ gnf ∂ 2 u ̲ ∂ y 2 − k 0in u ̲ ∂ 3 u ̲ ∂ x ∂ y 2 + v ̲ ∂ 3 u ̲ ∂ y 3 + ∂ u ̲ ∂ x ∂ 2 u ̲ ∂ y 2 − ∂ 2 u ̲ ∂ x ∂ y ∂ u ̲ ∂ y − σ gnf B x 2 ρ gnf u ̲ . (2)

The problem’s appropriate boundary criteria are as follows [40]: u ̲ = c x + β ∂ u ̲ ∂ y , v ̲ = 0 at y = 0 , u ̲ = 0 as y → ∞ , (3) where c is the stretching rate, β is the slip factor, ν _gnf is the kinematic viscosity, k _0in is the initial relaxation time distribution function, ρ _gnf is the density, B(x) is the magnetic parameter, u ̲ and v ̲ are the coordinates of the flowing fluid, μ _gnf is the dynamic viscosity, and α _gnf is the thermal diffusivity. The nanofluid is represented by the notation gnf. The thermal attributes are represented as follows [40]: α gnf = k gnf ρ c p gnf , μ gnf = μ f c 1 − ϕ 2.5 , ρ c p gnf = ϕ ρ c p gnf + 1 − ϕ ρ c p f c , ρ gnf = ϕ ρ gnf + 1 − ϕ ρ f c , ν gnf = μ gnf ρ gnf , k gnf = k gnf k f c + 2 k f c 2 − k f c 2 − k f c k gnf 2 ϕ k gnf + 2 k f c + ϕ k f c − k gnf , σ gnf = 3 ϕ σ f c σ gnf σ f c − 1 − ϕ σ gnf σ f c − 1 + σ gnf σ f c + 2 + σ f c . (4) In Eq. 4, ρ _fc is the density, ρ c p f c = is the effective heat capacity, ϕ is the particle volume ratio, k _fc is the thermal conductivity of the base liquid, and k _gnf is the nanoparticle thermal conductivity. The accompanying similarity variables have been proposed to non-dimensionalize the basic equations [40]: u ̲ = c x f ′ η , v ̲ = − ν c 1 / 2 f η , η = y c ν 1 / 2 , θ η = T − T ∞ T w − T ∞ . (5)

Applying Eqs 5, 2 and Eq. 3 become as: f ‴ − τ 2 τ 1 f ′ 2 − τ 1 M g n f ′ + τ 2 τ 1 f f ″ − τ 2 τ 1 K v e − f f i v + 2 f ′ f ‴ − f ″ 2 = 0 , (6) f η = 0 , f ′ η = 1 + β s f ″ 0 , at η = 0 , f ′ η → 0 as η → ∞ . (7)

Here, the slip parameter is β _s = β(c/ν)^1/2, τ ₁ = (1 − ϕ)^2.5, τ 2 = ρ gnf ρ f c ϕ − ϕ + 1 , K v e = c k 0in ν f c is the viscoelastic parameter, and M g n = 2 β s σ f c B 0 2 ρ is the Hartmann number. The solution of Eq. 6 in a closed form is as follows [41]: f η = γ 1 + γ 2 e − ϒ η . (8) Using Eq. 7, we determine the answer to Eq. 6, as shown as follows: f η = 1 β s ϒ 2 + ϒ − e − ϒ η β s ϒ 2 + ϒ . (9)

Applying Eq. 6 and Eq. 9 together, we gain ϒ = λ 1 + 12 3 τ 1 λ 2 − λ 3 β s − 8 3 6 β s + 2 λ 4 3 β s λ 5 + 12 3 τ 1 λ 6 − λ 7 β s − 24 3 + τ 1 τ 2 K ve − 1 3 β s , (10) where λ 1 = 8 τ 1 3 τ 2 3 K ve 3 + 36 τ 1 2 τ 2 K ve β s 2 M gn − 24 τ 1 2 τ 2 2 K ve 2 + 72 τ 1 M gn β s 2 + 108 τ 1 τ 2 β s 2 + 24 τ 1 τ 2 K ve , (11) λ 2 = − K ve 2 M gn 2 β s 2 τ 1 3 τ 2 2 + 4 K ve 3 M gn τ 1 3 τ 2 3 + 4 K ve 3 τ 1 3 τ 2 4 − 4 M gn 3 β s 4 τ 1 2 + 20 K ve M gn 2 β s 2 τ 1 2 τ 2 + 18 K ve M gn β s 2 τ 1 2 τ 2 2 , (12) λ 3 = − 12 K ve 2 M gn τ 1 2 τ 2 2 − 12 K ve 2 τ 1 2 τ 2 3 + 8 M gn 2 β s 2 τ 1 + 36 M gn β s 2 τ 1 τ 2 + 27 β s 2 τ 1 τ 2 2 + 12 K ve M gn τ 1 τ 2 + 12 K ve τ 1 τ 2 2 − 4 M gn − 4 τ 2 , (13) λ 4 = τ 1 2 τ 2 2 K ve 2 + 3 τ 1 M gn β s 2 − 2 τ 1 τ 2 K ve + 1 , (14) λ 5 = 8 τ 1 3 τ 2 3 K ve 3 + 36 τ 1 2 τ 2 K ve β s 2 M gn − 24 τ 1 2 τ 2 2 K ve 2 + 72 τ 1 M gn β s 2 + 108 τ 1 τ 2 β s 2 + 24 τ 1 τ 2 K ve , (15) λ 6 = − K ve 2 M gn 2 β s 2 τ 1 3 τ 2 2 + 4 K ve 3 M gn τ 1 3 τ 2 3 + 4 K ve 3 τ 1 3 τ 2 4 − 4 M gn 3 β s 4 τ 1 2 + 20 K ve M gn 2 β s 2 τ 1 2 τ 2 + 18 K ve M gn β s 2 τ 1 2 τ 2 2 − 12 K ve 2 M gn τ 1 2 τ 2 2 , (16) λ 7 = 12 K ve 2 τ 1 2 τ 2 3 + 8 M gn 2 β s 2 τ 1 + 36 M gn β s 2 τ 1 τ 2 + 27 β s 2 τ 1 τ 2 2 + 12 K ve M gn τ 1 τ 2 + 12 K ve τ 1 τ 2 2 − 4 M gn − 4 τ 2 , (17) where γ ₁, γ ₂, and ϒ are the constants.

This section describes the heat exchange investigation under the influence of Joule heating and thermal radiation phenomenon. The following is the elementary equation [42]: u ̲ ∂ T ∂ x + v ̲ ∂ T ∂ y = α gnf ∂ 2 T ∂ y 2 − 1 ρ c p gnf ∂ q rad ∂ y + σ gnf B x 2 ρ c p gnf u ̲ 2 , (18) where q rad = − σ * 3 k * ∂ T 4 ∂ y . (19) Putting Eq. 19 into Eq. 18, we acquire the following: u ̲ ∂ T ∂ x + v ̲ ∂ T ∂ y = α gnf ∂ 2 T ∂ y 2 + 4 2 σ * T ∞ 3 k * 3 ρ c p gnf ∂ 2 T ∂ y 2 + σ gnf B x 2 ρ c p gnf u ̲ 2 , (20) and the boundary restrictions are as follows: T = T w = T ∞ + T 0cr x / c 2 at y = 0 , T → T ∞ as y → ∞ . (21)

Here, T is the temperature field, c is the characteristic length, T _0cr is the constant reference temperature, α _gnf is the thermal diffusivity, c p gnf is the specific heat, T _w is the wall temperature, k* is the mass absorption coefficient, T _∞ = is the free stream temperature, and σ* is the Stefan−Boltzmann constant. After applying Eq. 5 and Eq. 21 together, the corresponding non-dimensional energy equation is obtained as follows: Ω Pr θ η η − 2 f η θ + f θ η + E crt M g n τ 4 f η 2 = 0 , (22) where E crt = u ̲ 2 Δ T c p , Pr = ν f c α f c , χ = k * k f c 2 2 σ * T ∞ 3 , Ω = τ 3 τ 4 3 χ τ 3 + 4 3 χ τ 3 τ 3 = k gnf + 2 k f c − 2 ϕ k f c − k gnf k gnf + 2 k f c + 2 ϕ k f c − k gnf , τ 4 = 1 − ϕ + ϕ ρ c p gnf ρ c p f c . (23)

Here, χ is the radiation entity, E _crt is the Eckert number, and Pr is the Prandtl number. The modified boundary constraints are as follows: θ η = 1 at η = 0 , θ η → 0 as η → ∞ . (24)

Consequently, it is simple to obtain by entering Eq. 9 into Eq. 24, expressed as Ω Pr θ η η − 2 e − ϒ η β s ϒ + 1 θ + 1 β s ϒ 2 + ϒ − e − ϒ η β s ϒ 2 + ϒ θ η + E crt M g n τ 4 e − ϒ η β s ϒ + 1 2 = 0 . (25) To convert Eq. 26 into Kummer’s ordinary differential equation, a latest expression is established: ζ = − Pr e − ϒ η Ω ϒ 2 β s ϒ + 1 . (26) As a consequence, Eq. 26 is converted into Kummer’s ordinary differential equation, which is as follows: ζ ∂ 2 θ ∂ ζ 2 + κ − ζ ∂ θ ∂ ζ + 2 θ = − E crt M g n τ 4 e − ϒ η β s ϒ + 1 2 , (27) where κ = (1 − κ ₁) and κ 1 = Pr Ω ϒ 2 ( β s ϒ + 1 ) . The updated boundary requirements are as follows: θ ζ = 1 , θ 0 = 0 . (28)

With respect to Kummer’s functions [43], the closed-form solution of Eq. 28 accompanying Eq. 29 is: θ ( ζ ) = − ( ζ Ψ 1 Ψ 2 E crt , ϒ 2 ) Ω 2 Pr τ 4 Ψ 4 Ψ 5 + − ζ Ψ 6 ϒ 2 Ω β s ϒ + 1 Pr Ψ 7 Ψ 2 Ψ 9 2 Ψ 6 Ψ 7 Ψ 3 Ω 2 ϒ 4 τ 4 Ψ 4 Ψ 5 + ζ M g n E crt ϒ 2 Ω 2 Pr τ 4 , (29) where Ψ 1 = 1 − Pr β s ϒ + 1 Ω ϒ 2 2 + 2 ζ + 1 1 − Pr β s ϒ + 1 Ω ϒ 2 + ζ 2 + 2 ζ , Ψ 2 = M − 2 + Pr β s ϒ + 1 Ω ϒ 2 , 1 + Pr β s ϒ + 1 Ω ϒ 2 , − ζ , Ψ 3 = M − 2 + Pr β s ϒ + 1 Ω ϒ 2 , 1 + Pr β s ϒ + 1 Ω ϒ 2 , − Pr β s ϒ + 1 Ω ϒ 2 , Ψ 4 = 1 − Pr β s ϒ + 1 Ω ϒ 2 , Ψ 5 = 2 − Pr β s ϒ + 1 Ω ϒ 2 , Ψ 6 = − Pr β s ϒ + 1 Ω ϒ 2 , Ψ 7 = Pr β s ϒ + 1 Ω ϒ 2 , Ψ 8 = − 2 ζ E crt M g n ϒ 4 Ω 2 Pr β s ϒ + 1 Pr , Ψ 9 = 2 Ω 2 ϒ 4 τ 4 Ψ 4 2 + 2 Ω 2 ϒ 4 τ 4 Ψ 4 + Ψ 4 Ψ 8 − Ψ 8 + ζ E crt M g n Pr ϒ 2 Ω β s ϒ + 1 2 .

Here, M is the confluent hypergeometric function. The solution of the energy equation is as follows: θ η = ω 1 Ψ 2 E crt e − ϒ η 2 β s ϒ + 2 τ 4 Ψ 4 Ψ 5 + Ψ 6 e − ϒ η Ψ 7 ω 2 ω 9 2 Ψ 6 Ψ 7 Ψ 3 Ω 2 ϒ 4 τ 4 Ψ 4 Ψ 5 − E crt M g n e − ϒ η 2 β s ϒ + 2 τ 4 , (30) where ω 1 = 1 − Pr β s ϒ + 1 Ω ϒ 2 2 + 2 Pr e − ϒ η β s ϒ + 1 Ω ϒ 2 + 1 1 − Pr β s ϒ + 1 Ω ϒ 2 + Pr 2 e − ϒ η 2 β s ϒ + 1 2 Ω 2 ϒ 4 + 2 Pr e − ϒ η β s ϒ + 1 Ω ϒ 2 , ω 2 = M − 2 + Pr β s ϒ + 1 Ω ϒ 2 , 1 + Pr β s ϒ + 1 Ω ϒ 2 , − Pr e − ϒ η β s ϒ + 1 Ω ϒ 2 , ω 8 = 2 E crt M g n e − ϒ η ϒ 2 Ω Pr β s ϒ + 1 2 , ω 9 = 2 Ω 2 ϒ 4 τ 4 Ψ 4 2 + 2 Ω 2 ϒ 4 τ 4 Ψ 4 + Ψ 4 ω 8 − ω 8 − E crt M g n e ϒ η Pr 2 β s ϒ + 1 3 . θ η ′ 0 = − A 3 M g n E crt 2 β s ϒ + 2 τ 4 A 1 A 2 + A 4 M g n E crt ϒ 2 β s ϒ + 2 τ 4 A 1 A 2 + A 5 Pr A 8 A 6 2 β s ϒ + 1 Ω 3 ϒ 5 A 5 A 8 τ 4 A 1 A 2 − A 10 , (31) where A 1 = 1 − Pr β s ϒ + 1 Ω ϒ 2 , A 2 = 2 − Pr β s ϒ + 1 Ω ϒ 2 , A 3 = − 2 p r A 1 ϒ β s ϒ + 1 Ω − 2 Pr 2 β s ϒ + 1 2 Ω 2 ϒ 3 − 2 Pr ϒ β s ϒ + 1 Ω , A 4 = A 1 2 + 2 Pr β s α + 1 Ω ϒ 2 + 1 A 1 + Pr 2 β s ϒ + 1 2 Ω 2 ϒ 4 + 2 p r β s ϒ + 1 Ω ϒ 2 , A 5 = − Pr β s ϒ + 1 Ω ϒ 2 1 − A 1 , A 6 = − 2 E crt M g n A 1 ϒ 2 Ω Pr β s ϒ + 1 2 + 2 Ω 2 ϒ 4 τ 4 A 1 2 − 2 E crt M g n ϒ 2 Ω Pr β s α + 1 2 + 2 Ω 2 ϒ 4 A 4 A 1 − E crt M g n Pr 2 β s ϒ + 1 3 , A 7 = 2 E crt M g n ϒ 3 A 1 Ω Pr β s ϒ + 1 2 + 2 E crt M g n ϒ 3 Ω Pr β s ϒ + 1 2 + E crt M g n ϒ e Pr 2 β s ϒ + 1 3 , A 8 = M − A 2 , 1 + Pr β s ϒ + 1 Ω ϒ , A 1 − 1 , A 9 = M − A 1 , 2 + Pr β s ϒ + 1 Ω ϒ 2 , A 1 − 1 , A 10 = − A 9 Pr A 6 ϒ 5 β s α + 1 Ω 3 A 8 τ 4 A 1 2 + 2 Pr β s ϒ + 1 Ω ϒ 2 − 1 − A 7 1 2 Ω 2 ϒ 4 τ 4 A 1 A 2 .

The expression for the local skin friction is as follows: C f = τ w ρ u ̲ w 2 = R e x − 1 / 2 τ 1 f ″ 0 , τ 1 C f R e x − 1 / 2 = f ″ 0 , (32) where τ w = μ gnf ∂ u ̲ ∂ y y = 0 and R e x = x u ̲ w ν are the stress at the wall and the Reynolds number, respectively.

The local Nusselt number is defined as follows: N u = − k gnf x k f c T w − T ∞ ∂ T ∂ y y = 0 = − k gnf k f c R e x 1 / 2 θ η 0 . (33)

It is obtained as follows in the present study: k f c k gnf N u x R e x − 1 / 2 = − θ η 0 . (34)

We have analyzed the fluid flow and energy transport of the MHD viscoelastic nanofluid past a stretching surface with the effects of velocity slip and thermal radiation entity. The influence of several emerging parameters such as ϕ, M _gn , β _s , K _ve , E _crt , Pr, and χ on f′(η), θ(η), −f′′(0), and −θ′(0) is shown. In this connection, we have drawn the graphs and the corresponding numerical tables. The physical perspective of the considered model is depicted in Figure 1. Figures 2–5 are created to discuss the influence of these parameters on the velocity field. In Figure 2, the variation in ϕ above the velocity field is illustrated. It is noted that the magnitude of the velocity distribution is enhanced by an increment in the values of ϕ. Figure 3 depicts the variation in M _gn with β _s , K _ve , and ϕ. The velocity field diminishes as the M _gn increases. Physically, the Lorentz forces are dominant, which causes the velocity field to de-escalate while increasing the magnitude of M _gn . The influence of β _s on f′(η) is portrayed in Figure 4. It is perceived that increasing the amount of β _s causes the velocity distribution to drop rapidly. Actually, when the velocity slip occurs, the stretching surface velocity is faster than the fluid velocity, resulting in a decrease in the velocity of the nanofluid. The influence of K _ve on the velocity field is delineated in Figure 5. Physically, viscoelasticity creates a tensile force that opposes fluid movement. Due to viscous and elastic factors, escalating values of K _ve cause a reduction in the flow velocity.

Figures 6–11 investigate the effects of ϕ, M _gn , β _s , E _crt , Pr, and χ on θ(η). Figure 6 shows the temperature field against the nanoparticle ratio. It implies that the temperature field is enhanced with an enlargement in the magnitude of ϕ. Physically, the thermal conductivity of the fluid is enhanced by adding nanoparticles. Consequently, the temperature decreases. The behavior of the temperature field against M _gn is studied in Figure 7. It is revealed that a rise in M _gn causes the fall in temperature distribution. It is due to friction between the fluid and nanoparticles by Lorentz forces that leads to an increase in the temperature field. The development of Pr on θ(η) is shown in Figure 8. The temperature field is observed to rise, increasing the magnitude of Pr. The Prandtl number, in terms of physics, is the fraction of momentum diffusivity to thermal diffusivity. Thermal diffusivity is enhanced by adding nanoparticles, which causes a decrement in temperature distribution. Figure 9 illustrates the impact of E _ert over θ(η). It is discovered that increasing E _ert increases the magnitude of θ(η). The influence of the thermal radiation parameter on θ(η) is depicted in Figure 10. It demonstrates that the magnitude of θ(η) de-escalates, while the magnitude of χ increases. Figure 11 shows how decreasing the amount of χ increases the value of θ(η).

Figures 12–14 investigate the change in M _gn , β _s , and ϕ over −f′′(0). The local skin friction factor is decreased with the growing magnitude of M _gn as shown in Figure 12. Physically, magnetic fields and electric forces produce Lorentz forces that create resistance forces because of which the skin friction on the wall increases. Figure 13 and Figure 14 show a plot of −f′′(0) versus β _s and ϕ, respectively. It is perceived that the local skin friction field is decreased with an augmentation in β _s , and an opposite trend is shown for ϕ. Figures 15–18 scrutinize the variation in M _gn , E _crt , Pr, and χ over −θ′(0). The impact of M _gn on −θ′(0) against ϕ is presented in Figure 15, which shows a decrement in the amount of −θ′(0). Figure 16 depicts the effects of the Eckert number on the dimensionless temperature variation, −θ′(0). It can be seen that when the Eckert number increases, −θ′(0) also increases. The Eckert number E _crt is expressed as the ratio of advective transmission to energy-dissipated capability in physical terms. Furthermore, a rise in the rate of heat transport suggests that the energy loss is accompanied by a fall in the Eckert number. Figure 17 and Figure 18 express the influence of Pr and χ on the heat exchange ratio. It is discovered that θ′(0) is augmented while the strengths of Pr and χ are increased.

The thermophysical characteristics of graphene and CMC−water are displayed in Table 1. Tables 2, 3 show the numerical values of local skin friction coefficient and local Nusselt number respectively. Figures 19–24 indicate the streamline characteristics. It has been shown that the boundary-layer thickness dramatically shrinks when K _ve is increased.

In this study, the analysis of the MHD graphene−CMC−water nanofluid past a stretchable wall with Joule heating and velocity slip impact was performed. We arrived at the following conclusions:

• The velocity profile is the growing function of M _gn , β _s , and K _ve , with the reversed mode shown in case of ϕ.