Abstract

This work aims to analyze the Casson thermal boundary layer flow over an expanding wedge in a porous medium with convective boundary conditions and ohmic heating. Moreover, the effects of porosity and viscous dissipation are studied in detail and included in the analysis. The importance of this study is due to its applications in biomedical engineering where the analysis of behavior of non-Newtonian blood flow in arteries and veins is desired. Within the context of blood flow, it is also applicable to many other fields, for instance, radiative therapy, MHD generators, soil machines, melt-spinning, and insulation processes. The modeled problem is a set of PDEs, which is nondimensionalized to derive a nonlinear boundary value problem (BVP). The obtained BVP is solved using the shooting technique, endowed with the order four Runge-Kutta and Newton methods. The impact of different parameters on the momentum and temperature fields is investigated along with two important parameters of physical significance, i.e., the Nusselt number and the surface drag force. Results are validated, and an excellent agreement is seen for the parameters of interest using MATLAB built-in function bvp4c. A significant finding is that by increasing the Casson liquid parameter, the velocity decreases as the wedge expands quicker than the free stream velocity at . However, the velocity increases for the case when . A decrease in the Darcy number increases the temperature profile. Furthermore, the convective parameter accelerated the heat transmission rate, and a rise in the Prandtl number thickens the thermal boundary layer. The findings of this investigation contribute to problems in fluid dynamics and heat transfer that involve studying the behavior of a non-Newtonian fluid with Casson rheological properties near a solid porous wedge surface.

1. Introduction

Mixed convection across an impermeable or permeable wedge is used in many manufacturing processes, including extruding molten polymers, fabricating plastic sheets for solar energy, and storing thermal energy. According to Shah et al. [1], mixed convection is a fluid flow phenomenon characterized by the simultaneous effect of forced convection, which is caused by external forces such as a pump or a fan, and natural convection, which is caused by buoyant forces caused by temperature gradients within the fluid. Systems such as technical equipment and natural phenomena such as the movement of temperature stratified masses of air and water on Earth exhibit mixed convection flow [2, 3]. It has uses in heat exchangers [4], thermal design, cooling electronic components, and other fields. It is also pertinent to heat transfer and fluid flow analysis. Karwe and Jaluria [5, 6] explained that the investigation is necessary to check the quality of resulting products and cooling rates. Ishak et al. [2] investigated the mixed convection magnetohydrodynamics boundary layer across a vertically expanding sheet while the wall temperature remained unchanged. The reader is referred to [1, 7, 8] for a good overview on the mixed convection MHD viscoelastic fluid flow analysis.

Some recent studies have also considered the mixed convective flow in analyzing the porous flow structures. In this connection, thermal analysis in a higher grade Forchheimer porous nanostructure bounded between nonisothermal plates is investigated by Saleem et al. [9]. They utilized a finite difference approach in their analysis. The transient problem of mixed convection over a vertically moving wedge was examined by Ravindran et al. [10]. Numerical results show that the buoyant force enhances the skin friction and Nusselt number. The thermal properties of magnetohydrodynamics of Falkner–Skan Casson flow sideways a moving wedge are analyzed by Ullah et al. [11]. They used finite difference and quasilinearization techniques to derive the problem’s solution. It is detected that an augmentation in the values of the Eckert number has a negligible impact on force convection flow. They discovered that the existence of the suction or injection changes the rates of skin friction, concentration, and heat transfer. Unsteady flow with mixed convection dusty flow across a vertical wedge was investigated by Hossain et al. [3].

Many of the above-described studies addressed the flow convection over a vertical porous wedge. It is worth mentioning that there is limited research about convection flow in the presence of suction/injection, viscous and ohmic dissipation upon a permeable/impermeable wedge. The numerical solution of thermal radiation influence on magnetohydrodynamics with forced convective flow adjacent to a wedge with a heat source/sink [12] was computed by Chamkha et al. [13] employing the implicit finite difference technique. Anwar Hossain et al. [14] have explored the problem of MHD flow via a wedge with variable surface temperature. Additionally, perturbation solutions for different values of dimensions of time are discovered. Su et al. [15] investigated the MHD mixed convective flow of heat transfer across a permeable stretched wedge using numerical and analytic (DTM) methods. Also, the authors demonstrated that as the domain becomes unbounded, the DTM-BF solutions diverge. The influence of chemical reaction on the mixed convection magnetohydrodynamics and heat transfer across a porous wedge was investigated by Kandasamy et al. [16]. This research concludes that in a mixed convective regime, there is an increase in the hydrodynamic layer. In contrast, the concentration and temperature boundary layers decrease by lowering the values of the suction parameter. The issue of MHD flow through a porous nonisothermal wedge was examined by Prasad et al. [17]. They observed that among all variables, the mixed convection parameter substantially impacts the surface drag force, thermal rate, normal velocity of the wall, and temperature profile of the MHD flow through a porous nonisothermal wedge. Transient Falkner–Skan liquid flow of Carreau microliquid via a still/moving wedge was demonstrated by Khan et al. [18].

Due to the numerous applications in radiative therapy, metal casting technology, and ceramic engineering, several studies [4, 19] have explored the influence of radiative effect at the boundary layer with the transfer of heat characteristics. In biomedical engineering, within the context of blood flow [20, 21], the Casson parameter has been considered in analyzing the behavior of non-Newtonian flow in arteries [2224]. For a detailed review of the use and applications of Casson nanofluid flow, the reader is also referred to the work of Nayak et al. [25] who have investigated the Casson nanofluid flow subjected to biomedicine applications and chemical reactions [2628]. For some other important applications related to heat transfer enhancement, the reader is referred to the work of [29] and references therein. Ullah et al. [30] have investigated the Casson liquid flow across a nonstationary wedge using the effect of thermophoresis and Brownian diffusion. The transportation of heat for homogeneous/heterogeneous Casson liquid in a porous medium was studied by Bilal et al. [31]. The velocity, temperature, and concentration of solutions for homogeneous and heterogeneous reactions are fundamental limitations of the basic subsidized flow parameters. Several researchers investigated the impact of heat radiation on different flow problems [3234] using finite element strategies.

Viscous dissipation can significantly affect the temperature distribution and overall energy balance of the system. Similarly, porous media have a large contact surface with fluids, which can significantly increase the heat transfer effect. The porous medium changes the flow field conditions, and the conduction heat transfer coefficient is usually higher than that of the fluid studied. The main aim of this study is to extend the work of Hussain et al. [35] by considering the following additional physical effects:(i)The Casson thermal boundary layer flow rheology is assumed(ii)The effects of porosity and viscous dissipation are taken into account

In this respect, two important parameters of physical interest arises are the viscous dissipation parameter and the porosity parameter . The specific goals and objectives of the study are to investigate the following research questions:(i)How does a porous wedge flow’s thermal reaction to Casson flow rheology behave?(ii)How does the porosity affect the velocity and thermal profiles in the presence of Casson rheology?(iii)How does the viscous dissipation affect the hydrodynamic profiles of the present rheology?(iv)How do the viscous dissipation and porosity of the Casson fluid affect the skin-friction coefficient and Nusselt number?

The next section presents the modeling aspects of the modeled problem which is first transformed to a boundary value problem in ordinary differential equations. The problem is then solved numerically by the shooting method, and obtained results are discussed.

2. Mathematical Modeling

This problem considers an electrically conducting thermal radiative 2D Casson fluid flow with mixed convection on a moving wedge as shown in Figure 1. The conservation equations for mass, momentum, and energy for viscous incompressible fluid in the presence of the viscous dissipation, Rosseland radiation, and porosity constitute the modeling description of the considered flow and are described respectively as [36]andwhere in equation (2), the symbol represents the stress tensor of the considered fluid, which is represented by

The rheological equation for the incompressible flow of the Casson fluid is given by [37]where are the components of the strain tensor defined as

In equation (5), is the product of the component of the deformation rate with itself, , is a critical value of this product based on the non-Newtonian model, and is the yield stress of the fluid. The Cartesian coordinates’ system is denoted by from which -direction is along the wedge and -direction is normal. The velocity of the wedge is represented by , and the velocity that is distant from the wedge is , which depends on , i.e., . The stretching wedge has suction and injection velocity, denoted by . During the heating of the wedge, a hot fluid with temperature produces a high heat transfer coefficient , where is the temperature of the fluid away from the wedge. The uniform magnetic field is denoted by , acting along -direction and equivalent to . The wedge angle is denoted by .

2.1. Dimensional Equations

Connected flow equations with boundary conditions are (see for reference [15])

The boundary conditions are described as follows:

For the physical description of these boundary conditions, the reader is referred to the work of Su et al. [15]. Here, , , , , , , and denote the kinematic viscosity, gravity field, coefficient of thermal expansion, conductivity of fluid due to electric current, fluid density, heat capacity, and thermal conductivity, respectively. The term here is radiative heat flux and is equal to , where and are known as the Stefan–Boltzman parameter and mean absorption coefficient, respectively. is the parameter of Casson fluid.

2.2. Similarity Transformations

Here, in this subsection, the similar similarity transformations are introduced. Similar transformations are employed [15] to simplify the set of governing equations of a complex physical problem by converting them dimensionless and making them similar in structure.

These are particularly useful for analyzing the flow and heat transfer phenomenon scales with changing parameters like Reynolds and Nusselt numbers. There are a lot of situations when one needs to introduce nonsimilar transformations [38] to convert this physical system into a nondimensional form. Such situations arise when nonstandard boundary conditions [3941] are introduced in the mathematical description of the flow problems.

2.3. Nondimensional Equations

Equations (8)–(13) are nondimensionalized using transformations in Section 2.2. The resulting equations take the form:along with the dimensionless BCs as

2.4. Dimensionless Parameters

The associated dimensionless parameters are

3. Quantities of Interest

3.1. Surface Drag Force

The dimensional form of the skin friction is as follows:where represents the viscosity of the liquid and is the liquid yield stress. The nondimensional form of surface drag force is given as

3.2. Local Nusselt Number

The local Nusselt parameter is described aswhere is the conductive heat flux from the surface of the wedge.

4. Solution Method through Shooting and bvp4c

The above nonlinear ODEs are now solved through the shooting method. This numerical technique is used for finding the solution to a boundary value problem. Using the following derivatives, the boundary value problem is transformed into a first-order initial value problem (see Table​ 1).

Using the above notations, equations (15) and (16) are transformed into the following fifteen ODEs:where and are the missing initial conditions. Choose the missing initial conditions that satisfy the above boundary conditions. The RK-4 method is applied to solve the first-order initial value problem. Newton’s numerical technique is used to refine the starting values. This iteration is repeated until the values meet the required accuracy, achieved by getting sufficiently close to the given boundary conditions. The stopping criteria for Newton’s method are given as follows:Here, .

The solution of the above model boundary value problem can also be calculated through the built-in MATLAB function bvp4c. This built-in function is suitable for solving higher-order boundary value problems in fluid mechanics. The criteria behind the bvp4c function are the three-stage Lobatto III-a formula, implemented by the finite difference algorithm [42]. This formula uses collocation, and the collocation polynomial offers a -continuous solution that is consistently fourth-order accurate across the integration interval of the problem. The residual of the continuous solution serves as the foundation for mesh selection and error control in the model problem. The integration interval is split into smaller intervals using the collocation technique and a mesh of points. The solver can arrive at a numerical solution by resolving a global system of algebraic equations caused by the boundary conditions and the collocation conditions placed on each subinterval. The solver then determines how much each subinterval’s numerical solution error is. The solver adjusts the mesh and restarts the procedure if the solution does not meet the tolerance conditions. In addition to an initial approximation of the solution at the mesh points, one must supply the initial mesh points.

5. Analysis of Results

This research is focused on analyzing the following aspects:(1)The steady nonlinear conducting radiative Casson flow through a penetrable wedge is examined with ohmic heating, thermal behavior, suction/injection, and convective BC.(2)The effects of porosity and viscous dissipation are considered.(3)Drag caused by the friction of a fluid against the surface and heat transfer due to pure conduction is examined by calculating the numerical values of the skin-friction coefficient and Nusselt number by varying the values of wedge angle and the magnetic parameter .(4)Thickness of the momentum boundary layer is examined for magnetic parameter , wedge range , Casson parameter , Prandtl number , and radiative heat transfer parameter .(5)Thermal boundary layer is analyzed through the numerical results by varying the Eckert number , Biot number Darcy number , Prandtl number , and radiative heat transfer parameter .

6. Results and Discussion

This section describes the outcomes of the model presented above through the discussed shooting technique. Table 2 shows the values of the heat rate parameter and drag force parameter for various values of the wedge angle and velocity ratio parameter . When is less than or more than 1, surface drag force increases as the angle increases. It is also seen that by an increase in , the temperature gradient coefficient of Nusselt number decreases. The results acquired by the shooting approach are verified by BVP4C as well. Here, the values of different parameters are fixed as , , , , , , , , , , and .

The numerical effects of drag force and Nusselt number for various values of the velocity ratio parameter R and the magnetic parameter M are shown in Table 3. The values of the Nusselt coefficient , surface drag force coefficient , and magnetic coefficient rise when is less than 1. When is more than 1, the skin-friction coefficient and both drop as the magnetic coefficient values increase. The outcomes obtained using the shooting approach are likewise met by bvp4c.

6.1. The Velocity Profiles

The impact of the Casson parameter on the velocity profile at is seen in Figure 2(a). Increases in steepen the velocity profile when is less than one, whereas increases in deplete the velocity profile as . For the case , Figure 2(b) shows an inverse relationship between the velocity profile and the magnetic parameter ; however, for , the relationship is the reverse. The Lorentz force that the magnetic field produces is what causes this phenomenon. For both situations of , Figure 2(c) shows an increase in velocity through an increase in the wedge’s angle. The velocity boundary layer thickens with an increase in the velocity ratio , as seen in Figure 3(a). The velocity boundary layer thickens with an increase in the velocity ratio , as seen in Figure 3. The effects of radiation parameter and Prandtl number on the velocity profile are shown by generating Figures 3(b) and 3(c) from the numerical results. The figures show that for both situations of R, there is a fall in the velocity profile due to cumulative values of the Prandtl number and decreasing values of the radiation parameter. Three key implications are shown in Figure 4. To demonstrate how the Biot number affects the velocity profile, Figure 4(a) is drawn. The ratio of the body’s internal conductive resistance to its surface-measured external convective resistance is known as the convection parameter . It is observed that temperature is improved when internal conductive resistance is increased rather than exterior convective resistance and thus raises the velocity profile . Figure 4(b) illustrates how Eckert number has a major influence. This figure shows that for , is a decreasing function of ; however, the behavior is reversed when . As can be seen in Figure 4(c), the depth of the momentum barrier layer decreases as the Darcy parameter decreases.

6.2. The Temperature Profile

The influence of the Casson coefficient on the temperature profile , at and , is shown in Figures 5(a) and 6. The temperature profile decreases in all scenarios when increases. Figure 5(b) illustrates how changing the magnetic parameter affects the momentum boundary layer’s behavior. The graph indicates that decreases for rising values of for and exhibits opposite behavior at . A direct relationship between the wedge angle and the temperature profile is seen in Figure 5(c). Three important effects on thermal profiles are depicted in Figure 7. Here, Figure 6(a) demonstrates how the thickness of the thermal boundary layer decreases as the radiation parameter decreases. The thermal profile decreases as the Prandtl number accumulates, as seen in Figure 6(b). The temperature profile decreases as the velocity ratio parameter increases, as seen in Figure 6(c). The temperature’s relationship with the Eckert number is seen in Figure 8(a). When rises, the temperature profile decreases. The effect of convective number is seen in Figure 8(b). Increasing levels increases the body’s internal conductive resistance. The behavior of temperature is steepened. Figure 8(c) illustrates how increases along with a drop in the Darcy number . The following fixed values apply to the selection of various parameters: , , , , , , , , , , and

6.3. Nusselt Number and Surface Drag Force Coefficient

The relationship between the Nusselt number and and is evident in Figure 7(a). On the other hand, Figure 7(b) shows that skin friction rises with increasing values but decreases as increases. The Nusselt number is an increasing function of , as Figure 7(c) illustrates, but for a constant value of , it decreases with decreasing values of the Casson parameter . Furthermore, Figure 7(d) shows that skin friction rises with the increasing values of for a given value of .

7. Conclusions

This paper considers magnetohydrodynamics thermal flow and boundary layer analysis using Casson fluid over an extending wedge by considering porosity, ohmic heating, convective boundary layer conditions, and viscous dissipation. The mathematical model of the problem is a system of nonlinear PDEs. This system is transformed into a set of ODEs with suitable transformations. The obtained boundary value problem is numerically solved by using the shooting method. Tables and graphs give a detailed picture of the different thermodynamic profiles.

The following observations are of notable significance:(i)As Casson parameter increases when , increases and increase in when , reduces the depth of momentum boundary layer.(ii)When rises with , increases and when , decreases.(iii)At and , the momentum boundary layer decreases by an increase in .(iv)As the velocity ratio parameter increases, also increases.(v) increases as rises at , and decreases when .(vi)At and , the thermal boundary’s thickness is directly related to the angle parameter .(vii)It is also detected that the temperature profile reduces with increasing values of the Eckert number .

With a rise in the Prandtl number , the depth of the thermal boundary layer decreases.(i)When the value of rises, decreases.(ii)Decrease in the Darcy number increases the temperature profile .(iii)Nusselt number and drag force are the increasing functions of Biot number .(iv)Increasing the Casson parameter enhances the skin friction and reduces the Nusselt number.

The findings of this investigation contribute to problems in fluid dynamics and heat transfer that involve studying the behavior of a non-Newtonian fluid with Casson rheological properties near a solid porous wedge surface. Such problems are common in various engineering and scientific applications, including chemical engineering, biomedical engineering, and materials processing. In biomedical engineering, within the context of blood flow, the Casson parameter considered in this analysis describes the behavior of non-Newtonian blood flow in arteries and veins. In the pharmaceutical industry, the Casson parameter characterizes the rheological properties of pharmaceutical suspensions and emulsions, thus playing a vital role in the formulation and manufacturing of drugs. In polymer suspensions, the Casson parameter is vital in extrusion and injection molding processes. The findings of this investigation and the data obtained through numerical calculations delineate the variational changes in the physical parameters of importance, i.e., the skin-friction coefficient and Nusselt number. Therefore, it is recommended to reinvestigate some related physical problems in the field of applications as mentioned earlier based on the findings of the present investigation.

Nomenclature

: and components of velocity
:Cartesian coordinates
:Temperature
:Temperature of hot fluid
:Velocity far away from wedge
:Velocity of expanding wedge
:Suction/injection velocity
:Uniform magnetic field
[]:Total angle of the wedge
[]:Wedge’s angle
:Heat transfer coefficient
:Gravity field
:Conductivity of fluid due to electric current
:Radiative heat flux
:Stream function temperature
:Specific heat capacity
[]:Casson parameter
:Drag coefficient
[]:Skin friction
[]:Nusselt number
[]:Magnetic parameter
:Mixed convective parameter
[]:Grashof number
[]:Reynold number
:Radiative heat transfer parameter
:Eckert number
[]:Prandtl number
[]:Biot number
[]:Velocity ratio parameter
[]:Darcy inverse
:Local inertia parameter
Greek Symbols
:Thermal conductivity
[]:Porosity parameter
[]:Thermal expansion parameter
:Mass density
:Kinematic viscosity
[]:Dimensionless temperature.

Data Availability

No data were used to support the findings of this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Dur-e-Shehwar Sagheer conducted conceptualization, methodology, handled software, formal analysis, and writing of the original draft. M. Sabeel Khan and M. Asif Memon performed writing of the original draft, collected resources, data curation, investigation, visualization, and validation. R. Shehzad and Amsalu Fenta conducted project administration, writing of the review, editing, and supervision. Mohammad Alqudah conducted validation, investigation, physical interpretation, writing of the review, and editing. Nawal A. Alshehri has contributed in the final drafting of the revised version of the paper by reviewing it critically for important intellectual content, editing, and proofreading. However, he has also addressed a few points of the revision.

Acknowledgments

The researchers would like to acknowledge the Deanship of Scientific Research, Taif University, Saudi Arabia, for funding this work.