Numerous formulations of finite volume schemes for the Euler and Navier-Stokes equations exist, but in the majority of cases they have been developed for structured and stationary meshes. In many applications, more flexible mesh geometries that can dynamically adjust to the problem at hand and move with the flow in a (quasi) Lagrangian fashion would, however, be highly desirable, as this can allow a significant reduction of advection errors and an accurate realization of curved and moving boundary conditions. Here we describe a novel formulation of viscous continuum hydrodynamics that solves the equations of motion on a Voronoi mesh created by a set of mesh-generating points. The points can move in an arbitrary manner, but the most natural motion is that given by the fluid velocity itself, such that the mesh dynamically adjusts to the flow. Owing to the mathematical properties of the Voronoi tessellation, pathological mesh-twisting effects are avoided. Our implementation considers the full Navier-Stokes equations and has been realized in the AREPO code both in 2D and 3D. We propose a new approach to compute accurate viscous fluxes for a dynamic Voronoi mesh, and use this to formulate a finite volume solver of the Navier-Stokes equations. Through a number of test problems, including circular Couette flow and flow past a cylindrical obstacle, we show that our new scheme combines good accuracy with geometric flexibility, and hence promises to be competitive with other highly refined Eulerian methods. This will in particular allow astrophysical applications of the AREPO code where physical viscosity is important, such as in the hot plasma in galaxy clusters, or for viscous accretion disk models.

Water-rich super-Earth exoplanets are expected to be common. We explore the effect of late giant impacts on the final bulk abundance of water in such planets. We present the results from smoothed particle hydrodynamics simulations of impacts between differentiated water(ice)-rock planets with masses between 0.5 and 5 M_Earth and projectile to target mass ratios from 1:1 to 1:4. We find that giant impacts between bodies of similar composition never decrease the bulk density of the target planet. If the commonly assumed maximum water fraction of 75wt% for bodies forming beyond the snow line is correct, giant impacts between similar composition bodies cannot serve as a mechanism for increasing the water fraction. Target planets either accrete materials in the same proportion, leaving the water fraction unchanged, or lose material from the water mantle, decreasing the water fraction. The criteria for catastrophic disruption of water-rock planets are similar to those found in previous work on super-Earths of terrestrial composition. Changes in bulk composition for giant impacts onto differentiated bodies of any composition (water-rock or rock-iron) are described by the same equations. These general laws can be incorporated into future N-body calculations of planet formation to track changes in composition from giant impacts.

Haumea, a rapidly rotating elongated dwarf planet (~ 1500 km in diameter), has two satellites and is associated with a "family" of several smaller Kuiper Belt objects (KBOs) in similar orbits. All members of the Haumea system share a water ice spectral feature that is distinct from all other KBOs. The relative velocities between the Haumea family members are too small to have formed by catastrophic disruption of a large precursor body, which is the process that formed families around much smaller asteroids in the Main Belt. Here we show that all of the unusual characteristics of the Haumea system are explained by a novel type of giant collision: a graze-and-merge impact between two comparably sized bodies. The grazing encounter imparted the high angular momentum that spun off fragments from the icy crust of the elongated merged body. The fragments became satellites and family members. Giant collision outcomes are extremely sensitive to the impact parameters. Compared to the Main Belt, the largest bodies in the Kuiper Belt are more massive and experience slower velocity collisions; hence, outcomes of giant collisions are dramatically different between the inner and outer solar system. The dwarf planets in the Kuiper Belt record an unexpectedly large number of giant collisions, requiring a special dynamical event at the end of solar system formation.

The detailed interior structure models of super-Earth planets show that there is degeneracy in the possible bulk compositions of a super-Earth at a given mass and radius, determined via radial velocity and transit measurements, respectively. In addition, the upper and lower envelopes in the mass--radius relationship, corresponding to pure ice planets and pure iron planets, respectively, are not astrophysically well motivated with regard to the physical processes involved in planet formation. Here we apply the results of numerical simulations of giant impacts to constrain the lower bound in the mass--radius diagram that could arise from collisional mantle stripping of differentiated rocky/iron planets. We provide a very conservative estimate for the minimum radius boundary for the entire mass range of large terrestrial planets. This envelope is a readily testable prediction for the population of planets to be discovered by the Kepler mission.

The final stage of planet formation is dominated by collisions between planetary embryos. The dynamics of this stage determine the orbital configuration and the mass and composition of planets in the system. In the solar system, late giant impacts have been proposed for Mercury, Earth, Mars, and Pluto. In the case of Mercury, this giant impact may have significantly altered the bulk composition of the planet. Here we present the results of smoothed particle hydrodynamics simulations of high-velocity (up to ~5 v_esc) collisions between 1 and 10 M_Earth planets of initially terrestrial composition to investigate the end stages of formation of extrasolar super-Earths. As found in previous simulations of collisions between smaller bodies, when collision energies exceed simple merging, giant impacts are divided into two regimes: (1) disruption and (2) hit-and-run (a grazing inelastic collision and projectile escape). Disruption occurs when the impact parameter is near zero, when the projectile mass is small compared to the target, or at extremely high velocities. In the disruption regime, we derive the criteria for catastrophic disruption (when half the total colliding mass is lost), the transition energy between accretion and erosion, and a scaling law for the change in bulk composition (iron-to-silicate ratio) resulting from collisional stripping of a mantle.