Numerical simulation and analysis of solid-liquid two-phase threedimensional unsteady flow in centrifugal slurry pump

Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.

A numerical simulation of unsteady flow in small diameter helical grafts

The application of small diameter arterial grafts is limited due to the fact of relatively poor long-time patency which is caused by thrombosis formation in the short term and intimal hyperplasia(IH) in the medium and long term.Thrombosis,obstructing the flow of blood

NUMERICAL SIMULATION OF UNSTEADY FLOW AROUND A CIRCULAR CYLINDER AT HIGH REYNOLDS NUMBERS

A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.

UNSTEADY/STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.

Numerical Simulation of Stall Flow Control Using a DBD Plasma Actuator in Pulse Mode

A numerical simulation method is employed to investigate the effects of the unsteady plasma body force over the stalled NACA 0015 airfoil at low Reynolds number flow conditions. The plasma body force created by a dielectric barrier discharge actuator is modeled with a phenomenological method for plasma simulation coupled with the compressible Navier-Stokes equations. The governing equations are solved using an efficient implicit finitevolume method. The responses of the separated flow field to the effects of an unsteady body force in various inter- pulses and duty cycles as well as different locations and magnitudes are studied. It is shown that the duty cycle and inter-pulse are key parameters for flow separation control. Additionally, it is concluded that the body force is able to attach the flow and can affect boundary layer grow that Mach number 0.1 and Reynolds number of 45000.

Unsteady Flow Simulations in a Three-lobe Positive Displacement Blower

To improve the performance of the positive displacement blower, it is imperative to understand the detailed internal flow characteristics or enable a visualization of flow status. However, the existing two-dimensional unsteady, three-dimensional steady or quasi-unsteady numerical simulation and theoretical analysis cannot provide the detailed flow information, which is unfavorable to improve the performance of positive displacement blower. Therefore, the unsteady flow characteristics in a three-lobe positive displacement blower are numerically investigated by solving the three-dimensional, unsteady, compressible Navier-Stokes equations coupled with RNG k-e turbulent model. In the numerical simulation, the dynamic mesh technique and overset mesh updating method are adopted. Due to the air being compressed in the process of the rotors rotating, the variation of the temperature field in the positive displacement blower is considered. By comparing the experimental measurements and the numerical results on the variation of flow rate with the outlet pressure, the maximum relative error of the flow rate is less than 2.15% even at the maximum outlet pressure condition, which means that the calculation model and numerical computational method used are effective. The numerical results show that in the intake region, the fluctuations of the inlet flow are greatly affected by the direction of the velocity vectors. In the exhaust region, the temperature changes significantly, which leads to the increase of the airflow pulsation. Through analysis on the velocity, pressure and temperature fields obtained from the numerical simulations, three-dimensional unsteady flow characteristics in the positive displacement blower are revealed. The studied results will provide useful reference for improving the performance and empirical correction in the design of the positive displacement blower.

Numerical simulation of a gas pipeline network using computational fluid dynamics simulators

This article describes numerical simulation of gas pipeline network operation using high-accuracy computational fluid dynamics (CFD) simulators of the modes of gas mixture transmission through long, multi-line pipeline systems (CFD-simulator). The approach used in CFD-simulators for modeling gas mixture transmission through long, branched, multi-section pipelines is based on tailoring the full system of fluid dynamics equations to conditions of unsteady, non-isothermal processes of the gas mixture flow. Identification, in a CFD-simulator, of safe parameters for gas transmission through compressor stations amounts to finding the interior points of admissible sets described by systems of nonlinear algebraic equalities and inequalities. Such systems of equalities and inequalities comprise a formal statement of technological, design, operational and other constraints to which operation of the network equipment is subject. To illustrate the practicability of the method of numerical simulation of a gas transmission network, we compare computation results and gas flow parameters measured on-site at the gas transmission enter-prise.

Lagrangian time scales and its relationship to Eulerian equivalents in turbulent channel flow

Lagrangian and Eulerian time scales were obtained from the direct numerical simulation of turbulent channel flow at two Reynolds numbers based on the friction velocity and channel half-height, Rer= 80, 100. The Lagrangian integral time scales and time microscales were compared to their Eulerian equivalents. It is found that the ratio of Lagrangian to TL Eulerian integral time scales is given by TE/TiE= 1 ＋ 0.1y＋ for y＋ ≤ 10, and that the ratios between the Lagrangian to theEulerian time microscales are almost the same irrespective of the components. Those increase with y＋ are approximated by ≈ 2.75 - 1.75 exp （-v＋/a） . These results also show that these expressions are independent of the Reynolds number.

Numerical Simulation of Wake Deflection Control around NACA0012 Airfoil Using Active Morphing Flaps

This study demonstrates an active flow control for deflecting a direction of wake vortex structures behind a NACA0012 airfoil using an active morphing flap. Two-dimensional direct numerical simulations are performed for flows at the chord Reynolds number of 10,000, and the vortex pattern in the controlled and noncontrolled wakes as well as the effect of an actuation frequency on the control ability are rigorously investigated. It is found that there is an optimum actuation-frequency regime at around F + = 2.00 which is normalized by the chord length and freestream velocity. The wake vortex pattern of the well-controlled case is classified as the 2P wake pattern according to the Williamson’s categorization [1] [2], where the forced oscillation frequency corresponds to the natural vortex shedding frequency without control. The present classification of wake vortex patterns and finding of the optimum frequency regime in the wake deflection control can lead to a more robust design suitable for vortex-induced-vibration (VIV) related engineering systems.

Numerical Simulation on Unsteady Flow Mechanism of a 1.5-Stage Axial Transonic Compressor

In this paper,a numerical simulation method is used to calculate a 1.5-stage axial transonic compressor to explore its unsteady flow mechanism.The performance curve is compared with the experimental data to verify the calculation method with a high numerical accuracy,which shows that the unsteady calculation has good reliability.According to the analysis of the data from the monitoring points under the near-stall condition,the unsteady disturbances originate from the tip region of blade and perform the strongest at the blade pressure surface with a broadband characteristic.Further analysis is conducted by combining with the characteristics of the transient flow field at the tip of blade.The results show that the unsteady pressure fluctuations are caused by the migration of the new vortex cores.These new vortex cores are generated by the breakdown of leakage vortex in the downstream,which is induced by the leakage vortex and shock wave interference.Moreover,the relationship between the unsteady flow characteristics and the working conditions is also studied.The leakage vortex intensity and the shock wave strength gradually increase with the decrease of flow rate.When the combination of the leakage vortex intensity and shock wave strength reaches the first threshold,a single frequency of unsteady disturbances appears at the blade tip.When the combination of the leakage vortex intensity and shock wave strength reaches the second threshold,the frequency of unsteady disturbances changes to a broadband.

Unsteady aerodynamic forces and power requirements of a bumblebee in forward flight

Aerodynamic forces and power requirements in forward flight in a bumblebee （Bombus terrestris） were studied using the method of computational fluid dynamics. Actual wing kinematic data of free flight were used in the study （the speed ranges from 0 m/s to 4.5 m/s; advance ratio ranges from 0-0.66）. The bumblebee employs the delayed stall mechanism and the fast pitching-up rotation mechanism to produce vertical force and thrust. The leading-edge vortex does not shed in the translatory phase of the half-strokes and is much more concentrated than that of the fruit fly in a previous study. At hovering and low-speed flight, the vertical force is produced by both the half-strokes and is contributed by wing lift; at medium and high speeds, the vertical force is mainly produced during the downstroke and is contributed by both wing lift and wing drag. At all speeds the thrust is mainly produced in the upstroke and is contributed by wing drag. The power requirement at low to medium speeds is not very different from that of hovering and is relatively large at the highest speed （advance ratio 0.66）, i.e. the power curve is Jshaped. Except at the highest flight speed, storing energy elastically can save power up to 20%-30%. At the highest speed, because of the large increase of aerodynamic torque and the slight decrease of inertial torque （due to the smaller stroke amplitude and stroke frequency used）, the power requirement is dominated by aerodynamic power and the effect of elastic storage of energy on power requirement is limited.

Numerical Investigation of the Aerodynamic Performance Affected by Spiral Inlet and Outlet in a Positive Displacement Blower

The flow in the positive displacement blower is very complex.The existing two-dimensional numerical simulation cannot provide the detailed flow information,especially flow characteristics along the axial direction,which is unfavorable to improve the performance of positive displacement blower.To investigate the effects of spiral inlet and outlet on the aerodynamic performance of positive displacement blower,three-dimensional unsteady flow characteristics in a three-lobe positive displacement blower with and without the spiral inlet and outlet are simulated by solving Navier-Stokes equations coupled with RNG k-ε turbulent model.In the numerical simulation,the dynamic mesh technique and overset mesh updating method are used.The computational results are compared with the experimental measurements on the variation of flow rate with the outlet pressure to verify the validity of the numerical method presented.The results show that the mass flow rate with the change of pressure is slightly affected by the application of spiral inlet and outlet,but the internal flow state is largely affected.In the exhaust region,the fluctuations of pressure,velocity and temperature as well as the average values of velocity are significantly reduced.This illustrates that the spiral outlet can effectively suppress the fluctuations of pressure,thus reducing reflux shock and energy dissipation.In the intake area,the average value of pressure,velocity and temperature are slightly declined,but the fluctuations of them are significantly reduced,indicating that the spiral inlet plays the role in making the flow more stable.The numerical results obtained reveal the three-dimensional flow characteristics of the positive displacement blower with spiral inlet and outlet,and provide useful reference to improve performance and empirical correction in the noise-reduction design of the positive displacement blowers.

Application of Factor Difference Scheme to Solving Discrete Flow Equations Based on Unstructured Grid

A second-order mixing difference scheme with a limiting factor is deduced with the reconstruction gradient method and applied to discretizing the Navier-Stokes equation in an unstructured grid.The transform of nonorthogonal diffusion items generated by the scheme in discrete equations is provided.The Delaunay triangulation method is improved to generate the unstructured grid.The computing program based on the SIMPLE algorithm in an unstructured grid is compiled and used to solve the discrete equations of two types of incompressible viscous flow.The numerical simulation results of the laminar flow driven by lid in cavity and flow behind a cylinder are compared with the theoretical solution and experimental data respectively.In the former case,a good agreement is achieved in the main velocity and drag coefficient curve.In the latter case,the numerical structure and development of vortex under several Reynolds numbers match well with that of the experiment.It is indicated that the factor difference scheme is of higher accuracy,and feasible to be applied to Navier-Stokes equation.

A LARGE EDDY SIMULATION OF THE NEAR WAKE OF A CIRCULAR CYLINDER

Viscous flow around a circular cylinder at a subcritical Reynolds number is investigated using a large eddy simulation (LES) coupled with the Smagorinsky subgrid-scale (SGS) model. A fractional-step method with a second-order in time and a combined finite-difference/spectral approximations are used to solve the filtered three-dimensional incompressible Navier-Stokes equations. Calculations have been performed with and without the SGS model. Turbulence statistical behaviors and flow structures in the near wake of the cylinder are studied. Some calculated results, including the lift and drag coefficients, shedding frequency, peak Reynolds stresses, and time-average velocity profile, are in good agreement with the experimental and computational data, which shows that the Smagorinsky model can reasonably predict the global features of the flow and some turbulent statistical behaviors.

Aerodynamic Coupling Analysis of Counter⁃rotating Propfan and Inlet Flow Field

Taking a propfan engine as the research object,the CFD method was used for 3D modeling and unsteady slip flow for numerical calculation.The propfan rotation domain and the nacelle outside flow domain were meshed by using the partition splicing grid technology.Used the Reynolds⁃averaged of N⁃S equation,the Reynolds stress term uses the RNG turbulence model;and based on the slip grid method,numerical calculation of the flow field with different Mach numbers,front and rear blade angles and engine state were carried out;and the change law of propeller fan characteristics and the influence of slip flow on the inlet flow field were analyzed.The blade angle was the key parameter of the propeller fan characteristic conditions.When the blade angle increases from 41°to 50°,the thrust coefficient increases by 31.2%,and the power coefficient increases by 33.4%;in the climbing state of the propeller fan,the maximum total pressure distortion at the inlet port of 6.8%;the cross section is less affected by the slip flow of the propfan;and the pressure distribution is relatively uniform,but the area of the flow channel is small.The research results can provide a solution for the matching of the counter⁃rotating propeller fan and the engine and the arrangement of the air inlet measuring rake.

超深层裂缝性致密砂岩气藏多尺度耦合流动数值模拟

塔里木盆地克拉苏气田白垩系气藏是罕见的超深层裂缝性致密砂岩气藏,该类气藏的储渗空间具有显著的多尺度特征,基质与多尺度裂缝、断层介质的渗透率级差相差5~6个数量级,常规渗流理论难以准确描述其流动规律和开发机理。为此,基于单裂缝流动物理实验结果及流体力学理论,结合多尺度裂缝几何信息,应用均化理论和体积平均尺度升级方法,将多尺度介质划分为5个流动系统,建立了考虑介质间的非稳态窜流多尺度耦合流动数学模型,并应用有限体积法对耦合流动模型进行了数值求解和数值试井分析。研究结果表明:(1)不同尺度裂缝中具有不同的流动特征,随缝宽增加流速加快,流动模态发生变化;(2)多尺度耦合流动模型与双重介质模型结果存在较大差异,导数曲线具有不同趋势特征;(3)应用所建立的多尺度耦合流动模型成功解释了气藏实际试井数据,模型能够反映实际地层中的流动过程。结论认为,超深层裂缝性致密砂岩气藏多尺度耦合流动模型揭示了多尺度裂缝以及致密基质间逐级动用、协同供气的开发机理,可为类似气藏制订合理开发技术政策及气藏提高采收率提供理论和技术支撑。

压气机内部旋转不稳定的研究综述

旋转不稳定是压气机工作在高负荷近失速工况时的一种常见现象。研究旋转不稳定在降低压气机工作噪声、减小流致振动以及保障航空发动机稳定工作等方面具有重要意义。首先对旋转不稳定现象进行了回顾,详细讨论了旋转不稳定的特征。其次,重点调研了旋转不稳定的起源和机理,将旋转不稳定产生的原因归纳为叶尖泄漏流、涡脱落以及流动剪切等类别。此外,回顾了模拟旋转不稳定的数值方法,讨论了多种流动控制手段对旋转不稳定的作用效果。最后,对旋转不稳定的研究现状进行了总结,对未来的研究趋势进行了展望。

立式轴流泵装置叶轮与导叶间非定常流场分析

为了分析立式轴流泵装置叶轮和导叶体的非定常流场特征,采用ANSYS CFX对立式轴流泵装置全流道进行三维非定常数值计算,预测泵装置的能量性能,并通过物理模型试验进行验证。结果表明,各流量工况时,叶轮进口圆周测线的轴向速度分布均呈类正弦分布,波峰、波谷数与叶轮的叶片数相一致。小流量工况时,叶轮进口面的速度分布受叶轮的影响明显,速度加权偏流角与导叶体出口的平均涡角均最大。随着流量的增大,肘形进水流道的能量损失逐渐增大;导叶体的能量损失随着流量的增大呈先减小后增大的趋势,在最优工况时导叶体的能量损失最小,在小流量工况时60°弯管和虹吸式出水流道的能量损失之和最大。研究结果有利于立式轴流泵装置的结构优化与能效提高。

喷针开度对冲击式水轮机转轮水力性能影响研究

冲击式水轮机是一种适用于中高水头段水力资源开发的流体机械,在我国西南地区高水头水力资源的开发中有着广阔的应用前景。本文以西藏紫霞水电站冲击式水轮机模型机转轮为研究对象,在5个不同的喷针开度下分别进行了转轮内部三维非定常多相流动数值模拟,分析了喷针开度对转轮水力性能与内部流动特性的的影响。预测结果显示转轮水力效率随开度上升而先上升后下降,在设计开度下水力效率最优。在0.75~1.60倍设计流量区间中,转轮水力效率的变化幅度小于1%。根据3个典型开度下转轮水斗内部非定常流动过程的对比分析,发现大开度工况时,水斗缺口处会发生水膜流的漏流现象,喷针开度越大时的泄漏量也越大,从水斗缺口处的漏流现象是其水力损失加剧的主要流动原因。研究为紫霞水电站230 MW大型冲击式水轮转轮水斗的水力优化设计提供了重要参考。

基于非恒定渐变流-急变流方程的变流量畦灌数值模型

传统畦灌模型多是基于非恒定渐变流方程建立的,在模拟变流量畦灌水流运动时的精度难以保障。本文综合分析了变流量畦灌过程中田面水流的运动状况,将其按照边界条件的不同划分为恒定流量进水阶段、变流量进水阶段、畦首消退阶段、田面消退第1阶段、田面消退第2阶段等5个阶段,基于非恒定渐变流方程和非恒定急变流方程构建了适用于变流量畦灌系统的渐变流-急变流数值模型,通过2组恒定流量畦灌、4组变流量畦灌的田间试验以及2组文献资料中的畦灌试验数据对模型进行了验证。结果表明,渐变流-急变流畦灌模型模拟值与现场实测结果吻合较好,模拟推进时间决定系数R^(2)均大于0.96、模拟消退时间R^(2)大于0.90。与目前常用的WinSRFR模型相比,渐变流-急变流畦灌数值模型在模拟恒定流量畦灌方面具有相似的精度,且在模拟变流量畦灌方面精度更高。渐变流-急变流畦灌模型可以较精准地模拟变流量畦灌的水流运动状况,可为分析变流量畦灌系统、优化变流量畦灌方案提供支撑。