可压缩气固混合层中离散相与连续相的相互作用研究

尽管已有许多文献采用数值模拟方法研究两相流问题,但主要是集中不可压流动方面.本文采用Eul-er-Lagrange颗粒-轨道双向耦合模型对时间模式下含有固粒的二维可压缩混合层流动进行了研究.气相流场采用非定常全Navier-Stokes方程描述,并应用具有空间三阶精度的WNND(Weighted Non-Oscillatory, Contai-ning No Free Parameters and Dissipative)格式进行数值高散.固相方程采用二阶单边三点差分离散.在考虑流场对固粒作用的同时,也计及颗粒对流场的反作用.主要研究混合层大尺度涡对颗粒扩散特性的影响及颗粒对流场结构的影响问题.在对流马赫数为0.5时,研究不同Stokes数颗粒在连续流场中的扩散特性,而在对流马赫数为0.8时研究了不同Stokes数颗粒对流场小激波结构的影响.     

The structure and dynamics of reacting plane mixing layers

The reacting two-dimensional plane mixing layer has been studied in two configurations: a rearward facing step and a two-stream mixing layer. Observations have been made of the steady state behavior, and the unsteady behavior when the flow was forced by a specific acoustic frequency. The steady behavior of the mean properties of the reacting flows is similar to that of non-reacting free shear flows except for the global effects of thermodynamic property changes. The structure of these flows is qualitatively similar to that of non-reacting flows. Vortices form by the two-dimensional Kelvin-Helmholtz instability and grow by subharmonic combination until the mixing layer interacts with the walls. Entrainment is dominated by the two-dimensional vortex motion. Three-dimensional instabilities give rise to secondary vortices which are coherent over several Kelvin-Helmholtz structures and dominate the fine scale mixing process. The mixing transition corresponds to a loss of coherence in the layer. Unsteady behavior occurs when there are resonant interactions with the Kelvin-Helmholtz instability or the instability associated with the recirculation vortex in the rearward facing step flow. Modeling efforts are reported which show promise of simulating the essential features of plane mixing layers.A version of this paper was presented at the ASME Winter Annual Meeting of 1984 and printed in AMD-Vol. 66     

DEM investigation of mixing indices in a ribbon mixer

Mixing index is an important parameter to understand and assess the mixing state in various mixers including ribbon mixers, the typical food processing devices. Many mixing indices based on either sample variance methods or non-sample variance methods have been proposed and used in the past, however, they were not well compared in the literature to evaluate their accuracy of assessing the final mixing state. In this study, discrete element method (DEM) modelling is used to investigate and compare the accuracy of these mixing indices for mixing of uniform particles in a horizontal cylindrical ribbon mixer. The sample variance methods for mixing indices are first compared both at particle- and macro-scale levels. In addition, non-sample variance methods, namely entropy and non-sampling indices are compared against the results from the sample variance methods. The simulation results indicate that, among the indices considered in this study, Lacey index shows the most accurate results. The Lacey index is regarded to be the most suitable mixing index to evaluate the steady-state mixing state of the ribbon mixer in the real-time (or without stopping the impeller) at both the particle- and macro-scale levels. The study is useful for the selection of a proper mixing index for a specific mixture in a given mixer.     

Quantifying macro-mixing and micro-mixing in a static mixer using two-tracer laser-induced fluorescence

An experimental procedure has been developed to quantify mixing at large scales (flow-induced) and at small scales (induced by molecular diffusion). It relies on the simultaneous imaging of two different fluorescent tracers using planar laser-induced fluorescence (PLIF). In order to quantify micro-mixing, a suitable neutralization reaction involving the fluorescent tracer uranine has been identified. Using PLIF, uranine is measured simultaneously with another fluorescent tracer, pyridine 2, employed to characterize macro-mixing. Since both tracers are quite inexpensive, this procedure allows an in-depth characterization of mixing properties even in large installations, by measuring the concentration fields of the involved tracers in a non-intrusive manner. This measurement procedure has been applied to a static mixer segment with geometrical features and dimensions similar to that found in practical applications. Laminar inflow conditions are employed. The flow and mixing analysis obtained by post-processing the measurement results is detailed in the present article.     

Direct numerical simulation of turbulent mixing of a passive scalar in pipe flow

Turbulent mixing of a passive scalar in fully developed turbulent pipe flow has been investigated by means of a Direct Numerical Simulation (DNS). The scalar is released from a point source located on the centreline of the pipe. The domain size of the concentration field has been chosen large enough to capture the different stages of turbulent mixing. Results are presented for mean concentration profiles, turbulent fluxes, concentration fluctuations, probability density functions and higher-order moments. To validate the numerical simulations the results are compared with experimental data on mixing in grid-turbulence that have been reported in the literature. The agreement between the experimental measurements and the computations is satisfactory. We have also considered the Probability Density Function (PDF). For small diffusion times and positions not on the plume centreline, our results lead to a PDF of an exponential form with a large peak at zero concentration. When the diffusion time increases, the PDF shifts from a exponential to a more Gaussian form.     

Large eddy simulation of hot and cold fluids mixing in a T-junction for predicting thermal fluctuations

Temperature fluctuations in a mixing T-junction have been simulated on the FLUENT platform using the large eddy simulation （LES） turbulent flow model and a sub-grid scale Smagorinsky-Lilly model. The normalized mean and root mean square temperatures for describing time-averaged temperature and temperature fluctuation intensity, and the velocity are obtained. The power spectrum densities of temperature fluctuations, which are key parameters for thermal fatigue analysis and lifetime evaluation, are analyzed. Simulation results are consistent with experimental data published in the literature, showing that the LES is reliable. Several mixing processes under different conditions are simulated in order to analyze the effects of varying Reynolds number and Richardson number on the mixing course and thermal fluctuations.     

Color schlieren deflectometry study of jet mixing: effect of buoyancy and perforation

Mixing of jets is crucial for optimal performance of many industrial applications and there is a need to optimize both nozzle geometry and flow conditions. The present study reports the influence of buoyancy and perforation on mixing between a jet and its environment. Optical techniques are ideal for the study of jet mixing due to their non-intrusive and inertia free properties. The present study gives an account of mixing between helium jet and the ambient fluid using a combination of color schlieren deflectometry and radial tomographic mathematics. Four different perforation sizes have been used and the experiments are performed for Reynolds numbers 21–676 and Richardson numbers 3.27–0.0015. Color schlieren images show distinct influence of perforation and flow conditions (Richardson number). Oxygen concentration and jet width quantify effectiveness of jet mixing. Buoyancy plays an important role in mixing at high Richardson number. Perforation improves jet mixing i.e. there is about 120% increase in jet width and the size of perforation plays an important role.     

Numerical simulation on micromixer based on synthetic jet

This paper studied a concept of micromixer with a synthetic jet placed at the bottom of a rectangular channel. Due to periodic ejections from and suctions into the channel, the fluids are mixed effectively. To study the effects of the inlet velocity, the jet intensity and frequency, and the jet location on the mixing efficiency, 3-D numerical simulations of the micromixer have been carried out. It has been found that when the jet intensity and the frequency are fixed, the mixing efficiency increases when Re 〈 50, and decreases when Re 〉 50 with the best mixing efficiency achieved at Re = 50. When the ratio of the jet velocity magnitude to the inlet velocity is taken as 10 and the jet frequency is 100 Hz, the mixing index reaches the highest value. It has also been found that to get better mixing efficiency, the orifice of the synthetic jet should be asymmetrically located away from the channel＇s centerline.     

Reducing mixing at the outlet bore of a cylinder

Summary In some technical applications mixing caused by contraction of cross section is to be avoided. An example is the outlet of a chromatographic column in preparative liquid chromatography (HPLC). First of all, it is desirable to reduce the volume of the cylinder within which the mixing occurs. This is done by placing a porous disk upstream of the contraction. Mixing within the contraction (or funnel) is avoided if a constant axial velocity component over every cross-section is achieved. Two such velocity fields have been found analytically. The funnel shapes have been taken to be the stream tubes of these fields. The mixing within these funnels, and indeed in other common funnels, is described by the breakthrough curves or residence-time distributions, which are numerically determined. An experiment has been set up to measure the breakthrough curves for different funnels. Received 6 September 1997; accepted for publication 26 January 1998     

Promotion of Turbulent Thermal Mixing of Hot and Cold Airflows in T-junction

An experimental study has been conducted on the promotion and control of turbulent thermal mixing of hot and cold airflows in a T-junction with rectangular cross sections, which simulates the HVAC unit for automobile air-conditioning system. In order to promote the turbulent thermal mixing, small jets have been blown into the main channel at the upstream edge of the T-junction in the direction of 45° against the main flow. Turbulence intensity in the upper part of the thermal mixing layer can be increased with these jets, and consequently the turbulent mixing of hot and cold airflows is promoted effectively. Moreover, it has been found that the degree of thermal mixing can be controlled by changing the jet velocity.     

Numerical Simulation of Particle Dispersion in a Spatially Developing Mixing Layer

Although there have been several numerical studies on particle dispersion in mixing layers, most of them have been conducted for temporally evolving mixing layers. In this study, numerical simulations of a spatially developing mixing layer are performed to investigate particle dispersion under various conditions. The full compressible Navier--Stokes equations are solved with a high-order compact finite difference scheme, along with high-order time-integration. Accurate non-reflecting boundary conditions for the fluid flow are used, and several methods for introducing particles into the computational domain are tested. The particles are traced using a Lagrangian approach assuming one-way coupling between the continuous and the dispersed phases. The study focuses on the roles of the large-scale vortex structures in particle dispersion at low, medium and high Stokes numbers, which highlights the important effects of interacting vortex structures in nearby regions in the spatially developing mixing layer. The effects of particles with randomly distributed sizes (or Stokes numbers) are also investigated. Both instantaneous flow fields and statistical quantities are analyzed, which reveals essential features of particle dispersion in spatially developing free shear flows, which are different from those observed in temporally developing flows. The inclusion of the gravity not only modifies the overall dispersion patterns, but also enhances stream-crossing by particles. Received 7 June 2001 and accepted 19 February 2002     

Modeling of Liquid Fuel Injection, Evaporation and Mixing in a Gas Turbine Burner Using Large Eddy Simulations

The interaction of turbulence, temperature fluctuation, liquid fuel transport, mixing and evaporation is studied by using Large Eddy Simulations (LES). To assess the accuracy of the different components of the methods we consider first isothermal, single phase flow in a straight duct. The results using different numerical methods incorporating dynamic Sub-Grid-Scale (SGS) models are compared with DNS and experimental data. The effects of the interactions among turbulence, temperature fluctuation, spray transport, evaporation and mixing of the gaseous fuel are studied by using different assumptions on the temperature field. It has been found that there are strong non-linear interactions among temperature-fluctuation, evaporation and turbulent mixing which require additional modeling if not full LES is used. The developed models and methods have been applied to a gas turbine burner into which liquid fuel is injected. The dispersion of the droplets in the burner is described. This revised version was published online in July 2006 with corrections to the Cover Date.     

Shock induced Richtmyer-Meshkov instability in the presence of a wall boundary layer

An experimental investigation on gaseous mixing zones originated from the Richtmyer-Meshkov instability has been undertaken in a square cross section shock tube. Mass concentration fields, of one of the two mixing constituents, have been determined within the mixing zone when the shock wave passes from the heavy gas to the light one, from one gas to an other of close density, and from the light gas to the heavy one. Results have been obtained before and after the coming back of the reflected shock wave. The diagnostic method is based on the infrared absorption of one of the two constituents of the mixing zone. It is shown that the mixing zone is strongly deformed by the wall boundary layer. The consequence is the presence of strong gradients of concentration in the direction perpendicular to the shock wave propagation. Finally, it is pointed out that the mixing goes more homogeneous when the Atwood number tends to zero.     

Simultaneous PIV and concentration measurements in a gas-turbine combustor model

Simultaneous velocity and concentration measurements have been performed in a gas-turbine combustor model. Particle image velocimetry (PIV) was used to acquire planar velocity information and to identify coherent flow structures. The Mie scattering technique, based on a slightly modified experimental setup, was used for concentration measurements in this mixing flow. The degree of mixing was assessed by examining local concentration measurements while inhomogeneously seeding the primary and secondary stream of the mixing layer. Connections between flow field and concentration distribution were highlighted using the proper orthogonal decomposition algorithm (POD). Uncertainties and systematic errors for the PIV measurements due to the suboptimal seeding are discussed using a comparison with a second test series at optimal seeding conditions. Results are presented for several flow parameters and at various lateral planes.     

Shock-induced mixing zone characterization: an attempt by hot-wire diagnostic

An attempt to extract velocity and molar fraction from a single hot-wire trace within a turbulent mixing zone induced by a shock accelerated gaseous interface has been proposed. Experiments have been conducted for negative and positive density jumps across the interface. The hot-wire signals clearly show interfaces between mixed and unmixed regions and the locations of incident and reflected shocks. With some hypotheses on the temperature, velocity and molar fraction profiles within the turbulent mixing zone have been obtained solving an inverse problem. Results show that if the molar fraction profiles follow physically coherent evolutions, those of the local velocity are strongly correlated with the choice of its variation range. So, we reasonably think that the results obtained from single wire have to remain limited to interface and shock locations. And it is only by coupling the present technique with the laser Doppler velocimetry, which we will be able to possibly obtain reliable estimates of the variations of quantities in the turbulent mixing zone.      

Passive and reactive scalar measurements in a transient high-Schmidt-number Rayleigh–Taylor mixing layer

An experimental study of mixing induced by Rayleigh?CTaylor (RT) instability at an Atwood number (A _t ) ~7.5?×?10^?4 and Schmidt number (Sc) ~1,000 has been performed. A new transient experimental facility developed on the working principles of the draw-tank facility at Cambridge (Dalziel et al. in J Fluid Mech, 399:1?C48, 1999) has been established and enhanced to observe a higher (2×) Reynolds number regime. Water and brine were used to produce the RT density stratification. The evolution of the instability was studied using passive and reactive scalar techniques and quantified using optical diagnostic methods. The data were combined to estimate local and global mixing metrics representative of the mixing mechanism across the mixing layer. In comparison with parameters reported from analogous experiments, the mixing phenomenon at a high Sc shows a strong dependency on the initial conditions prevailing at the onset of the instability and the evidence of a delay in the mixing transition. Values of global and integral mixing parameters did not reach late-time asymptotic values that have been reported previously from steady-state experiments (Texas A&M Water Channel) and may be attributed to the effect of the barrier pull and the overturning mechanism that is thought to hinder the progress of the mixing layer.     

Accuracy and resolution issues in NO/acetone PLIF measurements of gas-phase molecular mixing

Under-resolved passive-scalar measurements in liquid and gaseous shear layers have been known to over-predict mixed-fluid quantities as a result of sub-resolution stirring. In recent studies, simultaneous cold-chemistry (nitric oxide) and passive-scalar (acetone) planar laser-induced fluorescence have been employed for the measurement of sub-resolution molecular mixing. In the current work, the theory and experimental approach of this technique are reviewed, and new analyses of quenching corrections, differential diffusion errors, and imaging resolution are presented. Results from both experimental work and direct numerical simulations indicate that this technique provides a successful means of quantifying macro- and micro-scale turbulent molecular mixing in gaseous shear flows. Received: 20 December 2000/Accepted: 24 September 2001     

Determination of mixing time by colourimetric diagnosis—Application to a new mixing system

A method for determining mixing times in a transparent stirred vessel for highly viscous fluids has been developed. The method involves a chemical decolourization with two acid-base indicators and a colourimetric diagnosis of digital images, deduced from a video captured during the mixing process. The advantages of the method developed are:(1) there is no disturbance of the fluid flow, and (2) measuring and determination of the degree of homogeneity is relatively effortless and cheap since it requires only a standard colour camera and commercial software for digital image analysis. This paper describes protocols for carrying this method out. In addition, examples of experimental results obtained for a standard impeller (of the helical ribbon agitator type) and a new piece of mixing equipment are given and discussed.