Darcy–Brinkman Model for Ternary Dusty Nanofluid Flow across Stretching/Shrinking Surface with Suction/Injection

: Understanding of dusty fluids for different Brinkman numbers in porous media is limited. This study examines the Darcy–Brinkman model for two-dimensional magneto-hydrodynamic fluid flow across permeable stretching/shrinking surfaces with heat transfer. Water was considered as a conventional base fluid in which the copper ( Cu ), silver ( Ag ), and titanium dioxide ( TiO 2 ) nanoparticles were submerged in a preparation of a ternary dusty nanofluid. The governing nonlinear partial differential equations are converted to ordinary differential equations through suitable similarity conversions. Under radiation and mass transpiration, analytical solutions for stretching sheets/shrinking sheets are obtained. Several parameters are investigated, including the magnetic field, Darcy–Brinkman model, solution domain, and inverse Darcy number. The outcomes of the present article reveal that increasing the Brinkman number and inverse Darcy number decreases the velocity of the fluid and dusty phase. Increasing the magnetic field decreases the momentum of the boundary layer. Ternary dusty nanofluids have significantly improved the heat transmission process for manufacturing with applications in engineering, and biological and physical sciences. The findings of this study demonstrate that the ternary nanofluid phase’s heat and mass transpiration performance is better than the dusty phase’s performance.


Introduction
The Darcy-Brinkman model with magneto-hydrodynamics (MHD) has many applications in paint sparing, research on air pollution, capillary blood flow, and vehicle and factory emissions of smoke and other pollutants; dusty fluids contain fine dust particles.The laminar boundary layer flow of an electrically conducting dusty fluid across a stretching/shrinking sheet has drawn the attention of many researchers.For example, Saffaman et al. [1] studied stability solutions of dusty fluid flow.Modeling the flow and heat of a two-phase dusty fluid across deteriorating isothermal boundaries was investigated by Turkyilmazoglu et al. [2].Farooq et al. [3] investigated the influence of thermal radiation on the laminar fluid flow of dust-containing ternary nanofluid across permeable stretched sheets.Jalil et al. [4] and Datta et al. [5] explored an accurate solution to the dusty fluid flow problems in the MHD boundary layer across a shrunken sheet with fluid flowing across an infinitely flat surface.
Crane et al. [6] studied flow over a stretching sheet.Flow on a continuously stretched sheet with heat transfer was investigated by Carragher et al. [7].A reduction in temperature over a stretched sheet with constant heat flux was studied by Dutta et al. [8].In another study, Pradhan et al. [9] investigated the Marangoni convection and radiation effects on ternary nanofluid flow in a permeable media with mass transfer.An investigation on the heat transfer and viscous dissipation over stretched sheets and around the stagnation point was carried out by Aly et al. [10].Sheu et al. [11] explored a layer of porous material saturated with a nanofluid that exhibited thermal instability.
The Brinkman ratio is a crucial parameter in controlling heat conduction.Increasing the Brinkman parameter causes the temperature to rise because of the slower heat transfer in viscous dissipation.Many researchers have explored the effect of the Brinkman ratio on fluid flow; for example, Pop et al. [12] studied the Brinkman concept for flow across a circular cylinder contained in a permeable material.The thermal conveyance model for sinusoidally oscillating Brinkman-type nanofluids containing carbon nanotubes was examined by Jie et al. [13].Hsu et al. [14] studied a Brinkman model for spontaneous convection in a permeable material around a vertical flat surface with a semi-infinite height.This work was extended by Sachhin et al. [15] and Adun et al. [16], exploring the Brinkman model's effect on heat-transfer-enhanced porous stretched sheets with non-Newtonian fluid flow, and they also reviewed the synthesis stability and thermophysical characteristics.
Sarwar et al. [17] studied ternary nanoparticles in a magneto-hydrodynamics non-Newtonian fluid flow for heat and mass transpiration.The effects of radiation on a ternary nanofluid flow in a porous media with heat and mass transfer were examined by Mahabaleshwar et al. [18], while Sahoo et al. [19] and Adnan et al. [20] examined the novel correlation to measure the viscosity of ternary hybrid nanofluid and analyzed the heat transfer efficiency for TNF flow in a radiated channel under various conditions.Nagamgari et al. [21] studied the volume percentage of dusty particles with nanofluid flow in magneto-hydrodynamics across a stretched surface, with Sneha et al. [22] carrying out an investigation on the Darcy-Brinkman equation for heat transport and mass transpiration in a hybrid dusty nanofluid flow, while thermal management of curved solid conductive panels linked with various nanofluid cooling methods was investigated by Selimefendigil et al. [23].
Hayat et al. [24] and Roy et al. [25] studied copper and silver water nanofluids using nonlinear thermal radiation and magneto-hydrodynamics dusty hybrid nanofluid flows and heat transmission across a contracting sheet.Furthermore, the dynamics of a heated convective surface with magnetic flux density and partial slip in asymmetric flow over a stretched sheet was examined by Animasun et al. [26].Ariel et al. [27] and Jalili et al. [28] investigated the heat transmission in a polar system in a cylinder-shaped magnetic field.Nawaz et al. [29] studied the asymmetric MHD flow of a Casson fluid with varying thermal conductivity.Following the above work, Sachhin et al. [30] investigated the analytical solution of a Casson fluid flow across a stretched surface.Finally, several other studies regarding stretching surfaces under thermal radiation and mass transfer have been published [31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46].
Although the extrusion and cooling of dusty-fluid processes have various practical and scientific applications, the mass transpiration in the boundary layer of ternary dusty nanofluids has yet to be adequately studied as an analytical technique applicable to permeable media.The originality of this investigation concerns the influence of radiation and a heat source/sink on ternary dusty nanofluid water flow with heat transfer across permeable stretching/shrinking surfaces.The motivation is to combine dusty particles with the ternary nanofluid in the framework of momentum and temperature equations to derive an analytic solution.The analytic solution is important for supporting nanofluid computational research [47][48][49][50][51].The non-dimensional governing equations are obtained using suitable similarity transformations over various boundaries [52][53][54].The results of the current study impact several scientific and industrial applications.
Considering the findings and research gaps identified by previous studies and the associated impact on applications, the magnetic field and the Darcy-Brinkman model on ternary dusty nanofluid flow across permeable stretching/shrinking surfaces are investigated.The results of this study contribute to improving the understanding of ternary nanoliquids under MHD and different Brinkman ratios.A similarity transformation is used to reduce the governing equations of the problem into a system of nonlinear ordinary differential equations.The paper presents the analytical solution of the momentum equations and discusses the results of the porous medium, MHD, and volume fraction.
The findings of this paper are applicable in the development of efficient fuel cells and the polymer industry, concerning stretching/shrinking sheets.

Mathematical Formulation
Consider the steady two-dimensional, incompressible, viscous, laminar flow of a ternary dusty nanofluid across a permeable stretching/shrinking surface.The surface is parallel to the x-axis and normal to the y-axis (Figure 1); the stretching/shrinking sheet is parallel to the y-axis, and the momentum of the fluid flow is assumed as U w (x) = dax.The added ternary nanoparticles are copper (Cu), silver (Ag), and titanium dioxide (TiO 2 ), with water as the base fluid.The Darcy-Brinkman model is examined.

•
The fluid flow is generated by stretching/shrinking the sheet, and there is no pressure gradient affecting the fluid, i.e., ∇p = 0. • Finally, steady flow is considered, i.e., ∂u ∂t = 0. • A stretching velocity, u = dax, is also introduced.

A. Governing equations
The boundary layer flow [2,3,22] is described by the following equations: a. Fluid phase: • Continuity equation: • Momentum equation: • Temperature equation: Dusty phase: • Continuity equation: • Momentum equation: • Temperature equation: The imposed boundary conditions (BCs) include the fluid and dusty phases [2,22], given as Here, u, v, u p , and v p are the velocity components of a fluid and dusty fluid phase along the x-and y-directions, respectively; the dusty-and fluid-phase temperatures are Tp and T; µ tn f is the dynamic viscosity; ρ tn f is the effective density; κ tn f is the thermal conductivity; σ tn f is the electrical conductivity; and c p and c m are the specific heat coefficients; τ T is the heat equilibrium; L 1 is the Stokes's drag/resistance; ν tn f is the kinematic viscosity of nanoparticles N; ρCp tn f is the heat capacitance; k 1 is the flow permeability; and τ v = m L 1 is a relaxation time parameter, where m denotes the mass of dusty particles.

B. Similarity variables introduced to convert governing PDEs to ODEs [2,22]
where F and f denote stream functions of the dusty and fluid phases, η denotes the similarity variable, a > 0 is a constant, and the wall mass transfer velocity is given by v w = − √ ν f aS.We employ the Rosseland approximation for the value of the radiative heat flux, given by [33] where σ * denotes the Stefan-Boltzmann constant, and κ * is the mean absorption.We expand the above equation by a Taylor expansion, neglecting the higher-order terms, as a linear function of temperature [33]: Equations ( 9) and ( 10) yield Combining Equations ( 11) and ( 3) gives Using the similarity transformation, Equations ( 2), ( 3), (5), and ( 6) can be written as c.Fluid phase: d. Dusty phase: The modified boundary conditions are [2,3] where , and ρ denotes the mass number of dust particles; ρ p = Nm denotes the particle-phase density, fluid-particle interaction parameter denoted as β = 1 aτ v .The Prandtl number is given by Pr = ν p c p ν .S is the mass suction/injection parameter; the particular heat coefficient is γ = c m c p ; the temperature of the fluid interaction component is β T = 1 aτ T ; the Eckert number is denoted as Ec = a 2 bC p ; the inverse Darcy number is ; and the heat source/sink parameter is given by Ni = Q 0 aρ f , where d denotes a stretching/shrinking parameter.

Analytic Solution for Velocity Equations for Fluid-and Dusty-Phase Flow
Using the above boundary conditions, the exponential forms of solutions lead to exact solutions to Equations ( 13) and ( 14) as follows [2,3]: where λ > 0, and α denotes the stretching speed of dust particles.Differentiating Equation ( 18) for the fluid phase, Equation ( 13), yields The momentum equation for the dusty phase is Solving Equation (22) gives we substitute Equation (21) into Equation (19); the cubic equation is obtained: The roots of the above algebraic equation are for simplification −2ξ 3 2 + 9ξ 1 ξ 2 ξ 3 − 27ξ 2 1 ξ 4 is defined by z, where The wall shear can be calculated from Equation (20) with the values of λ and α [2,3]:

Analytic Solution for Temperature Equations for Fluid-and Dusty-Phase Flow
By using boundary conditions, we can write specific forms for solutions of the fluid and dusty phases of temperature, Equations ( 15) and ( 16), as follows [2,3]: differentiating the above equations and substituting into the energy Equations ( 15) and ( 16) leads to the polynomial equations After simplifying Equation ( 29), the following value of φ is obtained: The efficient heat exchange of the fluid and dusty phase are obtained as follows as [2,3] dθ dη η=0 = 2λ, and

Results and Discussion
The present study focuses on the Darcy-Brinkman model in conjunction with a magnetic field over a two-phase ternary nanofluid with dusty particles flowing across stretching/shrinking boundaries and under the influence of thermal radiation and a heat source/sink.The governing equations are converted to coupled ODEs via similarity transformations, and exact solutions for the momentum and temperature equations are obtained for the fluid and dusty phases.The Prandtl number for water, Pr, is considered to be 6.2.The range of parameters used are as follows: Brinkman number 0.1 ⩽ Λ ⩽ 5, mass suction/injection 0.1 ⩽ S ⩽ 5, stretching/shrinking parameter −2 ⩽ d ⩽ 2, magnetic field parameter 0.5 ⩽ M ⩽ 10, inverse Darcy number 0.5 ⩽ Da −1 ⩽ 10, 0.1 ⩽ l ⩽ 5, and volume fraction 0.1 ⩽ ϕ < 0.3.The application of the model concerns the parameters shown in Table 1.Note that past investigations served as the limiting case of the present work:

•
Absence of heat source/sink, and the presence of hybrid nanoparticles: limiting case is the results of Sneha et al. [22].

•
Absence of magnetic field, heat source/sink, Brinkman ratio: limiting case is the results of Farooq et al. [3].

•
Absence of magnetic field, heat source/sink, Brinkman ratio, volume fraction: limiting case is the results of Turkyilmazoglu et al. [2].
In the following sections, we will show the results of the analytic solution for a range of parameters.

e.
The ternary nanoparticles' thermophysical properties [17,37,39] Viscosity of the ternary nanofluids stated as [17]: Density of the ternary nanofluids stated as [17,39]: Thermal capacity of the ternary nanofluids stated as [37]: The nanoparticles' effects on the thermal conductivities in a ternary nanofluid [17] are described by the following set of equations: The electrical conductivities are given by f. Thermophysical properties of ternary nanofluid [17]: Figures 2 and 3 show the solution concerning mass transpiration for a magnetic field of 0.5, an inverse Darcy parameter of 0.1, a Brinkman number of 0.5, and various values of β.The results are obtained for two different values of mass number l: mass number l = 0 (dotted lines) and mass number l = 3 (solid lines).Black solid and dotted lines correspond to λ 1 ; red solid and dotted lines correspond to λ 2 ; purple solid and dotted lines correspond to λ 3 roots.All three solutions are in the positive x-axis.The free-stream velocity remains the same on the stretching surface and differs on the shrinking surface.The boundary layer increases as the β increases, which leads to the relaxation time becoming very small, and increasing β leads to a greater fluid-phase thickness than the thickness related to the dusty phase.7, which elucidates the effect of the higher value of a magnetic field on the solution β for both stretching/shrinking boundaries.The imposed magnetic field produces a kind of drag or resistive force and encourages the resistive Lorentz force which is present in an electrically conducting fluid.This force results in resistance to the fluid particle's momentum, which raises the fluid-phase stretching sheet, significantly retards the transport phenomenon, and consequently decreases the boundary layer thickness.Figures 8-10 show the dusty-and fluid-phase velocity profile plots concerning similarity variables, where the red lines denote the dusty phase and the black lines represent the fluid phase for a stretching boundary.Figure 8 corresponds to inverse Darcy numbers 0, 5.This condition will restrict the flow of both phases.Due to the permeable factor, which imposes frictional and drag forces on the liquid, the following values are used for the various parameters: the magnetic field is 1, the mass suction is 2, the Brinkman number is 0.1, the mass number is 2, and the fluid particle interaction parameter is 2. When increasing the inverse Darcy number, both fluid phases and the velocity profile decrease, and the fluid phase is more dominant than the dusty phase.The momentum of the ternary nanofluid phase and dust phase is decreased.Only one solution is obtained.As shown in Figure 9, for variation in the magnetic field M = 1, 5, the magnetic field encourages the resistive Lorentz force which is present in an electrically conducting fluid.This force resists the fluid particle's momentum, which lowers the fluid's velocity in the stretching sheet.When increasing the magnetic field, both the fluid-and dusty-phase velocity profiles decrease, and the fluid phase is more dominant than the dusty phase; it reduces the momentum of TNF and the dust phase.The Brinkman parameter is the relation between drag force and density.Therefore, the drag force increases when increasing the values of the Brinkman number, which slows down the fluid velocity.A rising Brinkman parameter reduces the ternary fluid velocity and decreases the dusty-phase velocity.16 show the temperatures for the fluid and dusty phases for different similarity variables.Figure 14 shows the temperature graph of the fluid phase while keeping the inverse Darcy number at 1, the magnetic field at 0.1, and the heat sink at −2.The fluid phase temperature increases for the stretching boundary compared to the dusty-phase temperature.Figure 15 shows the dusty-phase's temperature for a stretching boundary while keeping the inverse Darcy number at 1, Brinkman number at 2.1, magnetic field at 0.1, and heat sink at −2.The dusty-phase temperature boundary layer rapidly increases compared to the fluid-phase temperature boundary layer.Figure 16 shows the temperature of both the fluid-and dusty-phases for a stretching boundary for Darcy number of 1, Brinkman number of 2.1, magnetic field of 0.1, and heat source of 2.1.The fluid-phase temperature boundary layer significantly increases compared to the dusty-phase temperature boundary layer.The solution parameter and mass number enhance the thermal profile of both the fluid and dusty phases, and the temperature of the fluid phase increases with increasing the fluid-particle interaction parameter compared to the dusty-particle interaction.Figures 17 and 18 show the temperature and velocity graphs for the dusty phase concerning the similarity variable.Figure 17 shows the temperature of the dusty phase for a shrinking boundary for the inverse Darcy number of 1, magnetic field of 0.1, and heat sink of −2.The fluid-phase temperature increases dramatically, and the dusty fluid phase decreases the thermal boundary layer.Compared with the results of Figure 14, the dustyphase temperature drops for a shrinking boundary.Figure 16 shows the velocity of the dusty phase for the shrinking boundary with varying inverse Darcy number = 1, 5. It will resist the fluid flow (Figure 11).The porous media decreases the velocity profile.Furthermore, in the dusty fluid phase, the flow paths are restricted due to the improvement in the permeable factor.There is frictional and drag forces on the liquid for a Brinkman number of 2.1, magnetic field of 0.1, and heat sink of −2.The shrinking boundary fluid phase increases with decreasing the dusty phase, which reduces the velocity of the fluid flow by creating resistance for the flow by varying the porous media.Table 1.Comparison of the present study and related works by other authors [2,6,[47][48][49].

Previous Related Works Fluids Value of Momentum Solution
Turkyilmazoglu

Concluding Remarks
This analytical study examines the MHD's influence across a ternary dusty nanofluid with permeable media, radiation, and a heat source/sink also considered.Using suitable transformations, velocity and temperature equations are converted to a set of feasible ODEs and solved using an analytical method.Physical significance parameters like the solution domain, inverse Darcy number, and Brinkman parameter have been examined.Table 1 is the analytical validation, Tables 2 and 3 gives the numerical comparision with previous studies.The findings of this study are:

•
An analytical approach can be used to solve the magnetic field with the Darcy-Brinkman model.

•
The fluid-and dusty-phase velocity profiles decrease when increasing the inverse Darcy number.

•
The momentum boundary layer becomes thicker as the Brinkman number rises.

•
When increasing the solution domain thickness, the particle interaction parameters strength also increases.

•
The analysis yields a unique solution when considering mass flow suction.

•
Increasing the magnetic field decreases the velocity profile due to Lorentz's force.

•
The fluid-phase temperature in the boundary layer significantly increases compared to the dusty-phase temperature in the boundary layer.

Figure 1 .
Figure 1.Schematic diagram of the problem.Assumptions in the Mathematical Equations Describing the Physical Model • The fluid phase and nanoparticles are in the thermal equilibrium state.• Water is considered a conventional base fluid in which the copper (Cu), silver (Ag), and titanium dioxide (TiO 2 ) nanoparticles are submerged in the preparation of a ternary dusty nanofluid.• An inclined magnetic field is introduced.• No chemical reactions take place in the fluid layer.• There is negligible viscous dissipation.• The nanofluid is incompressible; Newtonian and laminar flow are considered.• The Darcy-Brinkman model is examined.• The fluid flow is generated by stretching/shrinking the sheet, and there is no pressure gradient affecting the fluid, i.e., ∇p = 0. • Finally, steady flow is considered, i.e., ∂u ∂t = 0. • A stretching velocity, u = dax, is also introduced.

5 d 5 Figure 2 .
Figure 2. Solution graph for stretching boundary with variation in β.The results are obtained for stretching the boundary (d = 2).

5 , 5 Figure 3 .
Figure 3. Solution graph for shrinking boundary with variation in β.The results are obtained for shrinking the boundary (d = −1.4).

Figures 4
Figures 4 and 5 show the solution concerning mass transpiration.The black solid and dotted lines correspond to λ 1 , red solid and dotted lines correspond to λ 2 , and purple solid and dotted lines correspond to λ 3 roots, with a magnetic field of 0.5, inverse Darcy parameter of 10, Brinkman number of 0.5, and variation in β.Figure 4 is obtained when stretching the boundary (d = 2), and Figure 5 when shrinking the boundary (d = −1.2) for two different values of mass number l, where the mass number is l = 0 (dotted lines) and l = 3 (solid lines); all three solutions are in the positive x-axis.Increasing the value of porous media and β leads to resisting the fluid flow and decreases the boundary layer thickness.The flow paths for the flow of the fluid are restricted due to improvement in the permeable factor, which exhibits frictional and drag forces on the fluid flow.Increasing the β values decreases the boundary layer thickness for stretching and shrinking boundaries.At specific points, the results reveal similar orientations for different values of the mass number.

Figure 4
Figures 4 and 5 show the solution concerning mass transpiration.The black solid and dotted lines correspond to λ 1 , red solid and dotted lines correspond to λ 2 , and purple solid and dotted lines correspond to λ 3 roots, with a magnetic field of 0.5, inverse Darcy parameter of 10, Brinkman number of 0.5, and variation in β.Figure 4 is obtained when stretching the boundary (d = 2), and Figure 5 when shrinking the boundary (d = −1.2) for two different values of mass number l, where the mass number is l = 0 (dotted lines) and l = 3 (solid lines); all three solutions are in the positive x-axis.Increasing the value of porous media and β leads to resisting the fluid flow and decreases the boundary layer thickness.The flow paths for the flow of the fluid are restricted due to improvement in the permeable factor, which exhibits frictional and drag forces on the fluid flow.Increasing the β values decreases the boundary layer thickness for stretching and shrinking boundaries.At specific points, the results reveal similar orientations for different values of the mass number.

Figures 6
Figures 6 and 7 show the solution concerning mass transpiration, a magnetic field of 10, inverse Darcy parameter of 1, Brinkman number of 0.5, and variation in β.Figure 6 concerns the stretching boundary, d = 2, and Figure 7 the shrinking, d = −1.2, for two different values of the mass number l: 0 (dotted lines) and 3 (solid lines).Similar notation to Figure 6 applies to Figure7, which elucidates the effect of the higher value of a magnetic field on the solution β for both stretching/shrinking boundaries.The imposed magnetic field produces a kind of drag or resistive force and encourages the resistive Lorentz force which is present in an electrically conducting fluid.This force results in resistance to the fluid particle's momentum, which raises the fluid-phase stretching sheet, significantly retards the transport phenomenon, and consequently decreases the boundary layer thickness.
Figures 6 and 7 show the solution concerning mass transpiration, a magnetic field of 10, inverse Darcy parameter of 1, Brinkman number of 0.5, and variation in β.Figure 6 concerns the stretching boundary, d = 2, and Figure 7 the shrinking, d = −1.2, for two different values of the mass number l: 0 (dotted lines) and 3 (solid lines).Similar notation to Figure 6 applies to Figure7, which elucidates the effect of the higher value of a magnetic field on the solution β for both stretching/shrinking boundaries.The imposed magnetic field produces a kind of drag or resistive force and encourages the resistive Lorentz force which is present in an electrically conducting fluid.This force results in resistance to the fluid particle's momentum, which raises the fluid-phase stretching sheet, significantly retards the transport phenomenon, and consequently decreases the boundary layer thickness.

Figure 6 .
Figure 6.Solution for stretching boundary for M = 10 with variation in β.

Figure 7 .
Figure 7. Solution for shrinking boundary for M = 10 with variation in β.

Figure 9 .
Figure 9. Velocity profile versus similarity variable M.

Figures 11 -
Figures 11-13 show the axial velocity profiles for the inverse Darcy number 1, 5 values that resist the fluid flow.The flow paths for the fluid in both phases are restricted due to improvement in the permeable factor, which exhibits frictional and drag forces on the liquid.The results include mass suction and injection values and a shrinking boundary while keeping the magnetic field at 0.5, mass transpiration values of 2, −2, and inverse Darcy number of 0.1.The momentum profile in both the fluid and dusty phases decreases with upsurges in the inverse Darcy number compared to the graph in Figure8, which shows almost similar phenomena for stretching and shrinking boundaries.Figure8shows nearly similar phenomena for both stretching and shrinking boundaries.Figure13shows the Brinkman parameter for the shrinking boundary and mass suction values.The drag force increases with increasing Brinkman number values, slowing the fluid velocity.Increasing the Brinkman parameter decreases the velocity of the ternary nanofluid phase, and the dusty phase slowly decays compared to the liquid phase.

Figure 8
Figures 11-13 show the axial velocity profiles for the inverse Darcy number 1, 5 values that resist the fluid flow.The flow paths for the fluid in both phases are restricted due to improvement in the permeable factor, which exhibits frictional and drag forces on the liquid.The results include mass suction and injection values and a shrinking boundary while keeping the magnetic field at 0.5, mass transpiration values of 2, −2, and inverse Darcy number of 0.1.The momentum profile in both the fluid and dusty phases decreases with upsurges in the inverse Darcy number compared to the graph in Figure8, which shows almost similar phenomena for stretching and shrinking boundaries.Figure8shows nearly similar phenomena for both stretching and shrinking boundaries.Figure13shows the Brinkman parameter for the shrinking boundary and mass suction values.The drag force increases with increasing Brinkman number values, slowing the fluid velocity.Increasing the Brinkman parameter decreases the velocity of the ternary nanofluid phase, and the dusty phase slowly decays compared to the liquid phase.
Figures 11-13 show the axial velocity profiles for the inverse Darcy number 1, 5 values that resist the fluid flow.The flow paths for the fluid in both phases are restricted due to improvement in the permeable factor, which exhibits frictional and drag forces on the liquid.The results include mass suction and injection values and a shrinking boundary while keeping the magnetic field at 0.5, mass transpiration values of 2, −2, and inverse Darcy number of 0.1.The momentum profile in both the fluid and dusty phases decreases with upsurges in the inverse Darcy number compared to the graph in Figure8, which shows almost similar phenomena for stretching and shrinking boundaries.Figure8shows nearly similar phenomena for both stretching and shrinking boundaries.Figure13shows the Brinkman parameter for the shrinking boundary and mass suction values.The drag force increases with increasing Brinkman number values, slowing the fluid velocity.Increasing the Brinkman parameter decreases the velocity of the ternary nanofluid phase, and the dusty phase slowly decays compared to the liquid phase.

Figures 14 -
show the temperatures for the fluid and dusty phases for different similarity variables.Figure14shows the temperature graph of the fluid phase while keeping the inverse Darcy number at 1, the magnetic field at 0.1, and the heat sink at −2.The fluid phase temperature increases for the stretching boundary compared to the dusty-phase temperature.Figure15shows the dusty-phase's temperature for a stretching boundary while keeping the inverse Darcy number at 1, Brinkman number at 2.1, magnetic field at 0.1, and heat sink at −2.The dusty-phase temperature boundary layer rapidly increases compared to the fluid-phase temperature boundary layer.Figure16shows the temperature of both the fluid-and dusty-phases for a stretching boundary for Darcy number of 1, Brinkman number of 2.1, magnetic field of 0.1, and heat source of 2.1.The fluid-phase temperature boundary layer significantly increases compared to the dusty-phase temperature boundary layer.The solution parameter and mass number enhance the thermal profile of both the fluid and dusty phases, and the temperature of the fluid phase increases with increasing the fluid-particle interaction parameter compared to the dusty-particle interaction.

Figure 14 .
Figure 14.Temperature profile for fluid phase versus similarity variable.

Figure 15 .
Figure 15.Temperature profile for dusty phase versus similarity variable.

Figure 16 .
Figure 16.Temperature profile for dusty and fluid phases versus similarity variable.

Figure 17 .
Figure 17.Temperature profile versus similarity variable for a shrinking boundary.