Volume 624, April 2019
|Number of page(s)||15|
|Section||Numerical methods and codes|
|Published online||18 April 2019|
A novel fourth-order WENO interpolation technique
A possible new tool designed for radiative transfer
Istituto Ricerche Solari Locarno (IRSOL), 6605 Locarno-Monti, Switzerland
2 Seminar for Applied Mathematics (SAM) ETHZ, D-MATH, 8093 Zurich, Switzerland
3 Kiepenheuer-Institut für Sonnenphysik (KIS), Schöneckstrasse 6, 79104 Freiburg i. Br., Germany
Accepted: 15 March 2019
Context. Several numerical problems require the interpolation of discrete data that present at the same time (i) complex smooth structures and (ii) various types of discontinuities. The radiative transfer in solar and stellar atmospheres is a typical example of such a problem. This calls for high-order well-behaved techniques that are able to interpolate both smooth and discontinuous data.
Aims. This article expands on different nonlinear interpolation techniques capable of guaranteeing high-order accuracy and handling discontinuities in an accurate and non-oscillatory fashion. The final aim is to propose new techniques which could be suitable for applications in the context of numerical radiative transfer.
Methods. We have proposed and tested two different techniques. Essentially non-oscillatory (ENO) techniques generate several candidate interpolations based on different substencils. The smoothest candidate interpolation is determined from a measure for the local smoothness, thereby enabling the essentially non-oscillatory property. Weighted ENO (WENO) techniques use a convex combination of all candidate substencils to obtain high-order accuracy in smooth regions while keeping the essentially non-oscillatory property. In particular, we have outlined and tested a novel well-performing fourth-order WENO interpolation technique for both uniform and nonuniform grids.
Results. Numerical tests prove that the fourth-order WENO interpolation guarantees fourth-order accuracy in smooth regions of the interpolated functions. In the presence of discontinuities, the fourth-order WENO interpolation enables the non-oscillatory property, avoiding oscillations. Unlike Bézier and monotonic high-order Hermite interpolations, it does not degenerate to a linear interpolation near smooth extrema of the interpolated function.
Conclusion. The novel fourth-order WENO interpolation guarantees high accuracy in smooth regions, while effectively handling discontinuities. This interpolation technique might be particularly suitable for several problems, including a number of radiative transfer applications such as multidimensional problems, multigrid methods, and formal solutions.
Key words: radiative transfer / methods: numerical
© ESO 2019
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.