Development of Methods of Asymptotic Analysis of Transition Layers in Reaction–Diffusion–Advection Equations: Theory and Applications

This work presents a review and analysis of modern asymptotic methods for analysis of singularly perturbed problems with interior and boundary layers. The central part of the work is a review of the work of the author and his colleagues and disciples. It highlights boundary and initial-boundary value problems for nonlinear elliptic and parabolic partial differential equations, as well as periodic parabolic problems, which are widely used in applications and are called reaction–diffusion and reaction–diffusion–advection equations. These problems can be interpreted as models in chemical kinetics, synergetics, astrophysics, biology, and other fields. The solutions of these problems often have both narrow boundary regions of rapid change and inner layers of various types (contrasting structures, moving interior layers: fronts), which leads to the need to develop new asymptotic methods in order to study them both formally and rigorously. A general scheme for a rigorous study of contrast structures in singularly perturbed problems for partial differential equations, based on the use of the asymptotic method of differential inequalities, is presented and illustrated on relevant problems. The main achievements of this line of studies of partial differential equations are reflected, and the key directions of its development are indicated.


INTRODUCTION
Presently, one of the topical problems in the singular perturbation theory is the study of nonlinear singularly perturbed partial differential equations, the solutions of which have boundary and interior layers. Such equations are of great interest both in the qualitative theory of differential equations and in many applied problems. In particular, such equations arise in mathematical models of processes of the reaction-diffusion and reaction-diffusion-advection type in chemical kinetics, synergetics, astrophysics, biology, and other fields, where the processes under study are characterized by narrow boundary regions of rapid changes in process parameters or by sharp interior layers of various types (contrast structures), stationary or moving (fronts). For an adequate mathematical description of such processes, it is necessary to develop new asymptotic methods for studying nonlinear singularly perturbed problems, the studies of which (initial value problems with small parameters at higher derivatives) began quite a long time ago. The direction of research presented in this work is the development of the deep ideas of the works of A.N. Tikhonov (see [1][2][3]), who laid the foundations of the modern singular perturbation theory. These works, which combined the singular perturbation theory with the qualitative theory of differential equations and Lyapunov's stability theory and were continued and developed in the works of A.B. Vasilieva, V.F. Butuzov, and their students and followers (see the book by Vasilieva and Butuzov [4]), to this day form the language of the singular perturbation theory, despite the significantly complicated problems under consideration.
Various areas of the theory of nonlinear singularly perturbed problems are being intensively developed both in our country and abroad. This is evidenced by a number of international conferences held in recent years and devoted to the singular perturbation theory, as well as conferences devoted to interior layers. In our country, these studies are carried out by the school of V.P. Maslov and S.Yu. Dobrokhotov; disciples of S.A. Lomov, S.A. Kashchenko, S.D. Glyzin, and many others. These works have been and are being actively pursued in the USA, European countries, Japan, and China (P. Fife, B. McLeod, P. Bates, N. Alikakos, Danielle Hilhorst, B. Fidler, L. Recke, K. Scheider, P. Szmolyan, H. Motano, H. Morita, K. Sakamoto, Vimura, Nishiura, and others). The results of these and other specialists are presented in the proceedings of the conference "Patterns of Dynamics 2016", held in Berlin, and in a special issue of the journal "Discrete and continuous dynamical systems," 37 (2), February 2017 (see [5]), dedicated to memory of Paul Fife.
The current state of the asymptotic analysis of singularly perturbed problems with boundary and inner layers can be found in the review [6]. A number of interesting results developing the studies of problems with interior and boundary layers were obtained in [7][8][9]. Currently, there is a significant increase in interest in singularly perturbed delay equations, which is associated with the consideration of more complex models important for applications (see, e.g., [10,11] and references therein). One of the most important applications of asymptotic analysis of singularly perturbed problems is the development of a theory of numerical methods for problems with transition layers. A significant contribution to the development of this direction was made in the works of A.M. Il'in, N.S. Bakhvalov, G.I. Shishkin, and N.V. Kopteva and Irish mathematicians M. Stynes, E. O'Riordan, and others. The results of these works can be found in the review [12]. A number of interesting results concerning theory and their application for the numerical solution of singularly perturbed problems with transition layers were obtained in [13][14][15][16][17][18]. Note that this actively developing area is based on the results of an asymptotic analysis of the works presented in this review.
In this paper, we present a development of the well-known method of boundary functions for constructing asymptotics for solutions of various classes of problems with interior layers. The main ideas are demonstrated on new classes of singularly perturbed problems, including problems on the motion of fronts and studies of stable stationary or periodic interior layers in multidimensional problems. These results are a further development of our studies of contrast structures, which were published in review articles [19][20][21][22].
Contrast structures are usually understood as solutions with interior layers of nonlinear singularly perturbed equations. Such solutions have been studied for a long time. The first results in this direction were obtained by A.B. Vasilieva back in the early 1970s for a two-point boundary value problem (1) It was proved that there exists a solution that, for small , is close to the root of the degenerate equation to the left of a certain point and close to the root to the right of . In this case, the point was not known in advance, but was determined in the course of constructing the asymptotics. Such solutions are now commonly referred to as step-type contrast structures. An example of nonlinearity for which solutions of the above type exist is . In this case, , , and is a simple root of the equation . These results were extended to some classes of quasilinear equations that are one-dimensional in the space variable. Methods for constructing the asymptotics of such classes of problems and the substantiation of the asymptotics constructed are contained in the book by A.B. Vasilieva and V.F. Butuzov (see [4]). A significant contribution to the development of asymptotic methods for analysis of nonlinear partial differential equations was made by the American mathematician P. Fife, communications with whom I recall with deep gratitude. He considered a multidimensional analogue of problem (1). The results of his work on boundary and interior layers for a multidimensional analogue of problem (1) (see [23,24]) played an essential role in the development of this line of studies. It should be noted that the existence of a boundary-layer solution and a solution with an interior transition layer in [23,24] was proved using methods similar to those used by A.B. Vasilieva for one-dimensional problems, while their realization in multidimensional problems is much more difficult. ε − , = , ∈ , , , ε = , , ε = .
If we consider a two-dimensional analogue of problem (1), where is a small parameter and is the Laplacian, the contrast structure can be defined as follows.
A step-type contrast structure is a solution of problem (2) that is close to two different solutions of the degenerate equation on opposite sides of some closed curve (the position of the curve is not known in advance, but is determined in the course of constructing the asymptotics). In the vicinity of the curve , a region of rapid change in the solution arises; the solution in this region is called an interior transition layer. In [25,26], a new method was proposed for proving the existence of solutions of singularly perturbed partial differential equations, called the asymptotic method of differential inequalities. This method turned out to be effective, and both the general scheme of this method was developed (see, e.g., [27][28][29]) and its application to other classes of singularly perturbed problems: initial-boundary value problems for parabolic equations in describing solutions with moving interior layers (fronts), periodic parabolic boundary value problems, boundary and initial-boundary value problems for some classes of integro-differential equations, and some classes of systems.
Another class of singularly perturbed problems discussed in this paper are problems with discontinuous nonlinearities (sources or advective coefficients), which simulate transition layers in the contact region of different media, as well as phase transitions. The paper presents new results concerning the existence of solutions of this class of problems, their asymptotics, and stability.
The study of the problems considered is based on the asymptotic method of differential inequalities. The main idea of this approach is the construction of lower and upper solutions of the problem using formal asymptotics. For all the problems considered, the existence of solutions is proved, the accuracy of the asymptotics is estimated, and their stability is analyzed.
The results presented in this paper form the basis for further development of methods for asymptotic analysis of new classes of nonlinear singularly perturbed problems for partial differential equations, which are widely used as mathematical models in many important applications. The development of asymptotic methods for studying complex nonlinear mathematical models with sharp transition layers will allowalong with solving important mathematical questions of existence, stability, and asymptotic approximation of solutions-to develop efficient numerical methods for solving these classes of problems, as well as new concepts of solving inverse problems based on asymptotic analysis, analyze conditions of destruction and sharpening of solutions, as well as the application of the results obtained in a number of applications: simulation of contrast structures in the studies of problems in biophysics, ecology, genetics, nonlinear wave theory, and other applications.
It should be noted that this review does not claim to be complete. This paper highlights the line of studies pursued by the scientific school to which the author belongs and the method to the development of which he is directly related.

BASIC IDEAS OF THE ASYMPTOTIC METHOD OF DIFFERENTIAL INEQUALITIES
In our works, a general scheme for a rigorous study of contrast structures in singularly perturbed problems for partial differential equations was developed based on the asymptotic method of differential inequalities. Below, for simplicity of presentation, we will explain this scheme on the example of problem (2), without specifying the requirements under which its implementation is possible. As noted in Introduction, this scheme consists of three main stages: constructing the formal asymptotics of the contrast structure, modifying this asymptotics to construct the upper and lower solutions of the problem, and using the lower and upper solutions constructed to prove the Lyapunov asymptotic stability of stationary or periodic solutions in the corresponding classes of problems. These steps will be outlined below in this section.

Construction of Formal Asymptotics
Step-type contrast structures in problem (2) and other similar problems are studied under the assumption that the nonlinearity is bistable, i.e., has two stable roots and an unstable root such ε Δ = , ,ε , ∈ ⊂ , , ε = , ∈ ∂ , for . The position of the transition curve between the roots and is usually determined by the condition (i.e., is defined as the projection of the curve of intersection between the solution and root onto the region ). Let be some closed sufficiently smooth curve lying in . In the neighborhood of , we introduce local coordinates and , where is the distance from a point from the neighborhood to the curve along the normal to the curve, taken with the "+" sign if the point lies inside the curve and with the " " sign if the point lies in the exterior of the curve, and is the coordinate of the point on the curve from which the specified normal issues. The equation of the curve is sought in the form The curve and functions are determined in the course of constructing the asymptotics.
The curve divides into regions and , internal and external to it, respectively.
In the region , we consider the problem We construct the asymptotics of problem (4) by the method of boundary functions (see [4,20]) in the form (5) where and denote the regular and boundary-layer (near ) parts of the asymptotics.
In the region , we consider the problem The asymptotics of problem (6) is constructed by analogy with (5) in the form (7) where denotes the regular part of the asymptotics and and are the boundary-layer parts of the asymptotics near the boundary and , respectively. The main terms of the regular parts of the asymptotics are chosen in the form , and the functions are used to describe the transition layer near .
By construction (see [4,20]), asymptotics (5) and (7) satisfy problems (4) and (6) with a residual up to , in particular, since all the terms in the sum, except the last one, are equal to zero by virtue of the equations for the coefficients of the asymptotics; in particular, and are nonlinear equations for determining the principal terms of asymptotics (5). The remaining coefficients of the asymptotics are deter- No. 12 2021 NEFEDOV mined from linear equations with the help of invertible operators: the operator , generating the regular part of the asymptotics, and the operator , generating the boundary-layer part of the asymptotics. A representation similar to (8) also holds for problem (6): (9) In this case, is a nonlinear equation for determining the leading term of the boundary-layer part of the asymptotics near the boundary and the subsequent boundary functions are determined from linear equations by virtue of the invertible operator . The rest of the terms in (9) are completely analogous to the corresponding terms in (8).
By virtue of the boundary conditions on the curve in problems (4) and (6), the asymptotics and are joined by continuity (the -matching condition is satisfied) on the curve . This curve (i.e., the curve ) and the coefficients in expansion (3) are determined from the condition of -matching of the asymptotics and on : (10) considered in the appropriate approximation in . In this case, is usually found from a nonlinear problem. To determine , the linear equations are obtained in the form (11) where is an invertible operator and is a function known at each step. This algorithm will be demonstrated in more detail with examples in the following sections.

Asymptotic Method of Differential Inequalities
We developed the asymptotic method of differential inequalities for those classes of singularly perturbed problems where the principle of comparison of elliptic boundary value problems, parabolic initialboundary value problems, and time-periodic parabolic boundary value problems, as well as problems for integro-differential equations, works. Basic comparison theorems for nonlinear equations, based on the principle of maximum or monotonicity of linearized operators, originate from the well-known S.A. Chaplygin's comparison theorems for initial-boundary value problems and Nagumo's theorems for boundary value problems and are one of the most important tools for studying nonlinear partial differential equations. These theorems can be found in a number of articles and monographs (see, e.g., [30] and references therein). A significant contribution to the development of the comparison principle was made by M.A. Krasnosel'skii, M.G. Krein, S.I. Pokhozhaev, H. Amann, P. Hess, and D.J. Sattinger (see [31][32][33][34][35][36][37]). These ideas were further developed in the works for special cases (see [38][39][40][41][42][43][44]).
The main idea of the asymptotic method of differential inequalities is the use of a modified formal asymptotics to construct comparison functions: lower and upper solutions; in this case, the main requirement for the implementation of this approach is the monotonicity of the asymptotics-generating operators linearized on the zeroth order. This approach was proposed in [25,26] and determined a strategy for proving the existence and estimating the accuracy of the asymptotics constructed. In [27], this method was extended to periodic parabolic boundary value problems. In that work, a fairly effective method was proposed for proving the asymptotic stability of periodic solutions, which was then transferred to the analysis of the asymptotic stability of stationary solutions of initial-boundary value problems for reaction-diffusion-type parabolic equations and then generalized to some more complex classes of reaction-diffusionadvection-type quasilinear equations (see [29,45]). In recent years, this approach has been extended to reaction-diffusion-advection problems with discontinuous nonlinearities and sources (see [46] and references therein). An essential development of the method was the results on the asymptotic analysis of boundary and interior layers in various problems for integro-differential equations (see [28], as well as [47] and references therein).
For the convenience of presentation, we recall the well-known concepts and facts concerning problem (2). Definition. Functions and are called lower and upper solutions of problem (2) if they satisfy the following conditions: It is known (see, e.g., [30]) that, if there exist ordered lower and upper solutions of problem (2), i.e., then problem (2) has a solution , satisfying the inequalities When constructing the formal asymptotics, the main requirement for the linear operators , , and is their invertibility, whereas, when constructing the lower and upper solutions, our requirements are more stringent: we require that inequalities , , and have positive solutions in the same classes of functions in which these operators are used to construct the corresponding expansions. Then, if we add and to the coefficients at the terms of the order of in the regular and boundarylayer parts of asymptotics (5) and (7) and add to from Eq. (3), we will obtain the upper solution . Similarly, by the addition of , , and , we construct the lower solution . Checking according to the standard scheme (see [26]) shows that and satisfy the relations (12) where in . In addition, by construction, It follows from (12) and (13) that and are ordered lower and upper solutions of problem (2) and, therefore, there exists a solution of this problem, satisfying the inequalities (14) Moreover, since, by construction, (15) it follows from (14) and (15) that the solution of problem (2) has an estimate Note that and satisfying (12) are called asymptotic lower and upper solutions of order .

Asymptotic Method of Differential Inequalities in the Analysis of Lyapunov Asymptotic Stability of Stationary and Periodic Solutions
As noted above, in [27], a fairly simple method was proposed for proving the Lyapunov asymptotic stability of periodic solutions of periodic parabolic problems and for determining the local region of stability and uniqueness of solutions with the asymptotics constructed. This method can also be used for analysis of the asymptotic stability of stationary solutions with boundary and interior layers of initial-boundary value problems for the reaction-diffusion and reaction-diffusion-advection equations. Note that the standard methods for analysis of stability, based on considering the corresponding spectral problems for equations linearized on the solution, are rather complicated and have not been significantly developed. The main difficulty of this approach is the presence of the so-called critical spectrum: eigenvalues that tend to zero as the small parameter tends to zero. This approach for partial differential equations, called the method of the singular limit of the eigenvalue problem, was first implemented in [48]. We also note No. 12 2021 NEFEDOV that this work has played an important role in the development of methods for constructing asymptotics in the so-called problems with balanced nonlinearity, as well as in problems with advection, where multidimensionality plays an essential role in the structure and stability properties of the solution. Let us illustrate the use of the method of differential inequalities on problem (2). The solution of problem (2) can be considered a stationary solution of the initial-boundary value problem (16) It can be shown that, if the upper and lower solutions of problem (2) are sufficiently accurate ( is large enough: see the previous paragraph), then the functions and (which are determined by the expressions (17) where is a sufficiently small positive number, and are the upper and lower solutions of problem (2), and is the solution of problem (2)) satisfy the corresponding differential inequalities for the parabolic initial-boundary value problem (16) (see, e.g., [30]). Then, for , a solution of problem (16) is , which satisfies the inequalities . These inequalities with representations (17) imply the asymptotic stability of the stationary solution of problem (16) with a stability region at least , . Moreover, is locally unique as a solution of problem (2) with the same uniqueness region.

SOME RELEVANT PROBLEMS WITH BOUNDARY AND INTERIOR LAYERS
In this section, we consider a number of topical problems of interest both for the singular perturbation theory and for applications. The consideration of these problems illustrates the main ideas and directions of development. One of the actively developed areas is the study of the existence and stability of stationary smooth solutions in one-dimensional and multidimensional reaction-diffusion problems with discontinuous sources and advective terms. These problems simulate transition layers at the interface between two media with different characteristics. Problems for the stationary reaction-diffusion equation on an interval and in a two-dimensional domain are considered. A specific feature of these problems is the discontinuity (of the first kind) of the reactive term in the one-dimensional case at some interior point and, in the two-dimensional case, on some smooth closed curve. This is a new previously unexplored class of problems. Methods for constructing asymptotic approximations of solutions are developed, the existence of smooth stationary solutions of reaction-diffusion equations with an interior transition layer in the vicinity of discontinuities in the one-dimensional and two-dimensional cases is proved, and conditions are obtained under which these solutions are locally unique and asymptotically Lyapunov-stable. In addition, the local region of attraction of stationary solutions is indicated. These studies are presented in [46,49,50] (see also the references therein).
Another important area of studies are new classes of problems: problems with modular nonlinearities. These are new classes of mathematical problems with discontinuous advection coefficients and/or discontinuous nonlinear source terms. The main feature of these problems is that discontinuities occur at some values of the sought-for function. These problems simulate phase transitions and are initiated by the studies in nonlinear wave theory (see [51][52][53]). In [54,55], initial-boundary value problems for Burgers-type equations were studied. The existence of new classes of solutions with a moving interior layer (front) was proved, conditions for a new effect in this class of problems-destruction (sharpening) of the front-were formulated, and the time for which the front moves from the initial point to the point of destruction was estimated. These works continued the study of problems for reaction-diffusion and reaction-diffusionadvection-type equations with moving interior layers (fronts). In [47,56] an asymptotic approximation of the fronts was obtained in one-dimensional parabolic and integro-parabolic reaction-diffusion-type problems. These results were then transferred to quasilinear reaction-diffusion-advection-type equations with a nonlinear advection coefficient in [57][58][59]. Note that the consideration of moving interior layers (fronts) in bounded volumes, taking into account the spatiotemporal inhomogeneity, differs signifi- cantly from the studies of solutions such as traveling waves and self-similar solutions, widely presented in the literature (see, e.g., [60][61][62] and references therein). Similar results were obtained in the study of periodic fronts in parabolic periodic boundary value problems, including Burgers-type equations. Various problems of this class are presented in [63][64][65][66][67][68]. In these works, periodic problems, multidimensional in the spatial variable, with an interior transition layer are considered and classes of new problems are singled out in which the multidimensionality leads to new previously unexplored conditions for the existence and stability of a solution with an interior transition layer. These are problems with a balance of nonlinearity, in which the results previously obtained for the reaction-diffusion problem are generalized to problems for the reaction-diffusion-advection equation. The relationships between advection and reaction are revealed, in which advection plays a significant role in the stability of the interior transition layers.
In [69][70][71], a new class of nonlinear singularly perturbed problems with singularly perturbed boundary conditions of the second and third kind was investigated. The existence of stationary and time-periodic boundary-layer-type solutions was proved, and an asymptotic approximation of such solutions over a small parameter was constructed. The set of boundary conditions under which such solutions exist was investigated, and the local uniqueness and asymptotic Lyapunov stability of such solutions was established. It was shown that, unlike the analogous Dirichlet problem, where such a solution is unique, the problem under consideration can have several such solutions with individual regions of stability and local uniqueness.
Among the actively developing areas, we should single out the studies of the existence and conditions of Lyapunov stability of stationary solutions of initial-boundary value problems and periodic solutions of systems of parabolic equations with fast sources of different velocities. The asymptotics of the solution was constructed, and the asymptotic stability of the periodic solution was investigated. Systems of equations with fast sources of different velocities in the case of different quasi-monotonicity of the right-hand sides (sources) have been studied. In [72,73], singularly perturbed systems of parabolic equations in onedimensional and two-dimensional spatial domains with periodic conditions in time were considered. The existence of a solution in the form of a periodic front was investigated. An algorithm for constructing an asymptotic approximation of the solution was presented, and a theorem of existence of a solution in the form of a periodic front, as well as its local uniqueness and asymptotic stability, was proved. A solution of this type can describe an abrupt change in the physical characteristics of some spatially inhomogeneous medium. Systems of this type are used to simulate transient processes in ecology, biophysics, chemical kinetics, semiconductor physics, and other fields. The aforementioned works are based on the development of classical methods for constructing asymptotics, proposed in the works of A.B. Vasilieva, V.F. Butuzov, and their disciples. In recent years, in the works of V.F. Butuzov, a new method was developed for constructing asymptotic approximations for the solution of singularly perturbed equations with multiple roots of degenerate equations (see [74][75][76][77][78][79][80][81][82][83] and references therein). In these works, it was shown that, in this case, the boundary and interior layers are multi-zone. The structure of the transition layers, in contrast to the cases of a simple layer, essentially depends on the type of boundary conditions. The study of the existence of solutions and their stability was based on the development of the asymptotic method of differential inequalities for this class of problems.
Among the important applications of asymptotic methods, we single out the development of the concept of asymptotic solution of inverse problems for nonlinear singularly perturbed reaction-diffusionadvection-type equations. A specific feature of mathematical models of this type is the presence in their solutions of stationary or moving narrow interior transition layers, which leads to instability in the numerical solution of direct problems and significantly complicates the further study of the corresponding inverse problems. The main idea of the concept of asymptotic solution of inverse problems is that asymptotic analysis makes it possible to reduce the original nonlinear singularly perturbed problem to a set of simpler problems that do not contain small parameters and have a lower spatial dimension or contain algebraic rather than differential equations. Thus, the use of asymptotic analysis makes it possible to establish simpler relationships between the observed parameters of the solution and the parameters of the inverse problem to be determined (coefficients in the equation, boundary and initial conditions, etc.), which can fundamentally change the approach to solving inverse problems. The presence of a small parameter in the original model provides in this case a sufficiently accurate qualitative and quantitative description of the solution. Recently, results have been obtained on the solution of inverse problems for a new class of nonlinear singularly perturbed Burgers-type equations with modular and quadratic advections with time-periodic coefficients, as well as some classes of initial-boundary value problems. It has been shown that, in such models, under certain conditions, the asymptotic approach makes it possible to reduce the coefficient inverse problems to linear algebraic equations relating, with a given accuracy, the observed position of the moving front with the input data of the model (coefficients in the equation and boundary condi-COMPUTATIONAL MATHEMATICS AND MATHEMATICAL PHYSICS Vol. 61 No. 12 2021 NEFEDOV tions). Thus, if there is a possibility to observe the trajectory of the front on a certain time interval, then the question of determining the unknown coefficient of the equation or the boundary regime reduces to a set of simple algebraic operations. An asymptotic solution was proposed for the problem of recovery of the source function in the Burgers equation with modular advection from known (exactly or approximately) information about the trajectory of the front on a given time interval (period). In a similar manner, the asymptotic solution of the boundary control problem, i.e., determining the boundary conditions under which the front moves according to a given law, was obtained. The approach developed was demonstrated by solving the inverse linear gain coefficient problem on the basis of asymptotic analysis and efficient numerical methods in the case when the linear gain is obtained by solving an ill-posed problem of differentiating the observed information. The high efficiency of the approach proposed in comparison with those previously developed for such classes of problems was demonstrated (see [84][85][86][87]).
In a number of works, a new approach has been developed for the analytico-numerical study of moving fronts in singularly perturbed reaction-diffusion-advection-type models. A method for generating a dynamically adapted mesh for efficient numerical solution of problems of this class was proposed. The method is based on a priori information about the motion and properties of the front, obtained as a result of a rigorous asymptotic analysis of a singularly perturbed parabolic problem. In particular, the essential parameters taken into account when constructing the mesh are estimates of the position of the transition layer and its width and structure. The analytico-numerical approach proposed can significantly save computing resources, reduce the computation time, and increase the stability of the computing process in comparison with the classical approaches. Examples demonstrating the main ideas and techniques of using the proposed approach were considered (see, e.g., [88] and references therein).
The abilities of the methods of asymptotic analysis in solving the problem of diagnosing the destruction in time and space of a solution of the Burgers equation with quadratic and modular nonlinearities and cubic amplification were demonstrated (see [54]). The result obtained by the methods of asymptotic analysis was refined using the numerical diagnostics of the destruction of the exact solution, based on calculating an a posteriori asymptotically exact estimate of the error obtained when calculating the approximate solution on a sequence of condensing meshes. Numerical experiments demonstrate that the asymptotic result gives a good estimate for the time of destruction of the solution and its localization in space for a fairly wide range of small parameters. Thus, it can be stated that a new method has been created and is being developed for the analytical study of the phenomena of destruction of solutions of singularly perturbed problems (for example, problems for a Burgers-type equation). An important difference between this method and those used by other researchers is that it not only can establish the fact of the destruction of the solution in time (as well as to obtain an asymptotic approximation of the destruction time), but also to determine the spatial localization of this phenomenon, which is an important result for many applied problems. Below, the ideas of the asymptotic method of differential inequalities are demonstrated in two classes of relevant problems. These sections are now being actively developed and transferred to new classes of problems of both theoretical and applied interest.

Multidimensional Interior Layers in a Problem with Balanced Nonlinearity
In this section, we will illustrate the general scheme of using the asymptotic method of differential inequalities on the example of the problem of a multidimensional stationary contrast structure. This problem arises in the mathematical modeling of various processes with interphase transitions.
Consider a spatially inhomogeneous reaction-diffusion equation with the Dirichlet boundary conditions: (18) and study the existence and Lyapunov stability of a stationary contrast structure. We will assume to be satisfied the following conditions: We introduce into consideration the function : and suppose that (A2) , .
Under condition (A2), the equation in problem (18) is called an equation with balanced nonlinearity.
(A3) Let the boundary condition belong to the influence domain of the root , i.e., the condition is satisfied for all .
We introduce the set of hypersurfaces The asymptotics of the solution with an interior transition layer localized near some hypersurface belonging to the set is constructed according to the scheme described in the previous section. In this case, conditions and ensure the solvability of the problems for determining the leading terms of asymptotic representations (5) and (7). The leading terms of the regular parts of these representations are determined from the degenerate equation, and those of the boundary-layer parts are determined from standard boundary-value problems for autonomous second-order ODEs, the solvability of which is established using analysis on the phase plane (see, e.g., [20]). The operators for . Therefore, they are invertible and, moreover, make it possible to construct lower and upper solutions of the required order by modifying the regular part of the asymptotics (see the previous section). The operator generating the boundary layer part of the asymptotic expansion, where is the stretched distance along the normal to the boundary , is well studied and, in particular, the corresponding inequality has a positive (exponentially decreasing) solution in the case of a monotonic boundary function , the addition of which to the boundary-layer part of the asymptotics makes it possible to achieve the ordering of the lower and upper solutions in the boundary-layer zone (see [20]).
We define a function for surfaces as follows: where is the mean curvature of the surface , ± ± ≡ ϕ , , > , ∈ .
No. 12 2021 NEFEDOV is the unit normal vector to the curve , directed to , is the solution of the boundary value problem describing the interior transition layer in the zeroth approximation: and is the stretched variable of the transition layer. It is known that this problem has a unique solution.
The -matching condition (10) is satisfied in the first approximation if we impose the following requirement (note that, in the case of balanced nonlinearity, condition (10) is satisfied in the zeroth approximation for any surface ): (A4) There is a surface such that As noted in the previous section, the -matching condition (10) in the following approximations leads to linear equations for determining the coefficients in the asymptotic representation (3) for the transition surface : where the operator is expressed via the original data as follows: where , , and denote, respectively, the Laplace-Beltrami, divergence, and gradient operators on the surface and ( ) are the principal curvatures of at a point . The function is defined as follows: The operator is self-adjoint, and all of its eigenvalues are real. Our next requirement is as follows.
(A5) The principal (maximum) eigenvalue of the operator satisfies the condition .
It ensures the unique solvability of the problems for determining , as well as the fulfillment of the condition for the existence of a positive solution of the inequality , which, as noted in Section 2.2, allows one to construct the lower and upper solutions, ordered in the zone of the transition layer. Thus, using the asymptotic method of differential inequalities, we obtain the following result for problem (18

. Periodic Fronts in the Problem with Modular Advection
Let us illustrate the development of an important line of research in problems with discontinuous coefficients and source terms by the following problem: (19) for , . The functions , , and are sufficiently smooth and -periodic in .
We define the regions , , and and give the following definition of the solution.

Definition. A function
is called a solution of problem (19) if it satisfies the equation in (19) in each and the boundary conditions.
Suppose that , , and are functions sufficiently smooth and -periodic in . Let and for .
For in (19), we obtain the so-called degenerate equation We consider Eq. (20) with one of the additional conditions of problem (19): Suppose to be satisfied the following conditions.
Problems (20), (21) and (20)  ). The formal asymptotics of the solution of each of these problems (23) is sought in the form (according to the scheme of Section 2.1; for details, see, e.g., [68]) (24) where and denote the regular and boundary-layer (near ) parts of asymptotics (19).
The functions are determined using the -matching conditions (25) in th order in . Using the standard procedure, we obtain a representation from which the coefficients of the regular part of asymptotics (24) are determined: Obviously, the leading terms of this representation, , are determined by condition and the higher orders terms are determined from the linear problems whose solution can be written out explicitly (it also implies the monotonicity of the operator generating the regular part of the asymptotics).
To construct the boundary layer part of asymptotics (24), we use the stretched variable , A.B. Vasilieva's representation of the nonlinearity of , and the representation of the differential operator , =: > ∈ .
For the coefficients of the boundary layer part of the asymptotics, we have the representation (26) From representation (26), we obtain equations , for which we use the standard condition on For and , we obtain the following problems: and (29) The solutions of problems (27), (28) and (29), (30) are written out explicitly: and . Obviously, is estimated as , , . The -matching condition (25) in the zeroth order, is satisfied, since it reduces to the problem which is solvable by condition . Note that this condition is analogous to the Rankine-Hugoniot condition for the class of problems under consideration.
From representation (26), one can obtain problems for and of the following orders, the solution of which can be written out explicitly. The -matching condition (25) in the first order in is reduced rather simply to the problem of finding a periodic solution of the equation where Φ 1 (t) are expressed via the known first-order asymptotic functions. The solution of this problem is guaranteed by condition . The functions of the next orders are determined from similar problems, and condition is the condition for the unique solvability and monotonicity of the operator generating the asymptotics of the position of the interior transition layer.
We define and : ,ε + εξ, ,ε + , ,ε + εξ, * * − , ε +ε ξ ,, ε, , ε +ε ξ , , * * , ,ε + ,ε , ,ε = . * We denote by partial sums of the th order of representations (24), and is replaced with . We also introduce the notation Consider the initial-boundary value problem corresponding to problem (19): The proof of the existence and stability theorem for a periodic solution is conducted using the asymptotic method of differential inequalities, which, in this problem, uses the following definition of the lower and upper solutions.
Definition. Functions and are called ordered lower and upper solutions of problem (19) for if they satisfy the following conditions: In the case when and are continuous in , but not smooth on some periodic curve , they must have an admissible jump of the derivative: Remark. It is known (see, e.g., [33]) that the existence of ordered lower and upper solutions implies the existence of a solution of problem (19), satisfying the inequalities In this problem, the lower and upper solutions are used with an admissible jump of the derivative on the modified curves of the formal asymptotics of the position of the interior layer and . In this case, there is the same result on existence (see, e.g., [36,40], as well as [39] and references therein). We have the following existence theorem for a solution of problem (19), its local uniqueness, and its asymptotic stability as a periodic solution of initial-boundary value problem (31 , ∈ := , : − < < , ∈ , − , ,ε = , , ,ε = ∈ , where is some positive constant independent of . This solution is asymptotically stable as a periodic solution of the initial-boundary value problem (31), with a stability region at least and locally unique in this region as a solution of problem (19).
To prove the existence of a solution and to estimate the accuracy of the asymptotic approximation of the solution of problem (19), we use the lower and upper solutions and constructed according to the scheme in Section 2.2. In this case, the modified asymptotics of the position of the transition layer, , has the form In these modifications, used to construct the upper and lower solutions, the shift is determined from the periodic problem introduced in condition : where is a periodic positive function. Therefore, due to the monotonicity of the operator of this problem. This ensures the necessary jump in the derivatives of and , as well as their ordering in the vicinity of the transition layer. The conditions of the definition are checked by the substitution into problem (19). The proof of the asymptotic stability of the periodic solution of problem (31) follows from the fact that the functions and defined at by the expressions where is a periodic solution of problem (19), and are the upper and lower solutions of problem (19), and and are the subdomains of lying to the right and left of the curves , are upper and lower solutions of problem (31). A similar proof scheme is used in other periodic problems cited above.

CONCLUSIONS
This review reflects the development of methods for asymptotic analysis of new classes of nonlinear singularly perturbed problems for parabolic and elliptic equations, including the so-called reaction-diffusion-advection equations with boundary and interior transition layers. These studies include both the development of methods for constructing formal asymptotic approximations (satisfying the problem in the residual) and a rigorous substantiation: the proof of the existence of a solution, an estimate of the accuracy of the asymptotic approximation constructed, and an analysis of the Lyapunov stability of the solutions. The novelty of each class of problems can be determined by new differential operators generating asymptotics. The rigorous mathematical results presented in the article are based on the use of the asymptotic method of differential inequalities. The essence of this method is to use asymptotics to construct comparison functions: lower and upper solutions. In this case, the properties of operators generating asymptotics-invertibility and monotonicity-are revealed. This method determines the research strategy. The transition to new classes of problems, including the generalization of a number of results to the multidimensional case, as well as to some classes of systems important for applications, requires the development of both methods for constructing asymptotics and the study of the invertibility of new classes of operators generating asymptotics, and the study of their monotonicity. In a number of cases, this also required the use and modification of basic comparison theorems that establish conditions for the existence of the corresponding classes of solutions.
The results of asymptotic analysis obtained are widely used and form a basis both for the development of the theory of numerical methods and for the creation of efficient numerical algorithms for the study of problems with transition layers. In the course of work in this direction, a new method for the analytical , ,ε NEFEDOV study of the phenomena of destruction of solutions was created and is being developed, the field of which is the class of singularly perturbed problems (e.g., problems for a Burgers-type equation). It is known that the dynamics of solutions of initial-boundary value problems for parabolic equations depends on the set of unstable stationary solutions. It is of interest to develop methods for proving the existence and finding instability conditions in the classes of problems presented above. Advances in this direction are presented in [89].
Among the important applications, we note the development of the concept of asymptotic solution of inverse problems for nonlinear singularly perturbed equations of the reaction-diffusion-advection type. It is known that the presence of stationary or moving narrow interior transition layers in their solutions leads to instability in the numerical solution of direct problems and significantly complicates the further study of the corresponding inverse problems. The main idea of the concept of asymptotic solution of inverse problems is that asymptotic analysis allows one to reduce the original nonlinear singularly perturbed problem to a set of simpler problems that do not contain small parameters and have a lower spatial dimension or contain algebraic rather than differential equations. Thus, the use of asymptotic analysis makes it possible to establish simpler relationships between the observed parameters of the solution and the parameters of the inverse problem to be determined (coefficients in the equation, boundary and initial conditions, etc.), which can fundamentally change the approach to solving inverse problems.
The results of asymptotic analysis obtained, the development of numerical methods and methods for solving inverse problems were applied to specific mathematical models of nonlinear wave theory, biology, magnetohydrodynamics, and urbo-ecological models. Among the important achievements, we note the creation of a model for the development of urban ecosystems, according to which the urban ecosystem can be regarded as an active environment with interacting natural and anthropogenic subsystems.