A Spectral Solver to Capture Unsteady Dynamics in the Aerospike Nozzle Wake

A robust, high fidelity spectral solver for supersonic and hypersonic flow over an aerospike nozzle is presented. Numerical stability is achieved with a novel hybrid artificial viscosity (AV) - finite volume subcell element (FVSE) shock capturing method. The AV method is computational efficient. The diffusive FVSE method is locally called in elements where the AV method cannot sufficiently smooth the solution gradients, serving as a "parachute" in the context of numerical stability. The efficacy of this hybrid scheme is demonstrated on perfectly expanded Mach 2 and under expanded Mach 2 flow over a truncated aerospike and under expanded Mach 2 flow through a thruster module - aerospike nozzle assembly. The power spectral density of the unsteady pressure forcing on the nozzle is computed for each case and primary oscillation frequencies are identified. These oscillations are linked to flow patterns observed in Mach and schlieren images of the temporally evolving flow fields. To the authors' knowledge this is the first successful computation with a spectral solver of the unsteady, expanding, hypersonic flow over an aerospike from the initial thruster discharge until a quasi-steady state is achieved, as well as the first numerical investigation of the thruster module - aerospike nozzle assembly.