We have conducted extensive lattice simulations to study the post-inflation dynamics of multifield models involving nonminimal couplings. We explore the parameter dependence of preheating in these models and describe the various timescales that control such nonlinear processes as energy transfer, re-scattering, and the approach to radiation-domination and thermalization. In the limit of large nonminimal couplings (ξ I ∼ 100), we find that efficient transfer of energy from the inflaton condensate to radiative degrees of freedom, emergence of a radiation-dominated equation of state, and the onset of thermalization each consistently occur within N reh 3 e-folds after the end of inflation, largely independent of the values of the other couplings in the models. The exception is the case of negative ellipticity, in which there is a misalignment between the dominant direction in field-space along which the system evolves and the larger of the nonminimal couplings ξ I. In those cases, the field-space-driven parametric resonance is effectively shut off. More generally, the competition between the scalar fields' potential and the field-space manifold structure can yield interesting phenomena such as two-stage resonances. Across many regions of parameter space, we find efficient rescattering between the distinct fields, leading to a partial memory loss of the shape of the initial fluctuation spectrum. Despite the explosive particle production, which can lead to a quick depletion of the background energy density, the nonlinear processes do not induce any super-horizon correlations after the end of inflation in these models, which keeps predictions for CMB observables unaffected by the late-time amplification of isocurvature fluctuations. Hence the excellent agreement between primordial observables and recent observations is preserved for this class of models, even when we consider post-inflation dynamics.