Wall-resolved Large Eddy Simulation Research Articles

Large-Eddy Simulation and Aeroacoustic Prediction of Supercritical Airfoil Side-Edge Noise

Airfoil self-noise is investigated on a cantilever wing with a supercritical profile with a compressible wall-resolved large-eddy simulation. The Reynolds number based on the chord is 620,000, and the angle of attack is 5°. The aerodynamic results reveal the development of a complex vortex system at the side edge, including primary, secondary, and tertiary vortices that govern aerodynamic noise production. The wall-pressure coefficient and root-mean-square pressure coefficient contours highlight the side-edge shear layer and flow impingement of the primary vortex at the pressure-side edge as important noise-generation mechanisms. Dynamic mode decomposition shows that the dominant pressure modes on the airfoil exist along primary and secondary vortex impingement lines. Far-field acoustic predictions based on both solid and porous Ffowcs-Willliams and Hawkings’ analogies demonstrate good agreement with experimental results between 1.5 and 9 kHz. Low-frequency tonal humps observed in the spectra result from duct acoustic modes excited by airfoil self-noise, a consequence of the installation effect. For this configuration, airfoil side-edge noise did not dominate over airfoil surface noise levels, suggesting that at low to moderate angles of attack, airfoil side-edge noise is not a significant contributor to the overall acoustic emissions.

Continuous Eddy Simulation (CES) of Transonic Shock-Induced Flow Separation

Reynolds-averaged Navier–Stokes (RANS), large eddy simulation (LES), and hybrid RANS-LES, first of all wall-modeled LES (WMLES) and detached eddy simulation (DES) methods, are regularly applied for wall-bounded turbulent flow simulations. Their characteristic advantages and disadvantages are well known: significant challenges arise from simulation performance, computational cost, and functionality issues. This paper describes the application of a new simulation approach: continuous eddy simulation (CES). CES is based on exact mathematics, and it is a minimal error method. Its functionality is different from currently applied simulation concepts. Knowledge of the actual amount of flow resolution enables the model to properly adjust to simulations by increasing or decreasing its contribution. The flow considered is a high Reynolds number complex flow, the Bachalo–Johnson axisymmetric transonic bump flow, which is often applied to evaluate the performance of turbulence models. A thorough analysis of simulation performance, computational cost, and functionality features of the CES model applied is presented in comparison with corresponding features of RANS, DES, WMLES, and wall-resolved LES (WRLES). We conclude that CES performs better than RANS, DES, WMLES, and even WRLES at a little fraction of computational cost applied for the latter methods. CES is independent of usual functionality requirements of other methods, which offers relevant additional advantages.

Sensitivity analysis of wall-modeled large-eddy simulation for separated turbulent flow

In this study, we conduct a parametric analysis to evaluate the sensitivities of wall-modeled large-eddy simulation (LES) with respect to subgrid-scale (SGS) models, mesh resolution, wall boundary conditions and mesh anisotropy. While such investigations have been conducted for attached/flat-plate flow configurations, systematic studies specifically targeting turbulent flows with separation are notably sparse. To bridge this gap, our study focuses on the flow over a two-dimensional Gaussian-shaped bump at a moderately high Reynolds number, which involves smooth-body separation of a turbulent boundary layer under pressure-gradient and surface-curvature effects. In the simulations, the no-slip condition at the wall is replaced by three different forms of boundary condition based on the thin boundary layer equations and the mean wall-shear stress from high-fidelity numerical simulation to avoid the additional complexity of modeling the wall-shear stress. Various statistics, including the mean separation bubble size, mean velocity profile, and dissipation from SGS model, are compared and analyzed. The results reveal that capturing the separation bubble strongly depends on the choice of SGS model. While simulations approach grid convergence with resolutions nearing those of wall-resolved LES meshes, above this limit, the LES predictions exhibit intricate sensitivities to mesh resolution. Furthermore, both wall boundary conditions and the anisotropy of mesh cells exert discernible impacts on the turbulent flow predictions, yet the magnitudes of these impacts vary based on the specific SGS model chosen for the simulation.

Improved self-adaptive turbulence eddy simulation for complex flows and stall prediction using high-order schemes

It is challenging to apply numerical simulations to accurately predict the stall behavior of aircraft equipped with high-lift devices. Simulations with Reynolds-Averaged NavierStokes (RANS) models suffer from lack of the reliability at high angles of attack with separated and reattached boundary layers, whereas wall-resolved Large Eddy Simulations (LES) of wall-bounded flows at high Reynolds numbers currently costs too much computational resources. A new unified hybrid turbulence modeling approach, denoted the Self-Adaptive Turbulence Eddy Simulation (SATES), is proposed and applied for complex turbulent flows combining with high-order numerical scheme of the Weighted Compact Nonlinear Scheme (WCNS) in the present study. It enables a seamless evolution from unsteady RANS to LES and finally approaches Direct Numerical Simulation (DNS) depending on the turbulent scales. In the framework of SATES, a new SATES-σ model with an adaptive model coefficient is developed by extending the underlying LES mode based on an enhanced sub-grid-scale model of the σ-model. The new SATES-σ is first examined in two benchmark cases of channel flow and flow past a square cylinder. Then, it is validated in supercritical flow past a circular cylinder to assess the performance of turbulent models. The results show significant improvements over the previous SATES and IDDES in the predictions of boundary layer flow. Finally, successful application is achieved in the accurate prediction of the stall of the MD-30P30N airfoil at a Reynolds number of 9×106 with wide angles of attack. The simulation results show good agreement with experimental results for surface pressure even for the challenging cases of 21 and 23 deg angles of attack. Again, the SATES-σ shows better results than the previous SATES and IDDES. The presented method has considerable potential for the challenging stall predictions.

Assessment of Machine Learning Wall Modeling Approaches for Large Eddy Simulation of Gas Turbine Film Cooling Flows: An a Priori Study

Abstract In this work, a priori analysis of machine learning (ML) strategies is carried out with the goal of data-driven wall modeling for large eddy simulation (LES) of gas turbine film cooling flows. High-fidelity flow datasets are extracted from wall-resolved LES (WRLES) of flow over a flat plate interacting with the coolant flow supplied by a single row of 7-7-7 shaped cooling holes inclined at 30 degrees with the flat plate at different blowing ratios (BR). The WRLES are performed using the high-order Nek5000 spectral element computational fluid dynamics (CFD) solver. Light gradient boosting machine (LightGBM) is employed as the ML algorithm for the data-driven wall model. Parametric tests are conducted to systematically assess the influence of a wide range of input flow features (velocity components, velocity gradients, pressure gradients, and fluid properties) on the accuracy of ML wall model with respect to prediction of wall shear stress. In addition, the use of spatial stencil and time delay is also explored within the ML wall modeling framework. It is shown that features associated with gradients of the streamwise and spanwise velocity components have a major impact on the prediction fidelity of wall model, while the effect of gradients of wall-normal velocity component is found to be negligible. Moreover, adding flow feature information from an x-y-z spatial stencil significantly improves the ML model accuracy and generalizability compared to just using local flow features from the matching location. Overall, highest prediction accuracy is achieved when both spatial stencil and time delay features are incorporated within the data-driven wall modeling paradigm.

Aerodynamic and Aeroacoustic Effects of Different Transition Mechanisms on an Airfoil

The paper delves into the detailed sound generation and propagation mechanisms associated with tripping-induced flow perturbations and trailing-edge scattering using wall-resolved large-eddy simulations. Two distinct boundary-layer tripping techniques, namely a geometrically resolved stair strip and an artificially modeled trip using suction and blowing, are investigated. To facilitate comparison, the natural boundary-layer transition is also simulated as a baseline scenario. The analysis takes into account a Reynolds number of 4×105, a Mach number of 0.058, and a nonzero angle of attack of 6.25° over a NACA 0012 airfoil configuration. The mechanisms for sound generation and propagation related to trailing-edge noise remain consistent across the three transition scenarios. However, boundary-layer tripping notably leads to intricate, scenario-specific noise generation: there is an interaction between the laminar separation bubble and tripping for the stair strip case, whereas laminar boundary-layer instability is evident for the suction and blowing scenario. Aerodynamic flowfields involving acoustic noise sources, their propagating natures near the wall, and far-field acoustics are cross-examined in detail. The comprehensive analysis of observed phenomena provides valuable insights for understanding nonlinear flow and acoustic interactions relevant to airfoil noise and designing new types of trips under adverse pressure gradient flows.

Calibrating sub-grid scale models for high-order wall-modeled large eddy simulation

AbstractHigh-order methods have demonstrated orders of magnitude reduction in computational cost for large eddy simulation (LES) over low-order methods in the past decade. Most such simulations are wall-resolved implicit LES (ILES) without an explicit sub-grid scale (SGS) model. The use of high-order ILES for severely under-resolved LES such as wall-modeled LES (WMLES) often runs into robustness and accuracy issues due to the low dissipation embedded in these methods. In the present study, we investigate the performance of several popular SGS models, the static Smagorinsky model, the wall-adapting local eddy-viscosity (WALE) model and the Vreman model, to improve the robustness and accuracy of under-resolved LES using high-order methods. The models are implemented in the high-order unstructured grid LES solver called hpMusic based on the discontinuous flux reconstruction method. The length scales in these SGS models are calibrated using the direct numerical simulation (DNS) database for the turbulent channel flow problem. The Vreman model has been found to produce the most accurate and consistent results with a proper choice of the length scale for WMLES.

Wall-resolved large eddy simulation of mixed-size sand-laden flow

Most of the existing numerical studies on wind-blown sand flow simplify sands into single-size particles, whereas natural wind-blown sand flow is a two-phase flow with mixed-size particles, thus, the simulation of mixed-size sand-laden flow is necessary. In the present work, wall-resolved large eddy simulations of mixed-size sand-laden flows are realized. Each sand in the wind field is tracked using the Lagrangian point-particle model. The transport characteristics of sand particles in mixed-size sand-laden flow are investigated under the premise of considering bed erosion. Considering the significant influence of sand-bed collision on simulation, the splash function is modified in the present simulation according to the previous experimental results. It reveals that in mixed-size sand-laden flow, the fraction of rebound sand particles in all the saltation particles is approximately 0.6, which is twice times of the ejected sand particles, and the modification of the sand rebound angle greatly affects the simulation results of mixed-size sand-laden flow. Meanwhile, the mean size of the saltation sand particles decreases with height and is 20% lower at the top of the saltation layer than that near the sand bed in the present simulation. Further analysis by grouping of sands with their size shows that the sand transport intensity of small sands decreased more rapidly with increasing height. The volume fraction and sand transport intensity of small sand particles exceed those of medium and large sand particles at heights y/δ = 0.05 and y/δ = 0.1.

A wall model for large-eddy simulation of highly compressible flows based on a new scaling of the law of the wall

Wall modelling in large-eddy simulation (LES) is of high importance to allow scale resolving simulations of industrial applications. Numerous models were developed and validated for incompressible flows, including a simple quasi-analytical model based on Reichardt's formula that approximates the law of the wall. In this paper, a scaling is proposed to generalize this wall model to highly compressible flows. First, the results of wall-resolved LES (wrLES) of adiabatic compressible turbulent channel flows at $Re_\tau = 1000$ and at centreline Mach numbers of $M_c= 0.76$ and $1.5$ are presented. Then, three potential scalings of the incompressible wall model are proposed, and their a priori performance is evaluated : (i) the Howarth–Stewartson scaling, (ii) an improved Van Driest scaling and (iii) a new scaling obtained from a blending of those two. The results of wall-modelled LES (wmLES) of compressible channel flows using these three models are compared with the reference wrLES data, showing the superior accuracy of the hybrid scaling. The consistency of the new wall model at low Mach numbers is also verified by comparing the results of a wmLES at $M_c= 0.25$ with those of reference incompressible DNS data at $Re_\tau = 1000$ and $5200$ . Finally, the proposed wall model is also applied to a turbulent channel flow at $M_c=1.5$ and $Re_\tau =5200$ .

Highly resolved LES and URANS computations of a differentially heated square cavity

This investigation focusses on the use of two wall-resolved LES approaches for the simulation of the buoyancy driven flow in a differentially heated square cavity, aiming to advance understanding of the flow phenomena involved. The study also assesses the predictive accuracy and cost-effectiveness of LES for natural convection flows. This assessment is pursued through comparisons of the resulting predictions, and also of their computational requirements, with those for DNS, produced in an earlier study, and those of a number of URANS models produced during the current investigation. URANS models tested here include both high- and low-Reynolds-number schemes. In the case of the high-Reynolds-number models a recently developed numerical variant of the Analytical Wall Function (AWF) has been employed to account for the effects of near-wall turbulence. Both LES approaches result in predictions that are close to those of the DNS but have considerably more modest resource requirements. The URANS approaches, which have the lowest resource requirements, for key parameters such as the local Nusselt number are not as accurate as the LES approaches and their performance is highly dependent on the turbulence model employed. Surprisingly, the simplest of the low-Re models employed, the Launder–Sharma two-equation model, is reasonably close to the DNS data. The approach that requires the least computational resource, the high-Re model with the recent version of the AWF, also offers a promising URANS route. Simulations were also conducted to investigate the flow behaviour at a Rayleigh number that is an order of magnitude higher than that of the available DNS. At such high Rayleigh number, at which the phenomena are of close relevance to cooling loops of nuclear reactors, comparisons are performed between the LES and the high-Re AWF combination. The findings show that as the Rayleigh number increases, the turbulence levels become stronger, but these are confined to more refined structures within the thinner near-wall boundary layers.

TOWARD AN INVARIANT-THRESHOLD TURBULENT-NON-TURBULENT DETECTOR FUNCTION IN VANE CASCADE TRANSITIONAL BOUNDARY LAYERS

This paper investigates the variations on the threshold values used for turbulent-nonturbulent discrimination in transitional boundary layers on the suction surface of a vane under varying levels of inflow turbulence intensities. The invariant-threshold here refers to a potential characteristic of a detector function if a single threshold value can be effectively used to distinguish turbulent and nonturbulent states at any stage of flow development or wall-normal distance. The analysis is based on velocity time series extracted at constant <i>d/δ</i> wall-normal distances obtained by wall-resolved large eddy simulations (LES), with the detector function being the sum of the time derivatives of the spanwise and wall-normal velocity signals, normalized by the boundary layer edge velocity. The detector function's probability density function (PDF) minima is chosen as the turbulent-nonturbulent discrimination threshold. The variations of the selected threshold value are examined at different stages of transition to turbulence and wall-normal distances. Findings suggest that a single threshold value is sufficient for turbulent-nonturbulent discrimination at <i>d/δ</i> > 0.2, throughout the depth of the flow, and for both levels of inflow turbulence intensities considered in this study. Although the appropriate threshold decreases closer to the wall for <i>d/δ</i> < 0.2, it remains invariant across different streamwise positions.

Coexistence of different mechanisms underlying the dynamics of supersonic turbulent flow over a compression ramp

Supersonic turbulent flow over a compression ramp is studied using wall-resolved large eddy simulation with a freestream Mach number of 2.95 and a Reynolds number [based on δ0: the thickness of incoming turbulent boundary layer (TBL)] of 63 560. The unsteady dynamics of the present shock wave/turbulent boundary layer interaction (STBLI) flow are investigated by using dynamic mode decomposition techniques, linear and nonlinear disambiguation optimization, local stability analysis (LSA), and global stability analysis (GSA). By analyzing the dynamic system for the STBLI flow, three dynamically important modes with characteristic spanwise wavelengths of 2δ0, 3δ0, and 6δ0 are captured. The 2δ0 mode approximates the spanwise scale of the Görtler-like vortices and Görtler mode of LSA, suggesting the presence of Görtler instability, which is believed to be related to the unsteady motion of streaks downstream of reattachment in the flow. The features of the 3δ0 mode are also observed in large-scale motions of the incoming TBL, implying the existence of a convective mechanism that is excited and maintained by such motions. Additionally, the GSA results show the most unstable mode features a spanwise wavelength of around 6δ0, indicating the existence of global instability that is believed to be related to the oscillating motion of separation shock. The coexistence of these three mechanisms is confirmed. Discussions on the above findings provide an interpretation for low-frequency unsteadiness that the unsteadiness of surface streaks results from the combined effects of the Görtler instability near flow reattachment and the convection of large-scale motions in the incoming boundary layer, while the low-frequency shock motion may be related to a global mode driven by upstream disturbances.

Shock-wave/turbulent boundary-layer interaction with a flexible panel

The dynamic coupling between a Mach 2.0 shock-wave/turbulent boundary-layer interaction (STBLI) and a flexible panel is investigated. Wall-resolved large-eddy simulations are performed for a baseline interaction over a flat-rigid wall, a coupled interaction with a flexible panel, and a third interaction over a rigid surface that is shaped according to the mean panel deflection of the coupled case. Results show that the flexible panel exhibits self-sustained oscillatory behavior over a broad frequency range, confirming the strong and complex fluid–structure interaction (FSI). The first three bending modes of the panel oscillation are found to contribute most to the unsteady panel response, at frequencies in close agreement with natural frequencies of the mean deformed panel rather than those for the unloaded flat panel. This highlights the importance of the mean panel deformation and the corresponding stiffening in the FSI dynamics. The time-averaged flow shows an enlarged reverse-flow region in the presence of mean surface deformations. The separation-shock unsteadiness is enhanced due to the panel motion, leading to higher wall-pressure fluctuations in the coupled interaction. Spectral analysis of the separation-shock location and bubble-volume signals shows that the STBLI flow strongly couples with the first bending mode of the panel oscillation. This is further confirmed by dynamic mode decomposition of the flow and displacement data, which reveals variations in the reverse-flow region that follow the panel bending motion and appear to drive the separation-shock unsteadiness. Low-frequency modes that are not associated with the fluid–structure coupling, in turn, are qualitatively similar to those obtained for the rigid-wall interactions, indicating that the characteristic low-frequency unsteadiness of STBLI coexists with the dynamics emerging from the fluid–structure coupling. Based on the present results, unsteady FSIs involving STBLIs and flexible panels are likely to accentuate rather than mitigate the undesirable features of STBLIs.

Passive flow control for the E423 airfoil utilizing an equilateral triangular trip

A passive flow control technique of utilizing an equilateral triangular trip close to the leading edge was developed and tested for a micro-scale Eppler E423 airfoil at a Reynolds number based on the chord of 40 000. The analysis was carried out via high-order wall-resolved large eddy simulation using the computational solver HpMusic. Angles of attack of 5° and 20° were tested. It was shown that at an angle of attack of 5°, the trip height of two times the local boundary layer thickness outperformed existing passive flow control techniques from the literature by almost a factor of five in terms of the lift-to-drag ratio. To understand the underlying physics which allowed the trip to provide this very significant performance benefit, metrics such as pressure coefficient profiles, oil flows, iso-surfaces of Q-criteria, and leading-edge flow behavior were examined. It was found that this trip configuration simultaneously removed the flow separation regions on both the suction and pressure sides of the wing.

Assessment of the DDES-γ Model for the Simulation of a Highly Loaded Turbine Cascade

Abstract The flow over the linear low-pressure turbine cascade MTU-T161 at Re = 90,000 is analyzed using delayed detached eddy simulations (DDES). At this operating point, the low Reynolds number and the high loading of the blade result in a separation bubble and a separation-induced transition of the flow over the suction side. The utilized DDES method is based on a vorticity-based formulation to calculate the subgrid length scales, and it incorporates the one-equation γ-transition model. The computational model of the MTU-T161 cascade consists of one blade passage, including the diverging viscous sidewalls. To reproduce realistic operating conditions and to mimic the experiments, synthetic turbulence is prescribed at the inlet of the computational domain. Several studies are performed to assess the accuracy and performance of the DDES one-equation γ-transition model against experimental data and a benchmark large eddy simulations (LES). The primary focus is on the prediction of the separation and the separation-induced transition mechanism. First of all, a systematic grid convergence study is conducted and grid criteria are derived in order to ensure a satisfactory agreement of the flow metrics, such as isentropic Mach number, friction coefficient distribution, and total pressure wake losses at mid-span with experimental data. Furthermore, a detailed analysis of the DDES model parameters, such as shielding function and subgrid length scale, is presented and the effect of these parameters on the prediction accuracy of the separation bubble region is analyzed. The analysis of the suction side boundary layer indicates that the turbulent kinetic energy should be resolved and modeled properly in order to represent the separation bubble correctly. In particular, the correct prediction of the separated shear layer above the separation bubble is of utmost importance. The results of the simulations reveal higher demands on grid resolution for such transitional flows than typically have been reported in the literature for turbulent boundary layers. This higher demand on grid resolution results in more expensive simulations than Reynolds-averaged Navier–Stokes (RANS). Nevertheless, DDES requires less computing time than wall-resolved LES. Additionally, the results of the transitional DDES model are compared to DDES without a transition model, an RANS eddy viscosity model, and a reference LES. The results show that the DDES approach needs to be coupled with a transition model, such as the one-equation γ-transition model, in order to capture the flow topology over a highly loaded turbine blade correctly. The benefit of the DDES one-equation γ-transition model becomes particularly evident when predicting the separated shear layer, the transition process, and the subsequent reattachment. The RANS eddy viscosity turbulence and transition models applied within our study are not able to predict the aforementioned mechanisms accurately. For highly loaded turbine blades in particular, the accurate prediction of flow separation and potential reattachment is crucial for the aerodynamic design of turbines, since large parts of the total pressure loss are generated in the separated region. For this reason, the DDES one-equation γ-transition model can be a good compromise in terms of predictive accuracy and computational costs.

Wake bi-modality: the effect of upstream boundary layer dynamics

The turbulent wake past a square-back Ahmed body in close proximity to the ground experiences random side-to-side switching between two asymmetric positions, a phenomenon known as bi-modality. It has been observed to be sensitive to the dynamics of the upstream boundary layers formed along the body surfaces. Close to the body fore end, these separate and reattach, with hairpin vortices emanating from the reattachment points and growing along the surfaces before breaking down upstream of the base. This study uses wall-resolved large eddy simulations to investigate the effect of using suction to suppress these upstream boundary layer separations on the wake bi-modality. It is seen that, in the unforced flow (in the absence of suction), the smaller top and side surface vortices resulting from breakdown interact as they convect downstream. Steady suction is confirmed to suppress the boundary layer separations on the different body surfaces. When the boundary layer separations on the two side (vertical) surfaces are suppressed, it is found that horizontal bi-modality is completely inhibited with weak vertical asymmetry preserved. The interaction of the small top/side surface vortices is interrupted, damping boundary layer fluctuations just upstream of the base. Applying suction on different combinations of side/top/bottom boundary layer separations is found to have different effects on the underbody flow and the wake vertical balance, with bi-modality suppression dependent on side surface suction. This confirms that bi-modality is triggered, at least in part, by boundary layer disturbances on the surfaces perpendicular to the switching direction.

Assessment of Reynolds number effects in supersonic turbulent boundary layers

Wall-resolved large-eddy simulations (LES) are performed to investigate Reynolds number effects in supersonic turbulent boundary layers (TBLs) at Mach 2.0. The resulting database covers more than a decade of friction Reynolds number Reτ, from 242 to 5554, which considerably extends the parameter range of current high-fidelity numerical studies. Reynolds number trends are identified on a variety of statistics for skin-friction, velocity and thermodynamic variables. The efficacy of recent scaling laws as well as compressibility effects are also assessed. In particular, we observe the breakdown of Morkovin’s hypothesis for third-order velocity statistics, in agreement with previous observations for variable-property flows at low Mach number. Special attention is also placed on the size and topology of the turbulent structures populating the TBL, with an emphasis on the outer-layer motions at high Reynolds number. The corresponding streamwise spectra of streamwise velocity fluctuations show a clear separation between inner and outer scales, where energetic peaks are found at streamwise wavelengths of λx+≈700 and λx/δ0≈6. The spanwise spacing of the outer-layer structures, in turn, is found to be insensitive to the Reynolds number and equal to ∼0.7δ0. It is also found that the integral length-scales in spanwise direction for the temperature, streamwise and spanwise velocity fields appear to progressively collapse with increasing Reynolds number. The modulating influence that the outer-layer structures exert on the near-wall turbulence is also clearly visible in many of the metrics discussed. In addition, the present LES data is further exploited to assess the Reτ-sensitivity of uniform momentum regions in the flow. We find that the resulting probability density function of the number of zones as well as its evolution with Reτ agrees well with incompressible data. This suggests that uniform zones, which have been associated with outer-layer dynamics, are not strongly influenced by compressibility at the considered Mach number.

Analysis of improved digital filter inflow generation methods for compressible turbulent boundary layers

We propose several enhancements to improve the accuracy and performance of the digital filter turbulent inflow generation technique and assess their efficacy in the context of wall-resolved large-eddy simulations of a compressible turbulent boundary layer. Improvements of accuracy include a more realistic correlation function for the transversal directions, target length scales that vary with wall-distance, and a counter-intuitive approach that involves the suppression of streamwise velocity fluctuations at the inflow. For improving the computational performance, we propose to generate the inflow data in parallel in single precision and at a prescribed time interval based on the turbulence time scale, and not at every time-step of the simulation. Based on the results of 7 wall-resolved large-eddy simulations, we find that the new correlation functions and the considered performance improvements are beneficial and therefore desired. Suppressing streamwise velocity fluctuations at the inflow leads to the fastest relaxation of the pressure fluctuations; however, this approach increases the adaptation length defined in terms of compliance with the von Kármán integral equation. The adaptation length can be shortened by artificially increasing the wall-normal Reynolds stresses, thereby preserving the desired turbulence kinetic energy level. A detailed inspection of the Reynolds stress transport budgets reveals that the observed spurious spatial transients are largely driven by pressure-related terms. For instance, increased values of u′p′¯ are found throughout the computational domain when a physical Reynolds stress distribution is prescribed at the inflow. Therefore, efforts to enhance digital filter techniques should aim at modeling pressure fluctuations as well as their correlation with the velocity components.

Large-eddy simulation of non-isothermal flow and heat transfer in an enclosed rotor–stator cavity

This paper presents a numerical study of non-isothermal flow and heat transfer in an enclosed rotor–stator disk cavity. Wall-resolved large-eddy simulation is implemented with attention focused on the Reynolds number effect (Reϕ=105 and 5×105) and non-isothermal effects, including different thermal Rossby numbers (βΔT=0.05–0.25) and heat convection types (axial heat convection, axial and radial heat convection). Verification and validation are implemented by comparing with available experimental data and examining the resolution of the present LES. Although the Batchelor flow type is unchanged, the non-isothermal effect is found to be responsible for the acceleration of flow in both disk boundary layers and cavity core. Temperature in the cavity is negatively correlated with the Rayleigh number. The radial inward temperature gradient is found to be responsible for destabilizing the flow near the cylinders, due to the effect of centrifugal buoyant force. Transition from laminar to turbulent is shown in the rotor boundary layer for Reϕ=5×105 conditions. The non-isothermal effect increases flow unsteadiness in both rotor and stator boundary layers, but it has little to no effect on the transition position, as well as the kinematic and thermal boundary layer thicknesses, which are mainly controlled by Reϕ. Spectral proper orthogonal decomposition is implemented to explore the Reynolds number effect. Analyzing the most energetic mode and its energy spectra provides some insights into the understanding of non-isothermal rotor–stator disk cavity flows.

Capturing the bi-stable asymmetry exchanges in the turbulent air wake of a simplified frigate using Large Eddy Simulation

Turbulent air wake of a Simplified Frigate Ship (SFS2) model is simulated using wall-resolved Large Eddy Simulation. The wake behind the superstructure presents an asymmetric configuration for an otherwise symmetric free stream condition. Bi-stable characteristics of the wake have also been simulated as evidenced by asymmetric state switching events. The wake behind the stern also exhibits bi-stable switching and adapts an opposite asymmetric orientation to the wake behind the superstructure, as also observed experimentally. The resulting time-averaged flow fields for each bi-stable state reveal a differential flux formation around the upstream sides of the superstructure. The instantaneous wake fields also suggest a possible influence of flow structures developed around the upstream end of the superstructure to generate a particular asymmetric wake state downstream.