We demonstrate that the critical magnetic Reynolds number $Rm_c$ for a turbulent non-helical dynamo in the low magnetic Prandtl number $Pm$ limit (i.e. $Pm = Rm/Re \ll 1$) can be significantly reduced if the flow is submitted to global rotation. Even for moderate rotation rates the required energy injection rate can be reduced by a factor more than $10^3$. This strong decrease of the onset is attributed to the reduction of the turbulent fluctuations that makes the flow to have a much larger cut-off length-scale compared to a non-rotating flow of the same Reynolds number. The dynamo thus behaves as if it is driven by laminar behaviour (i.e. high $Pm$ behaviour) even at high values of the Reynolds number (i.e. at low values of $Pm$). Our finding thus points into a new paradigm for the design of new liquid metal dynamo experiments.

At the heart of today's solar magnetic field evolution models lies the alpha dynamo description. In this work, we investigate the fate of alpha-dynamos as the magnetic Reynolds number $Rm$ is increased. Using Floquet theory, we are able to precisely quantify mean field effects like the alpha and beta effect (i) by rigorously distinguishing dynamo modes that involve large scale components from the ones that only involve small scales, and by (ii) providing a way to investigate arbitrary large scale separations with minimal computational cost. We apply this framework to helical and non-helical flows as well as to random flows with short correlation time. Our results determine that the alpha-description is valid for $Rm$ smaller than a critical value $Rm_c$ at which small scale dynamo instability starts. When $Rm$ is above $Rm_c$ the dynamo ceases to follow the mean field description and the growth rate of the large scale modes becomes independent of the scale separation while the energy in the large scale modes is inversely proportional to the square of the scale separation. The results in this second regime do not depend on the presence of helicity. Thus alpha-type modeling for solar and stellar models needs to be reevaluated and new directions for mean field modeling are proposed.

We consider the generation of magnetic field by the flow of a fluid for which the electrical conductivity is nonuniform. A new amplification mechanism is found which leads to dynamo action for flows much simpler than those considered so far. In particular, the fluctuations of the electrical conductivity provide a way to bypass anti-dynamo theorems. For astrophysical objects, we show through three-dimensional global numerical simulations that the temperature-driven fluctuations of the electrical conductivity can amplify an otherwise decaying large scale equatorial dipolar field. This effect could play a role for the generation of the unusually tilted magnetic field of the iced giants Neptune and Uranus.

Magnetohydrodynamic turbulent flows driven by random mechanical and electromagnetic external forces of zero helicities are investigated by means of direct numerical simulations. It is shown that despite the absence of helicities in the forcing, the system is attracted to self-organized helical states that exhibit laminar behaviour despite the large value of the Reynolds numbers examined. We demonstrate that the correlation time of the external forces is controlling the time spent on these states, i.e. for short correlation times the system remains in the turbulent state while as the correlation time is increased the system spends more and more time in the self-organised states. As a result, time averaged statistics can significantly be affected by the time spent on these states. These results have important theoretical implications for the understanding of the suppression of non-linearities in plasma fusion devises as well as in astrophysical observations.

We demonstrate that special correlations in the initial conditions of freely evolving, homogeneous magnetohydrodynamic (MHD) turbulence can lead to the formation of enormous current sheets. These coherent structures are observed at the peak of the energy dissipation rate and are the carriers of long-range correlations despite all the non-linear interactions during the formation of turbulence. Even though the largest scale separation has been achieved at this point, these structures are coherent with a size that spans our computational domain dominating the scaling of the energy spectrum, which follows a $E \propto k^{-2}$ power law. As Reynolds number increases curling of the current sheets, due to Kelvin-Helmhotlz type instabilities and reconnection, modifies the scaling of the energy spectrum from $k^{-2}$ towards $k^{-5/3}$. This transition occurs at the highest Reynolds numbers of direct numerical simulations with resolutions up to $2048^3$ grid points. Finite Reynolds number behaviour is observed due to the initial correlations without reaching a finite asymptote for the energy dissipation rate and with an unexpected scaling between the Taylor and the integral scale Reynolds numbers, i.e. $Re_\lambda \propto Re^{2/3}$. Our results, therefore, demonstrate that even state-of-the-art numerical simulations of the highest resolution can be influenced by the choice of initial conditions and consequently they are inadequate to deduce unequivocally the fate of universality in MHD turbulence.

We investigate the origins of $k^{-2}$ spectrum in a decaying Taylor-Green magnetohydrodynamic flow with zero large scale magnetic flux that was reported in Lee et al. (2010). A possible candidate for this scaling exponent has been the weak turbulence phenomenology. From our numerical simulations, we observe that current sheets in the magnetic Taylor-Green flow are formed in regions of magnetic discontinuities. Based on this observation and by studying the influence of the current sheets on the energy spectrum, using a filtering technique, we argue that the discontinuities are responsible for the $-2$ power law scaling of the energy spectra of this flow.

An investigation of the dynamo instability close to the threshold produced by an ABC forced flow is presented. We focus on the on-off intermittency behavior of the dynamo and the counter-effect of the Lorentz force in the non-linear stage of the dynamo. The Lorentz force drastically alters the statistics of the turbulent fluctuations of the flow and reduces their amplitude. As a result much longer burst (on-phases) are observed than what is expected based on the amplitude of the fluctuations in the kinematic regime of the dynamo. For large Reynolds numbers, the duration time of the ``On'' phase follows a power law distribution, while for smaller Reynolds numbers the Lorentz force completely kills the noise and the system transits from a chaotic state into a ``laminar'' time periodic flow. The behavior of the On-Off intermittency as the Reynolds number is increased is also examined. The connections with dynamo experiments and theoretical modeling are discussed.

Many classical nova ejecta are enriched in CNO and Ne. Rosner et al. recently suggested that the enrichment might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and gravity waves at the interface between the envelope and the underlying C/O white dwarf. The shear flow amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/O into the superincumbent H/He. Using two-dimensional simulations, we formulate a quantitative expression for the amount of C/O per unit area that can be entrained into the H/He at saturation. The fraction of the envelope that is enriched depends on the horizontal distribution of shear velocity and the density contrast between the C/O white dwarf and the H/He layer but is roughly independent of the vertical shape of the shear profile. Using this parameterization for the mixed mass, we then perform several one-dimensional Lagrangian calculations of an accreting white dwarf envelope and consider two scenarios: that the wave breaking and mixing is driven by the convective flows; and that the mixing occurs prior to the onset of convection. In the absence of enrichment prior to ignition, the base of the convective zone, as calculated from mixing-length theory with the Ledoux instability criterion, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more violent runaway. The envelope can be enriched by <25% of C/O by mass (consistent with that observed in some ejecta) for shear velocities, over the surface, with Mach numbers <0.4.

Observed novae abundances and explosion energies estimated from observations indicate that there must be significant mixing of the heavier material of the white dwarf (C+O) into the lighter accreted material (H+He). Accordingly, nova models must incorporate a mechanism that will dredge up the heavier white dwarf material, and fluid motions from an early convection phase is one proposed mechanism. We present results from two-dimensional simulations of classical nova precursor models that demonstrate the beginning of a convective phase during the `simmering' of a Nova precursor. We use a new hydrostatic equilibrium hydrodynamics module recently developed for the adaptive-mesh code FLASH. The two-dimensional models are based on the one-dimensional models of Ami Glasner, and were evolved with FLASH from a pre-convective state to the onset of convection.

We present the results of a simulation of a wind-driven non-linear gravity wave breaking on the surface of a white dwarf. The ``wind'' consists of H/He from an accreted envelope, and the simulation demonstrates that this breaking wave mechanism can produce a well-mixed layer of H/He with C/O from the white dwarf above the surface. Material from this mixed layer may then be transported throughout the accreted envelope by convection, which would enrich the C/O abundance of the envelope as is expected from observations of novae.

We examine the linear stability of fluid interfaces subjected to a shear flow. Our main object is to generalize previous work to arbitrary Atwood number, and to allow for surface tension and weak compressibility. The motivation derives from instances in astrophysical systems where mixing across material interfaces driven by shear flows may significantly affect the dynamical evolution of these systems.

Using the results of recent work in shear instabilities in stratified fluids, we show that the resonant interaction between large-scale flows in the accreted H/He envelope of white dwarf stars and interfacial gravity waves can mix with the star's envelope with the white dwarf's surface material, leading to the enhancement of the envelope's C/O abundance to levels required by extant models for nova outbursts.