Open access peer-reviewed chapter

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

Written By

Igor Andrianov and Anatoly Shatrov

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

DOI: 10.5772/intechopen.93084

From the Edited Volume

## Mathematical Theorems

Edited by Lyudmila Alexeyeva

Chapter metrics overview

View Full Metrics

## 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
• 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 fin 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.

“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. 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. 2. PAs gives meromorphic continuation of a given power series.

3. 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

1. 4. 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).

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

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

1. If c+df00, then

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

1. provided that there is fn/nε.

2. 6. 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

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

2. 7. 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 Rconfirms 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. 1. Functions F0(s) and F(s) are asymptotically equivalent at s→ ∞ and s→ 0, that is,

F0sFsats0ands.

1. 2. 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 Ais 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.

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

kcoefficients 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 aand bare 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, Ucan 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 xand yare the Cartesian coordinates such that the axis xis 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

Mis the Mach number, σis the Prandtl number,κis the adiabatic index, ψis the stream function, Tis the temperature, μis the viscosity coefficient, and xand yare 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 Mnumbers, 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 Ais 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 Bproved 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 Bis 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¨+μ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 nis 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.

## References

1. 1. Andrianov IV. Application of Padé approximants in perturbation methods. Advances in Mechancis. 1991;14(2):3-15
2. 2. Barantsev RG, Pashkevich DA, Shatrov AV. Heat transfer in a boundary layer of reactive gas. In: Heat/Mass Transfer MIF-2000 [in Russian]. Minsk; 2000. pp. 185-188
3. 3. Frost PA, Harper EY. Extended Padé procedure for constructing global approximations from asymptotic expansions: An explication with examples. SIAM Review. 1976;18:62-91. DOI: 10.1137/1018003
4. 4. Shatrov AV. Matching of interior and exterior asymptotics in boundary-value problems of mathematical physics. Journal of Mathematical Sciences. 2018;230(5):804-807. DOI: 10.1007/s10958-018-3794-0
5. 5. Andrianov IV, Awrejcewicz J. New trends in asymptotic approaches: Summation and interpolation methods. Applied Mechanics Reviews. 2001;54:69-92. DOI: 10.1115/1.3097289
6. 6. Brezinski C. Convergence acceleration during the 20th century. Journal of Computational and Applied Mathematics. 2000;122:1-21. DOI: 10.1016/S0377-0427(00)00360-5
7. 7. Bellman R. Perturbation Techniques in Mathematics, Physics and Engineering. London: Holt, Rinehart and Winston; 1964. 118 p
8. 8. Hinch EJ. Perturbation Methods. London: Cambridge University Press; 1991. 160 p. DOI: 10.1017/S0022112092220644
9. 9. Kublanovskaya VN. Application of analytic continuation by means of change of variables in numerical analysis. Trudy Matematicheskogo Instituta im. Steklov. 1959;53:145-185 (in Russian)
10. 10. Van Dyke M. Analysis and improvement of perturbation series. The Quarterly Journal of Mechanics and Applied Mathematics. 1974;27:423-450http://dx.doi.org/10.1093/qjmam/27.4.423
11. 11. Van Dyke M. Computer extension of perturbation series in fluid mechanics. SIAM Journal on Applied Mathematics. 1975;28:720-734. DOI: 10.1007/978-1-4757-9259-1_2
12. 12. Van Dyke M. Perturbation Methods in Fluid Mechanics. Stanford, CA: The Parabolic Press; 1975; 271 p. DOI: 10.1017/S000192400003390X
13. 13. Andrianov IV, Mikhlin Yu V, Tokarzhewsky S. Two-point Padé approximants and their applications to solving mechanical problems. Journal of Theoretical and Applied Mechanics. 1997;35(3):577-606
14. 14. Barantsev RG. Definition of asymptotics and system triades. In: Asymptotic Methods in the Theory of Systems [in Russian]. Irkutsk; 1980. pp. 70-81
15. 15. Barantsev RG, Pashkevich DA. Matching of asymptotics in transition layers. In: Asymptotic Methods in Problems of Aerodynamics and Design of Aircrafts [in Russian]. Irkutsk; 1994. pp. 67-70
16. 16. Gonchar AA. Some recent results in diagonal Padé approximants. In: Approximation Theory V. Boston: Academic Press; 1986. pp. 55-70
17. 17. Litvinov GL. Approximate construction of rational approximations and the effect of autocorrection error. Russian Journal of Mathematical Physics. 1994;1:313-352.arXiv:math/0101042v1
18. 18. Litvinov GL. Error autocorrection in rational approximation and interval estimation. Central European Journal of Mathematics. 2003;1:1-33. DOI: 10.2478/BF02475663
19. 19. Suetin SP. Padé approximants and efficient analytic continuation of a power series. Russian Mathematical Surveys. 2002;57:43-141. DOI: 10.4213/rm475
20. 20. Suetin SP. On interpolation properties of diagonal Padé approximants of elliptic functions. Russian Mathematical Surveys. 2004;59:201-202. DOI: 10.4213/rm767
21. 21. Alekseeva EA, Barantsev RG, Shatrov AV. Matching of thermal asymptotics in a boundary layer. Vestnik St. Petersburg University Ser. 1. 1996;8:96-99. ISSN: 00240850
23. 23. Shatrov AV. Combination of asymptotics in boundary-value problems of hydrodynamics. In: Proceedings of the International Conference (OFEA-2001); June 25–29, 2001; St. Petersburg. 2001. pp. 59-60
24. 24. Shatrov AV. Matching of interior and exterior asymptotics in transition layers of viscous liquids and gases. In: Proceedings of the VIII All-Russian Conference on Theoretical and Applied Mechanics; Perm, August 23–29, 2001; Ural Branch RAS, Ekaterinburg. 2001. pp. 602-603 (in Russian)
25. 25. Baker GA Jr. and Graves-Morris P. Padé Approximants. New York: Cambridge University Press, 1996. 673 p. DOI: 10.1017/CBO9780511530074
26. 26. Slepyan LI, Yakovlev Yu S. Integral Transforms in Non-Stationary Problems of Mechanics. Leningrad: Sudostroyenie; 1980; 343 p. (in Russian)
27. 27. Vyatchin AV. On convergence of Padé approximants. Moscow University Mathematics Bulletin. 1982;37(4):1-4
28. 28. Kashin AP, Kvaratskheliya TM, Maksimov MZ, Chikovani ZE.Higher approximations of the reduced method of matched asymptotic expansions and its convergence. Theoretical and Mathematical Physics. 1989;78:278-284. DOI: 10.1007/BF01017665
29. 29. Chalbaud E, Martin P. Two-point quasifractional approximant in physics: method improvement and application to Jn(x). Journal of Mathematical Physics. 1992;33:2483-2486. DOI: 10.1063/1.529618
30. 30. Gaunt DS, Guttman AJ. Series expansions: analysis of coefficients. In: Domb C, Green MS, editors. Phase Transitions and Critical Phenomena. Vol. 3. New York: Academic Press; 1974. pp. 181-243
31. 31. Kruskal MD. Asymptotology. In: Drobot S, Viebrock PA, editors. Proceedings of the Conference on Mathematical Models in Physical Sciences; April 1962; Notre Dame, Indiana. Englewood Cliffs, NJ: Prentice-Hall. 1963. pp. 17-48
32. 32. Martin P, Baker GA Jr. Two-point quasifractional approximant in physics. Truncation error. Journal of Mathematical Physics. 1991;32:313-328. DOI: 10.1063/1.529304
33. 33. Grasman J, Verhulst F, Shagi-Di S. The Lyapunov exponents of the Van der Pol oscillator. Mathematical Methods in the Applied Sciences. 2005;28:1131-1139. DOI: 10.1002/mma.606
34. 34. Andersen CM, Dadfar MB, Geer JF. Perturbation analysis of the limit cycle of the van der Pol equation. SIAM Journal on Applied Mathematics. 1984;44:881-895
35. 35. Bender SM, Orszag SA. Advanced Mathematical Methods for Scientists and Engineers. New York: McGraw-Hill; 1978. DOI: 10.1007/978-1-4757-3069-2
36. 36. Buonomo A. The periodic solution of Van der Pol’s equation. SIAM Journal of Applied Mathematics. 1999;59:156-171. DOI: 10.4236/jamp.2019.71001
37. 37. Miller PD. Applied Asymptotic Analysis. Providence, Rhode Island: AMS; 2006; 488 p. Available from:https://bookstore.ams.org/gsm-75/5
38. 38. Shanks D. Nonlinear transforms of divergent and slowly convergent sequences. Journal of Mathematical Physics. 1955;34:1-42
39. 39. Vinogradov VN, Gay EV, Rabotnov NC. Analytical Approximation of Data in Nuclear and Neutron Physics. Moscow: Energoatomizdat; 1987; 128 p. (in Russian)
40. 40. Weniger EI. Nonlinear sequence transformations for the acceleration of convergence and the summation of divergent series. Computer Physics Report. 1989;10:189-371.arXiv:math/0306302v1

Written By

Igor Andrianov and Anatoly Shatrov

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