Particle Simulation of Swirling Flows

The vortex particle method has been applied to the axisymmetric swirling flow of a viscous fluid. The formulation used yields two transport equations which have been solved within the lagrangian framework of particle method. The diffusion operator for both equations has been approximated by means of a Particle Strength Exchange scheme. Applications to the cases of one isolated vortex ring and two co-rotating vortex rings illustrate the interest of this new method. Special attention has been devoted to the vorticity production resulting from the interaction between the azimuthal components of vorticity and velocity. The generation of small eddies at the boundary of the vortex ring cross-section has been particularly investigated. This revised version was published online in July 2006 with corrections to the Cover Date.     

Analysis of High-order Explicit LES Dynamic Modeling Applied to Airfoil Flows

A high-order low dissipative numerical framework is discussed to tackle simultaneously the modeling of unresolved sub-grid scale flow turbulence and the capturing of shock waves. The flows around two different airfoil profiles are simulated using a Spectral Difference discretisation scheme. First, a transitional, almost incompressible, low Reynolds number flow over a Selig-Donovan 7003 airfoil. Second, a high Reynolds number flow over a RAE2822 airfoil under transonic conditions. These flows feature both laminar and turbulent flow physics and are thus particularly challenging for turbulence sub-grid scale modeling. The accuracy of the recently developed Spectral Element Dynamic Model, specifically capable of detecting spatial under-resolution in high-order flow simulations, is evaluated. Concerning the test in transonic conditions, the additional complexity due to the presence of shock waves has been handled using an artificial viscosity shock-capturing technique based on bulk viscosity. To mitigate the impact of the shock-capturing on turbulence dissipation, it was necessary to combine the high-order modal-type shock detection with a usual sensor measuring the local flow compressibility.

     

Fluid Dynamical Analysis of Atmospheric Reacting and Isothermal Swirling Flows

A series of unconfined swirling premixed natural gas/air flames was investigated. Reynolds-numbers spanned from 10,000 to 42,300. Respective isothermal flows were studied additionally to gain insight into changes of fluid dynamical features caused by combustion. Statistical moments, Reynolds-stresses, temporal time scales, spatial length scales, and power spectral densities were deduced from one- and two-point laser Doppler velocimetry (LDV) data. Properties of the turbulent flows and dependencies on Reynolds-number, swirl number, and chemical reactions are discussed. Most distinct differences between combusting and isothermal flows were precessing vortex cores (PVC) occurring only for the latter cases. The study is aimed to serve as a database of a generic flame geometry featuring important characteristics of industrial applications for validation of numerical simulations. Therefore, nozzle exit profiles as important inlet conditions to numerical simulations are thoroughly documented.     

Confined and Agitated Swirling Flows with Applications in Chemical Engineering

Turbulent, swirling flows are encountered frequently in chemical engineering practice. In this article, experiments and simulations on two classes of swirling flows, viz. agitated flows (stirred tanks), and confined swirling flows are discussed. Results of large-eddy simulations of stirred tank flow are compared with experimental data, mainly phase-resolved LDA data of the flow in the vicinity of the impeller. Next to the average velocity field, also the turbulent kinetic energy, and the anisotropy of the Reynolds stress tensor have been assessed. An important application of confined swirling flow is the cyclone separator (hydrocyclones for the separation of liquids, gas cyclones for gas-solid separation). The flow in a swirl tube geometry exhibiting many of the typical features of swirl flows (e.g. vortex breakdown) is discussed. Furthermore, a large-eddy simulation of the gas flow in a high-efficiency Stairmand cyclone separator is presented. Two examples of process modeling based on flow simulations are briefly treated: orthokinetic agglomeration of crystals in a stirred tank, and particle separation in a cyclone. This revised version was published online in July 2006 with corrections to the Cover Date.     

Large Eddy Simulation of Turbulent Confined Highly Swirling Annular Flows

The present paper discusses the Large Eddy Simulation of a confined non-reacting annular swirling jet. The configuration corresponds to a well investigated flow studied experimentally by Sheen (1993). The flow field is characterised by a high swirl number resulting in relatively complex features. The challenging behaviour of the flow is governed by the interaction of several recirculation zones. The central recirculation zone formed by the swirling jet is strongly affected by the cylindrical centre body which acts as a bluff body. The flow features coherent structures such as Precessing Vortex Cores (PVCs), which create regions with high velocity fluctuations. The simulations presented comprise a detailed investigation of the parameters controlling the inert flow and a thorough comparison with the experimental data.     

On the Use of Embedded Meshes in the LES of External Flows

In the present study, the possibility of performing cheap and accurate LES computations using unstructured grids on complex geometries has been examined. The test case considered is a turbulent flow around a thin controlled diffusion airfoil at low speed, with angle of attack of 8° to the incoming free-stream flow. The flow velocity is set to satisfy the chord Reynolds number of 1.2 × 10⁵. Two subgrid models are tested; the classical Smagorinsky model and the WALE model. The predictions are found to be in good agreement with the experimental data both in terms of flow dynamics and pressure spectra. The results are found to be independent of the LES physical models tested herein, although this conclusion is issued with the caveat that a careful grid generation procedure using (2:3 and 3:4 type) refinements was required to reduce the jump size across the different grid resolutions.     

The Decay of Swirling Flows in a Type of Cross-Section-Varying Pipes

ＩｎｔｒｏｄｕｃｔｉｏｎＳｗｉｒｌｉｎｇｆｌｏｗ (ｏｒｓｐｉｒａｌｆｌｏｗ )ｉｓｗｉｄｅｌｙｕｓｅｄｉｎｉｎｄｕｓｔｒｙ ,ｅ.ｇ .ｉｎｊｅｔｔｅｃｈｎｏｌｏｇｙ ;ｃｏｍｂｕｓｔｉｏｎｔｅｃｈｎｏｌｏｇｙ ;ｃｙｃｌｏｎｅ/ｈｙｄｒｏｃｙｃｌｏｎｅｓｅｐａｒａｔｏｒｓ;ｐｎｅｕｍａｔｉｃｔｒａｎｓｐｏｒｔａｔｉｏｎｅｔｃ .Ｓｗｉｒｌｉｓｃｏｍｍｏｎｌｙｉｍｐａｒｔｅｄｔｏｆｌｏｗｓｅｎｔｅｒｉｎｇｆｕｒｎａｃｅｓａｎｄｃｏｍｂｕｓｔｏｒｓｉｎｏｒｄｅｒｔｏｅｎｈａｎｃｅｆｌａｍｅｓｔａｂｉｌｉｔｙａｎ…     

Investigation of a Dynamic Hybrid RANS/LES Modelling Methodology for Finite-Volume CFD Simulations

This paper investigates a recently proposed dynamic hybrid RANS-LES framework using a general-purpose finite-volume flow solver. The new method is highly generalized, allowing coupling of any selected RANS model with any selected LES model and containing no explicit grid dependence in its formulation. Selected results are presented for three test cases: two-dimensional channel flow, backward facing step, and a nozzle flow relevant to biomedical applications. Comparison with experimental and DNS data, and with other hybrid RANS-LES approaches, highlights the advantages of the new method and suggests that further investigation is warranted.     

From Two-Point Closures of Isotropic Turbulence to LES of Shear Flows

We first recall the EDQNM two-point closure approach of three-dimensional isotropic turbulence. It allows in particular prediction of the infrared kinetic-energy dynamics (with ak ₄ backscatter) and the associated time-decay law of kinetic-energy, useful in particular for one-point closure modelling. Afterwards, we show how the spectral eddy viscosity concept may be used for large-eddy simulations: we introduce the plateau-peak model and the spectral-dynamic models. They are applied to decaying isotropic turbulence, and allow recovery of the EDQNM infrared energy dynamics. Anew infrared k ₂ law for the pressure spectrum, predicted by the closure, is also well verified. Assuming that subgrid scales are not too far from isotropy, the spectral-dynamic model is applied to the channel flow at h ₊= 390, with statistics in very good agreement with DNS, while reducing considerably the computational time. We study with the aid of DNS and LES the case of the channel rotating about an axis of spanwise direction. The calculations allow to recover the universal linear behaviour of the mean velocity profile, with a local Rossby number equal to −1. We present also LES (using the Grenoble Filtered Structure-Function Model), of a turbulent boundary layer passing over a cavity. Finally, we make some remarks on the future of LES for industrial applications. This revised version was published online in July 2006 with corrections to the Cover Date.     

Prediction of Global Extinction Conditions and Dynamics in Swirling Non-premixed Flames Using LES/CMC Modelling

The Large Eddy Simulation (LES)/three-dimensional Conditional Moment Closure (CMC) model with detailed chemistry is applied to predict the operating condition and dynamics of complete extinction (blow-off) in swirling non-premixed methane flames. Using model constants previously selected to provide relatively accurate predictions of the degree of local extinction in the piloted jet flames Sandia D ?F, the error in the blow-off air velocity predicted by LES/3D-CMC in short, recirculating flames with strong swirl for a range of fuel flow rates is within 25 % of the experimental value, which is considered a new and promising result for combustion LES that has not been applied before for the prediction of the whole blow-off curve in complex geometries. The results also show that during the blow-off transient, the total heat release gradually decreases over a duration that agrees well with experiment. The evolution of localized extinction, reactive scalars and scalar dissipation rate is analyzed. It has been observed that a consistent symptom for flames approaching blow-off is the appearance of high-frequency and high-magnitude fluctuations of the conditionally filtered stoichiometric scalar dissipation rate, resulting in an increased fraction of local extinction over the stoichiometric mixture fraction iso-surfaces. It is also shown that the blow-off time changes with the different blow-off conditions.     

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.     

基于正交设计的冲击荷载下型钢混凝土梁动力性能有限元分析

为研究冲击荷载下型钢混凝土梁的动力性能,取试验结果验证模型有效性后,使用有限元软件ABAQUS建立了型钢混凝土梁落锤冲击模型．在对比了型钢混凝土梁与钢筋混凝土梁的动力性能后,分析了落锤速度对型钢混凝土梁的冲击位移、冲击力、冲击力位移曲线、惯性效应、塑性耗能等性能的影响,然后基于正交设计,研究了落锤质量、落锤速度、混凝土强度等参数对型钢混凝土梁动力性能的影响．研究结果表明：型钢混凝土梁的抗冲击性能较好;落锤速度是影响构件动力性能的主要参数,落锤质量对动力性能的影响较大,提高混凝土强度有助于增强构件的抗冲击性能．     

A Mesh Adaptation Strategy to Predict Pressure Losses in LES of Swirled Flows

Large-Eddy Simulation (LES) has become a potent tool to investigate instabilities in swirl flows even for complex, industrial geometries. However, the accurate prediction of pressure losses on these complex flows remains difficult. The paper identifies localised near-wall resolution issues as an important factor to improve accuracy and proposes a solution with an adaptive mesh h-refinement strategy relying on the tetrahedral fully automatic MMG3D library of Dapogny et al. (J. Comput. Phys. 262, 358-378, 2014) using a novel sensor based on the dissipation of kinetic energy. Using a joint experimental and numerical LES study, the methodology is first validated on a simple diaphragm flow before to be applied on a swirler with two counter-rotating passages. The results demonstrate that the new sensor and adaptation approach can effectively produce the desired local mesh refinement to match the target losses, measured experimentally. Results shows that the accuracy of pressure losses prediction is mainly controlled by the mesh quality and density in the swirler passages. The refinement also improves the computed velocity and turbulence profiles at the swirler outlet, compared to PIV results. The significant improvement of results confirms that the sensor is able to identify the relevant physics of turbulent flows that is essential for the overall accuracy of LES. Finally, in the appendix, an additional comparison of the sensor fields on tetrahedral and hexahedral meshes demonstrates that the methodology is broadly applicable to all mesh types.     

Numerical Investigation of Unsteady Transitional Flows in Turbomachinery Components Based on a RANS Approach

This paper presents results of the numerical simulation of periodically unsteady flows with focus on turbomachinery applications. The unsteady CFD solver used for the simulations is based on the Reynolds averaged Navier–Stokes equations. The numerical scheme applies an extended version of the Spalart–Allmaras one-equation turbulence model coupled with a transition correlation. The first example of validation consists of boundary layer flow with separation bubble on a flat plate, both under steady and periodically unsteady main flow conditions. The investigation includes a variation of the major parameters Strouhal number, amplitude, and Reynolds number. The second, more complex test case consists of the flow through a cascade of turbine blades which is influenced by wakes periodically passing over the cascade. The computations were carried out for two different blade loadings. The results of the numerical simulations are discussed and compared with experimental data in detail. Special emphasis is given to the investigation of boundary layers with regard to transition, separation and reattachment under the influence of main flow unsteadiness. This revised version was published online in July 2006 with corrections to the Cover Date.     

The State of the Art of Hybrid RANS/LES Modeling for the Simulation of Turbulent Flows

This review presents the state of the art of hybrid RANS/LES modeling for the simulation of turbulent flows. After recalling the modeling used in RANS and LES methodologies, we propose in a first step a theoretical formalism developed in the spectral space that allows to unify the RANS and LES methods from a physical standpoint. In a second step, we discuss the principle of the hybrid RANS/LES methods capable of representing a RANS-type behavior in the vicinity of a solid boundary and an LES-type behavior far away from the wall boundary. Then, we analyze the principal hybrid RANS/LES methods usually used to perform numerical simulation of turbulent flows encountered in engineering applications. In particular, we investigate the very large eddy simulation (VLES), the detached eddy simulation (DES), the partially integrated transport modeling (PITM) method, the partially averaged Navier-Stokes (PANS) method, and the scale adaptive simulation (SAS) from a physical point of view. Finally, we establish the connection between these methods and more precisely, the link between PITM and PANS as well as DES and PITM showing that these methods that have been built by different ways, practical or theoretical manners have common points of comparison. It is the opinion of the author to consider that the most appropriate method for a particular application will depend on the expectations of the engineer and the computational resources the user is prepared to expend on the problem.     

Numerical Investigation of Off-Design Regimes of Flow Past Power-Law Waveriders Based on the Flows Behind Plane Shocks

The results of a numerical investigation of supersonic off-design flow past waveriders at the freestream Mach numbers M = 4 and 8 are presented. Flow regimes with M both greater and smaller than the design value M _d are analyzed. Configurations based on the flows behind plane shocks and described by power-law functions are considered. The results are obtained by the finite-volume solution of the Euler equations using higher-order TVD Runge-Kutta schemes.     

Performances of LES and RANS Models for Simulation of Complex Flows in a Coaxial Jet Mixer

Flow and mixing processes in a classical coaxial jet mixer have been investigated numerically. Calculations have been performed using three Large Eddy Simulation models and three unsteady RANS models. The time averaged mixture fraction and axial velocity, their rms values and energy spectra are compared with LIF and LDA measurements for both j- and r-modes of the jet mixer flow. A special attention is paid to the ability of different models to reproduce unsteady effects. The analysis demonstrates the superiority of the LES method with the dynamic mixed SGS model (DMM) with respect to other RANS and LES models.     

LES Investigation of the Flow Through an Effusion-Cooled Aeronautical Combustor Model

The present study is devoted to the analysis of the behaviour of the flow through an effusion-cooled aeronautical combustor model. High-fidelity calculations are performed on an experimental model of a combustion chamber multi-perforated wall and compared to experimental measurements. The effect of combustion instability on the effusion-cooling system is investigated by studying the interaction of an acoustic wave with the jets-in-crossflow issued from the cooling plate. It is shown that the mass-flow rate through the plate can be drastically reduced by the acoustic wave, which demonstrates the destructive effect that such instability may have on the cooling of an aeronautical combustion chamber.     

Numerical Investigation of Multi-Jet Channel Flows

The results of a numerical simulation of the three-dimensional outflow of a system of circular supersonic turbulent jets into a cocurrent supersonic (or subsonic) air flow in a partially bounded region are given. Solutions are obtained by the splitting method using a matrix sweep of the parabolized Navier-Stokes equations. Assuming that the flow is nonseparated in the boundary layer, features of the three-dimensional structure of the jet system are investigated as functions of the pressure ratio number and the jet and cocurrent flow Mach numbers.