Gradient estimates for Stokes and Navier-Stokes systems with piecewise DMO coefficients

We study stationary Stokes systems in divergence form with piecewise Dini mean oscillation coefficients and data in a bounded domain containing a finite number of subdomains with $C^{1,\mathrm{D}\mathrm{i}\mathrm{n}\mathrm{i}}$ boundaries. We prove that if $(u, p)$isaweaksolutionofthesystem,then$(Du, p)$ is bounded and piecewise continuous. The corresponding results for stationary Navier-Stokes systems are also established, from which the Lipschitz regularity of the stationary $H^1$-weak solution in dimensions $d=2,3,4$ is obtained.