振荡潜流带沉积层?水界面污染物输运的研究

潜流带中污染物质交换与输运特性是影响水资源环境的重要问题之一. 本文对底部为高渗透沉积层的三维槽道振荡流高Schmidt数传质问题进行了大涡模拟研究. 采用动力学亚格子模型来封闭滤波后的三维不可压缩Navier-Stokes方程以及污染物输运方程, 同时采用修正的Darcy-Brinkman-Forcheimer模型来描述沉淀有锌离子污染溶质的可渗透沉积层. 通过对沉积层内外流场和浓度场的统计特性以及瞬态结构的分析, 探究了上覆水体中振荡流驱动作用对污染物输运的动力学影响以及扩散率随振荡周期和振荡角的变化规律. 研究结果表明, 浓度通量中的湍流浓度分量在垂向物质交换中起主导作用, 流向、展向速度, 湍流强度和污染物浓度的波动跟随振荡驱动力呈现准周期变化, 同时发现沉积层?水交界面处的湍流浓度通量与法向湍流强度两者之间的变化具有明确的相关性. 并且在较大振荡角和低频振荡的情况下, 沉积层?水交界面处的有效扩散率增大, 这主要是来自于沉积层?水交界面处流体的猝发行为有效促进了湍流混合和物质交换, 进而增强了污染物的垂向输运.      

Effect of flow-obstruction geometry on pressure drops in horizontal air-water flow

Experiments were carried out to investigate the pressure drop due to flow obstructions in horizontal air-water flow. The axial pressure distribution along a 25.4 mm inside diameter tube, with and without flow obstructions was measured using multi-tube manometers. Various obstruction shapes and sizes were investigated. Pressure-loss coefficients and two-phase multipliers were derived for twelve different flow obstructions.Estimates of the kinetic energy and momentum of the flow were also obtained from radial void distribution measurements. The pressure drop through the obstructions in two-phase flow was found to depend strongly on the kinetic energy and momentum of the liquid intercepted by the flow obstruction. Buoyancy-induced flow stratification caused a strongly non-symmetrical effect, and resulted in large pressure drops for flow obstructions located in the bottom part of the channel.     

Laminar flow and mass transfer in channels with a porous bottom wall and with fins attached to the top wall

 The steady incompressible, viscous, two- dimensional flow of a solution in a channel was considered. The bottom wall was porous and the fins were attached to the top wall. Employing control volume approach, a computer program based on SIMPLE algorithm was developed. Computations were carried out to investigate the effects of the inlet Reynolds number, the fin length, the suction Reynolds number and the slip coefficient on the flow structure and the concentration distribution. It was observed that the thickness of concentration boundary layer increases in the flow direction. The concentration on the porous wall and the concentration boundary layer thickness decrease with increasing fin length, the slip coefficient and the inlet Reynolds number. These results show that fins attached to the upper wall of the channel can be utilized to reduce the concentration polarization and hence improve the effectiveness of the separation process. Received on 24 February 1999     

NUMERICAL ANALYSIS ON THE VELOCITY AND PRESSURE FIELDS INDUCED BYMULTI-OSCILLATIONS OF AN UNDERWATER EXPLOSION BUBBLE

假设水下爆炸气泡的内部气体在膨胀收缩过程中满足绝热条件,周围流体无黏无旋不可压缩. 基于势流理论,采用边界元法研究气泡动力学行为,重点关注气泡引起的流场脉动载荷以及滞后流特性,给出了相关的理论推导和数值计算方法. 通过将数值结果与解析解、实验值进行对比,数值模型的收敛性和有效性能够得到保证. 利用编写的程序进行计算和分析,发现在气泡加速膨胀阶段,流场压力在气泡径向不一定是逐渐衰减,还有可能以先增后减的规律变化;气泡射流后,为了能够继续描述环状气泡的运动以及流场特性,将此时的流场分为无旋场和一个布置在气泡内部涡环的叠加,计算过程中采用了一些数值技巧处理气泡的拓扑结构,得以连续模拟多个周期的气泡运动. 环状气泡具有相对较高的上浮迁移速度,而且在其顶部和底部附近分别形成两个高压区,顶部的高压区峰值相对较大,底部的高压区范围相对较大. 环状气泡中心轴上的流场速度会在气泡中心有一个加速过程,在气泡顶部附近又迅速减小.     

气泡多周期运动时引起的流场压力与速度

假设水下爆炸气泡的内部气体在膨胀收缩过程中满足绝热条件,周围流体无黏无旋不可压缩. 基于势流理论,采用边界元法研究气泡动力学行为,重点关注气泡引起的流场脉动载荷以及滞后流特性,给出了相关的理论推导和数值计算方法. 通过将数值结果与解析解、实验值进行对比,数值模型的收敛性和有效性能够得到保证. 利用编写的程序进行计算和分析,发现在气泡加速膨胀阶段,流场压力在气泡径向不一定是逐渐衰减,还有可能以先增后减的规律变化;气泡射流后,为了能够继续描述环状气泡的运动以及流场特性,将此时的流场分为无旋场和一个布置在气泡内部涡环的叠加,计算过程中采用了一些数值技巧处理气泡的拓扑结构,得以连续模拟多个周期的气泡运动. 环状气泡具有相对较高的上浮迁移速度,而且在其顶部和底部附近分别形成两个高压区,顶部的高压区峰值相对较大,底部的高压区范围相对较大. 环状气泡中心轴上的流场速度会在气泡中心有一个加速过程,在气泡顶部附近又迅速减小.      

A comparative study of flame-wall interaction and flame-cooling air interaction

This paper presents a three-dimensional (3D) direct numerical simulation (DNS) study of flame-wall interaction (FWI) and flame-cooling air interaction (FCAI). A preheated, methane/air mixture enters a channel with constant temperature walls, where the top wall is effusion cooled. An imposed vertical hot sheet near the inlet creates two flame branches interacting with the top and bottom walls. The flame is observed to be leaner in the region where it interacts with the effusion cooling jets. In this region, the flame is longer and features reduced CO mass fraction. The fluctuations in the heat release rate (HRR) and CO mass fraction are also relatively small near the top wall. Near the bottom wall, finger-like flame structures are formed due to the interaction of turbulent vortices with the flame surface. These flame structures initially move away from the wall as they propagate further downstream before eventually collapsing at the wall. This leads to the creation of regions of high wall heat flux and CO. While analysis of the CO thermochemical state shows a complex picture near the bottom wall, two-dimensional (2D) manifolds can be identified near the top wall. Therefore, a framework to estimate CO mass fraction due to FCAI based on 1D freely-propagating flame solutions is proposed showing a good agreement with the DNS results.     

Statistics of Particle Accumulation in Spatially Developing Turbulent Boundary Layers

We present the results of a Direct Numerical Simulation of a particle-laden spatially developing turbulent boundary layer up to Re _θ ?=?2500. Two main features differentiate the behavior of inertial particles in a zero-pressure-gradient turbulent boundary layer from the more commonly studied case of a parallel channel flow. The first is the variation along the streamwise direction of the local dimensionless parameters defining the fluid-particle interactions. The second is the coexistence of an irrotational free-stream and a near-wall rotational turbulent flow. As concerns the first issue, an inner and an outer Stokes number can be defined using inner and outer flow units. The inner Stokes number governs the near-wall behavior similarly to the case of channel flow. To understand the effect of a laminar-turbulent interface, we examine the behavior of particles initially released in the free stream and show that they present a distinct behavior with respect to those directly injected inside the boundary layer. A region of minimum concentration occurs inside the turbulent boundary layer at about one displacement thickness from the wall. Its formation is due to the competition between two transport mechanisms: a relatively slow turbulent diffusion towards the buffer layer and a fast turbophoretic drift towards the wall.     

Forced Convection Heat Transfer from a Finite-Height Cylinder

This paper presents a large eddy simulation of forced convection heat transfer in the flow around a surface-mounted finite-height circular cylinder. The study was carried out for a cylinder with height-to-diameter ratio of 2.5, a Reynolds number based on the cylinder diameter of 44 000 and a Prandtl number of 1. Only the surface of the cylinder is heated while the bottom wall and the inflow are kept at a lower fixed temperature. The approach flow boundary layer had a thickness of about 10% of the cylinder height. Local and averaged heat transfer coefficients are presented. The heat transfer coefficient is strongly affected by the free-end of the cylinder. As a result of the flow over the top being downwashed behind the cylinder, a vortex-shedding process does not occur in the upper part, leading to a lower value of the local heat transfer coefficient in that region. In the lower region, vortex-shedding takes place leading to higher values of the local heat transfer coefficient. The circumferentially averaged heat transfer coefficient is 20 % higher near the ground than near the top of the cylinder. The spreading and dilution of the mean temperature field in the wake of the cylinder are also discussed.     

Longitudinal diffusion of solid particle admixture in a flow through a tube

The propagation of solid particle admixture in a flow through a flat channel is studied.The processes of diffusion and convective transfer as well as solid particle deposition due to gravity result in varying admixture concentration both in depth and longtitudinally.The study of admixture longitudinal distribution is of great interest in a lot of applications, therefore this paper gives the derivation of longitudinal diffusion equation for a mean cross-section admixture concentration.The equation contains three effective parameters; i.e. convective tranfer velocity, longitudinal diffusion coefficient and particle deposition time. These parameters integrally reflect local processes of matter transfer as well as momentum.The proposed model is specific and differs from Taylor equation for longitudinal diffusion, since the fact of particle deposition and adhesion is taken into account. As a result of particle deposition a sediment layer is formed on the channel bottom which increases in thickness with time. To describe this process balance conditions for the whole flow mass and admixture mass on sediment sediment surface are formulated and a condition for matter movement towards the channel bottom is derived that is different from zero due to particle adhesion.     

Transient mixed convection flow arising due to thermal and mass diffusion over porous sensor surface inside squeezing horizontal channel

The double diffusion effect on the mixed convection flow over a horizontal porous sensor surface placed inside a horizontal channel is analyzed.With the appropriate transformations,the unsteady equations governing the flow are reduced to non-similar boundary layer equations which are solved numerically for the time-dependent mixed convection parameter.The asymptotic solutions are obtained for small and large values of the time-dependent mixed convection parameter.The results are discussed in terms of the skin friction,the heat transfer coefficient,the mass transfer coefficient,and the velocity,temperature,and concentration profiles for different values of the Prandtl number,the Schmidt number,the squeezing index,and the mixed convection parameter.     

Surface pressure on a cubic building exerted by conical vortices

Experiments were performed to study surface pressure on a cubic building underlying conical vortices, which are known to cause severe structural damage and failure. The focus is on the effects of turbulence in the incident flow. Three turbulent boundary layers were created in a boundary layer wind tunnel. A wall-mounted cube, i.e. a cube situated on the horizontal ground floor surface of the wind-tunnel test section, was used as an experimental model. The cube was subjected to the incidence flow at 40°. Steady and unsteady pressure measurements were performed on the cube surface. The analysis suggests that conical vortices developed above the top surface of the wall-mounted cube. A larger mean suction was observed on the top cube surface in the less turbulent boundary layer. With an increase in turbulence in the incoming flow, the strong suction zones decreased in size. The fluctuating pressure coefficient profiles retained their shape when the turbulence in the upstream flow of the cube increased. The fluctuating pressure coefficient was observed to be larger in more turbulent flows. The pressure fluctuations were larger on the cube surface underlying outer boundaries of the conical vortex. The fluctuating pressure coefficient under the conical vortex was three to four times larger than in the weak suction zone on the central area of the top cube surface. Close to the leading cube corner, the pressure spectra were dominated by a single low frequency peak. As the conical vortex developed, this primary peak weakened and a secondary peak emerged at a higher reduced frequency. There is a general trend of shifting the pressure spectra towards higher reduced frequencies when the turbulence in the undisturbed incident flow increases.     

Mass and momentum transport in the near field of swirling turbulent jets. Effect of swirl rate

Local transport of the flow momentum and scalar admixture in the near-field of turbulent swirling jets (Re = 5,000) has been investigated by using a combination of the particle image velocimetry and planar laser-induced fluorescence methods. Advection and turbulent and molecular diffusions are evaluated based on the measured distributions of the mean velocity and concentration and the Reynolds stresses and fluxes. As has been quantified from the data, the flow swirl intensifies the entrainment of the surrounding fluid and promotes mass and momentum exchange in the outer mixing layer. A superimposed swirl results in the appearance of a wake/recirculation region at the jet axis and, consequently, the formation of an inner shear layer. In contrast to the scalar admixture, the momentum exchange in the inner shear layer is found to be strongly intensified by the swirl. For the jet with the highest considered swirl rate, a substantial portion of the surrounding fluid is found to enter the unsteady central recirculation zone, where it mixes with the jet that is issued from the nozzle. The contribution of the coherent velocity fluctuations, which are induced by large-scale vortex structures, to the turbulent transport has been evaluated based on triple decomposition, which was based on proper orthogonal decomposition analysis of the velocity data sets. For the considered domain of the jet with the highest swirl rate and vortex breakdown, the contributions of detected helical vortex structures, inducing pressing vortex core, to the radial fluxes of the flow momentum and the scalar admixture are found to locally exceed 65% and 80%, respectively.     

An analysis of the chaotic motion of particles of different sizes in a gas fluidized bed

The dynamic behavior of individual particles during the mixing/segregation process of particle mixtures in a gas fluidized bed is analyzed. The analysis is based on the results generated from discrete particle simulation, with the focus on the trajectory of and forces acting on individual particles.Typical particles are selected representing three kinds of particle motion:a flotsam particle which is initially at the bottom part of the bed and finally fluidized at the top part of the bed; a jetsam particle which is initially at the top part of the bed and finally stays in the bottom de-fluidized layer of the bed; and a jetsam particle which is intermittently joining the top fluidized and bottom de-fluidized layers. The results show that the motion of a particle is chaotic at macroscopic or global scale, but can be well explained at a microscopic scale in terms of its interaction forces and contact conditions with other particles, particle-fluid interaction force, and local flow structure. They also highlight the need for establishing a suitable method to link the information generated and modeled at different time and length scales.     

Viscous effects in liquid encapsulated liquid bridges

An analytical derivation of the surface deflections and the streamfunctions for the flow inside a liquid encapsulated liquid bridge has been derived using an asymptotic expansion about a small capillary number. The model assumes an initially flat and cylindrical interface under the assumption that the densities of both fluids are equal. To simplify the analysis, the top and bottom walls are assumed to be stress-free and the Reynolds number is assumed to be negligible. Flow is generated either by a moving outer wall (shear-driven flow) or by applying a temperature difference across the top and bottom walls (Marangoni-driven flow). The resulting equations show that for the shear-driven flow, as the viscosity ratio increases, the surface deflections increase monotonically. For the Marangoni-driven flow there exist values of the viscosity ratio where the surface deflections reach a minimum and then switch signs. This investigation shows that it may be possible in more realistic systems to use an outer encapsulating liquid of the proper viscosity ratio to stabilize the liquid–liquid interface during float zone crystal growth.     

Analytical investigation of boundary layer growth and swirl intensity decay rate in a pipe

In this research, the developing turbulent swirling flow in the entrance region of a pipe is investigated analytically by using the boundary layer integral method. The governing equations are integrated through the boundary layer and obtained differential equations are solved with forth-order Adams predictor-corrector method. The general tangential velocity is applied at the inlet region to consider both free and forced vortex velocity profiles. The comparison between present model and available experimental data demonstrates the capability of the model in predicting boundary layer parameters (e.g. boundary layer growth, shear rate and swirl intensity decay rate). Analytical results showed that the free vortex velocity profile can better predict the boundary layer parameters in the entrance region than in the forced one. Also, effects of pressure gradient inside the boundary layer is investigated and showed that if pressure gradient is ignored inside the boundary layer, results deviate greatly from the experimental data.     

An analysis of the chaotic motion of particles of different sizes in a gas fluidized bed

The dynamic behavior of individual particles during the mixing/segregation process of particle mixtures in a gas fluidized bed is analyzed. The analysis is based on the results generated from discrete particle simulation, with the focus on the trajectory of and forces acting on individual particles. Typical particles are selected representing three kinds of particle motion： a flotsam particle which is initially at the bottom part of the bed and finally fluidized at the top part of the bed; a jetsam particle which is initially at the top part of the bed and finally stays in the bottom de-fluidized layer of the bed; and a jetsam particle which is intermittently joining the top fluidized and bottom de-fluidized layers. The results show that the motion of a particle is chaotic at macroscopic or global scale, but can be well explained at a microscopic scale in terms of its interaction forces and contact conditions with other particles, particle-fluid interaction force, and local flow structure. They also highlight the need for establishing a suitable method to link the information generated and modeled at different time and length scales.     

Mixing Layers in Sink Flow: Effect of Length of Flight on Mixing in a Channel Downstream

We have experimentally and analytically studied transport of a passive scalar in the wake of a thin flat plate located at the centerline of a planar contraction with flat walls. The constant Launder parameter in the contraction, K = 6.25 ×10^− 6, was twice the value required for a turbulent boundary layer to relaminarize. In addition to the mixing analysis inside the contraction, layer mixing is also investigated downstream, where the flow continues inside a constant cross-section channel. In order to generate the passive scalar, the airflow above the plate was heated and the temperature stratification in the wake was traced by measuring the temperature field using constant current anemometry. Using different plate lengths, we found that the degree of mixing, obtained at a given position in the straight channel, is a function of the distance from the plate trailing edge to the contraction outlet. For a plate which does not protrude into the straight channel, we demonstrate the existence of an optimal trailing edge-contraction outlet distance that results in the lowest possible degree of mixing at a given downstream position in the straight channel. This finding is also supported by a semi-empirical relationship based on our developed self-similar solution for mixing layers in planar contractions.     

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

Large-eddy simulation of turbulent channel flow at transcritical states

We present well-resolved large-eddy simulations (LES) of a channel flow solving the fully compressible Navier–Stokes equations in conservative form. An adaptive look-up table method is used for thermodynamic and transport properties. A physically consistent subgrid-scale turbulence model is incorporated, that is based on the Adaptive Local Deconvolution Method (ALDM) for implicit LES. The wall temperatures are set to enclose the pseudo-boiling temperature at a supercritical pressure, leading to strong property variations within the channel geometry. The hot wall at the top and the cold wall at the bottom produce asymmetric mean velocity and temperature profiles which result in different momentum and thermal boundary layer thicknesses. Different turbulent Prandtl number formulations and their components are discussed in context of strong property variations.     

Scalar transport from a point source in flows over wavy walls

Simultaneous measurements of the velocity and concentration field in fully developed turbulent flows over a wavy wall are described. The concentration field originates from a low-momentum plume of a passive tracer. PLIF and digital particle image velocimetry are used to make spatially resolved measurements of the structure of the scalar distribution and the velocity. The measurements are performed at three different Reynolds numbers of Re _b = 5,600, Re _b = 11,200 and Re _b = 22,400, respectively, based on the bulk velocity u _b and the total channel height 2h. The velocity field and the scalar field are investigated in a water channel with an aspect ratio of 12:1, where the bottom wall of the test section consists of a train of sinusoidal waves. The wavy wall is characterized by the amplitude to wavelength ratio α = 0.05 and the ratio β between the wave amplitude and the half channel height where β = 0.1. The scalar is released from a point source at the wave crest. For the concentration measurements, Rhodamine B is used as tracer dye. At low to moderate Reynolds number, the flow field is characterized through a recirculation zone which develops after the wave crest. The recirculation zone induces high intensities of the fluctuations of the streamwise velocity and wall-normal velocity. Furthermore, large-scale structures are apparent in the flow field. In previous investigations it has been shown that these large-scale structures meander laterally in flows over wavy bottom walls. The investigations show a strong effect of the wavy bottom wall on the scalar mixing. In the vicinity of the source, the scalar is transported by packets of fluid with a high scalar concentration. As they move downstream, these packets disintegrate into filament-like structures which are subject to strong gradients between the filaments and the surrounding fluid. The lateral scale of the turbulent plume is smaller than the lateral scale of the large-scale structures in the flow field and the plume dispersion is dominated by the structures in the flow field. Due to the lateral meandering of the large-scale structures of the flow field, also the scalar plume meanders laterally. Compared to turbulent plumes in plane channel flows, the wavy bottom wall enhances the mixing effect of the turbulent flow and the spreading rate of the scalar plume is increased.