Numerical non-LTE 3D radiative transfer using a multigrid method

Abstract

3D non-LTE radiative transfer problems are computationally demanding, and this sets limits on the size of the problems that can be solved. So far Multilevel Accelerated Lambda Iteration (MALI) has been to the method of choice to perform high-resolution computations in multidimensional problems. The disadvantage of MALI is that its computing time scales as $\mathcal{O}(n^2)$, with $n$ the number of grid points. When the grid gets finer, the computational cost increases quadratically. We aim to develop a 3D non-LTE radiative transfer code that is more efficient than MALI. We implement a non-linear multigrid, fast approximation storage scheme, into the existing Multi3D radiative transfer code. We verify our multigrid implementation by comparing with MALI computations. We show that multigrid can be employed in realistic problems with snapshots from 3D radiative-MHD simulations as input atmospheres. With multigrid, we obtain a factor 3.3-4.5 speedup compared to MALI. With full-multigrid the speed-up increases to a factor 6. The speedup is expected to increase for input atmospheres with more grid points and finer grid spacing. Solving 3D non-LTE radiative transfer problems using non-linear multigrid methods can be applied to realistic atmospheres with a substantial speed-up.

Accepted for publication by A&A

References in corpus (27)

• Line formation in solar granulation: IV. [O I], OI and OH lines and the photospheric O abundance
• The stellar atmosphere simulation code Bifrost
• The formation of the Halpha line in the solar chromosphere
• The formation of IRIS diagnostics II. The formation of the Mg II h&k lines in the solar atmosphere
• Approximations for radiative cooling and heating in the solar chromosphere
• A publicly available simulation of an enhanced network region of the Sun
• The formation of IRIS diagnostics I. A quintessential model atom of Mg II and general formation properties of the Mg II h&k lines
• Non-LTE oxygen line formation in 3D hydrodynamic model stellar atmospheres
• The Formation of IRIS diagnostics. IV. The Mg II triplet lines as a new diagnostic for lower chromospheric heating
• Multi-level 3D non-LTE computations of lithium lines in the metal-poor halo stars HD140283 and HD84937
• Three-dimensional non-LTE radiative transfer computation of the Ca 8542 infrared line from a radiation-MHD simulation
• 3D NLTE analysis of the most iron-deficient star, SMSS0313-6708
• PORTA: A three-dimensional multilevel radiative transfer code for modeling the intensity and polarization of spectral lines with massively parallel computers
• Oxygen lines in solar granulation. II. Centre-to-limb variation, NLTE line formation, blends and the solar oxygen abundance
• The formation of IRIS diagnostics V. A quintessential model atom of C II and general formation properties of the C II lines at 133.5 nm
• Millimeter radiation from a 3D model of the solar atmosphere I. Diagnosing chromospheric thermal structure
• Partial redistribution in 3D non-LTE radiative transfer in solar atmosphere models
• IRIS: A Generic Three-Dimensional Radiative Transfer Code
• The formation of IRIS diagnostics VI. The Diagnostic Potential of the C II Lines at 133.5 nm in the Solar Atmosphere
• Heating of the magnetic chromosphere: observational constraints from Ca II 8542 spectra
• 3-D non-LTE radiative transfer effects in Fe I lines: II. Line formation in 3-D radiation hydrodynamic simulations
• Three-dimensional radiative transfer simulations of the scattering polarization of the hydrogen Ly$\alpha$ line in a MHD model of the chromosphere-corona transition region
• The Scattering Polarization of the Sr I 4607 ÅLine at the Diffraction Limit Resolution of a 1-m Telescope
• 3-D non-LTE radiative transfer effects in Fe I lines: III. Line formation in magneto-hydrodynamic atmospheres
• Fast 2D non-LTE radiative modelling of prominences I. Numerical methods and benchmark results
• Photon mean free paths, scattering, and ever-increasing telescope resolution
• A conjugate gradient method for the solution of the non-LTE line radiation transfer problem

Cited by in corpus (9)

• Observational constraints on the origin of the elements. I. 3D NLTE formation of Mn lines in late-type stars
• An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates
• Study of the polarization produced by the Zeeman effect in the solar Mg I b lines
• The dynamic chromosphere: pushing the boundaries of observations and models
• Radiance Cascades: A Novel High-Resolution Formal Solution for Multidimensional Non-LTE Radiative Transfer
• Solar Flux Atlases: The new HARPS-N Quiet Sun Benchmark and Ca II H & K Lines Continuum Normalisation
• EMISSA: Exploring millimetre indicators of solar-stellar activity III. Comparison of Ca II indices and millimetre continua in a 3D model atmosphere
• Spherically symmetric model atmospheres using approximate lambda operators V. Static inhomogeneous atmospheres of hot dwarf stars
• Irregular grids for 3D NLTE radiative transfer in stellar atmospheres