BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jorrit Leenaarts

Jorrit Leenaarts contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

19 published item(s)

preprint2022arXiv

Radiative losses in the chromosphere during a C-class flare

Solar flares release an enormous amount of energy into the corona. A substantial fraction of this energy is transported to the lower atmosphere, which results in chromospheric heating. The mechanisms that transport energy to the lower solar atmosphere during a flare are still not fully understood. We aim to estimate the temporal evolution of the radiative losses in the chromosphere at the footpoints of a C-class flare, in order to set observational constraints on the electron beam parameters of a RADYN flare simulation. We estimated the radiative losses from hydrogen, and singly ionized Ca and Mg using semi-empirical model atmospheres. To estimate the integrated radiative losses in the chromosphere the net cooling rates were integrated between the temperature minimum and the height where the temperature reaches 10 kK. The stratification of the net cooling rate suggests that the Ca IR triplet lines are responsible for most of the radiative losses in the flaring atmosphere. During the flare peak time, the contribution from Ca II H & K and Mg II h & k lines are strong and comparable to the Ca IR triplet ($\sim$32 kW m$^{-2}$). Since our flare is a relatively weak event the chromosphere is not heated above 11 kK, which in turn yields a subdued Lyα contribution ($\sim$7 kW m$^{-2}$). The temporal evolution of total integrated radiative losses exhibits sharply-rising losses (0.4 kW m$^{-2}$ s$^{-1}$) and a relatively slow decay (0.23 kW~m$^{-2}$ s$^{-1}$). The maximum value of total radiative losses is reached around the flare peak time, and can go up to 175 kW m$^{-2}$ for a single pixel located at footpoint. After a small parameter study, we find the best model-data consistency in terms of the amplitude of radiative losses and the overall atmospheric structure with a RADYN flare simulation in the injected energy flux of $5\times10^{10}$ erg s$^{-1}$ cm$^{-2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mapping Solar Magnetic Fields from the Photosphere to the Base of the Corona

Routine ultraviolet imaging of the Sun's upper atmosphere shows the spectacular manifestation of solar activity; yet we remain blind to its main driver, the magnetic field. Here we report unprecedented spectropolarimetric observations of an active region plage and its surrounding enhanced network, showing circular polarization in ultraviolet (Mg II $h$ & $k$ and Mn I) and visible (Fe I) lines. We infer the longitudinal magnetic field from the photosphere to the very upper chromosphere. At the top of the plage chromosphere the field strengths reach more than 300 gauss, strongly correlated with the Mg II $k$ line core intensity and the electron pressure. This unique mapping shows how the magnetic field couples the different atmospheric layers and reveals the magnetic origin of the heating in the plage chromosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physical properties of bright Ca II K fibrils in the solar chromosphere

Broad-band images of the solar chromosphere in the Ca II H&K line cores around active regions are covered with fine bright elongated structures called bright fibrils. The mechanisms that form these structures and cause them to appear bright are still unknown. We aim to investigate the physical properties, such as temperature, line-of-sight velocity, and microturbulence, in the atmosphere that produces bright fibrils and to compare those to the properties of their surrounding atmosphere. We used simultaneous observations of a plage region in Fe I 6301-2 Å, Ca II 8542 Å, Ca II K, and H$α$ acquired by the CRISP and CHROMIS instruments on the Swedish 1-m Solar Telescope. We manually selected a sample of 282 Ca II K bright fibrils. We compared the appearance of the fibrils in our sample to the Ca II 8542 Å and H$α$ data. We performed non-local thermodynamic equilibrium (non-LTE) inversions using the inversion code STiC on the Fe I 6301-2 Å, Ca II 8542 Å, Ca II K lines to infer the physical properties of the atmosphere. The line profiles in bright fibrils have a higher intensity in their K$_2$ peaks compared to profiles formed in the surrounding atmosphere. The inversion results show that the atmosphere in fibrils is on average $100-200$~K hotter at an optical depth log$(τ) = -4.3$ compared to their surroundings. The line-of-sight velocity at chromospheric heights in the fibrils does not show any preference towards upflows or downflows. The microturbulence in the fibrils is on average 0.5 km s$^{-1}$ higher compared to their surroundings. Our results suggest that the fibrils have a limited extent in height, and they should be viewed as hot threads pervading the chromosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Radiation hydrodynamics in simulations of the solar atmosphere

Nearly all energy generated by fusion in the solar core is ultimately radiated away into space in the solar atmosphere, while the remaining energy is carried away in the form of neutrinos. The exchange of energy between the solar gas and the radiation field is thus an essential ingredient of atmospheric modeling. The equations describing these interactions are known, but their solution is so computationally expensive that they can only be solved in approximate form in multi-dimensional radiation-MHD modeling. In this review, I discuss the most commonly used approximations for energy exchange between gas and radiation in the photosphere, chromosphere, and corona.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Ion-neutral interactions and non-equilibrium ionization in the solar chromosphere

The thermal structure of the chromosphere is regulated through a complex interaction of various heating processes, radiative cooling, and the ionization degree of the plasma. Here we study the impact on the thermal properties of the chromosphere when including the combined action of non-equilibrium ionization (NEI) of hydrogen and helium and ion-neutral interaction effects. We have performed a 2.5D radiative magnetohydrodynamic simulation including ion-neutral interaction effects by solving the generalized Ohm's law (GOL) as well as NEI for hydrogen and helium using the Bifrost code. The GOL equation includes ambipolar diffusion and the Hall term. We compare this simulation with another simulation that computes the ionization in local thermodynamic equilibrium (LTE) including ion-neutral interaction effects. Our numerical models reveal substantial thermal differences in magneto-acoustic shocks, the wake behind the shocks, spicules, low-lying magnetic loops, and the transition region. In particular, we find that heating through ambipolar diffusion in shock wakes is substantially less efficient, while in the shock fronts themselves it is more efficient, under NEI conditions than when assuming LTE.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On fibrils and field lines: The nature of H$α$ fibrils in the solar chromosphere

Observations of the solar chromosphere in the line-core of the \Halpha\ line show dark elongated structures called fibrils that show swaying motion. We performed a 3D radiation-MHD simulation of a network region, and computed synthetic \Halpha\ images from this simulation to investigate the relation between fibrils and the magnetic field lines in the chromosphere. The periods, amplitudes and phase-speeds of the synthetic fibrils are consistent with those observed. We analyse the relation between the synthetic fibrils and the field lines threading through them, and find that some fibrils trace out the same field line along the fibril's length, but there are also fibrils that sample different field lines at different locations along their length. Fibrils sample the same field lines on a time scale of $\sim200$~s. This is shorter than their own lifetime. We analysed the evolution of the atmosphere along a number of field lines that thread through fibrils and find that they carry slow-mode waves that load mass into the field line, as well as transverse waves that propagate with the Alfvén speed. Transverse waves propagating in opposite directions cause an interference pattern with complex apparent phase speeds. The relationship between fibrils and field lines is complex. It is governed by constant migration and swaying of the field lines, their mass loading by slow modes and subsequent draining, and their actual visibility in \Halpha. Field lines are visible where they lie close to the optical depth unity surface. The location of the latter is governed by the height at which the column mass in the chromosphere reaches a certain value. We conclude that using the swaying motion of fibrils as a tracer of chromospheric transverse oscillations must be done with caution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Three-dimensional radiative transfer simulations of the scattering polarization of the hydrogen Ly$α$ line in a MHD model of the chromosphere-corona transition region

Probing the magnetism of the upper solar chromosphere requires measuring and modeling the scattering polarization produced by anisotropic radiation pumping in UV spectral lines. Here we apply PORTA (a novel radiative transfer code) to investigate the hydrogen Ly$α$ line in a 3D model of the solar atmosphere resulting from a state of the art MHD simulation. At full spatial resolution the linear polarization signals are very significant all over the solar disk, with a large fraction of the field of view showing line-center amplitudes well above the 1% level. Via the Hanle effect the line-center polarization signals are sensitive to the magnetic field of the model's transition region, even when its mean field strength is only 15 G. The breaking of the axial symmetry of the radiation field produces significant forward-scattering polarization in Ly$α$, without the need of an inclined magnetic field. Interestingly, the Hanle effect tends to decrease such forward-scattering polarization signals in most of the points of the field of view. When the spatial resolution is degraded, the line-center polarization of Ly$α$ drops below the 1% level, reaching values similar to those previously found in 1D semi-empirical models (i.e., up to about 0.5%). The center to limb variation of the spatially-averaged polarization signals is qualitatively similar to that found in 1D models, with the largest line-center amplitudes at $μ=\cosθ\approx 0.4$ ($θ$ being the heliocentric angle). These results are important, both for designing the needed space-based instrumentation and for a reliable interpretation of future observations of the Ly$α$ polarization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detailed and simplified non-equilibrium helium ionization in the solar atmosphere

Helium ionization plays an important role in the energy balance of the upper chromosphere and transition region. Helium spectral lines are also often used as diagnostics of these regions. We carry out 1D radiation-hydrodynamics simulations of the solar atmosphere and find that the helium ionization is mostly set by photoionization and direct collisional ionization, counteracted by radiative recombination cascades. By introducing an additional recombination rate mimicking the recombination cascades, we construct a simplified 3 level helium model atom consisting of only the ground states. This model atom is suitable for modeling non-equilibrium helium ionization in 3D numerical models. We perform a brief investigation of the formation of the He I 10830 and He II 304 spectral lines. Both lines show non-equilibrium features that are not recovered with statistical equilibrium models, and caution should therefore be exercised when such models are used as a basis in the interpretation of observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The effect of isotopic splitting on the bisector and inversions of the solar Ca II 854.2 nm line

The Ca II 854.2 nm spectral line is a common diagnostic of the solar chromosphere. The average line profile shows an asymmetric core, and its bisector shows a characteristic inverse-C shape. The line actually consists of six components with slightly different wavelengths depending on the isotope of calcium. This isotopic splitting of the line has been taken into account in studies of non-solar stars, but never for the Sun. We performed non-LTE radiative transfer computations from three models of the solar atmosphere and show that the asymmetric line-core and inverse C-shape of the bisector of the 854.2 nm line can be explained by isotopic splitting. We confirm this finding by analysing observations and showing that the line asymmetry is present irrespective of conditions in the solar atmosphere. Finally, we show that inversions based on the Ca II 854.2 nm line should take the isotopic splitting into account, otherwise the inferred atmospheres will contain erroneous velocity gradients and temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Detailed Comparison Between The Observed and Synthesized Properties of a Simulated Type ii Spicule

We performed a 3D radiative MHD simulation of the solar atmosphere. This simulation shows a jet-like feature that shows similarities to the type II spicules observed for the first time with Hinode. Rapid Blueshifted Events (RBEs) on the solar disk are associated with these spicules. Observational results suggest they may contribute significantly in supplying the corona with hot plasma. We perform a detailed comparison of the properties of the simulated jet with those of type II spicules (observed with Hinode) and RBEs (with ground-based instruments). We analyze variety of synthetic emission and absorption lines from the simulations including chromospheric Ca II and Ha to TR and coronal temperatures (10E4 to several 10E6K). We compare their synthetic intensities, line profiles, Doppler shifts, line widths and asymmetries with observations from Hinode/SOT and EIS, SOHO/SUMER, SST and SDO/AIA. Many properties of the synthetic observables resemble the observations, and we describe in detail the physical processes that lead to these observables. Detailed analysis of the synthetic observables provides insight into how observations should be analyzed to derive information about physical variables in such a dynamic event. For example, we find that LOS superposition in the optically thin atmosphere requires the combination of Doppler shifts and spectral line asymmetry to determine the velocity in the jet. Other properties differ from the observations, especially in the chromospheric lines. The mass density of the part of the spicule with a chromospheric temperature is too low to produce significant opacity in chromospheric lines. These and other discrepancies are described in detail, and we discuss which mechanisms and physical processes may need to be included in the MHD simulations to mimic the thermodynamic processes of the chromosphere and corona, in particular to reproduce type II spicules.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How realistic are solar model atmospheres?

Recently, new solar model atmospheres have been developed to replace classical 1D LTE hydrostatic models and used to for example derive the solar chemical composition. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D MHD model (with an imposed 10 mT vertical magnetic field), 1D models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & Müller. We confront the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLV), absolute continuum fluxes, and the wings of hydrogen lines. We also test the 3D models for the intensity distribution of the granulation and spectral line shapes. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & Müller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models are worse, given their steeper temperature gradients. For the continuum fluxes, predictions for most models agree well with the observations. No model fits all hydrogen lines perfectly, but again the 3D model comes ahead. The 3D model also reproduces the observed continuum intensity fluctuations and spectral line shapes very well. The excellent agreement of the 3D model with the observables reinforces the view that its temperature structure is realistic. It outperforms the MHD simulation in all diagnostics, implying that recent claims for revised abundances based on MHD modelling are premature. Several weaknesses in the 1D models are exposed. The differences between the PHOENIX LTE and NLTE models are small. We conclude that the 3D hydrodynamical model is superior to any of the tested 1D models, which gives further confidence in the solar abundance analyses based on it.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The formation of IRIS diagnostics. III. Near-ultraviolet Spectra and Images

The Mg II h&k lines are the prime chromospheric diagnostics of NASA's Interface Region Imaging Spectrograph (IRIS). In the previous papers of this series we used a realistic three-dimensional radiative magnetohydrodynamics model to calculate the h&k lines in detail and investigated how their spectral features relate to the underlying atmosphere. In this work, we employ the same approach to investigate how the h&k diagnostics fare when taking into account the finite resolution of IRIS and different noise levels. In addition, we investigate the diagnostic potential of several other photospheric lines and near-continuum regions present in the near-ultraviolet (NUV) window of IRIS and study the formation of the NUV slit-jaw images. We find that the instrumental resolution of IRIS has a small effect on the quality of the h&k diagnostics; the relations between the spectral features and atmospheric properties are mostly unchanged. The peak separation is the most affected diagnostic, but mainly due to limitations of the simulation. The effects of noise start to be noticeable at a signal-to-noise ratio (S/N) of 20, but we show that with noise filtering one can obtain reliable diagnostics at least down to a S/N of 5. The many photospheric lines present in the NUV window provide velocity information for at least eight distinct photospheric heights. Using line-free regions in the h&k far wings we derive good estimates of photospheric temperature for at least three heights. Both of these diagnostics, in particular the latter, can be obtained even at S/Ns as low as 5.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Approximations for radiative cooling and heating in the solar chromosphere

Context. The radiative energy balance in the solar chromosphere is dominated by strong spectral lines that are formed out of LTE. It is computationally prohibitive to solve the full equations of radiative transfer and statistical equilibrium in 3D time dependent MHD simulations. Aims. To find simple recipes to compute the radiative energy balance in the dominant lines under solar chromospheric conditions. Methods. We use detailed calculations in time-dependent and 2D MHD snapshots to derive empirical formulae for the radiative cooling and heating. Results. The radiative cooling in neutral hydrogen lines and the Lyman continuum, the H and K and intrared triplet lines of singly ionized calcium and the h and k lines of singly ionized magnesium can be written as a product of an optically thin emission (dependent on temperature), an escape probability (dependent on column mass) and an ionization fraction (dependent on temperature). In the cool pockets of the chromosphere the same transitions contribute to the heating of the gas and similar formulae can be derived for these processes. We finally derive a simple recipe for the radiative heating of the chromosphere from incoming coronal radiation. We compare our recipes with the detailed results and comment on the accuracy and applicability of the recipes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fast approximation of angle-dependent partial redistribution in moving atmospheres

Radiative transfer modeling of spectral lines including partial redistribution (PRD) effects requires the evaluation of the ratio of the emission to the absorption profile. This quantity requires a large amount of computational work if one employs the angle-dependent redistribution function, which prohibits its use in 3D radiative transfer computations with model atmospheres containing velocity fields. We aim to provide a method to compute the emission to absorption profile ratio that requires less computational work but retains the effect of angle-dependent scattering in the resulting line profiles. We present a method to compute the profile ratio that employs the angle-averaged redistribution function and wavelength transforms to and from the rest frame of the scattering particles. We compare the emergent line profiles of the \MgII\,k and \Lyalpha\ lines computed with angle-dependent PRD, angle-averaged PRD and our new method in two representative test atmospheres. The new method yields a good approximation of true angle-dependent profile ratio and the resulting emergent line profiles while keeping the computational speed and simplicity of angle-averaged PRD theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-local thermodynamic equilibrium inversions from a 3D MHD chromospheric model

The structure of the solar chromosphere is believed to be governed by magnetic fields, even in quiet-Sun regions that have a relatively weak photospheric field. During the past decade inversion methods have emerged as powerful tools for analyzing the chromosphere of active regions. The applicability of inversions to infer the stratification of the physical conditions in a dynamic 3D solar chromosphere has not yet been studied in detail. This study aims to establish the diagnostic capabilities of non-local thermodynamical equilibrium (NLTE) inversion techniques of Stokes profiles induced by the Zeeman effect in the Ca II 8542 line. We computed the Ca II atomic level populations in a snapshot from a 3D radiation-MHD simulation of the quiet solar atmosphere in non-LTE using the 3D radiative transfer code Multi3d. These populations were used to compute synthetic full-Stokes profiles in the Ca II 8542 line using 1.5D radiative transfer and the inversion code Nicole. The profiles were then spectrally degraded to account for finite filter width and Gaussian noise was added to account for finite photon flux. These profiles were inverted using Nicole and the results were compared with the original model atmosphere. Our NLTE inversions applied to quiet-Sun synthetic observations provide reasonably good estimates of the chromospheric magnetic field, line-of-sight velocities and somewhat less accurate, but still very useful, estimates of the temperature. Three dimensional scattering of photons cause cool pockets in the chromosphere to be invisible in the line profile and consequently they are also not recovered by the inversions. To successfully detect Stokes linear polarization in this quiet snapshot, a noise level below 10^{-3.5} is necessary.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The formation of the Halpha line in the solar chromosphere

We use state-of-the-art radiation-MHD simulations and 3D non-LTE radiative transfer computations to investigate \Halpha\ line formation in the solar chromosphere and apply the results of this investigation to develop the potential of \Halpha\ as diagnostic of the chromosphere. We show that one can accurately model \Halpha\ line formation assuming statistical equilibrium and complete frequency redistribution provided the computation of the model atmosphere included non-equilibrium ionization of hydrogen, and the Lyman-$α$ and Lyman-$β$ line profiles are described by Doppler profiles. We find that 3D radiative transfer is essential in modeling hydrogen lines due to the low photon destruction probability in \Halpha. The \Halpha\ opacity in the upper chromosphere is mainly sensitive to the mass density and only weakly sensitive to temperature. We find that the \Halpha\ line-core intensity is correlated with the average formation height: the larger the average formation height, the lower the intensity. The line-core width is a measure of the gas temperature in the line-forming region. The fibril-like dark structures seen in \Halpha\ line-core images computed from our model atmosphere are tracing magnetic field lines. These structures are caused by field-aligned ridges of enhanced chromospheric mass density that raise their average formation height, and therefore makes them appear dark against their deeper-formed surroundings. We compare with observations, and find that the simulated line-core widths are very similar to the observed ones, without the need for additional microturbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Hanle effect of Ly$α$ in an MHD model of the Solar Transition Region

In order to understand the heating of the solar corona it is crucial to obtain empirical information on the magnetic field in its lower boundary (the transition region). To this end, we need to measure and model the linear polarization produced by scattering processes in strong UV lines, such as the hydrogen Ly$α$ line. The interpretation of the observed Stokes profiles will require taking into account that the outer solar atmosphere is highly structured and dynamic, and that the height of the transition region may well vary from one place in the atmosphere to another. Here we report on the Ly$α$ scattering polarization signals we have calculated in a realistic model of an enhanced network region, resulting from a state-of-the-art radiation MHD simulation. This model is characterized by spatially complex variations of the physical quantities at transition region heights. The results of our investigation lead us to emphasize that scattering processes in the upper solar chromosphere should indeed produce measurable linear polarization in Ly$α$. More importantly, we show that via the Hanle effect the model's magnetic field produces significant changes in the emergent $Q/I$ and $U/I$ profiles. Therefore, we argue that by measuring the polarization signals produced by scattering processes and the Hanle effect in Ly$α$ and contrasting them with those computed in increasingly realistic atmospheric models, we should be able to decipher the magnetic, thermal and dynamic structure of the upper chromosphere and transition region of the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the minimum temperature of the quiet solar chromosphere

aims: We aim to provide an estimate of the minimum temperature of the quiet solar chromosphere. methods: We perform a 2D radiation-MHD simulation spanning the upper convection zone to the lower corona. The simulation includes non-LTE radiative transfer and {an equation-of-state that includes non-equilibrium ionization of hydrogen and non-equilibrium H_2 molecule formation}. We analyze the reliability of the various assumptions made in our model in order to assess the realism of the simulation. results: Our simulation contains pockets of cool gas with down to 1660 K from 1 Mm up to 3.2 Mm height. It overestimates the radiative heating, and contains non-physical heating below 1660 K. Therefore we conclude that cool pockets in the quiet solar chromosphere might have even lower temperatures than in the simulation, provided that there exist areas in the chromosphere without significant magnetic heating. We suggest off-limb molecular spectroscopy to look for such cool pockets and 3D simulations including a local dynamo and a magnetic carpet to investigate Joule heating in the quiet chromosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The stellar atmosphere simulation code Bifrost

Context: Numerical simulations of stellar convection and photospheres have been developed to the point where detailed shapes of observed spectral lines can be explained. Stellar atmospheres are very complex, and very different physical regimes are present in the convection zone, photosphere, chromosphere, transition region and corona. To understand the details of the atmosphere it is necessary to simulate the whole atmosphere since the different layers interact strongly. These physical regimes are very diverse and it takes a highly efficient massively parallel numerical code to solve the associated equations. Aims: The design, implementation and validation of the massively parallel numerical code Bifrost for simulating stellar atmospheres from the convection zone to the corona. Methods: The code is subjected to a number of validation tests, among them the Sod shock tube test, the Orzag-Tang colliding shock test, boundary condition tests and tests of how the code treats magnetic field advection, chromospheric radiation, radiative transfer in an isothermal scattering atmosphere, hydrogen ionization and thermal conduction. Results: Bifrost completes the tests with good results and shows near linear efficiency scaling to thousands of computing cores.

0 citations0 reviewsNo code linkedOpen access