1. Introduction
Singularly perturbed problems with an unstable spectrum of the limit operator for ordinary partial differential equations and partial differential equations have been studied by many authors [
1,
2,
3,
4,
5,
6,
7,
8]. The most difficult of them are problems with point instability, namely turning points.
For the first time, these problems arose, in particular, in quantum mechanics. One of the first methods of solution was the WKB method, a method of semiclassical calculation, which was developed in 1926 by G. Wentzel, H.A. Kramers, and L. Brillouin. At the same time, H. Jeffreys generalized the method of approximate solution of linear differential equations of the second order, including the solution of the Schrodinger equation. Methods for solving problems with spectral features were and are being developed by: the school of V.P. Maslov, the school of A.B. Vasilieva—V.F. Butuzov—N.N. Nefedov, and the school of S.A. Lomov, among others.
The regularization method defines three groups of turning points.
1. “Simple” turning point: The eigenvalues of the limit operator are isolated from each other, and one eigenvalue at separate points of t vanishes.
2. “Weak” turning point: At least one pair of eigenvalues intersect at separate points of t, but the limit operator preserves the diagonal structure up to the intersection points. The eigenvector basis remains smooth in t.
3. “Strong” turning point: At least one pair of eigenvalues intersect at separate points of t, but the limit operator changes the diagonal structure to Jordan at the intersection points. The basis at the intersection points loses its smoothness in t.
Classic turning points are of the third type. A feature of the problem with a “simple” turning point presented in the article is the pointwise irreversibility of the limit operator .
In the present paper, using the Lomov’s regularization method [
1], a regularized asymptotic solution of the singularly perturbed Cauchy problem on the entire segment
for a parabolic equation is constructed in the presence of a “simple” rational turning point of the limit operator.
The point
for the singularly perturbed Cauchy problem is singular in the sense that the classical existence theorems for the solution of the Cauchy problem do not hold at this point. Therefore, in the solution of singularly perturbed problems, essentially special singularities arise that describe the irregular dependence of the solution on
[
9].
One of the advantages of the regularization method is that it makes it possible to construct a global asymptotic solution over the entire domain of integration, and, under certain conditions on the coefficients of the equation, it gives an exact solution.
The fractional turning point of order
for an ordinary differential equation of the first order was considered by the method of boundary functions by [
6].
In this article, we continue on the “simple” rational turning point work carried out in [
10,
11].
2. Formulation of the Problem and Construction of an Asymptotic Solution
Consider the Cauchy problem
where
is a function depending on variables
x and
t and real parameters
,
and
.
Let the following conditions be satisfied:
These conditions ensure the existence and uniqueness of the bounded solution and the possibility of constructing the asymptotic series of the problem (
1) [
12].
Singularly perturbed problems arise when the domain of definition of the initial operator depending on at does not coincide with the domain of definition of the limit operator at .
Under the stability condition for the spectrum of the limiting operator, essentially singular singularities are described using exponentials of the form , , , where are the eigenvalues of the limit operator and are smooth (in general, complex) functions of a real variable t.
If the stability conditions for the limit operator are violated for at least one point of the spectrum of the limit operator, then new singularities arise in the solution of the inhomogeneous equation. When studying problems with a “simple” turning point, we are faced with a problem when the range of values of the original operator does not coincide with the range of values of the limit operator.
Special singularities of the problem (
1) have the form:
We look for a solution
in the form [
1,
9]:
We denote differentiation with respect to
t by “
” and differentiation with respect to
x by “
” or
. Substituting (
3) into (
1) we get the equation
By identifying the coefficients in the linear combination of
,
and 1, we obtain the system
We look for a solution (
4) in the form of power series in
:
Substituting (
5) into (
4) we get a series of iterative tasks:
To solve iterative problems, the solvability theorem is used (see
Section 4).
Consider the system (
6) with
:
From (
7) it follows that
The functions
are found from the solvability condition for the system (
6) for
:
Let us expand
by Maclaurin’s formula in
t:
From the condition of solvability of the equation for
at
, it follows that, if
,
,
, then
where
is the whole part of
; if
,
,
, then
As a result, we get the solution at the
iteration step:
Considering (
10) we can write the system (
9) as:
The right-hand side of the equation with respect to
of the system (
12), due to the choice of
, satisfies the solvability theorem (see
Section 4):
(
s is an integer). Notice, that
.
Solving (
12) we get
where
are determined from (
10).
at the iteration step
are unknown. The functions
are found from the condition of solvability for
of the iterative problem for
:
where
is the fractional part of
.
To determine
, we subordinate the right side of the equation to the conditions of point solvability. To do this, we expand
by Maclaurin’s formula in
t:
where
,
,
. Then,
Using this scheme, you can find a solution at any iteration step.
For
, we get
where
are determined from (
14) and (
17).
Thus, the main term of the asymptotics of the solution is written in the form
3. Remainder Estimate
Let
iterative problems be solved. Then, the solution to the Cauchy problem can be represented in the form
Substituting (
19) into (
1) and taking into account that
are solutions to iterative problems, we obtain the Cauchy problem for determining the remainder
:
where
The remainder estimate is based on the maximum principle for parabolic problems [
12]. This principle is used in the form of generality that we need to estimate the remainder. The classical solution to the problem (
1) is a function
, continuous in
, having continuous
,
,
in
, and satisfies Equation (
1) and the initial conditions at
at all points of
.
Theorem 1 (remainder estimate). Let the conditions be satisfied:
- (1)
Conditions (1)–(5) of the Cauchy problem (1); - (2)
, ;
- (3)
, , ; and
- (4)
, .
Then, , .
Proof of Theorem 1. We denote . We prove the theorem in two stages.
Stage 1. Consider the homogeneous Cauchy problem in the domain
:
In the proof of this theorem, ideas from [
12] are used to estimate the solution. We introduce the function
. Then,
(1) For
Take
, then
.
(2) For
Take
, then
.
Besides,
Hence, by the maximum theorem in a bounded domain, we have
in
, i.e.
Letting
, we get
.
Stage 2. Consider the inhomogeneous Cauchy problem
We introduce the function
Then,
From the result of Stage 1, it follows that
in
, i.e.
. Consequently,
,
. We write the remainder
Then,
,
. ▹ □
4. Appendix
Lemma 1 (on the solvability of iterative problems)
. Let the equation be givenand conditions are metThen, Equation (23) is solvable in the class of smooth functions if and only if Proof of Lemma 1. We expand
by Maclaurin’s formula in
t:
Then, Equation (
23) has the form
Necessity, Let Equation (
23) have a solution. For
, we have
:
Dividing the equation by
t if
:
For
, we get
.
Continuing this process to step , we get .
Adequacy: Let be
. Then, Equation (
23) has the form:
hence the decision
. □
5. Conclusions
The novelty of the article lies in the construction by the regularization method of an asymptotic solution of a singularly perturbed Cauchy problem for a parabolic equation in the presence of a “simple” rational turning point of the limit operator. The nature of this “simple” turning point affects the structure of the functions describing the singular dependence of the solution on the parameter . The asymptotic expansion of the constructed solution is justified using the maximum principle. This approach can be applied both in the study of applied problems containing turning points and in the construction of numerical algorithms for solving problems with spectral features.