184 research outputs found

    Coupled Kinetic-Fluid Simulations of Ganymede's Magnetosphere and Hybrid Parallelization of the Magnetohydrodynamics Model

    Full text link
    The largest moon in the solar system, Ganymede, is the only moon known to possess a strong intrinsic magnetic field. The interaction between the Jovian plasma and Ganymede's magnetic field creates a mini-magnetosphere with periodically varying upstream conditions, which creates a perfect laboratory in nature for studying magnetic reconnection and magnetospheric physics. Using the latest version of Space Weather Modeling Framework (SWMF), we study the upstream plasma interactions and dynamics in this subsonic, sub-Alfvénic system. We have developed a coupled fluid-kinetic Hall Magnetohydrodynamics with embedded Particle-in-Cell (MHD-EPIC) model for Ganymede's magnetosphere, with a self-consistently coupled resistive body representing the electrical properties of the moon's interior, improved inner boundary conditions, and high resolution charge and energy conserved PIC scheme. I reimplemented the boundary condition setup in SWMF for more versatile control and functionalities, and developed a new user module for Ganymede's simulation. Results from the models are validated with Galileo magnetometer data of all close encounters and compared with Plasma Subsystem (PLS) data. The energy fluxes associated with the upstream reconnection in the model is estimated to be about 10^-7 W/cm^2, which accounts for about 40% to the total peak auroral emissions observed by the Hubble Space Telescope. We find that under steady upstream conditions, magnetopause reconnection in our fluid-kinetic simulations occurs in a non-steady manner. Flux ropes with length of Ganymede's radius form on the magnetopause at a rate about 3/minute and create spatiotemporal variations in plasma and field properties. Upon reaching proper grid resolutions, the MHD-EPIC model can resolve both electron and ion kinetics at the magnetopause and show localized crescent shape distribution in both ion and electron phase space, non-gyrotropic and non-isotropic behavior inside the diffusion regions. The estimated global reconnection rate from the models is about 80 kV with 60% efficiency. There is weak evidence of sim1sim 1 minute periodicity in the temporal variations of the reconnection rate due to the dynamic reconnection process. The requirement of high fidelity results promotes the development of hybrid parallelized numerical model strategy and faster data processing techniques. The state-of-the-art finite volume/difference MHD code Block Adaptive Tree Solarwind Roe Upwind Scheme (BATS-R-US) was originally designed with pure MPI parallelization. The maximum problem size achievable was limited by the storage requirements of the block tree structure. To mitigate this limitation, we have added multithreaded OpenMP parallelization to the previous pure MPI implementation. We opt to use a coarse-grained approach by making the loops over grid blocks multithreaded and have succeeded in making BATS-R-US an efficient hybrid parallel code with modest changes in the source code while preserving the performance. Good weak scalings up to 50,0000 and 25,0000 of cores are achieved for the explicit and implicit time stepping schemes, respectively. This parallelization strategy greatly extends the possible simulation scale by an order of magnitude, and paves the way for future GPU-portable code development. To improve visualization and data processing, I have developed a whole new data processing workflow with the Julia programming language for efficient data analysis and visualization. As a summary, 1. I build a single fluid Hall MHD-EPIC model of Ganymede's magnetosphere; 2. I did detailed analysis of the upstream reconnection; 3. I developed a MPI+OpenMP parallel MHD model with BATSRUS; 4. I wrote a package for data analysis and visualization.PHDClimate and Space Sciences and EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/163032/1/hyzhou_1.pd

    Turbulent magnetic field amplification from spiral SASI modes in core-collapse supernovae

    Full text link
    We describe the initial implementation of magnetohydrodynamics (MHD) in our astrophysical simulation code \genasis. Then, we present MHD simulations exploring the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields by adding a weak magnetic field to an initially spherically symmetric fluid configuration that models a stalled shock in the post-bounce supernova environment. Upon perturbation and nonlinear SASI development, shear flows associated with the spiral SASI mode contributes to a widespread and turbulent field amplification mechanism. While the SASI may contribute to neutron star magnetization, these simulations do not show qualitatively new features in the global evolution of the shock as a result of SASI-induced magnetic field amplification.Comment: 15 pages, 7 figures, To appear in the Journal of Physics: Conference Series. Proceedings of the IUPAP Conference on Computational Physics (CCP2011

    On the Exchange of Kinetic and Magnetic Energy Between Clouds and the Interstellar Medium

    Get PDF
    We investigate, through 2D MHD numerical simulations, the interaction of a uniform magnetic field oblique to a moving interstellar cloud. In particular we explore the transformation of cloud kinetic energy into magnetic energy as a result of field line stretching. Some previous simulations have emphasized the possible dynamical importance of a ``magnetic shield'' formed around clouds when the magnetic field is perpendicular to the cloud motion (Jones et al. 1996, Miniati et al. 1998). It was not clear, however, how dependent those findings were to the assumed field configuration and cloud properties. To expand our understanding of this effect, we examine several new cases by varing the magnetic field orientation angle with respect to the cloud motion (\theta), the cloud-background density contrast, and the cloud Mach number. We show that in 2D and with \theta large enough, the magnetic field tension can become dominant in the dynamics of the motion of high density contrast, low Mach number clouds. In such cases a significant fraction of cloud kinetic energy can be transformed into magnetic energy with the magnetic pressure at the cloud nose exceeding the ram pressure of the impinging flow. We derive a characteristic timescale for this process of energy ``conversion''. We find also that unless the cloud motion is highly aligned to the magnetic field, reconnection through tearing mode instabilities in the cloud wake limit the formation of a strong flux rope feature following the cloud. Finally we attempt to interpret some observational properties of the magnetic field in view of our results.Comment: 24 pages in aaspp4 Latex and 7 figures. Accepted for publication in The Astrophysical Journa

    What Controls the Structure and Dynamics of Earth's Magnetosphere?

    Get PDF

    Modeling a Coronal Mass Ejection from an Extended Filament Channel. I. Eruption and Early Evolution

    Get PDF
    We present observations and modeling of the magnetic field configuration, morphology, and dynamics of a large-scale, high-latitude filament eruption observed by the Solar Dynamics Observatory. We analyze the 2015 July 9-10 filament eruption and the evolution of the resulting coronal mass ejection (CME) through the solar corona. The slow streamer-blowout CME leaves behind an elongated post-eruption arcade above the extended polarity inversion line that is only poorly visible in extreme ultraviolet (EUV) disk observations and does not resemble a typical bright flare-loop system. Magnetohydrodynamic (MHD) simulation results from our data-inspired modeling of this eruption compare favorably with the EUV and white-light coronagraph observations. We estimate the reconnection flux from the simulation's flare-arcade growth and examine the magnetic-field orientation and evolution of the erupting prominence, highlighting the transition from an erupting sheared-arcade filament channel into a streamer-blowout flux-rope CME. Our results represent the first numerical modeling of a global-scale filament eruption where multiple ambiguous and complex observational signatures in EUV and white light can be fully understood and explained with the MHD simulation. In this context, our findings also suggest that the so-called stealth CME classification, as a driver of unexpected or "problem" geomagnetic storms, belongs more to a continuum of observable/nonobservable signatures than to separate or distinct eruption processes.Peer reviewe
    corecore