The paper considers the asymptotic solution of two-point boundary value problems εy” + A(x)y’ = 0, 0 ≤ x ≤ 1, when 0 < ε ≪ 1, A(x) is smooth with isolated zeros, y(0) = 0 and y(1) = 1. By using perturbation method, the limit asymptotic solutions of various cases are obtained. We provide a reliable and direct method for solving similar problems. The limiting solutions are constants in this paper, except in narrow boundary and interior layers of nonuniform convergence. These provide simple examples of boundary layer resonance.
A typical turning point problem consists of the linear differential equation
ε y ″ − x y ′ + n y = 0 , (1)
for a nonnegative integer n on − 1 ≤ x ≤ 1 with prescribed boundary values y ( ± 1 ) and a small positive parameter ε , i.e., 0 < ε ≪ 1 . Limiting solutions, away from narrow so-called boundary and interior shock layers of rapid change, take the form
Y 0 ( x ) = x n C (2)
as ε → 0 for constants C, so satisfy the limiting reduced equation
x Y ′ 0 = n Y 0 . (3)
Determining constants C and the location of layers is a nontrivial task, the subject of boundary layer resonance [
A simple, but still rich, the related problem concerns the asymptotic solution of the two-point problem [
ε y ″ + A ( x ) y ′ = 0 , 0 ≤ x ≤ 1 , (4)
with the special boundary values
y ( 0 ) = 0 and y ( 1 ) = 1 (5)
and a smooth coefficient A ( x ) . Its unique exact solution is
y ( x , ε ) = ∫ 0 x e − 1 ε ∫ 0 s A ( t ) d t d s ∫ 0 1 e − 1 ε ∫ 0 s A ( t ) d t d s . (6)
Since
y ′ ( x , ε ) = e − 1 ε ∫ 0 x A ( t ) d t ∫ 0 1 e − 1 ε ∫ 0 s A ( t ) d t d s > 0 , (7)
the solution y will increase monotonically with x. The asymptotic value of ∫ 0 x I ( s , ε ) d s in (6) is the area under the curve
I ( s , ε ) ≡ e − 1 ε ∫ 0 s A ( t ) d t (8)
for ε → 0 . Sophisticated techniques to obtain the asymptotic evaluation of integrals can be found in Olver [
The variety of limiting behaviors to singularly perturbed linear two-point boundary value problems with turning points has not been clearly described. The first papers by Pearson in 1968 stressed a numerical approach. In the intervening fifty years, software has improved tremendously, though finding the limiting solution is extremely ill-conditioned as Trefethen recently observed. Due to the serious instability of direct numerical methods, the examples found in scattered literature are usually less detailed. Inspired by this, in this paper we consider the asymptotic solution of two-point boundary value problems (4)-(5). In our examples, we’ll find the constant “outer” limits 0, 0.5, and 1.
Case 1: A ( x ) > 0 .
Here, I ( s , ε ) decays exponentially as ε → 0 , so for any fixed x > 0 , the numerator and denominator of (6) are both O ( ε ) and the ratio (6) is asymptotically one. Since y ( 0, ε ) = 0 , there is an initial boundary layer region of O ( ε ) thickness involving nonuniform convergence of y. Here, we’re using the big O order symbol.
As an example, take A ( x ) ≡ 1 , ε = 0.1 , 0.01 and 0.001, and plot the solution (6). One gets
y ( x , ε ) = 1 ε ∫ 0 x e − s ε d s 1 ε ∫ 0 1 e − s ε d s = 1 − e − x ε 1 − e − 1 ε . (9)
The constant limiting solution Y 0 ( x ) = 1 for x > 0 as ε → 0 , satisfies the reduced equation Y ′ 0 = 0 away from x = 0 . We plot the solution for three small ε values in
The limit of y ( x , ε ) is discontinuous at x = 0 , signaling nonuniform convergence.
Case 2: A ( x ) < 0 .
Now I ( x , ε ) grows exponentially large as ε → 0 . This causes y to be asymptotically zero for any x < 1 and a terminal boundary layer of nonuniform convergence to occur near x = 1 .
As an example, take A ( x ) ≡ − 1 and plot
y ( x , ε ) = e − 1 ε ε ∫ 0 x e s ε d s e − 1 ε ε ∫ 0 1 e s ε d s = e x − 1 ε − e − 1 ε 1 − e − 1 ε . (10)
For x < 1 , the limiting solution e x − 1 ε as ε → 0 is trivial. The limiting terminal layer will have O ( ε ) thickness.
Case 3: A ( x ) = x − 0.5 .
Now there’s a simple turning point at x = 0.5 and
I ( s , ε ) = e − 1 2 ε [ ( s − 1 2 ) 2 − 1 4 ] = e 1 8 ε e − 1 2 ε ( s − 1 2 ) 2 .
We write the ratio (6) as
y ( x , ε ) = ∫ 0 x I ( s , ε ) d s ∫ 0 1 I ( s , ε ) d s . (11)
The exact solution is
y ( x , ε ) = 1 2 ( 1 + e r f ( x − 1 2 2 ε ) e r f ( 1 2 2 ε ) ) , (12)
where e r f ( z ) = 2 π ∫ 0 z e − t 2 d t is the error function [
Since the integrands of (11) peak at the turning point and are asymptotically negligible elsewhere, we will have
y ( x , ε ) ∼ ( 0 , for x < 1 2 1 , for x > 1 2 . (13)
The numerator and denominator of (11) are both O ( ε ) . Clearly, y ( 1 2 , ε ) = 1 2 , there’s antisymmetry about s = 1 2 , and an O ( ε ) thick region of nonuniform convergence about the midpoint.
Plotting the solution (12) for ε = 10 − 3 , we get
Case 4: A ( x ) = x − α , 0 < α < 1 .
We rescale I to get
e − α 2 2 ε I ( s , ε ) = e − 1 2 ε ( s − α ) 2 ,
a function that peaks in an O ( ε ) interval about s = α and is asymptotically negligible elsewhere. This implies that a shock layer of nonuniform convergence occurs about the turning point. The exact solution (6) is
y ( x , ε ) = e r f ( x − α 2 ε ) + e r f ( α 2 ε ) e r f ( 1 − α 2 ε ) + e r f ( α 2 ε ) . (14)
For α = 1 3 and ε = 10 − 3 , we get
Not surprisingly, the asymptotic solution is essentially a translation of that for α = 0.5 . For α = 0 ,
y ( x , ε ) = ∫ 0 x e − s 2 2 ε d s ∫ 0 1 e − s 2 2 ε d s , (15)
so we again get an O ( ε ) initial layer (see
For α = 1 , there is an analogous terminal layer. For α < 0 or α > 1 , the boundary layer is O ( ε ) , i.e. thinner.
Case 5: A ( x ) = ( x − α ) 3 , 0 < α < 1 .
We have a third order turning point at x = α . Again, the rescaled integral I ( s , ε ) peaks at s = α , causing y to jump there. The shock layer is now
O ( ε 1 4 ) thick. The exact solution is
y ( x , ε ) = Γ ( 1 4 , ( x − α ) 4 4 ε ) − Γ ( 1 4 , α 4 4 ε ) Γ ( 1 4 , ( 1 − α ) 4 4 ε ) − Γ ( 1 4 , α 4 4 ε ) , (16)
where Γ ( a , z ) = ∫ z + ∞ u a − 1 e − u d u is the incomplete gamma function.
Evaluating y ( x , ε ) for α = 1 3 and ε = 10 − 3 , we get
To steepen the shock layer, we must take ε much smaller.
We change the sign of A for the next three examples.
Case 6: A ( x ) = 1 2 − x .
Rewriting (6) as
y ( x , ε ) = ∫ 0 x e s 2 ε ( s − 1 ) d s ∫ 0 1 e s 2 ε ( s − 1 ) d s , (17)
the exact solution is
y ( x , ε ) = 1 2 ( 1 + ∫ 0 x − 1 2 2 ε e t 2 d t ∫ 0 1 2 2 ε e t 2 d t ) . (18)
We note that the solution could be expressed in terms of Dawson’s integral F ( z ) = e − z 2 ∫ 0 z e t 2 d t . The integrands in (17) peak symmetrically at s = 0 and 1,
being asymptotically negligible elsewhere. Moreover, y ( 1 2 ) = 1 2 . Indeed y ( x ) ∼ 1 2 for 0 < x < 1 , and twin O ( ε ) boundary layers occur near both endpoints. For ε = 10 − 3 , we have
Case 7: For A ( x ) = α − x , 0 < α < 1 2 , we have
y ( x , ε ) = ∫ 0 x e s 2 ε ( s − 2 α ) d s ∫ 0 1 e s 2 ε ( s − 2 α ) d s . (19)
Its exact solution is
y ( x , ε ) = ∫ 0 x − α 2 ε e t 2 d t + ∫ 0 α 2 ε e t 2 d t ∫ 0 1 − α 2 ε e t 2 d t + ∫ 0 α 2 ε e t 2 d t . (20)
The integrand of (19) is asymptotically negligible for s < 2 α , but asymptotically large for s > 2 α . This implies that
y ( x ) ∼ 0 for x < 1,
so there is as O ( ε ) -thick terminal layer.
As an example, consider
10 − 3 y ″ + ( 1 3 − x ) y ′ = 0.
We have
Case 8: For A ( x ) = α − x , α > 1 2 , the integrand in (19) decays for 0 < s < 2 α . Thus
y ( x ) ∼ 1 for x > 0
and there is an O ( ε ) thick initial layer (see
For appropriate α ∼ 1 2 + o ( 1 ) , we’d expect that the shock layer moves across the interval. We’re now using the little o order symbol, which admittedly isn’t very explicit.
Case 9: For A ( ( x ) = ( x − 1 4 ) ( x − 3 4 ) , we have simple turning points at 1 4 and 3 4 . Moreover, I ( x , ε ) = e − x 3 ε ( x − 3 4 ) 2 peaks at x = 0 and 3 4 and is asymptotically negligible elsewhere. The sizes of the contributions to the integral differ, however. The area under I near x = 0 is O ( ε ) , but that near x = 3 4 is O ( ε ) , i.e., larger. Thus, the ratio (6) is O ( ε ) for 0 < x < 3 4 and O ( 1 ) for 3 4 < x < 1 .
Computing for ε = 10 − 3 , we get
This relies on the following figures. We’ve increased ε in
We have not been exhaustive, but we have certainly demonstrated a wide variety of asymptotic solutions to turning point problems of the form (4) - (5). They mimic the asymptotics of the more general boundary layer resonance problem. When the problem of turning points becomes complicated, numerical methods will become unreliable. Finding the limiting solution is extremely ill conditioned as Trefethen recently observed. Due to the serious instability of direct numerical methods, the examples found in scattered literature are usually less detailed. In this paper, we only give asymptotic solutions for a class of singularly perturbed with a turning point. Indeed, the techniques developed here might be expected to apply to that problem. Readers are encouraged to study other limiting possibilities for (4) - (5).
This research is supported by the Natural Science Foundation of Shanghai Institute of Technology, Research Fund nos. ZQ2018-22, 391100190016027.
The author declares no conflicts of interest regarding the publication of this paper.
Wang, N. (2019) Simple Singular Perturbation Problems with Turning Points. Journal of Applied Mathematics and Physics, 7, 2979-2989. https://doi.org/10.4236/jamp.2019.712208