The linearized pressure Poisson equation for global instability analysis of incompressible flows

Abstract

The linearized pressure Poisson equation (LPPE) is used in two and three spatial dimensions in the respective matrix-forming solution of the BiGlobal and TriGlobal eigenvalue problem in primitive variables on collocated grids. It provides a disturbance pressure boundary condition which is compatible with the recovery of perturbation velocity components that satisfy exactly the linearized continuity equation. The LPPE is employed to analyze instability in wall-bounded flows and in the prototype open Blasius boundary layer flow. In the closed flows, excellent agreement is shown between results of the LPPE and those of global linear instability analyses based on the time-stepping nektar++, Semtex and nek5000 codes, as well as with those obtained from the FreeFEM++ matrix-forming code. In the flat plate boundary layer, solutions extracted from the two-dimensional LPPE eigenvector at constant streamwise locations are found to be in very good agreement with profiles delivered by the NOLOT/PSE space marching code. Benchmark eigenvalue data are provided in all flows analyzed. The performance of the LPPE is seen to be superior to that of the commonly used pressure compatibility (PC) boundary condition: at any given resolution, the discrete part of the LPPE eigenspectrum contains converged and not converged, but physically correct, eigenvalues. By contrast, the PC boundary closure delivers some of the LPPE eigenvalues and, in addition, physically wrong eigenmodes. It is concluded that the LPPE should be used in place of the PC pressure boundary closure, when BiGlobal or TriGlobal eigenvalue problems are solved in primitive variables by the matrix-forming approach on collocated grids.Keywords Global linear instability · Matrix-forming · Collocated grids · Pressure boundary conditions IntroductionGlobal linear instability theory is enjoying increasing acceptance in the fluid mechanics community, as also witnessed by the contents of this volume. The theory has permitted revealing previously unknown mechanisms responsible for laminar-turbulent flow transition in complex three-dimensional geometries with maximally one homogeneous spatial direction. Two main approaches are followed for the numerical work associated with modal and non-modal global linear instability analysis, namely matrix-forming vs. time-stepping, both of which have been discussed in a recent review [61]. The matrix-forming approach is interesting when dealing Communicated by Ati

Localised eigenmodes in a moving frame of reference representing convective instability

“…The only physically substantiated boundary condition for external incompressible flow problems is given at solid walls. There no-slip is imposed for ũ and ṽ and the Laplace equation for p, consistent with the conditions proposed by Theofilis [13].…”

Transient growth analysis of hypersonic flow over an elliptic cone

“…As shown in a number of studies in the incompressible limit (e.g. Alizard & Robinet 2007; Theofilis 2017), the global linear stability equations and the linear PSE lead to identical results when initialized with consistent inflow boundary conditions for the flow perturbations.…”

Optimal triggering of jet bifurcation: an example of optimal forcing applied to a time-periodic base flow

“…The stress-free conditions (2.19), (2.20), consistent with the Nek5000 boundary treatment (2.8), arise as the 'natural boundary conditions' (see § 10.2.4 in Dick 2009) in this finite-element formulation. Their quality as non-reflecting outflow conditions have been examined by Theofilis (2017) and by Lesshafft (2018).…”