Magnetohydrodynamic double-diffusive peristaltic flow of radiating fourth-grade nanofluid through a porous medium with viscous dissipation and heat generation/absorption

This article focuses on determining how to double diffusion affects the non-Newtonian fourth-grade nanofluids peristaltic motion within a symmetrical vertical elastic channel supported by a suitable porous medium as well as, concentrating on the impact of a few significant actual peculiarities on the development of the peristaltic liquid, such as rotation, initial pressure, non-linear thermal radiation, heat generation/absorption in the presence of viscous dissipation and joule heating with noting that the fluid inside the channel is subject to an externally induced magnetic field, giving it electromagnetic properties. Moreover, the constraints of the long-wavelength approximation and neglecting the wave number along with the low Reynolds number have been used to transform the nonlinear partial differential equations in two dimensions into a system of nonlinear ordinary differential equations in one dimension, which serve as the basic governing equations for fluid motion. The suitable numerical method for solving the new system of ordinary differential equations is the Runge–Kutta fourth-order numerical method with the shooting technique using the code MATLAB program. Using this code, a 2D and 3D graphical analysis was done to show how each physical parameter affected the distributions of axial velocity, temperature, nanoparticle volume fraction, solutal concentration, pressure gradients, induced magnetic field, magnetic forces, and finally the trapping phenomenon. Under the influence of rotation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega$$\end{document}Ω, heat Grashof number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Gr}_{t}$$\end{document}Grt, solutal Grashof number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Gr}_{t}$$\end{document}Grt, and initial stress \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P}^{*}$$\end{document}P∗, the axial velocity distribution \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$u$$\end{document}u changes from increasing to decreasing, according to some of the study’s findings. On the other hand, increasing values of nonlinear thermal radiation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document}R and temperature ratio \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\theta }_{w}$$\end{document}θw have a negative impact on the temperature distribution \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta$$\end{document}θ but a positive impact on the distributions of nanoparticle volume fraction \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi$$\end{document}Φ and solutal concentration \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma$$\end{document}γ. Darcy number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Da$$\end{document}Da and mean fluid rate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F$$\end{document}F also had a positive effect on the distribution of pressure gradients, making it an increasing function.


− → D
The displacement field D B The Brownian diffusion coefficient D T The thermophoresis diffusion D TF The Dufour diffusivity D F T The Soret diffusivity Da The Darcy number − → E The electric field intensity E The electric field strength Ec The Eckert number F The dimensionless time mean flow rate Gr t The heat Grashof number Gr p The nanoparticles Grashof number Gr c The solutal Grashof number g The gravitational acceleration − → H The magnetic field H o Constant magnetic field H + The total magnetic field Hr The Hartman number h x The induced magnetic field I The identity tensor − → J The current density K p The permeability.k * The mean absorption coefficient k The thermal conductivity.

J he
The Joule heating M The magnetic field parameter Nb The Brownian motion parameter Nt The thermophoresis parameter P The pressure N F T The Soret parameter N TF The Dufour parameter Pr The Prandtl number P * The Initial stress parameter Q o The constant of heat conduction − → q The velocity vector − → q r The radiation heat flux R The thermal radiation parameter

Re
The Reynolds number R m The magnetic Reynolds number S i The extra stress S xx , S xy , S yy The extra stress tensors S The Strommer number T The fluid temperature D s The solutal diffusivity T o The Initial temperature tr The trace t The time U, V , u, v The velocity components in the fixed and wave frames X, Y , x, y The fixed frame and the wave frame The wavelength µ e The magnetic permeability ς The magnetic diffusivity ε * The porosity µ app The apparent viscosity σ The Electrical conductivity ρ f The fluid density ρ p The density of nanoparticles τ The superscript is the matrix transpose − → The rotation ∇ The differential operator θ w The temperature ratio parameter δ The wave number α i

Material constants β i
Material constants γ i

Material constants d/dt
The material time derivative ∂P/∂x The pressure gradient α The amplitude ratio The heat generation/absorption parameter (ρc P ) f The heat capacity of the fluid µ The dynamic viscosity (ρc P ) p The effective heat capacity of the nanoparticles ω 3 The vortex γ The non-dimensional solutal concentration function ψ The stream function T The Cauchy stress tensor θ The non-dimensional temperature function ξ The time-mean flow ϑ The dimensionless time-mean flow β F The volumetric solutal expansion β T The thermal expansion coefficient ζ The wave frame φ The non-dimensional magnetic force function The non-dimensional concentration of nanoparticles function σ * The Stefan-Boltzmann constant y The step size F The solutal (species) concentration The initial solutal (species) concentration The component of peristaltic pumping or the peristaltic movement of liquids inside living frameworks is perhaps the main strategy that certainly stands out to analysts in the field of biomechanics for its extraordinary significance and the significant job it plays inside the body of the creature.Latham 1 was the first to study the transmission of fluids by the peristaltic pumping mechanism or the peristaltic wave method.Peristaltic pumping of fluids through flexible tubes is defined as the movement of fluids within these channels utilizing pressure gradients that push them making them move along these tubes, such as the movement of blood within the veins and arteries, the movement of food within the flexible esophageal channel, as well as the movement of food within the small and large intestines.On the other hand, urine moves from the kidney to the urinary bladder through the peristaltic pumping of urine through the flexible ureteral canal, as well as the movement of semen from the testicle to the outside through peristaltic pumping through the canal of the carrier vessel and many other forms of peristaltic movement within living systems.Akram and Razia 2 concentrated on the crossover impacts of warm and focus convection on the peristaltic stream of fourth-grade nanofluids in an inclined tapered channel.According to the findings of this study, as the Brownian motion coefficient values rise, so does the fluid's temperature distribution.Kothandapani et al. 3 investigated magnetic field impact on the peristaltic flow of a fourth-grade fluid in a tapered asymmetric channel.The impact of the magnetic field on the axial velocity distribution is nonuniform and changes among increment and reduction, as indicated by the aftereffects of this study.Abdulhadi and Ahmed 4 examined the peristaltic flow of Walters-B fluid through a porous medium in a tapered asymmetric channel in the present magnetic field.It was found through this study that the temperature distribution turned into a rising function affected by the thermophoresis coefficient and the Dufour coefficient while the solute concentration distribution became a decreasing function under the influence of the Brownian motion coefficient and the Soret coefficient.Asha and Deepa 5 analyzed the Joule heating and magnetic field effect on the peristaltic flow of a third-grade fluid in an asymmetric channel.Ranjit et al. 6 illustrated Joule heating and zeta potential effect on peristaltic blood flow through porous micro-vessels altered by electrohydrodynamics.It was found from this study that the effect of the Brinkman number on the temperature distribution was positive.
It is realized that non-Newtonian liquids are one of the kinds of liquids that have an extraordinary nature, which is that they consolidate consistency and flexibility in their properties together.The viscosity of non-Newtonian fluids, on the other hand, varies depending on other factors like pressure and temperature.As a result, the viscosity of each model of non-Newtonian fluid is unique because the relationship between stress and strain for these fluids is nonlinear.In actuality, non-Newtonian fluids respond to external factors in a significant way.At the point when this liquid is exposed to a bunch of outside drives, it abandons a fluid substance to an adaptable, rubbery strong and when these powers are debilitated, it gets back to what it was previously.Non-Newtonian liquids assume a significant part in most significant applications that rely upon their properties including non-Newtonian liquid erosion decrease, oil-pipeline grating decrease, surfactant applications to huge scope warming and cooling frameworks, and stream tracers.It is additionally utilized in numerous enterprises, for example, paints, toothpaste, glue, lubricants, car oils, and foodstuffs such as ketchup and others.As a result, numerous researchers attempt to investigate the flow issues of various non-Newtonian fluids; El-Dabe et al. 7 used a numerical approach to solve the MHD peristaltic flow of a non-Newtonian power-law nanofluid through a non-Darcy-porous-medium-inside-non-uniform inclined channel.This work found that the Darcy number and Reynolds number had a negative effect on the axial velocity distribution.Sayed et al. 8 studied the peristaltic transport of nano non-Newtonian fluid with heat transfer through an inclined asymmetric duct with slip and convective boundary conditions.It was observed that the axial velocity distribution changed irregularly between decreasing and increasing under the influence of heat Grashof number one of the outcomes of that study.Ali et al. 9 investigated the impact of the Hartmann boundary layer in peristaltic flow for viscoelastic fluid.Nadeem et al. 10 provided an elliptic cross-sectional illustration of the peristaltic flow of a heated Jeffrey fluid within a duct that study revealed that an increase in the axial velocity distribution was caused by the mean flow rate parameter.Akram and Saleem 11 discussed the analysis of heating effects and different waveforms on the peristaltic flow of Carreau fluid in a rectangular duct.It was found that the effect of the mean flow rate coefficient on the axial velocity distribution is positive, one of the outcomes of that study.Abbas et al. 12 debated the peristaltic flow of Casson fluid with heat and mass transfer in the presence of Soret and Dufour effects with Hall currents through a non-Darcy porous medium inside a vertical channel.According to the findings of this study, both the temperature distribution and the nanoparticles volume fraction distribution were affected negatively by the Grashof thermal number.A peristaltic blood flow of non-Newtonian Carreau fluid with heat transfer phenomenon through a curved channel was discussed by Tanveer et al. 13 .The results of this study indicate that the axial velocity distribution is positively influenced by the Darcy number and the heat Grashof number.He and Mostapha 14 agitated the significance of the Hall current and Joule heating impacts on Darcy-Forchheimer's peristaltic flow of a Rabinowitsch fluid through a tapered tube.
A substance with pores (voids) is known as a porous medium or porous material.The material's skeletal component is frequently referred to as the matrix or frame.The pores are normally loaded up with a liquid (fluid or gas).The skeletal material is typically strong, however, structures like froths are many times additionally helpfully investigated utilizing the idea of permeable media.A permeable medium is most frequently described by its porosity.Different properties of the medium (e.g., penetrability, rigidity, electrical conductivity) can in some cases get from the particular properties of its constituents (strong lattice and liquid) and the media porosity and pores structure, however, such a deduction is normally complicated.The idea of porous media, which is utilized in a lot of applied science and engineering fields: filtration, mechanics (acoustics, geomechanics, soil mechanics, rock mechanics), designing (oil designing, bio-remediation, development designing), geosciences (hydrogeology, petrol topography, geophysics), science and biophysics, material science, and so on.Liquid course through permeable media is a subject of most normal interest and has arisen in a different field of study.On the other hand, any natural and industrial areas are very interested in the flow through porous media.Models remember the progression of hydrocarbons for oil wells, groundwater courses through beds of rocks, as well as the vehicle of minerals and pollutants through the ground.One may likewise find this stream peculiarity straightforwardly significant in vaporous or fluid synergist and idle pressed bed reactors moves through separates channels, geothermal intensity the executives, softening, or cementing of double composites.Abou-zeid and Ouaf 15 demonstrated the hall currents effect on squeezing the flow of non-Newtonian nanofluid through a porous medium between two parallel plates.Mccash et al. 16 elucidated entropy analysis of the peristaltic flow of a hybrid nanofluid inside an elliptic duct with sinusoidally advancing boundaries.Ahmed et al. 17 explained the effects of nonlinear thermal radiation, convective boundary conditions, and heat generation/absorption on magnetohydrodynamic Maxwell nanofluid flow over a stretching surface through a porous medium.It was found from this study that the effects of the magnetic field coefficient and Darcy number on the velocity distribution is negative.Mohamed et al. 18 analyzed the nonlinear thermal radiation and heat generation/absorption effect on MHD Jeffrey nanofluid flow over a stretching sheet through a porous medium.This study included that the nanoparticles volume fraction distribution was an increasing function under the influence of the temperature ratio coefficient, while it became a decreasing function under the influence of the Brownian motion coefficient.The influences of thermal radiation and slip conditions on MHD Casson nanofluid flow over a stretching surface through a porous medium was studied by Mohamed et al. 19 .The temperature distribution in this study has become in a state of continuous increase under the influence of both the absorption and heat generation coefficient, as well as the thermal transfer coefficient.Bouslimi et al. 20 dissected MHD Williamson nanofluid flow over a stretching sheet through a porous medium under the effects of Joule heating, nonlinear thermal radiation, heat generation/ absorption, and chemical reaction.The velocity distribution decreases with an increase in the magnetic field coefficient, while the temperature distribution increases with an increase in this parameter, and this is one of the outcomes of that study.Rashed and Ahmed 21 studied the peristaltic flow of dusty nanofluids in curved channels.This study proved that the heat Grashof number affects positively on the temperature distribution and negatively on the nanoparticles volume fraction distribution, while the effect of the nanoparticles Grashof number is the opposite.Ahmed and Rashed 22 investigated magnetohydrodynamic dusty hybrid nanofluid peristaltic flow in curved channels.Mohamed et al. 23 analyzed the thermal radiation and MHD effects on the free convective flow of a polar fluid through a porous medium in the presence of internal heat generation and chemical reaction.Abd-Alla et al. 24 studied the effects of rotation and initial stress on peristaltic transport of fourth-grade fluid with heat transfer and induced magnetic field.The study shows that the axial velocity distribution decreases with an increase in the rotation coefficient, and increases with the enhancement in the heat generation and absorption coefficient.Abd-Alla et al. 25 investigated the influence of heat and mass transfer, initial stress, and radially varying magnetic fields on the peristaltic flow in an annulus with a gravity field.Abd-Alla et al. 26 illustrated the effect of rotation on the peristaltic flow of a micropolar fluid through a porous medium with an external magnetic field.Abd-Alla and Abo-Dahab 27 elaborated on a magnetic field and rotation effects on the peristaltic transport of a Jeffrey fluid in an asymmetric channel.Abd-Alla et al. 28 illustrated the peristaltic flow of a Jeffrey fluid under the effect of a radially varying magnetic field in a tube with an endoscope.
Nanofluids are a new class of fluids that was discovered by Choi and Eastman 29 .They are a basic fluid that is permeated with tiny metal particles called nanoparticles which are measured in nanometers.On the other hand, nanofluids play an important role in many modern technological applications as they are the most intelligent type of liquid.In the past few years, many researchers in the field of nanofluids have worked hard to understand the behavior of these fluids.From previous investigations, nanofluids have been found to possess enhanced thermophysical properties such as thermal conductivity, thermal diffusivity, viscosity, and convective heat transfer coefficients compared to those of base fluids like oil or water.Nanofluids can be used in a plethora of engineering applications ranging from use in the automotive industry to the medical arena to use in power plant cooling systems as well as computers.Ibrahim and Anbessa 30 examined the three-dimensional MHD mixed convection flow of Casson nanofluid with hall and ion slip effects.Hayat et al. 31 examined the impacts of slip in radiative MHD peristaltic flow of fourth-grade nanomaterial with chemical reaction.It was noted from the consequences of this study that the impact of the magnetic field coefficient on the axial velocity is negative, the temperature  32 analyzed the hybrid effects of thermal and concentration convection on the peristaltic flow of fourthgrade nanofluids in an inclined tapered channel.A portion of the consequences of this study incorporated that the impact of the magnetic field coefficient on the axial velocity distribution was negative, and the temperature dispersion was a rising function affected by the thermophoresis coefficient.Kothandapani et al. 33 discussed the effect of a magnetic field on the peristaltic flow of a fourth-grade fluid in a tapered asymmetric channel.It has been observed that the axial velocity distribution fluctuates between an increase and a decrease under the influence of the magnetic field coefficient.Haroun 34 explained the non-linear peristaltic flow of a fourth-grade fluid in an inclined asymmetric channel.Akram et al. 35 illustrated slip impact on double diffusion convection of magneto fourth-grade nanofluids with peristaltic propulsion through an inclined asymmetric channel.Hayat and Noreen 36 investigated the peristaltic transport of fourth-grade fluid with heat transfer and an induced magnetic field.A portion of the consequences of this study showed that the heat Grashof number affects irregularly the axial velocity distribution and the induced magnetic field distribution.Mohamed et al. 37 studied the MHD threedimensional flow of a couple of stress nanofluids over a stretching sheet through a porous medium in the presence of heat generation/absorption and nonlinear thermal radiation.Abo-Dahab et al. 38 discussed the Double-diffusive peristaltic MHD Sisko nanofluid flow through a porous medium in the presence of non-linear thermal radiation, heat generation/ absorption, and Joule heating.Abd-Alla et al. 39 examined the rotation and initial stress effect on MHD peristaltic flow of reacting radiating fourth-grade nanofluid with viscous dissipation and Joule heating.Some results of this study showed that the distribution of axial velocity changes its effect between increase and decrease under the influence of rotation, initial stress, and magnetic field coefficient.Mohamed et al. 40 analyzed the effects of nonlinear thermal radiation and heat generation/absorption on magnetohydrodynamic Carreau nanofluid flow on a nonlinear stretching surface through a porous medium.Raja et al. 41 elaborated the integrated intelligent computing application for the effectiveness of Au nanoparticles coated over MWCNTs with velocity slip in curved channel peristaltic flow.Abdelhafez et al. 42 illustrated the influence of an inclined magnetic field and heat and mass transfer on the peristaltic flow of blood in an asymmetric channel.Rafiq et al. 43 studied the theoretical exploration of thermal transportation with Lorentz force for a fourth-grade fluid model obeying peristaltic mechanism.Abbas et al. 44 investigated the rheology of peristaltic flow in a couple of stress fluids in an inclined tube.Abbas et al. 45 analyzed the peristaltic transport of a Casson fluid in a non-uniform inclined tube with Roseland approximation and wall properties.Abbas and Rafiq 46 examined the numerical simulation of thermal transportation with viscous dissipation for a peristaltic mechanism of micropolar-Casson fluid.Ghazanfari et al. 47 sicussed the numerical study on the thermal performance of the shell and tube heat exchanger using twisted tubes and Al 2 O 3 nanoparticles.Manh et al. 48elaborated the investigation of nanomaterial flow through non-parallel plates.Sheikholeslami et al. 49 illustrated hybrid nanoparticles dispersion into water inside a porous wavy tank involving magnetic force.Awais et al. 50examined the MHD effects on ciliary-induced peristaltic flow coatings with rheological hybrid nanofluid, Ali et al. 51 explained the entropy generation on MHD peristaltic flow of Cu-water nanofluid with slip conditions.Rooman et al. 52 investigated the entropy generation analysis of magnetized radiative Ellis (Cu-TiO 2 /Engine Oil) nanofluid flow using Cattaneo-Christov heat flux model with viscous dissipation and Joule heating effects.Jameel et al. 53 elaborated the statistical and entropy optimization modeling for radiative hybrid nanofluid flow with Hall effect over exponential stretching/shrinking plate.Asghar et al. 54 analyzed the dual solutions of convective rotating flow of three-dimensional hybrid nanofluid across the linear stretching/shrinking sheet.The double-diffusion of peristaltic flow of a non-Newtonian fourth-grade nanofluid through a vertically symmetric flexible channel through a porous medium was the subject of our investigation in this paper.The primary objective of this numerical study is to investigate the effects of rotation, initial stress, and induced magnetic field in the presence of non-linear thermal radiation, heat generation/absorption, and Joule heating with viscous dissipation.For the system of equations governing the peristaltic flow process, it is a set of partial differential equations that consist of the continuity equation, the momentum equation, the energy equation in addition to the solute concentration equation, the nanoparticles volume fraction equation, and the magnetic forces equation, and induced magnetic field equation.The limitations of the long-wavelength approximation were used to facilitate the solution of the aforementioned system of differential equations and thus the system of these differential equations was converted to the system of ordinary differential equations which was solved numerically using the fourth-order Runge-Kutta method with the shooting technique by MATLAB code.The effect of all important physical parameters resulting from the current study on important distributions such as velocity, temperature, solute concentration, nanoparticle volume fraction, stream function, and induced magnetic field was studied employing the graphs using the MATLAB program and the results are discussed with the physical meaning of the phenomenon.

Discerption of the problem and mathematical formulation
At the outset, it should be noted that the current numerical study is represented in the peristaltic motion of fluid inside the vertical symmetrical flexible channel called a fourth-grade nanofluid, A fourth-grade fluid is defined as a category of non-Newtonian fluids with very high viscosity and describe the shear thinning and shear thickening phenomena which cannot be expressed by the classical Navier-Stokes equations, and this study is a simulation of the movement of food and liquids from the mouth to the stomach through the flexible esophageal channel, food, and liquids work on the contraction and relaxation of the esophageal channel walls, which is what happens to the fourth-grade fluid movement inside the flexible channel used in the study additionally, the study of the model's non-Newtonian fluid peristaltic motion's natural properties in general.It is worth noting that the peristaltic movement of the fluid inside the channel occurs in the presence of a suitable porous medium with permeable properties that acts as a resistance force that affects the movement of the fluid inside the channel.The porous medium contains a group of holes or voids permeated by the fluid during its movement, In the simulation or applied study, it corresponds to the membrane lining the walls of the flexible esophageal channel, as it is a semi-permeable membrane, or food can be considered as a porous medium permeated with fluids also the size of the pores of the porous medium is appropriate to the size of the nanoparticles and that the mechanism of action of the porous medium is subject to Darcy's law, see 55 .Consider the magnetohydrodynamic (MHD) peristaltic steady flow of an incompressible fourth-grade nanofluid through a porous medium in a vertical twodimensional symmetric channel of uniform thickness 2a subjected to Joule heating, initial stress, induced mag- netic field, and rotation in presence of nonlinear thermal radiation, heat generation/absorption, and viscous dissipation, the fluid flow in the channel walls is produced, when the sinusoidal waves of the small amplitude b propagate the speed of the channel walls.c is the constant speed of the channel.In this current study, the Cartesian coordinate system will be used and will be selected X in the direction of wave propagation and Y transverse to it, see Fig. 1.A constant magnetic field of strength H o acting in the transverse direction results in an induced magnetic field and the total magnetic is The geometry of the symmetric channel can be written as follows 56 ; The velocity field − → q in the two-dimensional flow can be written as follows; The governing equations of the non-Newtonian nanofluid motion can be written as follows 24,39,56 ; (1) where the Eqs.( 5)- (10) in general are the continuity equation, the momentum equation, the energy equation, the solute concentration equation, the nanoparticles volume fraction equation, and the magnetic force equation respectively.
Cauchy stress tensor T for the non-Newtonian fourth-grade fluid can be written as follows 24 : where α i (i = 1, 2),β i (i = 1, 2, 3) and γ i (i = 1, 2, . . ., 8) , The Rivlin-Ericksen tensors A i (i = 1to4) are defined as follows 24,39,56 : In magnetic fluid mechanics Maxwell's equations in the absence of displacement current are presented as follows 24 ; Modified Darcy's law obeys 55 Now, x, y represent the wave frame moving with velocity c away from the fixed frame X, Y by the transformations 24 ; Now, the basic governing equations of MHD non-Newtonian fourth-grade nanofluid can be written mathematically as follows 24,39,48 ; (6) The stream function ψ , velocity fields, and the dimensionless axial induced magnetic field are defined as follows; Substituting Eq. ( 23) in Eqs. ( 18)-( 22) will produce (18)    It is important to note that the values of each S xx , S xy , and S yy can be changed in Eqs. ( 11)-( 13) respectively.The final form of the Eqs.( 26)-( 30) is as follows when using the long-wavelength approximation and disregarding the wave number and low Reynolds number in turn: The boundary conditions for the dimensionless in the wave frame are given as follows; Now, the extra stress tensor for fourth-grade fluid can be written as follows 31 where α 1 , β 1 , β * , γ 1 and γ * are material constants and µ is the coefficient of viscosity.Substituting Eq. ( 38) into Eq.( 31) results in that Thus, the final form of Eq. ( 39) becomes the form By Eqs. ( 31) and ( 40) we get

Rate of volume flow and boundary conditions
The dimensional rate of fluid flow in the laboratory frame is And in wave frame, the above equation is as follows By using Eqs.( 20), (33), and (34) it can be seen that The time-mean flow over a period T at a fixed position X is given as follows After using Eq. ( 36) the above expression becomes as follows The dimensionless time-mean flow ϑ and F can be defined in the laboratory and wave frame respectively, as follows; (38)

Numerical solution
The Runge-Kutta numerical method is one of the most important numerical analysis methods that are used in solving a system of ordinary differential equations, there are different formulas for the solution of the numerical method of Runge-Kutta method; for example, there is the Runge-Kutta method of the fourth order and also the Runge-Kutta method of the fifth order, and the most used method is the Runge-Kutta method of the fourth-order with shooting technique because it gives accurate results and is easy to use, its derivation depends on the Euler method, and the error rate in the solutions is very small and this is greatly required for the credibility of the results and also for their ease.Accordingly, the system of Eqs. ( 32)-( 35) and ( 41) is a new system of nonlinear ordinary differential equations that describe the study of the fourth-grade nanofluid inside the flexible channel in addition to the boundary conditions ( 36) and (37).This system of equations has been solved numerically using the fourth-order Runge-Kutta method with the shooting technique using MTLAB program, this method depends on the basic step size used in the calculation is y = 0.001 so that this value was reached after conducting many numerical experiments to reach the independence of the network and uses a variable step size to calculate the sharp changes in the variables in the areas dominated by viscosity inside the channel and where the limits the beginning and the end depend on the size of the basic step, so they are taken −1.5 ≤ y ≤ 1.5 it has been obtained by substituting in h 2 and h 1 in the boundary conditions ( 36) and ( 37) by ε = 0.5 and x = 1 .Our equations are dealt with in MATLAB because each nth-order equation is mutated into n of the first-order equations, then using the bvp4c function (in general, bvp4c is a finite difference code that implements the three-stage lobatto IIIa formula.This is a collocation formula and the collocation polynomial provides a C1-continuous solution that is fourth-order accurate uniformly.Mesh selection and error control are based on the residual of the continuous solution) to solve these first-order equations.The governing system of ordinary differential equations can be written in simplest form as follows; The system of differential Eqs. ( 48)-( 52) has been inserted into the MATLAB program as follows The boundary conditions can be written as follows

Results and discussion
Table 1 represents a comparison between previous works that considered some parameters used in the present estimation to show the differences between the related topics.The current study resulted in the emergence of a set of physical parameters that you, dear reader, can identify when looking at the symbol table.
The effect of each of the physical parameters resulting from the study on the distributions of axial velocity u , temperature θ , solute concentration γ , nanoparticles volume fraction , pressure gradients ∂P ∂x , magnetic force φ , and finally induced magnetic field h x has been studied taking into account, the study of the behavior of the streamlines responsible for the occurrence of the trapping phenomenon under the influence of some of these parameters.It should be noted that when studying the effect of any parameter on any of the aforementioned distributions, the values of the rest of the other parameters remain constant at the following values In some detail, the results of this study will be discussed as follows.

Axial velocity distribution.
The effects of the rotation 0.5 ≤ ≤ 3 and the initial stress 1 ≤ P * ≤ 6 on the axial velocity distribution u are depicted in Figs. 2 and 3 (2D and 3D).It was discovered that their impact was positive in the distance −1.5 ≤ y ≤ 0.25 and soon turned negative in the distance 0.25 ≤ y ≤ 1.5 , as they caused the axial velocity distribution u to exhibit irregular behavior.The fact that the direction of rotation is perpendicular to the axis of the flexible channel and the fluid is subjected to rotation around the axis of the tube when applied is the reason why rotation causes an irregular behavior in the axial fluid velocity on a physical level.A flow is rotational if its vortices vector is not zero in all of its regions while the fluid elements are rotating about their axis as it flows along streamlines.Also, the fluid that is moving inside the channel responds strongly to initial stress P * , causing turbulence between the fluid layers and making the fluid's velocity and movement irregular.In contrast, Figs. 4 and 5 (2D and 3D) can be used to explain the effects of heat Grashof number 1 ≤ Gr t ≤ 6 and solutal Grashof number 1 ≤ Gr c ≤ 6 on the axial velocity distribution u .You will notice that both numbers cause an enhancement in the axial velocity distribution in the interval −1.5 ≤ y ≤ 0.25 , while a reduction of this distribution occurs in the interval 0.25 ≤ y ≤ 1.5 .This indicates that both numbers also produce irregular behavior in the axial velocity distribution u .Physically, it is important to note that the ratio between the fluid's thermal buoyancy force and its viscous hydrodynamic force is not fixed.This means that an increase in the fluid's thermal buoyancy force directly causes a decrease in the fluid's viscosity, making it  14),       Temperature, solutal concentration, and nanoparticles volume fraction distributions.Take a look at Figs. 9 and 10 (2D and 3D) to learn how the non-linear thermal radiation parameter 0.1 ≤ R ≤ 1.7 and the temperature ratio parameter 0.1 ≤ θ w ≤ 0.6 affect the distributions of temperature, solutal concentration, and nanoparticles volume fraction.It was found that increasing the values of the previous two parameters caused a noticeable negative change in the distribution of temperature and an increase in the distributions of both solutal concentration and nanoparticles volume fraction.Physically, the inverse relationship between the thermal conductivity coefficient and these two parameters, on the one hand, and the occurrence of thermal diffusion far from the system, on the other, explains why there was a significant decrease in the temperature distribution under the influence of the various increasing values taken by the non-linear thermal radiation coefficient and also the temperature ratio coefficient it was assumed that there would be an enhancement in the temperature  distribution under the influence of the various increasing values taken by these two parameters.This also caused an increase in the interfacial distances between the basic fluid particles, with a concentration of nanoparticles higher inside the channel and a concentration of nanoparticles lower near the channel walls.This also explains the support in the distribution of the solvent concentration and the volume fraction of the nanoparticles.Figures 11 and 12 (2D and 3D) depict the effects of heat Grashof number 1 ≤ Gr t ≤ 6 and solutal Grashof number 1 ≤ Gr c ≤ 6 on the temperature θ , solutal concentration γ , and nanoparticles volume fraction distributions, respectively.It was found that these two numbers had a positive effect on the temperature distribution, while they had a negative effect on the distributions of both solutal concentration and nanoparticles volume fraction.Physically, The relative relationship between the thermal buoyancy force of the liquid and the viscous hydrodynamic force of the liquid is determined by the heat Grashof number 1 ≤ Gr t ≤ 6 , and accordingly, the enhance- ment in the temperature distribution resulting from the influence of the different increasing values of the thermal Grashof number was mainly caused by the dominance and control of the thermal buoyancy force over the hydrodynamic viscous force for the fluid, this means a higher temperature in exchange for lower viscosity, on the one hand.On the other hand, this also causes an increase in the movement of the basic fluid molecules, and thus a decrease in the interstitial distances between those molecules, and this means a decrease in the solutal concentration distribution, in addition to a decrease in the concentration of nanoparticles inside the fluid, and this is  www.nature.com/scientificreports/also a decrease in the nanoparticles volume fraction.It should be noted that the solutal Grashof number 1 ≤ Gr c ≤ 6 follows the same approach as the heat Grashof number 1 ≤ Gr t ≤ 6 .To concentrate on the impacts of the nanoparticle's Grashof number 1 ≤ Gr p ≤ 6 and the Brinkman number 1 ≤ Bn ≤ 6 on the way of behav- ing of the distributions of temperature θ , solutal focus γ , and nanoparticle volume division Φ, dear reader, you ought to take a gander at Figs. 13 and 14 (2D and 3D), you can see that these two numbers concur in their positive effect on the temperature distribution θ while they concur in their adverse consequence on the solutal con- centration distribution γ , and nanoparticles volume fraction distribution .Physically, the Grashof number of the nanoparticle,1 ≤ Gr p ≤ 6 has a significant impact on the behavior of the nanoparticle and its relationship to the temperature distribution.As a result, the positive effect of the temperature distribution under the influence of this number is to increase the thermal buoyancy of the nanoparticles within the fluid at the expense of the fluid viscosity.Additionally, this has a negative impact on the distribution of the solvent concentration because of the shorter distances between the underlying fluid particles.The activity of the thermal conductivity coefficient, which in turn works on the transfer of temperature from the channel walls to inside the fluid because the Brinkman number 1 ≤ Bn ≤ 6 expresses the heat transfer process from the wall inside the system and this also causes an increase in the movement of the basic fluid molecules and nanoparticles and a decrease in the interfa-  cial distances between them, is the primary reason for the enhancement of the temperature distribution in contrast to the decrease in the other two distributions under the influence of the various increasing values of the Brinkman number.Increasing the various values of the thermophoresis parameter 0.1 ≤ Nt ≤ 2.1 supports and enhances the temperature distribution θ in exchange for reducing the distributions of the solutal concentration γ and the nanoparticle volume fraction this behavior is depicted in Fig. 15 (2D and 3D).Physically, thermo- phoresis is a thermal phenomenon that occurs in a specific medium with a temperature difference.It causes matter particles to move from one region of the medium with a higher temperature to one with a lower temperature until the temperature is the same in all of the medium's regions.This phenomenon is controlled by a thermophoresis coefficient 0.1 ≤ Nt ≤ 2.1 .It is strengthening and increasing values have led to an increase in the fluid temperature as a result of the movement of fluid particles at a high temperature, the explanation for the reduction in the distributions of the solutal and the nanoparticles volume fraction is the scattering and precariousness that happens to the essential liquid atoms and nanoparticles inside the liquid when the upsides of the thermophoresis coefficient are improved.The enhancement in the values of the magnetic field 0.5 ≤ M ≤ 3 is straightforwardly trailed by an upgrade in the temperature distribution θ as well as, a diminishing in the distributions of solutal concentration γ, and nanoparticles volume fraction Φ, this is obvious to you, dear peruser in      Pressure gradient distribution.The pressure gradient plays an important role in the movement of the fluid within the elastic channel, as it is like a machine that works to pump the fluid and its flow at different speeds.Figure 21 explains the effect of rotation 0.5 ≤ ≤ 3 on the distribution of pressure gradients.It has been  observed that the pressure gradient increases significantly under the influence of the enhancement in rotation, and this works to support the movement of the fluid inside the channel due to rotation, which in turn resists the large viscosity that characterizes the fluid, it also works to increase the waves occurring in the channel walls, which move with the movement of the fluid inside the channel.Figure 22 shows the negative effect of the Darcy number 0.5 ≤ Da ≤ 3 on the distribution of pressure gradients so that the enhancement in the values of the Darcy number reduces the rate of pressure gradient, and the reason behind this is the strengthening of the resistive force resulting from the porous medium.On the other hand, it works to reduce the wave's incident in the walls of the flexible channel.The support that occurred in the distribution of pressure gradients under the effects of the magnetic field parameter 0.5 ≤ M ≤ 3 and the heat Grashof number 1 ≤ Gr ≤ 6 can be seen in Figs.Magnetic force, and induced magnetic field distributions.The magnetic force is one of the four known forces in nature and is considered the product of the electromagnetic force, or in other words, it is the result of the electromagnetic force, and the main cause of the magnetic force is the movement of charges.For example, when two bodies contain the same charge and have the same direction of motion, a magnetic attraction force is generated between them, while if the direction of motion is opposite, then a magnetic repulsion force is generated.Accordingly, Figs. 26 and 27 studies the impact of electric field strength coefficient 0.2 ≤ E ≤ 2.2 , and the electric Reynolds number 0.1 ≤ R m ≤ 1.1 on the distribution of the magnetic force φ .The distribution of the magnetic force that was directed on the movement of the fluid inside the channel under the influence of the two previous parameters becomes at its lowest value if you look at the two figures in order.This indicates that the fluid movement within the channel is subjected to slight values of magnetic force.In contrast, the induced  magnetic field is defined in electrical physics as the magnetic field produced by the passage of an electric current through a wire.It is directly responsible for the magnetic force produced when a magnet is brought close to a piece of metal, and its direction is perpendicular to the electric current's.In the present study, when exposing the fluid that flows inside the channel to an induced magnetic field, it may exhibit a behavior that is recognized by influencing that field with three physical factors, which are, in order, electric field strength coefficient 0    and 35a-c show the behavior of streamlines profiles, their density, and their relationship to the phenomenon of trapping under the influence of different values taken by each of the thermophoresis Nt ∈ {0.2, 0.5, 0.7} and the heat Grashof number Gr t ∈ {0.1, 1, 0.6} , respectively.The enhancement in the values of the two parameters led to a systematic increase in the density of the streamlines with the emergence of more than one bolus on both sides of the channel noting the occurrence of compressions on the boluses at the bottom of the channel as a result of the convergence of the streamlines compared to the boluses above the channel so that these boluses are trapped by ripples that move with it in the same speed.In general, it can be concluded that the phenomenon of trapping resulting from the state of the flow lines that describe the behavior of the fluid inside the channel from the bottom is different from the top so that the phenomenon of trapping increases significantly from the bottom to the top, causing irregular movement in general within the channel.The reader ought to focus on; the current work is due to Abdalla et al. 24 in the absence of a porous medium, nonlinear thermal radiation, heat generation/absorption, viscous dissipation, and Joule heating, and there is no solutal concentration equation and nanoparticle volume fraction equation with the difference in the geometric shape of the channel and the boundary conditions.In contrast, the graphs in Fig. 36a-c show a comparison between the current study and Abd-Alla et al. 's 24 investigation of the impact of the electric Reynolds number coefficient R m and the electric field strength coefficient E on the magnetic force distribution φ , as well as the impact of the strength magnetic field H o on the induced magnetic field distribution h x .This demonstrates the congruence and similarity between the two investigations, which lends considerable support to the numerical method that was used in obtaining the results of the current study.

Conclusion
The subject of this study is the double-diffusion peristaltic flow of a non-Newtonian fourth-grade nanofluid through a symmetric vertical elastic channel in the presence of a suitable porous medium.The flow process is also subject to the strength of the induced magnetic field and is affected by rotation, initial stress, nonlinear thermal radiation, heat generation/absorption, viscous dissipation, and Joule heating.The flow process is controlled by     www.nature.com/scientificreports/ a system of nonlinear partial differential equations with boundary conditions that have been transformed into a system of nonlinear ordinary differential equations with a new form of boundary conditions using the long wavelength approximation and ignoring the wave number and the lower Reynolds number.The new system of ordinary differential equations was mathematically solved numerically using the fourth-order Runge-Kutta strategy with the shooting technique in MATLAB code.Finally, graphs were made using MATLAB software to examine the influence of the relative polymorphism of the important physical parameters due to the constant focus on the main distributions namely axial velocity, temperature, solutal concentration, nanoparticle volume fraction, magnetic forces, induced magnetic field, and streamlines profiles.The following is a list of the most important findings of the study: • Under the influence of the rotation , the initial stress coefficient P * , the heat Grashof number Gr t , and the solutal Grashof number Gr c , the axial velocity distribution u changes from increasing to decreasing.• Under the influence of the Darcy number Da , the nanoparticles Grashof number Gr p , and the mean flow coefficient F , the axial velocity distribution u shifts from a decreasing function to an increasing function.• The heat Grashof number Gr t , the solutal Grashof number Gr c , the nanoparticles Grashof number Gr p , the Brinkman number Bn , the thermophoresis Nt , the magnetic field M , the heat generation/absorption β , the Brownian motion Nb , the Soret parameter N F T , and the Dufour parameter N TF all have a positive impact on the temperature distribution θ.• The temperature distribution is a decreasing function under effects of the nonlinear thermal radiation R and the temperature ratio θ w work on the decrease of this distribution.• The distribution of solutal concentration γ and nanoparticles volume fraction were negatively impacted by the heat Grashof number Gr t , the solutal Grashof number Gr c , the nanoparticles Grashof number Gr p , the Brinkman number Bn , the thermophoresis Nt , the magnetic field M , the heat generation/absorption β , the Brownian motion Nb , the Soret parameter N F T , and the Dufour parameter N TF .• The solutal concentration γ and nanoparticle volume fraction distributions were positively impacted by the nonlinear thermal radiation R and the temperature ratio coefficient θ w .• The pressure gradients distribution ∂P ∂x turned into a rising function affected by the rotation , the heat Grashof number Gr t , the magnetic field M , and the mean flow F , while it turned into a negative function affected by the Darcy number Da.
• The magnetic forces distribution φ decreases under the influence of the electric field strength E and the electric Reynolds number R m .• The behavior of the induced magnetic field distribution h x changes between decreasing and increasing under the influence of the electric field strength E , the electric Reynolds number R m , and the magnetic field strength H o .

Figure 1 .
Figure 1.The geometry of the problem.

Y a ( 12 )
= (�) y=−1.5 Y b (12) = (�) y=1.5 Y a (14) = (φ) y=−1.5 Y b (14) = (φ) y=1.5 easier for the fluid to move, and the other way around.According to what was said, the heat Grashof number Gr t describes the relative relationship between the fluid's thermal buoyancy force and its viscous hydrodynamic force.Therefore, an increase in the fluid velocity distribution is accompanied by a decrease in the fluid viscosity rate at a distance of −1.5 ≤ y ≤ 0.25 from the channel, while the opposite is true at a distance of 0.25 ≤ y ≤ 1.5 from the channel.The behavior of the solutal Grashof number Gr c is identical to that of the heat Grashof num- ber Gr t .Plots of the changes in the axial velocity distribution u caused by Darcy number 0.5 ≤ Da ≤ 3 and nanoparticles Grashof number 1 ≤ Gr p ≤ 6 are shown in Figs.6 and 7(2D and 3D).It was found that both numbers caused a decrease in the axial velocity distribution u in the range −1.5 ≤ y ≤ 0.25 of the channel, but that behavior quickly changed to an increase in the same distribution in the range 0.25 ≤ y ≤ 1.5 of the chan- nel.Physically, it is well known that a porous medium is a resistive medium that is permeated by many voids or

Figure 2 .
Figure 2. Effect of Ω on the axial velocity distribution u (2D and 3D).

Figure 3 .
Figure 3.Effect of P * on the axial velocity distribution u (2D and 3D).

Figure 4 .
Figure 4. Effect of Gr t on the axial velocity distribution u (2D and 3D).

Figure 5 .
Figure 5.Effect of Gr c on the axial velocity distribution u (2D and 3D).

Figure 6 .
Figure 6.Effect of Da on the axial velocity distribution u (2D and 3D).

Figure 7 .
Figure 7. Effect of Gr p on the axial velocity distribution u (2D and 3D).

Figure 8 .
Figure 8.Effect of F on the axial velocity distribution u (2D and 3D).

Figure 9 .
Figure 9.Effect of R on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 10 .
Figure 10.Effect of θ w on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 11 .
Figure 11.Effect of Gr t on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 12 .
Figure 12.Effect of Gr c on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 13 .
Figure 13.Effect of Gr p on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Fig. 16 (
Fig.16 (2D and 3D).Physically, the presence of a drag force that prevents the development of momentum and results in a decrease in the flow inside the channel causes friction between the layers of the viscous fluid and the walls of the channel, resulting in an increase in the fluid's temperature, this is the primary reason for the improvement in the temperature distribution under the influence of the magnetic field parameter 0.5 ≤ M ≤ 3 .Like- wise, the reaction of nanoparticles and natural liquid particles inside the liquid to the strength of the attractive field made their fixations decline.It turned out, dear reader, that the effect is positive on the temperature distribution θ and negative on the other two distributions γ , � when the heat generation/absorption parameter −0.4 ≤ β ≤ 0.1 was plotted in Fig.17(2D and 3D) specifically to explain its effect on the three main distribu- tions temperature θ , solutal concentration γ , and nanoparticles volume fraction .Physically, the expansion in the field of temperature dissemination is communicated from heat absorption (β < 0) to heat generation (β > 0) by the normal active energy of the liquid atoms coming about because of their retention of thermal power, which prompts an expansion in the development and speed of the fluid particles inside in light of the fact that the temperature decides the movement energy related with the development of liquid particles and nanoparticles which likewise causes little distances between liquid particles.The effects of Soret parameter 0.1 ≤ N F T ≤ 0.6 and Dufour parameter 0.1 ≤ N TF ≤ 0.6 on the three main distributions of temperature, solutal concentration, and Fig.16 (2D and 3D).Physically, the presence of a drag force that prevents the development of momentum and results in a decrease in the flow inside the channel causes friction between the layers of the viscous fluid and the walls of the channel, resulting in an increase in the fluid's temperature, this is the primary reason for the improvement in the temperature distribution under the influence of the magnetic field parameter 0.5 ≤ M ≤ 3 .Like- wise, the reaction of nanoparticles and natural liquid particles inside the liquid to the strength of the attractive field made their fixations decline.It turned out, dear reader, that the effect is positive on the temperature distribution θ and negative on the other two distributions γ , � when the heat generation/absorption parameter −0.4 ≤ β ≤ 0.1 was plotted in Fig.17(2D and 3D) specifically to explain its effect on the three main distribu- tions temperature θ , solutal concentration γ , and nanoparticles volume fraction .Physically, the expansion in the field of temperature dissemination is communicated from heat absorption (β < 0) to heat generation (β > 0) by the normal active energy of the liquid atoms coming about because of their retention of thermal power, which prompts an expansion in the development and speed of the fluid particles inside in light of the fact that the temperature decides the movement energy related with the development of liquid particles and nanoparticles which likewise causes little distances between liquid particles.The effects of Soret parameter 0.1 ≤ N F T ≤ 0.6 and Dufour parameter 0.1 ≤ N TF ≤ 0.6 on the three main distributions of temperature, solutal concentration, and

Figure 14 .
Figure 14.Effect of Bn on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 15 .
Figure 15.Effect of Nt on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 16 .
Figure 16.Effect of M on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 17 .
Figure 17.Effect of β on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

23 and 24 .
The large response to the pressure gradient contributed to the increase in the rate of fluid pumping into the channel.Finally, Fig.25presents the features of the effect of the mean flow coefficient on the distribution of pressure gradients.When looking at the figure, it was noticed that the large response range and the positive direct relationship between this coefficient and this distribution, and this is clear from the remarkable spacing between the curves, if this indicates something, it indicates that both of them help to support the fluid flow inside the flexible channel and help in the regularity of the fluid movement, especially in overcoming the narrow areas

Figure 18 .
Figure 18.Effect of N FT on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 19 .
Figure 19.Effect of N TF on the distributions of temperature, solutal concentration, and nanoparticles volume fraction (2D and 3D).

Figure 20 .
Figure 20.Effect of Nb on the distributions of temperature and nanoparticles volume fraction (2D and 3D).

Figure 21 .
Figure 21.Pressure gradient distribution for different values of Ω.
2 , electric Reynolds number 0.1 ≤ R m ≤ 1.1 and finally the strength magnetic field 0.1 ≤ H o ≤ 0.6 , this behavior appears in Figs. 28, 29, and 30, respectively, it becomes clear to you dear reader, from the first glance, that the distribution of the induced magnetic field h x becomes a decreasing function in the left side of the channel under the influence of the three parameters, respectively, in the interval −1.5 ≤ y ≤ 0 for E and R m parameters in Figs.28 and 29, and interval −1.5 ≤ y ≤ −0.5 concerning parameter H o in Fig. 30, and soon a rapid change occurs in the behavior of the induced magnetic field, it becomes an increasing function in the right side of the channel in distance 0 ≤ y ≤ 1.5 for the first and second parameters ( E and R m ) and in distance −0.5 ≤ y ≤ 1.0 for the third parameter (H o ).

Figure 22 .
Figure 22.Pressure gradient distribution for different values of Da.

Figure 23 .Figure 26 .
Figure 23.Pressure gradient distribution for different values of M.

Figure 27 .
Figure 27.Magnetic force distribution for different values of R m .

Figure 28 .
Figure 28.Induced magnetic field distribution for different values of E.

Figure 29 .
Figure 29.Induced magnetic field distribution for different values of R m .

Figure 30 .
Figure 30.Induced magnetic field for different values of H o .

Figure 31 .
Figure 31.Streamlines profiles for different values of Nb.

Figure 32 .
Figure 32.Streamlines profiles for different values of Ω.

Figure 33 .
Figure 33.Streamlines profiles for different values of Da.

Figure 34 .
Figure 34.Streamlines profiles for different values of Nt.

Figure 35 .
Figure 35.Streamlines profiles for different values of Gr t .

Table 1 .
A comparison between some previous works. Nano-