电弧等离子体一维非定常数值模拟

通过对毛细管内等离子体工作过程的分析给出了电弧等离子体的一维非定常数学模型 ,以热力学参量描述了非理想条件下等离子体的状态方程和电导率计算方法。利用差分法求解该数学方程组 ,给出了毛细管内电弧等离子体特征参量在毛细管内沿轴向随时间的变化规律。对数值模拟结果作相应的分析 ,并将部分模拟结果与实验数据作了比较。

非定常喷嘴射流形状预测研究

在对射流进行非定常CFD解析的基础上,提出了射流直径的计算公式和射流形状的统计预测方法,并由时间平均的射流形状确定了缩流位置、膨胀率等用于设计的重要参数。

直驱风力发电机翼型特性数值仿真

目前的发电机翼型特性数值仿真方法中,由于对翼身的抖振非定常数的计算不足,影响最终结果的准确性,因此提出直驱风力发电机特性数值仿真方法研究。首先建立风力空气动力模型,作为仿真动力基础。其次确定风翼参数,计算翼身的抖振非定常数,最后在数值仿真中将翼型实现实体化生成,并根据风力空气动力模型实现数值仿真。为了验证设计的数值仿真方法的可行性,利用某地直驱风力发电机的相关参数,使用不同对比数值仿真方法和设计方法进行翼型特性数值仿真。实验结果证明设计方法的特性系数更接近实际值,在仿真中准确性较高。

Pressure fluctuation and its influencing factors in circulating water pump

In order to investigate the effect of sampling frequency and time on pressure fluctuations, the three-dimensional unsteady numerical simulations were conducted in a circulating water pump. Through comparison of turbulence models with hydraulic performance experiment, SST k-co model was confirmed to study the rational determination of sampling frequency and time better. The Fast Fourier Transform （FFT） technology was then adopted to process those fluctuating pressure signals obtained. On these bases, the characteristics of pressure fluctuations acting on the tongue were discussed. It is found that aliasing errors decrease at higher sampling frequency of 17 640 Hz, but not at a lower sampling frequency of 1 764 Hz. Correspondingly, an output frequency range ten-times wider is obtained at 17 640 Hz. Compared with 8R, when the sampling time is shorter, the amplitudes may be overvalued, and the frequencies and amplitudes of low-frequency fluctuations can not be well predicted. The frequencies at the tongue are in good agreement with the values calculated by formula and the frequency compositions less than the blade passing frequency are accurately predicted.

Unsteady Numerical Simulation of Flow around 2-D Circular Cylinder for High Reynolds Numbers

In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54×106. The calculations were performed by means of solving the 2-D Unsteady Reynolds-Averaged Navier-Stokes （URANS） equations with a k-ε turbulence model. The calculated results, produced flow structure drag and lift coefficients, as well as Strouhal numbers. The findings were in good agreement with previous published data, which also supplied us with a good understanding of the flow across cylinders of different high Reynolds numbers. Meanwhile, an effective measure was presented to control the lift force on a cylinder, which points the way to decrease the vortex induced vibration of marine structure in future.

Experimental Study on Calculation of Hydro-Geological Parameters for Unsteady Flow

After the linear analytical method of unsteady flow theory is further improved,an innovative and faster algorithm is introduced.The water storage in a confined aquifer is derived from the water transmissivity coefficient and the water-pressure conductivity coefficient.The water transmissivity coefficient is approximated by a Taylorseries expansion of drawdown,and the water-pressure conductivity coefficient is obtained by the average drawdown.In this algorithm,the distance of the observation points from the pumping well must be short.When the distance is as short as the radius of the main pumping well,the data of the drawdown difference between the sidewall and the center of pumping well are difficult to measure,but the same results can be achieved based on the assumption that the drawdown difference approximates to the drawdown of the observation wells at a radial distance from the pumping well according to the algorithm.Without the help of charts,this algorithm is more concise and efficient,which has been verified by the test of water pumping project in Tianjin Binhai International Airport.

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.

Simulation of Hydrodynamic Performance of Drag and Double Reverse Propeller Podded Propulsors

The unsteady performance of drag and double reverse propeller podded propulsors in open water was numerically simulated using a computational fluid dynamics （CFD） method. A moving mesh method was used to more realistically simulate propulsor working conditions, and the thrust, torque, and lateral force coefficients of both propulsors were compared and analyzed. Forces acting on different parts of the propulsors along with the flow field distribution of steady and unsteady results at different advance coefficients were compared. Moreover, the change of the lateral force and the difference between the abovementioned two methods were mainly analyzed. It was shown that the thrust and torque results of both methods were similar, with the lateral force results having the highest deviation

Numerical simulation of the influence of ground effect on the performance of multi section wings

the establishment of multi-element airfoil in steady and unsteady ground effect N-S equation turbulence model, the S-A model of multi element airfoils during takeoff and landing high attack angle change numerical simulation analysis, the calculation results show that the lower altitude, lift and drag wing angle decreased; the greater the ground the effect is more obvious, the greater the loss of lift. The simulation results show that the lift coefficient is slightly less than that of unsteady numerical simulation, and the drag coefficient is slightly less than that of unsteady numerical simulation. The ground disturbance to the wing not only affects the steady state flow field, but also is closely related to the unsteady aerodynamic performance. The results of this study can provide a reference for the design and flight control of large aircraft wings.

Decomposition of the Unsteady Wave Patterns for Bessho form Translating-Pulsating Source Green Function

In order to interpret the physical feature of Bessho form translating-pulsating source Green function, the phase function is extracted from the integral representation and stationary-phase analysis is carried out in this paper. The complex characteristics of the integral variable and segmentation of the integral intervals are discussed in m complex plane. In θ space, the interval [-π/2+φ,-π/2+φ-iε] is dominant in the near-field flow, and there is a one-to-one correspondence between the real intervals in m space and the unsteady wave patterns in far field. If 4τ>1(τ is the Brard number), there are three kinds of propagation wave patterns such as ring-fan wave pattern, fan wave pattern and inner V wave pattern, and if 0<4τ<1, a ring wave pattern, an outer V and inner V wave pattern are presented in far field. The ring-fan or ring wave pattern corresponds to the interval [-π+α,-π/2+φ] for integral terms about k2, and the fan or outer V wave pattern and inner V wave pattern correspond to [-π+α,-π/2) and(-π/2,-π/2+φ] respectively for terms about k1. Numerical result shows that it is beneficial to decompose the unsteady wave patterns under the condition of τ≠0 by converting the integral variable θ to m. In addition, the constant-phase curve equations are derived when the source is performing only pulsating or translating.

A CFD Based Investigation of the Unsteady Hydrodynamic Coefficients of 3-D Fins in Viscous Flow

The motion of the fins and control surfaces of underwater vehicles in a fluid is an interesting and challenging research subject.Typically the effect of fin oscillations on the fluid flow around such a body is highly unsteady, generating vortices and requiring detailed analysis of fluid-structure interactions.An understanding of the complexities of such flows is of interest to engineers developing vehicles capable of high dynamic performance in their propulsion and maneuvering.In the present study, a CFD based RANS simulation of a 3-D fin body moving in a viscous fluid was developed.It investigated hydrodynamic performance by evaluating the hydrodynamic coefficients (lift, drag and moment) at two different oscillating frequencies.A parametric analysis of the factors that affect the hydrodynamic performance of the fin body was done, along with a comparison of results from experiments.The results of the simulation were found in close agreement with experimental results and this validated the simulation as an effective tool for evaluation of the unsteady hydrodynamic coefficients of 3-D fins.This work can be further be used for analysis of the stability and maneuverability of fin actuated underwater vehicles.

Analysis of unsteady supercavitating flow around a wedge

Supercavitating flow around a slender symmetric wedge moving at variable velocity in static fluid has been studied. Singular integral equation for the flow has been founded through distributing the sources and sinks on the symmetrical axis. The supereavity length at each moment is determined by solving the singular integral equation with finite difference method. The supercavity shape at each moment is obtained by solving the partial differential equation with variable coefficient. For the case that the wedge takes the impulse and uniformly variable motion, numerical results of time history of the supercavity length and shape are presented. The calculated results indicate that the shape and the length of the supercavity vary in a similar way to the case that the wedge takes variable motion, and there is a time lag in unsteady supercavitating flow induced by the variation of wedge velocity.

Numerical Tools for the Control of the Unsteady Heating of an Airfoil

This paper concerns the real time control of the boundary layer on an aircraft wing. This new approach consists in heating the surface in an unsteady regime using electrically resistant strips embedded in the wing skin. The control of the boundary layer＇s separation and transition point will provide a reduction in friction drag, and hence a reduction in fuel consumption. This new method consists in applying the required thermal power in the different strips in order to ensure the desired temperatures on the aircraft wing. We also have to determine the optimum size of these strips （length, width and distance between two strips）. This implies finding the best mathematical model corresponding to the physics enabling us to facilitate the calculation for any type of material used for the wings. Secondly, the heating being unsteady, and, as during a flight the flow conditions or the ambient temperatures vary, the thermal power needed changes and must be chosen as fast as possible in order to ensure optimal operating conditions.

Statistical analysis of pressure fluctuations during unsteady flow for low-specific-speed centrifugal pumps

A three-dimensional transient numerical simulation was conducted to study the pressure fluctuations in low-specific-speed centrifugal pumps. The characteristics of the inner flow were investigated using the SST k-ω turbulence model. The distributions of pressure fluctuations in the impeller and the volute were recorded, and the pressure fluctuation intensity was analyzed comprehensively, at the design condition, using statistical methods. The results show that the pressure fluctuation intensity increases along the impeller streamline from the leading edge to the trailing edge. In the impeller passage, the intensity near the shroud is much higher than that near the hub at thc inlet. However, the intensity at the middle passage is almost equal to the intensity at the outlet. The pressure fluctuation intensity is the highest at the trailing edge on the pressure side and near the tongue because of the rotor-stator interaction. The distribution of pressure fluctuation intensity is symmetrical in the axial cross sections of the volute channel. However, this intensity decreases with increasing radial distance. Hence, the pressure fluctuation intensity can be reduced by modifying the geometry of the leading edge in the impeller and the tongue in the volute.

Numerical Study of Unsteady Behavior of Partial Cavitation on Two Dimensional Hydrofoils

This paper deals with time dependent performance characteristics of cavitating hydrofoils, the flow around which has been simulated using pressure-based finite volume method. A bubble dynamics cavitation model was used to investigate the unsteady behavior of cavitating flow and describe the generation and evaporation of vapor phase. For choosing the turbulence model and mesh size a non cavitating study was conducted. Three turbulence models such as Spalart-Allmaras, Shear Stress Turbulence （SST） κ-ω model, Re-Normalization Group （RNG） κ-ε model with enhanced wall treatment are used to capture the turbulent boundary layer along the hydrofoil surface. The cavitating study presents an unsteady behavior of the partial cavity attached to the foil at different time steps for σ = 0.8 and σ = 0.4. Moreover, this study is focused on cavitation inception, the shape and general behavior of sheet cavitation, lift and drag forces for different cavitation numbers.

Application of 2D Numerical Model to Unsteady Performance Evaluation of Vertical-Axis Tidal Current Turbine

Abstract Tidal current energy is renewable and sustainable, which is a promising altemative energy resource for the future elec- tricity supply. The straight-bladed vertical-axis turbine is regarded as a useful tool to capture the tidal current energy especially under low-speed conditions. A 2D unsteady numerical model based on Ansys-Fluent 12.0 is established to conduct the numerical simulation, which is validated by the corresponding experimental data. For the unsteady calculations, the SST model, 2x 105 and 0.01 s are se- lected as the proper turbulence model, mesh number, and time step, respectively. Detailed contours of the velocity distributions around the rotor blade foils have been provided for a flow field analysis. The tip speed ratio （TSR） determines the azimuth angle of the appearance of the torque peak, which occurs once for a blade in a single revolution. It is also found that simply increasing the incident flow velocity could not improve the turbine performance accordingly. The peaks of the averaged power and torque coeffi- cients appear at TSRs of 2.1 and 1.8, respectively. Furthermore, several shapes of the duct augmentation are proposed to improve the turbine performance by contracting the flow path gradually from the open mouth of the duct to the rotor. The duct augmentation can significantly enhance the power and torque output. Furthermore, the elliptic shape enables the best performance of the turbine. The numerical results prove the capability of the present 2D model for the unsteady hydrodynamics and an operating performance analy- sis of the vertical tidal stream turbine.

Temporal Behaviour of a Corner Separation in a Radial Vaned Diffuser of a Centrifugal Compressor Operating near Surge

The temporal behaviour of a flow separation in the hub-suction side comer of a transonic diffuser is studied thanks to unsteady numerical simulations based on the phase-lagged approach. The validity of the numerical re- sults is confn＇med by comparison with experimental unsteady pressure measurements. An analysis of the instan- taneous skin-friction pattern and particles trajectories is presented. It highlights the topology of the separation and its temporal behaviour. The major result is that, despite of a highly time-dependent core flow, the separation is found to be a ＂fixed unsteady separation＂ characterized by a fixed location of the main saddle of the separation but an extent of the stall region modulated by the pressure waves induced by the impeller-diffuser interaction.

Investigation on Multiple Cylindrical Holes Casing Treatment for Transonic Axial Compressor Stability Enhancement

Casing treatment is one possible way of regaining axial compressor operating range. However, most of casing treatments extend the operating range with the cost of efficiency penalty. A new form of multiple cylindrical holes casing treatment(MHCT) with pre-swirl blowing for the NASA Rotor-37 has been designed based on profound understanding of the stall inception. Unsteady numerical simulations have been performed for Rotor-37 with and without MHCT. Parametric studies of the total extraction holes area and their axial locations show that the compressor performance deteriorates as the area ratio increases but the stall margin is extended and there is an optimum extraction holes axial location for stall margin extending. The better configuration of MHCT could extend the stall margin by 6.2% with only 0.23% peak efficiency reduction. Detailed analysis of the physical mechanism behind the stall margin improvement shows that the casing treatment could eliminate the passage blockage by suppressing breakup of tip leakage vortex and decrease the blade load in tip region, which both contribute to improve stall margin of transonic axial compressors.

Numerical simulation of unsteady cavitating flows around a transient pitching hydrofoil

The objective of this paper is to improve the understanding of the influence of multiphase flow on the turbulent closure model, the interplay between vorticity fields and cavity dynamics around a pitching hydrofoil. The effects of pitching rate on the sub- cavitating and cavitating response of the pitching hydrofoil are also investigated. In particular, we focus on the interactions between cavity inception, growth, and shedding and the vortex flow structures, and their impacts on the hydrofoil performance. The calculations are 2-D and performed by solving the incompressible, multiphase Unsteady Reynolds Averaged Navier Stokes （URANS） equations via the commercial CFD code CFX. The k-co SST （Shear Stress Transport） turbulence model is used along with the transport equation-based cavitation models. The density correction function is considered to reduce the eddy viscosity according to the computed local fluid mixture density. The calculation results are validated with experiments conducted by Ducoin et al. （see Computational and experimental investigation of flow over a transient pitching hydrofoil, Eur J Mech/B Fluids, 2009, 28：728-743 and An experimental analysis of fluid structure interaction of a flexible hydrofoil in vari- ous flow regimes including cavitating flow, Eur J Mech B/fluids, 2012, 36： 63-74）. Results are shown for a NACA66 hydro- foil subject to slow （quasi static, t2=6~/s, ＆＊ =0.18） and fast （dynamic, ＆=63~/s, dr＂ =1.89） pitching motions from a =0~ to a =15~. Both subcavitaing （or =8.0） and cavitating （cr=3.0） flows are considered. For subcavitating flow （or=8.0）, low frequency fluctuations have been observed when the leading edge vortex shedding occurs during stall, and delay of stall is ob- served with increasing pitching velocity. For cavitating flow （tr=3.0）, small leading edge cavities are observed with the slow pitching case, which significantly modified the vortex dynamics at high angles of attack, leading to high frequency fluctuations of the hydrodynamic coefficients and different stall behaviors compared to the subcavitating flow at the same pitching rate. On the other hand, for the fast pitching case at or=3.0, large-scale sheet/cloud cavitation is observed, the cavity behavior is un- steady and has a strong impact on the hydrodynamic response, which leads to high amplitude fluctuations of the hydrodynamic coefficients, as well as significant changes in the stall and post-stall behavior. The numerical results also show that the local density modification helps to reduce turbulent eddy viscosity in the cavitating region, which significantly modifies the cavity lengths and shedding frequencies, particularly for the fast pitching case. In general, compared with the experimental visualiza- tions, the numerical results with local density correction have been found to agree well with experimental measurements and observations for both slow and fast transient pitching cases.

Suction Performance and Internal Flow of a 2-Bladed Helical Inducer with Inlet Asymmetric Plate

It has been found in our past studies that the installation of asymmetric plate at the inlet of inducer is effective for the suppression of cavitation surge phenomenon.In the present study,the suction performance of 2-bladed helical inducer with an inlet asymmetric plate is experimentally investigated.It is observed that the suction performance in large flow rate conditions is not significantly influenced by the asymmetric plate,whereas the head of inducer with the asymmetric plate increases just before the head breakdown in partial flow conditions.To understand the mechanism of this additional head,the flow measurements and the numerical simulations are carried out.It is found that the circumferential component of absolute velocity at the exit of inducer slightly increases with the development of cavitation in both cases with and without the inlet asymmetric plate,indicating the increase of the theoretical head.The theoretical head increase with the inlet asymmetric plate is also confirmed by the unsteady numerical simulations,suggesting that the additional head is achieved through the increase of the theoretical head with the change of the exiting flow from the inducer associated with some amount of cavitation.