A Novel Numerical Procedure for Simulating Steady MHD Convective Flows of Radiative Casson Fluids over a Horizontal Stretching Sheet with Irregular Geometry under the Combined Influence of Temperature-Dependent Viscosity and Thermal Conductivity

A novel mathematical computing analysis for steady magnetohydrodynamic convective flows of radiative Casson fluids moving over a nonlinearly elongating elastic sheet with a nonuniform thickness is established successfully in this numerical exploration. Also, the significance of an externally applied magnetic field with space-dependent strength on the development of MHD convective flows of Casson viscoplastic fluids is evaluated thoroughly by including the momentous influence of linear thermal radiation along with the temperature-dependent viscosity and thermal conductivity effects. By combining the assumption of the low-inducing magnetic field with the boundary layer approximations, the governing partial differential equations monitoring the current flow model are transmuted accordingly into a set of nonlinear coupled ordinary differential equations by invoking appropriate similarity transformations. Moreover, these derived differential equations are resolved numerically by utilizing a new innovative GDQLLM algorithm integrating the local linearization technique with the generalized differential quadrature method. On the other hand, the behaviours of velocity and temperature fields are deliberated properly through various graphical illustrations and different sets of flow parameters. However, the accurate datasets generated for the skin friction coefficient and local Nusselt number are presented separately in tabular displays, whose physical insights are discussed comprehensively via the slope linear regression method (SLRM). As main results, it is demonstrated that the higher values of the Casson viscoplastic parameter reduce significantly the fluid velocity within the boundary layer region, while a partial reverse tendency is observed near the stretching sheet as long as the wall thickness parameter is increased. Besides the previously mentioned hydrodynamical features, it is also depicted that the thermal field throughout the medium is enhanced considerably with the elevating values of these parameters.


Introduction
Linear and nonlinear rheological features of complex fluidic mediums have currently attracted widespread interests and considerable physical insights in the scientific debates owing to their valuable applications and widespread benefits over the whole usual fluids. Additionally, the non-Newtonian fluids have been extensively investigated during the past few decades by many scientists to explore their distinguishing rheological properties and establish new realistic constitutive laws permitting detailed explanations for some anomalous experimental observations. Indeed, the dynamic of most fluids together with their diverse characteristics cannot be properly described by the classical Navier-Stokes equations because the most common fluids used in industries are of a non-Newtonian kind, whose Cauchy stress and strain tensors can be expressed mathematically by quasilinear or purely nonlinear relationships depending greatly on their rheological behaviours, such as shear-thinning, shearthickening, viscoelasticity, viscoplasticity, material memory, Weissenberg effect, yield stress, die-swelling, and relaxation and retardation times, which are not defined in the Newtonian model. For this reason, the non-Newtonian fluids are widely used in various branches of science and bioengineering, for instance in the modelling of biological structures, drilling muds, cement slurries, petroleum, coated sheets, polymer manufacturing, foodstuffs, lubrication fields, and biomedical flows. Generally, the rheological attitudes of non-Newtonian fluids are dissimilar from those of viscous fluids. Besides, the transport equations governing the dynamical evolution of non-Newtonian fluids are fully nonlinear and highly complicated from the handling point of view compared with the Newtonian fluids. Furthermore, an analysis of the literature survey shows that there is no single constitutive equation that reveals all properties of such fluids. So, several constitutive non-Newtonian models have been suggested by the scientific communities. ese models are principally classified into the differential, integral, and rate-type categories. In the rate-type rheological fluid models (e.g., Maxwell, Oldroyd-B, Jeffery, and Burgers models), the shear stress is an implicit relationship of the velocity gradient and its higher temporal derivatives. However, the differential rheological model (e.g., Leonov model) has been established mathematically to explain the existence of a set of discrete relaxation times. is phenomenological model is based on the presence of an internal thermodynamic equilibrium state, where the equilibrium elastic deformations have been memorized and any aberration kind from this state can cause a new nonequilibrium state. In contrast to the aforesaid models, the extra stress tensor for the integral constitutive model (e.g., K-BKZ and Wagner models) is expressed mathematically as the temporal integral of the memory function and the elastic potential distributions as well as the Cauchy-Green tensor and its inverse (i.e., Finger tensor). Moreover, this model shows an excellent capability of giving comparable upshots against the experimental rheological results for long dies towards the extrudate swell occurrence witnessed in some viscoelastic materials, such as acrylonitrile-butadiene-styrene nanocomposites and other polymer melts. is rheological phenomenon is due mainly to the elastic recovery of deformation, which reflects the relaxation effect on the material after leaving from the die state.
Numerous constitutive models have been proposed by pioneering researchers to describe complex convective boundary layer flows. For instance, Bilal et al. [1] adopted the tangent hyperbolic rheological model to examine the characteristics of an unsteady MHD convective flow of a non-Newtonian fluid moving over a nonlinear stretching porous sheet under the combined influence of convective heating and time-dependent magnetic field in the presence of magnetized dusty particles, heat generation/absorption, and variable thermal conductivity. Similarly, Bibi et al. [2] accomplished comprehensive computational estimations with the aid of the shooting technique along with the built-in MATLAB bvp4c package to handle approximatively an unsteady MHD convective flow of a tangent hyperbolic fluid conveying tiny magnetized particles over a linear stretching permeable sheet characterized by a uniform suction/injection velocity and affected by the significant effects of variable thermal conductivity, internal heat source/sink, convective heating, and time-dependent magnetic field. Otherwise, Mohebbi et al. [3] utilized the power-law model to evaluate the shear-thinning and shear-thickening features of a non-Newtonian fluid by conducting a numerical investigation on a fully developed forced heat transfer convection for a laminar non-Newtonian flow confined between two parallel plates containing partially porous media in the presence of a uniform internal arrangement of circular obstacles. In the same way, Liu et al. [4] developed an approximate approach to scrutinize the non-Newtonian rheological behaviour of fresh concrete based on the Herschel-Bulkley model [5] and various available experimental data sources. Sequel to the aforesaid background, Zhu et al. [6] investigated numerically the mutual influence of fluid yield stress and finite-size particles on a turbulent flow of a Bingham fluid [7] flowing inside a vertical channel under the effect of no-slip velocity boundary condition and the periodicity physical constraints in the streamwise and spanwise Cartesian directions. Based on the Carreau-Yasuda rheological model and Buongiorno's nonhomogeneous nanofluid model, Waqas et al. [8] proposed a realistic physical model to analyze the impacts of second-order velocity slip, zero nanoparticles mass flux condition, thermal radiation, and chemical reaction on the features of MHD convective flows over a convectively heated stretching surface in a porous medium by considering the coexistence of gyrotactic motile microorganisms and tiny solid particles and taken into account the thermophoresis and Brownian motion aspects. Keeping in mind the nonsimilarity characteristic of some flow problems causing by the non-Newtonian behaviour of the flowing fluids, Amanulla et al. [9] carried out an accurate numerical implementation to find nonsimilar solutions for the momentum and thermal boundary layers developed in a non-Darcy porous medium during the steady two-dimensional MHD free convection flow of a Prandtl-Eyring fluid past over an isothermal sphere under the effective influences of buoyancy forces and various slip conditions. Earlier, Muhammad et al. [10] exploited Buongiorno's mathematical formulation along with the generalized Fourier-Fick's law and Rosseland's approximation to establish a theoretical nanofluid model allowing to evaluate more efficiently the Cattaneo-Christov heat transfer for a three-dimensional MHD-EMHD convective flow of an Eyring-Powell nanofluid flowing horizontally over an elongated Riga plate by considering the effects of nonlinear thermal radiation, activation energy, velocity slip condition, and convective heating. To explore the impact of shear-thinning and shear-thickening rheological properties on different kinds of transportation phenomena in non-Newtonian nanofluids under the combined effect of Brownian motion and thermophoresis diffusion of solid nanomaterials, Hayat et al. [11] extended the Sisko rheological fluid model [12] along with Buongiorno's mathematical formulation [13] and the vanishing mass flux condition for nanoparticles [14,15] to simulate a threedimensional MHD forced convective nanofluid flow over a nonisothermal surface, which is taken to be materially impermeable, heated convectively, and stretched bidirectionally via linear velocities. In the objective of inspecting the surface drag forces in nanofluidic mediums, as well as the feasibility of enhancing their heat and mass transportations, Hayat et al. [16] addressed an imperative unsteady threedimensional MHD squeezing viscoelastic flow of a secondgrade nanofluid between two isothermal parallel disks, whose chemical species are controlled actively via Buongiorno's two-phase model under the assumption of small magnetic Reynolds number theory along with the blowing and suction effects. To foresee the impacts of shear-thinning and shear-thickening rheological properties, Sajid et al. [17] performed an approximative analytical investigation on the two-dimensional coating process for incompressible thirdgrade fluid flows in a thin slit located between a moving substrate and flexible/stiff coater by supposing that the coating layer thickness is greatly less than the blade length. Likewise, Arifuzzaman et al. [18] utilized the explicit finite difference technique to analyze the heat and mass transfer features of unsteady MHD radiative-convective flows of fourth-grade fluids moving perpendicularly from an isothermal permeable surface in a porous medium under the effects of magnetic field, suction/injection velocity, nonlinear order chemical reaction, heat generation, viscous dissipation, and buoyancy forces.
Due to the progress in the rheology of homogeneous and nonhomogeneous viscoelastic mixtures employed in widespread industrial and engineering areas, various constitutive rheological relationships have been developed for biological and polymeric rate-type liquids to reveal the memory and elasticity features of these liquids by introducing the relaxation and retardation times via the Oldroyd-B rheological model [19]. In the range of low shear rate model, the dynamic viscosity of Oldroyd-B fluids found to be increasingly varied, while this trend is reversed for higher shear rates indicating that these materials possess both solid and liquid characteristics. In particular, these properties can be constrained rheologically in some cases to deduce the Maxwell model [20] by disappearing the retardation time characterized by the solid tendency of viscoelastic fluids, in which the fluidic medium cannot respond instantaneously to any change in shear stress. In this direction, it is worth noting here that the Oldroyd-B model is incapable to explain the shear-thinning and shear-thickening rheological phenomena observed in some physiological and biological fluids (e.g., blood, serum, urine, saliva, and tears). Further, this model exhibits a quasilinear relationship between the shear stress and the shear strain rate, while this rheological dependency becomes a linear relationship for the Maxwell model. Among the most common rate-type model that having a time derivative in its linear rheological constitutive law and offering the possibility to investigate easily the viscoelastic behaviour of non-Newtonian fluids, we find the constitutive rheological Jeffrey's fluid model [21]. Bearing in mind the importance of this rheological model and its simplicity in simulating many realistic viscoelastic flows encountered widely in the bioscience sector and technology, this linear rheological law has lately obtained exceptional interest by researchers. In this regard, Hayat et al. [22] exploited the constitutive rheological equation of the Jeffrey model to examine the probable effects of an applied magnetic field, viscous dissipation, and internal heat source on the stagnation point flow of a non-Newtonian viscoelastic fluid towards a nonlinear stretching surface, which is characterized by a variable surface thickness and heated isothermally via an appropriate melting process. In another relevant investigation, Hayat et al. [23] inspected a three-dimensional MHD convective flow of an Oldroyd-B nanofluid moving horizontally over a bidirectionally stretching surface under the effects of zero nanoparticles mass flux condition, convective heating, Brownian motion, and thermophoresis. Because of the substantial interest given to these non-Newtonian models, Koteswara Reddy et al. [24] proposed an adjustable rheological model to analyze the features of MHD mixed convection viscoelastic flows of Maxwell/Jeffery/Oldroyd-B nanofluids moving vertically over a nonisothermal conical geometry by considering the existence of an internal heat source/sink along with the impact of thermophoresis and Brownian motion of solid nanoparticles. In a particular case, Awias et al. [25] carried out a numerical implementation via the bvp4c technique to study three-dimensional mixed convection flows of Maxwellian nanofluids driven by a vertically rotating and stretching isothermal cylinder, under the influences of various interactions including thermal/mass buoyancy, magnetic, thermophoretic, and Brownian forces. On the other hand, the existence of multiple relaxation times in the response of some materials (e.g., Earth's mantle and asphalt mastics) in creep and stress relaxation examinations can be adequately explained by another specified rate-type fluid model of higher order. For this purpose, many researchers prefer modelling numerous complex flow problems with the help of the Burgers rheological model [26] due to its capability of describing the coexistence of two different relaxation times. In the framework of this advanced rate-type fluid category, Hayat et al. [27] assumed the generalized Fourier-Fick's law along with Buongiorno's nanofluid model to analyze the heat and mass transportation phenomena occurred during the steady two-dimensional convective flows of Burgers nanofluids over a linearly stretching isothermal sheet.
Later on, various experimental investigations indicated that the Casson rheological fluid model [28] is still the best rheological way to understand the non-Newtonian behaviour of many viscoplastic fluids, like chocolate, blood, tomato sauce, paint, shampoo syrup, some polymer solutions, and other colloidal suspensions. is kind of fluids exhibits a shear-thinning rheological behaviour (i.e., pseudoplastic feature), whose viscosity decreases with the increasing shearing rate as long as the stress upsurges above a critical value recognized as the yield stress. More precisely, the presence of mineral ingredients in a liquid increases significantly the dynamic viscosity of mediums and modifies its rheological characteristic to a non-Newtonian behaviour with a rheological tendency to an increase in the fluidity, when the resulting mixture is sheared. However, the fluidic mediums show the capability to gel under certain situations, when it is not sheared. In this context, Kelessidis and Maglione [29] realized an exhaustive statistical analysis on the Casson and Robertson-Stiff rheological models [28,30] via various existing experimental viscometric data generated especially for water-Greek bentonite suspensions to choose the excellent rheological category permitting a more realistic description of the studied aqueous medium. After successful attempts, they concluded that both suggested models can describe well the addressed sample data, in which the true shear rates for the Casson and Robertson-Stiff models are found to be higher compared with the Newtonian shear rates, whose differences become lesser at higher shear rates.
ese findings indicate an assessment in the Newtonian behaviour at greater shear rates. In the framework of the practical use of non-Newtonian Casson fluids, Qasim and Noreen [31] exploited the stability analysis theory to illustrate more evidently the possibility of generating multiple solutions for dissipative boundary layer flows of Casson fluids moving over a porous shrinking surface. More recently, Bhatti et al. [32] studied semianalytically MHD radiative Casson fluid flows from a nonisothermal linear stretching surface by including the Hall and ion-slip effects in the momentum equation via the generalized form of Ohm's law and imposing local convective heating as a feasible thermal condition at the horizontal boundary for the energy conservation equation.
Building upon the aforementioned literature review, the authors tackled their flow problems by utilizing various numerical and semianalytical procedures, such as Runge-Kutta method (RKM), spectral local linearization method (SLLM), spectral quasilinearization method (SQLM), spectral relaxation method (SRM), Keller box method (KBM), homotopy perturbation method (HPM), and homotopy analysis method (HAM), which are extensively employed in other related investigations. Besides, the analysis of these previous research works reveals that there is a greater lack of innovative research studies on the numerical analysis topic dealing realistically with the heat transportation phenomenon via non-Newtonian flows over complex geometries (e.g., irregular sheet, needle, cone, sphere, and cylinder). Motivating by these facts, a new numerical method is established in this scrutiny to explore the rheological fluid flow and heat transfer characteristics of radiative Casson viscoplastic fluids having temperaturedependent thermophysical properties (i.e., dynamic viscosity and thermal conductivity), surrounding by a nonuniform applied magnetic field and moving horizontally over an isothermal stretching sheet with variable thickness. Moreover, the graphical and tabular highlights produced for every physical quantity showing up in the present numerical assessment (i.e., velocity, temperature, skin friction coefficient, and heat transfer rate) are discussed physically through plausible explanations. Furthermore, the flow patterns are further displayed for velocity boundary layers and streamlines to inspect the consequence of applying a nonuniform magnetic field on the fluid flow properties.

Problem Statement
Let us consider a steady two-dimensional MHD convective flow of a Casson viscoplastic fluid exhibiting a temperature dependency for the dynamic viscosity and thermal conductivity and moving over a nonlinearly slender sheet of variable thickness ε(x), which is heated uniformly via an imposed temperature T w and stretched horizontally in the positive sense of the x-axis with a spatially variable velocity u x (x) under the impacts of thermal radiation and variable magnetic field Here, the dimensional constants (a, b, c) are utilized in the expressions of ε(x) and u x (x) to characterize the slender sheet geometry and its dynamical behaviour, while the exponent n designates the velocity power parameter. Also, the constant c is selected properly in this physical problem, in such a way that the sheet is sufficiently thin compared to its width ε(x). As schematically portrayed in Figure 1, the mathematical characteristics of the studied flow geometry are presented in a Cartesian coordinate system (x, y, z), whose unit vectors are (e x , e y , e z ). In this configuration, the origin O is located at the leading edge of the transverse section centre of the sheet, in which the x-axis is positioned horizontally along the symmetrical axis of the elastic sheet, while the y-axis is taken perpendicularly to the flow direction. As assumptions, the radiative fluid flow is supposed laminar, incompressible, and having a weaker electrical conductivity σ and can be developed with a low Reynolds magnetic number, in which the induced magnetic field can be neglected reasonably.
By adopting the linearized form of the Rosseland radiative heat flux model [33] along with the boundary layer approximations, the various transport equations (i.e., mass, momentum, and energy conservation equations) governing the present flow problem are stated as follows: As emphasized above, u and v represent the velocity components in the x and y directions, respectively, ρ means the fluid density, C P corresponds to the specific heat of the Casson fluid, β indicates the Casson fluid parameter, μ T refers to the plastic dynamic viscosity of the medium, T symbolizes the fluid temperature, k T specifies the thermal conductivity of the fluid, σ e signifies the Stefan-Boltzmann constant, β R designates the mean absorption coefficient, σ defines the electrical conductivity of the medium, and B(x) shows the strength of the variable magnetic field. 4 Mathematical Problems in Engineering By assuming a linear temperature dependency for the thermophysical quantities (μ T , k T ), we get [34] Moreover, the thermophysical properties (μ, k) shown above designate the effective values of the variable quantities (μ T , k T ) at the stream temperature T ∞ , whereas dimensionless quantities Λ μ and Λ k represent the variable viscosity and thermal conductivity parameters, correspondingly.
Accordingly, equations (1)-(3) can be nondimensionalized by invoking the following authentic similarity transformations: In this state, the resulting nonlinear ordinary differential equations (ODEs) with their proper boundary conditions are given by Slender sheet Mathematical Problems in Engineering 5 As underlined above, the prime superscripts designate the derivation with respect to ζ. Also, the physical parameters M, R D , Pr, and χ involved in equations (8)-(10) are defined formally as follows: For further simplifications, we set the following feasible changes: erefore, the nonlinear differential system described by equations (8)-(11) is altered to 1 Here, the prime notations are utilized to refer to the derivation for the variable η.

Engineering Quantities of Interest
In this physical problem, the significance of the principal engineering quantities of major importance can be scrutinized appropriately by estimating locally the magnitude of dragging forces and heat transfer rate at the stretching surface in terms of the skin friction coefficient C fx and Nusselt number Nu x . ese physical quantities are defined as follows: , in which By making use of equations (5), (7), and (13), we get the following expressions: As mentioned above, the physical quantities C fr and Nu r represent the reduced forms of the local skin friction coefficient C fx and Nusselt number Nu x , which are given by Here, Re x denotes the local Reynolds number, where Re x � ((x + a)u w (x)/υ).

Numerical Modeling Strategy
A new efficient numerical procedure is presented exclusively in this investigation to solve the nonlinear differential system described in the fluid motion and heat transfer features of a non-Newtonian viscoplastic fluid possessing temperaturedependent viscosity and thermal conductivity and flowing over an isothermal stretching sheet with irregular and axisymmetric geometry under the influence of a nonuniform applied magnetic field with the presence of thermal radiation. e proposed numerical algorithm is built on a purely theoretical framework. In this regard, the resulting nonlinear differential system of equations (14)- (17) can be linearized and decoupled accordingly to reduce them into a sequence of dependent subsystems of linear ODEs. After several successful rearrangements, the constructed differential subsystems are discretized spatially to be tackled numerically by utilizing the generalized differential quadrature local linearization method (GDQLLM).
Based on the pioneering mathematical findings reported by Bellman and Kalaba [35], the present flow model can be handled iteratively by applying the local linearization technique (LLT) to the following generated system: As highlighted in equation (22), the value η ∞ corresponds to the estimated dimensionless boundary layer thickness assuring the asymptotical achievement of the stream boundary conditions with higher precision. Consequently, the linearized form of the nonlinear differential system (S NL ) is written as follows: Moreover, the solutions of the linear differential system (S L ) can be adjusted iteratively through the functions (A i (η), B i (η), R i (η))/1 ≤ i ≤ 3 and 0 ≤ η ≤ η ∞ , which are expressed as follows:

Mathematical Problems in Engineering
Here, the iterative subscripts r and (r + 1) are utilized in equations (23)- (26) to indicate the previous and current estimations, respectively. Furthermore, to boost the numerical procedure approved during the resolution steps, it is preferable to initialize the iterative process by testing solutions (f 0 (η), g 0 (η), θ 0 (η))/0 ≤ η ≤ η ∞ verifying all the boundary conditions imposed on the differential system (S L ). ese initial guessed solutions take the following expressions: By applying the generalized differential quadrature method (GDQM) along with the modified Gauss-Lobatto collocation points [36][37][38][39][40][41], the resulting linearized system (S L ) can be discretized spatially to give three decoupled linear algebraic subsystems (S θ ), (S g ), and (S f ), whose unknowns are the approximate discrete values of the functions θ r+1 (η), g r+1 (η), and f r+1 (η) in the computational ese subsystems can be written in the following matrix forms: Additionally, the square matrices A, B, C { } and vector columns T r+1 , G r+1 , F r+1 , R 1 , R 2 , R 3 are expressed as follows: 8 Mathematical Problems in Engineering As shown above, D (1) and D (2) � [(d (2) ij ) 1 ≤ i,j ≤ N ] symbolize the first-and second-order differentiation matrices, respectively, I denotes the unit matrix, the set η j /1 ≤ j ≤ N highlights the modified Gauss-Lobatto collocation points, and N represents the number of collocation nodal points η j .
According to Shu [42], the components of the modified differentiation matrices D (1) and D (2) can be computed from the following general system: Furthermore, it bears noting that, after discretizing the differential system (S L ), the derivative terms can be handled numerically in each collocation point η i as follows: where p is the order of derivation concerning the variable η. During the progression in the iterative procedure, the reduced quantities C fr and Nu r are computed as follows:

Accuracy and Convergence of the Proposed Numerical Algorithm
e convergence and constancy of the sequential iterative schemes established for the present boundary layer flow problem can be assessed fittingly from the nonlinear differential system (S NL ) by considering the following residual functions: After discretization, these residual functions take the following vector forms:

Mathematical Problems in Engineering
For the whole computational domain, we define the following infinite norms: By applying the above definitions on the constructed vectors (Res g , Res θ ) at any advanced iteration step (r + 1), the efficiency of the computed results can be authenticated persuasively for every set of the control parameters (i.e., η ∞ , N, n, χ, β, Pr, M, R D , Λ μ , and Λ k ) by utilizing the following convergence criteria: Based on the above considerations, the iterative numerical schemes converge strongly to higher accurate outcomes, when the systematic errors E d and E Res are less or equal than the specified convergence tolerance levels ε d and ε Res , respectively. On the other hand, the error values ε d and ε Res are generated adequately with magnitude orders of about 10 − 9 and 10 − 7 , correspondingly, to get more precise numerical results. Moreover, the computational cost corresponding to the GDQLLM code execution can be assessed automatically by measuring the CPU time via the MAT-LAB's Tic-Toc function. In parallel, various proper numerical implementations are carried out via this numerical code to estimate the approximate possible values of C fr , Nu r , ‖Res g ‖ ∞ , ‖Res θ ‖ ∞ , and CPU time in Table 1 against various iteration numbers, when η ∞ � 100, N � 250, n � 2, χ � 0.5, β � 0.5, Pr � 7, M � 2, R D � 1, Λ μ � 0.8, and Λ k � 1. From this tabular illustration, it is clear that the proposed numerical method GDQLLM exhibits a faster convergence response, higher flexibility, slighter computational CPU time, and higher robustness, in which the necessary convergence criteria can be attained excellently after just 10 repetitions with a CPU time less than 0.66 s.

Validation of Numerical Results
More importantly, the engineering applicability of GDQLLM as a new innovative numerical method can be realistically corroborated by making use of the aforementioned accuracy criteria and linking the present flow model with other existing related literature studies. For this purpose, the developed GDQLLM code is scrutinized extensively by comparing appropriately its particular results against those obtained previously by Fang et al. [43], Wang [44], and Mabood et al. [45] for the quantities −f ″ (0) and −θ ′ (0), when β − 1 � M � R D � Λ μ � Λ k � 0. Sequel to these multiple comparative tests enlisted in Tables 2 and 3, the accurateness of the GDQLLM results is strengthened more by the SLLM upshots generated by the spectral local linearization method [46]. A quantitative analysis of the tabular data proves that the implemented GDQLLM code exhibits very satisfactory numerical outcomes, which lends a high level of credibility to the present investigation.

Discussion of Results
After choosing the more appropriate values for the extra parameters (η ∞ , N), the GDQLLM results can be obtained successfully, when the number of iterations is achieved a sufficient value N it providing smaller residual errors with a higher order of precision. For this main purpose, it is preferable to select η ∞ � 100, N � 250, and N it � 100 as the main key parameters controlling the correctness of the present outcomes in all subsequent numerical implementations. By exploiting the accuracy numerical datasets (f N it (η i ), θ N it (η i ))/1 ≤ i ≤ N and 0 ≤ η i ≤ η ∞ outputted during the completion of the code execution for specified flow parameter values (i. e., n, χ, β, Pr, M, R D , Λ μ , and Λ k ), the fluid flow properties and heat transfer characteristics are revealed throughout the medium in terms of velocity and temperature fields, as well as locally via the skin friction coefficient and Nusselt number at the slender sheet. ese physical features are presented and discussed via numerous graphical and tabular displays by taking n � 2, χ � 0.5, β � 0.5, Pr � 7, M � 2, R D � 1, Λ μ � 0.8, and Λ k � 1 as default flow parameter values, unless otherwise indicated. Furthermore, to elucidate exhaustive physical insights from the graphical representations, the velocity and temperature profiles are clipped from 0 to 8 to show the impacts of every embedded parameter more visibly and avoid the overlap between the curves in each graphical illustration.
As exhibited below, the graphical and tabular representations displayed in Figures 2-10 and Tables 4 and 5 reflect completely the fluid motion and heat transfer characteristics (i.e., velocity profiles, temperature profiles, skin friction coefficients, heat transfer rates, velocity boundary

Existing results
Present results Fang et al. [43] GDQLLM SLLM  Table 3: Comparison of the present numerical findings for −θ′(0) with the latest existing literature works, when n � 1 and

Pr
Existing results Present results Wang [44] Mabood et al. [   layers, and streamlines) of the studied non-Newtonian fluid against sundry influencing parameters, namely, velocity power parameter n, wall thickness parameter χ, Casson fluid parameter β, Prandtl number Pr, magnetic parameter M, radiation parameter R D , variable viscosity parameter Λ μ , and variable thermal conductivity parameter Λ k . Generally, it is observed that the dimensionless profiles of the horizontal velocity distributions f ′ (η i ) � g(η i )/1 ≤ i ≤ N and 0 ≤ η i ≤ η ∞ are monotonically declining functions and tends asymptotically to the vanishing boundary condition f ′ (η ⟶ η ∞ ) � 0 as long as the similarity variable η exceeds a certain distance η � δ ( ≤ 8), which indicates the approximate value of the velocity boundary layer thickness. e probable evolution of the dimensionless velocity and temperature profiles (f ′ (η), θ(η)) in response to the increasing change in the values of the power index parameter n are delineated quite clearly in Figure 2. Graphically, it is shown that the dynamical parameter n exhibits a heightening effect towards the thermohydrodynamic profiles (f ′ (η), θ(η)) with a significant widening in the thermohydrodynamic boundary layer thicknesses. ese circumstances can be explained more by the higher increment in the magnitude of the nonlinear stretching velocity of the heated slender sheet as long as the velocity exponent n is upsurged, which leads to thicker thermohydrodynamic boundary layers.
Additionally, Figure 3 is plotted to elucidate the sensitive effect of the wall thickness parameter χ on the dimensionless velocity and temperature profiles (f ′ (η), θ(η)). Seemingly, the fluid flow shows a dual behaviour concerning the geometrical parameter χ depending on the nearness or the farness from the stretching sheet, whereas a thermal enhancement is disclosed throughout the fluidic medium. After analysing carefully the graphical results highlighted in Figure 3, it is found that the higher values of the wall thickness parameter χ cause a tenuous increase in the fluid velocity near the stretching sheet; however, a slower fluid motion is observed far from the horizontal boundary to get almost a thinner velocity boundary layer. Since the velocity exponent n is greater than the unit value in this problem (i.e., n � 2), the function f(η) takes negative values near the stretching sheet and decreases with the increasing values of the wall thickness parameter χ, due to the imposed boundary condition f(0) � ((1 − n)χ/(1 + n)). So, the rising trend revealed for the velocity profile f ′ (η) near the stretching sheet against the increasing values of the geometrical parameter χ can be explicated quite clearly by invoking the vanishing boundary condition for the dimensional transversal velocity v(x, y), considering the positive sign of the dimensionless horizontal velocity f ′ (η) and making use of the transformations defined in equations (7) and (13) for the expression of v(x, y). Furthermore, the thermal improvement noticed in the fluidic medium can be clarified by the blowing impact arising from the augmentation in the wall thickness parameter χ.
Rheologically, the Casson parameter β reflects the viscoplastic non-Newtonian behaviour of the present boundary layer flow. Moreover, the Casson viscoplastic fluid model can be reduced to the Newtonian viscous fluid case just by escalating the rheological parameter β to an extremely higher value (i.e., β ⟶ ∞). Also, an increase in the Casson parameter β rises the plastic dynamic viscosity of the Casson fluid and decreases the resulting yield stress. ereby, a viscous resistance is induced to the flow in the fluidic medium. For this reason, the fluid velocity lessens with the escalating values of the Casson viscoplastic parameter β, while an opposite effect is detected for the fluid temperature as portrayed in Figure 4. Consequently, the Casson parameter β has the proclivity to suppress the velocity boundary layer region and heighten the thermal boundary layer thickness. Via a physical standpoint, the temperature dependency involved in the dynamic viscosity and thermal conductivity expressions causes a fully mathematical coupling between the conservation equations (i.e., momentum and energy equations) and engenders an implicit effect of the Prandtl number Pr on the fluid flow and heat transfer features as depicted in Figure 5. Due to the augmenting values of the Prandtl number Pr, the temperature of the medium diminishes and the fluid slows down slightly near the isothermal horizontal boundary, while an accelerating impact is achieved after a certain distance from the sheet, which in turn provides an enhancement in the velocity boundary layer thickness. Besides, the Prandtl number Pr illustrates a direct link between the fluid motion inside the velocity boundary layer region and its corresponding temperature, since it represents the ratio between the kinematic viscosity υ (� μ/ρ) of the viscoplastic medium to its thermal diffusivity α (� k/(ρC P )).
erefore, an enhancement in the Prandtl number Pr weakens the thermal diffusive influence in front of the momentum diffusive impact. As a result of the viscous dominance, the thermal boundary layer becomes more thinner and the fluid motion shows a weak deceleration near the stretching sheet, whilst an opposite behaviour is perceived whenever we get away from the isothermal sheet because of the thermophysical interference noticed between the impacts of the variable properties arising from the temperature-dependent parameters (Λ μ , Λ k ) introduced in the governing equations.
On the one hand, the illustrative variations in the dimensionless velocity and temperature profiles (f ′ (η), θ(η)) against the varying values of the magnetic parameter M are established graphically in Figure 6. From this graphical visualization, it is evident that an increase in this flow parameter leads to a lessening in the fluid velocity. is slowing behaviour demonstrated for the fluid velocity towards the intensifying values of the magnetic parameter M can be attributed to the existence of Lorentz forces, which are happened in the electrically conducting fluidic medium due to the presence of an external magnetic field applied transversely to the fluid flow. Besides, these magnetic forces act meaningfully as a resistive factor by weakening the velocity boundary layer thickness as evidenced in Figure 7(a) and enlarging its corresponding thermal width. ese facts are occurred because of both the deceleration in the fluid motion and the elevation in its temperature. On the other hand, Figure 7(b) proves that the flow patterns are affected noticeably by the presence of an externally applied magnetic, in which the values of the dimensionless stream function decline notably, due to the slowing down of the fluid velocity within the momentum boundary layer. Moreover, the streamline shapes become wider and tend to be parallel as long as we go away from the stretching sheet. Energetically, the thermal radiation parameter R D (� (16σ e T 3 ∞ )/(3β R k)) represents the effective contribution rate between the radiative and conductive heat transports. So, an intensification in the radiation parameter R D causes an increase in the radiative manner of heat transfer, which in turn leads to a strengthening effect on the fluid motion near the hot slender sheet with an enhancement in both the temperature field and its boundary layer thickness as emphasized graphically in Figure 8. Also, it is witnessed that the fluid motion shows a slowing behaviour towards an upsurge in the values of the radiation parameter R D when going away from the stretching sheet, which implies a reduction in the velocity boundary layer thickness.
Additionally, the inevitable features of the variable viscosity parameter Λ μ ( > 0) on the dimensionless velocity and Figure 9. From this graphical representation, it is revealed that the rising values of the variable viscosity parameter Λ μ have a decelerating impact on the horizontal fluid motion, while a thermal enhancement is depicted in the medium due to this increasing variation in the control parameter Λ μ . As mentioned before, the linear temperature dependency adopted for the physical propriety μ T implies that the variable viscosity parameter Λ μ has a declining effect on the dynamic viscosity μ T , which results in turn to a shearthinning rheological characteristic for the studied non-Newtonian fluid. Indeed, the increasing values of the variable viscosity parameter Λ μ cause a stronger temperature difference (T w − T ∞ ) between the stretching sheet and the free-stream fluid. So, larger values of the variable viscosity parameter Λ μ increase the fluidity of the medium and strengthen the fluid motion in the vertically upward direction as compared with its horizontal flow. For this reason, a diminishing trend is seen for the hydrodynamic boundary layer thickness, while an opposite upshot is found for the thermal boundary layer thickness.
Unlike the aforediscussed results for the variable viscosity parameter Λ μ , the consequences of supposing a linear temperature dependency for the thermal conductivity k T on the dimensionless velocity and temperature profiles (f ′ (η), θ(η)) are underscored in Figure 10. Generally, an elevation in the variable thermal conductivity parameter Λ k leads to an important enrichment in the fluid thermal conductivity, which in turn improves the fluid temperature. However, a significant depletion in the dynamic viscosity μ T is inspected within the fluidic medium due to the augmentation in its liquidity. Keeping in mind that the flow and heat transport mechanisms are notably influenced by the thermal conductivity along the transverse direction, the fluid motion is reinforced near the heated starching sheet with the increasing values of the variable thermal conductivity parameter Λ k , while a reverse trend is seen away from the stretching sheet because of the vertical motion tendency. Consequently, the variable thermal conductivity parameter Λ k has a diminishing influence on the hydrodynamic boundary layer thickness, whilst an enlargement is noticed for the thermal boundary layer zone. By applying the statistical procedure SLRM [47][48][49] on the accurate GDQLLM and SLLM datasets generated for the present viscoplastic flow model, the fluid flow and heat transfer features against the rising values of each involved flow parameter (i.e., n, χ, β, Pr, M, R D , Λ μ , and Λ k ) are studied systematically based on the reduced physical quantities C fr and Nu r arranged in Tables 4 and 5. From these tabular demonstrations, it is found that the surface frictional effect arising from the amplifying values of the engineering quantity C fr can be minimized either by intensifying the values of wall thickness parameter χ, Casson fluid parameter β, radiation parameter R D , variable viscosity parameter Λ μ , and variable thermal conductivity parameter Λ k or by lowering the values of velocity power parameter n, Prandtl number Pr, and magnetic parameter M. However, an enhancement in the surface heat transfer rate can be attained either by mounting the values of the parameters Pr, R D , and Λ k or by lessening the values of the parameters n, χ, β, M, and Λ μ .

Concluding Remarks
Numerous interesting physical and mathematical insights have been realistically drawn for the present MHD radiativeconvective flow problem in the case of an enhanced Casson viscoplastic fluid having temperature-dependent thermophysical properties (i.e., dynamic viscosity and thermal conductivity) and moving horizontally over an irregular isothermal stretching sheet by employing an innovative numerical approach, called the generalized differential quadrature local linearization method (GDQLLM). In this framework, noteworthy GDQLLM outcomes have been derived from this numerical investigation, which can be itemized as follows: Data Availability e data of the present problem cannot be shared due to confidentiality requirements.

Conflicts of Interest
e authors declare that they have no conflicts of interest.