核主泵叶轮间隙的轴向振动特性研究

具有径向流的间隙结构广泛存在于轴承结构与旋转机械中,间隙中的流固耦合作用可能影响整体结构的运动稳定性。基于理论间隙模型和核主泵的实际结构,本文对径向间隙流引起的轴向振动进行了多方面的研究。当间隙的一个壁面产生轴向振动并处于旋转状态时,壁面受到由径向流引起的时变轴向力,因此间隙为叶轮提供附加的轴向刚度和阻尼。通过研究以水为介质的理论间隙模型,发现径向间隙流会引起负的等效轴向动力系数（刚度和阻尼）,并且流道形状是影响间隙轴向动力特性的重要因素。扩张流道和平行流道会产生负的轴向动力系数,特别是负阻尼会引起结构振动发散;而收缩流道间隙具有稳定的轴向动力特性。最后,对AP1000核主泵原型叶轮间隙模型进行分析,结果表明,间隙会引起轴向负刚度,并且在一定工况下出现负阻尼,此时系统轴向稳定性及结构安全运行将受到严重的不良影响。     

Modelling of solids distribution in stirred tanks: analysis of simulation strategies and comparison with experimental data

The predictive capabilities of CFD techniques as applied to solid–liquid stirred vessels are investigated. The distribution of solid particles was simulated in three baffled stirred tanks agitated with single and multiple impellers. Suspensions of glass beads of different diameters and average concentration up to 6?vol. % in water were studied. The simulations of solid–liquid suspensions in the stirred vessels were performed by using fully predictive approaches. Eulerian multiphase models were adopted for modelling the solid–liquid flow, coupled with three different extensions of the standard k-? model to the case of multiphase flows. The simulated particle axial concentration profiles are compared with experimental data and critically discussed. The most successful simulation strategy and one possible implementation are described.     

An experimental study of two-phase turbulent coaxial jets

The effect of solid particles on the flow structure of axisymmetric turbulent coaxial jets has been studied. A laser-Doppler anemometer was used to measure the mean and fluctuation velocities of both phases, and a Malvern laser diffraction instrument was applied to measure particle size and concentration. A series of velocity ratios and particle loading ratios were investigated, and the results were analysed for the effects of these ratios on the mixing characteristic and the similarity behavior of the jet. The effects of particle diameter and its distribution were also studied as well as their influence on the coaxial jet behavior.     

Dynamics of swirling jet flows

Experimental investigations of near-field structure of coaxial flows are presented for four different configurations: coaxial jets without rotation (reference case), outer flow rotating only (OFRO), inner-jet rotating only (IJRO) and corotating jets (CRJ). The investigations are performed in a cylindrical water tunnel, with an independent rotation of two coaxial flows. Laser tomography is used to document the flow field, and photographs are shown for different configurations. Time mean velocity profiles obtained by PIV, with and without swirl, are also presented. The dynamics of the swirling jets in the initial region (i.e. near the exit of the jets) is described. The effects of azimuthal velocity and axial velocity ratio variations on flow dynamics are examined. The appearance and growth of the first instabilities are presented and compared with some theoretical results, as is the influence of the rotation (inner or outer) on the dominating structures.     

Computational fluid dynamics simulations of aerosol behavior in a high-speed (Heubach) rotating drum dustiness tester

Potential exposure from hazardous dust may be assessed by evaluating the dustiness of the powders being handled. Dustiness is the tendency of a powder to aerosolize with a given input of energy. Previously we used computational fluid dynamics (CFD) to numerically investigate the flow inside the European Standard (EN15051) rotating drum dustiness tester during its operation. The present work extends those CFD studies to the widely used Heubach rotating drum. Air flow characteristics are investigated within the Abe-Kondoh-Nagano k-epsilon turbulence model; the aerosol is incorporated via a Euler-Lagrangian multiphase approach. The air flow inside these drums consists of a well-defined axial jet penetrating relatively quiescent air. The spreading of the Heubach jet results in a fraction of the jet recirculating as back-flow along the drum walls; at high rotation rates, the axial jet becomes unstable. This flow behavior qualitatively differs from the stable EN15051 flow pattern. The aerodynamic instability promotes efficient mixing within the Heubach drum, resulting in higher particle capture efficiencies for particle sizes d < 80 μm.     

Mixing and segregation of solid particles in a conical spouted bed: Effect of particle size and density

In this work, the mixing and segregation of binary mixtures of particles with different sizes and densities in a pseudo-2D spouted bed were studied experimentally. A binary mixture of solid particles including sand, gypsum, and polyurethane was used. To determine the particles mass fraction, and their mixing and segregation in the bed, an image-processing technique was developed and used. Important hydrodynamic parameters, such as the axial and radial segregation profiles of the solid particles, were measured. The effects of air velocity, particle size, and particle mass fraction were also evaluated. The flow regime in the spouted bed and the time required for reaching the equilibrium state of the solid particles were discussed. The results showed that the segregation of solid particles and the time to equilibrium both decreased when the air velocity increased to much larger than the minimum spouting velocity. The axial segregation increased with the diameter ratio of the particles. Upon completion of the test, coarse particles were concentrated mainly in the spout region, while fine particles were aggregated in the annulus region. Examination of the flow pattern in the spouted bed showed that the particles near the wall had longer flow paths, while those near the spout region had shorter flow paths.     

Effects of particle size and concentration on bubble coalescence and froth formation in a slurry bubble column

A new approach for simulating the formation of a froth layer in a slurry bubble column is proposed. Froth is considered a separate phase, comprised of a mixture of gas, liquid, and solid. The simulation was carried out using commercial flow simulation software (FIRE v2014) for particle sizes of 60–150 μm at solid concentrations of 0–40 vol%, and superficial gas velocities of 0.02–0.034 m/s in a slurry bubble column with a hydraulic diameter of 0.2 m and height of 1.2 m. Modelling calculations were conducted using a Eulerian–Eulerian multiphase approach with k–ε turbulence. The population balance equations for bubble breakup, bubble coalescence rate, and the interfacial exchange of mass and momentum were included in the computational fluid dynamics code by writing subroutines in Fortran to track the number density of different bubble sizes. Flow structure, radial gas holdup, and Sauter mean bubble diameter distributions at different column heights were predicted in the pulp zone, while froth volume fraction and density were predicted in the froth zone. The model was validated using available experimental data, and the predicted and experimental results showed reasonable agreement. To demonstrate the effect of increasing solid concentration on the coalescence rate, a solid-effect multiplier in the coalescence efficiency equation was used. The solid-effect multiplier decreased with increasing slurry concentration, causing an increase in bubble coalescence efficiency. A slight decrease in the coalescence efficiency was also observed owing to increasing particle size, which led to a decrease in Sauter mean bubble diameter. The froth volume fraction increased with solid concentration. These results provide an improved understanding of the dynamics of slurry bubble reactors in the presence of hydrophilic particles.     

Active Control of Swirling Coaxial Jet Mixing with Manipulation of Large-Scale Vortical Structures

Large-scale vortical structures and associated mixing in methane/air swirling coaxial jets are actively controlled by manipulating the outer shear layer of the outer swirling coaxial jet with miniature flap actuators. In order to investigate the control mechanisms, stereoscopic particle image verocimetry (stereo-PIV) and plannar laser-induced fluorescence (PLIF) techniques are employed. It is found that intense vortex rings are produced in the outer shear layer in phase with the periodic flap motion regardless of the swirl number examined. The vortical structures in the inner shear layer, however, are strongly dependent on the swirl rate. This is because the central methane jet is accelerated by the negative axial pressure gradient, of which strength is determined by the swirl. As a result, the inner vortex formation is significantly suppressed at a higher swirl rate. On the other hand, at a relatively low swirl rate, the inner vortices are shed continuously and the methane jet is pinched off. This particular mode promotes the mixing of methane and air, so that the flammable mixture can be formed at an earlier stage of the jet flow development. In addition, the evolution of secondary streamwise vortices is prompted by the combination of the periodic vortex ring shedding and the swirl. They also contribute to the mixing enhancement in the downstream region.     

Numerical simulation of the accumulation of heavy particles in a circular bounded vortex flow

In present work, an Eulerian–Lagrangian CFD model based on the discrete element method (DEM) and immersed boundary method (IBM) has been developed, validated and used to investigate the accumulation of heavy particles in a circular bounded viscous vortex flow. The inter-particle and particle-wall collisions are resolved by a hard-sphere model. Effects of one-way and two-way coupling, Reynolds number, and particle diameter are systematically explored. Results show that, in case of one-way coupling, the majority of particles will spiral into an accumulation point located near the stagnation point of the flow field. The accumulation point represents a stable equilibrium point as the drag created by the flow field balances the destabilizing centrifugal force on the particle. However, in case of two-way coupling, there does not exist a stable accumulation point due to the strong interaction between the particles and fluid dynamics. Instead most particles are expelled from the circular domain and accumulate on the confining wall. The percentage of accumulated particles on the wall increases with increasing Reynolds number and particle diameter. Moreover, influence of three well-known drag models is also studied and they give consistent results on the particle accumulation behavior, although small quantitative differences can still be discerned.     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

回转筒内沙石物料混合均匀性的数值模拟

采用离散单元法对回转筒内沙石物料的混合均匀性进行数值模拟研究,并选取回转筒转速、直径、提升条个数和石子填充率作为影响因素,分别设立3个水平,进行回转筒内物料混合的正交模拟,旨在对各因素的敏感性影响程度进行分析。基于试验结果选择物料混合时间和颗粒接触数率作为试验指标,对四个因素进行极差分析和方差分析。结果表明,回转筒转速、直径和石子填充率对物料混合时间的影响显著,提升条个数对混合时间的影响不显著;石子填充率、回转筒转速和提升条个数对物料混合后颗粒接触数率的影响显著,回转筒直径对物料混合后颗粒接触数率的影响不显著。     

CFD simulation of particle suspension in a stirred tank

Particle suspension characteristics are predicted computationally in a stirred tank driven by a Smith turbine. In order to verify the hydrodynamic model and numerical method, the predicted power number and flow pattern are compared with designed values and simulated results from the literature, respectively. The effects of particle density, particle diameter, liquid viscosity and initial solid loading on particle suspension behavior are investigated by using the Eulerian–Eulerian two-fluid model and the standard k–? turbulence model. The results indicate that solid concentration distribution depends on the flow field in the stirred tank. Higher particle density or larger particle size results in less homogenous distribution of solid particles in the tank. Increasing initial solid loading has an adverse impact on the homogeneous suspension of solid particles in a low-viscosity liquid, whilst more uniform particle distribution is found in a high-viscosity liquid.     

Analysis of minimum specific energy consumption and optimal transport concentration of slurry pipeline transport systems

Slurry pipeline transport is widely used in several industrial processes. Calculating the specific power consumption (SPC) and determining the best working conditions are important for the design and operation of transportation systems. Based on the Shanghai Jiao Tong University high-concentration multi-sized slurry pressure drop (SJTU-HMSPD) pipeline-resistance-calculation model, the SJTU-SPC model for calculating the power required to transport a unit volume of solid materials over a unit pipeline length is established for a slurry transport system. The said system demonstrates a uniformity coefficient in the 1.26–7.98 range, median particle size of 0.075–4 mm, particle volume concentration of 10–60%, and pipeline diameter of 0.203–0.8 m. The results obtained were successfully verified against existing experimental data. The influence of parameters, such as particle-gradation uniformity coefficient, median particle size, pipe diameter, and particle volume concentration, on the SPC were analysed. The results revealed that the greater is the uniformity coefficient, the smaller is the minimum specific energy consumption and the larger the optimal transport concentration for a constant, median particle size slurry. As observed, the optimal transport concentration for broad-graded sand equalled approximately 48%. These results supplement the conclusions of existing research, indicating that the optimal transport concentration is approximately 30% and provides theoretical support for high concentration transportation of broad graded slurry.     

Gas flow in rotary kilns

This paper describes an experimental investigation on the flow characteristics within a rotating cylinder containing a rolling bed of sand. The axis of the cylinder was horizontal and there was no axial bulk flow of particles. The velocity field of the gas flowing through the cylinder was measured by hot-wire anemometry. The measurements indicate that the velocity field is asymmetric with respect to a diameter perpendicular to the granular bed. CFD calculations confirm this finding. The gas velocity profiles are crucial in determining heat transfer from gas to solid.     

Numerical modelling of rotating tangential layers (jets) in shells under strong uniform magnetic field

A numerical study of tangential layers in steady‐state magnetohydrodynamic rotating flows is presented using CFD to solve the inductionless governing equations. The analysis considers two basic flow configurations. In the first, a fluid is enclosed in a cylinder with electrically perfect conducting walls and the flow is driven by a small rotating, conducting disk. In the second, a flow is considered in a spherical shell with an inner rotating sphere. The fluid in both cases is subjected to an external axial uniform magnetic field. The results show that these flows exhibit two different types of flow cores separated from each other by a tangential layer parallel to the axis of rotation. The inner core follows a solid‐body rotation while the outer is quasistagnant. A counter‐rotating jet is developed in the tangential layer between the cores. The characteristics of the tangential layer and the properties of the meridional motion are determined for a wide range of Hartmann numbers. Distributions of angular velocity of circumferential flow and electric potential are obtained and the results are compared with those of analytic methods. Copyright © 2009 John Wiley & Sons, Ltd.     

离心式渣浆泵叶轮磨损规律研究

在实验室条件下对离心式渣浆泵叶轮的磨损规律进行了研究，并对泵轮叶片间流道的颗粒运动轨迹进行数值模拟。探讨了磨粒浓度、泵转速、叶片几何参数、运行时间及磨粒粒径对泵轮磨损的影响。研究表明：叶片的磨损强度与泵轮转速呈５次方关系，叶片进口角对磨损强度影响明显，出口角对磨损强度影响不大，研究结果对渣浆泵的设计具有重要的参考价值。     

Circular Couette flow with pressure-driven axial flow and a porous inner cylinder

In a rotating filter separator a suspension is introduced at one end of the annulus between a rotating porous inner cylinder and a fixed impermeable outer cylinder. The filtrate is removed through the inner cylinder and the concentrate is removed from the opposite end of the annulus from which the suspension entered. The flow in a rotating filter separator is circular Couette flow with a pressure-driven axial flow and a suction boundary condition at the inner cylinder. Flow visualization was used to determine the effect of the Taylor number, axial Reynolds number, and radial Reynolds number on the types of flows present in the annulus. A rich variety of secondary vortical flows appear, depending upon the flow parameters. The radial inflow at the inner cylinder delays the appearance of supercritical circular Couette flow and prevents the appearance of certain flow regimes that have a helical vortex structure. Nevertheless, the average azimuthal velocity measured using laser Doppler velocimetry indicates that the velocity profile is nearly the same for all supercritical flow regimes.This work was supported by a grant from The Whitaker Foundation     

Surfactant-induced effects on turbulent swirling flows

 We employ digital particle imaging velocimetry (DPIV) to investigate the influence of a drag reducing cationic surfactant additive, cetyltrimethyl-ammonium chloride (CTAC), on turbulent swirling flows generated in a cylindrical vessel either by a rotating disk or a rotating disk fitted with vertical flat blades. The largest concentration of CTAC used in this study (0.05 ≤ C ≤ 0.5 mmol/l) is an order of magnitude smaller than those used in experimental investigations of surfactant induced drag reduction in turbulent pipe/channel flows. Even for such dilute systems, a number of dramatic and intriguing effects are observed. In the case of disk-driven flow, it is shown that the surfactant has a non-monotonic influence on turbulence intensity: both radial and axial root mean square velocity fluctuations first increase with increasing surfactant concentration C, reach a maximum and decrease upon further increase in C. Moreover, the maximum intensity is attained at a concentration that is practically independent of the angular frequency Ω of the disk. For the flow driven by bladed impeller, the introduction of the surfactant leads to flow reversal at the impeller plane for low concentrations. Enhancement in the radial and azimuthal mean velocities is also observed. For relatively larger concentrations (=0.5 mmol/l), a mean flow field that consists of multiple transient mixing pockets emerges as Ω exceeds a critical value. Plausible mechanisms are proposed to explain these observations. Received: 11 September 2000 Accepted: 10 April 2001     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics (CFD) modeling of the gas–solids two-phase flow in a circulating fluidized bed (CFB) riser is carried out. The Eularian–Eularian method with the kinetic theory of granular flow is used to solve the gas–solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory (Johnson & Jackson, 1987) with specularity coefficient and particle–wall restitution coefficient. The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operating conditions of the CFB riser.