For almost two decades the properties of 'dwarf' galaxies have challenged the cold dark matter (CDM) model of galaxy formation. Most observed dwarf galaxies consist of a rotating stellar disk embedded in a massive dark-matter halo with a near-constant-density core. Models based on the dominance of CDM, however, invariably form galaxies with dense spheroidal stellar bulges and steep central dark-matter profiles, because low-angular-momentum baryons and dark matter sink to the centres of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present-day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different dark-matter particle candidate. Here we report hydrodynamical simulations (in a framework assuming the presence of CDM and a cosmological constant) in which the inhomogeneous interstellar medium is resolved. Strong outflows from supernovae remove low-angular-momentum gas, which inhibits the formation of bulges and decreases the dark-matter density to less than half of what it would otherwise be within the central kiloparsec. The analogues of dwarf galaxies-bulgeless and with shallow central dark-matter profiles-arise naturally in these simulations.

We present an analysis of star formation and feedback recipes appropriate for galactic smoothed particle hydrodynamics simulations. Using an isolated Milky Way-like galaxy, we constrain these recipes based on well-established observational results. Our star formation recipe is based on that of Katz with the additional inclusion of physically motivated supernova feedback recipes. We propose a new feedback recipe in which Type II supernovae are modelled using an analytical treatment of blastwaves. With this feedback mechanism and a tuning of other star formation parameters, the star formation in our isolated Milky Way-like galaxy follows the slope and normalization of the observed Schmidt law. In addition, we reproduce the lowdensity cut-off and filamentary structure of star formation observed in disc galaxies. Our final recipe will enable better comparison of N-body simulations with observations.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.Comment: 15 pages, 13 figures, to be submitted to MNRAS. For high-resolution figures, please see http://www-theorie.physik.unizh.ch/~agertz

We used fully cosmological, high‐resolution N‐body + smooth particle hydrodynamic (SPH) simulations to follow the formation of disc galaxies with rotational velocities between 135 and 270 km s−1 in a Λ cold dark matter (CDM) universe. The simulations include gas cooling, star formation, the effects of a uniform ultraviolet (UV) background and a physically motivated description of feedback from supernovae (SNe). The host dark matter haloes have a spin and last major merger redshift typical of galaxy‐sized haloes as measured in recent large‐scale N‐body simulations. The simulated galaxies form rotationally supported discs with realistic exponential scalelengths and fall on both the I band and baryonic Tully–Fisher relations. An extended stellar disc forms inside the Milky Way (MW)‐sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a MW‐sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy's stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.

The key algorithms and features of the Gasoline code for parallel hydrodynamics with self-gravity are described. Gasoline is an extension of the efficient Pkdgrav parallel N -body code using smoothed particle hydrodynamics. Accuracy measurements, performance analysis and tests of the code are presented. Recent successful Gasoline applications are summarized. These cover a diverse set of areas in astrophysics including galaxy clusters, galaxy formation and gas-giant planets. Future directions for gasdynamical simulations in astrophysics and code development strategies for tackling cutting edge problems are discussed.