“…It is well known that the research of global regularity estimates in various function spaces for (non-)linear elliptic equations (or systems) in non-smooth domains is one of the most interesting and important topics in partial differential equations (see, for instance, [1,5,12,16,21,23,25,28,35,55] for the linear case, and [2,8,9,17,18,44,45,63,66] for the non-linear case). Moreover, the global regularity estimates for solutions to elliptic boundary problems depend not only on the structure of equations, the integrability of the right-hand side datum, and the properties of coefficients appearing in equations, but also on the smooth property or the geometric property of the boundary of domains (see, for instance, [1,5,16,24,25,28,36,41,44,62]).…”