Abstract

The Darcy-Forchheimer flow of a Williamson fluid over a Riga plate was analyzed in this paper. Energy and mass equations are modeled with Cattaneo-Christov theory and double stratifications. The governing PDE models are altered into ODE models. These models are numerically solved by MATLAB bvp4c and analytically solved by the homotopy analysis method. The impact of governing flow parameters on fluid velocity, fluid temperature, fluid concentration, skin-friction coefficient, local Nusselt number, and local Sherwood number is scrutinized via graphs and tables. We acknowledged that the speed of the fluid becomes diminishes for more presence of porosity parameter. Also, we noted that the thermal and solutal boundary layer thicknesses are waning due to their corresponding stratification parameters. In addition, the maximum decreasing percentage of skin friction is obtained when the suction/injection parameter varies from 0.0 to 0.4 for Williamson and viscous fluids. The maximum increasing percentage of local Nusselt number occurs when the suction/injection parameter varies from 0.4 to 0.8 for Williamson and viscous fluids.

1. Introduction

Non-Newtonian fluids are extensively implemented in diverse industrial processes such as petroleum drilling, drawing of plastic films, fibre spinning, and food production. The Williamson fluid model is one of the simplest non-Newtonian models to replicate the viscoelastic shear-thinning attributes, see Williamson [1]. The flow of thermally radiative Williamson fluid on a stretching sheet with chemical reaction was disclosed by Krishnamurthy et al. [2]. They proved the fluid temperature falling off due to the presence of the Williamson parameter. Khan et al. [3] demonstrated the impact of slip flow of Williamson nanofluid in a porous medium. They exposed that the surface drag force suppresses due to rising the Williamson fluid parameter. The 2D unsteady radiative Williamson fluid flow on a permeable stretching surface was deliberated by Hayat et al. [4]. They noticed that the fluid speed becomes slow when the Williamson parameter is high. Nadeem et al. [5] examined the Williamson fluid flow past a stretching sheet, and they found that the skin friction coefficient decreases with enhancing the Williamson parameter. Make use of the Keller box procedure to solve the problem of MHD flow of Williamson fluid over a stretching sheet by Salahuddin et al. [6]. Their outcome shows that the Williamson fluid parameter leads to suppress the fluid velocity. Few significant analysis for this area is seen in Refs. [7, 8].

Fluid flow over a porous medium is confronted in plentiful applications in industry. Few applications are wood drying, nuclear waste storage, food processing, oil purifying, drainage, and irrigation. Darcy’s principle is applied to analyze the flow behavior under the condition of small velocity and low porosity. When the quantity of Reynolds number overcomes unity, the Darcy principle was not applicable. Forchheimer [9] defeated this limitation by inserting the square velocity term in the momentum equation. After that, this is known as the Forchheimer number, which is applicable for working higher Reynolds number. Numerical analysis for a Darcy-Forchheimer flow of viscous fluid over a plate was inspected by Mukhopadhyay et al. [10]. They noted that the permeability parameter leads to a decrease in the warmth of the fluid. Hayat et al. [11] demonstrate the 3D Williamson nanomaterial flow on a Darcy-Forchheimer porous medium. They concluded that the surface shear stress diminishes for growing the Forchheimer number. The Darcy-Forchheimer flow of a viscous fluid with heterogeneous-homogeneous chemical reactions was portrayed by Khan et al. [12]. Their results clearly show that the fluid speed becomes slowdown due to the availability of Darcy number. Haider et al. [13] scrutinize the Darcy-Forchheimer and slip flow of hybrid nanofluid on a rotating disk. They proved that the larger estimation of Forchheimer enhances the fluid temperature. Steady 3D Darcy-Forchheimer flow of carbon nanotubes on a rotating disk was revealed by Sadiq et al. [14]. Some important studies for these concepts are collected in Refs. [1518].

The magnetic field plays a significant role in the development of fluid thermophysical traits. The demeanour of broadly used fluids like liquid metals, plasma, and electrolytes has a low conductor of electricity. Therefore, an external agent is required to boost up the heat transfer attributes through superior conductivity and thermophysical traits. A magnetic bar with permanently fixed magnets and alternate electronics, known as a Riga plate, can be acted as an external agent to improve fluid electricity. This plate was introduced by Gailitis and Lielausis [19]. Nanofluid flow over a Riga plate was deliberated by Ahmad et al. [20]. Nazeer et al. [21] inspected the chemically reacting Eyring-Powell nanofluid on a Riga plate. They proved that the fluid speed enhances when enhancing the modified Hartmann number. Chemically reacting Prandtl fluid on a Riga plate was addressed by Gireesha et al. [22]. Their results show that the velocity boundary layer thickens due to the more presence of the modified Hartmann number. Mehmood et al. [23] performed the impact of Soret and Dufour effects of a Casson fluid flow on a Riga plate with chemical reaction. Ayub et al. [24], Nayak et al. [25], and Rasool et al. [26] are few essential studies of fluid flow over a Riga plate.

Stratification is a natural process that combines two or more fluids with different densities, temperatures, and concentrations. The double stratification occurs due to both the heat and mass transfer differences. Cheng [27] examined the mass and heat transfer analysis of a power-law fluid in a stratified medium. He noticed that the heat transfer gradient declines for escalating the thermal stratification parameter. The radiative flow of a hyperbolic tangent fluid with chemical reaction and dual stratification’s in a porous medium was elucidated by Khan et al. [28]. They found that the fluid concentration downturns for the high magnitude of the stratified thermal parameter. Rehman et al. [29] evaluated the problem of a chemically reacting Williamson fluid with dual stratification, and they have seen that the rate of heat transfer rate is declined for the presence of thermal stratification parameter. The impact of solutal and thermal stratification of a Williamson nanofluid was deliberated by Khan et al. [30]. It is noticed that the horizontal velocity suppresses when the higher magnitude of the thermal stratification parameter. Mallawi et al. [31] derived the series solution of thermally radiative non-Newtonian fluid flow with double stratifications. They have seen that the fluid concentration depresses for enhancing the solutal stratification parameter. Time-dependent MHD nanofluid flow with dual stratifications was performed by Hayat et al. [32]. They proved that surface shear stress enriches for higher values of thermal and solutal stratification parameters.

The aforementioned inspection announces that most of the researchers are involved in revealing the nature of the Darcy-Forchheimer flow with Cattaneo-Christov theory through prescribed wall temperature but not analyzed dual stratifications on a Riga plate. Therefore, our key argument is to fulfill this gap. So, our study elucidates the outcome of the Darcy-Forchheimer flow of a Williamson fluid in the presence of double stratifications, thermal radiation, and chemical reaction on a Riga plate. These types of outcomes will be definitely helpful for a thermal engineer to modeling the thermal systems. Here, the heat and mass transfer phenomena are illustrated by the Cattaneo-Christov dual flux model and the Riga plate is used to control the fluid flow.

2. Mathematical Formulation

Let us consider the 2D Darcy-Forchheimer flow of a Williamson fluid on a Riga plate. Here, the surface temperature and the concentration are denoted by and which are always larger than the free stream temperature and the free stream concentration , respectively, see Figures 1(a) and 1(b). The thermal radiation and first-order chemical reaction are taking into account. Flow situation is manifested with double stratifications. The fluid phase is heat consumption/generation. In addition, the heat and mass transfer phenomenon is inspected through Cattaneo-Christov dual models. The governing equations are modeled as follows: where where is the velocity in and directions, is the kinematic viscosity, is the time constant, is the current density, is the density of the fluid, is the magnetization of the magnet, is the width of the magnet and the electrodes, is the drag coefficient, is the permeability of porous medium, is the fluid temperature, is the relaxation time of heat flux, is the thermal diffusivity, is the specific heat, is the Stefan-Boltzmann constant, is the mean absorption coefficient, is the heat generation/absorption coefficient, is the fluid concentration, is the relaxation time of mass flux, is the mass diffusivity, and is the chemical reaction parameter.

The boundary conditions are

Now, we consider the following dimensionless variables:

By using (7), we can modify equations (2)–(4) as follows:

is the Weissenberg number, is the modified Hartmann number, is the dimensionless parameter, is the local porosity parameter, is the Forchheimer number, is the Prandtl number, is the radiation parameter, is the thermal stratification parameter, is the heat absorption/generation parameter, is the heat relaxation time parameter, is the Schmidt number, is the chemical reaction parameter, is the solutal stratification parameter, is the mass relaxation time parameter, and is the suction/injection parameter.

The corresponding boundary conditions are

The skin friction coefficient, local Nusselt number, and local Sherwood number are expressed as here, the wall shear stress, heat, and mass flux are as follows:

The dimensionless form of the above parameters are expressed as

3. Solutions

3.1. Numerical Solution

In this section, the bvp4c solver has been used for gaining the solution. In order to solve the problem, equations (8)–(10) are commuted into a system of first-order differential equations with the boundary conditions also modified in the same manner. For this, let us take

The system of equations are

With the boundary conditions

The above set of equations are numerically solved by MATLAB built-in function bvp4c.

3.2. HAM Solution

The obtained ODE’s (8)–(10) with conditions (11) are analytically solved by applying the HAM scheme. Because this method is powerful tool for solving nonlinear problems, see Sarwar and Rashidi [33]. Let the initial approximations are chosen as , , and , and linear operators are , , and , where is the differential operator and the linear property is , where are constants.

After implementing the order HAM technique, we found the following: Here, , and are the particular solutions.

These HAM techniques have the parameters (, and ), and these are responsible for the convergence of solutions, see Refs. [3437]. Figure 2 portrays the range value of , and are , , and . We assign for better convergency.

4. Results and Discussion

Here, we revealed the results by graphs and tables which describes the shift in velocity, temperature, concentration, skin friction coefficient, local Nusselt number, and local Sherwood number concerning the disparate values of the parameters, such as Weissenberg number , local porosity parameter , Forchheimer number , modified Hartmann number , thermal radiation parameter , thermal stratification parameter , heat generation/absorption parameter , heat relaxation time parameter , chemical reaction parameter , solute stratification parameter , mass relaxation time parameter , and the suction/injection parameter . The numerically obtained values are compared with the results fetched by the analytical approach by HAM.

Table 1 provides the HAM order and CPU timings. From this table, we concluded that order is sufficient for all profiles. Table 2 delineates the changes of skin friction coefficient, local Nusselt number, and local Sherwood number for the distinct values of , and . We noted that the surface shear stress upsurges when heightening the and values and it declines for enhancing the , , and values. The local Nusselt and Sherwood numbers reduce for raising the , , and , and it rises for increasing the and . Table 3 describes the influence of , , , and over the heat flux. The heat transfer gradient decimates when developing the and values, and it grows when growing the and values. Table 4 helps to figure out the shift of mass flux for the various values of , and . The mass transfer rate escalates for the enriching values of and , and it suppresses for increasing values. Also, we proved that our numerical and analytical results are almost same.

Figures 3(a) and 3(b) establish the impact of (a) and (b) on velocity profile for DFRP and NDFRP. We uncovered that the fluid speed aggravates due to more presence of Hartmann number and quite the opposite behavior is obtained for the parameter. The MBLT is high in NDFRP compared to DFRP for both parameters. The outcomes for disparate values of and on temperature profile are illustrated in Figures 4(a) and 4(b) for DFRP and NDFRP (a) and heat consumption/generation (b). We ascertained that the fluid temperature dwindles because of the high quantity of Hartmann number. However, it is enhanced for raising the radiation parameter. Figures 5(a) and 5(b) explains the changes of fluid temperature for distinct values of (a) and (b) for heat generation/consumption fluid. We revealed that the fluid warmness becomes subsides for hike values of and . In addition, the thermal boundary layer thickness is large in the heat generation case compared to the heat consumption case. The significance of and on concentration profile on destructive chemical reaction and the generative chemical reaction is plotted in Figures 6(a) and 6(b). We ascertained that the fluid concentration lessens for large values of and .

The variations of SFC for distinct combination of and (a) for DF flow (lower plate) and NDF flow (upper plate) and (b) for DF flow of Williamson fluid (upper plate) and viscous fluid (lower plate) are portrayed in Figures 7(a) and 7(b). We found that the plate shear stress depresses for enhancing values of , and it raises for rising the values of Hartmann number. Figures 8(a) and 8(b) disclose the changes of LNN for disparate combination of and (a) for DF flow (lower plate) and NDF flow (upper plate) and (b) for DF flow of Williamson fluid (lower plate) and viscous fluid (upper plate). We noted that heat transfer gradient escalates for increasing the values of and . The alternations of LNN for various combination of (a) and with (upper plate) and (lower plate) and (b) and with (upper plate) and (lower plate) are plotted in Figures 9(a) and 9(b). We noticed that the heat transfer gradient upgrades for the available of , and it weakens for high quantity of and . In addition, the larger heat transfer gradient occurs in heat consumption case compared to heat generation case. Figures 10(a) and 10(b) express the LSN variation with respect to the distinct combination of and (a) for DF flow (upper plate) and NDF flow (lower plate) and (b) for DF flow of Williamson fluid (upper plate) and viscous fluid (lower plate). We proved that the mass transfer gradient slashes due to larger values of and . The deviation of LSN for different combination of (a) and with generative chemical reaction (upper plate) and destructive chemical reaction (lower plate) and (b) and with generative chemical reaction (upper plate) and destructive chemical reaction (lower plate) is shown in Figures 11(a) and 11(b). We concluded that the LSN declines for upgrading the values of , , and .

The decrement percentage of SFC for various values of on Williamson fluid and viscous fluid is plotted in Figures 12(a) and 12(b) and observed that the maximum decreasing percentage of surface shear stress is obtained when varies from 0.0 to 0.4 for both fluids. Figures 13(a) and 13(b) give the increment percentage of LNN for various values of on Williamson fluid and viscous fluid. The maximum increment percentage has occurred when vary from 0.4 to 0.8 for both fluids. The decrement percentage of LSN for various values of on Williamson fluid and viscous fluid is illustrated in Figures 14(a) and 14(b), and we have seen that the maximum decrement percentage occurred when varies from 0.4 to 0.8 for both fluids. Figures 15(a) and 15(b) display the decrement/increment percentage of LNN for various values of and LSN for various values of . We concluded that the heat transfer gradient enhances in heat generation case and maximum increment percentage is obtained when varies from 0 to 0.2 and it suppresses in heat consumption case and maximum decrement percentage is obtained when varies from 0.8 to 1. The mass transfer gradient enhances in destructive chemical reaction case and maximum increment percentage is obtained when varies from 0.8 to 1, and it suppresses in the generative case and maximum decrement percentage is obtained when varies from 0 to 0.2.

5. Conclusion

The current study figures out the Darcy-Forchheimer flow of the Williamson fluid over a Riga plate with Cattaneo-Christov double diffusion and double stratification. The solutions are fetched numerically by a bvp4c solver in MATLAB and analytically by HAM. The outcomes are employed as follows: (i)The fluid velocity drops when the suction/injection parameter enriches(ii)The larger values of the thermal radiation parameter boost up the temperature(iii)The high temperature is noticed in the heat generation case, and the low temperature is produced in the heat absorption case(iv)The generative chemical reaction case upturns the fluid concentration while the destructive chemical reaction case declines(v)The high wall shear stress is produced when enlarging the modified Hartmann number

Abbreviations

CC:Cattaneo-Christov
DC:Destructive chemical reaction
DF:Darcy-Forchheimer
FM:Fourier model
GC:Generative chemical reaction
HA:Heat absorption
HG:Heat generation
LNN:Local Nusselt number
LSN:Local Sherwood number
MBLT:Momentum boundary layer thickness
NDF:Non-Darcy-Forchheimer
RP:Riga plate
SFC:Skin friction coefficient
SP:Stationary plate.

Data Availability

The raw data supporting the conclusions of this article will be made available by the authors without undue reservation.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

All authors contributed equally to this work. And all the authors have read and approved the final version of the manuscript.

Acknowledgments

This research was funded by the Deanship of Scientific Research at Princess Nourah Bint Abdulrahman University through the Fast-track Research Funding Program.