Comparison of very high order (O3–18) flux and Crowley flux, and (O3–17) WENO flux schemes with a 2D nonlinear test problem

Two‐dimensional, nonlinear diffusion‐limited colliding plumes simulations were used to demonstrate the improved solution accuracy with very high order O9–18 flux schemes, including upwind‐biased and even‐centred constant grid flux and Crowley constant grid flux schemes, and odd order weighted essentially non‐oscillatory (WENO) flux schemes, along with variations and hybrids of these. All schemes were coupled with comparably high order even‐centred Lagrangian interpolations and pressure gradient/divergence approximations, and O18 spatial filtering. Subgrid‐scale (SGS) turbulent flux calculations, with a constant eddy‐mixing coefficient, were made with O2 spatial approximations (O4–20 accurate SGS turbulent fluxes had little impact). Using a range of resolutions from Δx = Δz = 25–166.66… m for all schemes comparisons against an O17 flux, 25 m resolution reference solution showed solutions made with ≥O9 fluxes produced (often substantially) improved solutions, both visually and usually objectively, compared to solutions produced with lower order (