Search

In this work, a non-equilibrium wall model is proposed for the prediction of turbulent flows experiencing adverse pressure gradients, including separated flow regimes. The mean-flow nonequilibrium is identified by comparing two characteristic velocities: the friction velocity (u_tau) and the viscous-pressure gradient velocity (up). In regions where the pressure gradient velocity is comparable to the friction velocity (up \sim u_tau, the near-wall turbulent closure is modified to include the effect of the pressure-gradient and convective terms. The performance of this wall model is evaluated in two canonical flows experiencing smooth body separation: the NASA-Boeing speed bump and the Bachalo-Johnson bump. Improvements in the predictive capabilities of the proposed model for the conventional equilibrium wall model are theorized and then demonstrated through numerical experiments. In particular, the proposed wall model can capture the onset of boundary layer separation observed in experiments or DNS calculations at resolutions where the equilibrium wall model fails to separate., Comment: Under review for publication in Physical Review Fluids. 20 pages, 14 figures

This article utilizes the Large-Eddy Simulation (LES) paradigm with a physics-based turbulence modeling approach, including a dynamic subgrid-scale model and an equilibrium wall model, to examine the flow over the NASA transonic Common Research Model (CRM), a flow configuration that has been the focus of several AIAA Drag PredictionWorkshops (DPWs). The current work explores sensitivities to laminar-to-turbulent transition, wind tunnel mounting system, grid resolution, and grid topology and suggests current best practices in the context of large-eddy simulations of transonic aircraft flows. It is found that promoting the flow transition to turbulence via an array of cylindrical trip dots, including the sting mounting system in the simulations, and leveraging stranded boundary layer grids all tend to improve the quality of the LES solutions. Non-monotonic grid convergence in the LES calculations is observed to be strongly sensitive to grid topology, and stranded meshes rectify this issue relative to their hexagonal close-packed (HCP) counterparts. The details of the boundary layer profiles both at the leading edge of the wing and within the shock-induced separation bubble are studied, with thicknesses and integral measures reported, providing details about the boundary layer characteristics to turbulence modelers not typically available from complex aircraft flows. Finally, an assessment of the initial buffet prediction capabilities of LES is made in the context of a simpler NACA 0012 flow, with computational predictions showing reasonable agreement with available experimental data for the angle of attack at initial buffet onset and shock oscillation frequency associated with sustained buffet., Comment: Under review for publication consideration in the Journal of Aircraft

This article proposes a Reynolds number scaling of the required grid points to perform wall-modeled LES of turbulent flows encountering separation off a solid surface. Based on comparisons between the various time scales in a non-equilibrium (due to the action of an external pressure gradient) turbulent boundary layer, a simple definition of the near-wall ``under-equilibrium" and ``out-of-equilibrium" scales is put forward (where ``under-equilibrium" refers to scales governed by a quasi-balance between the viscous and the pressure gradient terms). It is shown that the former length scale varies with Reynolds number as lp Re^(-2/3). The same scaling is obtained from a simplified Green's function solution of the Poisson equation in the vicinity of the separation point. A-priori analysis demonstrates that the resolution required to reasonably predict the wall-shear stress (for example, errors lower than approximately 10-15% in the entire domain) in several nonequilibrium flows is at least O(10) lp irrespective of the Reynolds number and the Clauser parameter. Further, a series of a-posteriori validation studies are performed to determine the accuracy of this scaling including the flow over the Boeing speed bump, Song-Eaton diffuser, Notre-Dame Ramp, and the backward-facing step. The results suggest that for these flows, scaling the computational grids () such that / lp is independent of the Reynolds number results in accurate predictions of flow separation at the same ``nominal" grid resolution across different Reynolds numbers. Finally, it is suggested that in the vicinity of the separation and reattachment points, the grid-point requirements for wall-modeled large eddy simulations may scale as Re^4/3, which is more restrictive than the previously proposed flat-plate boundary layer-based estimates (Re1) of Choi and Moin (Phys. Fluids, 2012) and Yang and Griffin (Phys. Fluids, 2021)., Comment: 21 pages, 17 figures. Submitted to AIAA Journal for publication consideration

Thwaites (1949) developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of a flow parameter, $m = \theta^2/\nu \frac{dU_e}{ds}$, thus reducing the von Karman momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. However, for turbulent flows, a similar approximation for turbulent boundary layers subjected to pressure gradients does not yet exist. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed. It is shown that the method provides good estimates of the momentum thickness, when compared to available high-fidelity simulation data, for multiple boundary layers including both favorable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model which is shown to be in agreement with prior empirical observations (Alber, \textit{$9^{th}$ Aerospace Sciences Meeting, 1971}). The sensitivity of the separation location with respect to upstream perturbations is also analyzed through this model for the NASA/Boeing speed bump and the transonic Bachalo-Johnson bump, Comment: 21 pages, 13 figures. Under consideration for publication in J. Fluid Mech

In this study, five ice shapes generated at NASA Glenn's Icing Research Tunnel (IRT) are simulated at multiple angles of attack (Broeren et al., J. of Aircraft, 2018). These geometries target different icing environments, both early-time and longer-duration glaze and rime ice exposure events, including a geometry that results from using a thermal ice-protection system. Using the laser-scanned geometries, detailed representations of the three-dimensional ice geometries are resolved on the grid and simulated using wall-modeled LES. Integrated loads (lift, drag, and moment coefficients) and pressure distributions are compared against experimental measurements in both clean and iced conditions for several angles of attack in both pre-and post-stall regions. The relevant comparisons to the experimental results show that qualitative and acceptable quantitative agreement with the data is observed across all geometries. Glaze ice formations exhibit larger and highly nonuniform ice features, such as `horns', in contrast to rime ice formations characterized by smaller, uniformly distributed roughness elements. In wall-modeled LES, it was observed that larger roughness scales in the glaze ice that trigger transition can be accurately resolved. Therefore, it is possible for WMLES to accurately capture the aerodynamics of glaze ice shapes without the need for additional modeling. In contrast, rime ice geometries required additional resolution to accurately represent the aerodynamic loads. This study demonstrates the effectiveness of the wall-modeled LES technique in simulating the complex aerodynamic effects of iced airfoils, providing valuable insights for aircraft design in icing environments and highlighting the importance of accurately representing ice geometries and roughness scales in simulations.

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data has shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [Bardina et al., AIAA 1980; Meneveau and Liu, Ann. Rev. Fluid Mech. 2002]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g.: Bardina et al., AIAA 1980; Clark et al. J. Fluid Mech., 1979; Stolz and Adams, Phys. of Fluids, 1999]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for LES, vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence (HIT), and wall-modeled turbulent channel flow at high Reynolds numbers; the results show favorable agreement with DNS data. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to the existing static coefficient (Vreman) and dynamic Smagorinsky model., Comment: Revised Manuscript, under Review, Physical Review Fluids

We investigate the error scaling and computational cost of wall-modeled large-eddy simulation (WMLES) for external aerodynamic applications. The NASA Juncture Flow is used as representative of an aircraft with trailing-edge smooth-body separation. Two gridding strategies are examined: i) constant-size grid, in which the near-wall grid size has a constant value and ii) boundary-layer-conforming grid (BL-conforming grid), in which the grid size varies to accommodate the growth of the boundary-layer thickness. Our results are accompanied by a theoretical analysis of the cost and expected error scaling for the mean pressure coefficient ($C_p$) and mean velocity profiles. The prediction of $C_p$ is within less than $5\%$ error for all the grids studied, even when the boundary layers are marginally resolved. The high accuracy in the prediction of $C_p$ is attributed to the outer-layer nature of the mean pressure in attached flows. The errors in the predicted mean velocity profiles exhibit a large variability depending on the location considered, namely, fuselage, wing-body juncture, or separated trailing-edge. WMLES performs as expected in regions where the flow resembles a zero-pressure-gradient turbulent boundary layer such as the fuselage ($<5\%$ error). However, there is a decline in accuracy of WMLES predictions of mean velocities in the vicinity of wing-body junctions and, more acutely, in separated zones. The impact of the propagation of errors from the underresolved wing leading-edge is also investigated. It is shown that BL-conforming grids enable a higher accuracy in wing-body junctions and separated regions due to the more effective distribution of grid points, which in turn diminishes the streamwise propagation of errors., Comment: arXiv admin note: text overlap with arXiv:2101.00331

We investigate the performance of wall-modeled LES for external aerodynamics in the NASA Juncture Flow. We characterize the errors in the prediction of mean velocity profiles and pressure coefficient for three different locations over the aircraft: the upstream region of the fuselage, the wing-body juncture, and the wing-body juncture close to the trailing-edge. The last two locations are characterized by strong mean-flow three-dimensionality and separation. The message conveyed by our error analysis is that WMLES performs as expected in regions where the flow resembles a zero-pressure-gradient flat plate boundary layer. However, there is a clear decline of the current models in the presence of wing-body junctions and, more acutely, in separated zones. The slow convergence to the solution in these regions renders the brute-force grid-refinement approach to improve the accuracy of the solution unfeasible. The results reported above pertain to the mean velocity profile predicted using the typical grid resolution for external aerodynamics applications, i.e., 5--20 points per boundary-layer thickness. The prediction of the pressure coefficient is below 5\% error for all grid sizes considered, even when boundary layers were marginally resolved. The latter accuracy can be attributed to the outer-layer nature of the mean pressure, which becomes less sensitive to flow details within the turbulent boundary layer. Finally, we show that boundary-layer-conforming grids (i.e., grids maintaining a constant number of points per boundary-layer thickness) allow for a more efficient distribution of grid points and smaller errors. Our results suggest that novel modeling venues encompassing physical insights, together with numerical and gridding advancements, must be exploited to attain predictions within the tolerance required for Certification by Analysis.

The interaction between an incident shock wave and a Mach-6 undisturbed hypersonic laminar boundary layer over a cold wall is addressed using direct numerical simulations (DNS) and wall-modeled large-eddy simulations (WMLES) at different angles of incidence. At sufficiently high shock-incidence angles, the boundary layer transitions to turbulence via breakdown of near-wall streaks shortly downstream of the shock impingement, without the need of any inflow free-stream disturbances. The transition causes a localized significant increase in the Stanton number and skin-friction coefficient, with high incidence angles augmenting the peak thermomechanical loads in an approximately linear way. Statistical analyses of the boundary layer downstream of the interaction for each case are provided that quantify streamwise spatial variations of the Reynolds analogy factors and indicate a breakdown of the Morkovin's hypothesis near the wall, where velocity and temperature become correlated. A modified strong Reynolds analogy with a fixed turbulent Prandtl number is observed to perform best. Conventional transformations fail at collapsing the mean velocity profiles on the incompressible log law. The WMLES prompts transition and peak heating, delays separation, and advances reattachment, thereby shortening the separation bubble. When the shock leads to transition, WMLES provides predictions of DNS peak thermomechanical loads within $\pm 10\%$ at a computational cost lower than DNS by two orders of magnitude. Downstream of the interaction, in the turbulent boundary layer, WMLES agrees well with DNS results for the Reynolds analogy factor, the mean profiles of velocity and temperature, including the temperature peak, and the temperature/velocity correlation., Comment: 48 pages, 36 figures

The accurate prediction of aerothermal surface loading is of paramount importance for the design of high speed flight vehicles. In this work, we consider the numerical solution of hypersonic flow over a double-finned geometry, representative of the inlet of an air-breathing flight vehicle, characterized by three-dimensional intersecting shock-wave/turbulent boundary-layer interaction at Mach 8.3. High Reynolds numbers ($Re_L \approx 11.6 \times 10^6$ based on free-stream conditions) and the presence of cold walls ($T_w/T_o \approx 0.27$) leading to large near-wall temperature gradients necessitate the use of wall-modeled large-eddy simulation (WMLES) in order to make calculations computationally tractable. The comparison of the WMLES results with experimental measurements shows good agreement in the time-averaged surface heat flux and wall pressure distributions, and the WMLES predictions show reduced errors with respect to the experimental measurements than prior RANS calculations. The favorable comparisons are obtained using an LES wall model based on equilibrium boundary layer approximations despite the presence of numerous non-equilibrium conditions including three dimensionality, shock-boundary layer interactions, and flow separation. Lastly, it is also demonstrated that the use of semi-local eddy viscosity scaling (in lieu of the commonly used van Driest scaling) in the LES wall model is necessary to accurately predict the surface pressure loading and heat fluxes., Comment: 27 pages, 27 figures

Wall modelling in large-eddy simulation (LES) is necessary to overcome the prohibitive near-wall resolution requirements in high-Reynolds-number turbulent flows. Most existing wall models rely on assumptions about the state of the boundary layer and require a priori prescription of tunable coefficients. They also impose the predicted wall stress by replacing the no-slip boundary condition at the wall with a Neumann boundary condition in the wall-parallel directions while maintaining the no-transpiration condition in the wall-normal direction. In the present study, we first motivate and analyse the Robin (slip) boundary condition with transpiration (nonzero wall-normal velocity) in the context of wall-modelled LES. The effect of the slip boundary condition on the one-point statistics of the flow is investigated in LES of turbulent channel flow and flat-plate turbulent boundary layer. It is shown that the slip condition provides a framework to compensate for the deficit or excess of mean momentum at the wall. Moreover, the resulting nonzero stress at the wall alleviates the well-known problem of the wall-stress under-estimation by current subgrid-scale (SGS) models. Secondly, we discuss the requirements for the slip condition to be used in conjunction with wall models and derive the equation that connects the slip boundary condition with the stress at the wall. Finally, a dynamic procedure for the slip coefficients is formulated, providing a dynamic slip wall model free of a priori specified coefficients. The performance of the proposed dynamic wall model is tested in a series of LES of turbulent channel flow, non-equilibrium three-dimensional channel flow, and flat-plate turbulent boundary layer. The results show that the dynamic wall model is able to accurately predict one-point turbulence statistics for various flow configurations, Reynolds numbers, and grid resolutions.

One of the greatest challenges to using large-eddy simulations (LES) in engineering applications is the large number of grid points required near walls. To mitigate this issue, researchers often couple LES with a simplified model of the flow close to the wall, known as the wall model. One feature common to most wall models is that the first few (about 3) grid points must be located below the inviscid log-layer, and the grid must have near-isotropic resolution near the wall. Hence, wall-modeled LES may still require a large number of grid points in both the wall-normal and span-wise directions. Because of these requirements, wall-modeled LES is still unfeasible in many applications. We present a new formulation of wall-modeled LES that is being developed to address this issue. In this formulation, LES is used to solve only for the features of the velocity field that can be adequately represented on the LES grid. The effects of the unresolved features are captured by imposing an integral balance of kinetic energy in the near-wall region. This integral energy balance translates into a dynamic partial differential equation defined on the walls, which is coupled to the LES equations. We discuss details of the new formulation and present results obtained in laminar channel flows., Comment: 10 pages, 7 figures, final report of the 16th biennial summer program of the Center for Turbulence Research at Stanford University

We describe and characterize a method for estimating the pressure field corresponding to velocity field measurements, such as those obtained by using particle image velocimetry. The pressure gradient is estimated from a time series of velocity fields for unsteady calculations or from a single velocity field for quasi-steady calculations. The corresponding pressure field is determined based on median polling of several integration paths through the pressure gradient field in order to reduce the effect of measurement errors that accumulate along individual integration paths. Integration paths are restricted to the nodes of the measured velocity field, thereby eliminating the need for measurement interpolation during this step and significantly reducing the computational cost of the algorithm relative to previous approaches. The method is validated by using numerically-simulated flow past a stationary, two-dimensional bluff body and a computational model of a three-dimensional, self-propelled anguilliform swimmer to study the effects of spatial and temporal resolution, domain size, signal-to-noise ratio, and out-of-plane effects. Particle image velocimetry measurements of a freely-swimming jellyfish medusa and a freely-swimming lamprey are analyzed using the method to demonstrate the efficacy of the approach when applied to empirical data., Comment: A free MATLAB implementation of this algorithm is available at http://dabiri.caltech.edu/software.html