Source author record

P. Charbonneau

P. Charbonneau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR cond-mat.stat-mech cond-mat.dis-nn

Catalog footprint

What is connected

7works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Characterisation of the turbulent electromotive force and its magnetically-mediated quenching in a global EULAG-MHD simulation of solar convection

We perform a mean-field analysis of the EULAG-MHD millenium simulation of global magnetohydrodynamical convection presented in Passos et al. 2014. The turbulent electromotive force operating in the simulation is assumed to be linearly related to the cyclic axisymmetric mean magnetic field and its first spatial derivatives. At every grid point in the simulation's meridional plane, this assumed relationship involves 27 independent tensorial coefficients. Expanding on Racine et al. 2011, we extract these coefficients from the simulation data through a least-squares minimization procedure based on singular value decomposition. The reconstructed alpha-tensor shows good agreement with that obtained by Racine et al. 2011, who did not include derivatives of the mean-field in their fit, as well as with the alpha-tensor extracted by Augustson et al. 2015 from a distinct ASH MHD simulation. The isotropic part of the turbulent magnetic diffusivity tensor beta is positive definite and reaches values of 5.0x10^7 m2s-1 in the middle of the convecting fluid layers. The spatial variations of both alpha_phiphi and beta_phiphi component are well reproduced by expressions obtained under the SOCA, with a good matching of amplitude requiring a turbulent correlation time about five times smaller than the estimated turnover time of the small-scale turbulent flow. We find the magnetic quenching of the alpha-effect to be driven primarily by a reduction of the small-scale flow's kinetic helicity, with variations of the current helicity playing a lesser role in most locations in the simulation domain. Our measurements of turbulent diffusivity quenching are restricted to the beta_phiphi component, but indicate a weaker quenching, by a factor of 1.36, than of the alpha effect, which in our simulation drops by a factor of three between the minimum and maximum phases of the magnetic cycle.

preprint2016arXiv

Modelling turbulent stellar convection zones: sub-grid scales effects

The impressive development of global numerical simulations of turbulent stellar interiors unveiled a variety of possible differential rotation (solar or anti-solar), meridional circulation (single or multi-cellular), and dynamo states (stable large scale toroidal field or periodically reversing magnetic fields). Various numerical schemes, based on the so-called anelastic set of equations, were used to obtain these results. It appears today mandatory to assess their robustness with respect to the details of the numerics, and in particular to the treatment of turbulent sub-grid scales. We report on an ongoing comparison between two global models, the ASH and EULAG codes. In EULAG the sub-grid scales are treated implicitly by the numerical scheme, while in ASH their effect is generally modelled by using enhanced dissipation coefficients. We characterize the sub-grid scales effect in a turbulent convection simulation with EULAG. We assess their effect at each resolved scale with a detailed energy budget. We derive equivalent eddy-diffusion coefficients and use the derived diffusivities in twin ASH numerical simulations. We find a good agreement between the large-scale flows developing in the two codes in the hydrodynamic regime, which encourages further investigation in the magnetohydrodynamic regime for various dynamo solutions.

preprint2015arXiv

Evidence of active MHD instability in EULAG-MHD simulations of solar convection

We investigate the possible development of magnetohydrodynamical instabilities in the EULAG-MHD "millenium simulation" of Passos & Charbonneau (2014). This simulation sustains a large-scale magnetic cycle characterized by solar-like polarity reversals taking place on a regular multidecadal cadence, and in which zonally-oriented bands of strong magnetic field accumulate below the convective layers, in response to turbulent pumping from above in successive magnetic half-cycles. Key aspects of this simulation include low numerical dissipation and a strongly subadiabatic fluid layer underlying the convectively unstable layers corresponding to the modeled solar convection zone. These properties are conducive to the growth and development of two-dimensional instabilities otherwise suppressed by stronger dissipation. We find evidence for the action of a non-axisymmetric magnetoshear instability operating in the upper portions of the stably stratified fluid layers. We also investigate the possibility that the Tayler instability may be contributing to the destabilization of the large-scale axisymmetric magnetic component at high latitudes. On the basis of our analyses, we propose a global dynamo scenario whereby the magnetic cycle is driven primarily by turbulent dynamo action in the convecting layers, but MHD instabilities accelerate the dissipation of the magnetic field pumped down into the overshoot and stable layers, thus perhaps significantly influencing the magnetic cycle period. Support for this scenario is found in the distinct global dynamo behaviors observed in an otherwise identical EULAG-MHD simulations, using a different degree of subadiabaticity in the stable fluid layers underlying the convection zone.

preprint2015arXiv

Jamming Criticality Revealed by Removing Localized Buckling Excitations

Recent theoretical advances offer an exact, first-principle theory of jamming criticality in infinite dimension as well as universal scaling relations between critical exponents in all dimensions. For packings of frictionless spheres near the jamming transition, these advances predict that nontrivial power-law exponents characterize the critical distribution of (i) small inter-particle gaps and (ii) weak contact forces, both of which are crucial for mechanical stability. The scaling of the inter-particle gaps is known to be constant in all spatial dimensions $d$ -- including the physically relevant $d=2$ and 3, but the value of the weak force exponent remains the object of debate and confusion. Here, we resolve this ambiguity by numerical simulations. We construct isostatic jammed packings with extremely high accuracy, and introduce a simple criterion to separate the contribution of particles that give rise to localized buckling excitations, i.e., bucklers, from the others. This analysis reveals the remarkable dimensional robustness of mean-field marginality and its associated criticality.

preprint2012arXiv

Torsional Oscillations in a Global Solar Dynamo

We characterize and analyze rotational torsional oscillations developing in a large-eddy magnetohydrodynamical simulation of solar convection (Ghizaru, Charbonneau, and Smolarkiewicz, Astrophys. J. Lett., 715, L133 (2010); Racine et al., Astrophys. J., 735, 46 (2011)) producing an axisymmetric large-scale magnetic field undergoing periodic polarity reversals. Motivated by the many solar-like features exhibited by these oscillations, we carry out an analysis of the large-scale zonal dynamics. We demonstrate that simulated torsional oscillations are not driven primarily by the periodically-varying large-scale magnetic torque, as one might have expected, but rather via the magnetic modulation of angular-momentum transport by the large-scale meridional flow. This result is confirmed by a straightforward energy analysis. We also detect a fairly sharp transition in rotational dynamics taking place as one moves from the base of the convecting layers to the base of the thin tachocline-like shear layer formed in the stably stratified fluid layers immediately below. We conclude by discussing the implications of our analyses with regards to the mechanism of amplitude saturation in the global dynamo operating in the simulation, and speculate on the possible precursor value of torsional oscillations for the forecast of solar cycle characteristics.

preprint2010arXiv

Numerical and Theoretical Study of a Monodisperse Hard-Sphere Glass Former

There exists a variety of theories of the glass transition and many more numerical models. But because the models need built-in complexity to prevent crystallization, comparisons with theory can be difficult. We study the dynamics of a deeply supersaturated \emph{monodisperse} four-dimensional (4D) hard-sphere fluid, which has no such complexity, but whose strong intrinsic geometrical frustration inhibits crystallization, even when deeply supersaturated. As an application, we compare its behavior to the mode-coupling theory (MCT) of glass formation. We find MCT to describe this system better than any other structural glass formers in lower dimensions. The reduction in dynamical heterogeneity in 4D suggested by a milder violation of the Stokes-Einstein relation could explain the agreement. These results are consistent with a mean-field scenario of the glass transition.

preprint2007arXiv

Dynamical Heterogeneity and Nonlinear Susceptibility in Short-Ranged Attractive Supercooled Liquids

Recent work has demonstrated the strong qualitative differences between the dynamics near a glass transition driven by short-ranged repulsion and one governed by short-ranged attraction. Here, we study in detail the behavior of non-linear, higher-order correlation functions that measure the growth of length scales associated with dynamical heterogeneity in both types of systems. We find that this measure is qualitatively different in the repulsive and attractive cases with regards to the wave vector dependence as well as the time dependence of the standard non-linear four-point dynamical susceptibility. We discuss the implications of these results for the general understanding of dynamical heterogeneity in glass-forming liquids.