EXTENSION OF THE MURAM RADIATIVE MHD CODE FOR CORONAL SIMULATIONS

We present a new version of the MURaM radiative magnetohydrodynamics (MHD) code that allows for simulations spanning from the upper convection zone into the solar corona. We implement the relevant coronal physics in terms of optically thin radiative loss, field aligned heat conduction, and an equilibrium ionization equation of state. We artificially limit the coronal Alfven and heat conduction speeds to computationally manageable values using an approximation to semi-relativistic MHD with an artificially reduced speed of light (Boris correction). We present example solutions ranging from quiet to active Sun in order to verify the validity of our approach. We quantify the role of numerical diffusivity for the effective coronal heating. We find that the (numerical) magnetic Prandtl number determines the ratio of resistive to viscous heating and that owing to the very large magnetic Prandtl number of the solar corona, heating is expected to happen predominantly through viscous dissipation. We find that reasonable solutions can be obtained with values of the reduced speed of light just marginally larger than the maximum sound speed. Overall this leads to a fully explicit code that can compute the time evolution of the solar corona in response to photospheric driving using numerical time steps not much smaller than 0.1 s. Numerical simulations of the coronal response to flux emergence covering a time span of a few days are well within reach using this approach.

[1]  R. Canfield,et al.  Flare loop radiative hydrodynamics. V - Response to thick-target heating. VI - Chromospheric evaporation due to heating by nonthermal electrons. VII - Dynamics of the thick-target heated chromosphere , 1985 .

[2]  M. Linton,et al.  GAS-DYNAMIC SHOCK HEATING OF POST-FLARE LOOPS DUE TO RETRACTION FOLLOWING LOCALIZED, IMPULSIVE RECONNECTION , 2008, 0810.3661.

[3]  M. Rempel.  NUMERICAL SIMULATIONS OF QUIET SUN MAGNETISM: ON THE CONTRIBUTION FROM A SMALL-SCALE DYNAMO , 2014, 1405.6814.

[4]  John Lyon,et al.  The Lyon-Fedder-Mobarry (LFM) global MHD magnetospheric simulation code , 2004 .

[5]  H. Peter,et al.  Using coronal seismology to estimate the magnetic field strength in a realistic coronal model , 2015, 1508.00593.

[6]  B. Gudiksen,et al.  An Ab Initio Approach to the Solar Coronal Heating Problem , 2004, astro-ph/0407266.

[7]  T. Gombosi,et al.  ALFVÉN WAVE SOLAR MODEL (AWSoM): CORONAL HEATING , 2013, 1311.4093.

[8]  M. Cheung,et al.  A model for the formation of the active region corona driven by magnetic flux emergence , 2014, 1402.5343.

[9]  M. Knoelker,et al.  RADIATIVE MAGNETOHYDRODYNAMIC SIMULATION OF SUNSPOT STRUCTURE , 2008, 0808.3294.

[10]  C. Schrijver,et al.  THERMAL DIAGNOSTICS WITH THE ATMOSPHERIC IMAGING ASSEMBLY ON BOARD THE SOLAR DYNAMICS OBSERVATORY: A VALIDATED METHOD FOR DIFFERENTIAL EMISSION MEASURE INVERSIONS , 2015, 1504.03258.

[11]  A. Brandenburg.  MAGNETIC PRANDTL NUMBER DEPENDENCE OF THE KINETIC-TO-MAGNETIC DISSIPATION RATIO , 2014, 1404.6964.

[12]  T. Emonet,et al.  Simulations of magneto-convection in the solar photosphere Equations, methods, and results of the MURaM code , 2005 .

[13]  Coronal Heating Through Braiding of Magnetic Field Lines Synthesized Coronal EUV Emission and Magnetic Structure , 2004, astro-ph/0409504.

[14]  Non-Fickian diffusion and tau approximation from numerical turbulence , 2003, astro-ph/0306521.

[15]  Forward Modeling of the Corona of the Sun and Solar-like Stars: From a Three-dimensional Magnetohydrodynamic Model to Synthetic Extreme-Ultraviolet Spectra , 2006 .

[16]  S. Bingert,et al.  Observationally driven 3D magnetohydrodynamics model of the solar corona above an active region , 2013, 1305.5693.

[17]  J. Linker,et al.  The Formation of Coronal Loops by Thermal Instability in Three Dimensions , 2008 .

[18]  Sven Bingert,et al.  Intermittent heating in the solar corona employing a 3D MHD model , 2011, 1103.6042.

[19]  T. Aslam,et al.  A second-order accurate Super TimeStepping formulation for anisotropic thermal conduction , 2012 .

[20]  V. Hansteen,et al.  A 16‐moment solar wind model: From the chromosphere to 1 AU , 2001 .

[21]  H. Mason,et al.  CHIANTI—AN ATOMIC DATABASE FOR EMISSION LINES. XII. VERSION 7 OF THE DATABASE , 2012 .

[22]  B. Gudiksen,et al.  An Ab Initio Approach to Solar Coronal Loops , 2004, astro-ph/0407267.

[23]  A new approach to turbulent transport of a mean scalar , 2003, astro-ph/0302013.

[24]  W. Abbett,et al.  The Magnetic Connection between the Convection Zone and Corona in the Quiet Sun , 2007 .

[25]  M. Rempel.  NUMERICAL SIMULATIONS OF SUNSPOT DECAY: ON THE PENUMBRA–EVERSHED FLOW–MOAT FLOW CONNECTION , 2015, 1511.01410.

[26]  D. D. Zeeuw,et al.  Semirelativistic Magnetohydrodynamics and Physics-Based Convergence Acceleration , 2002 .

[27]  P. Cargill,et al.  THE INFLUENCE OF NUMERICAL RESOLUTION ON CORONAL DENSITY IN HYDRODYNAMIC MODELS OF IMPULSIVE HEATING , 2013, 1305.1902.

[28]  D. Longcope,et al.  SHOCKS AND THERMAL CONDUCTION FRONTS IN RETRACTING RECONNECTED FLUX TUBES , 2010, 1006.4398.

[29]  Tamas I. Gombosi,et al.  The telegraph equation in charged particle transport , 1993 .

[30]  H. Peter,et al.  Forward modeling of stellar coronae: from a 3D MHD model to synthetic EUV spectra , 2005, astro-ph/0503342.

[31]  E. Parker.  Nanoflares and the solar X-ray corona , 1988 .

[32]  Å. Nordlund,et al.  Bulk Heating and Slender Magnetic Loops in the Solar Corona , 2002 .

[33]  J. Linker,et al.  Calculating the Thermal Structure of Solar Active Regions in Three Dimensions , 2005 .

[34]  L. Spitzer.  Physics of fully ionized gases , 1956 .

[35]  H. Peter,et al.  Nanoflare statistics in an active region 3D MHD coronal model , 2012, 1211.6417.

[36]  M. Berger,et al.  ON THE RELATIONSHIP BETWEEN PHOTOSPHERIC FOOTPOINT MOTIONS AND CORONAL HEATING IN SOLAR ACTIVE REGIONS , 2014 .

[37]  E. Parker.  Topological dissipation and the small-scale fields in turbulent gases. , 1972 .

[38]  C. Munz,et al.  Hyperbolic divergence cleaning for the MHD equations , 2002 .

[39]  E. N. Parker,et al.  Magnetic Neutral Sheets in Evolving Fields - Part Two - Formation of the Solar Corona , 1983 .

[40]  James A. Klimchuk,et al.  On Solving the Coronal Heating Problem , 2006 .

[41]  R. Schunk.  Transport equations for aeronomy , 1975 .

[42]  Forrest J. Rogers,et al.  Opal equation-of-state tables for astrophysical applications , 1996 .

[43]  A. Brandenburg.  NONLINEAR SMALL-SCALE DYNAMOS AT LOW MAGNETIC PRANDTL NUMBERS , 2011, 1106.5777.

[44]  M. Cheung,et al.  Magnetic jam in the corona of the Sun , 2015, Nature Physics.

[45]  George L. Withbroe,et al.  Mass and Energy Flow in the Solar Chromosphere and Corona , 1977 .