Classical dynamics calculations are performed for O((3)P) + H2O((1)A1) collisions from 2 to 10 km s(-1) (4.1-101.3 kcal mol(-1)), focusing on product internal energies. Several methods are used to produce ro-vibrationally state-resolved product cross sections and to enforce zero-point maintenance from analysis of the classical trajectories. Two potential energy surfaces are used: (1) a recently developed set of global reactive surfaces for the three lowest triplet states which model OH formation, H elimination to make H + OOH, O-atom exchange, and collisional excitation and (2) a non-reactive surface used in past classical and quantum collision studies. Comparisons to these previous studies suggest that for H2O vibrational excitation, classical dynamics which include gaussian binning procedures and/or selected zero-point maintenance algorithms can produce results which approximate quantum scattering cross sections fairly well. Without these procedures, the classical cross sections can be many orders of magnitude greater than the quantum cross sections for exciting the bending vibration of H2O, especially near threshold. The classical cross section over-estimate is due to energy borrowing from stretching modes which dip below zero-point values. For results on the reactive surfaces, the present calculations show that at higher velocities there is an unusually large amount of product internal excitation. For OOH, where 40% of available collision energy goes into internal motion, the excited product vibrational and rotational energy distributions are relatively flat and values of the OOH rotational angular momentum exceed J = 100. Other product channel distributions show an exponential fall-off with energy consistent with an energy gap law. The present detailed distributions and cross sections can serve as a guide for future hyperthermal measurements of this system.