Nonlinear method of angular boundary functions for singularly perturbed parabolic 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. $$
Research is carried out under the assumption that at the corner points $(k,0)$ of the rectangle $\Omega$, where $k=0$ or $1$, the function $F(u)=F(u,k,0,0)$ is cubic and has the form
$$ F(u)=(u-\alpha(k))(u-\beta(k))(u-\bar u_0(k)), \text{ where } \alpha(k)\leq\beta(k)<\bar u_0(k). $$

The nonlinear method of angular boundary functions is used, which combines the (linear) method of angular boundary functions, the method of upper and lower solutions (barriers), and the method of differential inequalities. Under the condition $\varphi(k)>\bar u_0(k)$, a complete asymptotic expansion of the solution for $\varepsilon\rightarrow 0$ is constructed and its uniformity in a closed rectangle is substantiated.
Previously, the following cases of cubic nonlinearities were considered:
$$ F(u)=u^3-\bar u^3_0, \text{ where } \bar u_0=\bar u_0(k)>0, $$
under the assumption that the boundary value $\varphi( k)>\bar u_0(k)$, as well as the case
$$ F(u)=u^3-\bar u^3_0, \text{ where } \bar u_0=\bar u_0(k)< 0, $$
under the assumption that the boundary value $\varphi(k)$ is contained in the interval
$$ \bar u_0<\varphi(k)<\frac{\bar u_0}{2}< 0. $$