We present Direct Numerical Simulations of the transport of heat and heavy elements across a double-diffusive interface or a double-diffusive staircase, in conditions that are close to those one may expect to find near the boundary between the heavy-element rich core and the hydrogen-helium envelope of giant planets such as Jupiter. We find that the non-dimensional ratio of the buoyancy flux associated with heavy element transport to the buoyancy flux associated with heat transport lies roughly between 0.5 and 1, which is much larger than previous estimates derived by analogy with geophysical double-diffusive convection. Using these results in combination with a core-erosion model proposed by Guillot et al. (2004), we find that the entire core of Jupiter would be eroded within less than 1Myr assuming that the core-envelope boundary is composed of a single interface. We also propose an alternative model that is more appropriate in the presence of a well-established double-diffusive staircase, and find that in this limit a large fraction of the core could be preserved. These findings are interesting in the context of Juno's recent results, but call for further modeling efforts to better understand the process of core erosion from first principles.

Oscillatory double-diffusive convection (ODDC, more traditionally called semiconvection) is a form of linear double-diffusive instability that occurs in fluids that are unstably stratified in temperature (Schwarzschild unstable), but stably stratified in chemical composition (Ledoux stable). This scenario is thought to be quite common in the interiors of stars and giant planets, and understanding the transport of heat and chemical species by ODDC is of great importance to stellar and planetary evolution models. Fluids unstable to ODDC have a tendency to form convective thermo-compositional layers which significantly enhance the fluxes of temperature and chemical composition compared with microscopic diffusion. Although a number of recent studies have focused on studying properties of both layered and non-layered ODDC, few have addressed how additional physical processes such as global rotation affect its dynamics. In this work we study first how rotation affects the linear stability properties of rotating ODDC. Using direct numerical simulations we then analyze the effect of rotation on properties of layered and non-layered ODDC, and study how the angle of the rotation axis with respect to the direction of gravity affects layering. We find that rotating systems can be broadly grouped into two categories, based on the strength of rotation. Qualitative behavior in the more weakly rotating group is similar to non-rotating ODDC, but strongly rotating systems become dominated by vortices that are invariant in the direction of the rotation vector and strongly influence transport. We find that whenever layers form, rotation always acts to reduce thermal and compositional transport.

Oscillatory double-diffusive convection (ODDC) (also known as semi- convection) refers to a type of double diffusive instability that occurs in regions of planetary and stellar interiors which have a destabilizing thermal stratification and a stabilizing mean molecular weight stratification. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and chemical species by ODDC. Rosenblum et al. (2011) first showed that ODDC can either spontaneously form layers, which significantly enhance the transport of heat and chemical species compared to mi- croscopic transport, or remain in a state dominated by large scale gravity waves, in which there is a more modest enhancement of the turbulent transport rates. Subsequent studies in this series have focused on identifying under what condi- tions layers form (Mirouh et al. 2012), and quantifying transport through layered systems (Wood et al. 2013). Here we proceed to characterize transport through systems that are unstable to the ODDC instability, but do not undergo spon- taneous layer formation. We measure the thermal and compositional fluxes in non-layered ODDC from both 2D and 3D numerical simulations and show that 3D simulations are well approximated by similar simulations in a 2D domain. We find that the turbulent mixing rate in this regime is weak and can, to a first level approximation, be neglected. We conclude by summarizing the findings of papers I through III into a single prescription for transport by ODDC.

Merging carbon-oxygen (CO) white dwarfs are a promising progenitor system for Type Ia supernovae (SN Ia), but the underlying physics and timing of the detonation are still debated. If an explosion occurs after the secondary star is fully disrupted, the exploding primary will expand into a dense CO medium that may still have a disk-like structure. This interaction will decelerate and distort the ejecta. Here we carry out multi-dimensional simulations of ``tamped" SN Ia models, using both particle and grid-based codes to study the merger and explosion dynamics, and a radiative transfer code to calculate synthetic spectra and light curves. We find that post-merger explosions exhibit an hourglass-shaped asymmetry, leading to strong variations in the light curves with viewing angle. The two most important factors affecting the outcome are the scale-height of the disk, which depends sensitively on the binary mass ratio, and the total ${}^{56}$Ni yield, which is governed by the central density of the remnant core. The synthetic broadband light curves rise and decline very slowly, and the spectra generally look peculiar, with weak features from intermediate mass elements but relatively strong carbon absorption. We also consider the effects of the viscous evolution of the remnant, and show that a longer time delay between merger and explosion probably leads to larger ${}^{56}$Ni yields and more symmetrical remnants. We discuss the relevance of this class of aspherical ``tamped" SN Ia for explaining the class of ``super-Chandrasekhar'' SN Ia.

Merging white dwarfs are a possible progenitor of Type Ia supernovae (SNe Ia). While it is not entirely clear if and when an explosion is triggered in such systems, numerical models suggest that a detonation might be initiated before the stars have coalesced to form a single compact object. Here we study such "peri-merger" detonations by means of numerical simulations, modeling the disruption and nucleosynthesis of the stars until the ejecta reach the coasting phase. Synthetic light curves and spectra are generated for comparison with observations. Three models are considered with primary masses 0.96 Msun, 1.06 Msun, and 1.20 Msun. Of these, the 0.96 Msun dwarf merging with an 0.81 Msun companion, with a Ni56 yield of 0.58 Msun, is the most promising candidate for reproducing common SNe Ia. The more massive mergers produce unusually luminous SNe Ia with peak luminosities approaching those attributed to "super-Chandrasekhar" mass SNe Ia. While the synthetic light curves and spectra of some of the models resemble observed SNe Ia, the significant asymmetry of the ejecta leads to large orientation effects. The peak bolometric luminosity varies by more than a factor of 2 with the viewing angle, and the velocities of the spectral absorption features are lower when observed from angles where the light curve is brightest. The largest orientation effects are seen in the ultraviolet, where the flux varies by more than an order of magnitude. Despite the large variation with viewing angle, the set of three models roughly obeys a width-luminosity relation, with the brighter light curves declining more slowly in the B-band. Spectral features due to unburned carbon from the secondary star are also seen in some cases.

Using 2D and 3D simulation, we study the "robustness" of the double detonation scenario for Type Ia supernovae, in which a detonation in the helium shell of a carbon-oxygen white dwarf induces a secondary detonation in the underlying core. We find that a helium detonation cannot easily descend into the core unless it commences (artificially) well above the hottest layer calculated for the helium shell in current presupernova models. Compressional waves induced by the sliding helium detonation, however, robustly generate hot spots which trigger a detonation in the core. Our simulations show that this is true even for non-axisymmetric initial conditions. If the helium is ignited at multiple points, the internal waves can pass through one another or be reflected, but this added complexity does not defeat the generation of the hot spot. The ignition of very low-mass helium shells depends on whether a thermonuclear runaway can simultaneously commence in a sufficiently large region.

The launching process of a magnetically driven outflow from an accretion disk is investigated in a local, shearing box model which allows a study of the feedback between accretion and angular momentum loss. The mass-flux instability found in previous linear analyses of this problem is recovered in a series of 2D (axisymmetric) simulations in the MRI-stable (high magnetic field strength) regime. At low field strengths that are still sufficient to suppress MRI, the instability develops on a short radial length scale and saturates at a modest amplitude. At high field strengths, a long-wavelength "clump" instability of large amplitude is observed, with growth times of a few orbits. As speculated before, the unstable connection between disk and outflow may be relevant for the time dependence observed in jet-producing disks. The success of the simulations is due in a large part to the implementation of an effective wave-transmitting upper boundary condition.

Aims: We study the differences between non-magnetic and magnetic regions in the flow and thermal structure of the upper solar photosphere. Methods: Radiative MHD simulations representing a quiet region and a plage region, respectively, which extend into the layers around the temperature minimum, are analyzed. Results: The flow structure in the upper photospheric layers of the two simulations is considerably different: the non-magnetic simulation is dominated by a pattern of moving shock fronts while the magnetic simulation shows vertically extended vortices associated with magnetic flux concentrations. Both kinds of structures induce substantial local heating. The resulting average temperature profiles are characterized by a steep rise above the temperature minimum due to shock heating in the non-magnetic case and by a flat photospheric temperature gradient mainly caused by Ohmic dissipation in the magnetic run. Conclusions: Shocks in the quiet Sun and vortices in the strongly magnetized regions represent the dominant flow structures in the layers around the temperature minimum. They are closely connected with dissipation processes providing localized heating.

We report on the occurrence of small-scale vortices in simulations of the convective solar surface. Using an eigenanalysis of the velocity gradient tensor, we find the subset of high vorticity regions in which the plasma is swirling. The swirling regions form an unsteady, tangled network of filaments in the turbulent downflow lanes. Near-surface vertical vortices are underdense and cause a local depression of the optical surface. They are potentially observable as bright points in the dark intergranular lanes. Vortex features typically exist for a few minutes, during which they are moved and twisted by the motion of the ambient plasma. The bigger vortices found in the simulations are possibly, but not necessarily, related to observations of granular-scale spiraling pathlines in "cork animations" or feature tracking.

We quantify possible differences between turbulent dynamo action in the Sun and the dynamo action studied in idealized simulations. For this purpose we compare Fourier-space shell-to-shell energy transfer rates of three incrementally more complex dynamo simulations: an incompressible, periodic simulation driven by random flow, a simulation of Boussinesq convection, and a simulation of fully compressible convection that includes physics relevant to the near-surface layers of the Sun. For each of the simulations studied, we find that the dynamo mechanism is universal in the kinematic regime because energy is transferred from the turbulent flow to the magnetic field from wavenumbers in the inertial range of the energy spectrum. The addition of physical effects relevant to the solar near-surface layers, including stratification, compressibility, partial ionization, and radiative energy transport, does not appear to affect the nature of the dynamo mechanism. The role of inertial-range shear stresses in magnetic field amplification is independent from outer-scale circumstances, including forcing and stratification. Although the shell-to-shell energy transfer functions have similar properties to those seen in mean-flow driven dynamos in each simulation studied, the saturated states of these simulations are not universal because the flow at the driving wavenumbers is a significant source of energy for the magnetic field.

We consider the conditions under which a rotating magnetic object can produce a magnetically powered outflow in an initially unmagnetized medium stratified under gravity. 3D MHD simulations are presented in which the footpoints of localized, arcade-shaped magnetic fields are put into rotation. It is shown how the effectiveness in producing a collimated magnetically powered outflow depends on the rotation rate, the strength and the geometry of the field. The flows produced by uniformly rotating, non-axisymmetric fields are found to consist mainly of buoyant plumes heated by dissipation of rotational energy. Collimated magnetically powered flows are formed if the field and the rotating surface are arranged such that a toroidal magnetic field is produced. This requires a differential rotation of the arcades' footpoints. Such jets are well-collimated; we follow their propagation through the stratified atmosphere over 100 times the source size. The magnetic field is tightly wound and its propagation is dominated by the development of non-axisymmetric instabilities. We observe a Poynting flux conversion efficiency of over 75% in the longest simulations. Applications to the collapsar model and protostellar jets are discussed.

A 3D simulation of a non-relativistic, magnetically driven jet propagating in a stratified atmosphere is presented, covering about three decades in distance and two decades in sideways expansion. The simulation captures the jet acceleration through the critical surfaces and the development of (kink-)instabilities driven by the free energy in the toroidal magnetic field component. The instabilities destroy the ordered helical structure of the magnetic field, dissipating the toroidal field energy on a length scale of about 2-15 times the Alfven distance. We compare the results with a 2.5D (axisymmetric) simulation, which does not become unstable. The acceleration of the flow is found to be quite similar in both cases, but the mechanisms of acceleration differ. In the 2.5D case approximately 20% of the Poynting flux remains in the flow, in the 3D case this fraction is largely dissipated internally. Half of the dissipated energy is available for light emission; the resulting radiation would produce structures resembling those seen in protostellar jets.

We have performed 2.5D and 3D simulations of conical jets driven by the rotation of an ordered, large-scale magnetic field in a stratified atmosphere. The simulations cover about three orders of magnitude in distance to capture the centrifugal acceleration as well as the evolution past the Alfven surface. We find that the jets develop kink instabilities, the characteristics of which depend on the velocity profile imposed at the base of the flow. The instabilities are especially pronounced with a rigid rotation profile, which induces a shearless magnetic field. The jet's expansion appears to be limiting the growth of Alfven mode instabilities.

We assess the importance of AGN outflows with respect to the metal enrichment of the intracluster medium (ICM) in galaxy clusters. We use combined N-body and hydrodynamic simulations, along with a semi-numerical galaxy formation and evolution model. Using assumptions based on observations, we attribute outflows of metal-rich gas initiated by AGN activity to a certain fraction of our model galaxies. The gas is added to the model ICM, where the evolution of the metallicity distribution is calculated by the hydrodynamic simulations. For the parameters describing the AGN content of clusters and their outflow properties, we use the observationally most favorable values. We find that AGNs have the potential to contribute significantly to the metal content of the ICM or even explain the complete abundance, which is typically ~0.5 Z_sun in core regions. Furthermore, the metals end up being inhomogeneously distributed, in accordance with observations.