Search Results

Now showing 1 - 2 of 2
  • Item
    Numerical solutions to linear transfer problems of polarized radiation III. Parallel preconditioned Krylov solver tailored for modeling PRD effects
    (Les Ulis : EDP Sciences, 2022) Benedusi, Pietro; Janett, Gioele; Riva, Simone; Krause, Rolf; Belluzzi, Luca
    Context. The polarization signals produced by the scattering of anistropic radiation in strong resonance lines encode important information about the elusive magnetic fields in the outer layers of the solar atmosphere. An accurate modeling of these signals is a very challenging problem from the computational point of view, in particular when partial frequency redistribution (PRD) effects in scattering processes are accounted for with a general angle-dependent treatment. Aims. We aim at solving the radiative transfer problem for polarized radiation in nonlocal thermodynamic equilibrium conditions, taking angle-dependent PRD effects into account. The problem is formulated for a two-level atomic model in the presence of arbitrary magnetic and bulk velocity fields. The polarization produced by scattering processes and the Zeeman effect is considered. Methods. The proposed solution strategy is based on an algebraic formulation of the problem and relies on a convenient physical assumption, which allows its linearization. We applied a nested matrix-free GMRES iterative method. Effective preconditioning is obtained in a multifidelity framework by considering the light-weight description of scattering processes in the limit of complete frequency redistribution (CRD). Results. Numerical experiments for a one-dimensional (1D) atmospheric model show near optimal strong and weak scaling of the proposed CRD-preconditioned GMRES method, which converges in few iterations, independently of the discretization parameters. A suitable parallelization strategy and high-performance computing tools lead to competitive run times, providing accurate solutions in a few minutes. Conclusions. The proposed solution strategy allows the fast systematic modeling of the scattering polarization signals of strong resonance lines, taking angle-dependent PRD effects into account together with the impact of arbitrary magnetic and bulk velocity fields. Almost optimal strong and weak scaling results suggest that this strategy is applicable to realistic 3D models. Moreover, the proposed strategy is general, and applications to more complex atomic models are possible.
  • Item
    Interpretation of Radio Wave Scintillation Observed through LOFAR Radio Telescopes
    (London : Institute of Physics Publ., 2022) Forte, Biagio; Fallows, Richard A.; Bisi, Mario M.; Zhang, Jinge; Krankowski, Andrzej; Dabrowski, Bartosz; Rothkaehl, Hanna; Vocks, Christian
    Radio waves propagating through a medium containing irregularities in the spatial distribution of the electron density develop fluctuations in their intensities and phases. In the case of radio waves emitted from astronomical objects, they propagate through electron density irregularities in the interstellar medium, the interplanetary medium, and Earth’s ionosphere. The LOFAR radio telescope, with stations across Europe, can measure intensity across the VHF radio band and thus intensity scintillation on the signals received from compact astronomical objects. Modeling intensity scintillation allows the estimate of various parameters of the propagation medium, for example, its drift velocity and its turbulent power spectrum. However, these estimates are based on the assumptions of ergodicity of the observed intensity fluctuations and, typically, of weak scattering. A case study of single-station LOFAR observations of the strong astronomical source Cassiopeia A in the VHF range is utilized to illustrate deviations from ergodicity, as well as the presence of both weak and strong scattering. Here it is demonstrated how these aspects can lead to misleading estimates of the propagation medium properties, for example, in the solar wind. This analysis provides a method to model errors in these estimates, which can be used in the characterization of both the interplanetary medium and Earth’s ionosphere. Although the discussion is limited to the case of the interplanetary medium and Earth’s ionosphere, its ideas are also applicable to the case of the interstellar medium.