Compressibility profoundly affects many aspects of turbulence in high-speed flows - most notably stability characteristics, anisotropy, kinetic-potential energy interchange and spectral cascade rate. Many of the features observed in compressible flows are due to the changing nature of pressure. Whereas for incompressible flows pressure merely serves to enforce incompressibility, in compressible flows pressure becomes a thermodynamic variable that introduces a strong coupling between energy, state, and momentum equations. Closure models that attempt to address compressibility effects must begin their development from sound first-principles related to the changing nature of pressure as a flow goes from incompressible to compressible regime. In this thesis, a unified framework is developed for modeling pressure-related compressibility effects by characterizing the role and action of pressure at different speed regimes. Rapid distortion theory is used to examine the physical connection between the various compressibility effects leading to model form suggestions for the pressure-strain correlation, pressure-dilatation and dissipation evolution equation. The pressure-strain correlation closure coefficients are established using fixed point analysis by requiring consistency between model and direct numerical simulation asymptotic behavior in compressible homogeneous shear flow. The closure models are employed to compute high-speed mixing-layers and boundary layers in a differential Reynolds stress modeling solver. The self-similar mixing-layer profile, increased Reynolds stress anisotropy and diminished mixing-layer growth rates with increasing relative Mach number are all well captured. High-speed boundary layer results are also adequately replicated even without the use of advanced thermal-flux models or low Reynolds number corrections. To reduce the computational burden required for differential Reynolds stress calculations, the present compressible pressure-strain correlation model is incorporated into the algebraic modeling framework. The resulting closure is fully explicit, physically realizable, and is a function of mean flow strain rate, rotation rate, turbulent kinetic energy, dissipation rate, and gradient Mach number. The new algebraic model is validated with direct numerical simulations of homogeneous shear flow and experimental data of high-speed mixing-layers. Homogeneous shear flow calculations show that the model captures the asymptotic behavior of direct numerical simulations quite well. Calculations of plane supersonic mixing-layers are performed and comparison with experimental data shows good agreement. Therefore the algebraic model may serve as a surrogate for the more computationally expensive differential Reynolds stress model for flows that permit the weak-equilibrium simplification.

In this study, the behavior of pressure in the Rapid Distortion Limit, along with its concomitant modeling, are addressed. In the first part of the work, the role of pressure in the initiation, propagation and suppression of flow instabilities for quadratic flows is analyzed. The paradigm of analysis considers the Reynolds stress transport equations to govern the evolution of a dynamical system, in a state space composed of the Reynolds stress tensor components. This dynamical system is scrutinized via the identification of the invariant sets and the bifurcation analysis. The changing role of pressure in quadratic flows, viz. hyperbolic, shear and elliptic, is established mathematically and the underlying physics is explained. Along the maxim of "understanding before prediction", this allows for a deeper insight into the behavior of pressure, thus aiding in its modeling. The second part of this work deals with Rapid Pressure Strain Correlation modeling in earnest. Based on the comprehension developed in the preceding section, the classical pressure strain correlation modeling approaches are revisited. Their shortcomings, along with their successes, are articulated and explained, mathematically and from the viewpoint of the governing physics. Some of the salient issues addressed include, but are not limited to, the requisite nature of the model, viz. a linear or a nonlinear structure, the success of the extant models for hyperbolic flows, their inability to capture elliptic flows and the use of RDT simulations to validate models. Through this analysis, the schism between mathematical and physical guidelines and the engineering approach, at present, is substantiated. Subsequently, a model is developed that adheres to the classical modeling framework and shows excellent agreement with the RDT simulations. The performance of this model is compared to that of other nominations prevalent in engineering simulations. The work concludes with a summary, pertinent observations and recommendations for future research in the germane field.

Temporally developing direct numerical simulations (T-DNS) are performed and validated for bypass transition of a zero pressure gradient flat plate boundary layer to understand the interplay between pressure-strain terms and flow instability mechanisms, and to propose and validate a phenomenological hypothesis for the identification of a robust transition onset marker for use in transition-sensitive Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations. Results show that transition initiates at a location where the slow pressure-strain term becomes more dominant than the rapid term in the pre-transitional boundary layer region. The slow pressure strain term is responsible for the transfer of turbulence energy from the streamwise component to other components while the rapid pressure strain term counteracts with the slow term in the pre-transitional regime before transition onset akin to a shear sheltering like effect. The relative magnitudes of the slow and rapid terms thus provide a basis for the development of physically meaningful large-scale parameters that can be used as a transition onset marker for Reynolds averaged Navier-Stokes RANS simulations.

Three different topics of turbulence research that cover modeling, theory and model computation categories are selected and studied in depth. In the first topic, "velocity gradient dynamics in turbulence" (modeling), the Lagrangian linear diffusion model that accounts for the viscous-effect is proposed to make the existing restricted-Euler velocity gradient dynamics model quantitatively useful. Results show good agreement with DNS data. In the second topic, "pressure-strain correlation in homogeneous anisotropic turbulence subject to rapid strain-dominated distortion" (theory), extensive rapid distortion calculation is performed for various anisotropic initial turbulence conditions in strain-dominated mean flows. The behavior of the rapid pressure-strain correlation is investigated and constraining criteria for the rapid pressure-strain correlation models are developed. In the last topic, "unsteady computation of turbulent flow past a square cylinder using partially-averaged Navier-Stokes method" (model computation), the basic philosophy of the PANS method is reviewed and a practical problem of flow past a square cylinder is computed for various levels of physical resolution. It is revealed that the PANS method can capture many important unsteady flow features at an affordable computational effort.

Three different topics of turbulence research that cover modeling, theory and model computation categories are selected and studied in depth. In the first topic, "velocity gradient dynamics in turbulence" (modeling), the Lagrangian linear diffusion model that accounts for the viscous-effect is proposed to make the existing restricted-Euler velocity gradient dynamics model quantitatively useful. Results show good agreement with DNS data. In the second topic, "pressure-strain correlation in homogeneous anisotropic turbulence subject to rapid strain-dominated distortion" (theory), extensive rapid distortion calculation is performed for various anisotropic initial turbulence conditions in strain-dominated mean flows. The behavior of the rapid pressure-strain correlation is investigated and constraining criteria for the rapid pressure-strain correlation models are developed. In the last topic, "unsteady computation of turbulent flow past a square cylinder using partially-averaged Navier-Stokes method" (model computation), the basic philosophy of the PANS method is reviewed and a practical problem of flow past a square cylinder is computed for various levels of physical resolution. It is revealed that the PANS method can capture many important unsteady flow features at an affordable computational effort.