Irregular double-phase evolution problem: existence and global regularity

We investigate the homogeneous Dirichlet problem for the irregular double-phase evolution equation \[ u_t-\operatorname{div} \left( a(z)|\nabla u|^{p(z)-2} \nabla u + b(z)|\nabla u|^{q(z)-2} \nabla u\right)=f(z),\quad z=(x,t)\in Q_T:=Ω\times (0,T), \] where $Ω\subset \mathbb{R}^N$, $N \geq 2$ is a bounded domain, $T>0$, The non-differentiable coefficients $a(z)$, $b(z)$, the free term $f$, and the variable exponents $p$, $q$ are given functions. The coefficients $a$ and $b$ are nonnegative, bounded, satisfy the inequality \[ a(z)+b(z)\geq α\quad \text{in} \ Q_T, \quad \text{and} \quad |\nabla a|, |\nabla b|, a_t, b_t \in L^d(Q_T) \] for some constant $α>0$, and with $d>2$ depending on $\sup p(z)$, $\sup q(z)$, $N$, and the regularity of initial data $u(x,0)$. The free term $f$ and initial data $u(x,0)$ satisfy \[ f\in L^σ(Q_T) \ \text{with} \ σ>2 \quad \text{and} \quad |\nabla u(x,0)|\in L^{r}(Ω) \ \text{with} \ r\geq \max \bigg\{2,\sup_{Q_T}p(z),\sup_{Q_T}q(z)\bigg\}. \] The variable exponents $p,q \in C^{0,1}(\overline{Q}_T)$ satisfy the balance condition \[ \frac{2N}{N+2} < p(z), q(z)< +\infty \ \text{in} \ \overline Q_T \quad \text{and} \quad \max\limits_{\overline Q_T}|p(z)-q(z)|< \dfrac{2}{N+2}. \] Under the above assumptions, we establish the existence of a solution, which is obtained as the limit of classical solutions to a family of regularized problems and preserves initial temporal integrability: \[ |\nabla u(\cdot, t)| \in L^r(Ω) \ \text{for a.e.} \ t \in (0,T), \] gains global higher integrability: \[ |\nabla u|^{\min\{p(z), q(z)\} + s +r} \in L^1(Q_T) \ \text{for any} \ s \in \left(0, \frac{4}{N+2}\right), \] and attains second-order regularity: \[ a(z) |\nabla u|^{\frac{p+r-2}{2}}+b(z) |\nabla u|^{\frac{q+r-2}{2}}\in L^2(0,T;W^{1,2}(Ω)). \]