Next Article in Journal
Fractional Whitham–Broer–Kaup Equations within Modified Analytical Approaches
Next Article in Special Issue
Axiomatic Approach in the Analytic Theory of Singular Perturbations
Previous Article in Journal
Advanced Control Strategies to Improve Nonlinear Automotive Dynamical Systems Consumption
Previous Article in Special Issue
Complete Controllability Conditions for Linear Singularly Perturbed Time-Invariant Systems with Multiple Delays via Chang-Type Transformation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Regularized Solution of Singularly Perturbed Cauchy Problem in the Presence of Rational “Simple” Turning Point in Two-Dimensional Case

by
Alexander Eliseev
†,‡ and
Tatjana Ratnikova
*,‡
Research Group on Macro-Structural Modeling of the Russian Economy, National Research University ”MPEI”, 111250 Moscow, Russia
*
Author to whom correspondence should be addressed.
Current address: National Research University ”MPEI”, Kracnokazarmennaya st., 14, 111250 Moscow, Russia.
These authors contributed equally to this work.
Axioms 2019, 8(4), 124; https://doi.org/10.3390/axioms8040124
Submission received: 17 September 2019 / Revised: 28 October 2019 / Accepted: 30 October 2019 / Published: 1 November 2019

Abstract

:
By Lomov’s S.A. regularization method, we constructed an asymptotic solution of the singularly perturbed Cauchy problem in a two-dimensional case in the case of violation of stability conditions of the limit-operator spectrum. In particular, the problem with a ”simple” turning point was considered, i.e., one eigenvalue vanishes for t = 0 and has the form t m / n a ( t ) (limit operator is discretely irreversible). The regularization method allows us to construct an asymptotic solution that is uniform over the entire segment [ 0 , T ] , and under additional conditions on the parameters of the singularly perturbed problem and its right-hand side, the exact solution.

1. Introduction

This work consists of five parts. The first part is an introduction. The second part is nomenclature. The third part presents the formulation of the Cauchy problem in the two-dimensional case if stability conditions for the spectrum of the limit operator are violated (the spectrum-stability condition means that eigenvalues of the operator A ( τ ) satisfy conditions λ 1 ( τ ) λ 2 ( τ ) , τ [ 0 , T ] and λ i 0 , i = 1 , 2 ).
A ”simple” pivot point of a limit operator (matrix A ( τ ) ) is understood when one eigenvalue vanishes at one point (i.e., matrix A ( τ ) is irreversible at this point). In [1], the case was considered of when one of the eigenvalues that had the form τ n a ( τ ) , a ( τ ) 0 , n was natural; in [2] the features of the solution were identified and described for a rational ”simple” turning point in the one-dimensional case (when the eigenvalue had the form τ m / n a ( τ ) , a ( τ ) 0 ).
In this article, we consider the case with a ”simple” turning point when one of the two eigenvalues of the operator vanishes at τ = 0 and has the form τ m / n a ( τ ) , a ( τ ) 0 .
The fourth part describes the formalism of the Lomov regularization method [1,3,4] that allows one to construct an asymptotic solution uniform over the entire segment [ 0 , T ] , under additional conditions on the parameters of a singularly perturbed problem, and its right side is the exact solution. The idea of this paper goes back to [1], in which methods were developed for solving a singularly perturbed Cauchy problem in the case of a ”simple” turning point of a limit operator with a natural exponent. A lemma is given on the estimation of basic singular functions, a theorem on the point solvability of iterative problems is proved, and the leading term of the asymptotic behavior of a singularly perturbed Cauchy problem is written out.
In the fifth part of the paper, we prove a theorem on the asymptotic behavior of a regularized series and a theorem on the passage to the limit as a small parameter tends to zero. For a parabolic equation, an example of solving a singularly perturbed Cauchy problem with a fractional turning point λ ( τ ) = τ 1 / 2 is given.
The sixth part is the conclusion.

2. Problem Formulation

Consider the Cauchy problem:
ε ¯ u ˙ ( τ ) = A ( τ ) u ( τ ) + h ( τ ) , u ( 0 , ε ) = u 0 ,
where
(1)
τ is a variable, τ [ 0 , T ] ;
(2)
u ( τ ) is a function, u ( τ ) C [ 0 , T ] ;
(3)
A ( τ ) is a matrix of size ( 2 × 2 ) , A ( τ ) C ( 0 , T ] ;
(4)
h ( τ ) is a function, h ( τ ) C [ 0 , T ] ;
(5)
λ 1 ( τ ) , λ 2 ( τ ) are eigenvalues of matrix A ( τ ) ; λ 1 ( τ ) λ 2 ( τ ) , τ [ 0 , T ] ; λ 2 ( τ ) = τ m / n a ( τ ) , where a ( τ ) < 0 , τ [ 0 , T ] , a ( τ ) C [ 0 , T ] ;
(6)
m, n are natural numbers;
(7)
Re λ 1 ( τ ) 0 ;
(8)
A ( t ) C [ 0 , T ] , where t = τ 1 / n ;
(9)
ε ¯ , ε = ε ¯ / n R there is a small parameter of the problem.
We make the change of variables in Problem (1): t = τ 1 / n . Then τ m / n = t m and
d u d τ = d u d t · d t d τ = u ˙ ( t ) 1 n τ ( 1 n ) / n = u ˙ ( t ) 1 n t 1 n .
Equation (1) takes the form:
ε ¯ n u ˙ ( t ) t 1 n = A ( t n ) u ( t ) + h ( t n )
or
ε ¯ n u ˙ ( t ) = t n 1 A ( t n ) u ( t ) + t n 1 h ( t n ) .
Denote ε ¯ / n = ε , t n 1 A ( t n ) = B ( t ) . Task (1) takes the form:
ε u ˙ ( t ) = B ( t ) u ( t ) + t n 1 h ( t n ) , u ( 0 , ε ) = u 0 .
Operator B ( t ) has eigenvalues λ ¯ 1 ( t ) = t n 1 λ 1 ( t n ) , λ ¯ 2 ( t ) = t p a ( t n ) , where p = m + n 1 , and corresponding vectors e ¯ 1 ( t ) = e 1 ( t n ) , e ¯ 2 ( t ) = e 2 ( t n ) , where e 1 ( τ ) , e 2 ( τ ) are eigenvectors of operator A ( τ ) , i.e.,
B ( t ) e ¯ 1 ( t ) = λ ¯ 1 ( t ) e ¯ 1 ( t ) = t n 1 λ 1 ( t n ) e 1 ( t n ) ; B ( t ) e ¯ 2 ( t ) = λ ¯ 2 ( t ) e ¯ 2 ( t ) = t p a ( t n ) e 1 ( t n ) .
Methods for solving the Cauchy problem (2) are described in [1]. Basic singularities (2) have the form:
e φ i ( t ) / ε , i = 1 , 2 ; σ i ( t , ε ) = e φ 2 ( t ) / ε 0 t e φ 2 ( s ) / ε s i d s , i = 0 , p 1 ¯ ;
where φ 1 ( t ) = 0 t s n 1 λ 1 ( s n ) d s ,   φ 2 ( t ) = 0 t s p a ( s n ) d s .
Singularities (3) in the source variables have the form
e φ 1 ( τ ) / ε ¯ , e φ 2 ( τ ) / ε ¯ , σ i ( τ , ε ) = e φ 2 ( τ ) / ε ¯ 0 τ e φ 2 ( s ) / ε s ( i + 1 n ) / n d s , i = 0 , p 1 ¯ ;
where φ 1 ( τ ) = 0 τ λ 1 ( s ) d s ,   φ 2 ( τ ) = 0 τ a ( s ) s m / n d s .

3. Formalism of Regularization Method

Point ε = 0 for Problem (1) is special in the sense that classical existence theorems for the solution of the Cauchy problem do not take place. Therefore, in solving this problem, essentially singular singularities arise. When the stability condition for spectrum A ( t ) is satisfied, singular singularities are described using exponentials of the form:
e φ i ( t ) / ε , φ i ( t ) = 0 t λ i ( s ) d s , i = 1 , 2 ¯ , λ 1 ( t ) λ 2 ( t ) , λ i ( t ) 0 , t [ 0 , T ] ,
where φ i ( t ) is a smooth function (in the general case, complex) of a real variable t. To solve linear homogeneous equations, such singularities have been described by Liouville [5,6,7,8].
If stability conditions are violated for at least one point of the spectrum of operator A ( t ) , then besides exponentially essentially singularities in the solution of the inhomogeneous equation, singularities of the following form also appear:
σ i = e φ 1 ( t ) / ε 0 t e φ 1 ( s ) / ε s i d s , i = 0 , k 1 ¯ ,
(k is the extreme zero of λ 1 ( t ) ), which, for ε 0 , has a power character of decreasing under the corresponding restrictions on λ 1 ( t ) , while it is assumed that the remaining points of the spectrum do not vanish at t = 0 .
Singularly perturbed problems arise in cases when the domain of definition of the initial operator, depending on ε with ε 0 , does not coincide with the domain of definition of the limit operator with ε = 0 . When studying problems with a ”simple” turning point, additional conditions arise when the domain of values of the original operator does not coincide with the domain of values of the limit operator.
Further, we need estimates of functions describing the basic singularities.
Lemma 1.
Let the conditions on the spectrum of operator A ( t )   5 ) ÷ 7 ) be satisfied. Then, the estimates hold:
(a) if t [ 0 , T ]   Re λ i ( t ) 0 , i = 1 , 2 , then
| e 1 ε 0 t λ ¯ i ( s ) d s | C , | σ k ( t , ε ) | C ,
where C is a constant, k = 0 , p 1 ¯ , p = m + n 1 ;
(b) if Re λ 1 ( t ) α < 0 , Re a ( t ) α < 0 , then
| e 1 ε 0 t λ ¯ 1 ( s ) d s | e α t n ε n , | e 1 ε 0 t λ ¯ 2 ( s ) d s | e α t p + 1 ε ( p + 1 ) , | σ k ( t , ε ) | C ε k + 1 p + 1 , k = 0 , p 1 ¯ , p = m + n 1 .
Proof of Lemma 1. 
(a) In this case, estimates are obvious.
(b) | e 1 ε 0 t λ ¯ 1 ( s ) d s | e α ε 0 t s n 1 d s = e α t n ε n , | e 1 ε 0 t λ ¯ 2 ( s ) d s | e α ε 0 t s p d s = e α t p + 1 ε ( p + 1 ) ,
| σ k ( t , s ) | = | 0 t e 1 ε s t λ ¯ 2 ( s ) d s s k d s | 0 t e α ε s t s 1 p d s 1 s k d s = 0 t e α ( t p + 1 s p + 1 ) ε ( p + 1 ) s k d s = 0 t e α s p + 1 ε ( p + 1 ) s k d s e α t p + 1 ε ( p + 1 ) =
= ε k + 1 p + 1 0 t / ε 1 / ( p + 1 ) e α ξ p + 1 p + 1 ξ k d ξ e α t p + 1 ε ( p + 1 ) .
Denote τ = t ε 1 / ( p + 1 ) . Consider a fraction when τ ; then, we have
0 τ e α ξ p + 1 p + 1 ξ k d ξ e α τ p + 1 p + 1 τ k α τ p τ 0 , as k < p .
Consequently, σ k ( t , ε ) = O ̲ ( ε k + 1 p + 1 ) . □
Remark 1.
Estimates in the source variables have the form:
| e 1 ε 0 t λ 1 ( s ) d s | e α t ε , | e 1 ε 0 t λ 2 ( s ) d s | e α t m / n + 1 ε ( m / n + 1 ) , σ k ( t , ε ) = O ̲ ( ε k + 1 m / n + 1 ) .
According to the regularization method, we seek a solution of Problem (2) in the form
u ( t , ε ) = x ( t , ε ) e φ 1 ( t ) / ε + y ( t , ε ) e φ 2 ( t ) / ε + i = 0 p 1 z i ( t , ε ) σ i ( t , ε ) + W ( t , ε ) ,
where x ( t , ε ) , y ( t , ε ) , W ( t , ε ) , z i ( t , ε ) , i = 0 , p 1 ¯ are smooth with respect to t functions that depend on power on ε . Substituting Problem (4) into Problem (2), we get system
( B ( t ) λ ¯ 1 ( t ) ) x ( t , ε ) = ε x ˙ ( t , ε ) , ( B ( t ) λ ¯ 2 ( t ) ) y ( t , ε ) = ε y ˙ ( t , ε ) , ( B ( t ) λ ¯ 2 ( t ) ) z i ( t , ε ) = ε z ˙ i ( t , ε ) , i = 0 , p 1 ¯ , B ( t ) W ( t , ε ) = ε W ˙ ( t , ε ) t n 1 h ( t n ) + i = 0 p 1 t i z i ( t , ε ) , x ( 0 , ε ) + y ( 0 , ε ) + W ( 0 , ε ) = u 0 .
Decomposing the unknown vector functions in a series in powers of ε , we obtain a series of iterative problems:
( B ( t ) λ ¯ 1 ( t ) ) x k ( t ) = x ˙ k 1 ( t ) , ( B ( t ) λ ¯ 2 ( t ) ) y k ( t ) = y ˙ k 1 ( t ) , ( B ( t ) λ ¯ 2 ( t ) ) z k i ( t , ε ) = z ˙ k 1 i ( t ) , i = 0 , p 1 ¯ , B ( t ) W k ( t ) = W ˙ k 1 ( t ) δ 0 k t n 1 h ( t n ) + i = 0 p 1 t i z k 1 i ( t ) , x k ( 0 ) + y k ( 0 ) + W k ( 0 ) = δ k 0 u 0 .
To solve iterative Problems (6), we formulate a point-solvability theorem.
Theorem 1.
Let the following equation be given:
B ( t ) u ( t ) t n 1 A ( t n ) = t n s h ( t n ) , 0 s n 1
and let the following conditions are met:
(1) 
B ( t ) has eigenvalues λ ¯ 1 ( t ) = t n 1 λ 1 ( t n ) , λ ¯ 2 ( t ) = t p a ( t n ) and eigenvectors e ¯ 1 ( t ) , e ¯ 2 ( t ) ;
(2) 
h ( t n ) C [ 0 , T ] .
Then, Problem (7) is solvable if and only if
(a) 
h 1 ( 0 ) = 0 , s = 2 , n 1 ¯ ;
(b) 
h 2 ( k ) = 0 , k = 0 , [ m + s 1 n ] ¯ , s = 0 , n 1 ¯ ,
where h 1 ( t n ) , h 2 ( t n ) are the components of decomposition h ( t ) on the basis of eigenvectors of operator B ( t ) ; u 1 ( t n ) , u 2 ( t n ) are the components of the expansion of u ( t ) on the basis of eigenvectors of operator B ( t ) .
Proof of Theorem 1. 
Let us prove the need. Let system
t n 1 λ 1 ( t n ) u 1 ( t ) = t n s h 1 ( t n ) , t p a ( t n ) u 2 ( t ) = t n s h 2 ( t n )
have a solution. Then,
(1)
the first equation of System (8) is solvable:
(a)
if s = 0 , 1 , then u 1 ( t ) = t 1 s h 1 ( t n ) λ 1 ( t n ) ,
(b)
if s = 2 , n 1 ¯ , then h 1 ( 0 ) = 0 and u 1 ( t ) = t n + 1 s h ¯ 1 ( t n ) λ 1 ( t n ) , where h 1 ( t n ) = t n h ¯ 1 ( t n ) ;
(2)
the second equation of System (8) is solvable if ( k + 1 ) n s p < ( k + 2 ) n s , which is equivalent to h 2 ( k ) ( 0 ) = 0 , k = 0 , [ m + s 1 n ] ¯ and u 2 ( t ) = t n j h ¯ 2 ( t n ) a ( t n ) , 0 j n 1 .
Sufficiency is obvious. □
Consider Problem (6) as k = 1 :
( B ( t ) λ ¯ 1 ( t ) ) x 1 ( t ) = 0 , ( B ( t ) λ ¯ 2 ( t ) ) y 1 ( t ) = 0 , ( B ( t ) λ ¯ 2 ( t ) ) z 1 i ( t ) = 0 , i = 0 , p 1 ¯ , B ( t ) W 1 ( t ) = 0 , x 1 ( 0 ) + y 1 ( 0 ) + W 1 ( 0 ) = 0 .
Solution (9) has the form
x 1 ( t ) = α 1 1 ( t ) e ¯ 1 ( t ) , y 1 ( t ) = β 1 2 ( t ) e ¯ 2 ( t ) , z 1 i ( t ) = γ 1 i , 2 ( t ) e ¯ 2 ( t ) , W 1 ( t ) 0 , α 1 1 ( 0 ) = 0 , β 1 2 ( 0 ) = 0 .
Functions x 1 ( t ) , y 1 ( t ) , z 1 i ( t ) are determined at the next iteration step k = 0 from the solvability conditions:
( B ( t ) λ ¯ 1 ( t ) ) x 0 ( t ) = x ˙ 1 ( t ) , ( B ( t ) λ ¯ 2 ( t ) ) y 0 ( t ) = y ˙ 1 ( t ) , ( B ( t ) λ ¯ 2 ( t ) ) z 0 i ( t ) = z ˙ 1 i ( t ) , i = 0 , p 1 ¯ , B ( t ) W 0 ( t ) = t n h ( t n ) + i = 0 p 1 t i z 1 i ( t ) , x 0 ( 0 ) + y 0 ( 0 ) + W 0 ( 0 ) = u 0 .
Let be
e ¯ ˙ i ( t ) = e ˙ i ( t n ) = n t n 1 C i 1 ( t n ) e ¯ 1 ( t ) + n t n 1 C i 2 ( t n ) e ¯ 2 ( t ) , i = 1 , 2 .
Denote by C ¯ i j ( t ) = n t n 1 C i j ( t n ) , i , j = 1 , 2 . Then
e ¯ ˙ i ( t ) = j = 1 2 C ¯ i j ( t ) e ¯ j ( t ) , i = 1 , 2 .
System (10) takes the form:
( B ( t ) λ ¯ 1 ( t ) ) x 0 ( t ) = ( α ˙ 1 1 ( t ) + C ¯ 1 1 ( t ) α 1 1 ( t ) ) e ¯ 1 ( t ) + α 1 1 ( t ) C ¯ 1 2 ( t ) e ¯ 2 ( t ) , ( B ( t ) λ ¯ 2 ( t ) ) y 0 ( t ) = ( β ˙ 1 2 ( t ) + C ¯ 2 2 ( t ) β 1 2 ( t ) ) e ¯ 2 ( t ) + β 1 2 ( t ) C ¯ 2 1 ( t ) e ¯ 1 ( t ) , ( B ( t ) λ ¯ 2 ( t ) ) z 0 i ( t ) = ( γ ˙ 1 i , 2 ( t ) + C ¯ 2 2 ( t ) γ 1 i , 2 ( t ) ) e ¯ 2 ( t ) + γ 1 i , 2 ( t ) C ¯ 2 1 ( t ) e ¯ 1 ( t ) , i = 0 , p 1 ¯ , B ( t ) W 0 ( t ) = t n h ( t n ) + i = 0 p 1 t i z 1 i ( t ) , x 0 ( 0 ) + y 0 ( 0 ) + W 0 ( 0 ) = u 0 .
The conditions for the solvability of System (11) and the initial conditions at the k = 1 step imply that α 1 1 ( t ) 0 , β 1 2 ( t ) 0 . To determine z 1 i ( t ) , we wrote by coordinate the equation for the W 0 ( t ) of System (11):
λ ¯ 1 ( t ) W 0 1 ( t ) t n 1 λ 1 ( t n ) W 0 1 ( t ) = t n 1 h 1 ( t n ) , λ ¯ 2 ( t ) W 0 2 ( t ) t p a ( t n ) W 0 2 ( t ) = t n 1 h 2 ( t n ) + i = 0 p 1 t i γ 1 i , 2 ( t ) .
Then, W 0 1 ( t ) = h 1 ( t n ) λ 1 ( t n ) . On the basis of the point-solvability theorem, we obtained:
γ 1 n 1 , 2 ( 0 ) = h 2 ( 0 ) , γ 1 2 n 1 , 2 ( 0 ) = h ˙ 2 ( 0 ) , , γ 1 ( k + 1 ) n 1 , 2 ( 0 ) = h 2 ( k ) ( 0 ) k ! e 0 t C ¯ 2 2 ( s ) d s ,
where k = [ m / n ] is the integer part, so when order ord ( t i ) is equal to order ord ( t ( j + 1 ) n 1 ) in the expansion of t n 1 h ( t n ) in Taylor–Maclaurin series, other γ 1 i , 2 ( 0 ) = 0 . Thus, the solution is determined at step k = 1 :
u 1 ( t , ε ) = i = 0 [ m / n ] γ 1 ( i + 1 ) n 1 , 2 ( t ) e ¯ 2 ( t ) σ ( i + 1 ) n 1 ( t , ε ) ,
where γ 1 ( i + 1 ) n 1 , 2 ( t ) = h 2 ( i ) ( 0 ) i ! e 0 t C ¯ 2 2 ( s ) d s .
The solution at zero step k = 0 is written in the form
x 0 ( t ) = α 0 1 ( t ) e ¯ 1 ( t ) , y 0 ( t ) = β 0 2 ( t ) e ¯ 2 ( t ) , z 0 i ( t ) = γ 0 i , 2 ( t ) e ¯ 2 ( t ) , ord ( t i ) ord ( t ( j + 1 ) n 1 ) , γ 0 i , 2 ( t ) e ¯ 2 ( t ) γ 1 i , 2 ( t ) C ¯ 2 1 ( t ) λ ¯ 1 ( t ) λ ¯ 2 ( t ) e ¯ 1 ( t ) , ord ( t i ) = ord ( t ( j + 1 ) n 1 ) , i = 0 , p 1 ¯ , j = 0 , [ m / n ] ¯ , W 0 ( t ) = h 1 ( t n ) λ 1 ( t n ) e ¯ 1 ( t ) + t n s H 0 ( t n ) e ¯ 2 ( t ) ,
where
(a)
s = { m n } n is the remainder of dividing m by n;
(b)
t n s H 0 ( t n ) t n 1 h ( t n ) + i = 0 [ m / n ] t ( i + 1 ) n 1 γ 1 ( i + 1 ) n 1 , 2 ( t ) t p a ( t n ) .
Arbitrary functions α 0 1 ( t ) , β 0 2 ( t ) , γ 0 i , 2 ( t ) are determined from the conditions for the solvability of the system at step k = 1 :
( B ( t ) λ ¯ 1 ( t ) ) x 1 ( t ) = ( α ˙ 0 1 ( t ) + C ¯ 1 1 ( t ) α 0 1 ( t ) ) e ¯ 1 ( t ) + α 0 1 ( t ) C ¯ 1 2 ( t ) e ¯ 2 ( t ) , ( B ( t ) λ ¯ 2 ( t ) ) y 1 ( t ) = ( β ˙ 0 2 ( t ) + C ¯ 2 2 ( t ) β 0 2 ( t ) ) e ¯ 2 ( t ) + β 0 2 ( t ) C ¯ 2 1 ( t ) e ¯ 1 ( t ) , ( B ( t ) λ ¯ 1 ( t ) ) z 1 i ( t ) = ( γ ˙ 0 i , 2 ( t ) + C ¯ 2 2 ( t ) γ 0 i , 2 ( t ) ) e ¯ 2 ( t ) + γ 0 i , 2 ( t ) C ¯ 2 1 ( t ) e ¯ 1 ( t ) , ord ( t i ) ord ( t ( j + 1 ) n 1 ) , ( B ( t ) λ ¯ 2 ( t ) ) z 1 i ( t ) = ( γ ˙ 0 i , 2 ( t ) + C ¯ 2 2 ( t ) γ 0 i , 2 ( t ) γ 1 i , 2 ( t ) C ¯ 2 1 ( t ) C ¯ 1 2 ( t ) λ ¯ 1 ( t ) λ ¯ 2 ( t ) ) e ¯ 2 ( t ) + + ( γ ˙ 0 i , 2 ( t ) C ¯ 2 1 ( t ) ( γ 1 i , 2 ( t ) C ¯ 2 1 ( t ) λ ¯ 1 ( t ) λ ¯ 2 ( t ) ) · γ 1 i , 2 ( t ) C ¯ 2 1 ( t ) C ¯ 1 1 ( t ) λ ¯ 1 ( t ) λ ¯ 2 ( t ) ) e ¯ 1 ( t ) , ord ( t i ) = ord ( t ( j + 1 ) n 1 ) , i = 0 , p 1 ¯ , j = 0 , [ m / n ] ¯ , B ( t ) W 1 ( t ) = W ˙ 0 ( t ) + i = 0 p 1 t i z 0 i ( t ) , x 1 ( 0 ) + y 1 ( 0 ) + W 1 ( 0 ) = 0 .
The solvability theorem of System (15) gives
α 0 1 ( t ) = ( u 1 0 + h 1 ( 0 ) λ 1 ( 0 ) ) e 0 t C ¯ 1 1 ( s ) d s , β 0 1 ( t ) u 2 0 e 0 t C ¯ 2 2 ( s ) d s .
Consider the equation for W 1 ( t ) . Given the expression for C ¯ i j ( t ) = n t n 1 C i j ( t n ) , this equation can be written as follows:
B ( t ) W 1 ( t ) = t n 1 ( W ˙ 0 ) 1 ( t ) e ¯ 1 ( t ) + t n 1 s ( W ˙ 0 ) 2 ( t ) e ¯ 2 ( t ) + i = 0 p 1 t i z 0 i ( t ) .
Consider Equation (16) component-wise:
λ ¯ 1 ( t ) W 1 1 ( t ) = t n 1 ( W ˙ 0 ) 1 ( t ) i = 0 [ m / n ] γ 1 ( i + 1 ) n 1 , 2 ( t ) C ¯ 2 1 ( t ) t ( i + 1 ) n 1 λ ¯ 1 ( t ) λ ¯ 2 ( t ) , λ ¯ 2 ( t ) W 1 2 ( t ) = t n 1 s ( W ˙ 0 ) 2 ( t ) + i = 0 p 1 t i γ 0 i , 2 ( t ) .
Solution of the first equation of System (17) is written as follows:
W 1 1 ( t ) = ( W ˙ 0 ) 1 ( t ) λ 1 ( t n ) i = 0 [ m / n ] γ 1 ( i + 1 ) n 1 , 2 ( t ) n C 2 1 ( t n ) t i n λ 1 ( t n ) ( λ 1 ( t n ) t m a ( t n ) ) .
For the solvability of the second equation of System (17), it is necessary and sufficient that
γ 0 ( i + 1 ) n 1 s , 2 ( 0 ) = ( W ˙ 0 ) 2 ( i ) ( 0 ) i ! , i = 0 , [ m + s n ] ¯ ,
here [ m + s n ] = [ [ m n ] + 2 s n ] = [ m n ] + [ 2 s n ] .
The other γ 0 j , 2 ( 0 ) = 0 , j ( i + 1 ) n 1 s , j = 0 , p 1 ¯ . Defining γ 0 i , 2 ( 0 ) , we can write the expression for z 0 i ( t ) :
(a)
if j = ( i + 1 ) n 1 s , i = 0 , [ m n ] + [ 2 s n ] ¯ , then
γ 0 j , 2 ( t ) = ( W ˙ 0 ) 2 ( i ) ( 0 ) i ! e 0 t C ¯ 2 2 ( s ) d s , z 0 j ( t ) = γ 0 j , 2 ( t ) C ¯ 2 ( t ) ;
(b)
if j ( i + 1 ) n 1 s , j = ( i + 1 ) n 1 , then
γ 0 j , 2 ( t ) = e 0 t C ¯ 2 2 ( s ) d s 0 t e 0 s C ¯ 2 2 ( z ) d s γ 1 j , 2 ( s ) C ¯ 2 1 ( s ) C ¯ 1 2 ( s ) λ ¯ 1 ( s ) λ ¯ 2 ( s ) d s ,
z 0 j ( t ) = γ 0 j , 2 ( t ) e ¯ 2 ( t ) γ 1 j , 2 ( t ) C ¯ 2 1 ( t ) λ ¯ 1 ( t ) λ ¯ 2 ( t ) e ¯ 1 ( t ) ;
(c)
if j ( i + 1 ) n 1 s , j ( i + 1 ) n 1 , then
γ 0 j , 2 ( t ) 0 , z 0 j , 2 ( t ) 0 .
The solution of the second equation of System (17) is written as follows:
W 1 2 ( t ) = t n 1 2 s n H 1 ( t n ) ,
where H 1 ( t n ) = t n s 1 ( W ˙ 0 ) 2 i = 0 p 1 t i γ 0 i , 2 ( t ) t p a ( t n ) .
Thus, the solution is determined at the zero iterative step:
u 0 ( t , ε ) = α 0 1 ( t ) e ¯ 1 ( t ) e φ 1 ( t ) / ε + β 0 2 ( t ) e ¯ 2 ( t ) e φ 2 ( t ) / ε + i = 0 m n + 2 s n z 0 ( i + 1 ) n 1 s ( t ) σ ( i + 1 ) n 1 s ( t , ε ) + i = 0 m n z 0 ( i + 1 ) n 1 ( t ) σ ( i + 1 ) n 1 ( t , ε ) h 1 ( t n ) λ 1 ( t n ) e ¯ 1 ( t ) + t n s H 0 ( t n ) e ¯ 2 ( t ) .
Similarly, according to this scheme, the solutions of subsequent iteration problems are determined. Thus, we can get an expression for any member of a regularized series.
We write the main term of the asymptotics of Problem (2):
u main = 1 ε u 1 ( t , ε ) + u 0 ( t , ε ) .

4. Limit-Transition Theorem

To prove the asymptoticity of a regularized series, we prove a theorem on estimating the remainder term for ε 0 .
Let be u ( t , ε ) = k = 1 n ε k u k ( t , ε ) + ε n + 1 R n ( t , ε ) , where
u k ( t , ε ) = x k ( t ) e φ 1 ( t ) / ε + y k ( t ) e φ 2 ( t ) / ε + i = 0 p 1 z k i ( t ) σ i ( t , ε ) + W k ( t ) .
Substituting Problem (18) into Problem (1), we obtain the Cauchy problem for the remainder R n ( t , ε ) :
ε R ˙ n ( t , ε ) = B ( t ) R n ( t , ε ) + H ( t , ε ) , R ( 0 , ε ) = 0 ,
where
H ( t , ε ) = x ˙ n ( t ) e φ 1 ( t ) / ε + y ˙ n ( t ) e φ 2 ( t ) / ε + i = 0 p 1 z ˙ n ( t ) σ i ( t , ε ) + W ˙ n ( t ) + i = 0 p 1 t i z n i ( t ) ,
in this case, it is assumed that H ( t , ε ) satisfies the conditions of the solvability theorem.
Theorem 2.
Let Cauchy Problem (1) be given and Conditions 1–9 be satisfied. Then, the estimate is correct
u ( t , ε ) k = 1 n ε k u k ( t , ε ) C ε n + 1 ,
where C > 0 in the norm C [ 0 , T ] for any ( t , ε ) [ 0 , T ] × ( 0 , ε 0 ] , x ( t ) C [ 0 , T ] = max t [ 0 , T ] | x ( t ) | .
Proof of Theorem 2. 
Solution (19) is written as follows:
R n ( t , ε ) = 1 ε 0 t U ε ( t , s ) H ( s , ε ) d s ,
where U ε ( t , s ) is resolving operator (fundamental solution system) satisfying system
ε U ˙ ε ( t , s ) = B ( t ) U ε ( t , s ) , U ε ( t , s ) | s = t = I .
Let S ( t ) be a matrix of eigenvectors e ¯ 1 ( t ) , e ¯ 2 ( t ) of operator B ( t ) . Then, System (21) is equivalent to system
ε V ˙ ε ( t , s ) = Λ ( t ) V ε ( t , s ) ε S 1 ( t ) S ˙ ( t ) V ε ( t , ε ) , V ε ( t , s ) | s = t = S 1 ( 0 ) ,
here, Λ ( t ) = λ ¯ 1 ( t ) 0 0 λ ¯ 2 ( t ) , V ε ( t , s ) = S 1 ( t ) U ε ( t , s ) . We reduce System (22) to an integral equation
V ε ( t , s ) = e 1 ε s t Λ ( s 1 ) d s 1 S 1 ( 0 ) s t e 1 ε s 1 t Λ ( s 2 ) d s 2 S 1 ( s 1 ) S ˙ ( s 1 ) V ε ( s 1 , ε ) d s 1 .
Let us estimate Equation (23) at the norm C [ 0 , T ] . Using the conditions on the spectrum of operator B ( t ) , we obtain
V ε ( t , s ) C 1 S 1 ( 0 ) + C 2 s t V ε ( s 1 , s ) d s 1 .
Using the Bellman–Gronuola inequality, we obtain U ε ( t , s ) C on [ 0 , T ] . To estimate the remaining term, it is important to take into account that operator B ( t ) is invertible on vector functions that satisfy the conditions of the solvability theorem. Then, integrating over parts of Solution (20), we obtain chain of equalities
R n ( t , ε ) = 1 ε 0 t U ε ( t , s ) H ( s , ε ) d s = 1 ε 0 t U ε ( t , s ) B ( s ) B 1 ( s ) H ( s , ε ) d s = = U ε ( t , s ) B 1 ( s ) H ( s , ε ) | 0 t + 0 t U ε ( t , s ) d d s B 1 ( s ) H ( s , ε ) d s = = B 1 ( t ) H ( t , ε ) + U ε ( t , s ) B 1 ( s ) H ( s , ε ) | s = 0 + 0 t U ε ( t , s ) d d s B 1 ( s ) H ( s , ε ) d s .
Since, by virtue of Conditions 1–9, H ( t , ε ) admits estimate H ( t , ε ) C 1 in norm C [ 0 , T ] , then remainder R n ( t , ε ) satisfies estimate
R n ( t , ε ) C 2 ( t , ε ) [ 0 , T ] × ( 0 , ε 0 ] .
Therefore, the asymptoticity of series k = 1 ε k u k ( t , ε ) is proved. □
Theorem 3 (The limit theorem). 
Let Cauchy Problem (1) be given and have satisfied the conditions:
(1) 
Conditions 1–9;
(2) 
h 2 ( i ) ( 0 ) = 0 , i = 0 , [ m / n ] ¯ , where h 2 ( t ) is the second coordinate in the expansion of h ( t ) = h 1 ( t ) e 1 ( t ) + h 2 ( t ) e 2 ( t ) in eigenvectors of the original matrix.
Then,
(1) 
for any δ > 0   t [ δ , T ] , Re λ i ( t ) α < 0
lim ε 0 u ( t , ε ) = A 1 ( t ) h ( t ) ;
(2) 
if Re λ i ( t ) = 0 , then
u ( t , ε ) ε 0 w e a k A 1 ( t ) h ( t ) i n a w e a k s e n s e .
Proof of Theorem 3 
(1) Conditions h 2 ( i ) ( 0 ) = 0 , i = 0 , [ m / n ] ¯ cause u 1 ( t , ε ) = 0 . Then,
u main ( t ) = u 0 ( t , ε ) .
By virtue of the singularity estimates described in the lemma, it follows that for any δ > 0   t [ δ , T ]
lim ε 0 u 0 ( t , ε ) = B 1 ( t ) t n 1 h ( t n ) ,
equivalent in source variables lim ε 0 u 0 ( t , ε ) = A 1 ( t ) h ( t ) .
(2) If Re λ i ( t ) 0 , i = 1 , 2 , then singularities are rapidly oscillating exponents as ε 0 . From here, according to Lebesgue’s lemma, for any φ ( t ) C ( 0 , T )
0 T ( u 0 ( t , ε ) + A 1 ( t ) h ( t ) ) φ ( s ) d t ε 0 0 .
 □
Example 1.
Consider the Cauchy problem for a parabolic equation
ε u t ε 2 2 u x 2 = t u + h ( x , t ) , u ( x , 0 ) = φ ( x ) , < x < ,
where φ ( x ) , h ( x , t ) C 0 ( , ) are smooth functions with compact support.Using the technique of the regularization method outlined above, we obtain the principal term of the asymptotics of the solution:
u ( x , t ) = 1 ε h ( x , 0 ) e 2 3 ε t 3 / 2 0 t e 2 s 3 / 2 3 ε d s + φ ( x ) e 2 t 3 / 2 3 ε h ˙ ( x , 0 ) 2 e 2 t 3 / 2 3 ε 0 t e 2 s 3 / 2 3 ε d s s +
+ t h ( x , 0 ) e 2 t 3 / 2 3 ε 0 t e 2 s 3 / 2 3 ε d s + h ( x , t ) h ( x , 0 ) t .

5. Conclusions

In this paper, the regularization method was developed into the class of singularly perturbed Cauchy problems in the case of a simple rational turning point for the limit operator (for ε = 0 ). The main singularities of the solution are highlighted:
e φ 1 ( t ) / ε , e φ 2 ( t ) / ε , σ i ( t , ε ) = e φ 2 ( t ) / ε 0 t e φ 2 ( s ) / ε s ( i + 1 n ) / n d s , i = 0 , p 1 ¯ ,
which allowed us to present the solution in the form:
u ( t , ε ) = x ( t , ε ) e φ 1 ( t ) / ε + y ( t , ε ) e φ 2 ( t ) / ε + i = 0 p 1 z i ( t , ε ) σ i ( t , ε ) + W ( t , ε ) ,
where x ( t , ε ) , y ( t , ε ) , W ( t , ε ) , z i ( t , ε ) , i = 0 , p 1 ¯ are t smooth functions that depend on power ε .
Estimates of the main singularities for ε 0 were given, and theorems on the solvability of iterative problems were proved. A theorem on the asymptotic convergence of the solution of the problem was proved, and conditions on the right-hand side of h ( t ) were described, under which the passage to the limit theorem is valid. An example of solving the Cauchy problem for a parabolic equation with a fractional turning point λ ( τ ) = τ 1 / 2 was given.

Author Contributions

Conceptualization and writing—review and editing by A.E. and T.R.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

u ( t ) , x ( t ) , y ( t ) , z i ( t ) , i = 0 , m + n 1 ¯ , w ( t ) , h ( t ) a vector of a function of a real variable
A ( t ) , B ( t ) matrices of order 2 × 2
λ 1 , λ 2 eigenvalues of matrix A
λ ¯ 1 , λ ¯ 2 eigenvalues of matrix B
e 1 , e 2 eigenvectors of matrix A
e ¯ 1 , e ¯ 2 eigenvectors of matrix B
ε a small task parameter
S ( t ) a matrix of eigenvectors e ¯ 1 , e ¯ 2
Λ ( t ) a matrix of eigenvalues of matrix B
R n the remainder term of the asymptotic series

References

  1. Eliseev, A.G.; Lomov, S.A. The theory of singular perturbations in the case of spectral singularities of the limit operator. Sb. Math. 1986, 131, 544–557. [Google Scholar] [CrossRef]
  2. Eliseev, A.G.; Ratnikova, T.A. A singularly perturbed Cauchy problem in the presence of a rational ”simple” turning point for the limit operator. Differ. Equ. Control Process. 2019, 3, 63–73. [Google Scholar]
  3. Lomov, S.A. Introduction to the General Theory of Singular Perturbations; Nauka: Moscow, Russia, 1981. [Google Scholar]
  4. Safonov, V.F.; Bobodzhanov, A.A. Course of higher mathematics. In Singularly Perturbed Equations and the Regularization Method: Textbook; Izdatelstvo MPEI: Moscow, Russia, 2012. [Google Scholar]
  5. Liouville, J. Second memoire sur le development des fonctions en series dont divers termes sont assujettis, a une meme equation. J. Math. Pure Appl. 1837, 2, 16–35. [Google Scholar]
  6. Lomov, S.A.; Eliseev, A.G. Asymptotic integration of singularly perturbed problems. Russ. Math. Surv. 1988, 43, 1–63. [Google Scholar] [CrossRef]
  7. Eliseev, A.G. Singular perturbation theory for systems of differential equations in the case of multiple spectrum of the limit operator. III. Math. USSR-Izv. 1985, 25, 475–500. [Google Scholar] [CrossRef]
  8. Tursunov, D.A.; Kozhbekov, K.G. Asymptotics of the solution of singularly perturbed differential equations with a fractional turning point. News Irkutsk State Univ. 2017, 21, 108–121. [Google Scholar]

Share and Cite

MDPI and ACS Style

Eliseev, A.; Ratnikova, T. Regularized Solution of Singularly Perturbed Cauchy Problem in the Presence of Rational “Simple” Turning Point in Two-Dimensional Case. Axioms 2019, 8, 124. https://doi.org/10.3390/axioms8040124

AMA Style

Eliseev A, Ratnikova T. Regularized Solution of Singularly Perturbed Cauchy Problem in the Presence of Rational “Simple” Turning Point in Two-Dimensional Case. Axioms. 2019; 8(4):124. https://doi.org/10.3390/axioms8040124

Chicago/Turabian Style

Eliseev, Alexander, and Tatjana Ratnikova. 2019. "Regularized Solution of Singularly Perturbed Cauchy Problem in the Presence of Rational “Simple” Turning Point in Two-Dimensional Case" Axioms 8, no. 4: 124. https://doi.org/10.3390/axioms8040124

APA Style

Eliseev, A., & Ratnikova, T. (2019). Regularized Solution of Singularly Perturbed Cauchy Problem in the Presence of Rational “Simple” Turning Point in Two-Dimensional Case. Axioms, 8(4), 124. https://doi.org/10.3390/axioms8040124

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

Article Metrics

Back to TopTop