纳米粒子在方管中传输与沉降的研究

采用有限体积法离散并应用Simple方法对方截面弯曲管道内的纳米粒子传输和沉降进行了数值计算,结果表明Reynolds数和Schmidt数是影响纳米粒子传输和沉降的重要参数。粒子较小时,弯管中轴向速度较大的区域就是粒子的高浓度区域,沉降增强因子最大值出现在外弯侧的中心位置;粒子较大时,截面浓度的梯度值降低,沉降增强因子趋于平均,此时整个截面的粒子平均沉降。弯曲作用对于粒子较大且Dean数也较大时的影响更加明显。     

Solution of general dynamic equation for nanoparticles in turbulent flow considering fluctuating coagulation

A new averaged general dynamic equation (GDE) for nanoparticles in the turbulent flow is derived by considering the combined effect of convection, Brownian diffusion, turbulent diffusion, turbulent coagulation, and fluctuating coagulation. The equation is solved with the Taylor-series expansion moment method in a turbulent pipe flow. The experiments are performed. The numerical results of particle size distribution correlate well with the experimental data. The results show that, for a turbulent nanoparticulate flow, a fluctuating coagulation term should be included in the averaged particle GDE. The larger the Schmidt number is and the lower the Reynolds number is, the smaller the value of ratio of particle diameter at the outlet to that at the inlet is. At the outlet, the particle number concentration increases from the near-wall region to the near-center region. The larger the Schmidt number is and the higher the Reynolds number is, the larger the difference in particle number concentration between the near-wall region and near-center region is. Particle polydispersity increases from the near-center region to the near-wall region. The particles with a smaller Schmidt number and the flow with a higher Reynolds number show a higher polydispersity. The degree of particle polydispersity is higher considering fluctuating coagulation than that without considering fluctuating coagulation.     

Transport and deposition of neutral particles in magnetohydrodynamic turbulent channel flows at low magnetic Reynolds numbers

The effect of Lorentz force on particle transport and deposition is studied by using direct numerical simulation of turbulent channel flow of electrically conducting fluids combined with discrete particle simulation of the trajectories of uncharged, spherical particles. The magnetohydrodynamic equations for fluid flows at low magnetic Reynolds numbers are adopted. The particle motion is determined by the drag, added mass, and pressure gradient forces. Results are obtained for flows with particle ensembles of various densities and diameters in the presence of streamwise, wall-normal or spanwise magnetic fields. It is found that the particle dispersion in the wall-normal and spanwise directions is decreased due to the changes of the underlying fluid turbulence by the Lorentz force, while it is increased in the streamwise direction. The particle accumulation in the near-wall region is diminished in the magnetohydrodynamic flows. In addition, the tendency of small inertia particles to concentrate preferentially in the low-speed streaks near the walls is strengthened with increasing Hartmann number. The particle transport by turbophoretic drift and turbulent diffusion is damped by the magnetic field and, consequently, particle deposition is reduced.     

Hydrodynamic force and torque models for a particle moving near a wall at finite particle Reynolds numbers

This research work is aimed at proposing models for the hydrodynamic force and torque experienced by a spherical particle moving near a solid wall in a viscous fluid at finite particle Reynolds numbers. Conventional lubrication theory was developed based on the theory of Stokes flow around the particle at vanishing particle Reynolds number. In order to account for the effects of finite particle Reynolds number on the models for hydrodynamic force and torque near a wall, we use four types of simple motions at different particle Reynolds numbers. Using the lattice Boltzmann method and considering the moving boundary conditions, we fully resolve the flow field near the particle and obtain the models for hydrodynamic force and torque as functions of particle Reynolds number and the dimensionless gap between the particle and the wall. The resolution is up to 50 grids per particle diameter. After comparing numerical results of the coefficients with conventional results based on Stokes flow, we propose new models for hydrodynamic force and torque at different particle Reynolds numbers. It is shown that the particle Reynolds number has a significant impact on the models for hydrodynamic force and torque. Furthermore, the models are validated against general motions of a particle and available modeling results from literature. The proposed models could be used as sub-grid scale models where the flows between particle and wall can not be fully resolved, or be used in Lagrangian simulations of particle-laden flows when particles are close to a wall instead of the currently used models for an isolated particle.     

Lattice Boltzmann simulations for flow and heat/mass transfer problems in a three‐dimensional porous structure

The lattice Boltzmann method (LBM) for a binary miscible fluid mixture is applied to problems of transport phenomena in a three‐dimensional porous structure. Boundary conditions for the particle distribution function of a diffusing component are described in detail. Flow characteristics and concentration profiles of diffusing species at a pore scale in the structure are obtained at various Reynolds numbers. At high Reynolds numbers, the concentration profiles are highly affected by the flow convection and become completely different from those at low Reynolds numbers. The Sherwood numbers are calculated and compared in good agreement with available experimental data. The results indicate that the present method is useful for the investigation of transport phenomena in porous structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Impingement of buoyancy-driven flows at a stratified interface

Laser-induced fluorescence (LIF) and particle-image velocimetry (PIV) are used to study both thermals and plumes impinging on a stratified interface. Data are obtained for a central slice of the flow near the stratified interface. Both the thermal and plume are generated by releasing fresh water at the bottom of a tank filled with two layers of salt water of different densities. Thermals and plumes are studied at Reynolds numbers ranging from 3,000 to 8,000, above the value for the mixing transition, a Schmidt number of about 600, and Richardson numbers from 1 to 22. The Richardson and Reynolds numbers are based on the thermal or plume characteristics (size and vertical velocity) before impingement and the initial density difference across the interface. Laser-induced fluorescence (LIF) is used to determine the maximum penetration height, rebound distance and lateral spreading velocity. The vorticity results obtained from the PIV data reveal the vortical structure near impingement. When the thermal impinges upon the stratified interface, a baroclinic eddy generated at the interface appears to merge with eddies comprising the thermal itself to form a vortex ring. This ring remains near the interface, moving mainly along the lateral or horizontal direction away from the region of impingement. These results suggest that lateral transport is significant for thermals impinging on stratified interfaces, and that ignoring such transport may greatly underestimate overall transport and mixing in such flows.     

DNS of passive scalar transport in turbulent channel flow at high Schmidt numbers

We perform DNS of passive scalar transport in low Reynolds number turbulent channel flow at Schmidt numbers up to Sc = 49. The high resolutions required to resolve the scalar concentration fields at such Schmidt numbers are achieved by a hierarchical algorithm in which only the scalar fields are solved on the grid dictated by the Batchelor scale. The velocity fields are solved on coarser grids and prolonged by a conservative interpolation to the fine-grid.

The trends observed so far at lower Schmidt numbers Sc  10 are confirmed, i.e. the mean scalar gradient steepens at the wall with increasing Schmidt number, the peaks of turbulent quantities increase and move towards the wall. The instantaneous scalar fields show a dramatic change. Observable structures get longer and thinner which is connected with the occurrence of steeper gradients, but the wall concentrations penetrate less deeply into the plateau in the core of the channel.

Our data shows that the thickness of the conductive sublayer, as defined by the intersection point of the linear with the logarithmic asymptote scales with Sc^−0.29. With this information it is possible to derive an expression for the dimensionless transfer coefficient K⁺ which is only dependent on Sc and Re_τ. This expression is in full accordance to previous results which demonstrates that the thickness of the conductive sublayer is the dominating quantity for the mean scalar profile.     

Turbulent heat and mass transfer over a rotating disk for the Prandtl or Schmidt numbers much larger than unity: an integral method

Turbulent heat and mass transfer of a rotating disk for Prandtl and Schmidt numbers much larger than unity was modeled using an integral method validated against empirical equations of different authors for Sherwood numbers. As shown, decrease in relative thickness of thermal/diffusion boundary layers with increasing local radii entails additional increase of the exponent at the Reynolds number in expressions for Nusselt and Sherwood numbers in comparison with air flows.     

On the fluid–structure interaction of a splitter plate: vibration modes and Reynolds number effects

Previous work (Eloranta et al. in Exp Fluids 39:841–855, 2005) has shown that flow separation from the trailing edge of a splitter plate in a convergent channel involves a fluid–structure interaction (FSI), which modifies the fundamental instability related to vortex shedding. Under certain conditions, the FSI induces cellular vortex shedding from the trailing edge. This paper reports detailed measurements of the plate vibration mode and studies the effect of the Reynolds number on the FSI. Experimental techniques including laser vibrometer and digital imaging are used to measure the response of the plate and particle image velocimetry is used to measure the flow field in the near wake. Combining data from these techniques, the development of the vibration frequency and mode can be addressed together with the imprint of the vibration mode in the flow. The results show that over most of the Reynolds numbers measured, the plate vibrates in a distinct mode characterized by a spanwise standing wave along the plate trailing edge. The vibration frequency and the spacing between the nodes of the standing wave depend on the Reynolds number. As the Reynolds number is increased, the frequency of the dominant vibration mode does not increase linearly. The plot of the vibration frequency as a function of the Reynolds number shows that the vibration tends to lock to a rather constant frequency over of range of Reynolds numbers. After certain Reynolds number if threshold is exceeded, the frequency jumps to a new level, which also involves a new vibration mode.     

Dynamics of dual-particles settling under gravity

As a first step towards understanding particle–particle interaction in fluid flows, the motion of two spherical particles settling in close proximity under gravity in Newtonian fluids was investigated experimentally for particle Reynolds numbers ranging from 0.01 to 2000. It was observed that particles repel each other for Re>0.1 and that the separation distance of settling particles is Reynolds number dependent. At lower Reynolds numbers, i.e. for Re<0.1, particles settling under gravity do not separate.The orientation preference of two spherical particles was found to be Reynolds number dependent. At higher Reynolds numbers, the line connecting the centres of the two particles is always horizontal, regardless of the way the two particles are launched. At lower Reynolds numbers, however, the particle centreline tends to tilt to an arbitrary angle, even of the two particles are launched in the horizontal plane. Because of the tilt, a side migration of the two particles was found to exist. A linear theory was developed to estimate the side migration velocity. It was found that the maximum side migration velocity is approximately 6% of the vertical settling velocity, in good agreement with the experimental results.Counter-rotating spinning of the two particles was observed and measured in the range of Re=0–10. Using the linear model, it is possible to estimate the influence of the tilt angle on the rate of rotation at low Reynolds numbers. Dual particles settle faster than a single particle at small Reynolds numbers but not at higher Reynolds numbers, because of particle separation. The variation of particle settling velocity with Reynolds number is presented. An equation which can be used to estimate the influence of tilt angle on particle settling velocity at low Reynolds number is also derived.     

Effect of thermophoresis on fog droplet deposition on low pressure steam turbine guide blades

Deposition of droplets in the sub-micron size range, when discouraged by thermophoresis forces, has been studied by a simulation method employing uranine particles entrained in air and flowing through a cascade of full-sized low-pressure steam turbine blades. The particles, generated by a Collison atomiser, had a mass-median diameter range of 0.05–0.25 μm. The test blades were internally heated using hot air and had an output of 600 W/m² of swept surface. Circumferential tape around their surfaces provided a reception medium for the particles, the deposition density variation of which was found by fluorimetry. Prediction of the deposition on an unheated blade using the method of C.N. Davies showed acceptable agreement with the experimental results. The simulation method was validated by showing that the relevant non-dimensional numbers. Schmidt and Reynolds, can be acceptably matched for the air/particle and the steam/droplet cases. Tests on the heated blade showed that the deposition was reduced by 30–90% of the corresponding value for the unheated blade. The extent of the reduction decreased with particle size decrease.     

基于LBM-LES方法翼型纯音噪声数值研究

基于LBM-LES方法,对中低雷诺数下的NACA0012翼型纯音噪声进行了直接计算,研究了不同迎角和雷诺数对纯音噪声的影响。计算结果表明,翼型的声源主要位于翼型的分离区和后缘处,在不同迎角和雷诺数下的声辐射特征均具有偶极子声场的特点;迎角的增大将引起较大的旋涡尺度和湍流强度,吸力面声源区域前移。声压级频谱分析表明,随着迎角的增大,纯音噪声逐渐消失,噪声谱最终呈现宽频特征;随着雷诺数的增大,后缘压力脉动增大。声压级频谱中,主频频率随着雷诺数的增大而增大,且符合Paterson公式的幂律关系。此外,声压级频谱特性随着雷诺数的增大表现出由离散特性向宽频特性转变的趋势。     

Evaporation phenomenon past a rotating hydrocarbon droplet of ternary components

A numerical study of heat and mass transfer from an evaporating fuel droplet rotating around its vertical axis was performed in forced convection only on the side opposite to the flow. The flow was assumed to be laminar, and the droplet was assumed to maintain its spherical shape during its lifetime. Based on the abovementioned assumption, the conservation equations in a general curvilinear coordinate were solved numerically. The behavior of rotating droplet evaporation in the forced convection flow can be investigated by analyzing the effects of the rotation of the droplet on the evaporation process of multi-component hydrocarbons droplet. The droplet is simulated to behave as a hard sphere. The transfer equations are discretized using an implicit finite difference method. Thomas algorithm is used to solve the system of algebraic equations. Moreover, dimensionless parameters of heat and mass transfer phenomena around a rotating hydrocarbon droplet were determined. The thickness of the boundary layer is unknown for this model and therefore, it was determined in function of time. Additionally, the study concerns “Dgheim dimensionless number” which is the ratio of the rotation forces over the viscosity forces. Dgheim dimensionless number is correlated to Nusselt and Sherwood numbers for multi-component hydrocarbon droplets in evaporation by taking into account the effect of heat and mass Spalding, Prandtl and Schmidt numbers respectively. Also, correlations for Nusselt and Sherwood numbers in terms of Reynolds, Prandtl and Schmidt numbers are proposed. These correlations consider the rotation phenomenon and advance the variation of the thermophysical and transport properties in the vapor phase of multi-component blends.     

Mass transfer to rotating cones

The mass transfer of oxygen to a rotating cone has been measured in a quiescent fluid. Cones with vertex angles of 19.2°, 30°, 60°, 90°, 120° and 180° (disks) were used in both laminar and turbulent flow at a Schmidt number of 400. The laminar flow mass transfer data is in agreement with an equation derived for cones at high Schmidt numbers. The turbulent flow mass transfer data for cones with vertex angles greater than 19.2° is in agreement with the high Schmidt number expression of Deissler. The 19.2° cone angle, however, results in experimental values for the mass transfer which are greater than the predicted values. Lastly, the values of the Reynolds number for transition from laminar to turbulent flow, as determined directly from the mass transfer data, are in agreement with previous results.     

Überlegungen zur Scherviskosität von Suspensionen aufgrund einer Dimensionsanalyse

Dimensional analysis of the motion of solid particles suspended in a fluid phase shows that the macroscopic relative shear viscosity of suspensions generally depends not only on the volume concentration and particle shape but also on two Reynolds numbers and a dimensionless sedimentation number. These dimensionless numbers are formed using parameters characterizing the structure and motion of the suspension at the microscopic level. The analysis was based on the assumptions that the dispersed particles are rigid and sufficiently large that Brownian motion may be neglected, that the continuous fluid phase is Newtonian and that the interactions between particles and between particles and fluid phase are only hydrodynamic. The Reynolds numbers describe the influence of the inertial forces at the microscopic level, and the sedimentation number the influence of gravity. The dimensionless numbers can be neglected if their values are much smaller than one. For each of the dimensionless numbers both the shear rate and the particle size influence the shear viscosity. Thus sedimentation number is large for low shear rates, whereas the Reynolds numbers are large for high shear rates. The viscosity function for one suspension can be transformed into the viscosity function for another suspension with geometrically similar particles but of a different size. The scale-up rules are derived from the requirement that the relevant dimensionless numbers must be constant. The influence of non-hydrodynamic effects at the microscopic level on the shear viscosity can be detected by deviations from the derived scale-up rules.     

An experimental study of drag on a single tube and on a tube in an array under two-phase cross flow

In this work, the drag coefficient and the void fraction around a tube subjected to two-phase cross flow were studied for a single tube and for a tube placed in an array. The drag coefficients were determined by measuring the pressure distribution around the perimeter of the tube. Single tube drag data were taken when the tube was held both rigidly and flexibly. The test tube was made of acrylic and was 2.2 cm in diameter and 20 cm in length. In the experiments, liquid Reynolds number ranged from 430 to 21,900 for the single tube and liquid gap Reynolds number ranged from 32,900 and 61,600 for the tube placed in a triangular array. Free stream void fraction was varied from 0 to 0.4. At low Reynolds numbers, the ratio of two-phase to single-phase drag coefficient is found to be a strong function of εGr/Re2. However, at high Reynolds numbers only void fraction is the important parameter. Empirical correlations have been developed for the ratio of two-phase drag on a single tube and on a tube placed in an array.     

低雷诺数俯仰振荡翼型等离子体流动控制

针对低雷诺数翼型气动性能差的特点, 通过介质阻挡放电(dielectric barrier discharge, DBD)等离子体激励控制的方法, 提高翼型低雷诺数下的气动特性,改善其流场结构. 采用二维准直接数值模拟方法求解非定常不可压Navier-Stokes方程,对具有俯仰运动的NACA0012翼型的低雷诺数流动展开数值模拟.同时将介质阻挡放电激励对流动的作用以彻体力源项的形式加入Navier-Stokes方程,通过数值模拟探究稳态DBD等离子体激励对俯仰振荡NACA0012翼型气动特性和流场特性的影响.为了进行流动控制, 分别在上下表面的前缘和后缘处安装DBD等离子体激励器,并提出四种激励器的开环控制策略,通过对比研究了这些控制策略在不同雷诺数、不同减缩频率以及激励位置下的控制效果.通过流场结构和动态压强分析了等离子体进行流场控制的机理. 结果表明,前缘DBD控制中控制策略B(负攻角时开启上表面激励器,正攻角时开启下表面激励器)效果最好,后缘DBD控制中控制策略C(逆时针旋转时开启上表面激励器,顺时针旋转时开启下表面激励器)效果最好,前缘DBD控制效果会随着减缩频率的增大而下降, 同时会导致阻力增大.而后缘DBD控制可以减小压差阻力, 优于前缘DBD控制,对于计算的所有减缩频率(5.01~11.82)都有较好的增升减阻效果.在不同雷诺数下, DBD控制的增升效果较为稳定, 而减阻效果随着雷诺数的降低而变差,这是由流体黏性效应增强导致的.      

Laboratory investigation and commercial test for rotors-assembled strand applied in smooth tube

In the commercial test for smooth tube inserted with rotors-assembled strand comparing with non-inserted ones on condensers in electric power plant, using water as working fluid, the single-phase pressure drop and heat transfer were measured. It was found difficult to receive reliable and accurate enough data through commercial test. Meanwhile, the single-phase pressure drop and heat transfer in a rotors-assembled strand inserted tube were measured in laboratory, with the tube side Prandtl numbers varying from 5.67 to 5.80 and the tube side Reynolds numbers varying from 21,300 to 72,200. Before that, a validation experiment based on the same smooth tube was carried out to testify the experimental system and the data reduction method, in which fixed mounts were employed to eliminate entrance effects. The Prandtl numbers varied from 5.64 to 5.76 and the Reynolds numbers varied from 19,000 to 56,000 in the tube. The annular side Reynolds numbers remained nearly constant at the value of around 50,000 for all experiments, with the annular side Prandtl numbers varying from 8.02 to 8.22. The experimental results of smooth tube show that employment of fixed mounts leads to a visible bias of friction factor at relative low Reynolds numbers while it hardly affects the Nusselt numbers. On the other hand, experiment for the tube inserted with rotors-assembled strand show remarkable improvement for heat transfer with the Nusselt number increased by 9.764–11.87% and the overall heat transfer coefficient increased by 7.08–7.49% within the range of Reynolds number from about 21,300 to 55,500. Meanwhile, friction factor increases inevitably by 278.1–353.9% within the same range of Reynolds number. Based on through multivariant linear normal regression method, the Reynolds number and Prandtl number dependencies of the Nusselt number and friction factor were determined to be Nu = 0.0031Re^0.9Pr^1.0849 and f = 0.993Re^−0.22.     

Simulations of Three-Dimensional Turbulent Mixing for Schmidt Numbers of the Order 1000

We report basic results from new numerical simulations of passive scalar mixing at Schmidt numbers (Sc) of the order of 1000 in isotropic turbulence. The required high grid-resolution is made possible by simulating turbulence at very low Reynolds numbers, which nevertheless possesses universality in dissipative scales of motion. The results obtained are qualitatively consistent with those based on another study (Yeung et al., Phys. Fluids 14 (2002) 4178-4191) with a less extended Schmidt number range and a higher Reynolds number. In the stationary state maintained by a uniform mean scalar gradient, the scalar variance increases slightly with Sc but scalar dissipation is nearly constant. As the Schmidt number increases, there is an increasing trend towards k ^?1 scaling predicted by Batchelor (Batchelor, J. Fluid Mech. 5 (1959) 113-133) for the viscous-convective range of the scalar spectrum; the scalar gradient skewness approaches zero; and the intermittency measured by the scalar gradient flatness approaches its asymptotic state. However, the value of Sc needed for the asymptotic behavior to emerge appears to increase with decreasing Reynolds number of the turbulence. In the viscous-diffusive range, the scalar spectrum is in better agreement with Kraichnan's (Kraichnan., Phys. Fluids 11 (1968) 945-953) result than with Batchelor's.     

Sensitivity analyses in a buoyancy-driven closed system with high resolution CFD using Boussinesq approximation and variable density models

The influences of the Schmidt number and fluid models of the Boussinesq approximation and variable density (low-Mach formulation) are investigated using large eddy simulation (LES) where fluid naturally advects between two tanks connected by a pipe. The geometry resembles a reactor-like vessel which comprises of a cold leg-downcomer region. The closed system consists of two miscible fluids with different densities and viscosities separated by a valve. Analysis with Reynolds-Averaged Navier–Stokes (RANS) is used for calculating Taylor length scales to establish minimum grid resolution for LES. Using particle image velocimetry (PIV) data, validation is conducted for four different magnitudes of Schmidt numbers. Evidence indicates that capturing the effects from the highest Schmidt number is essential for properly predicting fluctuations and instabilities in the flow. Although downcomer comparisons with PIV consist of larger errors, they are explained through additional sensitivity analyses involving initial conditions, geometry, and residual energy. Lastly, the influence of larger density differences is demonstrated using simulations with a higher Atwood number, temporally scaled using the fluid front velocity. Recommendations are made regarding future experiments and simulations for design of similar buoyant mixing problems.