二维平板间驻波声流的MRT-LBM数值研究

采用多松弛时间格子玻尔兹曼方法（Multiple Relaxation Time Lattice Boltzmann Method,MRT-LBM）对二维平板间的驻波声流进行数值模拟,模拟结果与Rayleigh流近似解析解相符,研究黏度和板间宽度对驻波声流的影响,得到不同黏度下x=L/4截面无量纲水平速度分布和x=L/2截面无量纲竖直速度分布,板间宽度对边界层内声流区域厚度的影响及驻波声流的形成过程,结果表明MRT-LBM模型能有效模拟驻波声流效应.     

Acoustic Streaming under Thermal Boundary Conditions of the Third Kind

A numerical study of acoustic streaming in a cylindrical cavity subjected to vibrational action with a small vibration amplitude is performed. The dependence of the streaming character on the intensity of heat exchange with the surrounding medium is studied. A change in the forms of the acoustic streaming vortices is shown for a smooth transition from adiabatic to isothermal boundary conditions, which occurs via variation in the heat transfer coefficient. The values of the dimensionless heat transfer coefficient are determined for which the acoustic streaming pattern is close to the limiting cases corresponding to adiabatic and isothermal boundary conditions. For limiting cases of thermal boundary conditions, comparison with an analytical solution is performed.     

Acoustic streaming in micro-scale cylindrical channels

The focus of this work is to extend the theory of boundary layer induced acoustic streaming to include cylindrical geometries and to highlight the effects of boundary layer induced streaming on flow velocities in micro-scale channels. The work presented here includes the development of a model for streaming in a cylindrical channel by a method of successive approximations. The validity of this model is established by comparison with a well-established model for streaming between parallel plates of infinite extent. This is followed by a discussion on the importance of employing a cylindrical solution including boundary layer induced streaming for the analysis of streaming in micro-scale channels.     

Acoustic force model for the fluid flow under standing waves

An acoustic Strouhal number is introduced to demonstrate that the viscosity of fluid can be ignored in the process of sound propagation and acoustic streaming is independent on the frequency of the acoustic wave. Furthermore, acoustic force based on the periodic velocity fluctuation caused by standing acoustic wave is introduced into Navier–Stokes equation in order to describe the fluid flow in the acoustic boundary layer. The numerical results show that the predicted results are consistent with the analytic solution. And the effect of the nonlinear terms cannot be ignored so the analytic solution derived by boundary-velocity condition is only an approximation for acoustic streaming.     

声波作用下球形颗粒外声流分布的数值模拟

综合考虑声学边界层内的热损失和黏性损失,建立处于平面驻波声压波节位置二维球形颗粒外声流计算模型,利用分离时间尺度的数值方法对颗粒外声流流场特征进行模拟.将模拟结果与相应的解析解和实验结果对比,验证了数值模拟的可靠性.在此基础上,研究了雷诺数Re和斯特劳哈尔数Sr对球形颗粒声学边界层内二阶声流流场结构、涡流强度及范围的影...     

Acoustic streaming in an ultrasonic air pump with three-dimensional finite-difference time-domain analysis and comparison to the measurement

The direct finite-difference fluid simulation of acoustic streaming on a fine-meshed three-dimensional model using a graphics processing unit (GPU)-based calculation array is discussed. Airflows are induced by an acoustic traveling wave when an intense sound field is generated in a gap between a bending transducer and a reflector. The calculation results showed good agreement with measurements in a pressure distribution. Several flow vortices were observed near the boundary layer of the reflector and the transducer, which have often been observed near the boundary of acoustic tubes, but have not been observed in previous calculations for this type of ultrasonic air pump.     

Generation of vorticity motion by sound in a chemically reacting gas and inversion of acoustic streaming in the non-equilibrium regime

Nonlinear stimulation of the vorticity mode caused by losses in the momentum of sound in a chemically reacting gas is considered. The instantaneous dynamic equation for the vorticity mode is derived. It includes a quadratic nonlinear acoustic source, which reflects the fact that the reason for the interaction between sound and the vorticity mode is nonlinear. Both periodic and aperiodic sound may be considered as the origin of the vorticity flow. The equation governing the mean flow (the acoustic streaming) in the field of periodic sound is also derived. In the non-equilibrium regime of a chemical reaction, there may exist streaming vortices whose direction of rotation is opposite to that of the vortices in the standard thermoviscous flows. For periodic sound, this is illustrated by an example. The theory and the example describe both equilibrium and non-equilibrium chemical reactions.     

Localization of Rayleigh waves in microfluidic channels with quadratic profile

Localization of Rayleigh waves due to mechanical loading of the substrate surface by a liquid layer with quadratic thickness variation in the transverse direction was theoretically studied. Under conditions typical of acoustic biochips, strong localization of the wave field was revealed. This makes it possible to calculate the lowest waveguide modes of the microfluidic channel, neglecting its finite width.     

A piecewise linear mean flow model for studying stability in a lined channel

Acoustic liners are used to reduce sound emission by turbofan engines. Under grazing flow they may sustain hydrodynamic instabilities and these are studied using a stability analysis, based on a simplified model: the liner is a mass–spring–damper system, the base channel flow is piecewise linear, and the inviscid, incompressible Rayleigh equation is used. The model is an extension to the channel case of a boundary layer model by Rienstra and Darau. The piecewise linear profile introduces a finite boundary layer thickness which ensures well-posedness, allowing an initial value problem to be conducted to investigate absolute stability. For typical values in aeronautics the flow above the liner is unstable. Absolute instability is obtained for somewhat extreme values of the mean flow (tiny boundary layer thickness), and under realistic conditions the flow is convectively unstable. The effect of finite channel height is investigated in both cases. In particular, for large boundary layer thicknesses associated with convective instability the channel height has little effect on the unstable mode. Favorable outcomes and failures of the model are shown by comparison to a published experimental work.     

Characterization of acoustic streaming in water and aluminum melt during ultrasonic irradiation

It is well known that ultrasonic cavitation causes a steady flow termed acoustic streaming. In the present study, the velocity of acoustic streaming in water and molten aluminum is measured. The method is based on the measurement of oscillation frequency of Karman vortices around a cylinder immersed into liquid. For the case of acoustic streaming in molten metal, such measurements were performed for the first time. Four types of experiments were conducted in the present study: (1) Particle Image Velocimetry (PIV) measurement in a water bath to measure the acoustic streaming velocity visually, (2) frequency measurement of Karman vortices generated around a cylinder in water, and (3) in aluminum melt, and (4) cavitation intensity measurements in molten aluminum. Based on the measurement results (1) and (2), the Strouhal number for acoustic streaming was determined. Then, using the same Strouhal number and measuring oscillation frequency of Karman vortices in aluminum melt, the acoustic streaming velocity was measured. The velocity of acoustic streaming was found to be independent of amplitude of sonotrode tip oscillation both in water and aluminum melt. This can be explained by the effect of acoustic shielding and liquid density.     

Ion acoustic solitons in a relativistic warm plasma with densitygradient

Modified Korteweg-deVries equation (mK-dV), which governs the behavior of ion acoustic solitons in a relativistic warm plasma with density gradient, is derived. The electron inertia is also taken into account which is important when the streaming ions are present in the plasma. A solution of the mK-dV equation is obtained for the constant density gradient. When the ion acoustic soliton propagates into the lower plasma density region, its amplitude and energy increase, but the width decreases; the same is the case for the stronger density gradients     

Analytic model of acoustic streaming in thermoacoustic waveguides with slowly varying cross-section

An analytic model of acoustic streaming generated in two-dimensional thermoa-coustic waveguides with slowly varying cross-section was developed for more general applications. The analytical solutions of acoustic streaming characteristics in the closed straight tube and the annular tube are given based on the model.The solution for the closed straight tube can be applied to the case with any transverse scale.The solution for the annular tube is obtained under the assumption that the width of the varying cross-section part is much larger than the viscous and thermal penetration depths.The effects of cross-section variation,time-averaged temperature distribution and components of sound field are reflected in the analytic solutions. The magnitude and distribution of acoustic streaming velocity vary with the characteristic scale of the waveguides.The analytic model of acoustic streaming can be applied in research under thermoacoustic and other physical backgrounds.     

渐变截面热声波导管内的声流解析模型

建立了二维渐变截面热声波导管内的声流模型,分别考虑了封闭直管和环形回路两种不同结构,获得了更为普适的解析结果。封闭直管结构的声流结果可应用于任意宽度的波导管,环路结构的结果考虑了渐变截面管段宽度远大于热、黏穿透深度的情形。研究结果表明,渐变截面热声波导管内的声流主要受声场结构、截面变化及轴向时均温度分布的影响,在其它参量不变时声流量值及分布随波导管特征尺度的不同而变化。该解析模型可应用于热声及其它物理背景下的声流分析。      

Propagation of shock waves in a medium with the Rayleigh energy release mechanism

The effect of the familiar Rayleigh mechanism of energy release in an elastic medium (which plays an important role, in particular, in gas discharge plasma) on the structure of a running shock wave (SW) is considered in the general case in the 1D approximation. The equation describing the propagation of the SW in this case is derived. An analytic solution to this equation is obtained for small values of the parameter characterizing the properties of the medium. The type of the solution for different signs of this parameter and for its values modulo equal to unity is analyzed. It is found that, for positive values of this parameter, a SW in the form of a step is suppressed in such a medium and degenerated into a perturbation in the form of a hump. On the contrary, for negative values of the parameter, the SW is enhanced. It is found that a stationary solution exists in the system of coordinates associated with the SW propagation in a medium with the Rayleigh energy release mechanism only if the boundary of the medium lies downstream from the shock layer. The position of this boundary corresponds to the so-called critical energy supply and the local Mach number is equal to unity at this point. For a positive value of the parameter of the medium with the Rayleigh energy release mechanism, the equation of propagation has no stationary solution for any position of the boundary of the medium upstream from the shock layer when the value of the parameter exceeds a certain limiting value. The results make it possible to analyze the features of SW propagation in a weakly ionized gas discharge plasma.     

Acoustic streaming in a rotating fluid

Acoustic streaming theory is derived that is applicable to a fluid that is slow moving in a reference frame that rotates with a constant angular velocity omega. A simplified streaming equation is obtained for the special case in which the acoustic angular frequency omega is large relative to omega, and the change in fluid density due to rotation alone is negligible. For this special case it is shown that the "driving force" for the acoustic streaming is independent of omega. Thus, if no acoustic streaming is present in a fluid system that is stationary, then no steady-state acoustic streaming is predicted for a similar system that rotates with constant angular velocity. For a system in which acoustic streaming is present, the flow behavior depends on the relative magnitudes of the Coriolis forces and the viscous forces. If the Ekman number is large (that is, the viscous force dominates) then the predicted flow is identical to that which would exist in a stationary system. If, on the other hand, the Ekman number is small then the Coriolis force dominates and the component of flow in the direction of the axis of rotation can be much smaller in the rotating system than in a similar system at rest.     

Measurements of inner and outer streaming vortices in a standing waveguide using laser doppler velocimetry

Measurements of the axial streaming velocity are performed by means of laser doppler velocimetry in an experimental apparatus consisting of a waveguide having loudspeakers at each end for high intensity sound levels. Streaming is characterized by an appropriate Reynolds number Re(NL), the case Re(NL)<1 corresponding to the so-called slow streaming and the case Re(NL)>/=1 being referred to as fast streaming. The variation of axial streaming velocity with respect to the transverse coordinate is compared to the available slow streaming theory. Streaming fluid flow is measured both in the core region and in the near wall region. Streaming velocity in the center of the guide agrees reasonably well with the slow streaming theory for small Re(NL) but deviates significantly from such predictions for Re(NL)>20 and its evolution for further increasing Re(NL) is discussed. Then streaming behavior in the near wall region is particularly studied. For Re(NL)<70, two vortices are present across the guide section as predicted by slow streaming theory. Then it appears that, when the Reynolds number is increased, two other vortices become visible in the near wall region. Different stages for the generation and evolution of these inner streaming vortices are presented.     

Analysis of the flame thickness of turbulent flamelets in the thin reaction zones regime

The thickness of the instantaneous flamelets in a turbulent flame brush on a weak-swirl burner burning in the thin reaction zones regime has been analysed experimentally, theoretically, and numerically. The experimental flame thickness has been measured correlating two simultaneous Rayleigh images and one OH-image from two closely spaced cross sections in the flame. It appears that the low temperature edge of the flame is thickened by turbulent eddies but that these structures cannot penetrate far enough into the flame front to distort the inner layer for the moderate Karlovitz numbers used. The flame front based on the temperature gradient at the inner layer becomes thinner for lean flames and thicker for rich methane–air flames. This has been explained theoretically and numerically by studying the influence of flame stretch and preferential diffusion on the flame thickness. It appears that the flame front thickness at the inner layer (and mass burning rate) is not influenced by turbulent mixing processes, and it seems that eddies of the size of the inner layer have to be used to change this picture. Experiments closer to the boundary of the broken reaction zones regime have to confirm this in the future.     

Compressibility effects on steady streaming from a noncompact rigid sphere

The problem of steady streaming around a rigid isolated sphere in a plane standing acoustic field is considered. Existing results in the literature have been generalized to allow for noncompactness of the sphere, and the influence of fluid compressibility on the streaming behavior has been included. It is found that in the high-frequency limit of interest for which the streaming is strongest, the effective steady slip velocity at the edge of the inner boundary layer region that is responsible for driving the steady streaming in the bulk of the fluid in the outer region, has a complex variation over the surface of the sphere that depends on (i) the sphere position (with respect to the node/antinode of the acoustic field), (ii) the extent of sphere compactness, and (iii) on a well-defined function (representing compressibility effects) of the fluid Prandtl number and its ratio of specific heats. Not surprisingly, the contribution from this function is negligible when the host fluid is a liquid. The steady streaming behavior around the sphere is demonstrated with the help of flow streamlines for various cases in the diffusive limit of weak outer flow for low streaming Reynolds numbers.     

Controlling acoustic streaming in an ultrasonic heptagonal tweezers with application to cell manipulation

Acoustic radiation force has been demonstrated as a method for manipulating micron-scale particles, but is frequently affected by unwanted streaming. In this paper the streaming in a multi-transducer quasi-standing wave acoustic particle manipulation device is assessed, and found to be dominated by a form of Eckart streaming. The experimentally observed streaming takes the form of two main vortices that have their highest velocity in the region where the standing wave is established. A finite element model is developed that agrees well with experimental results, and shows that the Reynolds stresses that give rise to the fluid motion are strongest in the high velocity region. A technical solution to reduce the streaming is explored that entails the introduction of a biocompatible agar gel layer at the bottom of the chamber so as to reduce the fluid depth and volume. By this means, we reduce the region of fluid that experiences the Reynolds stresses; the viscous drag per unit volume of fluid is also increased. Particle Image Velocimetry data is used to observe the streaming as a function of agar-modified cavity depth. It was found that, in an optimised structure, Eckart streaming could be reduced to negligible levels so that we could make a sonotweezers device with a large working area of up to 13 mm × 13 mm.