On the resolution of critical flow regions in inviscid linear and nonlinear instability calculations

Abstract

Numerical methods for tackling the inviscid instability problem are discussed. Convergence is demonstrated to be a necessary, but not a sufficient condition for accuracy. Inviscid flow physics set requirements regarding grid-point distribution in order for physically accurate results to be obtained. These requirements are relevant to the viscous problem also and are shown to be related to the resolution of the critical layers. In this respect, high-resolution nonlinear calculations based on the inviscid initial-boundary-value problem are presented for a model shear-layer flow, aiming at identification of the regions that require attention in the course of high-Reynolds-number viscous calculations. The results bear a remarkable resemblances with those pertinent to viscous flow, with a cascade of high-shear regions being shed towards the vortex-core centre as time progresses. In parallel, numerical instability related to the finite-time singularity of the nonlinear equations solved globally contaminates and eventually destroys the simulations, irrespective of resolution.