Large-eddy simulation study of an air flow past a heated square cylinder

A large-eddy simulation (LES) technique has been used for the calculation of an air flow past a heated square cylinder. The LES method is a conventional one with the Smagorinsky eddy-viscosity model, and the computational grid is small enough to be handled by workstations. The computed turbulent flow field quantities agree well with the experimental results reported by Lyn et al. (J Fluid Mech 304:285–319, 1995). Particular attention has been spent to the convective heat transfer and the prediction of the thermal fluctuations in case of a fluid with Pr = 0.71. The LES results agree well with the empirical correlations proposed by Hilpert and Sparrow et al. for the average Nusselt number.     

Smagorinsky constant in LES modeling of anisotropic MHD turbulence

Turbulent fluctuations in magnetohydrodynamic (MHD) flows can become strongly anisotropic or even quasi-2D under the action of an applied magnetic field. We investigate this phenomenon in the case of low magnetic Reynolds numbers. It has been found in earlier DNS and LES of homogeneous turbulence that the degree of anisotropy is predominantly determined by the value of the magnetic interaction parameter and only slightly depends on the Reynolds number, type of large-scale dynamics, and the length scale. Furthermore, it has been demonstrated that the dynamic Smagorinsky model is capable of self-adjustment to the effects of anisotropy. In this paper, we capitalize on these results and propose a simple and effective generalization of the traditional non-dynamic Smagorinsky model to the case of anisotropic MHD turbulence.      

Validation of large eddy simulation method for study of flatness and skewness of decaying compressible magnetohydrodynamic turbulence

In the present work we study potential applicability of large eddy simulation (LES) method for prediction of flatness and skewness of compressible magnetohydrodynamic (MHD) turbulence. The knowledge of these quantities characterizes non-Gaussian properties of turbulence and can be used for verification of hypothesis on Gaussianity for the turbulent flow under consideration. Prediction accuracy of these quantities by means of LES method directly determines efficiency of reconstruction of probability density function (PDF) that depends on used subgrid-scale (SGS) parameterizations. Applicability of LES approach for studying of PDF properties of turbulent compressible magnetic fluid flow is investigated and potential feasibilities of five SGS parameterizations by means of comparison with direct numerical simulation results are explored. The skewness and the flatness of the velocity and the magnetic field components under various hydrodynamic Reynolds numbers, sonic Mach numbers, and magnetic Reynolds numbers are studied. It is shown that various SGS closures demonstrate the best results depending on change of similarity numbers of turbulent MHD flow. The case without any subgrid modeling yields sufficiently good results as well. This indicates that the energy pile-up at the small scales that is characteristic for the model without any subgrid closure, does not significantly influence on determination of PDF. It is shown that, among the subgrid models, the best results for studying of the flatness and the skewness of velocity and magnetic field components are demonstrated by the Smagorinsky model for MHD turbulence and the model based on cross-helicity for MHD case. It is visible from the numerical results that the influence of a choice subgrid parametrization for the flatness and the skewness of velocity is more essential than for the same characteristics of magnetic field.     

基于PIV测量的超声波流量计内流场特性研究

采用PIV(Particle Image Velocimetry)测量手段,考察了小口径超声波流量计的流动特性。首先针对前端安装直管段时,不同流量条件下的流场特性建立基本认识,实验结果表明,在低流量条件下,流量计内流场存在明显的不稳定演变和非定常流动特征。进一步以上游前端安装球阀为典型案例,考察了安装条件对超声波流量计响应特性和测量偏差的影响。结合直管段的实验观测结果,发现此种结构超声波流量计的适应性与其流场非定常性的关系具有很好的一致性,即流场结构稳定则适应性强。此外,综合多参数的实验结果表明,雷诺数是判断小口径超声波流量计测量准确性的重要无量纲参数。     

Turbulent Flow and Dispersion of Inertial Particles in a Confined Jet Issued by a Long Cylindrical Pipe

In this work we examine first the flow field of a confined jet produced by a turbulent flow in a long cylindrical pipe issuing in an abrupt angle diffuser. Second, we examine the dispersion of inertial micro-particles entrained by the turbulent flow. Specifically, we examine how the particle dispersion field evolves in the multiscale flow generated by the interactions between the large-scale structures, which are geometry dependent, with the smaller turbulent scales issued by the pipe which are advected downstream. We use Large-Eddy-Simulation (LES) for the flow field and Lagrangian tracking for particle dispersion. The complex shape of the domain is modelled using the immersed-boundaries method. Fully developed turbulence inlet conditions are derived from an independent LES of a spatially periodic cylindrical pipe flow. The flow field is analyzed in terms of local velocity signals to determine spatial coherence and decay rate of the coherent K–H vortices and to make quantitative comparisons with experimental data on free jets. Particle dispersion is analyzed in terms of statistical quantities and also with reference to the dynamics of the coherent structures. Results show that the particle dynamics is initially dominated by the Kelvin–Helmholtz (K–H) rolls which form at the expansion and only eventually by the advected smaller turbulence scales.     

A numerical study of flows driven by a rotating magnetic field in a square container

The magnetically induced fluid flow in a square container is investigated by means of numerical simulations. Low frequency/ low induction conditions are assumed. The effect of the rotating magnetic field gives rise to a time-independent magnetic body force, computed via the electrical potential equation and Ohm's law and a time-dependent part that is neglected due to the low-interaction parameter. The magnetic body force calculation is verified successfully by comparison with the exact solution. The behavior of the fluid flow in the square container reveals similar features to the flow in the cylindrical container, for instance, in the dependence on the intensities of the magnetic field. However, we did find differences in the velocity field distribution. Particularly, in the finite as well as infinite geometry, the velocity field is influenced by the corner of the container and remains non-axisymmetric in a wide range of Taylor numbers.     

Spectral-based simulations of particle-laden turbulent flows

In this paper, we discuss the application of spectral-based methods to simulation of particle-laden turbulent flows. The primary focus of the article is on the past and ongoing works by the authors. The particles are tracked in Lagrangian framework, while direct numerical simulation (DNS) or large-eddy simulation (LES) is used to describe the carrier-phase flow field. Two different spectral methods are considered, namely Fourier pseudo-spectral method and Chebyshev multidomain spectral method. The pseudo-spectral method is used for the simulation of homogeneous turbulence. DNS of both incompressible and compressible flows with one- and two-way couplings are reported. For LES of particle-laden flows, two new models, developed by the authors, account for the effect of sub-grid fluctuations on the dispersed phase. The Chebyshev multidomain method is employed for the works on inhomogeneous flows. A number of canonical flows are discussed, including flow past a square cylinder, channel flow and flow over backward-facing step. Ongoing research on particle-laden LES of inhomogeneous flows is briefly reported.     

Thermal instability in natural convection flow over a boundary layer subject to external electrical and magnetic fields

This paper presents a numerical study of external electrical and magnetic effects on the formation of longitudinal vortices in natural convection flow over a heated horizontal plate. The criterion on the position marking on the onset of longitudinal vortices is defined in the present paper. The onset position characterized by the Grashof number depends on the Prandtl number, the wave number, the electric field parameter, and the magnetic field parameter. The flow is found more stable as the value of the magnetic field parameter increases. Moreover, the stabilizing effect is also found on the flow when the positive electric field parameter Ec is applied. The results of the present numerical prediction show reasonable agreement with the experimental data with zero magnetic field and electric field parameters in literature.     

Direct measurement of the stream-function in a quasi-two-dimensional liquid metal flow

 Flows of electrically conducting fluids become quasi-two-dimensional (2D) under the influence of a homogeneous magnetic field. In this regime the electrical potential induced in the flow is directly proportional to the 2D stream function. We utilize this property in order to perform simultaneous multi-channel measurements of the stream function in an electromagnetically driven annular liquid metal flow in the weakly nonlinear regime. Spatio-temporal maps of the stream functions of a quasi-periodic unstable flow are observed and used as a basis for a systematic decomposition of the flow into spatial modes belonging to different frequencies. The approach can be generalized to flows with complex time-dependence. Received: 3 March 1997/Accepted: 27 September 1997     

Solution of the direct problem of the mixed flow of a conducting gas in a laval nozzle with a meridional magnetic field

We consider the direct problem in the theory of the axisymmetric Laval nozzle (including sonic transition) for the steady flow of an inviscid and nonheat-conducting gas of finite electrical conductivity. The problem is solved by numerical integration of the equations of unsteady gas flow using an explicit difference scheme that was proposed by Godunov [1,2], and was used to calculate steady and unsteady flows of a nonconducting gas in nozzles by Ivanov and Kraiko [3]. The subsonic and the supersonic flows of a conducting gas in an axisymmetric channel when there is no external electric field, the magnetic field is meridional, and the magnetic Reynolds numbers are small have previously been completely investigated. Thus, Kheins, Ioller and Élers [4] investigated experimentally and theoretically the flow of a conducting gas in a cylindrical pipe when there is interaction between the flow and the magnetic field of a loop current that is coaxial with the pipe. Two different approaches were used in the theoretical analysis in [4]: linearization with respect to the parameter S of the magnetogasdynamic interaction and numerical calculation by the method of characteristics. The first approach was used for weakly perturbed subsonic and supersonic flows and the solutions obtained in analytic form hold only for small S. This is the approach used by Bam-Zelikovich [5] to investigate subsonic and supersonic jet flows through a current loop. The numerical calculations of supersonic flows in a cylindrical pipe in [4] were restricted to comparatively small values of S since, as S increases, shock waves and subsonic waves appear in the flow. Katskova and Chushkin [6] used the method of characteristics to calculate the flow of the type in the supersonic part of an axisymmetric nozzle with a point of inflection. The flow at the entrance to the section of the nozzle under consideration was supersonic and uniform, while the magnetic field was assumed to be constant and parallel to the axis of symmetry. The plane case was also studied in [6]. The solution of the direct problem is the subject of a paper by Brushlinskii, Gerlakh, and Morozov [7], who considered the flow of an electrically conducting gas between two coaxial electrodes of given shape. There was no applied magnetic field, and the induced magnetic field was in the direction perpendicular to the meridional plane. The problem was solved numerically in [7] using a standard process. However, the boundary conditions adopted, which were chosen largely to simplify the calculations, and the accuracy achieved only allowed the authors [7] to make reliable judgments about the qualitative features of the flow. Recently, in addition to [7], several papers have been published [8–10] in which the authors used a similar approach to solve the direct problem in the theory of the Laval nozzle (in the case of a nonconducting gas).Translated from Izvestiya Akademiya Nauk SSSR, Mekhanika Zhidkosti i Gaza., No. 5, pp. 14–20, September–October, 1971.In conclusion the author wishes to thank M. Ya. Ivanov, who kindly made available his program for calculating the flow of a conducting gas, and also A. B. Vatazhin and A. N. Kraiko for useful advice.     

Turbulent heat transfer characteristics of water flow in a rotating pipe

Large-Eddy-Simulation of turbulent heat transfer for water flow in rotating pipe is performed, for various rotation ratios (0 ≤ N ≤ 14). The value of the Reynolds number, based on the bulk velocity and pipe diameter, is Re = 5,500. The aim of this study is to examine the effect of the rotating pipe on the turbulent heat transfer for water flow, as well as the reliability of the LES approach for predicting turbulent heat transfer in water flow. Some predictions for the case of non-rotating pipe are compared to the available results of literature for validation. To depict the influence of the rotation ratio on turbulent heat transfer, many statistical quantities are analyzed (distributions of mean temperature, rms of fluctuating temperature, turbulent heat fluxes, higher-order statistics). Some contours of instantaneous temperature fluctuations are examined.     

Vortex mechanism of heat transfer enhancement in a channel with spherical and oval dimples

Vortex mechanism of heat transfer enhancement in a narrow channel with dimples has been investigated numerically using LES and URANS methods. The flow separation results in a formation of vortex structures which significantly enhance heat transfer on dimpled surfaces leading to a small increase in pressure loss. The heat transfer can be significantly increased by rounding the dimple edge and use of oval dimples. To get a deep insight into flow physics LES is performed for single phase flow in a channel with a spherical dimple. The instantaneous vortex formation and separation are investigated in and around the dimple area. Considered are Reynolds numbers (based on dimple print diameter) Re_D = 20,000 and Re_D = 40,000 the depth to print diameter ratio of Δ = 0.26. Frequency analysis of LES data revealed the presence of dominating frequencies in unsteady flow oscillations. Direct analysis of the flow field revealed the presence of coherent vortex structure inclined to the mean flow. The structure changes its orientation in time causing the long period oscillations with opposite-of-phase motion. Three dimensional proper orthogonal decomposition (POD) analysis is carried out on LES pressure and velocity fields to identify spatio-temporal structures hidden in the random fluctuations. Tornado-like spatial POD structures have been determined inside dimples.     

Euler–Euler large eddy simulations for dispersed turbulent bubbly flows

In this paper we present detailed Euler–Euler Large Eddy Simulations (LES) of dispersed bubbly flow in a rectangular bubble column. The motivation of this study is to investigate the potential of this approach for the prediction of bubbly flows, in terms of mean quantities. The physical models describing the momentum exchange between the phases including drag, lift and wall force were chosen according to previous experiences of the authors. Experimental data, Euler–Lagrange LES and unsteady Euler–Euler Reynolds-Averaged Navier–Stokes results are used for comparison. It is found that the present model combination provides good agreement with experimental data for the mean flow and liquid velocity fluctuations. The energy spectrum obtained from the resolved velocity of the Euler–Euler LES is presented as well.     

Thermal instability of a boundary layer subjected to external electric and magnetic fields

This paper presents a numerical study of the electrical and magnetic fields on thermal instability in a boundary layer. The criterion on the position marking on the onset of longitudinal vortices is defined in the present paper. The results show that the onset position characterized by the Grashof number depends on the Prandtl number, the Reynolds number, the wave number, the electric field parameter, and the Hartmann number. The flow becomes more stable as the magnetic field increases. However, the destabilizing effect is found on the flow when the negative electric field parameter is applied. The results of the present numerical prediction show reasonable agreement with the experimental data in the case of zero Hartmann number and zero electric field parameter in the literature. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical Evaluation of Combustion Regimes in a GDI Engine

There is significant interest in the gasoline direct-injection engine due to its potential for improvements in fuel consumption but it still remains an area of active research due to a number of challenges including the effect of cycle-by-cycle variations. The current paper presents the use of a 3D-CFD model using both the RANS and LES turbulence modelling approaches, and a Lagrangian DDM to model an early fuel injection event, to evaluate the regimes of combustion in a gasoline direct-injection engine. The velocity fluctuations were investigated as an average value across the cylinder and in the region between the spark plug electrodes. The velocity fluctuations near the spark plug electrodes were seen to be of lower magnitude than the globally averaged fluctuations but exhibited higher levels of cyclic variation due to the influence of the spark plug electrode and the pent-roof geometry on the in-cylinder flow field. Differences in the predicted flame structure due to differences in the predicted velocity fluctuations between RANS and LES modelling approaches were seen as a consequence of the inherently higher dissipation levels present in the RANS methodology. The increased cyclic variation in velocity fluctuations near the spark plug electrodes in the LES predictions suggested significant variation in the relative strength of the in-cylinder turbulence and that may subsequently result in a thickening of the propagating flame front from cycle-to-cycle in this region. Throughout this paper, the numerical results were validated against published experimental data of the same engine geometry under investigation.     

The electrohydrodynamics of superimposed fluids subjected to a nonuniform transverse electric field

This study aims to investigate electrohydrodynamics of two superimposed fluids that are confined between a pair of two-dimensional flat plates and are exposed to a sinusoidal electric field in zero gravity. The goal is to identify the parameters that affect the flow structure and interface deformation using a simple closed form solution. The governing electrohydrodynamic equations are solved analytically for Newtonian and immiscible fluids in the framework of leaky-dielectric theory and in the limit of small electric field and fluid inertia. A detailed analysis of the electric and flow fields is presented and it is shown that the electric field induces sinusoidal electrical stresses at the interface, which lead to periodic convection cells. The parameters affecting the sense of flow circulation and strength are investigated and it is shown that the former depends on the relative magnitude of the electric permittivity and conductivity ratios while the latter is controlled by the relative thicknesses of the fluid layers and the ratio of the electric conductivities and viscosities of the fluids. The maximum flow strength is achieved at a relative thickness that is set by the competition between the electric and hydrodynamic effects. For small deformation, the distortion of the interface is examined using a normal stress balance at the interface, and it is shown that the degree of interface deformation scales with the square of the amplitude of the electric potential nonuniformity, while its wavenumber is twice that of the imposed potential nonuniformity. Furthermore, a zero-deformation curve is found, which delineates the region in the permittivity-conductivity space according to the sense of interface deformation. The results show that for certain ranges of fluid layer thicknesses and permittivity ratios, the interface will remain flat, despite the action of the nonuiform field.     

Numerical Investigation of Active Flow Control Around a Generic Truck A-Pillar

Large Eddy Simulations (LES) are conducted to study the actuated flow field around a bluff body. The model is a simplified section of a truck. The aim of the work is to model the separation of the flow acting at the front rounded corners, the so called A-pillars, and to minimize the separation of the flow by means of Zero Net Mass Flux synthetic jets. LES data show the interaction of the flow main structures, the separation mechanism and the effects of the actuation on the flow field. The flow is post processed using modal and frequency decompositions. Relevant results in terms of drag reduction were observed for the actuated flow. The principle flow mechanisms are discussed and an optimal actuation frequency, in terms of induced fluctuations and drag reduction, is identified.     

Error of large-eddy simulation in the wall pressure fluctuation of a turbulent channel flow

We analyze the error of large-eddy simulation(LES) in wall pressure fluctuation of a turbulent channel flow. To separate different sources of the error, we conduct both direct numerical simulations(DNS)and LES, and apply an explicit filter on DNS data to obtain filtered DNS(FDNS) data. The error of LES is consequently decomposed into two parts: The first part is the error of FDNS with respect to DNS,which quantifies the influence of the filter operation. The second part is the difference between LES and FDNS induced by the error of LES in velocity field. By comparing the root-mean-square value and the wavenumber-frequency spectrum of the wall pressure fluctuation, it is found that the inaccuracy of the velocity fluctuations is the dominant source that induces the error of LES in the wall pressure fluctuation.The present study provides a basis on future LES studies of the wall pressure fluctuation.     

Large Eddy Simulation of the sound field of a round turbulent jet

In this paper we will use Large Eddy Simulation (LES) to obtain the flow field of a turbulent round jet at a Reynolds number based on the jet orifice velocity of 11000. In the simulations it is assumed that the flow field is incompressible. The acoustic field of the jet is calculated with help of the Lighthill acoustic analogy. The coupling between the flow solver and the acoustic solver is discussed in detail. The Mach number used in the acoustic calculation was equal to 0.6. It is shown that the decay of the jet centerline velocity and centerline rms are in good agreement with experimental data of [12]. Furthermore, it is shown that the influence of the LES modeling on the acoustic field is very small, if the dynamic subgrid model is used.     

Hybrid LES-RANS using synthesized turbulent fluctuations for forcing in the interface region

The main bottleneck in using Large Eddy Simulations at high Reynolds number is the requirement of very fine meshes near walls. One of the main reasons why hybrid LES-RANS was invented was to eliminate this limitation. In this method unsteady RANS (URANS) is used near walls and LES is used away from walls. The present paper evaluates a method for improving standard LES-RANS. The improvement consists of adding instantaneous turbulent fluctuations (forcing conditions) at the matching plane between the LES and URANS regions in order to trigger the equations to resolve turbulence. The turbulent fluctuations are taken from synthesized homogeneous turbulence assuming a modified von Kármán spectrum. Both isotropic and non-isotropic fluctuations are evaluated. The new approach is applied to fully developed channel flow and it is shown that the imposed fluctuations considerably improve the predictions. It is found that increasing the prescribed turbulent length scale of the synthesized turbulence provides excellent agreement with the classical log-law.