Direct Assessment and Investigation of Nonlinear and Nonlocal Turbulent Constitutive Relations in Three-Dimensional Boundary Layer Flow

Three-dimensional (3D) turbulent boundary layers (TBLs) play a crucial role in determining the aerodynamic properties of most aero-mechanical devices. However, accurately predicting these flows remains a challenge due to the complex nonlinear and nonlocal physics involved, which makes it difficult to develop universally applicable models. This limitation is particularly significant as the industry increasingly relies on simulations to make decisions in high-consequence environments, such as the certification or aircraft, and high-fidelity simulation methods that don't rely on modeling are prohibitively expensive. To address this challenge, it is essential to gain a better understanding of the physics underlying 3D TBLs. This research aims to improve the predictive accuracy of turbulence models in 3D TBLs by examining the impact of model assumptions underpinning turbulent constitutive relations, which are fundamental building blocks of every turbulence model. Specifically, the study focuses on the relevance and necessity of nonlinear and nonlocal model assumptions for accurately predicting 3D TBLs. The study considers the attached 3D boundary layer flow over the textbf{Be}nchmark textbf{V}alidation textbf{E}xperiment for textbf{R}ANS/textbf{L}ES textbf{I}nvestiagtions (BeVERLI) Hill as a test case and corresponding particle image velocimetry data for the investigation. In a first step, the BeVERLI Hill experiment is comprehensively described, and the important characteristics of the flow over the BeVERLI Hill are elucidated, including complex symmetry breaking characteristics of this flow. Reynolds-averaged Navier-Stokes simulations of the case using standard eddy viscosity models are then presented to establish the baseline behavior of local and linear constitutive relations, i.e., the standard Boussinesq approximation. The tested eddy viscosity models fail in the highly accelerated hill top region of the BeVERLI hill and near separation. In a further step, several nonlinear and nonlocal turbulent constitutive relations, including the QCR model, the model by Gatski and Speziale, and the difference-quotient model by Egolf are used as metrics to gauge the impact of nonlinearities and nonlocalities for the modeling of 3D TBLs. It is shown that nonlinear and nonlocal approaches are essential for effective 3D TBL modeling. However, simplified reduced-order models could accurately predict 3D TBLs without high computational costs. A constitutive relation with local second-order nonlinear mean strain relations and simplified nonlocal terms may provide such a minimal model. In a final step, the structure and response of non-equilibrium turbulence to continuous straining are studied to reveal new scaling laws and structural models.