Spectral and modal analysis of the flow in a helical coil steam generator experiment with large eddy simulation

This paper presents an analysis of flow in the shell side of a 24-degree cut helical coil steam generator (HCSG) with large eddy simulation (LES). In the study, the open source, high-order spectral element CFD code Nek5000 was employed for numerical simulation. LES was performed at a Reynolds number of 9000 based on the inlet velocity and tube diameter. Spectral analysis was carried out along the tube bundle to investigate vortex shedding and to reveal shedding frequencies in HCSGs. Wavelet analysis was carried out to determine the time-frequency information of the fluctuating velocities in the domain. Flow was decomposed into modes with proper orthogonal decomposition, which revealed coherent structures in the flow and evolution of coherent structures in time. For the validation of Nek5000 in complex geometries, results from the numerical simulation were compared with experiments conducted at Texas A&M University. Overall, the LES results showed good agreement with experimental data.     

Numerical Analysis of the Impact of the Interior Nozzle Geometry on Low Mach Number Jet Acoustics

The flow and acoustic fields of subsonic turbulent hot jets exhausting from three divergent nozzles at a Mach number M=0.12 based on the nozzle exit velocity are conducted using a hybrid CFD-CAA method. The flow field is computed by highly resolved large-eddy simulations (LES) and the acoustic field is computed by solving the acoustic perturbation equations (APE) whose acoustic source terms are determined by the LES. The LES of the computational domain includes the interior of the nozzle geometry. Synthetic turbulence is prescribed at the inlet of the nozzle to mimic the exit conditions downstream of the last turbine stage. The LES is based on hierarchically refined Cartesian meshes, where the nozzle wall boundaries are resolved by a conservative cut-cell method. The APE solution is determined on a block structured mesh. Three nozzle geometries of increasing complexity are considered, i.e., the flow and acoustic fields of a clean geometry without any built-in components, a nozzle with a centerbody, and a nozzle with a centerbody plus struts are computed. Spectral distributions of the LES based turbulent fluctuated quantities inside the nozzle and further downstream are analyzed in detail. The noise sources in the near field are noticeably influenced by the nozzle built-in components. The centerbody nozzle increases the overall sound pressure level (OASPL) in the near field with respect to the clean nozzle and the centerbody-plus-strut nozzle reduces it compared to the centerbody nozzle due to the increased turbulent mixing. The centerbody perturbed nozzle configurations generate a remarkable spectral peak at S t=0.56 which also occurs in the APE findings in the near field region. This tone is generated by large scale vortical structures shed from the centerbody. The analysis of the individual noise sources shows that the entropy term possesses the highest acoustic contribution in the sideline direction whereas the vortex sound source dominates the downstream acoustics.     

Large-Scale Turbulent Vortical Structures inside a Sudden Expansion Cylinder Chamber

A large eddy simulation (LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model and the Lagrangian dynamic subgrid-scale model are employed and tested. The calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structures behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent behavior behind the flame holder is analyzed by visualizing the sectional views of vortical structure. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effects of free‐stream turbulence on large‐scale coherent structures of separated boundary layer transition

It has been well established that large‐scale structures, usually called coherent structures, exist in many transitional and turbulent flows. The topology and range of scales of those large‐scale structures vary from flow to flow such as counter‐rotating vortices in wake flows, streaks and hairpin vortices in turbulent boundary layer. There has been relatively little study of large‐scale structures in separated and reattached transitional flows. Large‐eddy simulation (LES) is employed in the current study to investigate a separated boundary layer transition under 2% free‐stream turbulence on a flat plate with a blunt leading edge. The Reynolds number based on the inlet free stream velocity and the plate thickness is 6500. A dynamic subgrid‐scale model is employed to compute the subgrid‐scale stresses more accurately in the current transitional flow case. Flow visualization has shown that the Kelvin–Helmholtz rolls, which have been so clearly visible under no free‐stream turbulence (NFST) are not as apparent in the present study. The Lambda‐shaped vortical structures which can be clearly seen in the NFST case can hardly be identified in the free‐stream turbulence (FST) case. Generally speaking, the effects of free‐stream turbulence have led to an early breakdown of the boundary layer, and hence increased the randomization in the vortical structures, degraded the spanwise coherence of those large‐scale structures. Copyright © 2005 John Wiley & Sons, Ltd.     

Assessment of LES Subgrid-scale Models and Investigation of Hydrodynamic Behaviour for an Axisymmetrical Bluff Body Flow

This work is concerned with the investigation of fluid-mechanical behaviour and the performance of different subgrid-scale models for LES in the numerical prediction of a confined axisymmetrical bluff-body flow. Four subgrid-scale turbulence models comprising the Smagorinsky model, Dynamic Smagorinsky model, WALE model and subgrid turbulent kinetic energy model, are validated and compared directly against the experimental data. Two different mesh counts are used for the LES studies, one with a higher mesh resolution in the shear layer than the other. It is found that increasing the mesh resolution improves the time-averaged fluctuating velocity profiles, but has less effect on the time-averaged filtered velocity profiles. A comparison against experiment shows that the recirculation zone length is well predicted using LES. The accuracy of the four different subgrid scale models is then assessed by comparing the LES results using the dense mesh with the experiment. Comparisons with the time-averaged axial and radial velocity profiles demonstrate that LES displays good agreement with the experimental data, with the essential flow features captured both qualitative and quantitatively. The subgrid velocity also matches well with the experimental results, but a slight underprediction of the inner shear layer is observed for all subgrid models. In general, it is found that the Smagorinsky and WALE models are more dissipative than the Dynamic Smagorinsky model and subgrid TKE model. Comparison of the spectra against the experiment shows that LES can capture dominant features of the turbulent flow with reasonable accuracy, and weak spectral peaks related to the Kevin-Helmholtz instability and helical vortex shedding are present.     

Wave–mean flow interaction in an MHD wake behind bluff body

Experimental studies of the effect of constant magnetic field on the process of mean velocity profile stabilization in a wake behind a bluff body are described. To interpret the obtained results, a theoretical model is proposed explaining the scheme of wave-mean flow interaction. We assume that the stabilization process is based on the injection of energy of respective turbulent modes into vortical structures through non-attenuating inertial waves generated in these vortical structures. We also take into account helical character of turbulence in the conditions under study, with weak energy dissipation and its accumulation in large scales.     

Hybrid LES–RANS technique based on a one-equation near-wall model

In order to reduce the high computational effort of wall-resolved large-eddy simulations (LES), the present paper suggests a hybrid LES–RANS approach which splits up the simulation into a near-wall RANS part and an outer LES part. Generally, RANS is adequate for attached boundary layers requiring reasonable CPU-time and memory, where LES can also be applied but demands extremely large resources. Contrarily, RANS often fails in flows with massive separation or large-scale vortical structures. Here, LES is without a doubt the best choice. The basic concept of hybrid methods is to combine the advantages of both approaches yielding a prediction method, which, on the one hand, assures reliable results for complex turbulent flows, including large-scale flow phenomena and massive separation, but, on the other hand, consumes much fewer resources than LES, especially for high Reynolds number flows encountered in technical applications. In the present study, a non-zonal hybrid technique is considered (according to the signification retained by the authors concerning the terms zonal and non-zonal), which leads to an approach where the suitable simulation technique is chosen more or less automatically. For this purpose the hybrid approach proposed relies on a unique modeling concept. In the LES mode a subgrid-scale model based on a one-equation model for the subgrid-scale turbulent kinetic energy is applied, where the length scale is defined by the filter width. For the viscosity-affected near-wall RANS mode the one-equation model proposed by Rodi et al. (J Fluids Eng 115:196–205, 1993) is used, which is based on the wall-normal velocity fluctuations as the velocity scale and algebraic relations for the length scales. Although the idea of combined LES–RANS methods is not new, a variety of open questions still has to be answered. This includes, in particular, the demand for appropriate coupling techniques between LES and RANS, adaptive control mechanisms, and proper subgrid-scale and RANS models. Here, in addition to the study on the behavior of the suggested hybrid LES–RANS approach, special emphasis is put on the investigation of suitable interface criteria and the adjustment of the RANS model. To investigate these issues, two different test cases are considered. Besides the standard plane channel flow test case, the flow over a periodic arrangement of hills is studied in detail. This test case includes a pressure-induced flow separation and subsequent reattachment. In comparison with a wall-resolved LES prediction encouraging results are achieved.      

湍流边界层拟序结构的实验研究

20世纪60年代后, 先后从流动显示发现了快慢斑、猝发、上升流、下扫流和多种涡结构等湍流边界层的拟序结构. 它们对湍流边界层的摩阻、传热传质和湍动能的产生等特性有重要影响. 涡结构是上述拟序结构的核心, 它影响其它拟序结构的发展和演变. 发卡涡通常被认为是基本涡结构. 发卡涡等涡结构的再生, 是湍流边界层拟序结构能够自持续的必要的因素.壁面低速流上升产生猝发, 是湍流边界层湍能的主要来源; 条件采样是测量猝发频率和其它拟序结构出现频率的重要手段. 流动显示对湍流边界层拟序结构作了大量定性观察, 有许多减阻和增加传热率等应用性研究在此基础上发展起来. 80年代后, 出现了测量湍流边界层的瞬时流速矢量场的多热线法和PIV技术, 三维PIV技术可望将来为湍流边界层的实验研究带来重大进展. 本文评述了流动显示法、多热线法和PIV技术的优点和不足之处, 以及它们在对湍流边界层拟序结构的研究中的贡献.      

Evolution analysis of the main mechanisms of the jet/vortex interaction

A new analysis method of the jet/vortex interaction is presented to better understand the effect of each vortical structure on the others composing the jet and even on itself. Vortical structures are isolated from each other and the velocity field they produce is calculated by means of the Biot–Savart law. The structure influence on each another is then distinguished from the others by using its velocity field in the calculation of the vorticity momentum equations. Based on 3D‐datafields provided by a large‐eddy‐simulation (LES) of the interaction between a Lamb–Oseen vortex and a turbulent round jet, a post‐processing tool educes the vortical structures and quantifies the different vorticity momentum equations terms. The study has been limited to the first steps of the interaction regime when there are still few vortical structures so that the subsequent analysis of the post‐process results remains simple enough to be performed. The present method confirmed the intensification of azimuthal structures as they roll up around the main vortex. Core disturbance of the latter was found to be initially caused by those structures and not by core dynamic instabilities, which become dominant with the appearance of non‐linear mechanisms. Copyright © 2010 John Wiley & Sons, Ltd.     

Large eddy simulation of fuel injection and mixing process in a diesel engine

The large eddy simulation(LES) approach implemented in the KIVA-3V code and based on one-equation sub-grid turbulent kinetic energy model are employed for numerical computation of diesel sprays in a constant volume vessel and in a Caterpillar 3400 series diesel engine.Computational results are compared with those obtained by an RANS(RNG k-ε) model as well as with experimental data.The sensitivity of the LES results to mesh resolution is also discussed.The results show that LES generally provides flow and spray characteristics in better agreement with experimental data than RANS;and that small-scale random vortical structures of the in-cylinder turbulent spray field can be captured by LES.Furthermore,the penetrations of fuel droplets and vapors calculated by LES are larger than the RANS result,and the sub-grid turbulent kinetic energy and sub-grid turbulent viscosity provided by the LES model are evidently less than those calculated by the RANS model.Finally,it is found that the initial swirl significantly affects the spray penetration and the distribution of fuel vapor within the combustion chamber.     

Shear‐improved Smagorinsky modeling of turbulent channel flow using generalized Lattice Boltzmann equation

Generalized Lattice Boltzmann equation (GLBE) was used for computation of turbulent channel flow for which large eddy simulation (LES) was employed as a turbulence model. The subgrid‐scale turbulence effects were simulated through a shear‐improved Smagorinsky model (SISM), which is capable of predicting turbulent near wall region accurately without any wall function. Computations were done for a relatively coarse grid with shear Reynolds number of 180 in a parallelized code. Good numerical stability was observed for this computational framework. The results of mean velocity distribution across the channel showed good correspondence with direct numerical simulation (DNS) data. Negligible discrepancies were observed between the present computations and those reported from DNS for the computed turbulent statistics. Three‐dimensional instantaneous vorticity contours showed complex vortical structures that appeared in such flow geometries. It was concluded that such a framework is capable of predicting accurate results for turbulent channel flow without adding significant complications and the computational cost to the standard Smagorinsky model. As this modeling was entirely local in space it was therefore adapted for parallelization. Copyright © 2010 John Wiley & Sons, Ltd.     

Experimental and high-order LES analysis of the flow in near-wall region of a square cylinder

A coupled experimental/numerical analysis of turbulent flow past a square cylinder is performed at the ERCOFTAC Reynolds number Re = U_∞D/ν = 21,400, where U_∞ is the inflow velocity and D the cylinder height. Complementary Laser Doppler Velocimetry (LDV) and high-order large-eddy simulations (LES) approaches, based on a spectral vanishing technique (SVV-LES), provide a comprehensive data base including both instantaneous data and post-processed statistics. Beyond these results, an achievement of the paper is to investigate the coherent structures developing on the sides and in the wake of the cylinder with a special focus on the flow features in the near-wall region. The flow is found to separate at the leading edge of the cylinder with the occurence of three-dimensional Kelvin-Helmholtz (KH) pairings localized in the separating shear layer. The interaction between these KH vortical structures and Von Kármán vortex shedding (VK) in the near wake is discussed based on both visualisations and frequency analysis. In particular, signatures of VK and KH vortical structures are found on velocity time samples.     

Perturbed vortical layers and shear sheltering

New theoretical results and physical interpretations are presented concerning the interactions between different types of velocity fields that are separated by thin interfacial layers, where there are dynamically significant variations of vorticity across the layers and, in some cases within them. It is shown how, in different types of complex engineering and environmental flow, the strengths of these interactions vary from the weakest kind of superposition to those where they determine the flow structure, for example by mutual exclusion of velocity fields from the other region across the interface, or by local resonance near the interface. We focus here on the excluding kinds of interactions between, on the one hand, elongated and compact regions containing vortical flows and large variations in velocity, and on the other hand various kinds of weak perturbation in the surrounding external flow region: rotational, irrotational; time-varying, steady; large, small; coplanar, non-coplanar; non-diffusive, diffusive. It is shown how all these kinds of external disturbances can be wholly, or partially, ‘blocked’ at the interface with the vortical region, so that beyond a certain sheltering distance into the interior of this region the fluctuations can be very small. For the special case of quasi-parallel co-planar external straining motions outside non-directional shear flows, weak sheltering occurs if the mean velocity of the shear flow increases – otherwise the perturbations are amplified. For non-parallel flows, the sheltering effect can be greater when the vorticity is distributed in thin vortex sheets. The mechanism whereby the vortical flow induces ‘blocking’ and ‘shear-sheltering’ effects can be quantitatively explained in terms of the small adjustments of the vorticity in the vortical layers, and in some cases by the change in impulse of these layers. If the vorticity in the outer part of the vortical region is weak, it can be ‘stripped away’ by the external disturbances until the remaining vorticity is strong enough to ‘block’ the disturbances and shelter the inner flow of the vortical region. The mechanisms presented here appear to explain on the one hand some aspects of the observed robustness of vortical structures and jet or plume like shear flows in turbulent and geophysical flows, and on the other hand the levels of external perturbation needed to erode or breakdown turbulent shear flows.     

PIV analysis of near-wake behind a sphere at a subcritical Reynolds number

The vortical structure of near-wake behind a sphere is investigated using a PIV technique in a circulating water channel at Re = 11,000. The measured velocity fields show a detailed vortical structure in the recirculation region such as recirculation vortices, reversed velocity zone, and out-of-plane vorticity distribution. The vorticity distribution of the sphere wake shows waviness in cross-sectional planes. The time-averaged turbulent structures are consistent with the visualized flow showing the onset of shear layer instability. The spatial distributions of turbulent intensities provide turbulent statistics for validating numerical predictions.     

Numerical prediction of turbulent flow structure generated by a synthetic cross‐jet into a turbulent boundary layer

A detailed numerical study using large‐eddy simulation (LES) and unsteady Reynolds‐averaged Navier–Stokes (URANS) was undertaken to investigate physical processes that are engendered in the injection of a circular synthetic (zero‐net mass flux) jet in a zero pressure gradient turbulent boundary layer. A complementary study was carried out and was verified by comparisons with the available experimental data that were obtained at corresponding conditions with the aim of achieving an improved understanding of fluid dynamics of the studied processes. The computations were conducted by OpenFOAM C++, and the physical realism of the incoming turbulent boundary layer was secured by employing random field generation algorithm. The cavity was computed with a sinusoidal transpiration boundary condition on its floor. The results from URANS computation and LES were compared and described qualitatively and quantitatively. There is a particular interest for acquiring the turbulent structures from the present numerical data. The numerical methods can capture vortical structures including a hairpin (primary) vortex and secondary structures. However, the present computations confirmed that URANS and LES are capable of predicting current flow field with a more detailed structure presented by LES data as expected. Copyright © 2011 John Wiley & Sons, Ltd.     

On the performance of WENO/TENO schemes to resolve turbulence in DNS/LES of high‐speed compressible flows

High‐speed compressible turbulent flows typically contain discontinuities and have been widely modeled using Weighted Essentially Non‐Oscillatory (WENO) schemes due to their high‐order accuracy and sharp shock capturing capability. However, such schemes may damp the small scales of turbulence and result in inaccurate solutions in the context of turbulence‐resolving simulations. In this connection, the recently developed Targeted Essentially Non‐Oscillatory (TENO) schemes, including adaptive variants, may offer significant improvements. The present study aims to quantify the potential of these new schemes for a fully turbulent supersonic flow. Specifically, DNS of a compressible turbulent channel flow with M = 1.5 and Re_τ = 222 is conducted using OpenSBLI, a high‐order finite difference computational fluid dynamics framework. This flow configuration is chosen to decouple the effect of flow discontinuities and turbulence and focus on the capability of the aforementioned high‐order schemes to resolve turbulent structures. The effect of the spatial resolution in different directions and coarse grid implicit LES are also evaluated against the WALE LES model. The TENO schemes are found to exhibit significant performance improvements over the WENO schemes in terms of the accuracy of the statistics and the resolution of the three‐dimensional vortical structures. The sixth‐order adaptive TENO scheme is found to produce comparable results to those obtained with nondissipative fourth‐ and sixth‐order central schemes and reference data obtained with spectral methods. Although the most computationally expensive scheme, it is shown that this adaptive scheme can produce satisfactory results if used as an implicit LES model.     

Large-eddy simulation of counter-rotating Taylor–Couette flow: The effects of angular velocity and eccentricity

In the present work, turbulent flow in the annulus of a counter-rotating Taylor-Couette (CRTC) system is studied using large-eddy simulation. The numerical methodology employed is validated, for both the mean and second-order statistics, with the direct numerical simulation (DNS) data available in the literature, for a range of Reynolds numbers from 500 to 4000. Thereafter, turbulent flow occurring in this system at Reynolds numbers of 8000 and 16000 are studied, and the results obtained are analyzed using mean and second-order statistics, vortical structures, velocity vector plots and power energy spectra. Further, the spatio-temporal variation of azimuthal velocity, extracted near the inner cylinder, shows the existence of herringbone like patterns similar to that observed in the previous studies. The effect of eccentricity of the inner cylinder with respect to the outer cylinder is studied, on the turbulent flow in the CRTC system, for two different eccentricity ratios of 0.2 and 0.5 and for two different Reynolds numbers of 1500 and 4000. The results of the eccentric CRTC are analyzed using contours of pressure, mean and second-order statistics, velocity vectors, vortical structures, and turbulence anisotropy maps. It is observed from the eccentric CRTC simulations that the smaller-gap region seems to contain higher amplitude fluctuations and more vortical structures when compared with the larger-gap region. The mean turbulent kinetic energy contours do not change qualitatively with the Reynolds number, however, quantitatively a higher turbulent kinetic energy is observed in the higher Reynolds number case of 4000.     

A DNS study on the effects of convex streamwise curvature on coherent structures in a temporally-developing turbulent boundary layer with supercritical water

Direct numerical simulation (DNS) results are used to establish the effect of convex streamwise curvature on the development of turbulent boundary layers, and the effect of such curvature on the forced-convection heat transfer variations observed at certain supercritical thermodynamic states. The results illustrate the stabilizing effects of this flow geometry through modification of the structure and distribution of hairpin-like vortical flow structures in the boundary layer. Furthermore, enhancement of convective heat transfer realized at a particular heat flux-to-mass flux ratio with the working fluid at a supercritical state is observed to be reduced by the stabilizing effect of convex surface curvature.     

Numerical investigation of turbulent flow,heat transfer and entropy generation in a helical coiled tube with larger curvature ratio

Three-dimensional turbulent forced convective heat transfer and flow characteristics, and the non-dimensional entropy generation number in a helical coiled tube subjected to uniform wall temperature are simulated using the k–ε standard turbulence model. A finite volume method is employed to solve the governing equations. The effects of Reynolds number, curvature ratio, and coil pitch on the average friction factor and Nusselt number are discussed. The results presented in this paper cover a Reynolds number range of 2 × 10⁴ to 6 × 10⁴, a pitch range of 0.1–0.2 and a curvature ratio range of 0.1–0.3. The results show that the coil pitch, curvature ratio and Reynolds number have different effects on the average friction factor and Nusselt number at different cross-sections. In addition, the flow and heat transfer characteristics in a helical coiled tube with a larger curvature ratio for turbulent flow are different from that of smaller curvature ratio for laminar and turbulent flow in certain ways. Some new features that are not obtained in previous researches are revealed. Moreover, the effects of Reynolds number, curvature ratio, and coil pitch on the non-dimensional entropy generation number of turbulent forced convection in a helical coiled tube are also discussed.     

Anisotropic Turbulence and Coherent Structures in Eccentric Annular Channels

Eccentric annular pipe flows represent an ideal model for investigating inhomogeneous turbulent shear flows, where conditions of turbulence production and transport vary significantly within the cross-section. Moreover recent works have proven that in geometries characterized by the presence of a narrow gap, large-scale coherent structures are present. The eccentric annular channel represents, in the opinion of the present authors, the prototype of these geometries. The aim of the present work is to verify the capability of a numerical methodology to fully reproduce the main features of the flow field in this geometry, to verify and characterize the presence of large-scale coherent structures, to examine their behavior at different Reynolds numbers and eccentricities and to analyze the anisotropy associated to these structures. The numerical approach is based upon LES, boundary fitted coordinates and a fractional step algorithm. A dynamic Sub Grid Scale (SGS) model suited for this numerical environment has been implemented and tested. An additional interest of this work is therefore in the approach employed itself, considering it as a step into the development of an effective LES methodology for flows in complex channel geometries. Agreement with previous experimental and DNS results has been found good overall for the streamwise velocity, shear stress and the rms of the velocity components. The instantaneous flow field presented large-scale coherent structures in the streamwise direction at low Reynolds numbers, while these are absent or less dominant at higher Reynolds and low eccentricity. After Reynolds averaging is performed over a long integration time the existence of secondary flows in the cross session is proven. Their shape is found to be constant over the Reynolds range surveyed, and dependent on the geometric parameters. The effect of secondary flows on anisotropy is studied over an extensive Reynolds range through invariant analysis. Additional insight on the mechanics of turbulence in this geometry is obtained.