Evaluation of bubble-induced turbulence using direct numerical simulation

The presented research evaluates the interaction between a single bubble and homogeneous turbulent flow using direct numerical simulation (DNS) approach. The homogeneous single-phase turbulence is numerically generated by passing a uniform flow through grid planes. The turbulence decay rate is compared with experiment-based correlation. The single phase turbulence is then used as an inflow boundary condition for a set of single bubble studies. By estimating the turbulent field around the fully resolved bubble, the effects of bubble deformability, turbulent intensity and relative velocity on the bubble-induced turbulence are investigated. The existence of bubble creates new vortices in the wake region and the enhancement of turbulence is observed in the region behind the bubble. The results show that the magnitude of the turbulence enhancement would increase as the bubble encounters larger liquid turbulent intensity or higher relative velocity. Set of bubble Weber numbers from 0.34 to 3.39 are used to investigate the effect of bubble deformability. The more deformable bubble is the higher the increase in the magnitude of the turbulence enhancement behind the bubble. This research provides systematic insight on the bubble-induced turbulence (BIT) mechanism and is important for multiphase computational fluid dynamics (M-CFD) closure model development.     

Effects of sub-millimeter-bubble injection on transition to turbulence in natural convection boundary layer along a vertical plate in water

Temperature and velocity measurements are performed to clarify the effects of sub-millimeter-bubble injection on the transition to turbulence in the natural convection boundary layer along a vertical plate in water. In particular, we focus on the relationship between the bubble injection position L and the transition to turbulence in the natural convection boundary layer. The bubble injection positions used in our experiments are L = 1.6 and 3.6 mm. Bubble injection at L = 1.6 mm delays the transition to turbulence in the natural convection boundary layer, while that at L = 3.6 mm accelerates the transition to turbulence in the boundary layer. In the case of L = 1.6 mm, the appearance region of the liquid velocity fluctuation in the bubble-induced upward flow in the upstream unheated section is restricted to near the wall, although the peak of the liquid velocity fluctuation is high. In contrast, in the case of L = 3.6 mm, the relatively large liquid velocity fluctuation is distributed widely over the laminar boundary layer width. These results suggest that the effect of the liquid velocity fluctuation on the laminar boundary layer is quite different between L = 1.6 and 3.6 mm. It is therefore expected that the transition to turbulence in the natural convection boundary layer for the case with bubble injection is dependent on the magnitude and appearance region of the liquid velocity fluctuation in the bubble-induced upward flow in the upstream unheated section.     

A general model for unconfined bubble plumes from extended sources

A model for unconfined bubble plumes is developed without assumptions as to the form of the velocity and gas-fraction profiles or as to the fraction of the momentum flux associated with the mean flow. Dimensionless solutions for axisymmetric and plane-symmetric extended sources indicate an initial contraction followed by an almost linear expansion which closely resembles the single-phase case. A second contraction is predicted near the surface of deep bodies of liquid. Gas and liquid velocity measurements in laboratory-scale plumes are presented, providing information on entrainment coefficients, velocity profiles, bubble velocities and the contribution of velocity fluctutations to the total momentum flux. The latter effect is considerable, providing an explanation of the “too-low” plume velocities found by other investigators.     

Experimental investigation of horizontal air–water bubbly-to-plug and bubbly-to-slug transition flows in a 3.81 cm ID pipe

The present study seeks to investigate horizontal bubbly-to-plug and bubbly-to-slug transition flows. The two-phase flow structures and transition mechanisms in these transition flows are studied based on experimental database established using the local four-sensor conductivity probe in a 3.81 cm inner diameter pipe. While slug flow needs to be distinguished from plug flow due to the presence of large number of small bubbles (and thus, large interfacial area concentration), both differences and similarities are observed in the evolution of interfacial structures in bubbly-to-plug and bubbly-to-slug transitions. The bubbly-to-plug transition is studied by decreasing the liquid flow rate at a fixed gas flow rate. It is found that as the liquid flow rate is lowered, bubbles pack near the top wall of the pipe due to the diminished role of turbulent mixing. As the flow rate is lowered further, bubbles begin to coalesce and form the large bubbles characteristic of plug flow. Bubble size increases while bubble velocity decreases as liquid flow rate decreases, and the profile of the bubble velocity changes its shape due to the changing interfacial structure. The bubbly-to-slug transition is investigated by increasing the gas flow rate at a fixed liquid flow rate. In this transition, gas phase becomes more uniformly distributed throughout the cross-section due to the formation of large bubbles and the increasing bubble-induced turbulence. The size of small bubbles decreases while bubble velocity increases as gas flow rate increases. The distributions of bubble size and bubble velocity become more symmetric in this transition. While differences are observed in these two transitions, similarities are also noticed. As bubbly-to-plug or bubbly-to-slug transition occurs, the formation of large elongated bubbles is observed not in the uppermost region of bubble layer, but in a lower region. At the beginning of transitions, relative differences in phase velocities near the top of the pipe cross-section to those near the pipe center become larger for both gas and liquid phases, because more densely packed bubbles introduce more resistance to both phases.     

Experimental study on turbulent natural convection heat transfer in water with sub-millimeter-bubble injection

Using thermocouples and a particle tracking velocimetry technique, temperature and velocity measurements are conducted to investigate flow and heat transfer characteristics of turbulent natural convection from a vertical heated plate in water with sub-millimeter-bubble injection. Hydrogen-bubbles generated by the electrolysis of water are used as the sub-millimeter-bubbles. In the turbulent region, the heat transfer deterioration occurs for a bubble flow rate Q = 33 mm³/s, while the heat transfer enhancement occurs for Q = 56 mm³/s. Temperature and velocity measurements suggest that the former is caused by a delay of the transition due to the bubble-induced upward flow. On the other hand, the latter is mainly due to two factors: one is the enhancement of the rotation of eddies in the outer layer, and the other is the increase in the gradient of the streamwise liquid velocity at the heated wall. These are caused by bubbles, which are located in the inner layer, rising at high speed.     

An advanced LIF-PLV system for analysing the hydrodynamics in a laboratory bubble column at higher void fractions

Bubble columns are widely used in the chemical industry and biotechnology. Flow and turbulence in such an apparatus are induced by the bubble rise, and the bubble behaviour is strongly affected by swarm effects (i.e. the interaction between bubbles). For analysing the bubble swarm behaviour and simultaneously evaluating the flow structure and bubble-induced turbulence, a bubble column of 140 mm diameter and a height of 650 mm or 1,400 mm (initial water level) were considered. The bubble column was aerated with relatively fine bubbles having a mean size between about 0.5 and 4.0 mm. The gas hold-up was varied in the range between 0.5 and 19%. A two-phase pulsed-light velocimetry (PLV) system was developed to evaluate instantaneous flow fields of both rising bubbles and the continuous phase. The measurement of the liquid velocities in the bubble swarm was achieved by adding fluorescing seed particles. Images of bubbles and fluorescing tracer particles were acquired by two CCD cameras. Hence, the images from tracers and bubbles were easily separated by optical interference filters with a bandwidth corresponding to the emitting wavelength of the fluorescing tracer particles and the wavelength of the applied Nd-YAG pulsed laser, respectively. To improve the phase separation of the system, the CCD cameras were additionally placed in a non-perpendicular arrangement with respect to the light sheet. The acquired images were evaluated with the minimum-quadratic-difference algorithm. The potential of this technique for the analysis of bubbly flows with higher void fraction was explored. In order to obtain averaged velocity maps of bubble and fluid within the entire column, about 1,000 image pairs were recorded and evaluated for each phase. In addition, turbulence intensities of the fluid were deduced from the measurements. The turbulence properties were used to characterise bubble-induced turbulence for various bubble mean diameters and gas hold-ups. Moreover, the determination of the average bubble slip velocity within the bubble swarm was possible.     

Dynamic responses of a solid wall in contact with a bubbly liquid excited by thermal shock loading

Several MW-class spallation neutron sources are being developed in the world. Specifically, intensive and high energy protons are injected into heavy liquid metals (mercury, lead or lead-bismuth eutectic) to induce the spallation reaction that produces neutrons. At the moment when the proton beams are injected, thermal shock occurs in the liquid metal, causing pressure waves to propagate in the liquid metal, collide against the container and damage it.It is proposed that microbubbles are injected into the liquid metal to mitigate the impulsive pressure waves by means of absorption and attenuation effects. These effects are dependent on the relationship between bubble size and the rate of pressure increase. In the present experiment, a very rapid rise in pressure in the order of MPa/μs, equivalent to the rise in pressure due to proton beam injection, was simulated by the electric discharge method in a water loop test to investigate the impulsive pressure mitigation effect of injected microbubbles. The solid wall response was measured using an accelerometer, and the dynamic responses of microbubbles were observed using an ultra-high-speed camera filming at 5 × 10⁵ frame/s. The sound velocity in bubbly water was estimated using a differential image technique. It was confirmed from the experimental results that microbubbles are effective in reducing impulsive pressure waves and to suppressing the impact vibration of the solid wall in contact with the liquid.     

Hydrodynamic and mass transfer characteristics of slug flow in a vertical pipe with and without dispersed small bubbles

In this work, we present a numerical study to investigate the hydrodynamic characteristics of slug flow and the mechanism of slug flow induced CO₂ corrosion with and without dispersed small bubbles. The simulations are performed using the coupled model put forward by the authors in previous paper, which can deal with the multiphase flow with the gas–liquid interfaces of different length scales. A quasi slug flow, where two hypotheses are imposed, is built to approximate real slug flow. In the region ahead of the Taylor bubble and the liquid film region, the presence of dispersed small bubbles has less impacts on velocity field, because there are no non-regular intensive disturbance forces or centrifugal forces breaking the balance of the liquid and the dispersed small bubbles. In the liquid slug region, the strong centrifugal forces generated by the recirculation below the Taylor bubble lead to the effect of heterogeneity, which makes the profile of the radial liquid velocity component sharper with higher volume fraction of dispersed small bubbles. The volume fraction has a maximum value in the range of r/R = 0.5–0.6. Meanwhile, it is usually higher than 0.35, which means that larger dispersed bubbles can be formed by coalescences in this region. These calculated results are in good agreement with experimental results. The wall shear stress and the mass transfer coefficient with dispersed small bubbles are higher than those without dispersed small bubbles due to enhanced fluctuations. For short Taylor bubble length, the average mass transfer coefficient is increased when the gas or liquid superficial velocity is increased. However, there may be an inflection point at low mixture superficial velocities. For the slug with dispersed small bubbles, the product scales still cannot be damaged directly despite higher wall shear stress. In fact, the alternate wall shear stress and the pressure fluctuations perpendicular to the pipe wall with high frequency are the main cause for breaking the product scales.     

Effects of bubble size on heat transfer enhancement by sub-millimeter bubbles for laminar natural convection along a vertical plate

Injection of sub-millimeter bubbles is considered a promising technique for enhancing natural convection heat transfer for liquids. So far, we have experimentally investigated heat transfer characteristics of laminar natural convection flows with sub-millimeter bubbles. However, the effects of the bubble size on the heat transfer have not yet been understood. The purpose of this study is to clarify the effects of the bubble size on the heat transfer enhancement for the laminar natural convection of water along a vertical heated plate with uniform heat flux. Temperature and velocity measurements, in which thermocouples and a particle tracking velocimetry technique are, respectively used, are conducted to investigate heat transfer and flow characteristics for different bubble sizes. Moreover, two-dimensional numerical simulations are performed to comprehensively understand the effects of bubble injection on the flow near the heated plate. The result shows that the ratio of the heat transfer coefficient with sub-millimeter-bubble injection to that without injection ranges from 1.3 to 2.2. The result also shows that for a constant bubble flow rate, the heat transfer coefficient ratio increases with a decrease in the mean bubble diameter. It is expected from our estimation based on both experimental data and simulation results that this increase results from an increase in the advection effect due to bubbles.     

Gas–liquid mass transfer in the gas–liquid–solid mini fluidized beds

The gas–liquid–solid mini fluidized bed (GLSMFB) combines the advantages of fluidized bed and micro-reactor, and meets the requirements for safety and efficiency of green development of process industry. However, there are few studies on its flow performance and no studies on its mass and heat transfer performance. In this paper, the characteristics of gas–liquid mass transfer in a GLSMFB were studied in order to provide basic guidance for the study of GLSMFB reaction performance and application. Using CO₂ absorption by NaOH as the model process, the gas–liquid mass transfer performance of GLSMFB was investigated. The results show that the liquid volumetric mass transfer coefficient and the gas–liquid interfacial area both increase with the increase of the superficial gas velocity within the experimental parameter range under the same given superficial liquid velocity. At the same ratio of superficial gas to liquid velocity, the liquid volumetric mass transfer coefficient increases with the increase of the superficial liquid velocity. Fluidized solid particles strengthen the liquid mass transfer process, and the liquid volumetric mass transfer coefficient is about 13% higher than that of gas–liquid mini bubble column.     

Two-phase flow characteristics in gas–liquid microreactors

Multiphase chemical microreactors require a detailed knowledge of the flow conditions inside the reaction system. This paper reports flow visualization measurements of the two-phase gas–liquid flow pattern and the liquid velocity distribution inside liquid plugs of an intermittent flow. Rectangular cross-section silicon microchannels with hydraulic diameters between 187.5 and 218 μm are fabricated. Laser Induced Fluorescence (LIF) is used to determine the flow pattern. To analyze the influence of the liquid properties and the channel diameter on the two-phase flow pattern, we present flow regime maps using different channel geometries and fluids. A universal flow pattern map based on dimensional analysis is presented. In contrast to microchannel flows, a great number of correlations for flow characteristics for multiphase flow in (round) pipes with diameters >1 mm exist. We compare our experimental results from optical flow visualizations in microreactors with common flow correlations and regime maps for macro- and microchannels. The recirculation motion in the liquid segments of an intermittent gas–liquid flow is analyzed using micron-resolution particle image velocimetry (μPIV). The velocity distribution influences the mixing and the mass transport towards the reactive phase interface dealing with two-phase chemical reactions. For straight microchannels hardly any mass transport over the center line is quantified. For enhanced mixing geometrical adaptations are suggested.     

Macro-to-microchannel transition in two-phase flow: Part 2 - Flow boiling heat transfer and critical heat flux

This part of the paper presents the current experimental flow boiling heat transfer and CHF data acquired for R134a, R236fa and R245fa in single, horizontal channels of 1.03, 2.20 and 3.04 mm diameters over a range of experimental conditions. The aim of this study is to investigate the effects of channel confinement, heat flux, flow pattern, saturation temperature, subcooling and working fluid properties on the two-phase heat transfer and CHF. Experimentally, it was observed that the flow boiling heat transfer coefficients are a significant function of the type of two-phase flow pattern. Furthermore, the monotonically increasing heat transfer coefficients at higher vapor qualities, corresponding to annular flow, signifies convective boiling as the dominant heat transfer mechanism in these small scale channels. The decreasing heat transfer trend at low vapor qualities in the slug flow (coalescing bubble dominated regime) was indicative of thin film evaporation with intermittent dry patch formation and rewetting at these conditions. The coalescing bubble flow heat transfer data were well predicted by the three-zone model when setting the dryout thickness to the measured surface roughness, indicating for the first time a roughness effect on the flow boiling heat transfer coefficient in this regime. The CHF data acquired during the experimental campaign indicated the influence of saturation temperature, mass velocity, channel confinement and fluid properties on CHF but no influence of inlet subcooling for the conditions tested. When globally comparing the CHF values for R134a in the 0.51-3.04 mm diameter channels, a peak in CHF peak was observed lying in between the 0.79 (Co ≈ 0.99) and 1.03 (Co ≈ 0.78) mm channels. A new CHF correlation has been proposed involving the confinement number, Co that is able to predict CHF for R134a, R236fa and R245fa in single-circular channels, rectangular multichannels and split flow rectangular multichannels. In summary, the present flow boiling and CHF trends point to a macro-to-microscale transition as indicated by the results presented in Ong and Thome (2011) [1].     

Interaction between Taylor bubbles rising in stagnant non-Newtonian fluids

The interaction between Taylor bubbles rising in stagnant non-Newtonian solutions was studied. Aqueous solutions of carboxymethylcellulose (CMC) and polyacrylamide (PAA) polymers were used to study the effect of different rheological properties: shear viscosity and viscoelasticity. The solutions studied covered a range of Reynolds numbers between 10 and 714, and Deborah numbers up to 14. The study was performed with pairs of Taylor bubbles rising in a vertical column (0.032 m internal diameter) filled with stagnant liquid. The velocities of the leading and trailing bubbles were measured by sets of laser diodes/photocells placed along the column. The velocity of the trailing bubble was analysed together with the liquid velocity profile in the wake of a single rising bubble (Particle Image Velocimetry data obtained from the literature). For the less concentrated CMC solutions, with moderate shear viscosity and low viscoelasticity, the interaction between Taylor bubbles was similar to that found in Newtonian fluids. For the most concentrated CMC solution, which has high shear viscosity and moderate viscoelasticity, a negative wake forms behind the Taylor bubbles, inhibiting coalescence since the bubbles maintain a minimum distance of about 1D between them. For the PAA solutions, with moderate shear viscosity but higher viscoelasticity than the CMC solutions, longer wake lengths are seen, which are responsible for trailing bubble acceleration at greater distances from the leading bubble. Also in the PAA solutions, the long time needed for the fluid to recover its initial shear viscosity after the passage of the first bubble makes the fluid less resistant to the trailing bubble flow. Hence, the trailing bubble can travel at a higher velocity than the leading bubble, even at distances above 90D.     

An investigation on near wall transport characteristics in an adiabatic upward gas–liquid two-phase slug flow

The near-wall transport characteristics, inclusive of mass transfer coefficient and wall shear stress, which have a great effect on gas–liquid two-phase flow induced internal corrosion of low alloy pipelines in vertical upward oil and gas mixing transport, have been both mechanistically and experimentally investigated in this paper. Based on the analyses on the hydrodynamic characteristics of an upward slug unit, the mass transfer in the near wall can be divided into four zones, Taylor bubble nose zone, falling liquid film zone, Taylor bubble wake zone and the remaining liquid slug zone; the wall shear stress can be divided into two zones, the positive wall shear stress zone associated with the falling liquid film and the negative wall shear stress zone associated with the liquid slug. Based on the conventional mass transfer and wall shear stress characteristics formulas of single phase liquid full-pipe turbulent flow, corrected normalized mass transfer coefficient formula and wall shear stress formula are proposed. The calculated results are in good agreement with the experimental data. The shear stress and the mass transfer coefficient in the near wall zone are increased with the increase of superficial gas velocity and decreased with the increase of superficial liquid velocity. The mass transfer coefficients in the falling liquid film zone and the wake zone of leading Taylor bubble are lager than those in the Taylor bubble nose zone and the remaining liquid slug zone, and the wall shear stress associated falling liquid film is larger than that associated the liquid slug. The mass transfer coefficient is within 10⁻³ m/s, and the wall shear stress below 10³ Pa. It can be concluded that the alternate wall shear stress due to upward gas–liquid slug flow is considered to be the major cause of the corrosion production film fatigue cracking.     

Experimental investigation of vapor bubble growth during flow boiling in a microchannel

Experiments were conducted to analyze flow boiling characteristics of water in a single brass microchannel of 25 mm length, 201 μm width, and 266 μm depth. Different heat flux conditions were tested for each of two different mass flow rates over three different values of inlet fluid temperature. Temporal and spatial surface temperature profiles were analyzed to show the relative effect of axial heat conduction on temperature rise along the channel length and the effect of flow regime transition on local surface temperature oscillation. Vapor bubble growth rate increased with increasing wall superheat. The slower a bubble grew, the further it was carried downstream by the moving liquid. Bubble growth was suppressed for increased mass flux while the vapor bubble was less than the channel diameter. The pressure spike of an elongating vapor bubble was shown to suppress the growth of a neighboring bubble by more than 50% of its volume. An upstream progression of the Onset of Bubble Elongation (OBE) was observed that began at the channel exit and progressed upstream. The effects of conjugate heat transfer were observed when different flow regime transitions produced different rates of progression for the elongation sequence. Instability was observed at lower heat fluxes for this single channel experiment than for similar studies with multiple channels.     

Heat transfer to air–water annular flow in a horizontal pipe

Two-phase air–water flow and heat transfer in a 25 mm internal diameter horizontal pipe were investigated experimentally. The water superficial velocity varied from 24.2 m/s to 41.5 m/s and the air superficial velocity varied from 0.02 m/s to 0.09 m/s. The aim of the study was to determine the heat transfer coefficient and its connection to flow pattern and liquid film thickness. The flow patterns were visualized using a high speed video camera, and the film thickness was measured by the conductive tomography technique. The heat transfer coefficient was calculated from the temperature measurements using the infrared thermography method. It was found that the heat transfer coefficient at the bottom of the pipe is up to three times higher than that at the top, and becomes more uniform around the pipe for higher air flow-rates. Correlations on local and average Nusselt number were obtained and compared to results reported in the literature. The behavior of local heat transfer coefficient was analyzed and the role of film thickness and flow pattern was clarified.     

Flow regime development analysis in adiabatic upward two-phase flow in a vertical annulus

In this work radial and axial flow regime development in adiabatic upward air-water two-phase flow in a vertical annulus has been investigated. Local flow regimes have been identified using conductivity probes and neural networks techniques. The inner and outer diameters of the annulus are 19.1 mm and 38.1 mm, respectively. The equivalent hydraulic diameter of the flow channel, D_H, is 19.0 mm and the total length is 4.37 m. The flow regime map includes 1080 local flow regimes identifications in 72 flow conditions within a range of 0.01 m/s < 〈j_g〉 < 30 m/s and 0.2 m/s < 〈j_f〉 < 3.5 m/s where 〈j_g〉 and 〈j_f〉 are, respectively, superficial gas and liquid velocities. The local flow regime has been classified into four categories: bubbly, cap-slug, churn-turbulent and annular flows. In order to study the radial and axial development of flow regime the measurements have been performed at five radial locations. The three axial positions correspond to z/D_H = 52, 149 and 230, where z represents the axial position. The flow regime indicator has been chosen as some statistical parameters of local bubble chord length distributions and self-organized neural networks have been used as mapping system. This information has been also used to compare the results given by the existing flow regime transition models. The local flow regime is characterized basically by the void fraction and bubble chord length. The radial development of flow regime shows partial and complete local flow regime combinations. The radial development is controlled by axial location and superficial liquid velocity. The radial flow regime transition is always initiated in the center of the flow channel and it is propagated towards the channel boundaries. The axial development of flow regime is observed in all the flow maps and it is governed by superficial liquid velocity and radial location. The prediction results of the models are compared for each flow regime transition.     

An updated three-zone heat transfer model for slug flow boiling in microchannels

This work proposes a novel physics-based model for the fluid mechanics and heat transfer associated with slug flow boiling in horizontal circular microchannels to update the widely used three-zone model of Thome et al. (2004). The heat transfer model has a convective boiling nature and predicts the time-dependent variation of the local heat transfer coefficient during the cyclic passage of a liquid slug, an evaporating elongated bubble and a vapor plug. The capillary flow theory, extended to incorporate evaporation effects, is applied to estimate the bubble velocity along the channel. A liquid film thickness prediction method also considering bubble proximity effects, which may limit the radial extension of the film, is included. The minimum liquid film thickness at dryout is set to the channel wall roughness. Theoretical heat transfer models accounting for the thermal inertia of the liquid film and for the recirculating flow within the liquid slug are utilized. The heat transfer model is compared to experimental data taken from three independent studies. The 833 slug flow boiling data points cover the fluids R134a, R245fa and R236fa, and channel diameters below 1 mm. The proposed evaporation model predicts more than 80% of the database to within ±30%. It demonstrates a stronger contribution to heat transfer by the liquid slugs and correspondingly less by the thin film evaporation process compared to the original three-zone model. This model represents a new step towards a complete physics-based modelling of the bubble dynamics and heat transfer within microchannels under evaporating flow conditions.     

煤油液滴直径对两相旋转爆轰发动机流场的影响

为探究煤油液滴不同初始直径对气液两相旋转爆轰发动机流场的影响,假设初始注入的煤油液滴具有均匀直径,考虑雾化破碎、蒸发等过程,建立了非定常两相爆轰的Eulerian-Lagrangian模型,进行了液态煤油/高温空气爆轰的非预混二维数值模拟。结果表明：在初始液滴直径为1～70μm的工况范围,燃烧室内均形成了单个稳定传播的旋转爆轰波;全局当量比为1时,爆轰波前的空气区域大于液滴煤油的蒸气区域,导致波前燃料空气混合不均匀,波前均存在富油区和贫油区,两相速度差导致分离出的空气形成低温条带;当煤油液滴的初始直径较小时,波前的反应物混合过程主要受蒸发的影响,爆轰波可稳定传播;当直径减小至1μm时,煤油液滴在入口处即蒸发,旋转爆轰波表现为气相传播的特性,爆轰波结构平整;当煤油液滴的初始直径较大时,波前的反应物混合过程主要受液滴破碎的影响;对于相同的燃料质量流量,在不同初始煤油液滴直径工况下,煤油液滴最大的停留时间均占爆轰波传播时间尺度的80%以上;爆轰波前燃料预蒸发为气相的占比越高,爆轰波的传播速度越高;初始液滴直径为10～70μm的工况范围内,爆轰波的速度随初始直径的增大先升高后降低。     

Flow structure of microbubble-laden turbulent channel flow measured by PIV combined with the shadow image technique

The turbulence structure of a horizontal channel flow with microbubbles is experimentally investigated using combined particle image velocimetry (PIV) in order to clarify the mechanism of drag reduction caused by microbubbles. A new system which simultaneously measures the liquid phase and the dispersed bubbles is proposed, based on a combination of particle tracking velocimetry (PTV), laser-induced fluorescence (LIF) and the shadow image technique (SIT). To accurately obtain the velocity of the liquid phase, tracer particles which overlap with the bubble shadow images are almost entirely eliminated in the post-processing. Finally, the turbulence characteristics of the flow field are presented, including measurements for both phases, and the bubble effect on the turbulence is quantified.