LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.We acknowledge “Red Española de Surpercomputación” (RES) for awarding us access to the MareNostrum III machine based in Barcelona, Spain (Ref. FI-2015-2-0026 and FI-2015-3-0011). We also acknowledge PRACE for awarding us access to Fermi and Marconi Supercomputers at Cineca, Italy (Ref. 2015133120). Oriol Lehmkuhl acknowledges a PDJ 2014 Grant by AGAUR (Generalitat de Catalunya). Ugo Piomelli acknowledges the support of the Natural Sciences and Engineering Research Council (NSERC) of Canada under the Discovery Grant Programme (Grant No. RGPIN-2016-04391). Ricard Borrell acknowledges a Juan de la Cierva postdoctoral grant (IJCI-2014-21034). Ivette Rodriguez, Oriol Lehmkuhl, Ricard Borrell and Assensi Oliva acknowledge Ministerio de Economía y Competitividad, Secretaría de Estado de Investigación, Desarrollo e Innovación, Spain (ref. ENE2014-60577-R).Peer ReviewedPostprint (author's final draft

Using LES to Study Reacting Flows and Instabilities in Annular Combustion Chambers

Great prominence is put on the design of aeronautical gas turbines due to increasingly stringent regulations and the need to tackle rising fuel prices. This drive towards innovation has resulted sometimes in new concepts being prone to combustion instabilities. In the particular field of annular combustion chambers, these instabilities often take the form of azimuthal modes. To predict these modes, one must compute the full combustion chamber, which remained out of reach until very recently and the development of massively parallel computers. Since one of the most limiting factors in performing Large Eddy Simulation (LES) of real combustors is estimating the adequate grid, the effects of mesh resolution are investigated by computing full annular LES of a realistic helicopter combustion chamber on three grids, respectively made of 38, 93 and 336 million elements. Results are compared in terms of mean and fluctuating fields. LES captures self-established azimuthal modes. The presence and structure of the modes is discussed. This study therefore highlights the potential of LES for studying combustion instabilities in annular gas turbine combustors

Heat transfer of a buoyancy-aided turbulent flow in a trapezoidal annulus

The objective of this paper is to report a numerical investigation into the heat transfer of a buoyancy-aided flow in a rod-bundle-like channel. The flow field is simulated using large eddy simulation (LES) with WALE SGS model and the buoyant force is taken into account using the Boussinesq approximation. The general trend of the effect of buoyancy on the overall heat transfer is similar to that in a pipe flow, but the effect on the regional heat transfer varies greatly. This has resulted from a number of interplaying factors, including, the redistribution of the mass flow in the various sub-channels, the non-uniform buoyancy effects on turbulence in different regions of the domain and the behaviour of the large flow structures in the flow channel. These factors together make the effect of buoyancy on heat transfer in the considered flow channel really complicated, while the last factor has been found to have the most pronounced effect in most cases studied

SPLEND1D, a reduced one-dimensional model to investigate the physics of plasma detachment

Studying the process of divertor detachment and the associated complex interplay of plasma dynamics and atomic physics processes is of utmost importance for future fusion reactors. Whilst simplified analytical models exist to interpret the general features of detachment, they are limited in their predictive power, and complex 2D or even 3D codes are generally required to provide a self-consistent picture of the divertor. As an intermediate step, 1D models of the Scrape-Off Layer (SOL) can be particularly insightful as the dynamics are greatly simplified, while still self-consistently including various source and sink terms at play, as well as additional important effects such as flows. These codes can be used to shed light on the physics at play, to perform fast parameter scans, or to interpret experiments. In this paper, we introduce the SPLEND1D (Simulator of PLasma ENabling Detachment in 1D) code: a fast and versatile 1D SOL model. We present in detail the model that is implemented in SPLEND1D. We then employ the code to explore various elements of detachment physics for parameters typical of the Tokamak \`a Configuration Variable (TCV), including the atomic physics and other processes behind power and momentum losses, and explore the various hypotheses and free parameters of the model

Numerical and analytical investigation of the indirect combustion noise in a nozzle

International audienceAnalytical and numerical assessments of the indirect noise generated through a nozzle are presented. The configuration corresponds to an experimental setup operated at DLR by Bake et al. (2008) where an entropy wave is generated upstream of the nozzle by means of an electrical heating device. Both 3-D and 2-D axisymmetric simulations are performed to demonstrate that the experiment is mostly driven by linear acoustic phenomena, including pressure wave reflection at the outlet and entropy-to-acoustic conversion in the accelerated regions. Results show that the acoustic impedance downstream of the nozzle must be accounted for appropriately in order to recover the experimental pressure signal. A good agreement is also obtained with a purely analytical assessment based on the Marble and Candel compact nozzle approximation

Crashworthiness of foam-filled and reinforced honeycomb crash absorbers in transverse direction

Honeycomb crash absorbers have been widely studied as energy absorption devices for use in automotive industries. However, none of these investigations have studied the side impact of empty and foam-filled honeycomb absorbers and adding stiffeners between the different layers of the corrugated sheets which are composing the honeycomb structure to analyse the structure under transverse (L-direction) impacts. In this paper, the foam-filled and reinforced honeycomb crash absorbers are investigated under axial (T) and transverse (L) loading directions. Experimental results for both empty and foam-filled specimens under quasi-static and impact loads were implemented to validate the developed finite element model. Finite element analysis (FEA) was performed to find out the crashworthiness behaviour of the structure under axial and transverse impacts according to road conditions. Finally, a new design of stiffened honeycomb crash absorber was developed and investigated to reduce the level of acceleration experienced by the passengers during the crash event. In this regard, it is concluded that all the requirements related to the energy absorption capabilities and generated deceleration under impact loading can be met by introducing an advanced method to reinforce honeycomb absorbers using stiffeners. It is also proven that the thickness of these stiffeners will not significantly influence the force levels. Due to increase of wall thickness from 1 to 3 mm, the mean crushing force increased from 129 kN to 148 kN. This growth is not sufficient as the goal is to obtain a mean crushing force of 300 kN. Thickening the stiffeners would lead to a loss of efficiency of the structure, as the small increase in mean force would not make up for the gain in mass. Thus, increasing the corrugated sheet’ thickness becomes necessary.Thailand Metal and Materials Technology Center (MTEC); Bangkok Expressway; Metro Public Company Ltd; Royal Academy of Engineering (RAE) through the Engineering X Transforming Systems through Partnership programm

Constructing Physically Consistent Subgrid-Scale Models for Large-Eddy Simulation of Incompressible Turbulent Flows

We studied the construction of subgrid-scale models for large-eddy simulationof incompressible turbulent flows, focusing on consistency with importantmathematical and physical properties. In particular, we considered the symmetriesof the Navier-Stokes equations, and the near-wall scaling and dissipation behaviorof the turbulent stresses. After showing that existing models do not all satisfy thedesired properties, we discussed a general class of subgrid-scale models based onthe local filtered velocity gradient. We provided examples of models from this classthat preserve several of the symmetries of the Navier-Stokes equations and exhibitthe same near-wall scaling behavior as the turbulent stresses. Furthermore, thesemodels are capable of describing nondissipative effects

Numerical analysis of the linear and nonlinear vortex-sound interaction in a T-junction

T-junctions correspond to a classical academic configuration employed to unravel the vortexsound interaction leading to self-sustained oscillations. It is composed of a closed deep cavity exposed to a low-Mach grazing flow, in which an unstable shear layer can develop. Many studies usually consider this hydrodynamic instability either as a “flapping shear layer", or as a “discrete vortex shedding", which then couples with the acoustic field. This paper follows the idea that these two descriptions are related to the linear and non-linear regimes of the shear layer response to acoustic waves, and thus intends to further analyze these regimes and their transition. To do so, a typical T-junction turbulent flow is computed by forced Large Eddy Simulation (LES) where acoustic waves are injected at several amplitudes. The flow response is extracted, and exhibits a linear regime as well as two distinct non-linear regimes where a partial saturation of the response occurs. The post-processing of the flow field in the three situations reveals that a flapping mechanism exists at low wave amplitudes, whereas a vortex shedding appears for highest acoustic levels. For moderate wave amplitudes, the behavior of the shear layer lies between these two classical views. This suggests that during self-sustained oscillations, a transition between these scenarios occurs, starting from a flapping motion followed by a vortex shedding

Reduction in benefits of total flux expansion on divertor detachment due to parallel flows

The Super-X divertor (SXD) is an alternative divertor configuration leveraging total flux expansion at the outer strike point (OSP). Key features for the attractiveness of the SXD are facilitated detachment access and control, as predicted by the extended 2-point model (2PM). However, parallel flows are not consistently included in the 2PM. In this work, the 2PM is refined to overcome this limitation: the role of total flux expansion on the pressure balance is made explicit, by including the effect of parallel flows. In consequence, the effect of total flux expansion on detachment access and control is weakened, compared to predictions of the 2PM. This new model partially explains discrepancies between the 2PM and experiments performed on TCV, in ohmic L-mode scenarios, where in core density ramps in lower single-null (SN) configuration, the impact of the OSP major radius Rt on the CIII emission front movement in the divertor outer leg - used as a proxy for the plasma temperature - is substantially weaker than 2PM predictions; and in OSP sweeps in lower and upper SN configurations, with a constant core density, the peak parallel particle flux density at the OSP is almost independent of Rt, while the 2PM predicts a linear dependence. Finally, analytical and numerical modelling of parallel flows in the divertor is presented, to support the argument. It is shown that an increase in total flux expansion can favour supersonic flows at the OSP. Parallel flows are also shown to be relevant by analysing SOLPS-ITER simulations of TCV