Singularly Perturbed Cauchy Problem for a Parabolic Equation with a Rational “Simple” Turning Point

: The aim of the research is to develop the regularization method. By Lomov’s regularization method, we constructed a uniform asymptotic solution of the singularly perturbed Cauchy problem for a parabolic equation in the case of violation of stability conditions of the limit-operator spectrum. The problem with a “simple” turning point is considered in the case, when the eigenvalue vanishes at t = 0 and has the form t m / n a ( t ) . The asymptotic convergence of the regularized series is proved.

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 t m/n .
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 [0, T] for a parabolic equation is constructed in the presence of a "simple" rational turning point of the limit operator.
The point ε = 0 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 1/2 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].

Formulation of the Problem and Construction of an Asymptotic Solution
Consider the Cauchy problem where u = u(x, t, ε) is a function depending on variables x and t and real parameters ε, ε > 0 and ε << 1. Let the following conditions be satisfied: , the function h(x, t) and all its derivatives are bounded , the function f (x) and all its derivatives are bounded to R; (4) a(t) ∈ C ∞ ([0, T]), a(t) > 0; and (5) m, n ∈ N, p = m + n − 1, m/n = r is a fractional number. (2) 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 ε = 0 does not coincide with the domain of definition of the limit operator at ε = 0.
Under the stability condition for the spectrum of the limiting operator, essentially singular singularities are described using exponentials of the form where λ i (t) are the eigenvalues of the limit operator and ϕ i (t) 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 u(x, t, ε) in the form [1,9]: We denote differentiation with respect to t by " · " and differentiation with respect to x by " " or ∂ ∂x . Substituting (3) into (1) we get the equation By identifying the coefficients in the linear combination of e −ϕ(t)/ε , σ i 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 k = −1: The functions Z i −1 (x, 0) are found from the solvability condition for the system (6) for k = 0: Let us expand h(x, t) by Maclaurin's formula in t: From the condition of solvability of the equation for W 0 at t = 0, it follows that, if i+1−n n = j, where m n is the whole part of m/n; if i+1−n n = j, i = 0, (p − 1), 0 ≤ j ≤ m n , then As a result, we get the solution at the (−1) iteration step: Considering (10) we can write the system (9) as: The right-hand side of the equation with respect to W 0 of the system (12), due to the choice of Z i −1 (x, 0), satisfies the solvability theorem (see Section 4): (s is an integer). Notice, that W 0 (x, 0) = 0.
Solving (12) we get where Z i −1 (x, 0) are determined from (10). Z i 0 (x, 0) at the iteration step k = 0 are unknown. The functions Z i 0 (x, 0) are found from the condition of solvability for t = 0 of the iterative problem for k = 1: where m n is the fractional part of m/n. To determine W 1 (x, t), we subordinate the right side of the equation to the conditions of point solvability. To do this, we expand h 1 (x, t) by Maclaurin's formula in t: where k = m n + m n , m n = s n , 1 ≤ s ≤ n − 1. Then, Using this scheme, you can find a solution at any iteration step. For u 0 (x, t, ε), we get where Z i 0 (x, t) are determined from (14) and (17). Thus, the main term of the asymptotics of the solution is written in the form

Remainder Estimate
Let (N + 1) 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 u k (x, t, ε) are solutions to iterative problems, we obtain the Cauchy problem for determining the remainder R N (x, t, ε): 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 R(x, t, ε), continuous in Q T = (−∞, +∞) × [0, T] × (0, ε], having continuous ∂R ∂t , ∂R ∂x , ∂ 2 R ∂x 2 in Q T , and satisfies Equation (1) and the initial conditions at t = 0 at all points of Q T . Theorem 1 (remainder estimate). Let the conditions be satisfied: (1) Conditions (1)-(5) of the Cauchy problem (1); Proof of Theorem 1. We denote t m/n a(t) = q(t). We prove the theorem in two stages.
Letting L → +∞, we get u ≥ 0. Stage 2. Consider the inhomogeneous Cauchy problem We introduce the function Then, From the result of Stage 1, it follows that We write the remainder

Appendix
Lemma 1 (on the solvability of iterative problems). Let the equation be given and conditions are met

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.