On the breakup of bubbles at high Reynolds numbers and subcritical Weber numbers

In this paper we propose a simplified two-dimensional model to describe some aspects of the turbulent breakage of bubbles at subcritical Weber numbers. In particular we focus on the breakup of bubbles owing to their interaction with an array of successive eddies, modeled by a train of straining flows. Our simulations show that, under certain conditions, a bubble accumulates energy due to its interaction with a sequence of turbulent structures until it eventually breaks, even if none of the eddies is sufficiently energetic to split the particle by itself. It is also shown that the different strain directions of the eddies acting on the surface of the bubble, and the resonance effect between their characteristic frequency and the natural oscillation frequency of the bubble immersed into the straining flow are the two key factors in the bubble deformation, and subsequent breakup mechanism. Moreover, the breakup patterns obtained from our simulations seem to agree qualitatively well with the experimental observations.     

Communication between the buffer layer and the wall in a turbulent channel flow

Direct numerical simulations obtained in large computational domains of a fully developed turbulent channel flow up to the Karman number 1100 are analyzed to determine the scaling of the spanwise correlation coefficients and the effect of the outer eddies. The local fluctuating velocity field is narrow-band-pass and low-pass filtered along the streamwise wavenumber. The spanwise correlations of the narrow-band passed signals in the low buffer layer adequately provide length scales and signatures of the active structures. The low-pass filtering is used to investigate the relative role of the outer eddies. The impact of the active and passive eddies on the wall is analyzed separately through the cross-correlations of the filtered velocity field with the wall shear stress fluctuations. Characteristic length-scales resulting from the analysis of the velocity field differ depending on the quantity and some are related to the conventional streak spacing but not all. The quasi-streamwise vortex paradigm, for the most part, allows the interpretation of these characteristics, but fails in some cases.     

电磁力诱导二次流对微流体混合效果影响的研究

设计了一种微流道中液体混合的方法, 基于电磁力的作用, 诱导微流道中二次流动 发生, 从而产生Dean涡, 使流道中流体出现往复运动及流体界面的弯曲延伸, 使不 同流体的接触面积大大地增加, 从而提高混合效率. 对3种结构的工况进行了比较. 建立了微流道混合系统的理论模型并进行了数值模拟,对流体混合过程进行了分析并对混合效率进行了评价. 进行了验证性的试验研究, 对模拟结果与试验结果进行了比较, 结果相吻合.     

The nonlinear galerkin method: A multiscale method applied to the simulation of homogeneous turbulent flows

Using results of DNS in the case of two-dimensional homogeneous isotropic flows, we first analyze in detail the behavior of the small and large scales of Kolmogorov-like flows at moderate Reynolds numbers. We derive several estimates on the time variations of the small eddies and the nonlinear interaction terms; these terms play the role of the Reynolds stress tensor in the case of LES. Since the time step of a numerical scheme is determined as a function of the energy-containing eddies of the flow, the variations of the small scales and of the nonlinear interaction terms over one iteration can become negligible by comparison with the accuracy of the computation. Based on this remark, we propose a multilevel scheme which treats the small and the large eddies differently. Using mathematical developments, we derive estimates of all the parameters involved in the algorithm, which then becomes a completely self-adaptative procedure. Finally, we perform realistic simulations of (Kolmogorov-like) flows over several eddy-turnover times. The results are analyzed in detail and a parametric study of the nonlinear Galerkin method is performed.This work was partially supported by NFS Grant DMS 9024769 and DOE DE-FG02-92ER25120. The computations were made on the Cray-2 of the National Center for Supercomputing Applications (NCSA, Illinois). Part of the work was done while T. Dubois was visiting the Institute for Computer and Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681-0001, U.S.A.     

Numerical simulations of vortices in mixing layers and plane jets

Direct numerical simulations of two-dimensional mixing layers and jets are presented first within the temporal approximation, and then with spatial calculations. The evolution of a temperature field as a passive scalar is investigated simultaneously. Besides usual statistical quantities, the results are presented through visualisations of coloured vorticity- and temperature-contours in order to describe the evolutions of the fields with time. Concentrated vortices develop and the evolution of the layer is dominated by the interactions between these eddies which can be viewed as the coherent structures observed in laboratory experiments.     

Wavelet spectrum analysis on energy transfer of multi-scale structures in wall turbulence

The streamwise velocity components at different vertical heights in wall turbulence were measured. Wavelet transform was used to study the turbulent energy spectra, indicating that the global spectrum results from the weighted average of Fourier spectrum based on wavelet scales. W＇avelet transform with more vanishing moments can express the declining of turbulent spectrum. The local wavelet spectrum shows that the physical phenomena such as deformation position in the boundary layer, and the or breakup of eddies are related to the vertical energy-containing eddies exist in a multi-scale form. Moreover, the size of these eddies increases with the measured points moving out of the wall. In the buffer region, the small scale energy-containing eddies with higher frequency are excited. In the outer region, the maximal energy is concentrated in the low-frequency large-scale eddies, and the frequency domain of energy-containing eddies becomes narrower.     

Turbulence Modeling for the Stable Atmospheric Boundary Layer and Implications for Wind Energy

The near-surface structure of atmospheric turbulence affects the design and operation of wind turbines and is especially difficult to predict under stably-stratified conditions. This study uses large-eddy simulation (LES) to explore properties of the stable boundary layer (SBL) using an explicit filtering and reconstruction turbulence modeling approach. Simulations of the atmospheric boundary layer over flat terrain, under both moderately and strongly stable conditions are performed. Results from high-resolution simulations are used to investigate SBL flow structures including mean profiles and turbulence statistics, which are relevant to wind energy applications. The applicability of power-law relations and empirical similarity formulations for predicting wind speed depend on the strength of stratification and are shown to be inadequate. Low-level jets form in the simulations. Under strong stability, vertical wind shear below the jet triggers intermittent turbulence. The associated sporadic “bursting” events are extremely energetic and last longer than the time scale of the largest eddies. Such phenomena can have adverse effects on turbine lifetime and performance.     

Numerical experiments on stability and transition in plane channel flow

Various secondary and tertiary instabilities in plane channel flow are explored via time-dependent numerical simulations using the incompressible Navier-Stokes equations. Comparisons are made between transitional flows at Reynolds numbers 1500, 5000, and 8000. The lambda vortex, detached shear layer, and inverted vortex regions are identified and the origin of the latter is explained. The laminar breakdown of the Re=1500 flow is computed with high resolution and the nature of its ensuing hairpin eddies is clarified by numerical particle paths. The potential of center-mode rather than wall-mode transitions is proposed and the resulting flow structure is described.     

Statistical Structure at the Wall of the High Reynolds Number Turbulent Boundary Layer

The one and two-point statistical structure of very high Reynolds number turbulence in the surface layer near a rigid `wall' is analysed. The essential mechanisms for turbulent eddies impinging on the wall are studied using linearised rapid distortion theory, which show how the mean shear and blocking actions of the surface act first independently and then, over the life time of the eddy, interactively. Previous analytical results are reinterpreted and some new results are derived to show how the integral length scales, cross correlations and spectra of the different components of the turbulence are distorted depending on the form of the spectra of eddies above the surface layer and how they are related to motions of characteristic eddy structures near the surface. These results are applied to derive some quantitative and qualitative predictions in the surface layers (SL), where the eddies are affected by local shear dynamics, and in the `eddy surface layer' (ESL) where quasi independents loping elongated eddies interact directly with the wall, and where there is a large range of wave number within which the spectra of the horizontal velocity components are proportional to k ⁻¹. The longest eddies in the boundary layer occur near the wall. Field experiments agree with the theoretical model predictions for the quite different forms for the spectra, cospectra and cross correlations for the vertical and horizontal components of the velocity field. By showing that in SL the energy exchange between the large and small scale eddies is local(`staircase') energy cascade, whereas in ESL there is a direct nonlocal (`elevator-like')energy transfer to the small scales, it is shown why the thickness of the ESL increases over rougher surfaces and as the Reynolds number decreases. This revised version was published online in August 2006 with corrections to the Cover Date.     

The size effect on the dispersion of a particle in a homogeneous isotropic turbulence

The mechanism of the response motion of a suspended particle to turbulent motion of its surrounding fluid is different according to size of turbulent eddies. The particle is dragged by the viscous force of large eddies, and meanwhile driven randomly by small eddies. Based on this understanding, the dispersion of a particle with finite size in a homogeneous isotropic turbulence is calculated in this study. Results show that there are two competing effects: when enhanced by the inertia of a particle, the long-term particle diffusivity is reduced by the finite size of the particle.     

Studies on stability and dynamics of a swirling jet

The temporal instability and nonlinear evolution of the swirling jet near a nozzle exit are studied by both normal-mode method and three dimensional direct numerical simulation (3D DNS). It is found that the swirl enhances the maximum linear growth rates for negative helical modes, while decreases the growth rate for axisymmetric mode. Numerical simulations show that the evolution in early stage is compared well with the linear stability theory. In nonlinear stage, the swirl promotes the breakup of 3D large scale organized structures in the flow into small eddies. The project supported by the National Natural Science Foundation of China (19772052)     

Analysis of Momentum Transfer in a Lid-Driven Cavity Containing a Brinkman–Forchheimer Medium

The CE/SE (the space-time conservation element and solution element method) scheme with the second-order accuracy has been proposed. And the pretreatment method has been introduced to convert the parabolic equations to the hyperbolic equations, which are accurately solved by the CE/SE method. The lid-driven rectangular cavity containing a porous Brinkman–Forchheimer medium is studied in this numerical investigation. The Brinkman–Forchheimer equation is used such that both the inertial and viscous effects are incorporated. The governing equations are solved by the improved CE/SE approach. The characteristics of the flow are analyzed with emphasis on the influence of the Darcy number and the cavity depth. It is found that the porous medium effect decreases both the strength and the number of eddies, especially for deep cavities.     

Numerical Simulations of the Sandia Flame D Using the Eddy Dissipation Concept

A turbulent piloted methane/air diffusion flame (Sandia Flame D) is calculated using both compressible Reynolds-averaged and large-eddy simulations (RAS and LES, respectively). The Eddy Dissipation Concept (EDC) is used for the turbulence-chemistry interaction, which assumes that molecular mixing and the subsequent combustion occur in the fine structures (smaller dissipative eddies, which are close to the Kolmogorov length scales). Assuming the full turbulence energy cascade, the characteristic length and velocity scales of the fine structures are evaluated using a standard k- ?? turbulence model for RAS and a one-equation eddy-viscosity sub-grid scale model for LES. Finite-rate chemical kinetics are taken into account by treating the fine structures as constant pressure and adiabatic homogeneous reactors (calculated as a system of ordinary-differential equations (ODEs)) described by a Perfectly Stirred Reactor (PSR) concept. A robust implicit Runge-Kutta method (RADAU5) is used for integrating stiff ODEs to evaluate reaction rates. The radiation heat transfer is treated by the P1-approximation. The assumed β-PDF approach is applied to assess the influence of modeling of the turbulence-chemistry interaction. Numerical results are compared with available experimental data. In general, there is good agreement between present simulations and measurements both for RAS and LES, which gives a good indication on the adequacy and accuracy of the method and its further application for turbulent combustion simulations.     

提取壁湍流相干结构的数字滤波法

以三丝热线探头测得平板湍流边界层的数据为对象，提出用数字滤波技术将湍流信号分解为接近各向同性的小尺度涡和非各向同性的大尺度涡的方法。并用条件采样技术从大涡信号中提取相干结构。     

Dynamics of large-scale eddies in turbulent flows

A model equation based on the equipartition of the turbulent dissipation is proposed for describing the dynamics of large-scale eddies in turbulent flows. The equation is reducible to the equation of motion of an inviscid fluid, so that the motion of the large-scale eddies can be described in terms of inviscid fluid dynamics. It is found that the large-scale eddies are always weakened by the background turbulence and their evolution is slowed down compared with the corresponding inviscid motion. In the case of turbulent mixing layer, its linear growth in downstream direction is accounted for by the exponential growth in time of the perturbation in an inviscid plane vortex sheet.     

Multiscale simulations and experiments on water jet atomization

Numerical simulation of primary atomization at high Reynolds number is still a challenging problem. In this work a multiscale approach for the numerical simulation of liquid jet primary atomization is applied, using an Eulerian-Lagrangian coupling. In this approach, an Eulerian volume of fluid (VOF) method, where the Reynolds stresses are closed by a Reynolds stress model is applied to model the global spreading of the liquid jet. The formation of the micro-scale droplets, which are usually smaller than the grid spacing in the computational domain, is modelled by an energy-based sub-grid model. Where the disruptive forces (turbulence and surface pressure) of turbulent eddies near the surface of the jet overcome the capillary forces, droplets are released with the local properties of the corresponding eddies. The dynamics of the generated droplets are modelled using Lagrangian particle tracking (LPT). A numerical coupling between the Eulerian and Lagrangian frames is then established via source terms in conservation equations. As a follow-up study to our investigation in Saeedipour et al. (2016a), the present paper aims at modelling drop formation from liquid jets at high Reynolds numbers in the atomization regime and validating the simulation results against in-house experiments. For this purpose, phase-Doppler anemometry (PDA) was used to measure the droplet size and velocity distributions in sprays produced by water jet breakup at different Reynolds numbers in the atomization regime. The spray properties, such as droplet size spectra, local and global Sauter-mean drop sizes and velocity distributions obtained from the simulations are compared with experiment at various locations with very good agreement.     

KOLMOGOROV BEHAVIOR OF NEAR-WALL TURBULENCE AND ITS APPLICATION IN TURBULENCE MODELING

The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder¹ and followers^2,3,4,5. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterised by Kolmogorov microscales. According t o this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become “Kolmogorov” eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As a n example, the concept is incorporated in the standard κ - εmodel which is then applied t o channel and boundary layer flows. Using appropriate boundary conditions (based on Kolmogorov behaviour of near-wall turbulence), there is no need for any wall-modification to the κ - ε equations (including model constants). Results compare very well with the DNS and experimental data.     

Joint RANS/LES Approach to Premixed Flame Modelling in the Context of the TFC Combustion Model

We present an original timesaving joint RANS/LES approach to simulate turbulent premixed combustion. It is intended mainly for industrial applications where LES may not be practical. It is based on successive RANS/LES numerical modelling, where turbulent characteristics determined from RANS simulations are used in LES equations for estimation of the subgrid chemical source and viscosity. This approach has been developed using our TFC premixed combustion model, which is based on a generalization of the Kolmogorov’s ideas. We assume existence of small-scale statistically equilibrium structures not only of turbulence but also of the reaction zones. At the same time, non-equilibrium large-scale structures of reaction sheets and turbulent eddies are described statistically by model combustion and turbulence equations in RANS simulations or follow directly without modelling in LES. Assumption of small-scale equilibrium gives an opportunity to express the mean combustion rate (controlled by small-scale coupling of turbulence and chemistry) in the RANS and LES sub-problems in terms of integral or subgrid parameters of turbulence and the chemical time, i.e. the definition of the reaction rate is similar to that of the mean dissipation rate in turbulence models where it is expressed in terms of integral or subgrid turbulent parameters. Our approach therefore renders compatible the combustion and turbulent parts of the RANS and LES sub-problems and yields reasonable agreement between the RANS and averaged LES results. Combining RANS simulations of averaged fields with LES method (and especially coupled and acoustic codes) for simulation of corresponding nonstationary process (and unsteady combustion regimes) is a promising strategy for industrial applications. In this work we present results of simulations carried out employing the joint RANS/LES approach for three examples: High velocity premixed combustion in a channel, combustion in the shear flow behind an obstacle and the impinging flame (a premixed flame attached to an obstacle).     

A note on a new theory of turbulence

In this paper,the creation and annihilaton of turbulent eddies are described aselementary particles in the quantum field theory.An elementary particle may be consideredas a solid entity as it exists in quantum theory,but a turbulent eddy is often changed in sizeand shape with time due to its energy dissipation in a turbulent field.Therefore,in order toapply the method of the quantum field theory to the turbulent field by analogy,the entity ofthe same eddy should be defined firstly.According to the linearized theory,the turbulenteddies with the similarity character in lime duration may be considered as the entity of thesame eddy,and the creation and annihilation of turbulent eddies without the similarcharacters are related to the interaction term φ_2 in equation(2.6).Then,the creationoperator and annihilation operator similar to those in the quantum field theory are used todescribe the state of turbulent eddy field.Finally,a“Schrǒdinger”equation of turbulenteddies is formulated based upon the nonline