Open access peer-reviewed chapter

Perturbation Expansion to the Solution of Differential Equations

By Jugal Mohapatra

Submitted: July 21st 2020Reviewed: September 24th 2020Published: October 29th 2020

DOI: 10.5772/intechopen.94173

Downloaded: 188


The main purpose of this chapter is to describe the application of perturbation expansion techniques to the solution of differential equations. Approximate expressions are generated in the form of asymptotic series. These may not and often do not converge but in a truncated form of only two or three terms, provide a useful approximation to the original problem. These analytical techniques provide an alternative to the direct computer solution. Before attempting to solve these problems numerically, one should have an awareness of the perturbation approach.


  • perturbation methods
  • asymptotic expansion
  • boundary layer
  • principle of least degeneracy
  • inner and outer expansion

1. Introduction

The governing equations of physical, biological and economical models often involve features which make it impossible to obtain their exact solution. For instance, problems where we observe “a complicated algebraic equations”, “the occurrence of a complicated integral”, in case of differential equations (DE), “a varying coefficients or nonlinear term” sometimes problems with an awkwardly shaped boundary are tough to solve with the limited methods for finding analytical solutions. The main purpose of this chapter is to describe the application of perturbation expansion techniques to the solution of DE. Approximate expressions are generated in the form of asymptotic series. These may not and often do not converge but in a truncated form of only two or three terms, provide a useful approximation to the original problem. These analytical techniques provide an alternative to the direct computer solution. Before attempting to solve these DE numerically, one should have an awareness of the perturbation approach. An example of this occurs in boundary layer problems where there are regions of rapid change of quantities such as fluid velocity, temperature or concentration. Appropriate scaling of the boundary layer dimension is required before a numerical solution can be generated which will capture the behavior in the rapidly changing region.

When a large or small parameter occurs in a mathematical model of a process there are various methods of constructing perturbation expansions for the solution of the governing equations. Often the terms in the perturbation expansions are governed by simpler equations for which the exact solution techniques are available. Even if exact solutions cannot be obtained, the numerical methods used to solve the perturbation equations approximately are often easier to construct than the numerical approximation for the original governing equation.

First, we consider a model problem for which an exact solution is available against which the perturbation expansion can be compared. A feature of the perturbation expansions is that they often form divergence series. The concept of an asymptotic expansion will be introduced and the value of a truncated divergent series will be demonstrated.


2. Projectile motion

This example studies the effect of small damping on the motion of a particle. Consider a particle of mass Mwhich is projected vertically upward with an initial speed U0. Let Udenote the speed at some general time T. If air resistance is neglected then the only force acting on the particle is gravity, Mg, where gis the acceleration due to gravity and the minus sign occurs because the upward direction is chosen to be the positive direction. Newton’s second law governs the motion of the projectile, i.e.,


Integrating (1), we obtain the solution U=CgT. The constant of integration is determined from the initial condition U0=U0, so that


On defining the non-dimensional velocity v, and time t, by v=U/U0and t=gT/U0, the governing equation becomes


with the solution vt=1t.

Taking account of the air resistance, and is included in the Newton’s second law as a force dependent on the velocity in a linear way, we obtain the following linear equation


where the drag constant Kis the dimensions of masa/time. In the non-dimensional variables, it becomes


Let us denote the dimensionless drag constant by ε, then the governing equation is


where ε>0is a “small” parameter and the disturbances are very “small”. The damping constant Kin (4) is small, since Khas the dimensions of mass/time and a small quantity in units of kilograms per second.

2.1 Perturbation expansion

It is possible to solve (6) exactly since it is of variables separable form. Here, we solve by an iterative process, known as perturbation expansion for the solution.

Let videnotes the ith iterate, which is obtained from the equation


The justification for this iterative scheme is that the term εvinvolves the small multiplying coefficient ε, and so the term itself may be expected to be small. Thus, the term εviwhich should appear on the RHS of (7) to make it exact, may be replaced by εvi1with an error which is expected to be small.

The first iterate is obtained by neglecting the perturbation, thus


This is known as the unperturbed problem, and direct integration yields


The next iterate v1, satisfies


and integration yields


Similarly, v2satisfies


Direct integration yields the solution


Rearranging the terms in these iterates in ascending powers of ε, we obtain


Clearly as the iteration proceeds the expressions are refined by terms which involve increasing powers of ε. These terms become progressively smaller since εis a small parameter. This is an example of a perturbation expansion. It will often be the case that perturbation expansions involve ascending integer powers of the small parameter, i.e., ε0ε1ε2. Such a sequence is called an asymptotic sequence. Although this is the most common sequence which we shall meet, it is by no means unique. Examples of other asymptotic sequences are ε1/2εε3/2ε2and ε0ε2ε4. In each case the essential feature is that subsequent terms tend to zero faster than previous terms as ε0.

An alternative procedure to that of developing the expansion by iteration is to assume the form of the expansion at the outset. Thus, if we assume that the perturbation expansion involves the standard asymptotic sequence ε0ε1ε2, then the solution v, which depends on the variable t, and the parameter ε, is expressed in the form


The coefficients v0t,v1t,of powers of εare functions of tonly. Substituting expansion (9) in the governing Eq. (6) yields the following


Thus, the coefficients of powers of εcan be equated on the left– and right–hand sides of (10):


The proof of validity of this fundamental procedure can be developed by first setting ε=0in (10) which yields the first equation of (11). This result allows the first member of the left– and right–hand side of Eq. (10) to be removed. Then, after dividing the remaining terms by εwe obtain the equation


This is valid for all nonzero values of εso that on taking the limit as ε0we obtain the second equation of (11). Repeating the procedure, we obtain the other equations.

Integrating the equations in (11), we obtain


Using these values in (9), we obtain that


This is the same as the expansion (8) which is generated by iteration.

The IVP (6) can be solved exactly as


The perturbation expansion can be obtained from (12) by replacing the exponential function by its Maclaurin expansion, i.e.,


This is the same as the expansion (12). Thus, the perturbation expansion approach is justified in this case. One can refer the books [1, 2].

3. Asymptotics

The letters Oand oare order symbols. They are used to describe the rate at which functions approach limit values. We will consider the types of limit values, namely zero, a finite number but nonzero and infinite.

If a function fxapproaches a limiting value at the same rate of another function gxas xx0, then we write


The functions are said to be of the same order as xx0. The test for this is the limit of the ratio. Thus, if limxx0fxgx=C,where Cis finite, then we say (15) holds.

For example, we have the following functions:


The expression


means that limxx0fxgx=0.This is a stronger assertion that the corresponding O–formula. The relation (16) implies the relation (15), as convergence implies boundedness from a certain point onwards.

We have the following functions satisfy the o–relation:


3.1 Asymptotic expansions

Consider the expansion


is an asymptotic expansion as x, if, for any N,


The following expansion is used when (17) and (18) hold,


Here, limnRN=0, for any value of N.

The sequence 11/x1/x2is an asymptotic sequenceas x. The characteristic feature of such sequences is that each member is dominated by the previous member. In constructing examples it is easier to deal with the limit zero than any other. Thus, for the case x, we let ε=1/x, which for xx0, we let ε=xx0so that without loss of generality we may confirm our attention to the limit ε0. The standard asymptotic sequence is 1εε2as ε0. If we let δnεrepresent members of an asymptotic sequence δ0εδ1εas ε0, then the following condition must hold


Some examples of asymptotic sequences are

  1. 1sinεsinε2sinε3, here we have


  • 1ln1+εln1+ε2ln1+ε3, with δ0=1,δn=ln1+εnn1, we have

  • limε0δ1δ0=limε0ln1+ε=0,limε0δn+1δn=limε0ln1+εn+1ln1+εn=limε0εn+1+Oε2n+2εn+Oε2n=0.

    The general expression for an asymptotic expansion of a function fε, in terms of an asymptotic sequence δnεis


    where the coefficients anare independent of ε. The expression (20) involving the symbol , means that for all N,




    If a function possesses an asymptotic expansion involving the sequence δ0εδ1εthen the coefficients anof the expansion (21) given by the expression (24) are unique. However, another function may share the same set of coefficients. Thus, while functions have unique expansions, an expansion does not correspond to a unique function.

    Consider a function fxε, which depends on both an independent variable x, and a small parameter ε. Suppose that fxεis expanded using an asymptotic sequence δnε,


    The coefficients of the gauge functions δnεare functions of x, and the remainder after Nterms is a function of both xand ε. For this to be an asymptotic expansion, we require


    Refer [3, 4] for more details. For (24) to be a uniform asymptotic expansion the ultimate proportionality between RNand δN+1must be bounded by a number independent of x, i.e.,


    for εin the neighborhood near zero, where Kis a fixed constant.

    An example of a uniform asymptotic expansion is fxε=11εsinx.

    An example of a nonuniform expansion is


    Here, one cannot find a fixed Kwhich satisfy RNKεN+1,because for any choice of K, xcan be chosen so that xN+1exceeds this value.

    3.2 Nonuniformity

    The expansion (27) becomes nonuniform when subsequent terms are no longer small corrections to previous terms. This occurs when subsequent terms are of the same order or of dominant order than previous terms. Subsequent terms dominate previous terms for larger x, for example, when x=O1/ε2. The expansion is valid for x=O1since then subsequent terms decrease by a factor of ε. The expansion remains valid for large x, provided xis not as large as 1/ε. For instance, the expansion is valid for x=O1/ε, as ε0.

    The critical case is such that subsequent terms are of the same order. This determines the region of nonuniformity. In (27), the region of nonuniformity occurs when εx=O1, i.e., x=Oε1, as ε0.

    3.2.1 Sources of nonuniformity

    There are two common reasons for nonuniformities in asymptotic expansions, they are

    1. Infinite domains which allow long-term effects of small perturbations to accumulate.

    2. Singularities in governing equations which lead to localized regions of rapid change.

    Consider the nonlinear Duffing equation


    Suppose the solution may be expanded using the standard asymptotic sequence


    On substituting this in (28) and in the initial conditions, we get


    Equating like of powers of εon both sides, we get




    Solving Eqs. (30) and (31), we obtain


    The term tsintin the expansion (32) is called a secular term. It is an oscillating term of growing amplitude. All other terms are oscillating of fixed amplitude. The secular term leads to a nonuniformity for large t. The region of nonuniformity is obtained by equating the order of the first and second terms,


    The trigonometric functions are treated as O1terms. Thus, the region of nonuniformity is t=O1/ε, as ε0.

    The second common source of nonuniformities is associated with the presence of singularities. Consider, the following initial-value problem:


    where ε>0is a small parameter. Suppose yhas the expansion


    Substituting (34) in (33), we have


    Equating coefficients of like powers of εon both sides, we get


    Clearly, y0cannot satisfy the boundary condition y00=2as no constant of integration is available because the equation determining y0is an algebraic equation not a differential equation, and no additional conditions are required. Thus, we have obtain the expression


    but the initial condition y0=2has not been satisfied.

    The unperturbed problem, obtained by setting ε=0is not a DE, but an algebraic equation y=ex. This cannot satisfy an arbitrarily imposed condition at x=0. For any nonzero value of ε, (33) becomes a first-order DE which can satisfy an initial condition. This is an example of a singular perturbation problem (SPP), where the behavior of the perturbed problem is very different from that of the unperturbed problem.

    Thus, the perturbation expansion (36) is a good approximation of the exact solution away from the region x=0. To see this, let us compare (36) with the following exact solution:


    The perturbation expansion (36) generates the second member of (37), but not the first member. The coefficient ex/ε)is a rapidly varying function which takes the value of unity at x=0, and rapidly decays to zero for x>0. Clearly, y0provides a good approximation away from the region x=0. The region near x=0is called the boundary layer. These regions usually occur when the highest order derivative of a DE is multiplied by a small parameter. The unperturbed problem, obtained by setting ε=0is of lower order and consequently cannot satisfy all the boundary conditions. This leads to boundary layer regions where the solution varies rapidly in order to satisfy the boundary condition.

    Boundary layers are regions of nonuniformity in perturbation expansions of the form (36).

    4. Boundary layer

    Boundary layers are regions in which a rapid change occurs in the value of a variable. Some physical examples include “the fluid velocity near a solid wall”, “the velocity at the edge of a jet of fluid”, “the temperature of a fluid near a solid wall.” Ludwig Prandtl pioneered the subject of boundary layer theory in his explanation of how a quantity as small as the viscosity of common fluids such as water and air could nevertheless play a crucial role in determining their flow. The viscosity of many fluids is very small and yet taking account of this small quantity is vital. The essential point is that the viscous term involves higher order derivatives so that its omission necessitates the loss of a boundary condition. The ideal flow solution allow slip to occur between a solid and fluid. In reality the tangential velocity of a fluid relative to a solid is zero. The fluid is brought to rest by the action of a tangential stress resulting from the viscous force.

    Mathematically the occurrence of boundary layers is associated with the presence of a small parameter multiplying the highest derivative in the governing equation of a process. A straightforward perturbation expansion using an asymptotic sequence in the small parameter leads to differential equations of lower order than the original governing equation. In consequence not all of the boundary and initial conditions can be satisfied by the perturbation expansion. This is an example of what is commonly referred to as a singular perturbation problem. The technique for overcoming the difficulty is to combine the straightforward expansion, which is valid away from the layer adjacent to the boundary. The straightforward expansion is referred to as the outer expansion. The inner expansionassociated with the boundary layer region is expressed in terms of a stretched variable, rather than the original independent variable, which takes due account of the scale of certain derivative terms. The inner and outer expansions are matched over a region located at the edge of the boundary layer. The technique is called the method of matched asymptotic expansions.

    Consider the following two-point boundary value problem:


    where ε>0is a small parameter. If we assume that upossesses a straightforward expansion in powers of ε,


    then the equations associated with powers of εleads to


    and the boundary conditions require


    which leads to


    Equation (42) require that each unxsatisfy two boundary conditions. This is in general impossible since Eqs. (41) and (42) governing each unare of first-order. Now the question is which one of the boundary condition has to be taken into account. We will find out that the boundary condition at x=0must be abandoned and consequently the expansion (39) is invalid near x=0.

    The general solution of (42) is u0x=x2+x+C, using the boundary condition u01=4, we obtain


    From (42), we obtain the equations


    and its solutions are


    Therefore, the outer expansion is


    where ‘out’ label is used to indicate that the solution is valid away from the region near x=0. Clearly uoutfails to satisfy the boundary condition at x=0. The reason why the outer solution is of use is that it closely follows the exact solution of the problem except in a narrow region near x=0, where the exact solution changes rapidly in order to satisfy the boundary condition.

    The exact solution of the BVP (38) can be obtained as


    The constants Aand Bare determined from the boundary conditions:


    We know that e1/ε=oεN, as ε0, for all N. This means that the exponential term tends to zero faster than any power of ε, as ε0. It is called a transcendentally small term(T.S.T.) and can always be neglected since its contribution is asymptotically always less than any power of ε. Thus, (45) gives


    and the exact solution is


    after rearranging the terms in asymptotic order, we obtain


    Comparing the exact solution with the outer expansion shows that the terms involving ex/εare absent. The effect of these terms is negligible when x=O1. But, when x=Oε, then e=O1. It is clear that as ε0the region in which the outer solution departs from the exact solution becomes arbitrarily close to x=0with a thickness Oε. This region is called the boundary layer.

    The behavior of the exact solution and the zeroth-order term of the outer expansion are plotted in Figure 1 for various values of ε.

    Figure 1.

    Exact solution of(38)for various values ofε.

    By differentiating the leading order term u0ex, of the exact solution, we have


    Outside the boundary layer, i.e., for x=O1, we have ex/ε=oεN,N, so ε1ex/εand ε2ex/εare also transcendentally small. Within the boundary layer when x=Oε, we have ex/ε=O1. The order of u0exand its derivatives are given below:


    This indicates that xis the appropriate independent variable outside the boundary layer where u0exand its derivatives are of O1quantities. However, within the boundary layer the appropriately scaled independent variable is s=x/ε, then


    so that within the boundary layer


    The variable s=x/εis called a stretched variable. The differential equations becomes


    We assume a boundary layer expansion, called the inner expansionof the form


    The inner expansion will satisfy the boundary condition at x=s=0namely v0s=0=1giving v00=1, and vn0=0,n=1,2,. Substituting (49) into the DE (48), we obtain the following set of equations:


    with solutions


    The boundary condition at x=1cannot be used to determine the constants appearing in these solutions because the DEs (50) are only valid in the boundary layer. The constants in (51) are determined by matching the inner and outer expansions. We shall first restrict our attention to matching the leading order expansions u0and v0. The method which we shall apply is Prandtl’s matching condition.

    The leading order terms in the ‘inner’ and ‘outer’ expansions are to be matched at the ‘edge of the boundary layer’. Of course there is no precise edge of the boundary layer, we simply know that it has thickness of order Oε. A plausible matching procedure would be to equate u0and v0at a value of xsuch that the region of rapid change has passed. We might choose to equate the terms at the point x=5ε. The leading order expansions are


    Equating at x=5εgives the following:


    If, instead we choose to match at x=6ε, then we obtain


    These two expressions differ in the argument of the exponential and differ algebraically with 5εreplaced by 6ε. The exponential functions are approaching transcendentally small values so that their contribution can be neglected. The algebraic difference is of Oε. Thus, the arbitrariness in the decision of the point at which we choose to equate the expansions leads to a difference of Oε. But we are only dealing with leading order expansions anyway. The difference between the exact solution and the leading order expansions will of Oεso that an arbitrariness in v0and u0of Oεis immaterial. Rather than choose between, for example, 5εand 6εas the value of xto evaluate u0we may take the value at x=0, since


    where the remainder is uniformly Oεsince the gradient of u0is O1. For the inner expansion we are to ensure that the rapidly varying function has achieved its asymptotic value at the edge of the boundary layer. This means that the term ex/εshould be replaced by zero. This can be achieved by taking the limit s. Thus, rather than choosing a specific point to equate the inner and outer terms er are led to the following Prandtl’s matching condition:


    The limit smay appear rather dangerous since although it certainly removes the exponential term it could lead to an algebraically unbounded term. For example, if v0=As+1Aes, then the first member would be unbounded as s. This possibility can be eliminated since the inner expansion must be of a form which varies rapidly for x=Oεbut not for x=O1, i.e., not for s. In practice, if the boundary layer has been properly located and the correct inner variable is used then Prandtl’s matching condition is valid and elegantly avoids the need to choose an arbitrary ‘edge’ of the boundary layer.

    Applying these conditions to the current example leads to


    which yields A=2. Thus the leading order terms in the expansion solutions are


    To prove that these are valid leading terms we consider uex:


    We conclude that the matching condition has correctly predicted the leading order terms.

    4.1 Composite expansion

    As single composite expression for these leading order terms can be constructed using the combination


    where u0matchis given by (52). Then,


    For the current example, u0match=2, so the composite expansion is


    Prandtl’s matching condition can only be used for the leading order terms in the asymptotic expansions.

    The outer, inner and composite expansions of the BVP (38) are presented in Figures 2 and 3 for different values of ε. From these figures, one can easily identify the need and efficiency of the composite expansion.

    Figure 2.

    Outer, inner and composite expansions. (a) For ε = 0.2; (b) For ε = 0.1.

    Figure 3.

    Outer, inner and composite expansions. (a) For ε = 0.05; (b) For ε = 0.025.

    4.2 Boundary layer location

    Consider the following linear DE


    The following general statements can be made about the boundary layer location and the nature of the inner expansion.

    Case I.If ax>0throughout x1x2, then the boundary layer will occur at x=x1. The stretching transformation will be s=xx1/ε, and the one-term inner expansion will satisfy


    The solution of this equation is


    where A+B=ux=x1. The other condition to determine the constants Aand Bis obtained by matching with the value of the outer expansion at x=x1.

    Case II.If ax<0throughout x1x2, then the boundary layer will occur at x=x2. The stretching transformation will be s=x2x/ε, and the one-term inner expansion will involve the rapidly decaying function eax2x2x/ε.

    Case III.If axchanges sign in the interval x1<x<x2, then a boundary layer occurs at an interior point x0, where ax0=0and boundary layers may also occur at both ends x1and x2.

    4.3 Boundary layer thickness and the principle of least degeneracy

    The boundary layers which we have met so far have all had thickness Oε. By this we mean that a variation of Oεin the independent variable will encompass the region of rapid change in the dependent variable. The associated stretched independent variable s, appropriate for the boundary layer is related to xby a linear transformation involving division by ε.

    There are practical situations where the boundary layer thickness will be of Oεp. This means that if the boundary layer is located at x=x0, then the appropriate stretching transformation is s=xx0/εp. More generally, the choice of the function δεto use in the stretching transformation s=xx0/δεis determined by the need to represent the region of rapid change correctly. We must ensure that the boundary layer solution contains rapidly varying functions. The form of the governing equation in the boundary layer region must have sufficient structure to allow such solutions.

    Consider the example


    Since the signs of the first and second derivatives are the same, and the boundary layer will occur at x=0. We are not going to assume at the outset that the boundary layer thickness is Oε. Our intension is to deduce that the appropriate stretching variable is s=x/ε.

    The one-term outer expansion u0satisfies du0dx+u0=x,u01=2The solution is


    To determine the inner expansion we first wrongly assume that the boundary layer thickness is Oε1/2. The stretching transformation s=x/ε1/2changes the original DE (56) into the following one:


    If the appropriate stretching transformation has been used for the boundary layer then dv/dsand d2v/ds2will be of O1within it. The leading order expansion v0will satisfy the dominant part of (58), i.e., the component of Oε1/2


    The solution is v0s=1. This of course does not have the rapidly varying behavior which we anticipate in the boundary layer. Prandtl’s matching condition cannot be satisfied since


    Thus, we reject the assumption of a boundary layer of thickness Oε1/2.

    Next, suppose that the boundary layer thickness is Oε2and again we will discover that this is incorrect because the corresponding inner expansion cannot be matched to the outer expansion. Proceeding with the analysis we introduce the stretching transformation s=x/ε2which leads to the equation


    Again we argue that if the appropriate stretching has been used then all derivatives are of O1so that the governing equation for the leading term in Oε3, namely


    The solution is v0s=1+As, where the constant Ais to be determined from matching. This solution is rapidly varying but the rapidity does not decay at the edge of the boundary layer (i.e., as s). Indeed, we cannot match v0to the outer expansion because the term Asbecomes arbitrarily large as s.

    The correct choice of stretching transformation is s=x/εshowing that the boundary layer thickness is Oε. The boundary layer equation becomes


    The dominant equation satisfied by v0is O1/ε, namely


    The solution is v0s=1A+Aes. The last member provides the necessary rapid decay away from the point x=s=0. Prandtl’s matching condition requirest


    which leads to A=22e, and


    The one-term composite expansion is


    The leading order boundary layer equation associated with the stretching transformation s=x/ε, (61) involves more terms than (59), associated with s=x/ε1/2, and (60) associated with s=x/ε2. The extra term in (61) allows sufficient structure in the solution to produce the required boundary layer behavior. An aid for choosing the boundary layer thickness is to seek a stretching transformation which retains the largest number of terms in the dominant equation governing v0. This referred to as the principle of least degeneracyby Van Dyke.

    The composite expansion (62) can be verified by comparing with the exact solution of (56). The general solution of (56) is




    We expand 14εusing the binomial series, 14ε=12ε+Oε2, then


    so that


    Using the boundary conditions and by neglecting the transcendentally small term e1/ε, we have C1=2e,C2=21e. Then, (63) becomes


    There is an apparent discrepancy between (64) and the composite expansion (62) in the coefficient of the ex/εterm. There is an extra term only contributes in the boundary layer where x=Oεso that the coefficient exmay to leading order, be replaced by unity. Thus, the leading order composite expansion and the leading order term in the exact solution are in complete agreement.

    4.4 Boundary layer of thickness of Oε

    Consider the following two-point BVP:


    We seek a one-term composite expansion for the above BVP. We will tentatively assume that a boundary layer occurs at x=0although the vanishing of the coefficient of the first derivative suggests the possibility of nonstandard behavior.

    The one term outer expansion satisfies


    Its exact solution is u0x=2e11/x.

    Let us assume that the boundary layer thickness is of Oεp, where pis to be determined from the principle of least degeneracy. The stretched variable is s=x/εp, and (65) becomes


    The second-term is always dominated by the third, so the principle of degeneracy requires the first term to be of the same order as the third term (i.e., O1). Thus, p=1/2, and the one-term inner expansion satisfies


    The solution of the above problem is v0s=Aes+1Aes. Prandtl’s matching condition requires


    which yields A=0. This example is rather special in that Awill be zero for all boundary conditions.

    The on-term composite expansion is


    We conclude this example with the observation that a choice for the value of the index pother than p=1/2leads to boundary layer equations with insufficient structure to generate the required rapidly decaying behavior.

    Thus, if p>1/2, the dominant equation becomes


    which gives v0s=1+As. It is obvious that Prandtl’s matching condition cannot be used to determine A. Whereas, if p<1/2the dominant equation degenerates to v0s=0which does not satisfy the boundary condition at s=0.

    4.5 Interior layer

    Consider the BVP:


    The coefficient of the first derivative (convective term) is positive in 01which indicates the occurrence of a boundary layer at the left hand limit x=0. While the corresponding coefficient is negative in the range 1<x<0indicates a boundary layer located at the right-hand limit which is again is x=0. Thus, we are led to expr = ect two outer expansions for positive and negative xrespectively and an inner expansion in the interior layer located at x=0. We denote the leading term in the outer expansion for positive xby u0+, it satisfies


    with the solution u0+x=2ex.

    The outer expansion for negative x, u0satisfies


    with the solution u0x=ex.

    We suppose the boundary layer at x=0has thickness Oεpand determine the index pusing the principle of least degeneracy. Let s=x/εpso that the DE becomes


    The third term is dominated by the second term. The first term has the same order as the second term if p=1/2. For this choice of pthe leading term of the inner expansion v0satisfies


    Its solution can be given by


    Prandtl;s matching condition applied to the region x>0is


    and corresponding for x<0, we have


    Using the limiting values erf±=±1yields v00=1.5and B=0.5. The leading order terms over the whole region are


    A composite expansion cannot be formed in the standard way when there is more than one outer solution. However, the behavior of v0for x>Oεis as follows:


    Utilizing this enables a uniformly valid one-term composite expansion to be constructed which yields the correct coefficient of exoutside the boundary layer and the correct leading order behavior within the boundary layer. It is


    4.6 Nonlinear differential equation

    Consider the following semilinear


    The coefficient of the first and second order derivatives have the same sign, so the boundary layer will occur at the left boundary x=0. The one-term outer expansion satisfies


    and the solution is u0x=1/1+x. The stretching transformation for the inner region will be s=x/εand therefore, the inner expansion satisfies


    The one-term inner expansion v0satisfies the dominant part of this equation, i.e.,


    which gives v0s=A+2Aes. Prandtl’s matching condition yields A=1, and the composite one-term uniformly valid expansion is


    Next, consider the quasilinear problem


    The nonlinearity is associated with the first derivative term. The location of the boundary layer depends on the relative sign of the first and second derivative coefficients. If we assume that the dependent variable is nonnegative throughout the interval 0<x<1, then the boundary layer will occur at x=0. The one-term outer expansion satisfies


    with the solution u0x=2x+2.

    Assuming that the boundary layer thickness is Oε, therefore, the dominant-order equation for the one-term inner expansion becomes


    Its solution is v0s=atanhas. Prandtl’s matching condition yields a=2. Thus, v0s=2tanh2s, and the uniformly valid one-term composite expansion is


    Application of perturbation techniques to partial differential equations, and other types of problems can be seen in the books [5, 6].

    Conflict of interest

    The authors declare no conflict of interest.


    Notes/thanks/other declarations

    I owe a great debt to my mentor Prof. S Natesan, Department of Mathematics, IIT Guwahati, who introduced me to this topic. The chapter was discussed during my stay at IIT Guwahati.

    © 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

    How to cite and reference

    Link to this chapter Copy to clipboard

    Cite this chapter Copy to clipboard

    Jugal Mohapatra (October 29th 2020). Perturbation Expansion to the Solution of Differential Equations, A Collection of Papers on Chaos Theory and Its Applications, Paul Bracken and Dimo I. Uzunov, IntechOpen, DOI: 10.5772/intechopen.94173. Available from:

    chapter statistics

    188total chapter downloads

    More statistics for editors and authors

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

    Access personal reporting

    Related Content

    This Book

    A Collection of Papers on Chaos Theory and Its Applications

    Edited by Paul Bracken

    Next chapter

    Application of Perturbation Theory in Heat Flow Analysis

    By Neelam Gupta and Neel Kanth

    Related Book

    First chapter

    Introductory Chapter: Dynamical Symmetries and Quantum Chaos

    By Paul Bracken

    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