Stable and Efficient One-Way Modelling of Convective Disturbances in Laminar Boundary Layers: OWNS-Summation

Abstract

One-way spatial marching methods separate upstream- from downstream-propagating disturbances using a rational approximation of a spectral projector. Among existing one-way Navier--Stokes (OWNS) formulations, the recursive variant, OWNS-R, is the most economical, evaluating the approximation as a product of $N$ resolvent factors. This product amplifies rounding errors multiplicatively, imposing a flow-dependent upper limit on the approximation order that cannot be known in advance. We reformulate the same approximation as an additive partial-fraction sum (OWNS-Summation, OWNS-S). In exact arithmetic, the recursive and summation evaluations are equivalent when supplied with identical weights and poles; in floating-point arithmetic, they are not. Each of the $N$ resolvent solves acts on the same input state and contributes independently to a weighted sum, preventing multiplicative error amplification. The approximation order $N$ therefore becomes a pure convergence parameter, and the solves can run in parallel. Using the same auxiliary poles as OWNS-R, OWNS-S remains accurate in every configuration tested, showing that the recursive evaluation, rather than the poles, is the dominant source of instability. A paired greedy parameter-selection procedure is also introduced, with candidates drawn from analytic estimates of the upstream and downstream spectral regions, avoiding eigen-decomposition during the numerical march. OWNS-S is validated on incompressible, hypersonic and transonic boundary layers. In the transonic case, the disturbance spectrum reorganises from a subsonic to a supersonic topology during the march. The one-way computation proceeds continuously through this transition at a streamwise resolution unattainable by the parabolised stability equations.