Open access peer-reviewed chapter - ONLINE FIRST

Padé Approximation to Solve the Problems of Aerodynamics and Heat Transfer in the Boundary Layer

By Igor Andrianov and Anatoly Shatrov

Submitted: March 13th 2020Reviewed: June 1st 2020Published: September 4th 2020

DOI: 10.5772/intechopen.93084

Downloaded: 20

Abstract

In this chapter, we describe the applications of asymptotic methods to the problems of mathematical physics and mechanics, primarily, to the solution of nonlinear singular perturbed problems. We also discuss the applications of Padé approximations for the transformation of asymptotic expansions to rational or quasi-fractional functions. The applications of the method of matching of internal and external asymptotics in the problem of boundary layer of viscous gas by means of Padé approximation are considered.

Keywords

  • asymptotic methods
  • Padé approximation
  • boundary-value problem of mathematical physics
  • boundary layer

1. Introduction

An important drawback of asymptotic methods is the local character of solutions obtained [1, 2, 3, 4]. Since the constructed series are often asymptotic, a simple increase in the number of terms does not remove this drawback. Essence of the problem consists of divergence of obtained series. There exist a lot of approaches to these problems [5, 6]. The method of analytic continuation (e.g., the Euler transform or generalized Euler transform [7, 8, 9, 10, 11, 12]) requires a priori information about the singularities of the searched function in the complex domain [4, 9]. These methods are useful if a large number of terms of the series are known. In this case, it is possible to use the Domb-Sykes plot [5, 8]. But usually only a few terms of asymptotic series are known, and to get information from them, the method of Padé approximations (PAs) is useful [1, 2, 5, 13, 14, 15]. PAs yield meromorphic continuations of functions defined by power series and can be used even in cases where analytic continuations are inapplicable. If a PAs converges to the given function, then roots of the denominator tend to points of singularities. One-point PAs give possibilities to improve convergence of series [16, 17, 18, 19, 20]. Two-point PAs (TPPAs) allow matching asymptotics in transition zones and are widely used in mechanics and physics [1, 2, 4, 14, 21, 22, 23, 24]. Overcoming the mentioned limitations of asymptotic methods for practically important problem is the purpose of this chapter. We consider at the beginning (Section 2) the mathematical bases of asymptotic methods and the use of Padé approximants for the summations of the asymptotic series. Section 3 discusses the method of combining of internal and external asymptotics (matching method) by means of Padé approximants. In the Section 4, the methods of solving specific problems of mathematical physics and mechanics of fluid and gas are demonstrated. Section 5 presents a discussion of the obtained.

2. Mathematical background: summation of asymptotic series

2.1 Analysis of power series

We suppose that by the result of the asymptotic study, one obtains the following series:

fεn=0Cnεnfor        ε0.E1

As is known, the radius of convergence ε0series (1) is determined by the distance to the nearest singularity of the function f(ε) on the complex plane. To define ε0, the Domb-Sykes plot may be useful [8, 10]. In many cases, one can effectively use the conformal mapping of the series, a fairly complete catalog of which is given in [9]. In particular, it sometimes turns out to be a successful Euler transformation [8, 10], based on the introduction of a new variable:

ε˜=ε1ε/ε0.E2

Recast the function f in terms of ε˜, fn=0dnε˜n, transfer the singularity at the point ε˜=.

A natural generalization of Euler transformation looks as follows:

ε˜=ε1ε/ε0α,

where αis the certain number.

2.2 Padé approximants

“The coefficients of the Taylor series in the aggregate have a lot more information about the values of features than its partial sums. It is only necessary to be able to retrieve it, and some of the ways to do this is to construct a Padé approximant” [11]. Padé approximants (PAs) allow us to transform of power series to a fractional-rational function. Let us define PAs, following Baker and Graves-Morris [25].

Suppose we are given the power series:

fε=i=1ciεi,E3

PAs can be written as the following expression:

fn/mε=a0+a1ε++anεn1+b1ε++bmεm,E4

whose coefficients are determined from the condition

1+b1ε++bmεmc0+c1ε+c2ε2+=a0+a1ε++anεn+Oεn+m+1E5

Equating coefficients near the same powers ε, one obtains a system of linear algebraic equations. In the case where this system is solvable, one can obtain the Padé coefficients of the numerator and denominator of the PAs.

We note some properties of the PAs [5, 13, 19]. If the PAs at the chosen m and n exists, then it is unique.

  1. If the PAs sequence converges to some function, the roots of its denominator tend to the poles of the function. This allows for a sufficiently large number of terms to determine the pole and then perform an analytical continuation.

  2. PAs gives meromorphic continuation of a given power series.

  3. PAs of the inverse function is treated as the PAs function inverse itself. This property is called duality and is more exactly formulated as follows. Let

    qε=f1εandf00,thenqn/mε=fn/m1εE6

  • Diagonal PAs are invariant under fractional linear transformations of the argument. Suppose that the function is given by their expansion (3). Consider the linear fractional transformation that preserves the origin W=/1+and the function q(W) =fε. Then qn/nW=fn/nε, provided that one of these approximations exist. In particular, the diagonal PAs is invariant concerning Euler transformation (2).

  • Diagonal PAs are invariant under fractional linear transformations of functions. Let us analyze a function (3). Let

    qε=a+bfεc+dfε.

    If c+df00, then

    qn/nε=a+bfn/nεc+dfn/nε

    provided that there is fn/nε.

  • PAs can get the upper and lower bounds for fn/nε. For the diagonal PAs, one has the following estimate:

    fn/n1εfn/nεfn/n+1ε.E7

    Typically, this estimate is valid for the function itself, that is, fn/nεin Eq. (7) can be replaced by fε.

  • Diagonal and close to them a sequence of PAs often possesses the property of autocorrection [17, 18]. It consists of the following. To determine the coefficients of the numerator and denominator of PAs, we have to solve systems of linear algebraic equations. This is an ill-posed procedure, so the coefficients of PAs can be determined with large errors. However, these errors in a certain sense are of self-consistent, so the accuracy of PAs is high. This is the radical difference the PAs from the Taylor series, the calculation error of which only increases with increasing number of terms.

  • Autocorrection property is verified for a number of special functions. At the same time, even for elliptic functions, the so-called Froissart doublets phenomenon arises [26]. Thus, in general, having no information about the location of the poles of the PAs, but relying solely on the very PAs (computed exactly as you wish), we cannot say that you have found a good approximated function. Now consider the question: In what sense the available mathematical results on the convergence of the PAs can facilitate the solution of practical problems? Gonchar’s theorem [16] states: If none of the diagonal PAs fn/nεhas poles in the circle of radius R, then the sequence fn/nεis uniformly convergent in the circle to the original function fε. Moreover, the absence of poles of the sequence of the fn/nεin a circle of radius R confirms convergence of the Taylor series in the circle. Since the diagonal PAs is invariant under fractional linear maps εε/+b, the theorem is true for any open circle containing the point of decomposition, and for any area, which is the union of these circles. A significant drawback in practice is the need to check all diagonal PAs. The fact is that if a circle of radius R has no poles only for a subsequence of the sequence of diagonal PAs, then the uniform convergence to its original holomorphic in the disk is guaranteed only with r<r0, where 0.583<r0<0.584[27]. How can we use these results? Suppose that there are a few terms of the perturbation series and one wants to estimate its radius of convergence R. Consider the interval [0,ε0], where the truncated perturbation series and the diagonal PAs of the maximal possible order differ by no more than 5% (adopted in the engineering accuracy of the calculations). If none of the previous diagonal PAs does not have in a circle of radius ε0poles, then it is a high level of confidence to assert that Rε0.

    3. Matching of limiting asymptotic expansions

    3.1 Method of asymptotically equivalent functions

    This method was originally proposed by Slepyan and Yakovlev for the inversion of the integral transformations. Here is a description of this method, following [26].

    Suppose that the Laplace transform of a function of a real variable f(t) is

    Fs=0ftestds.

    To obtain an approximate expression for the inverse transform, it is necessary to clarify the behavior of the transform to the vicinity of the points s = 0 and s = ∞ and to determine whether the nature and location of its singular points are on the exact boundary of the regularity or near it. Then the transform F(s) is replaced by the function F0(s), approximated the exact inversion and satisfying the following conditions:

    1. Functions F0(s) and F(s) are asymptotically equivalent at s → ∞ and s → 0, that is,

      F0sFsats0ands.

  • Singular points of functions F0(s) and F(s), located on the exact boundary of the regularity, coincide.

  • The free parameters of the function F0(s) are chosen so as to satisfy the conditions of the good approximation of F(s) in the sense of minimum relative error for all real values s ≥ 0:

    maxF0sα1α2αkFs1minE8

    Condition (8) is achieved by variation of free parameters αk. Often the implementation of equalities

    0F0sds=0Fsds

    or F0sFsat s → 0 leads to a rather precise fulfillment of the requirements (8).

    Constructed in such a way function F0(s) is called asymptotically equivalent function for F(s) (AEF). Let’s dwell on the terminology. In the following sections, we will use the symbols of ordinal relations. We will give strict definitions of these concepts.

    Let’s consider the function f(x). To describe the ordinal relationships with respect to another function φ(x), enter the following definitions:

    Definition 1. Let us say that f(x) is a value of order φ(x) at xx0, that is,

    fx=Oφx

    if δ>0A:xx0<δfxAφx.

    Definition 2. Let us say that f(x) is a value of order less than φ(x) at xx0, that is,

    fx=oφx

    if δ>0ε:xx0<δfxεφx.

    Here A is a finite number, and ε, δ are infinitely small.

    Definition 3. Let us say that f (x) is asymptotically equal to φ(x) at xx0, that is,

    fxφx if fxφx1.

    Here we use the term “asymptotically equivalent function.” Other terms (“reduced method of matched asymptotic expansions” [28], “quasi-fractional approximants” (QAs) [29], and “mimic function” [30]) are also used.

    3.2 Two-point Padé approximants

    The analysis of numerous examples confirms “complementarity principle”: if for ε0, one can construct a physically meaningful asymptotics, there is a nontrivial asymptotics and for ε. The most difficult in the asymptotic approach is the intermediate case of ε1. In this domain, typically numerical methods work well; however, if the task is to investigate the solution depending on the parameter ε, then it is inconvenient to use different solutions in different areas. Construction of a unified solution on the basis of limiting asymptotics is not a trivial task, and for this purpose, one can use a two-point Padé approximants (TPPAs). We give the definition following [25]. Let

    Fε=i=0ciεiatε0,E9
    Fε=i=0diεiatεE10

    TPPA is a rational function of the form:

    fn/mε=a0+a1ε++anεn1+b1ε++bmεm,E11

    k coefficients which are determined from the condition

    1+b1ε++bmεmc0+c1ε+c2ε2+=a0+a1ε++anεn+Oεn+m+1E12

    and the remaining coefficients from a similar condition for ε1.

    4. Application of Padé approximants

    4.1 Using of TPPAs in boundary-value problems

    For boundary-value problems, we assume that there exist two asymptotics for limit values of the parameter. In this case, the method of matching of asymptotic expansions is usually used [4]. However, for correct application of the matching method, it is necessary to know the matching point or, at least, the domain of overlapping of asymptotics. An exact description of the transition layer 0 <ε< ∞ exists only in the cases where solutions with different behaviors on opposite sides of the layer can be matched by a special function (e.g., the Airy function).

    For the matching of nonoverlapping asymptotics, a method based on TPPAs has recently been developed. In [15, 21, 23], this method was applied for the construction of thermal profiles in a boundary layer of gas. In [2, 6], this method allowed one to examine the heat exchange in hypersonic boundary layers.

    Two-point Padé approximations (TPPAs) are defined in Section 3.2 [see formulas (2)(4)]. As an example of application of TPPAs, we consider the Airy boundary-value problem [4, 10, 31]:

    yλ2xy=gxyasλE13

    with boundary conditions

    y0=1,y=0E14

    This boundary-value problem has the form in terms of Airy function U(s):

    UsU=0,U0=1,U=0E15

    The asymptotic solution for problems (13) and (14) has the form:

    yx=Us1+Oλ1ass=xλ2/3.E16

    The interior asymptotic (s → 0) has the form of a power function:

    Ui=1as+16s3+Os4E17

    The exterior asymptotic has the form of an exponential function:

    Ue=bs14exp23s121548s32+Os3E18

    as a0.7290, b0.7946.

    The transition layer is defined by the domain, where x=Oλ23

    Airy function approaches with TPPA:

    Ua=1as+23s3223as52+325as41+325abs124exp23s32E19

    The TPPA (19) preserves three terms of the asymptotics at both ends and provides accuracy with relative error:

    Δ=UUaU1.5%

    Parameters a and b are obtained from the integral equations (relations). The relations (20) and (21) can be obtained by multiplying Eq. (18) by 1, s, s2,and then by integrating from 0 to ∞.

    U=sU0Uds=0sUds0Uds=0sUdsUo=0sUds

    This is the first integral relation.

    0sUds=aE20
    sU=s2U0sUds=0s2Udss=t,dt=dsUds=dV,V=U'U00Uds=0s2Uds0s2Uds=1

    This is the next integral relation.

    0s2Uds=1E21

    Substituting in Eqs. (20) and (21) instead of U (4) interpolation Ua (7), calculate using quadrature integration formulas a = 0.7287 and b = 0.7922.

    In the same manner, integral relations with weights U, U’ can be obtained by part integration. Multiplying Eq. (18) by U,U,U, we get after integration from 0 to ∞,

    UU=sU20UUds=0sU2dsU=t,dt=dUUds=dV,V=UUU200U2ds=0sU2dsa0U2ds=0sU2ds0U2+sU2ds=a

    This is the first integral relation for the second method of producing it:

    0U2+sU2ds=aE22
    UU=sUU0UUds=0sUUdsU=t,dt=UdsUds=dV,V=UU200UUds=0sUUdsa2=20sUUdss=t,dt=dsUUds=dV,V=U22a2=2sU2200U22ds0sU2ds0U2ds=a2

    And this is the next integral relation for the second method of producing it:

    0U2ds=a2E23

    Using Eq. (19), from Eqs. (22) and (23), we calculate а = 0.7277, and b = 0.7966.

    From the given example, it follows that the features of the asymptotic connection method are the ambiguity of the algorithm, the freedom to choose both the form of TPPAs, integral relations, and methods for calculating the parameters of the TPPAs. The question of choosing integral relations is, in fact, a question of controlling the asymptotic approximation using weights selected to obtain integral relations. Choosing the weight allows you to achieve acceptable accuracy in a particular area of the boundary layer: a weight equal to 1 means that the uniform influence of the entire layer is taken into account; a weight equal to 1,s,s2,increases the influence of the outer region of the layer; and if the desired solution U,U,Uis chosen as the weight, then its inhomogeneity increases the influence of the local region where the inhomogeneity is concentrated.

    4.2 Quasi-fractional Padé approximants (modification of TPPA)

    In the illustrated example (5), Eq. (18) TPPA represents a modified (quasi-fractional) two-point Padé approximant (10) by an exponential weight function, the choice of which is dictated by a kind of exterior asymptotics. Evidently, the TPPAs are not panacea. For example, one of the “bottlenecks” of the TPPAs method is related to the presence of logarithmic components in numerous asymptotic expansions. This problem is the most essential for the TPPAs, because, as a rule, one of the limits ε0orεfor a real mechanical problem gives expansions with logarithmic terms or other complicated functions. It is worth noting that in some cases these obstacles may be overcome by using an approximate method of TPPAs’ construction by tacking as limit points not ε=0andε=, but some small and large values. On the other hand, Martin and Baker [32] proposed the so-called quasi-fractional approximants (QAs). Let us suppose that we have a perturbation approach in powers of εfor ε0and asymptotic expansions Fεcontaining, for example, logarithm for ε. By definition, QA is a ratio R with unknown coefficients ai,bi, containing both powers of εand Fε. We give this modification of TPPA [2, 14, 15]. Let the series give for Eq. (5). Then the modification of TPPA is represented by the irrational function:

    Fε=k=0makεkk=0nbkεkexpk=0lckεk,E24

    where k+1coefficients ck,k=012are determined by means of l+1integral equations for function from Eqs. (20) and (21). We notice that exponential terms [multiplier in expressions (17) and (18)] give for ε=0andε=coincidence with TPPA (19). When considering the computational aspects of the connection method, it should first be assumed that the system of equations for determining the TPPA parameters is substantially nonlinear. To solve it, we developed a modification of the method of solving nonlinear algebraic systems [4, 23, 24].

    4.3 Application of TPPAs in problems of incompressible liquid and gas mechanics

    Consider the Blasius equation (45), which describes laminar boundary layers on a flat plate:

    φ+φφ=0;φ0=φ0=0;φ=2E25

    where φζ=ψ/x, ψyis the stream function, ζ=y2Rexis the automodel variable, and x and y are the Cartesian coordinates such that the axis x is directed along the flow. The interior asymptotic (ζ → 0) has the form:

    φ=a2ζ2a2230ζ5+Oζ8E26

    The procedure for obtaining external asymptotics is nontrivial due to the presence of logarithmic components in the main elements. We describe in detail the mechanism for obtaining and evaluating both primary and secondary members of asymptotic. From Eq. (25) follows:

    φφ=φE27

    After integration of Eq. (27) by the coordinate ζ follows:

    lnφζ0ζ=0ζφdζlnφζln2a2=0ζφdζφζ=2a2exp0ζφdζE28

    After reintegration of Eq. (28) by the coordinate

    0ζφζ1dζ1=0ζ2a2exp0ζ4φdζ2dζ1

    follows:

    φζ10ζ=0ζ2a2exp0ζ4φdζ2dζ1

    subject to boundary conditions

    φζ1=0ζ2a2exp0ζ1φdζ2dζ1φζ1=0ζ2a2φζ1exp0ζ1φζ2dζ2φζ1dζ1φζ1=2a20ζ1φζ1dexp0ζ1φζ2dζ2

    Let us make a limit transition ζin the last equation and represent the integration interval as

    0=0ζζfollows:

    φζ=2+2a2ζ1φζ1dexp0ζ1φζ2dζ2

    We use the mean theorem in the last equation

    φζ=2+2a21φζexp0ζ1φζ2dζ2E29

    In the resulting equation, the first compound is the principal member of the external asymptotics. To obtain the following members of the asymptotic, we will present the function as

    φ=2ζc+z

    where z0, if ζ. Given the last expression of the function φ, Eq. (29) is obtained as follows:

    φζ=2+2a22ζc+zexp0ζ2ζ1c+zdζ1

    If z=φ2ζ+c, then

    φζ=2+2a22ζc+zexpζ2+cζexp0ζφ2ζ1+cdζ1

    In the external domain, where ζand z0, let us receive an exterior asymptotic:

    φζ=2+2a2D2ζcexpζ2+cζ+o1ζ2E30

    where D=exp0φ2ζ+c.

    To calculate parameter a2, use the procedure of Section 4.1 [see formula (20)], and using weight equal to 1:

    a2=120φa2φE31

    At that, in external domain, ζ

    φ=2ζc,z0

    Therefore,

    c=02φE32

    and

    D=exp0φ2ζ+cE33

    Type of generalized and normalized TPPA of order (4,4):

    φaζ=211+α1ζ+α2ζ2+α3ζ3+α4ζ4expζ2+1+β1ζ+β2ζ2+β3ζ3+β4ζ4E34

    Species of TPPA taking into account four nontrivial parameters:

    α3,β1,β2,β4

    Therefore,

    φaζ=211+α3ζ3expζ2+1+β1ζ+β2ζ2+β4ζ4E35

    Parameter values are determined using local asymptotic and TPPA in the respective domain. Taking into account the decomposition of the exponent in the internal domain, we will write down the local equality:

    2a2ζa226ζ4=211+α3ζ31+ζ2+ζ222+1+β1ζ+β2ζ2+β4ζ4E36

    Taking into account Eq. (33) in the external domain, we will write down the second local equality:

    22a2D2ζcexpζ2+=211+α3ζ3expζ2+1+β1ζ+β2ζ2+β4ζ4a2D1+β1ζ+β2ζ2+β4ζ4=2ζc1+α3ζ3E37

    Equalizing the coefficients in Eqs. (36) and (37) at the same degrees ζ, we get

    α3=a2D,β1=a2+c,β2=c22+a2a2+c,β4=2.

    Therefore, the TPPA has the form:

    φaζ=211+a2Dζ3expζ2+1+a2+cζ+a22+a2c+c221ζ2+2ζ4E38

    After systems (31)–(33) are solved, we will obtain

    a2=0.6641,c=1.7308,D=0.3357E39

    By substituting (39) in (38), we get an explicit expression for the TPPA.

    4.4 Combining method of interior and exterior asymptotics for boundary layer of supersonic flow in compressed viscous gas by TPPA

    We consider the boundary layer in hypersonic flow of viscous gas and solve a model problem which reduces to ordinary differential equations with appropriate boundary conditions. The TPPAs parameters are calculated and relevant questions are discussed. The equations of laminar boundary layer near a semi-infinite plate in the supersonic flow of viscous perfect gas, as it is known [2, 7], can be reduced to the form:

    φμT+φφ=0,E40
    μTT+σφT+μTφ2=0E41

    where

    φ=ψx=φζ,T=Tζ,ζ=η2x,η=0ydyT,a=14M2κ1

    M is the Mach number, σis the Prandtl number,κis the adiabatic index, ψis the stream function, Tis the temperature, μis the viscosity coefficient, and x and y are the Cartesian coordinates.

    The boundary conditions at the wall are

    φ0=φ0=0,T0=TsE42

    At external boundary of layer is

    φ=2,T=1.E43

    Interior asymptotic expansions are for μ=Tn

    φ=2a2ζn1a2T1Tsζ2+Oζ3,T=Ts+T1ζ2a22+n12T12Tsζ2+Oζ3E44

    where two constants a2and T1remain undefined.

    Exterior asymptotics for ς

    lnφ=c2++lnA+o1,lnT=σζ2+σcζ+lnB+o1E45

    where three constants are unknown: c, A, and B.

    We solve boundary problems (40) and (41) approximately by connecting asymptotics (44) and (45) TPPA

    φaζ=211+Aζ3expζ2+1+α1ζ+α2ζ2+α4ζ4E46
    Taζ=ζmζβ0+β1ζexpσζ2+E47

    Boundary conditions (45) and (46) are satisfied if to put

    φaζ=0φaζ,Taζ=Ts+0TaζE48

    We complement the last equalities (50) and (51) with a normalizing condition:

    1=Ts+0TaζE49

    Following the procedure of the previous section, we will calculate the coefficients at ζ and ζ2 in asymptotic expansions (44) and, equating them with the corresponding expressions from Eqs. (46) and (47), we will obtain equalities, from which values α1,α2,α4,β0,β1are expressed through a1,c,T1,ζm:

    α1=a2+c,α2=12n1a2T1Ts1+c22+a2a2+c,α4=4,β0=ζmT1,β1=σcζmT11T1+n1ζmT12E50

    Three parameters in asymptotics (44) are defined in the outer region if the following condition is met:

    β1=1BE51

    A priori at large M numbers, it is known that the temperature profile is non-monotonic and has a maximum within the layer at point ςmat which, as can be seen from the second equation of the systems (40) and (41), the following condition is used:

    Tζm=φ2ζmE52

    From the convexity condition of the temperature profile in the vicinity of the point ςm, the following equality is used:

    β0+β1ζmφa2ζm=expσζm2ζmE53

    Let us add the received equations with the integrated ratios received on the basis of coincidence of TPPAs (46) and (47); in this case, three members in asymptotic decompositions (50) and (51), the initial system of Eqs. (40) and (41), with boundary conditions (42) and (43), by using the technique stated in the previous sections.

    a2=120φa2φc=02φE54

    The integral relation for parameter A is obtained by multiplying Eq. (40) by

    expζ2

    and integrating from 0 to ∞ taking into account Eq. (48):

    A=2a2μTsTs0φμTT2ζcφexpζ2E55

    Similarly, from Eq. (41), we get

    B=σ0TφμTT2ζc+a2μTTφ2expσζ2E56

    Thus, the integral relations (52) and (55)–(47) form a nonlinear system of equations for determining the following parameters:

    T1,a2,c,A,B.

    Integrals of the systems (37) and (42)–(44) solution were approximated using Simpson quadrature formulas. The behavior of magnitude B proved to be highly dependent on the behavior of the exponent at large, so the integral relation had to be replaced by the local condition (52), besides controlling the behavior of the TPPA near the maximum is more important than the weight of the exponent away from the wall. Thus, instead of the value of B, we include the value among the parameters sought, and the value of B is expressed from Eqs. (50) and (51).

    5. Results

    As an example of TPPA (see Section 3.2) used for matching of limiting asymptotics, consider the paper by Grasman et al. [33]. They dealt with Lyapunov exponents which characterize the dynamics of a system near its attractor. For the Van der Pol oscillator:

    x¨+μẋx21+x=0E57

    Similar to the asymptotic approximation of amplitude and period, expressions are derived for the nonzero Lyapunov exponent λ2for both small and large parameter μvalues:

    λ2=μ116μ3+26318432μ5+,μ0,E58
    λ2=3+4ln2232ln2μ+,μ.E59

    The overlap of these series does not take place. The authors of [33] remark: “Such an overlap comes within reach if in the regular expansion a large number of terms is included.” This is not correct, because the obtained series is asymptotic; so, with increasing of number of terms, the results will be worst. So, one needs a summation procedure. Some authors [34] proposed to use PAs, but in this case one needs hundreds of perturbation series terms. That is why we use TPPA. Using two terms from expansion (58) and one term from expansion (59), one obtains

    Lλ2μ=1+0.14μ21+0.079μ2E60

    Expression (60) has a pole at μ=12.66. Below, one can see some numerical results.

    In Table 1, the second column is made by calculation results by formula (4), the third column is made by paper data [33]. One can see that TPPA gives good result for any value of used parameter.

    μL (4)L (NR)
    11.0571.0648
    51.5131.4724
    101.6851.6358
    251.7591.7398
    501.7681.7691

    Table 1.

    Comparison for L of numerical results (NR) from paper by [33] with TPPA formulate (60).

    In Section 4.4, the problem was solved for several variants of the Mach number and the heating temperature: M=5;10;15,Ts=3;5;7of the streamlined flat plate, with constant Prandtl number values σ=0.76, adiabatic index κ=1.4, and two values of dynamical viscosity index μ=Tn:n=1;0.76. When the first equation of the systems (43) and (44) is solved, it becomes independent of the second equation and can be compared with the known Blasius solution (see Section 3), which was used as a test when compared to our method [35, 36, 37, 38, 39, 40]. Thus, the value of the parameters according to the exact solution is equal: a2=0.664; c=1.72. Our decision gives a2=0.6641;c=1.7308. Of course, such a good match is due to the fact that these parameters are largely determined by local internal asymptotics, more precisely, derived from the function on the wall. But also within the transition area, the deviation from the exact solution does not exceed 1÷2%(for φand T, respectively). Design values of parameters for determining approximations (37) and (38) for n = 1 are given in Table 2.

    MTsςmT1
    530.4261.340
    1030.7449.127
    1050.6377.929
    1070.5316.676
    1530.80622.00
    1550.75620.82
    1570.71219.75

    Table 2.

    TPPAs parameters for different Mach numbers M, temperature TS, and n = 1 values.

    If n=0.76, this value corresponds to the physical characteristics of the air, and the constant calculation results for the approximation formulas (49) and (50) are shown in Table 3.

    M51010
    Ts335
    a20.170.740.56
    ςm0.7290.7140.798
    c1.451.421.38
    T10.568.698.16

    Table 3.

    TPPAs parameters for different Mach numbers M, temperature TS, and n = 0.76 values.

    6. Conclusion

    The procedure of constructing the PA is much less labor-intensive than the construction of higher approximations of perturbation theory. PA can be applied to power series but also to the series of orthogonal polynomials. PA is locally the best rational approximation of a given power series. They are constructed directly and allow for efficient analytic continuation of the series outside its circle of convergence, and their poles in a certain sense localize the singular points (including the poles and their multiplicities) of the function at the corresponding region of convergence and on its boundary. PA is fundamentally different from rational approximations with (fully or partially) fixed poles, including the polynomial approximation, when all the poles are fixed in infinity. That is the above property of PA—effectively solving the problem of analytic continuation of power series—lies at the basis of their many successful applications in the analysis and the study of applied problems. Currently, the PA method is one of the most promising nonlinear methods of summation of power series and the localization of its singular points. Including the reason why the theory of the PA turned into a completely independent section of approximation theory, and these approximations have found a variety of applications both directly in the theory of rational approximations, and in perturbation theory.

    Thus, the main advantages of PA compared with the Taylor series are as follows:

    1. Typically, the rate of convergence of rational approximations greatly exceeds the rate of convergence of polynomial approximation. For example, the function eεin the circle of convergence approximated by rational polynomials Pn(ε)/Qn(ε) in 4п times better than an algebraic polynomial of degree 2n. More tangible, it is property for functions of limited smoothness. Thus, the function |ε| on the interval [−1,1] cannot be approximated by algebraic polynomials so that the order of approximation was better than 1/n, where n is the degree of polynomial. PA gives the rate of convergence ∼exp2n.

    2. Typically, the radius of convergence of rational approximation is large compared with the power series. Thus, for the function arctanx, Taylor polynomials converge only if ε1, and PA is everywhere in С\((- i∞, - i] [i, i∞)).

    3. PA can establish the position of singularities of the function.

    TPPA allows to overcome the locality of asymptotic expansions, using only a few terms of asymptotics. Unfortunately, the situations when both asymptotic limits have the form of power expansions are rarely encountered in practice, so we have to resort to other methods of AEFs construction, for example, the method quasirational approximation which is described in [23]. The method of combination (combining method) of asymptotics by using TPPA is alternative to the well-known matching method [6]; it is useful in local domains of transition layers where asymptotics are not uniform. This method was tested on well-known problems of mathematical physics, in particular, problems of fluid dynamics. The main advantage of the method is that it has an analytic form.

    How to cite and reference

    Link to this chapter Copy to clipboard

    Cite this chapter Copy to clipboard

    Igor Andrianov and Anatoly Shatrov (September 4th 2020). Padé Approximation to Solve the Problems of Aerodynamics and Heat Transfer in the Boundary Layer [Online First], IntechOpen, DOI: 10.5772/intechopen.93084. Available from:

    chapter statistics

    20total chapter downloads

    More statistics for editors and authors

    Login to your personal dashboard for more detailed statistics on your publications.

    Access personal reporting

    We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

    More About Us