A novel fourth-order WENO interpolation technique

Abstract

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.

Intraoperative Beat-to-Beat Pulse Transit Time (PTT) Monitoring via Non-Invasive Piezoelectric/Piezocapacitive Peripheral Sensors Can Predict Changes in Invasively Acquired Blood Pressure in High-Risk Surgical Patients

“…(4) This procedure resulted in a PTT time series with missing values that were interpolated using the 4th order weighted essentially non-oscillatory (WENO4) technique [ 15 ].…”

The Lightweaver Framework for Nonlocal Thermal Equilibrium Radiative Transfer in Python

“…For interpolation duties, other than those in the formal solver and calculation of the PRD terms, we adopt the rapid, but robust fourth-order weighted essentially non-oscillatory approach presented in Janett et al (2019). While this technique does not guarantee monotonicity around discontinuities, the over-and under-shoots remain very small, with no ringing artifacts, and we feel that the high quality of the solution in smooth regions makes it worthwhile.…”