离心泵内流场空化特性的数值模拟研究

针对离心泵内流场特性分析困难的问题,对离心泵流场数值模拟的几何模型建立、模型网格划分和边界条件设定进行了研究,采用计算流体力学方法,获取了在敞水性能条件下离心泵的扬程-流量、效率-流量的变化关系;结合Zwart空化模型,重点对不同有效汽蚀余量时离心泵的空化流场进行了数值模拟,得到了离心泵的内部流线和空泡分布的情况,并与该离心泵机组进行了性能测试实验,最后在此基础之上进行了对比分析。研究结果表明,所采用的数值模拟方法和空化模型合理有效,此结果可为进一步开展离心泵空化监测技术研究提供借鉴。     

Numerical and experimental investigation of cavitation flows in a multistage centrifugal pump

Cavitation behavior is very important in pumps for long time operation. However, there is difficulty in predicting the cavitation phenomena of pumps by Computational fluid dynamics (CFD). In order to accurately ascertain cavitation behavior, a comparison between CFD and experimental data is a significant and essential process. The purpose of this study is to analyze cavitating behavior in multistage centrifugal pumps numerically and experimentally. For this investigation an experimental set up was used to obtain cavitation performance results. The CFD method was used to investigate the multistage centrifugal pump performance under developed cavitating conditions. The Reynolds-averaged Navier-Stokes (RANS) equations were discretized by the finite volume method. The two-equation SST turbulence model was adopted to account for turbulent flows. Numerical data were validated with experimental data and a good comparison of results was achieved. Numerically, cavitation performances were obtained for different pump stages and the effects on cavitation were described according to different NPSH (Net positive suction head). The occurrence of cavitation was also described according to NPSH3% in the head drop lines and water vapor volume fraction on the impeller blade. The rapid drop in head at low NPSH was captured for different flow conditions. It was found that for stage to stage performance, the head drop changes could be related to losses inside the pump. It was also shown that the simulation results can truly represent the development of the attached sheet cavitation in the impeller.     

Critical flow venturi with a step-wall diffuser

The study investigates the critical flow venturi (CVN) with step-wall diffusers for constant volume gas flow rate in critical flow mode. An analytical method for the estimation of the pressure ratio required for the onset of critical flow from these CVN is given. Numerical simulation and experiments have both been employed to compare the onset of critical flow in two types of CVN: ISO 9300 [1] CVN and a CVN with a geometrically similar inlet and a step-wall diffuser with 4 cylindrical sections replacing the conical diffuser. Good agreement between the analytical, numerical and experimental data has been obtained. The CVN with step-wall diffusers have been demonstrated to provide stable flow rate in the critical flow operation mode, however somewhat higher pressure ratio is required to ensure the critical flow in this case.     

Numerical simulation of flows around axisymmetric inlet with bleed regions

A numerical simulation of flows in an axisymmetric supersonic inlet with bleed regions is performed. An existing code which solves the Reynolds Averaged Navier-Stokes equations and the two-equation turbulence model equations is converted into an axisymmetric code. In addition, a bleed boundary condition model has been applied to the code. In this paper, the modified code is validated by comparing numerical results against experimental data and other computational results for flows on a bump and over an oblique shock with bleed region. Using the code, numerical simulation is performed for the flows in an inlet with multiple bleed regions.     

Numerical investigations on swirl intensity decay rate for turbulent swirling flow in a fixed pipe

Due to the importance of predicting the SIDR¹ associated with engineering problems such as combustion chambers, draft tube of the Francis and Kaplan turbines, heat exchanger tubes, separators and so forth, in this research the trend of SIDR and its affecting factors, through a turbulent swirl decay pipe flow have been investigated. The swirling flow is created by means of a rotating honeycomb which produces solid body rotation at the inlet of a fixed pipe. First of all, turbulent swirling decay flow has been numerically surveyed using different flow conditions at the pipe inlet. The numerical results have been validated and compared with the existing experimental data and mathematical relations, showing satisfactory coincide. The obtained results show that, the SIDR depends mainly on the Reynolds number of the passing flow. On this basis, correlations have been proposed in order to improve the predictions of swirl intensity decay rate at upstream regions as well as those with high swirl intensity. In addition, conducted analyses demonstrates (analyses have been made to demonstrate) that the previous developed correlations for predicting swirl intensity decay rate, agree with those provided in this study only for regions far enough from downstream having the low swirl intensity. This implies that the swirl intensity decay rate should be a function of the type of swirl generator at the pipe inlet.     

A multi-tube pressure probe calibration method for measurements of mean flow parameters in swirling flows

A simple and sensitive measurement method of multi-tube pressure probes applicable in swirling flow fields is presented in this paper. Determination of flow direction associated with local pitch, , and yaw, ψ, angles, magnitudes of local static and dynamic pressures can be achieved through a calibration method used with 3- and 5-tube versions of a multi-tube pressure probe. The 5-tube probe was tested in a conventional air cyclone where a strong rotational flow prevailed while the 3-tube probe was used in a low swirl flow field in a pipe internally fitted with a helical coiled wire insert. The method is based on the rotational sensitivity of the pressure probe handled through non-dimensional calibration parameters.     

Transitional behavior of a supersonic flow in a two-dimensional diffuser

Two-dimensional blow-down type supersonic wind tunnel was designed and built to investigate the transient behavior of the startup of a supersonic flow from rest. The contour of the divergent part of the nozzle was determined by the MOC calculation. The converging part of the nozzle, upstream of the throat was contoured to make the flow profile uniform at the throat. The flow characteristics of the steady supersonic condition were visualized using the highspeed schlieren photography. The Mach number was evaluated from the oblique shock wave angle on a sharp wedge with half angle of 5 degree. The measured Mach number was 2.4 and was slightly less than the value predicted by the design calculation. The initial transient behavior of the nozzle was recorded by a high-speed digital video camera with schlieren technique. The measured transition time from standstill to a steady supersonic flow was estimated by analyzing the serial images. Typical transition time was approximately O.1sec.     

Numerical analysis of flow in a Francis turbine on an equal critical cavitation coefficient line

Numerical simulation and model test were applied to study the cavitation flow in Francis turbines. The SST k-ω turbulence model and the mixture model were used to simulate the cavitating flow in the Francis turbine. An equal critical cavitation coefficient line was calculated and the flow in the Francis turbine was analyzed. Simulation results show reasonable agreement with the experimental data. It is confirmed that these cavitation model and numerical method is a useful way to study the two-phase cavitation flow in Francis turbines. On the equal critical cavitation coefficient line, the energy loss in the turbine may be caused by the rotating of vortex rope in the draft tube or flow separation in the runner. The study of equal critical cavitation coefficient line can provide a basic guidance for industry practice.     

An investigation of swirling flow in a cylindrical tube

An experimental study was performed for measuring velocity and turbulence intensity in a circular tube for Re=10,000, 15,000 and 20,000, with swirl and without swirling flow. The velocity fields were measured using PIV techniques and swirl motion was produced by a tangential inlet condition. Some preliminary measurements indicated that over the first 4 diameter, two regions of flow reversal were set up (the so called 2-cell structure). At the highest Reynolds numbers, the maximum values of the measured axial velocity components had moved toward the test tube wall and produce more flow reversal at the corner of the tube. As the Reynolds number increased, the turbulence intensity of swirling flow at the tube inlet also increased.     

超微粉碎过程中气/固两相流场的数值模拟

为认知超微粉碎机粉碎腔内的流场结构和流动状态,采用离散相模型对粉碎机粉碎腔内的气/固两相三维定常流场进行了数值模拟.气相采用重整化群k-ε双方程模型,颗粒相采用随机颗粒轨道模型,同时综合考虑了颗粒受力和湍流扩散对流场状态的影响.模拟结果直观显示了粉碎腔内的速度、总压和流线等特征,研究结果对进行超微粉碎设备关键结构优化设计提供了一定的理论指导.     

Numerical simulation of water flow in an axial flow pump with adjustable guide vanes

A new adjustable guide vane (AGV) is proposed in this paper. This vane can reduce hydraulic losses and improve the performance of an axial flow pump. The formula of AGV adjustment was obtained after theoretical analysis. The fluid flow inside the axial flow pump with a fixed guide vane and adjustable guide vane was simulated. The calculated Q-H curves for the fixed guide vane agreed well with the experimental ones. The results show that the attack angle and flow separation have an important contribution to the vortices which create hydraulic losses in the guide vane channel. The AGV can decrease hydraulic losses and significantly enhance the pump head and efficiency by changing the guide vane angle.     

Numerical simulation of a binary gas flow inside a rotating cylinder

A new approach to calculate the axially symmetric binary gas flow is proposed. Dalton’s law for partial pressures contributed by each species of a binary gas mixture (argon and helium) is incorporated into numerical simulation of rarefied axially symmetric flow inside a rotating cylinder by using the time relaxed Monte-Carlo (TRMC) scheme and the direct simulation Monte-Carlo (DSMC) method. The results of flow simulations are compared with the analytical solution and results obtained by Bird [1]. The results of the flow simulations show better agreement than the results obtained by Bird [1] in comparison with the analytical solutions. However, the results of the flow simulations using the TRMC scheme show better agreement than those obtained using the DSMC method in comparison with the analytical solutions.     

Numerical simulation of multiphase flows with material interface on an unstructured grid system

Two-dimensional multiphase flows with material interface due to density difference are numerically simulated on an unstructured grid system by a Navier-Stokes solver developed by Myong and Kim (2006), since numerical computation for these flows is still known to be difficult, especially if the interface separates fluids of large different densities. This solver employs an unstructured cell-centered method based on a conservative pressure-based finite volume method, since the unstructured grid approach makes the solver very flexible in dealing with complex boundaries, and adopts a high resolution method (CICSAM) in a volume of fluid (VOF) scheme for the accurate phase interface capturing. The test cases are the Rayleigh-Taylor instability (density ratio of 2), the oil bubble rising in a partially filled container (density ratio of 2), the air bubble rising in a fully filled container with bubble shedding (density ratio of 100) and the droplet splash (density ratio of about 1000), which are typical benchmark problems among multiphase flows with material interface due to density difference. The present results are compared with other numerical solutions found in the literature. The present method (solver) efficiently and accurately simulates complex interface flows such as multiphase flows with material interface due to both density difference and instability.     

基于STAR-CCM+的某汽车外流场的数值模拟

以某自主品牌车型为案例,对汽车外流场CFD计算分析过程与结果做了通用性的描述,通过试验和设计数据的对比,得到了比较合理的结果.在此基础上,模拟出汽车外流场尾部的分离流动,从压力场、速度场和流动迹线几个方面分析了模型的数值模拟结果.对汽车虚拟技术平台CFD分析应用流程做了初步的归纳与总结,给相关设计流程以参考性的建议.     

Numerical simulation of fluid flow and heat transfer in a plasma spray gun

A computational fluid dynamics (CFD) simulation for analyzing fluid flow patterns in a plasma spray gun is presented in this study. It is coupled with a heat transfer simulation of the plasma spray gun. Based on CFD and heat transfer theory, the numerical model of the nozzle in the plasma spray gun is developed, and the coupled simulation of the flow fluid and heat transfer is carried out with the semi-implicit method for pressure-linked equations (SIMPLE) method. Local turbulence, which will lead to appearance of a static-water region, is found at the front corner of the cooling channel in the nozzle. The locations insufficiently cooled are found in the wall near the heat source and in the gasket in the rear of the nozzle. Then, cooling processes with different parameters of cooling water are analyzed. The optimal velocity and direction of cooling water, which efficiently cool the nozzle and improve the service life of the plasma jet, are obtained .     

Numerical simulation of fluid flow and heat transfer in a plasma spray gun

A computational fluid dynamics (CFD) simulation for analyzing fluid flow patterns in a plasma spray gun is presented in this study. It is coupled with a heat transfer simulation of the plasma spray gun. Based on CFD and heat transfer theory, the numerical model of the nozzle in the plasma spray gun is developed, and the coupled simulation of the flow fluid and heat transfer is carried out with the semi-implicit method for pressure-linked equations (SIMPLE) method. Local turbulence, which will lead to appearance of a static-water region, is found at the front corner of the cooling channel in the nozzle. The locations insufficiently cooled are found in the wall near the heat source and in the gasket in the rear of the nozzle. Then, cooling processes with different parameters of cooling water are analyzed. The optimal velocity and direction of cooling water, which efficiently cool the nozzle and improve the service life of the plasma jet, are obtained .     

Numerical simulation of flows of non-Newtonian fluids in the stenotic and bifurcated tubes

Steady flows of Newtonian and non-Newtonian fluids in the stenotic and bifurcated tubes are numerically simulated. Four rheologically different fluids such as water, aqueous sugar solution, aqueous Carbopol solution and blood are selected for the namerical simulation and the modified power-law model is used for the numerical simulation of non-Newtonian fluids in the stenotic and bifurcated tubes. Apparent viscosity of a non-Newtonian fluid in the modified power-law model is expressed as a function of the shear rate. Flows in the circular tube with sudden contraction-sudden expansion and gradual contraction-gradual expansion are studied numerically. Analyses in the stenotic tubes are concentrated on the effects of rheological properties, the stenotic geometry and Reynolds number. Flow characteristics of Carbopol solution in the stenotic tubes are compared with those of blood. Effects of the bifurcation geometry on the flow behaviors of Newtonian and non-Newtonian fluids are numerically investigated. Numerical analyses are focused on the flow patterns in the branch tubes of which angles are 30°, 60° and 90° and on the diameter ratios for Newtonian and non-Newtonian fluids. Variations of the axial velocity and pressure drop along the bifurcated tubes for various flow parameters are presented for Newtonian and non-Newtonian fluids.     

水滴迷宫式调节阀空化特性计算与结构改进

针对超(超)临界机组中水滴迷宫式调节阀在高温高压工况下引起的严重气蚀问题,基于计算流体力学理论和空化机理,选用标准k-ε湍流模型、Mixture模型和Schnerr-Sauer空化模型,比较了改进前后调节阀在典型开度下的压力、速度、气相体积分数等结构性能.计算结果表明:原始碟片结构压降大,最大可至19.95 MPa,流...     

Numerical simulation of HVAB rotor in hover using a mixed-mesh flow solver

Numerical simulations of the HVAB were carried out by using a Reynolds-averaged Navier-Stokes computational fluid dynamics flow solver. For the computation, mixed meshes of unstructured/Cartesian grids were used. An improved laminar-turbulent crossflow transition model γ − Re_θt − CF⁺ was used to predict laminar-turbulent transition phenomena. In addition, to achieve high resolute flow solutions, an improved scheme ESWENO-P was employed when calculating inviscid fluxes. To find the turbulence intensity for the simulations, a parametric study was first conducted with the PSP rotor configuration. The influence of the facility walls of the National Full-Scale Aerodynamics Complex was also investigated. Then, the HVAB rotor performance was predicted under conditions obtained by simulations of the PSP rotor. The HVAB rotor performance of the rotor, including thrust/torque coefficients and figure of merit, were analyzed. In addition, flow characteristics, such as laminar-turbulent locations and tip vortex trajectories, were also investigated.

     

微管道中压力驱动液体流动特性的数值模拟研究

MEMS系统的飞速发展让人们对微尺度领域的研究产生了极大的兴趣.对压力驱动下蒸馏水流过直径20μm微管道的流量-压力、流速-压力特性进行了数值模拟研究.结果表明,在该仿真条件下,蒸馏水的流动规律基本符合宏观条件下的Navier-Stokes和Hagen-Poiseuille方程.在仿真研究中还考虑了重力对流动特性的影响,结果发现,在微米尺度下,重力对流动的流量与速度的影响很小,几乎可以忽略不计.