#### Abstract

In this paper, the Darcy–Forchheimer laminar thin film flow with MHD and heat transfer on an unsteady horizontal stretched surface is investigated. The impact of thermal radiation and viscous dissipation is also considered for thin film flow. The analysis of heat source with thermal radiation in a boundary layer flow can play a great role in manufacturing engineering procedures such as the production of electric power, solar energy modernization, and astrophysical flow. By using similarity transformation, the system of PDEs is converted to ODEs. Then HAM is applied for the solution of the problem. Moreover, the velocity and temperature profile for various embedded variables are discussed through graphs while the numerical solution of concerned physical quantities such as skin friction and Nusselt number are discussed through tables. The analysis shows that velocity profile is reducing function of , , , and . The temperature profile is an increasing function and while a reducing function of and . The overall behavior of the proposed system is highlighted analytically.

#### 1. Introduction

Heat source, fluid flow within liquid film due to an unsteady stretching sheet, plays a very crucial role in many utilization engineering processes and technologies; for example, polymer extrusion, continuous casting, manufacturing plastic film, tinning of copper wires, crystal growth, a polymer sheet, artificial fiber, and wire and fiber coating. In addition to these, many production processes involve heat transmission in various forms such as cooling of a machine tool, pasteurization of food, heat treatment of a part in industrial furnaces, or the temperature control for initiating a chemical reaction. Heat transfer is accomplished in most of these applications using heat transfer instruments such as heat swaps, desiccation, condensers, boilers, and heat sinks. Mathematics plays an important role in different areas of research and has effective outcomes [1–4]. The importance of mathematical tools in the investigation of fluid flow cannot be ignored. Mathematical tools help the researchers to visualize the dynamical behavior of a system and to conceptualize the importance of the input factors of a system. Numerous researchers workout to understand the flow of fluid with different assumptions and in different circumstances. The nanofluid flow over a stretching surface has been originally considered by Wang [5]; they also reduce unstable Navier–Stokes equations using similarity transformation in nonlinear ODEs, which are then calculated by multiple shooting methods. The solutions based on the analysis of a liquid thin film are obtained by Wang [6], who also validated the correctness of the homotopy analysis technique. Sufian et al [7] study the impact of heat source and fluid flow over a stretch cylinder as a function of time. Zhang et al. [8] demonstrated mass diffusion and solutal Marangoni impact on heat transmission and fluid flow in a film toward a stretched sheet. Kalyani et al. [9] investigated the effects of Lorentz force, heat generation, and mass transfer of a Casson incompressible fluid flowing on a surface. They also discussed the impact of thermal radiation on the mass and heat transfer property of flow in a Casson liquid film [10]. The unstable mixed convection flow of the stretch sheet is considered by Devi et al. [11]. Jhankal et al. [12] studied heat source, the forced convective flow of an incompressible fluid past plate inserted in a Darcy–Forchheimer porous medium of the two-dimensional boundary layer MHD flow. The effects of thermal radiation in an upper-convective Maxwell thin film, unsteady stretched sheet inserted in a porous medium are conceptualized by Waheed [13]. Ahmad et al. [14] illustrated a solution for the heat source, MHD flow toward an exponentially stretched surface in the absence of radiation impacts and Darcy’s resistance. On the basis of boundary layer hypothesis and similarity transformation, Sparrow and Gregg [15] were the first to examine the Laminar-film condensing difficulty on a vertical plate. Dandapat and Ray explored the flow of a thin boundary layer across a spinning horizontal disc [16, 17]. Nath and Kumari [18] investigated the time-dependent MHD problem of a fluid dispersed by a spinning incompressible disc. Magnetohydrodynamic and heat flow toward a stretch sheet has been studied by Mohd Nasir et al. [19] in which we consider both the stretching and shrinking cases. Zaidi et al. [20] studied the impacts of thermal radiation, heat creation, and induced magnetic field on hydromagnetic free convection flow of pair stress fluid in an isoflux-isothermal vertical channel. Kishore et al. [21] studied the influence of thermal radiation and viscous dissipation on MHD mass, heat transfer flow past an oscillating vertical plate inserted in a porous medium with unknown surface situations. Elbashbeshy et al. [22] studied the influence of thermal radiation and heat transmission on an unstable extended sheet inserted in a porous medium in the existence of a heat sink. Thermal radioactivity and heat transfer of Sisko fluid on an uneven stretched sheet with a constant magnetic field were investigated by Khan et al. [23]. The authors of [24] discussed mass transfer and unstable boundary layer MHD flow of a UCM fluid in the absence of higher-level chemical reaction. The authors of [25] considered heat transfer, MHD flow in a liquid film over an unsteady stretched sheet. Numerical solution of MHD boundary layer flow of a nanofluid over a moving sheet in the absence of thermal radiation has been carried out by Shateyi and Prakash [26]. The MHD laminar boundary layer flow with heat and mass transport of an electrically conducting water-based nanofluid across a nonlinear stretching sheet with viscous dissipation is studied by Mabood and Khan in [27]. It has been noticed that further investigation is needed to conceptualize the impact of heat transfer, thermal radiation in a fluid film of an unsteady stretching surface inserted in a porous medium in the absence of viscous dissipation and magnetic field impact. Therefore, in the current work, we opt to investigate the impact of heat transfer, thermal radiation in a fluid film of an unsteady stretching surface inserted in a porous medium in the absence of viscous dissipation and magnetic field impact.

The aim of the present work is to study Darcy–Forchheimer's laminar thin film flow on an unsteady horizontal stretched surface. The laminar thin-film flow is investigated with the impact of MHD, thermal radiation, and viscous dissipation. The HAM is used to obtain analytic approximations for the solution of nonlinear ODEs. The comparison of HAM and numerical solutions are displayed. The results have been illustrated in graphical or tabular forms followed by a brief discussion and physical interpretations.

#### 2. Mathematical Expressions of the Problem

##### 2.1. Basic Equations and Boundary Constraints

Assume a thin liquid film extending on a stretched surface with a laminar unsteady flow which arises from a narrow slit where x-coordinates are parallel to the slit and y-coordinates are perpendicular to the sheet as exhibited in Figure 1. Along *x*-axis, the stretching sheet is stretched by two equal and opposite forces. Further, the latent heat effect is neglected due to evaporation by assuming the liquid to be nonvolatile. Moreover, due to a relatively thin film, the buoyancy force is also neglected.

The continuous sheet velocity is given by [25], where and are constants having per-time dimension. The term is localized Reynolds number on the sheet velocity , and the stretched rate reduces with time as . The nomenclature of this research work has been presented.

The stretched sheet’s surface temperature is determined as [25]where is the slit temperature and is the temperature with reference, which might be a constant temperature with a reference or a constant temperature difference. In this work, will be interpreted as . Normally, a transverse magnetic field is applied to the extended surface, which is believed to be changeable in nature and chosen as

The temperature and velocity influence of a thin film obeys the following governing time-dependent equations with the boundary conditions [25].where *u* and denote the velocity terms in the directions of *x*-axis and *y*-axis, respectively, is the dynamic viscosity, is the magnetic field, is the density, is the temperature, is the thermal diffusivity, is the kinematic viscosity, is the time, is the electrical conductivity, and is the uniform thickness of the film sheet. In the x- and y-directions, the two terms of equation (3) reflect the net outflow (outflow minus inflow) of mass, for an unsteady flow the sum of these terms must be zero. The net rate at which x-momentum leaves due to fluid motion across its limits is signified by the 2nd and 3rd components on the left side of equation (4) while the net force owing to viscous shear stresses is shown off by the 1st component on the right side. The net rate at which thermal energy escapes due to fluid motion is represented by the 2nd and 3rd terms on the left side of equation (5). The 1st component on the right side of equation (5) shows the net input of thermal energy owing to *y*-direction, while the 2nd term represents the viscous dissipation that remains.

In this case, when the influence of surface tension is insignificant at a distance of from the surface, shear stresses operating in planes parallel to the fluid velocity are defined aswhere is a fluid property and is known as the dynamic viscosity and is the rate of shear stress.

The flow of heat at is defined byIn the above equation, the constant is called thermal conductivity. Equation (8) tells us that if temperature decreases with then will be positive and it will flow in the *y*-direction, if temperature increases with then will be negative and it will flow opposite in the *y*-direction.

##### 2.2. Similarity Transformations

Here, we use the stream function to specify the velocity terms:satisfying continuity equation (3). New dependent and independent variables, , , and are, respectively, introduced as

Substitute the similarity transformation equations (9)–(11) into equations (4) and (5) to obtain a set of nonlinear ODEs in the form of

Following are the corresponding boundary constraints of the differential equations:where the prime denotes derivative in terms of , is the unsteadiness parameter, is the magnetic parameter, is the porosity parameter, is the Eckert number, is the Prandtl number, and is thermal radiation. Despite the fact that the dimensionless film thickness is constant and only relies on , the current film thickness *h* is only depending on time . So the free surface denotes the value of similarity variable by such that

By differentiating with respect to , we have

##### 2.3. Physical Interest of the Problem

In this subsection, we calculated the skin friction coefficient and local Nusselt number as

And

From equations (7) and (8), at the surface *y* = 0, the surface shear stress and surface heat flux are defined as

And the Rosseland estimation is given by [23]where is the Stefan–Boltzman constant and has the numerical value

Inserting equations (19)–(21) into equations (17) and (18) to obtain and in terms of ,where is the local Reynolds number and is the radiation parameter.

#### 3. Solution by HAM

HAM is applied for the solutions of equations (12) and (13) with boundary constraints (12). Mathematica is used for this purpose. The modeled equations are solved by HAM, the basic derivation is identified as

Linear operators , and are specified by

The nonlinear operators and are defined as follows:

The order system for equations (12) and (13) is presented aswhile BCs are as follows

and are utilized to standard for the solution confluence whereas the embedding constraint is ; when and , we get the following:

Through Taylor’s series, and are expanded for :while BCs are

Nowwhile

#### 4. Result and Discussion

The present work analyzes the Darcy–Forchheimer laminar thin film flow with MHD and heat transfer on an unsteady horizontal stretched surface in the absence of viscous dissipation, thermal radiation, and magnetic field of impact. The numerical, as well as analytical solutions, is given. Using similarity transformations, the governing partial differential equations of heat transfer flow are transformed into a nonlinear ordinary differential equation. Geometric representations of the problem are plotted in Figure 1 and the impact of physical parametric quantity on and curves are depicted in Figures 2–5. Figures 2 and 3 show the impact of , , , and on . The velocity profile reduces when () increases, as shown in Figure 2(a). This is due to the reducing nature of the Lorentz forces which reduces the motion of thin-film fluid in the boundary layer and raises its temperature. Figures 2(b) and 3(b) convey that the fluid velocity drop by rising values of () and (). This is because the increasing hydromagnetic boundary layer increases the strength of porosity. Figure 3(a) illustrates the distribution of velocity gradient in the film decreases by rising values of unsteadiness parameter. As the value of () in the system rises, the thickness of the boundary layer declines, preventing the establishment of a transition from laminar to turbulent flow. The effect of embedded parameters , , , and on temperature profiles is exhibited in Figures 4 and 5. Figure 4(a) displays the variation of temperature reduces, as the Prandtl number increases. Being the ratio of molecular diffusivity of momentum to the molecular diffusivity of heat, it demonstrates the relative thickness of the momentum boundary layer to the thermal boundary layer. The rise in Prandtl number () decreases the temperature of the fluid. Increasing the Prandtl number, the thermal conductivity of the fluid decreases and the rate of heat transfer through stretch surface reduces; consequently, the temperature and thermal boundary layer reduce. Figure 4(b) demonstrates the temperature gradient for different values of (). As the values of () upsurge, the temperature of thin film declines. Physically, when the radiation parameter increases, the Rosseland radiative absorption coefficient declines as defined by the expression for . Therefore, the heat flux diminishes. That is why the thin film temperature reduces. Figure 5(a) discusses the impact of unsteadiness parameter () on temperature profiles. As the values of () rise, the thermal boundary layer grows as a result, and the temperature field becomes increment for the nanofluid. Figure 5(b) demonstrates that increases in the value of () produce more temperature in the thin film. Eckert number is the ratio among kinetic energy and enthalpy, which is detected in the dissipation function of the energy equation, also known as the dissipation parameter. = 0 means there is no viscous dissipation. Figure 5(b) displays that the behavior of () is directly proportional to the temperature, as () upsurge viscous dissipation rises, i.e., kinetic energy transformed into heat energy which in turn raises the temperature.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

The physical significance of embedded parameters are displayed in Table 1 and the numerical values of skin friction and local Nusselt number are displayed in Tables 2 and 3. Table 2 demonstrates that the skin friction coefficient increases by rising various parameters, namely, , , , and . As these parameters increase, the unsteady parameter, the skin friction enhances due to the exertion of a drag force in the fluid thin film. The Lorentz force is a rising magnetic field that increases the resistive force on fluid movement. That is why more values of () enhance . Table 3 lists the numerical values of for a few physical variables, that is, , , , and . The Eckert number is normally made up of the nanofluid thermal conductivity components, which is used to raise the temperature field, which has the reverse effect on cooling operations, and hence the Nusselt number reduced for high values of (). With an increase in (), the boundary layer temperature of fluid raises, increasing the heat source of the thin film fluid; as a result, boosts. In this table, we may observe that the Nusselt number upsurges as () grows, but drops as () grows. Table 4 displays excellent agreement between the HAM and ND solve solution through computer-based package Mathematica 11.1.0.

#### 5. Concluding Remarks

The boundary layer two-dimensional Darcy–Forchheimer laminar thin film flow over a horizontal stretched sheet with the impact of thermal radiation and viscous dissipation are considered. The cause of heat source and thermal radiation over an unsteady stretched surface inserted in a porous medium in the absence of a magnetic field has been solved by HAM. The important key findings of this research work are given below:(i)Velocity profile reduces for raise in values of and (ii)When the values of and increase, the velocity profile reduces(iii)Temperature field increases with an increase in the values of and (iv)With increasing values of and , the temperature field reduces(v)Skin friction coefficient rises as upsurges the values of physical parameters that is , , , and (vi)By increasing the values of and , the local Nusselt number upsurge while it declines for more values of and

#### Nomenclature

Symbols: | Description |

: | Horizontal velocity component (ms^{−1}) |

: | Vertical velocity component (ms^{−1}) |

: | Horizontal coordinate (m) |

: | Vertical coordinate (m) |

: | Stretching rate (s^{−1}) |

: | Sheet velocity (ms^{−1}) |

: | Temperature (K) |

: | Film thickness (m) |

: | Dimensionless stream function, equation (9) |

: | Skin friction coefficient, shown in Table 1 |

: | Time (s) |

: | Unsteadiness parameter, |

: | Specific heat (JKg^{−1}K^{−1}) |

: | Heat flux, (Js^{−1}m^{−1}) |

: | Magnetic parameter, |

: | Porosity parameter, |

: | Local Reynolds number, shown in Table 1 |

: | Prandtl number, shown in Table 1 |

: | Local Nusselt number, shown in Table 1 |

: | Eckert number, shown in Table 1 |

Greek Symbols: | Descriptions: |

: | Dimensionless film thickness, |

: | Dimensionless temperature, equation (9) |

: | Constant (s^{−1}) |

: | Thermal conductivity (m^{2}s^{−1}) |

: | Kinematic viscosity (m^{2}s^{−1}) |

: | Dynamic viscosity (Kgm^{−1}s^{−1}) |

: | Similarity variable, equation (9) |

: | Shear stress, (Kgm^{−1}s^{−2}) |

: | Stream function (m^{2}s^{−1}) |

: | Density (Kgm^{−3}) |

Subscripts: | Descriptions: |

o: | Origin |

Ref: | Reference value |

: | Sheet |

: | Localized value |

Superscripts: | Descriptions: |

: | 1nd derivative |

: | 2nd derivative |

: | 3rd derivative |

#### Data Availability

No data are available regarding this manuscript.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.