A Numerical Study on Newtonian Heating Effect on Heat Absorbing MHD Casson Flow of Dissipative Fluid past an Oscillating Vertical Porous Plate

The numerical study of Newtonian heating eﬀect on unsteady free convection MHD Casson ﬂow of radiating and chemically reacting ﬂuid past an oscillating vertical porous plate embedded in a porous medium was conducted by considering the eﬀects of heat sink and viscous dissipation. The ﬂuid motion is persuaded due to the periodic oscillations of the plate along its length. This phenomenon is represented as nonlinear PDEs with initial and boundary conditions. By introducing suitable nondimensional variable and parameters, the leading equation with initial and boundary conditions is converted into dimensionless form, which are then solved numerically using a ﬁnite diﬀerence method. The eﬀects of several relevant parameters on the velocity, temperature, and concentration are displayed graphically, whilst the eﬀects of these parameters on the skin friction and Nusselt and Sherwood numbers are exhibited in tabular format and then discussed in detail. The ﬁnal outcomes divulge that the radiation parameter and Eckert number have an increasing eﬀect on the velocity and temperature, whilst reverse tendency is detected with increasing Prandtl number and heat absorption parameter. Newtonian heating parameter and thermal and mass buoyancy forces boost ﬂuid velocity, whilst Schmidt number and chemical reaction have the opposite impact. It is noteworthy to point out in this study that the velocity boundary layer thickness for the Newtonian ﬂuid is larger than the Casson ﬂuid.


Introduction
e study of non-Newtonian fluid flow models has fascinated investigators owing to its significant applications in chemical, pharmaceutical, and cosmic and food industries such as in the production of several chemicals, foams, oil, gas, paint, syrup, juice, chocolate, and many other food products. However, the study of non-Newtonian fluids is not as straightforward as Newtonian fluids, since there does not exist a single constructive equation for the non-Newtonian fluids that can be used to explicate them. In order to analyze their characteristics, a number of researchers introduced different constitutive equations or models. e foremost investigations [1][2][3][4][5][6][7][8][9] considered, respectively, are power law, second-grade, Jeffery, Maxwell, viscoplastic, Bingham plastic, Brinkman type, Oldroyd-B, and Walters-B fluid models. e trendiest Casson model was first introduced for the prediction of the flow narration of pigment oil suspensions of the printing ink type by Casson [10]. Subsequently, Mustafa et al. [11] studied unsteady boundary layer flow of a Casson fluid due to an impulsively started moving flat plate. Rao et al. [12] analyzed heat transfer in a Casson rheological fluid from a semi-infinite vertical plate with partial slip conditions. Khalid et al. [13] studied unsteady MHD free convection flow of Casson fluid past over an oscillating vertical plate embedded in a porous medium. Kataria and Patel [14] examined radiation and chemical reaction effects on MHD Casson fluid flow past an oscillating vertical plate embedded in a porous medium. Kataria and Patel [15] studied heat and mass transfer in MHD Casson fluid flow past over an oscillating vertical plate embedded in a porous medium with ramped wall temyperature. Mahanthesh et al. [16] studied the two-phase flow of dusty Casson fluid Cattaneo-Christov heat flux and heat source past a cone, wedge, and plate. Shamshuddin et al. [17] analyzed chemical reaction effects on unsteady free convection radiated Casson fluid flow over an inclined porous plate numerically. Shamshuddin et al. [18] presented chemical reaction effects on thermosolutal micropolar fluid past a permeable stretching porous sheet. Shamshuddin et al. [19] examined radiation, heat source/sink effects on micropolar fluid flow due to a permeable stretching sheet with partial slip, and surface heat flux boundary conditions. Rajput et al. [20] studied numerically radiated Casson thermosolutal convective fluid flow over a vertical plate propagated by Arrhenius kinetics with heat source/sink.
In all these investigations, the conditions at the wall surface are commonly ramped, constant or variable wall temperature. However, there were several problems of physical interest where the heat is transported to the fluid via a bounding surface with a finite heat capacity. In this case, the above thermal boundary conditions fail to work; consequently, the Newtonian heating conditions are highly required.
e applications of Newtonian heating are significant in heat exchangers, petroleum industry, turbines, and also in convective flows set up when the bounding surface absorbs heat by solar radiation. Merkin [21] first introduced the Newtonian heating condition in natural convection boundary layer flow over a vertical surface. Lesnic et al. [22,23] studied free convection boundary layer flows along vertical and horizontal surfaces in a porous medium generated by Newtonian heating. Chaudhary and Jain [24] studied unsteady free convection boundary layer flow past an impulsively started vertical plate with Newtonian heating. Salleh et al. [25] studied forced convection boundary layer flow near the forward stagnation point of an infinite plane wall generated by Newtonian heating. Salleh et al. [26] discussed the impact of Newtonian heating on boundary layer flow and heat transfer over a stretching sheet. Hussanan et al. [27] performed an exact analysis of heat and mass transfer flow past a vertical plate in the presence of Newtonian heating. Hussanan et al. [28,29] evaluated the influence of Newtonian heating on unsteady boundary layer flow by considering non-Newtonian and Newtonian fluids through the porous medium. Unsteady hydromagnetic natural convection flow past an impulsively moving vertical plate with Newtonian heating in a rotating framework was reported by Seth et al. [30]. Das et al. [31] analyzed the effects of Newtonian heating on MHD heat and mass transfer unsteady Casson fluid flow past a flat plate. Hussanan et al. [32] scrutinized the impact of Newtonian heating on heat transfer in the magneto-hydrodynamic flow of a Casson fluid through the porous medium. Recently, Prabhakar Reddy [33] evaluated the effect of Newtonian heating on a radiating hydromagnetic flow past an impulsively moving infinite vertical plate in the presence of Hall current.
In most of the studies mentioned above, the effects of viscous dissipation on the flow are ignored. On the other hand, in some cases, the fluid viscosity may experience a significant change with variation in the temperature and also, for the extreme size of flow field or at low temperature, the impact of viscous dissipation is essential. is type of flow takes place in many industrial applications such as the aerodynamic extrusion of plastic sheets, material handling conveyors, cooling or drying of papers and textiles, and glass fibre production. Due to these essential applications, the study of viscous dissipation on the flow becomes a special curiosity among many investigators. Mansour et al. [34] and Suneetha et al. [35] considered in their study the effect of viscous dissipation in MHD natural convection flow under different situations. Kumar [36] and Singh and Singh [37] scrutinized the impact of viscous dissipation in hydromagnetic flow through the stretching surface. Kishore et al. [38] analyzed the effects of thermal radiation and viscous dissipation on MHD heat and mass diffusion flow past an oscillating plate embedded in a porous medium with variable surface conditions. Pal et al. [39] studied convectiveradiation effects on stagnation point flow of nanofluid over a stretching surface with viscous dissipation. Pal and Biswas [40] examined the viscous dissipation effects on mixed convective MHD oscillatory flow of Casson fluid in the presence of thermal radiation, chemical reaction, and Soret effect. Prabhakar Reddy [41] discussed viscous dissipation and mass transfer effects on unsteady MHD free convective flow of an incompressible past an infinite vertical porous plate. Zigta [42] analyzed the influence of viscous dissipation on hydromagnetic flow in between a pair of infinite vertical Couette channel walls in the presence of thermal radiation and chemical reaction. Shamshuddin et al. [43] studied dissipative radiated micropolar convective fluid from a rotating vertical plate with oscillatory plate velocity adjacent to a permeable medium numerically. Shamshuddin et al. [44] studied numerically rotating unsteady multiphysicochemical magneto-micropolar flow in a porous medium from a vertical plate in the presence of Hall current. Shamshuddin and Tirupathi [45] carried out a numerical investigation for dissipative micropolar fluid flow past an inclined porous plate with heat source/sink. Ferdows et al. [46] studied micropolar fluid flow past a horizontal stretching sheet through a porous medium in the presence of Soret-Dufour effect and viscous dissipation. Prabhakar Reddy and Muthucumaraswamy [47] analyzed numerically the impact of viscous dissipation on MHD flow past an oscillating vertical porous plate with variable surface conditions in the presence of thermal radiation and chemical reaction. e above comprehensive literature analysis motivated us to determine the impact of Newtonian heating effect on unsteady MHD Casson fluid past an oscillating vertical porous plate in the existence of thermal radiation, chemical reaction, heat sink, and viscous dissipation effects as this problem was not being discussed. is kind of exploration has several medical and industrial applications such as in the study of blood flow in the cardiovascular system, glass manufacturing, paper production, and purification of crude oil. e leading coupled nonlinear PDEs of the model have been solved numerically by adopting the finite difference method. is numerical approach coupled with the model problem, which in essence extends the previously published work of Das et al. [31] by including heat sink, surface porosity, and viscous dissipation effect, which comprise the novelty of the present study. Further, to validate the present 2 International Journal of Mathematics and Mathematical Sciences numerical procedure, a comparison was executed in the limiting conditions between the present results and earlier published works in the literature. It has been found significant conformity between the results, which provides the accuracy of the present results.

Mathematical Model
We consider an unsteady heat-absorbing Casson MHD free convection heat and mass transfer flow of viscous, incompressible electrically conducting fluid past a vertical porous plate oscillating in its own plane fixed at y ′ � 0, where y ′ is the coordinate axis normal to the plate as shown in Figure 1. e uniform magnetic field of strength B 0 applied in a transverse direction to the flow. e induced magnetic field formed by the motion of the fluid is neglected in comparison to the applied one as the magnetic Reynolds number is small enough. Initially, at time t ′ � 0, both the fluid and plate are at uniform temperature T ∞ ′ and the concentration near to the plate C ∞ ′ . At time t ′ > 0, the plate starts to oscillate in the vertical direction in its own plane with velocity u ′ � U 0 cos(ω ′ t ′ ) against the gravitational field, where U 0 is the amplitude and ω ′ is the frequency of the plate oscillations. At the same time, between the local surface temperature T ′ and heat transfer from the plate, there is a direct proportionality relation and the concentration level near to the plate surface raised from C ∞ ′ to C w ′ . ere are first-order chemical reactions existing between the diffusing species and the fluid. e viscous dissipation term is taken into account in the energy equation. e rheological equation of the state for the isotropic and incompressible flow of a Casson fluid following Khalid et al. [13] and Hussanan et al. [28] is given as follows: or where τ is the shear stress, τ 0 is the Casson yield stress, μ is the dynamic viscosity, α * is the shear rate, π � e ij e ij and e ij is the (i, j) th component of the deformation rate, π is the product of the component of deformation rate with itself, π c is the critical value of this product based on the non-Newtonian fluid model, μ B is the plastic dynamic viscosity of the non-Newtonian fluid, and p y is the yield stress of the fluid. By considering these mentioned assumptions, the governing system of PDEs of the model following Das et al. [31] and Hussanan et al. [32] is as follows: where u ′ is the velocity of the fluid, g is the acceleration due gravity, β T is the volumetric coefficient of thermal expansion, β C is the volumetric coefficient of concentration expansion, T ′ is the temperature of the fluid, C ′ is the concentration of the fluid, k is the thermal conductivity, υ is the kinematic viscosity, μ is the fluid viscosity, ρ is the fluid density, σ is the electrical conductivity, c p is the specific heat at constant pressure, q r is the radiation heat flux, Q 0 is the heat absorption constant, k r ′ is the chemical reaction constant, D m chemical molecular diffusivity, φ is the porosity, k 1 is the permeability of the medium, and t ′ time. e associated initial and boundary conditions of the illustrated model are as follows: where h s is the heat transfer coefficient. e radiative heat flux q r under the Rosseland approximation is given by the following equation: where σ * and k * denote, respectively, the Stefan-Boltzmann constant and absorption coefficient. According to Das et al. [31], after simplifying T ′ 4 , by employing the Taylor series, we obtain the following equation: Substituting equations (5) and (6) into equation (2), the energy equation reduces to Introducing the following nondimensional variables and parameters, we obtain the following equation: Substituting equations (1) and (3) into (6), we obtain the dimensionless governing system of PDEs of the model as follows: 4 International Journal of Mathematics and Mathematical Sciences e analogous initial and boundary conditions become

Numerical Procedure
e modeled systems of coupled nonlinear PDEs in equations (8)-(10) are solved with connected initial and boundary conditions of equation (11) In general, finding the exact solutions of these equations is not possible. erefore, DuFort-Frankel's explicit finite difference technique is utilized to solve the governing system of equations (8)-(10) with related initial and boundary conditions of equation (11), which is quite easy, accurate, efficient, and unconditionally stable. e comprehensive elucidation of this numerical procedure was given in Jain et al. [48]. On applying the DuFort-Frankel's procedure to equations (8)-(10), it results in International Journal of Mathematics and Mathematical Sciences 5 e initial and boundary conditions of equation (11) can be expressed as follows: θ(0, j)), Here, the suffix i corresponds to space ξ and the suffix j corresponds to time τ, Δξ is the mesh size along ξ− direction, and Δτ is the mesh size along τ− direction. e physical domain of the problem is considered as finite restricted as the rectangle of finite dimensions for computational purpose. e finite difference equations (12) and (13) under the initial and boundary conditions (16) at every internal node on a particular j level constitute a tri-diagonal system of equations, which are solved by using omas algorithm. e free-stream boundary condition ξ ⟶ ∞ has been fixed by ξ max � 4, the value greater than 4 does not result in any significant change in the numerical values of u(ξ, τ), θ(ξ, τ), and C(ξ, τ) under various key parameters. By fixing the mesh size Δτ � 0.01, the developed code for the present numerical scheme was run with various mesh sizes Δξ � 0.001, 0.002 and 0.009, where no substantial variation was detected in the numerical values of velocity u, temperature θ, and concentration C. us, the mesh size Δξ � 0.001 has been fixed for the complete computational process. Some essential physical quantities connected to the heat and mass transfer characteristics of the transport, such as the skin friction and Nusselt and Sherwood numbers, are defined in dimensionless form by the following equation:

Validation of the Present Numerical Procedure.
In order to exhibit the validation and accuracy of the present numerical procedure, a comparison of the skin friction calculated in the present work and earlier published works of Das et al. [31] and Hussanan et al. [32] obtained analytically is provided in Table 1 when α ⟶ ∞, G m � 0, M � 0, K ⟶ ∞, N r � 0, Q � 0, E c � 0, and C � 0. It is perceived from Table 1 that there is an outstanding agreement among the results of Das et al. [31] and Hussanan et al. [32] obtained analytically and the present results obtained numerically by the developed code in this work. is entrusts the accuracy of the present numerical procedure.

Graphical Results and Discussion
To investigate the physical significance of the modeled problem, the influence of various key parameters involved in the problem such as magnetic parameter M, heat absorption parameter Q, radiation parameter N r , Eckert number E c , Newtonian heating parameter c, Prandtl number P r , Schmidt number S c , chemical reaction rate k r , porosity parameter K, phase angle ωτ, Casson parameter α, thermal Grashof number G r , and mass Grashof number G m on the fluid velocity, temperature, and concentration are evaluated numerically and executed graphically, whilst the skin friction and Nusselt and Sherwood numbers are put in tabular format. e numerical calculations have been done by adopting the default values of several involved parameters in the problem as P r � 0.71, N r � 0.71, Q � 1, E c � 0.01, c � 0.4, S c � 0.6, k r � 0.5, G r � 5, G m � 4, M � 2, K � 0.5, α � 1, and ωτ � π/3 until otherwise specified particularly. e variation of velocity and temperature distributions against ξ under the influence of Prandtl number P r is shown in Figure 2 for mixture of noble gas (P r � 0.5), air (P r � 0.71), gas (P r � 1), and water (P r � 7). It is observed a decreasing effect on both velocity and temperature in the boundary layer as the Prandtl number increased. Physically, an increase in P r leads to an increase of fluid viscosity, so that fluid becomes thick, and consequently, fluid velocity and temperature decrease. e effect of heat absorption parameter Q on fluid velocity and temperature against ξ is depicted in Figure 3. It is clearly seen that increasing values of heat absorption Q give lessening fluid temperature and velocity. Physically, when heat absorption is present in the energy equation, it absorbs energy and as a result declines the fluid temperature and consequently reduces fluid motion. Figure 4 exhibits the effect of radiation N r on fluid velocity and temperature against ξ. It is seen that growing values of N r improve the thickness of thermal and momentum boundary layers, which gives an increasing effect on both fluid motion and temperature. Figure 5 illustrates the variation of velocity and temperature distributions against ξ under the effect of Eckert number E c . We found a rise in both velocity and temperature of the fluid with enhancement in E c . Physically, increasing Eckert number E c converts the kinetic energy into the internal energy, which causes to enhance both fluid velocity and temperature. e effect of Newtonian heating on fluid velocity and temperature against ξ is shown in Figure 6. A significant effect on both fluid velocity and temperature within the boundary layer is noticed. Increasing values of Newtonian heating 6 International Journal of Mathematics and Mathematical Sciences   Figure 7 depicts the effect of Schmidt number S c on the velocity and concentration distributions against ξ. By the definition, Schmidt number S c is the relative influence of momentum diffusion on species diffusion. e momentum diffusion is quicker than species when S c > 1, and opposite is true when S c < 1. When S c � 1, both momentum and species boundary layers are of the same order of magnitude. As expected, when values of S c are increased, both fluid concentration and velocity fall down due to increasing momentum being diffused at a slighter rate than species. Figure 8 is plotted to show the variation of velocity and concentration distributions against ξ for various values of chemical reaction rate k r . It is observed that a decreasing effect on both fluid velocity and concentration with enhancing values of chemical reaction rate. is is because increasing values of k r reduce the chemical molecular diffusivity, and as a result, fluid concentration lessens therefore fluid velocity decreases. Figure 9 displays the effect of thermal Grashof number G r and mass Grashof number G m on the velocity distribution against ξ. It is clearly seen a remarkable increasing effect on the fluid velocity near to the plate surface and then starts to decline and receive the stationary position in the free stream with an increase in either G r or G m . By the definition, G r is the ratio of thermal buoyancy force to the viscous hydrodynamic force and G m is the ratio of species buoyancy force to the viscous hydrodynamic force acting on the fluid transport. Increasing in either G r or G m causes to decrease in viscous hydrodynamic forces and consequently, the momentum of the fluid is higher. e effects of magnetic parameter M and permeability parameter K on the fluid velocity against ξ are illustrated in Figure 10. It is seen from Figure 10(a) that a decreasing effect on the fluid velocity with raising values of M and the reverse effect is noted on the fluid velocity with increasing values of K in Figure 10(b). When the applied magnetic field is active against the fluid motion, the strength of the Lorentz force is maximum and consequently, fluid motion decreases. An increase in the values of K furnishes reduction in the resistance of the porous medium as a result accelerates the fluid motion. Figure 11 elucidates the effects of Casson parameter α and phase angle ωτ on the velocity against ξ. It is observed from Figure 11(a) that an increase in Casson fluid parameter α tend to decrease in the fluid motion. Physically, an increase in α causes to increase of plasticity of the fluid and consequently, fluid motion decreases. Figure 11(b) exposed that the fluid behavior is We conclude from these results that the velocity boundary layer thickness for the Newtonian fluid is larger than the Casson fluid.
e computed numerical results of the skin friction and Nusselt and Sherwood numbers are listed in Tables 2-4 respectively. It is clearly seen from Table 2 that there is a rise in the skin friction with an enhancement in values of radiation parameter, Eckert number, Newtonian heating parameter, porosity parameters, and thermal and mass Grashof numbers, whilst reverse tendency is perceived when the values of Prandtl number, heat absorption parameter, Schmidt number, chemical reaction rate, magnetic parameter, phase angle,  Table 3 lists that the Nusselt number increases with growing values of radiation parameter, Eckert number, Newtonian heating parameter, and time, whilst opposite behavior is noticed when the Prandtl number and heat absorption parameter are increased.
It is professed from Table 4 that an increase in Schmidt number and chemical reaction rate tends to ascend Sherwood number, whilst overturn is noticed with time progression.
Bold values represent the variation of the physical parameter. International Journal of Mathematics and Mathematical Sciences 13 the present work and previously published work of Das et al. [31] computed analytically, which provides the exactness of our present results.

Conclusions
In this article, the governing system of PDEs of the model is analyzed numerically to determine the impact of Newtonian heating on unsteady heat-absorbing free convection MHD Casson dissipative fluid flow of radiating and reacting past an oscillating vertical porous plate. Lastly, this study can be further extended to include the Soret-Dufour effect in the presence of different boundary conditions, which is momentous when the density variation exists in the transport. Also, such analysis is essential in the context of geosciences, oceanography, and chemical engineering. e major concluding observations of this analysis can be summarized as follows: 1. e velocity and temperature increase with rising values of N r , E c , and c, while an increase in P r and Q have opposite effect. 2. An increase in S c causes to decline velocity and concentration, while reverse effect is noticed when k r is increased. 3. Increasing K, G r , and G m escalates the velocity, while increasing M, α, and ωτ have opposite effect. 4. Skin friction decreases with growing values of P r , Q, S c , k r , M, α, and ωτ, while it raises with increasing values of N r , E c , c, K, G r , and G m . 5. Nusselt number increases with increasing N r , E c , c, and τ, while it reduces with increasing P r and Q. 6. Increasing S c and k r enhances Sherwood number, while it lessens with progression of τ.

Data Availability
e data used to support the findings of this research are available within the article.
is kind of investigation has several medical and industrial applications such as in the study of blood flow in the cardiovascular system, glass manufacturing, paper production, and purification of crude oil.

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