Nonlinear method of angular boundary functions in problems with cubic nonlinearities

In the rectangle $\Omega =\{(x,t) | 0<x<1, 0<t<T\}$ we consider an initial-boundary value problem for a singularly perturbed parabolic equation
$$ \varepsilon^2\left(a^2\frac{\partial^2 u}{\partial x^2}-\frac{\partial u}{\partial t}\right)=F(u,x,t,\varepsilon), (x,t)\in \Omega, $$

$$ u(x,0,\varepsilon)=\varphi(x), 0\le x\le 1, $$

$$ u(0,t,\varepsilon)=\psi_1(t), u(1,t,\varepsilon)=\psi_2(t), 0\le t\le T. $$
It is assumed that at the corner points of the rectangle the function $F$ with respect to the variable $u$ is cubic. To construct the asymptotics of the solution to the problem, the nonlinear method of angular boundary functions is used, which involves the following steps:
1) splitting the area into parts;
2) construction in each subdomain of lower and upper solutions of the problem;
3) continuous joining of the lower and upper solutions on the common boundaries of the subdomains;
4) subsequent smoothing of piecewise continuous lower and upper solutions.
In the present work, we succeeded in constructing barrier functions suitable for the entire region at once. The form of barrier functions is determined using boundary-layer functions that are solutions of ordinary differential equations, as well as taking into account the necessary properties of the desired solutions. As a result, a complete asymptotic expansion of the solution for $\varepsilon\rightarrow 0$ is constructed and its uniformity in a closed rectangle is justified.