This simulation shows the first 20 milliseconds in the life of a neutron star which is formed in a Type II supernova. After an initial collapse phase, the neutron star becomes unstable to convection. The resulting convective motions destroy the spherical symmetry of the star and rapidly mix the inner regions. In addition, the neutrino flux from the neutron star will be non-spherical and will be significantly enhanced by the convective motions. This may have major implications for the Type II supernova mechanism. The calculation was performed using the Piecewise-Parabolic Method for hydrodynamics. The computational grid contained 300 zones in radius and 200 zones in angle. The inner 200 zones in radius were uniformly spaced, ranging from the inner boundary at 25 km to 175 km. The outer 100 zones were non-uniformly spaced and stretched to 2000 km. Only the inner 200 zones are plotted. The inner boundary was treated as a hard sphere. At the outer boundary, zero gradients for all the variables were assumed. Periodic boundary conditions were used along the sides of the grid. The following sequence shows the mixing of composition which results from the convective motions. The variable plotted is the electron fraction Ye, which ranges from 0.2 to 0.5. Note: Evolution of the electron fraction during the first 20 milliseconds in the life of a neutron star formed from a Type II supernova. Animator: Bruce Fryxell (NASA/GSFC), Pamela ONeil (NASA). Scientist: Bruce Fryxell (NASA/GSFC), Adam Burrows (University of Arizona).