Optimal second-order regularity in the space variables is established for solutions to Cauchy-Dirichlet problems for nonlinear parabolic equations and systems of p-Laplacian type, with square-integrable right-hand sides and initial data in a Sobolev space. As a consequence, generalized solutions are shown to be strong solutions. Minimal regularity on the boundary of the domain is required, though the results are new even for smooth domains. In particular, they hold in arbitrary bounded convex domains.Here, p > 1, Ω is an open set in R n , n ≥ 1, with finite Lebesgue measure |Ω|, and T > 0. Moreover, Ω T = Ω × (0, T ), the functions f : Ω T → R N and ψ : Ω → R N , N ≥ 1 are given, and u : Ω → R N is the unknown. According to usage, ∇u stands for the gradient of u with respect to the space variables x ∈ Ω, and u t for its derivative in time t ∈ (0, T ).Mathematics Subject Classifications: 35K20, 35K65, 35B65.