Next Article in Journal
The Theoretical Relationship between the CCR Model and the Two-Stage DEA Model with an Application in the Efficiency Analysis of the Financial Industry
Next Article in Special Issue
Two-Phase Biofluid Flow Model for Magnetic Drug Targeting
Previous Article in Journal
Feature-Level Fusion of Finger Vein and Fingerprint Based on a Single Finger Image: The Use of Incompletely Closed Near-Infrared Equipment
Previous Article in Special Issue
MHD Bioconvection Flow and Heat Transfer of Nanofluid through an Exponentially Stretchable Sheet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving the Nonlinear Boundary Layer Flow Equations with Pressure Gradient and Radiation

by
Michalis A. Xenos
1,†,
Eugenia N. Petropoulou
2,*,†,
Anastasios Siokis
1,‡ and
U. S. Mahabaleshwar
3
1
Department of Mathematics, University of Ioannina, 45110 Ioannina, Greece
2
Department of Civil Engineering, Geotechnical Engineering Laboratory, University of Patras, 26504 Patras, Greece
3
Department of Mathematics, Davangere University, Shivagangotri, Davangere 577 007, India
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, 38124 Braunschweig, Germany.
Symmetry 2020, 12(5), 710; https://doi.org/10.3390/sym12050710
Submission received: 20 February 2020 / Revised: 3 April 2020 / Accepted: 14 April 2020 / Published: 2 May 2020
(This article belongs to the Special Issue Fluid Mechanics Physical Problems and Symmetry)

Abstract

:
The physical problem under consideration is the boundary layer problem of an incompressible, laminar flow, taking place over a flat plate in the presence of a pressure gradient and radiation. For the mathematical formulation of the problem, the partial differential equations of continuity, energy, and momentum are taken into consideration with the boundary layer simplifications. Using the dimensionless Falkner–Skan transformation, a nonlinear, nonhomogeneous, coupled system of partial differential equations (PDEs) is obtained, which is solved via the homotopy analysis method. The obtained analytical solution describes radiation and pressure gradient effects on the boundary layer flow. These analytical results reveal that the adverse or favorable pressure gradient influences the dimensionless velocity and the dimensionless temperature of the boundary layer. An adverse pressure gradient causes significant changes on the dimensionless wall shear parameter and the dimensionless wall heat-transfer parameter. Thermal radiation influences the thermal boundary layer. The analytical results are in very good agreement with the corresponding numerical ones obtained using a modification of the Keller’s-box method.

Graphical Abstract

1. Introduction

The mathematical description of the laminar, incompressible boundary layer flow remains a challenging problem with several aspects still unexplored. The first to describe the steady two-dimensional boundary layer flow over a semi-infinite flat plate was Blasius [1]. The adverse pressure gradient on the laminar boundary layer flow was numerically examined by Howarth [2]. Other studies on the effects of adverse pressure gradients, heat transfer, or thermal radiation on the fluid flow include, although not limited to, [3,4,5,6,7,8]. The effect of radiation in the presence of a magnetic field or with temperature dependent viscosity was studied in [9,10].
At high temperatures, thermal radiation has a significant effect on the flow field, which has important applications in many areas of engineering. This effect has been reported in several works where various flows in the presence of radiation have been studied, such as the free convective laminar flow in [11,12] and [13], the flow of an optically thin gray fluid in [14,15], and [10] or the flow past a moving vertical plate or a vertical oscillating plate in [16,17].
In [18], a computational model was presented for convective and radiative heat transfer of very high temperature gas cooled and gaseous core reactors, with results that are in very good agreement with experimentally based correlations. In [19,20], the interaction of flow with radiation was studied in hypersonic boundary layer flows. In [21,22], a comprehensive numerical analysis of heat transfer in a rectangular enclosure was performed, and several effects on the fluid flow and heat transfer have been extensively studied, such as the effects of the Rayleigh number, thermal conductivity ratio, as well as internal surface emissivity. These analyses revealed a significant cooling effect of the internal volume together with an increase in the thermal conductivity ratio [22]. In [23], it was shown that the thermal behavior of a compressible turbulent flow was influenced by radiation and that the radiative heat flux plays a role in the recirculating zone. More recent studies show the impact of thermal radiation and the viscous dissipation effect on the dynamic three-dimensional rotating flow of carbon nanotubes [24]. In [25], the slip effects on magnetohydrodynamic (MHD) axisymmetric buoyant nanofluid flow over a stretching sheet in the presence of radiation and a chemical reaction were studied. In [26], the mixed convection stagnation-point flow of a nanofluid past a permeable stretching or shrinking sheet was numerically examined with suction, radiation, and a heat source/sink.
In this study, the standard formulation of the boundary layer equations is used and analytically solved until the point of flow separation, for the steady state, two-dimensional, laminar, boundary layer flow of an incompressible fluid over a flat plate under radiation, and pressure gradient effects. The nondimensional form of the corresponding equations with pressure gradient and thermal radiation are obtained in Section 2, utilizing the Falkner–Skan transformation [8]. The resulting nondimensional equations form a nonlinear, coupled system of partial differential equations (PDEs), accompanied by adequate boundary conditions. Such kinds of systems have been studied in the literature using mostly numerical methods since finding an analytical solution is very difficult, if at all possible. However, instead of using a numerical method in order to solve this system of PDEs, one may choose an approximate technique instead. The main advantage when using an approximate technique instead of a numerical is that a closed form approximate formula for the solution is obtained. This may lead to useful formulae, even if they are approximations, for quantities of physical interest. In addition, the implementation of approximate techniques is highly facilitated nowadays by the use of symbolic computation software, such as Mathematica or Maple.
Approximate methods have been successfully applied in the past in order to study fluid mechanics problems. The most well known approximate methods are the perturbation techniques [27,28]. For recent applications of perturbation techniques to fluid mechanics problems, one may refer to [29,30,31,32,33,34]. Other approximate techniques have also been applied in the study of various physical problems, including boundary layer flows. The references of the present paper are confined to few indicative papers of the last decade. The variational iteration method was utilized in [35] for solving the Blasius equation and in [36] for studying the magnetohydrodynamic flow over a nonlinear stretching sheet. The Adomian decomposition method was used not only for the solution of the Blasius equation in [37] but also for the study of micropolar flow in a porous channel [38]. Also, the differential transform method was used in [39] and [40] for the study of specific wedge flow and magnetohydrodynamic flow problems. It is worth mentioning that the differential transform method, combined with the method of steps, has been successfully applied for solving partial differential equations [41], as well as delayed differential equations and systems [42,43].
Maybe a breakthrough for the approximate methods was Liao’s homotopy analysis method (HAM), which was introduced in [44]. Since then, it has been extensively used for a variety of physical problems, such as beam problems, oscillation problems, wave propagation problems, and of course, fluid mechanics problems. A thorough description of HAM can be found in the books [45,46], where it is also shown how other approximate methods can be unified by HAM. Probably the first papers where HAM was used for the study of boundary layer problems were [47,48,49], where the Blasius and the Falkner Skan equations were solved. Ever since, HAM has been successfully used, literally in hundreds of papers, for the study of boundary layer flows. Indicatively, it is mentioned that recently HAM has been used for the study of flows of nanofluids [24,50,51,52,53,54,55] and flows in porous media [56,57,58,59].
In the vast majority of the papers regarding boundary layer flows where HAM was used, the nonlinear equations under consideration were homogeneous with boundary conditions at . There are, however, also papers of this kind where the solved nonlinear equations are nonhomogeneous, but in most cases, the nonhomogeneous term is a constant. Even less are the papers of this kind, where the boundary conditions are not placed at but at some finite point.
In this paper, the steady state, two-dimensional, laminar, boundary layer flow of an incompressible fluid over a flat plate under radiation and pressure gradient effects, is solved for the first time using HAM, an analytical technique. This problem is mathematically described by a nonlinear, nonhomogeneous, and coupled system of PDEs. The nonhomogeneous terms in the momentum and energy equations are nonconstants, and in addition, the boundary conditions are placed not at but at the edge of the boundary layer. The appearance of such nonhomogeneous terms as well the specific type of boundary conditions, increase the difficulty of applying HAM. To the best of the authors’ knowledge, only in [60] were similar equations and boundary conditions solved via HAM, where the two-dimensional, steady, boundary layer flow of a heat-conducting perfect gas was studied. However, no radiation effects were taken into consideration in the aforementioned paper.
The rest of the paper is organized as follows: In Section 2, the physical problem is given and the mathematical equations that describe it are derived. In Section 3, a short description of HAM is given together with its specific application for the physical problem studied in this paper. Finally, in Section 4, a thorough discussion of the obtained results is given with respect not only to the effect of a pressure gradient and radiation but also to numerical results concerning similar but not identical physical problems [3,5,6,19,20,61,62]. The presented results are novel, especially concerning the radiation effect on the boundary layer flow under the effect of a pressure gradient. More precisely, the obtained analytical results, reveal that the adverse and favorable pressure gradients substantially influence the dimensionless velocity of the boundary layer. Pressure gradient also plays a significant role in the changes of the dimensionless temperature of the boundary layer. Furthermore, the adverse pressure gradient causes significant changes on the dimensionless wall shear parameter and the dimensionless wall heat-transfer parameter, leading to boundary layer separation. Thermal radiation significantly influences the thermal boundary layer.

2. Mathematical Formulation

Consider the steady state, two-dimensional laminar flow of an incompressible fluid in the boundary layer that forms over a flat plate parallel to the free stream velocity, U . In a Cartesian coordinate system, O x y , the surface is located at y = 0 , 0 x L , as shown in Figure 1.
The governing equations, after the boundary layer simplifications, form a parabolic system that constitutes of the continuity equation:
u x + v y = 0 ,
the x-momentum equation:
u u x + v u y = 1 ρ d p d x + ν 2 u y 2
and the energy equation:
ρ C p u T x + v T y = u d p d x + k 2 T y 2 + μ u y 2 q r y ,
with boundary conditions defined as follows:
y = 0 : u = v = 0 , T = T w , y δ : u = u e , T = T e ,
where u, v are the velocity components, T the temperature, p the pressure, ρ the density, ν the kinematic viscosity, μ the viscosity, C p the specific temperature of the fluid under constant pressure, k the coefficient of thermal conductivity, q r y the local radiant absorption, T w the temperature of the flat plate, u e , T e the velocity and the temperature at the edge of the boundary layer, respectively, and δ the boundary layer thickness.
Assuming that the temperature differences within the flow are sufficiently small, using the Taylor series of T 4 about T e and neglecting higher-order terms, T 4 can be expressed as a linear function of the temperature of the form T 4 4 T e 3 T 3 T e 4 [10]. Then, the local radiant absorption is given by the expression:
q r y = 16 α σ T e 3 ( T e T ) ,
where α is the absorption coefficient and σ is the Stefan–Boltzmann constant [4].
At the edge of the boundary layer, Howarth’s flow is considered, where the external velocity varies linearly with x, and at the edge, the velocity is u e = U 1 x L , where U is the free stream velocity. By using Bernoulli’s equation, the flow at the edge of the boundary layer and the pressure drop are associated:
d p d x = ρ u e d u e d x .
Moreover, introducing the dimensionless Falkner–Skan transformation and the dimensionless stream function [3]:
η = y u e ν x , Ψ = u e ν x f ( x , η ) ,
as well as the dimensionless quantity
θ = T w T T w T e T = T w θ ( T w T e ) ,
the following system of PDEs is obtained:
f + 1 2 1 + x u e d u e d x f f + x u e d u e d x 1 ( f ) 2 = x f f x f f x ,
1 P r θ + 1 2 1 + x u e d u e d x f θ + 16 x α σ T e 3 ρ C p u e ( 1 θ ) = x f θ x f x θ x f C p ( T w T e ) u e d u e d x + u e 2 C p ( T w T e ) ( f ) 2 ,
accompanied by the boundary conditions
η = 0 : f ( x , 0 ) = f ( x , 0 ) = 0 , θ ( x , 0 ) = 0 , η δ ˜ : f ( x , η ) = 1 , θ ( x , η ) = 1 ,
where the prime denotes partial differentiation with respect to η , P r is the Prandtl number, f = f ( x , η ) , θ = θ ( x , η ) and δ ˜ is again the boundary layer thickness, resulting now from Equation (7).
Equations (8)–(10) are still dimensional. Thus by setting
x = L x ¯ , u e = U u ¯ e ,
and dropping the · ¯ for reasons of simplicity, the following dimensionless equations are obtained
f + 1 2 1 + x u e d u e d x f f + x u e d u e d x 1 ( f ) 2 = x f f x f f x ,
1 P r θ + 1 2 1 + x u e d u e d x f θ + 16 x α σ L T e 3 ρ C p u e U ( 1 θ ) = x f θ x f x θ U 2 x f C p ( T w T e ) u e d u e d x + U 2 u e 2 C p ( T w T e ) ( f ) 2 ,
accompanied by the boundary conditions
η = 0 : f ( x , 0 ) = f ( x , 0 ) = 0 , θ ( x , 0 ) = 0 , η δ ˜ : f ( x , η ) = 1 , θ ( x , η ) = 1 .

3. Application of Homotopy Analysis Method

3.1. A Short Description of HAM

As already mentioned, one of the most powerful approximate methods for solving nonlinear differential equations (DEs) is the Homotopy Analysis Method (HAM) as proposed by Liao, see [44,45,46]. According to HAM, instead of solving the original nonlinear DE, a series of linear DEs is solved and the solution of the nonlinear DE is expressed as a series of the solutions of these linear DEs. In the following, a short description of HAM is given, whereas in Section 3.2, it is demonstrated how it can be applied to System (12)–(14).
In order to keep things simple, HAM will be demonstrated by applying it to the nonlinear ordinary differential equation of n th order
A [ y ( t ) ] = 0 ,
accompanied by the initial conditions
y ( 0 ) = a 0 , y ( 0 ) = a 1 , y ( n 1 ) ( 0 ) = a n 1 ,
where y ( t ) is the unknown function and a i , i = 0 , 1 , , n 1 are known constants.
According to HAM, an initial approximation y 0 ( t ) of y ( t ) should be constructed together with a homotopy H , so that
  • y 0 ( t ) is a solution of equation
    H [ Φ ( t ; q ) , q ] q = 0 = 0 ,
  • and y ( t ) , i.e., the solution of Equations (15) and (16), is a solution of equation
H [ Φ ( t ; q ) , q ] q = 1 = 0 ,
where q [ 0 , 1 ] is a parameter, often referred to as embedding parameter. Liao proposed the homotopy [45]:
H [ Φ ( t ; q ) , q ] = ( 1 q ) L [ Φ ( t ; q ) y 0 ( t ) ] q H ( t ) A [ Φ ( t ; q ) ] ,
where L is an auxiliary linear operator “inspired” by Equation (15) and , H ( t ) are nonzero auxiliary quantites, called the convergence control parameter and auxiliary function, respectively. It is obvious that
H [ Φ ( t ; q ) , q ] q = 0 = 0 L [ Φ ( t ; q ) y 0 ( t ) ] q = 0 = 0 Φ ( t ; 0 ) = y 0 ( t ) ,
due to the linearity of L and
H [ Φ ( t ; q ) , q ] q = 1 = 0 A [ Φ ( t ; q ) ] q = 1 = 0 ,
which is Equation (15). This means that as q varies from 0 to 1, the solution Φ ( t ; q ) of H [ Φ ( t ; q ) , q ] = 0 , varies from y 0 ( t ) to y ( t ) .
Since Φ ( t ; q ) depends not only on t but also on q, it can be expressed in the form of a series as
Φ ( t ; q ) = y 0 ( t ) + k = 1 y k ( t ) q k ,
where y k ( t ) = 1 k ! k [ Φ ( t ; q ) ] q k q = 0 satisfy the linear ODEs
L [ y k ( t ) χ k y k 1 ( t ) ] = H ( t ) R k [ y 0 ( t ) , y k 1 ( t ) , t ] , k 1 ,
called high-order deformation equations, where
χ k = 0 , k = 1 1 , k = 2 , 3 ,
and
R k [ y 0 ( t ) , , y k 1 ( t ) , t ] = 1 ( k 1 ) ! k 1 [ A ( Φ ( t ; q ) ] q k 1 q = 0 .
Thus, if all y k ( t ) can be found explicitly and the series on the right hand side of Equation (20) converges for q = 1 , the solution of Equation (15) is given in the form of the so-called homotopy series
y ( t ) = y 0 ( t ) + k = 1 y k ( t ) .
In practice, K terms of the homotopy series are calculated, in order to obtain an adequate approximation of y ( t ) . The corresponding series
y a p ( t ) = y 0 ( t ) + k = 1 K y k ( t )
is called the approximate solution of K order.
The HAM gives great freedom in choosing the initial approximation y 0 ( t ) , the convergence control parameter , the linear operator L , and the auxiliary function H ( t ) . In order to choose these functions and operators, a set of base functions, for example
S B = { e i ( t ) , i = 0 , 1 , 2 , } ,
is first chosen. This set S B is supposed to adequately describe the expected solution of Equations (15) and (16). Then, y 0 ( t ) should be chosen in such a way so that it a) can be expressed only in terms of S B and b) should satisfy the initial conditions (16), i.e.,
y 0 ( 0 ) = a 0 , y 0 ( 0 ) = a 1 , y 0 ( n 1 ) ( 0 ) = a n 1 .
The linear ODEs of Equation (21) are accompanied by the homogeneous initial conditions
y k ( 0 ) = 0 , y k ( 0 ) = 0 , y k ( n 1 ) ( 0 ) = 0 , k = 1 , 2 , .
In this way, the homotopy Series (24) satisfies Equation (16). Then L and H ( t ) are chosen in such a way so that every y k ( t ) , k = 1 , 2 , is expressed as a series of e i ( t ) , i.e.,
y k ( t ) = i = 0 c i , k e i ( t ) ,
where c i , k , i = 0 , 1 , 2 , coefficients. In this way, y ( t ) is also expressed in the form
y ( t ) = i = 0 c i e i ( t ) ,
where c i , i = 0 , 1 , 2 , coefficients.
Moreover, is chosen in such a way so that the series on the right hand side of Equation (24) is convergent. One way to choose , is to find the optimum for which the so-called “discrete squared residual" error defined by
E K 1 N + 1 j = 0 N A y a p ( t j )
is minimized, where t j = j Δ t and Δ t is the step used to discretize the interval of t using N + 1 . This has been proven to be a very good approximation of the total squared residual error; see, for example, [50] and [46] (Chapter 3).
Another way to choose is to use the so-called curves, which were proposed in [45] and have been used in a great amount of research papers. The curves are graphs of quantities of y a p ( t ) calculated at specific values of t, which usually have a physical meaning. For example, y a p ( n + 1 ) ( 0 ) may be used. If the Series (24) is convergent, y a p ( n + 1 ) ( 0 ) should always have the same value, which of course depends on . This means that in the graph of y a p ( n + 1 ) ( 0 ) against , a horizontal line segment will be apparent, which is called the valid region of . If a value of is chosen from this region, it is quite certain that the Series (24) converges. Of course, maybe more than one important quantities with physical meaning can exist. Thus, more than one curves can be used.

3.2. Application of HAM to System (12)–(14)

In order to apply HAM to System (12)–(14), an appropriate set of base functions is
S = η m e k η , m , k = 0 , 1 , 2 , .
As initial approximations, the following functions
f 0 ( x , η ) = 1 1 e δ ˜ 1 + η + e η , θ 0 ( x , η ) = 1 1 e δ ˜ 1 e η
may be considered. It should be noted that the multiplicative term 1 1 e δ ˜ , as well as the term η in f 0 ( x , η ) , are due to the boundary conditions at δ ˜ .
The corresponding linear and nonlinear operators are chosen to be
L 1 [ f ( x , η ) ] = f + f , L 2 [ θ ( x , η ) ] = θ + θ ,
N 1 [ f ( x , η ) , θ ( x , η ) ] = f + 1 2 1 + x u e d u e d x f f +
+ x u e d u e d x 1 ( f ) 2 x f f x f f x ,
N 2 [ f ( x , η ) , θ ( x , η ) ] = 1 P r θ + 1 2 1 + x u e d u e d x f θ + 16 x α σ L T e 3 ρ C p u e U ( 1 θ )
x f θ x f x θ + U 2 x f C p ( T w T e ) u e d u e d x U 2 u e 2 C p ( T w T e ) ( f ) 2 ,
where the prime denotes partial differentiation with respect to η .
Then, the higher-order deformation equations take the form
L 1 [ f k ( x , η ) χ k f k 1 ( x , η ) ] = f H f ( η ) R k f
and
L 2 [ θ k ( x , η ) χ k θ k 1 ( x , η ) ] = θ H θ ( η ) R k θ ,
accompanied by the boundary conditions
f k ( x , 0 ) = 0 . f k ( x , 0 ) = 0 , f k ( x , δ ˜ ) = 0 ,
and
θ k ( x , 0 ) = 0 . θ k ( x , δ ˜ ) = 0 ,
where k = 1 , 2 , 3 , , χ k is defined by Equation (22), the auxiliary functions H f ( η ) , H θ ( η ) are both chosen to be equal to 1, and the terms R k f and R k θ are determined via Equation (23) and are found to be:
R k f = f k 1 ( x , η ) + x u e d u e d x ( 1 χ k ) + m = 0 k 1 1 2 1 + x u e d u e d x f k 1 m ( x , η ) f m ( x , η )
x u e d u e d x f m ( x , η ) f k 1 m ( x , η ) x f m ( x , η ) f k 1 m ( x , η ) x f k 1 m ( x , η ) f m ( x , η ) x
and
R k θ = 1 P r θ k 1 ( x , η ) 16 x α σ L T e 3 ρ C p u e U θ k 1 ( x , η ) + x U 2 C p ( T w T e ) u e d u e d x f k 1 ( x , η )
+ 16 x α σ L T e 3 ρ C p u e U ( 1 χ k ) + m = 0 k 1 1 2 1 + x u e d u e d x f k 1 m ( x , η ) θ m ( x , η )
U 2 u e 2 C p ( T w T e ) f m ( x , η ) f k 1 m ( x , η ) x f m ( x , η ) θ k 1 m ( x , η ) x θ m ( x , η ) f k 1 m ( x , η ) x .
Thus, beginning with f 0 ( x , η ) , θ 0 ( x , η ) and solving Equations (31), (33) and (32), (34), several approximations of f and θ can be found. It should be mentioned that the nonhomogeneous terms x u e d u e d x and 16 x α σ L T e 3 ρ C p u e U , appearing in (12) and (13), respectively, result in the appearance of the terms x u e d u e d x ( 1 χ k ) and 16 x α σ L T e 3 ρ C p u e U ( 1 χ k ) in R k f and R k θ , respectively. These terms greatly increase the complexity of Equations (31) and (32), and as a consequence, the computational time needed for obtaining several approximations of f and θ is also considerably increased. For example, in the case where both homogeneous terms for the pressure gradient and the radiation are considered, the computational time was 50% more compared to the homogeneous case, i.e., in the absence of the pressure gradient and/or radiation effect.
Regarding the convergence of the method, the following is true
Theorem 1.
Consider the series
f ( x , η ) = f 0 ( x , η ) + k = 1 f k ( x , η ) , θ ( x , η ) = θ 0 ( x , η ) + k = 1 θ k ( x , η ) ,
where f 0 ( x , η ) , θ 0 ( x , η ) are the initial approximations previously defined and f k ( x , η ) , θ k ( x , η ) satisfy Equations (31), (33) and (32), (34), respectively. If the Series (35) converge, then they are solutions of Equations (12)–(14).
The proof of this theorem is omitted since it is similar to the proof of the general convergence theorem given in [45].
The values of the convergence control parameters f and θ are chosen in such a way so that the “discrete square residual” errors defined by
E K f 1 N + 1 j = 0 N N 1 i = 0 K f i ( η j ) ,
E K θ 1 N + 1 j = 0 N N 2 i = 0 K θ i ( η j ) ,
are minimized for both f and θ , where η j = j Δ η , Δ η = 0.16 and N = 50 . Due to the fact that θ k does not appear in Equation (31), the optimum value for f can be found independently and then used in order to find the optimum value for θ .
Table 1 shows the optimum values of f and θ , as well as the corresponding errors for several values of the order K from 2 up to 12, in the case where L = 8 , U = 1 , x = 0.5 , P r = 0.7 , T e = 300 K , a = 1 / 2 , σ = 5.67 × 10 8 , ρ = 1 , c p = 1004 and T w T e = 10 K . As the order increases, the error reduces for both equations, as expected.
The values for the convergence control parameters are in accordance with the -curves of the 15th order approximations, as shown in Figure 2. In this study, all the figures presented in the next section show results obtained utilizing the 15th order approximation for both dimensionless quantities, velocity, f , and temperature, θ . It is worth mentioning that in the case of no pressure gradient and/or no radiation effect, i.e., in the absence of nonhomogeneous terms, it suffices to consider only seven terms for f and only six terms for θ in order to obtain a good approximation of the solution.

4. Results and Discussion

To show the combined effect of a pressure gradient (adverse or favorable) and thermal radiation on the boundary layer flow, it is important to examine the dimensionless quantities of the problem under consideration, such as the dimensionless velocity, f ( η ) , and the dimensionless temperature, θ ( η ) . The shape of the boundary layer under adverse or favorable pressure gradient is clearly described by these quantities, and it turns out that it is very different from the shape of the boundary layer in the case with no pressure gradient.
In this study, the adverse pressure effect is achieved by changing the length, L ( L = 5 , 6 , 8 , 12 ) and retaining the free stream velocity, U constant. The free-stream temperature was considered to be T e = 300 K , P r = 0.7 , and c p = 1004 . The temperature difference between the plate and the free stream was considered to be, T w T e = 10 K . The radiation parameter, R, is given by the expression, R = 16 x a σ L T e 3 ρ c p u e U , where a is the absorption coefficient ( a = 1 / 2 ), σ is the Stefan–Boltzman constant ( σ = 5.67 × 10 8 ), L is a characteristic length (the length of the flat plate), u e is the velocity at the edge of the boundary layer, U is the free stream velocity, c p is the specific heat under constant pressure, and ρ is the density of the fluid [10,61,62] (see Equation (13)).
The pressure gradient effect plays a significant role on the flow, as indicated in Figure 3. The dimensionless velocity, f ( η ) , is presented in Figure 3a. Adverse pressure gradient (APG) decreases the dimensionless flow velocity, while a favorable pressure gradient (FPG) increases the dimensionless velocity compared to the case with no pressure gradient (PG). The same effect is observed for the dimensionless temperature, θ ( η ) , as shown in Figure 3b. The thermal boundary layer follows the behavior of the boundary layer. Numerically a similar problem has been studied without the effect of radiation [3,5,6].
A quantitative comparison of the analytical results obtained from HAM, with numerical results for the problem with the adverse pressure gradient, is presented in Figure 4. The numerical approach utilized is a modification of the Keller’s-box method [5], where the governing equations are reduced to a first order system and then are approximated using central differences of second order accuracy. A multi-dimensional Newton–Raphson iteration scheme is then implemented for the solution of the resulting non-linear algebraic system, and the corresponding Jacobian matrix is not prescribed by the programmer but computed numerically. It can be observed that the numerical solution is in very good agreement with the analytical one, utilizing only 20 terms from HAM.
Additionally, the most important parameters for engineering applications are the skin friction coefficient, C f x and the local Stanton number, S t x , defined per unit width of the plate, for the cases of a heating or cooling wall [5,6]. These parameters are functions of the dimensionless wall shear parameter, f w = f ( x , 0 ) , and the dimensionless wall heat-transfer parameter, θ w = θ ( x , 0 ) .
Figure 5 shows the effect of adverse pressure gradient on the skin friction coefficient, f w with the dimensionless x-distance. It is observed that with different adverse pressure gradients, the dimensionless f w provides substantially different behavior. More precisely, as the APG is increased, the flow can separate faster concerning the x-direction, as indicated in Figure 5 for different realizations of the edge fluid velocity, u e = U ( 1 x / L ) , representing the adverse pressure effect (ADV). Separation of the boundary layer, as shown in the A D V 4 case, could lead to adverse effects and boundary layer separation.
Figure 6 shows the effect of APG on the dimensionless wall heat-transfer parameter, θ w with the dimensionless x-distance. It is observed that the dimensionless wall heat-transfer parameter θ w is appreciably influenced by the adverse pressure gradient, following the boundary layer behavior. This effect is more pronounced close to the separation region, as indicated in Figure 6.
Finally, Figure 7 shows the effect of radiation on the thermal boundary layer. The flow is under the influence of an adverse pressure gradient. Close to the separation area ( x = 0.8 ), a significant cooling effect of the thermal boundary layer is observed, which is mainly due to the influence of thermal radiation. This phenomenon was also observed in previous numerical and theoretical studies regarding compressible and turbulent flows [19,20,61,62].

5. Conclusions

The main goal of this work was to study the boundary layer problem of the incompressible, laminar flow past a flat plate in the presence of a pressure gradient and radiation. The physical problem is described by the continuity, the momentum, and the energy equations with the boundary layer simplifications. Under specific assumptions on the radiation effect and using the Bernoulli equation for the pressure drop and the Falkner–Skan transformation, a nonlinear, nonhomogeneous, and coupled set of PDEs is obtained, accompanied by boundary conditions at the edge of the boundary layer.
This system is solved using the HAM, although the specific form of the boundary conditions, increase the complexity of the initial approximations f 0 ( x , η ) and θ 0 ( x , η ) , which together with the existence of the nonhomogeneous terms, increase the complexity of the higher-order deformation Equations (31) and (32). As a consequence, the computational time needed for obtaining several approximations of f and θ is considerably increased. It is indicative that in one of the cases, where both homogeneous terms for the pressure gradient and the radiation are considered, the computational time was 50% more compared to the homogeneous case. Moreover, the nonhomogeneous terms also affect the number of terms required in order to obtain a good approximation, which also increases the computational time. In some cases, only half terms are needed in the homogenous case.
To the best of the authors’ knowledge, in only one previous study were similar equations and boundary conditions solved via HAM, highlighting the novelty of the current work. Additionally, no radiation effects were taken into consideration in the aforementioned paper. A short description of this mathematical method is given together with its specific application for the problem studied.
The obtained analytical results show that the pressure gradient plays a significant role in the changes of the dimensionless velocity and temperature of the boundary layer. An adverse pressure gradient causes significant alterations on important engineering parameters, such as the dimensionless wall shear parameter, f w , and the dimensionless wall heat-transfer parameter, θ w , leading to adverse effects and flow separation. An essential cooling of the boundary layer was found, revealing that thermal behavior is influenced by radiation. This finding is important because it could be utilized in many engineering applications, such as the high temperature, gas-cooled, and gaseous fuel nuclear reactors, as discussed in the introduction section. Thermal radiation could influence the flow field and thermal boundary layer rendering the above application as a possible control technique.

Author Contributions

Conceptualization by M.A.X., E.N.P.; investigation by M.A.X., E.N.P., A.S., U.S.M.; methodology by M.A.X., E.N.P., A.S., U.S.M.; supervision by M.A.X., E.N.P.; validation by M.A.X., E.N.P., A.S.; visualization by M.A.X., E.N.P., A.S.; writing—original draft preparation by M.A.X., E.N.P., A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Mahabaleshwar (USM) wishes to thank Professor Pradeep G. Siddheshwar, F.N.A.Sc., F.I.M.A.(UK), Department of Mathematics, Bangalore University, Bangalore, India. We dedicate this work to him on the occasion of his retirement. The authors would also like to thank the reviewers for their valuable feedback.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PDE(s)Partial Differential Equation(s)
MHDMagnetohydrodynamic
HAMHomotopy Analysis Method
DE(s)Differential Equation(s)
ODE(s)Ordinary Differential Equation(s)
APGAdverse Pressure Gradient
FPGFavorable Pressure Gradient
PGPressure Gradient
ADVAdverse Pressure Effect
RADRadiation
eEdge
O x y coordinate system
x, ycoordinates
U free stream velocity
Llength of the plate
δ dimensional boundary layer thickness
u, vdimensional velocity components of the fluid
ρ density of the fluid
ppressure of the fluid
ν kinematic viscosity of the fluid
μ dynamic viscosity of the fluid
C p specific temperature of the fluid under constant pressure
Ttemperature of the fluid
kcoefficient of thermal conductivity
q r y local radiant absorption
q r radiative energy flux
T w temperature of the flat plate
T e temperature of the fluid at the edge of the boundary layer
u e velocity of the fluid at the edge of the boundary layer
α absorption coefficient
σ Stefan–Boltzmann constant
η ( x , y ) dimensionless coordinate variable
Ψ stream function
fdimensionless stream function
θ dimensionless temperature of the fluid
P r Prandtl number
x ¯ dimensionless coordinate variable
u ¯ e dimensionless velocity of the fluid at the edge of the boundary layer
δ ¯ dimensionless boundary layer thickness
H homotopy
qembedding parameter
L auxiliary linear operator
N auxiliary nonlinear operator
H ( t ) , H f ( η ) , H θ ( η ) auxiliary functions
, f , θ convergence control parameters
S B set of base functions
E K , E K f , E K θ “discrete square residual” errors
Rradiation parameter
C f x local skin friction coefficient
S t x local Stanton number

References

  1. Blasius, H. Grenzschichten in flüssigkeiten mit kleiner reibung. Z. Angew. Math. Phys. 1908, 56, 1–37. [Google Scholar]
  2. Howarth, L. On the solution of the laminar boundary layer equations. Proc. R. Soc. Lond. Ser. A 1938, 164, 547–579. [Google Scholar] [CrossRef] [Green Version]
  3. Cebeci, T.; Bradshaw, P. Physical and Computational Aspects of Convective Heat Transfer; Springer: New York, NY, USA, 1984. [Google Scholar]
  4. Howell, J.R.; Siegel, R.; Mengüç, M.P. Thermal Radiation, Heat Transfer, 5th ed.; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  5. Kafoussias, N.; Karabis, A.; Xenos, M. Numerical study of two dimensional laminar boundary layer compressible flow with pressure gradient and heat and mass transfer. Int. J. Eng. Sci. 1999, 37, 1795–1812. [Google Scholar] [CrossRef] [Green Version]
  6. Kafoussias, N.; Xenos, M. Numerical study of two-dimensional turbulent boundary layer compressible flow with adverse pressure gradient and heat and mass transfer. Acta Mech. 2000, 141, 201–223. [Google Scholar] [CrossRef]
  7. Raptis, A.; Perdikis, C.; Takhar, H.S. Effect of thermal radiation on MHD flow. Appl. Math. Comput. 2004, 153, 645–649. [Google Scholar] [CrossRef]
  8. Schlichting, H.; Gersten, K. Boundary–Layer Theory, 8th ed.; Springer: Berlin, Germany, 2000. [Google Scholar]
  9. Raptis, A.; Perdikis, C.; Leontitsis, A. Effects of radiation in an optically thin gray gas flowing past a vertical infinite plate in the presence of magnetic field. Heat Mass Transf. 2003, 39, 771–773. [Google Scholar] [CrossRef]
  10. Xenos, M. Radiation Effects on flow past a stretching plate with temperature dependent viscosity. Appl. Math. 2013, 4, 1–5. [Google Scholar] [CrossRef] [Green Version]
  11. Ali, M.M.; Chen, T.S.; Armaly, B.F. Natural convection-radiation interaction in boundary-layer flow over horizontal surfaces. AIAA J. 1984, 22, 1797–1803. [Google Scholar] [CrossRef]
  12. Raptis, A.; Toki, C.J. Thermal radiation in the presence of free convective flow past a moving vertical porous plate: An analytical solution. Int. J. Appl. Mech. Eng. 2009, 14, 1115–1126. [Google Scholar]
  13. Raptis, A.; Perdikis, C. Radiation and free convection flow past a moving plate. Appl. Mech. Eng. 1999, 4, 817–821. [Google Scholar]
  14. Raptis, A.; Perdikis, C. Thermal radiation of an optically thin gray gas. Int. J. Appl. Mech. Eng. 2003, 8, 131–134. [Google Scholar]
  15. Seddeek, M.A.; Abdelmeguid, M.S. Effects of radiation and thermal diffusivity on heat transfer over a stretching surface with variable heat flux. Phys. Lett. A 2006, 348, 172–179. [Google Scholar] [CrossRef]
  16. Muthucumaraswamy, R. The interaction of thermal radiation on vertical oscillating plate with variable temperature and mass diffusion. Theor. Appl. Mech. 2006, 33, 107–121. [Google Scholar] [CrossRef] [Green Version]
  17. Muthucumaraswamy, R.; Chandrakala, P. Effects of thermal radiation on moving vertical plate in the presence of an optically thin gray gas. Forschung im Ingenieurwesen 2005, 69, 205–208. [Google Scholar] [CrossRef]
  18. Anghaie, S.; Chen, G. Application of computational fluid dynamics for thermal analysis of high temperature gas cooled and gaseous core reactors. Nucl. Sci. Eng. 1998, 130, 361–373. [Google Scholar] [CrossRef]
  19. Duan, L.; Martin, M.P.; Feldick, A.M.; Modest, M.F.; Levin, D.A. Study of turbulence-radiation interaction in hypersonic turbulent boundary layers. AIAA J. 2012, 50, 447–453. [Google Scholar] [CrossRef]
  20. Duan, L.; Martin, M.P.; Sohn, I.; Levin, D.A.; Modest, M.F. Study of emission turbulence-radiation interaction in hypersonic boundary layers. AIAA J. 2011, 49, 340–348. [Google Scholar] [CrossRef]
  21. Miroshnichenko, I.V.; Sheremet, M.A. Numerical simulation of turbulent natural convection combined with surface thermal radiation in a square cavity. Int. J. Numer. Methods Heat Fluid Flow 2015, 25, 1600–1618. [Google Scholar] [CrossRef]
  22. Miroshnichenko, I.V.; Sheremet, M.A.; Mohamad, A.A. Numerical simulation of a conjugate turbulent natural convection combined with surface thermal radiation in an enclosure with a heat source. Int. J. Therm. Sci. 2016, 109, 172–181. [Google Scholar] [CrossRef]
  23. Kim, S.S.; Baek, S.W. Radiation affected compressible turbulent flow over a backward facing step. Int. J. Heat Mass Transfer 1996, 39, 3325–3332. [Google Scholar] [CrossRef]
  24. Jawad, M.; Shah, Z.; Islam, S.; Majdoubi, J.; Tlili, I.; Khan, W.; Khan, I. Impact of nonlinear thermal radiation and the viscous dissipation effect on the unsteady three-dimensional rotating flow of single-wall carbon nanotubes with aqueous suspensions. Symmetry 2019, 11, 207. [Google Scholar] [CrossRef] [Green Version]
  25. Khan, S.A.; Nie, Y.; Ali, B. Multiple slip effects on magnetohydrodynamic axisymmetric buoyant nanofluid flow above a stretching sheet with radiation and chemical reaction. Symmetry 2019, 11, 1171. [Google Scholar] [CrossRef] [Green Version]
  26. Jamaludin, A.; Nazar, R.; Pop, I. Mixed convection stagnation-point flow of a nanofluid past a permeable stretching/shrinking sheet in the presence of thermal radiation and heat source/sink. Energies 2019, 12, 788. [Google Scholar] [CrossRef] [Green Version]
  27. Nayfeh, A.H. Perturbation Methods; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  28. Dyke, M.V. Perturbation Methods in Fluid Mechanics; Parabolic Press: Stanford, CA, USA, 1975. [Google Scholar]
  29. Housiadas, K.D. Steady sedimentation of a spherical particle under constant rotation. Phys. Rev. Fluids 2019, 4, 103301. [Google Scholar] [CrossRef]
  30. Housiadas, K.D.; Ioannou, I.; Georgiou, G.C. Lubrication solution of the axisymmetric Poiseuille flow of a Bingham fluid with pressure-dependent rheological parameters. J. Non-Newton. Fluid Mech. 2018, 260, 76–86. [Google Scholar] [CrossRef]
  31. Housiadas, K.D.; Tanner, R.I. A high-order perturbation solution for the steady sedimentation of a sphere in a viscoelastic fluid. J. Non-Newton. Fluid Mech. 2016, 233, 166–180. [Google Scholar] [CrossRef]
  32. Housiadas, K.D.; Tanner, R.I. Viscoelastic shear flow past an infinitely long and freely rotating cylinder. Phys. Fluids 2018, 30, 073101. [Google Scholar] [CrossRef]
  33. Si, X.; Yuan, L.; Cao, L.; Zheng, L.; Shen, Y.; Li, L. Perturbation solutions for a micropolar fluid flow in a semi-infinite expanding or contracting pipe with large injection or suction through porous wall. Open Phys. 2016, 14, 231–238. [Google Scholar] [CrossRef]
  34. Zhang, Y.; Lin, P.; Si, X.-H. Perturbation solutions for asymmetric laminar flow in porous channel with expanding and contracting walls. Appl. Math. Mech.-Engl. Ed. 2014, 35, 203–220. [Google Scholar] [CrossRef]
  35. Liu, Y.; Kurra, S.N. Solution of Blasius equation by variational iteration. Appl. Math. 2011, 1, 24–27. [Google Scholar] [CrossRef] [Green Version]
  36. Xu, L.; Lee, E.W.M. Variational iteration method for the magnetohydrodynamic flow over a nonlinear stretching sheet. Abstr. Appl. Anal. 2013, 2013, 573782. [Google Scholar] [CrossRef]
  37. Abbasbandy, S. A numerical solution of Blasuis equation by Adomians decomposition method and comparison with homotopy perturbation method. Chaos Solitons Fractals 2007, 31, 257–260. [Google Scholar] [CrossRef]
  38. Aski, F.S.; Nasirkhani, S.J.; Mohammadiand, E.; Asgari, A. Application of Adomian decomposition method for micropolar flow in a porous channel. Prop. Power Res. 2014, 3, 15–21. [Google Scholar] [CrossRef] [Green Version]
  39. Kuo, B.-L. Application of the differential transformation method to the solutions of Falkner-Skan wedge flow. Acta Mech. 2003, 164, 161–174. [Google Scholar] [CrossRef]
  40. Thiagarajan, M.; Senthilkumar, K. DTM-Padé approximants for MHD Flow with suction/blowing. J. Appl. Fluid Mech. 2012, 6, 724–728. [Google Scholar]
  41. Rebenda, J.; Šmarda, Z. Convergence analysis of an iterative scheme for solving initial value problem for multidimensional partial differential equations. Comp. Math. Appl. 2015, 70, 1772–1780. [Google Scholar] [CrossRef]
  42. Rebenda, J.; Šmarda, Z. A differential transformation approach for solving functional differential equations with multiple delays. Commun. Nonlinear Sci. Numer. Simul. 2017, 48, 246–257. [Google Scholar] [CrossRef] [Green Version]
  43. Rebenda, J.; Šmarda, Z. Numerical algorithm for nonlinear delayed differential systems of nth order. Adv. Differ. Equ. 2019, 2019, 26. [Google Scholar] [CrossRef] [Green Version]
  44. Liao, S.J. The Proposed Homotopy Analysis Technique for the Solution of Nonlinear Problems. Ph.D. Thesis, Shanghai Jiao Tong University, Shanghai, China, 1992. [Google Scholar]
  45. Liao, S.J. Beyond Perturbation: Introduction to the Homotoy Analysis Method; Chapman & Hall/CRC, CRC Press LLC: Boca Raton, FL, USA, 2004. [Google Scholar]
  46. Liao, S.J. Homotopy Analysis Method in Nonlinear Differential Equations; Springer: Berlin, Germany, 2011. [Google Scholar]
  47. Liao, S.J. A kind of approximate solution technique which does not depend upon small parameters-II. An application in fluid mechanics. Int. J. Nonlin. Mech. 1997, 32, 815–822. [Google Scholar] [CrossRef]
  48. Liao, S.J. An explicit totally analytic solution of laminar viscous flow over a semi-infinite flat plate. Commun. Nonlinear Sci. Numer. Simul. 1999, 3, 53–57. [Google Scholar] [CrossRef]
  49. Liao, S.J. A uniformly valid analytic solution of 2D viscous flow past a semi-infinite flat plate. J. Fluid Mech. 1999, 385, 101–128. [Google Scholar] [CrossRef] [Green Version]
  50. Farooq, U.; Hayat, T.; Alsaedi, A.; Liao, S.J. Series solutions of non-similarity boundary layer flows of nano-fluids over stretching surfaces. Numer. Algor. 2015, 70, 43–59. [Google Scholar]
  51. Jawad, M.; Shah, Z.; Kahn, A.; Kumam, P.; Islam, S. Entropy generation and heat transfer analysis in MHD unsteady rotating flow for aqueous suspensions of carbon nanotubes with nonlinear thermal radiation and viscous dissipation effect. Entropy 2019, 21, 492. [Google Scholar] [CrossRef] [Green Version]
  52. Naz, R.; Noor, M.; Hayat, T.; Javed, M.; Alsaed, A. Dynamism of magnetohydrodynamic cross nanofluid with particulars of entropy generation and gyrotactic motile microorganisms. Int. Commun. Heat Mass Transf. 2020, 2020 110, 104431. [Google Scholar] [CrossRef]
  53. Ullah, A.; Alzahrani, E.O.; Shah, Z.; Ayaz, M.; Islam, S. Nanofluids thin film flow of Reiner-Philippoff fluid over an unstable stretching surface with Brownian motion and thermophoresis effects. Coatings 2019, 9, 21. [Google Scholar] [CrossRef] [Green Version]
  54. Ullah, A.; Shah, Z.; Kumam, P.; Ayaz, M.; Islam, S.; Jameel, M. Viscoelastic MHD nanofluid thin film flow over an unsteady vertical stretching sheet with entropy generation. Processes 2019, 7, 262. [Google Scholar] [CrossRef] [Green Version]
  55. Waqas, M.; Dogonchi, A.S.; Shehzad, S.A.; Khan, M.I.; Hayat, T.; Alsaedi, A. Nonlinear convection and joule heating impacts in magneto-thixotropic nanofluid stratified flow by convectively heated variable thicked surface. J. Molecular Liquids 2020, 300, 1111945. [Google Scholar] [CrossRef]
  56. Hussain, M.; Ashraf, M.; Nadeem, S.; Khan, M. Radiation effects on the thermal boundary layer flow of a micropolar fluid towards a permeable stretching sheet. J. Franklin Inst. 2013, 350, 194–210. [Google Scholar] [CrossRef]
  57. Sharma, R.P.; Seshadri, R.; Mishra, S.R.; Munjam, S.R. Effect of thermal radiation on magnetohydrodynamic three-dimensional motion of a nanofluid past a shrinking surface under the influence of a heat source. Heat Transfer-Asian Res. 2019, 48, 2105–2121. [Google Scholar] [CrossRef]
  58. Sohail, M.; Naz, R.; Abdelsalam, S.I. Application of non-Fourier double diffusions theories to the boundary-layer flow of a yield stress exhibiting fluid model. Physica A 2020, 537, 122573. [Google Scholar] [CrossRef]
  59. Sravanthi, C.S. Effect of heat source/sink on MHD flow due to porous rotating disk with variable thickness in the presence of cross diffusion. Heat Transf.-Asian Res. 2019, 48, 4016–4032. [Google Scholar] [CrossRef]
  60. Nargund, A.L.; Madhusudhan, R.; Sathyanarayana, S.B. Study of compressible fluid flow in boundary layer region by homotopy analysis method. Int. J. Latest Trend Eng. Technol. 2017, 9, 28–39. [Google Scholar]
  61. Xenos, M. Thermal radiation effect on the MHD turbulent compressible boundary layer flow with adverse pressure gradient, heat transfer and local suction. Open J. Fluid Dyn. 2017, 7, 1–14. [Google Scholar] [CrossRef] [Green Version]
  62. Xenos, M.; Pop, I. Radiation effect on the turbulent compressible boundary layer flow with adverse pressure gradient. Appl. Math. Comput. 2017, 299, 153–164. [Google Scholar] [CrossRef]
Figure 1. Flow configuration and coordinate system.
Figure 1. Flow configuration and coordinate system.
Symmetry 12 00710 g001
Figure 2. (a) curve for the dimensionless velocity, f ( η ) for 15 terms. (b) curve for the dimensionless temperature, θ ( η ) for 15 terms.
Figure 2. (a) curve for the dimensionless velocity, f ( η ) for 15 terms. (b) curve for the dimensionless temperature, θ ( η ) for 15 terms.
Symmetry 12 00710 g002
Figure 3. (a) The dimensionless velocity, f ( η ) with respect to η . (b) The dimensionless temperature, θ ( η ) with respect to η , for the case where L = 8 , for adverse pressure gradient (APG) and favorable pressure gradient (FPG) cases.
Figure 3. (a) The dimensionless velocity, f ( η ) with respect to η . (b) The dimensionless temperature, θ ( η ) with respect to η , for the case where L = 8 , for adverse pressure gradient (APG) and favorable pressure gradient (FPG) cases.
Symmetry 12 00710 g003
Figure 4. Comparison of the dimensionless velocity, f ( η ) with respect to the dimensionless distance η for the analytical (HAM, 20 terms), and the numerical solution obtained with the numerical approach introduced in [5].
Figure 4. Comparison of the dimensionless velocity, f ( η ) with respect to the dimensionless distance η for the analytical (HAM, 20 terms), and the numerical solution obtained with the numerical approach introduced in [5].
Symmetry 12 00710 g004
Figure 5. Skin friction coefficient, f w with respect to the dimensionless x-direction.
Figure 5. Skin friction coefficient, f w with respect to the dimensionless x-direction.
Symmetry 12 00710 g005
Figure 6. Heat transfer coefficient, θ w with respect to the dimensionless x-direction.
Figure 6. Heat transfer coefficient, θ w with respect to the dimensionless x-direction.
Symmetry 12 00710 g006
Figure 7. Effect of radiation on the dimensionless temperature, θ ( η ) , depending on the parameter, R, and the adverse pressure gradient, with respect to η .
Figure 7. Effect of radiation on the dimensionless temperature, θ ( η ) , depending on the parameter, R, and the adverse pressure gradient, with respect to η .
Symmetry 12 00710 g007
Table 1. Optimum values of f , θ , and corresponding errors.
Table 1. Optimum values of f , θ , and corresponding errors.
Order f -error f θ -error θ
20.00239722−0.4854580.0021427−0.637823
40.000427796−0.8032830.000710738−0.596724
60.0000874967−0.8246880.000489771−0.63417
80.0000507706−0.8921260.00045183−0.653759
100.0000637655−0.874670.000441161−0.681927
120.0000661033−0.6510430.000436139−0.68495

Share and Cite

MDPI and ACS Style

Xenos, M.A.; Petropoulou, E.N.; Siokis, A.; Mahabaleshwar, U.S. Solving the Nonlinear Boundary Layer Flow Equations with Pressure Gradient and Radiation. Symmetry 2020, 12, 710. https://doi.org/10.3390/sym12050710

AMA Style

Xenos MA, Petropoulou EN, Siokis A, Mahabaleshwar US. Solving the Nonlinear Boundary Layer Flow Equations with Pressure Gradient and Radiation. Symmetry. 2020; 12(5):710. https://doi.org/10.3390/sym12050710

Chicago/Turabian Style

Xenos, Michalis A., Eugenia N. Petropoulou, Anastasios Siokis, and U. S. Mahabaleshwar. 2020. "Solving the Nonlinear Boundary Layer Flow Equations with Pressure Gradient and Radiation" Symmetry 12, no. 5: 710. https://doi.org/10.3390/sym12050710

APA Style

Xenos, M. A., Petropoulou, E. N., Siokis, A., & Mahabaleshwar, U. S. (2020). Solving the Nonlinear Boundary Layer Flow Equations with Pressure Gradient and Radiation. Symmetry, 12(5), 710. https://doi.org/10.3390/sym12050710

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop