In this study, the biglobal stability approach, which is frequently used to characterize the hydrodynamic instability waves associated with parietal vortex shedding, is applied to a simulated rocket chamber with variable headwall injection. The analysis is motivated by the need to assess the viability of three dissimilar outflow boundary conditions that are relevant to the modeling of cylindrically-shaped solid and hybrid rocket chambers with headwall mass injection. To this end, three sets of boundary conditions are considered and discussed: first, an acoustically-closed chamber is assumed; second, an adjoint-compliant method is relied upon; and, third, a continual outflow requirement is imposed. The first set of boundary conditions is ubiquitously used in conjunction with the Local Nonparallel (LNP) formulation in rocket-flow problems that entail a choked chamber endwall (e.g., a nozzle throat). By suppressing velocity fluctuations, using axial symmetry, and setting the normal stress tensor at the exit boundary equal to zero, one arrives at a second unique set of conditions. These so-called adjoint-compliant requirements lead to a well-posed adjointconforming problem, which has been shown to be essential for simulations based on the finite element method. The final set of constraints is precipitated by the need to secure flow continuation at the exit, which may be achieved through linear extrapolation. As with previous implementations of the biglobal stability framework, we do not limit ourselves to a stream function formulation. Instead we consider the linearized Navier-Stokes equations and solve the resulting eigenvalue problem using pseudo-spectral methods. Although the stability spectrum for a solid rocket motor configuration with a non-injecting headwall may be recovered from our formulation as a special case, our main discussion concentrates on the hybrid rocket flow configuration with an order of magnitude difference between the headwall and sidewall injection speeds. Lastly, comparisons to previous investigations that rely on the stream function approach and the LNP technique are provided and discussed.
Nomenclatureij A = operator matrix a = chamber radius ij B = the right-hand-side coefficient matrix of a matrix pencil N D = Chebyshev pseudo-spectral derivative matrix of size N d = weight coefficients for pseudo-spectral derivative matrices N I = identity matrix of size N L = chamber length M = base flow component M = instantaneous flow component m = general amplitude function m = acoustic fluctuation 2 m = unsteady hydrodynamic fluctuation m = vortical fluctuation = Landau order symbol P = base flow pressure p = pressure amplitude function p = hydrodynamic pressure fluctuation q = azimuthal integer wave number , rz = normalized radial and axial coordinates Re = Reynolds number N T = Chebyshev polynomial of the first kind U = base flow velocity, r U ,U , z U u = velocity amplitude function u = hydrodynamic velocity fluctuation Greek = gradient operator = 1/ Re = streamwise Chebyshev variables m...