Numerical Investigation of the Unsteady Flow in a Transonic Compressor with Curved Rotors

The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.

GRIDLESS METHOD FOR UNSTEADY VISCOUS FLOWS

Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.

UNSTEADY FLOWS OF A GENERALIZED SECOND GRADE FLUID WITH THE FRACTIONAL DERIVATIVE MODEL BETWEEN TWO PARALLEL PLATES

The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.

Studies of the unsteady supersonic base flows around three afterbodies

Unsteady supersonic base flows around three afterbodies, cylindrical （Cy）, boattailed （BT） and three-step （MS）, are investigated in this paper. Reynolds-averaged Navier-Stokes （RANS） and two RANS/LES （large-eddy simulation） hybrid methods, detached eddy simulation （DES） and delayed-DES （DDES）, are used to predict the base flow characteristics around the baseline Cy afterbody. All the RANS and hybrid methods are based on the two-equation SST （shear-stress transport） model with compressible corrections （CC）. According to the comparison of measurements, both DES and DDES can produce more satisfactory results than RANS. RANS can only present the ＂stable＂ flow pat- terns, while the hybrid methods can demonstrate unsteady flow structures. DDES and DES results are little different from one another although the latter exhibits better agreement with the experiment. DES is taken to investigate the 5° BT and three-step afterbodies. The mean flow data and the instantaneous turbulent coherent structures are compared against available measurements.

Generation of Dynamic Grids and Computation of Unsteady Transonic Flows around Assemblies

Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.

Unsteady rotating flows of a viscoelastic fluid with the fractional Maxwell model between coaxial cylinders

The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases are solved and the exact solutions are obtained by using the Weber transform and the Laplace transform for fractional calculus.

Unsteady Flow and Structural Behaviors of Centrifugal Pump under Cavitation Conditions

Cavitation has a significant e ect on the flow fields and structural behaviors of a centrifugal pump. In this study, the unsteady flow and structural behaviors of a centrifugal pump are investigated numerically under di erent cavitation conditions. A strong two-way coupling fluid-structure interaction simulation is applied to obtain interior views of the e ects of cavitating bubbles on the flow and structural dynamics of a pump. The renormalization-group k-ε turbulence model and the Zwart–Gerbe–Belamri cavitation model are solved for the fluid side, while a transient structural dynamic analysis is employed for the structure side. The di erent cavitation states are mapped in the head-net positive suction head(H-NPSH) curves and flow field features inside the impeller are fully revealed. Results indicate that cavitating bubbles grow and expand rapidly with decreasing NPSH. In addition, the pressure fluctuations, both in the impeller and volute, are quantitatively analyzed and associated with the cavitation states. It is shown that influence of the cavitation on the flow field is critical, specifically in the super-cavitation state. The e ect of cavitation on the unsteady radial force and blade loads is also discussed. The results indicate that the averaged radial force increased from 8.5 N to 54.4 N in the transition progress from an onset cavitation state to a super-cavitation state. Furthermore, the structural behaviors, including blade deformation, stress, and natural frequencies, corresponding to the cavitation states are discussed. A large volume of cavitating bubbles weakens the fluid forces on the blade and decreases the natural frequencies of the rotor system. This study could enhance the understanding of the e ects of cavitation on pump flow and structural behaviors.

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.

An application of interacting shear flows theory: exact solution for unsteady oblique stagnation point flow

An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS equations and physical analysis and relevant discussions are then presented.

Analysis of Unsteady Flow over Offshore Wind Turbines in Combination with Different types of Foundations

Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.

An SOR Implicit Time-accurate Scheme in Calculating Unsteady Flows

A new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation （SOR） iteration method and Jacobian matrix diagonally dominant splitting method of LUSGS. One advantage of this algorithm is the second-order accuracy because of no factorization error. Another advantage is the low computational cost because the Jacobian matrices and fluxes are only calculated once in each physical time step. And, the SOR algorithm has better convergence property than Gauss-Seidel. To investigate its accuracy and convergency, several unsteady flow computational tests are carded out by using the proposed SOR algorithm. Roe＇s FDS scheme is used to discritize the inviscid flux terms. Unsteady computational results of SOR are compared with the experiment results and those of Gauss-Seidel, Results reveal that the numerical results agree well with the experimental data and the second-order accuracy can be obtained as the Gauss-Seidel for unsteady flow computations. The impact of SOR factor is investigated for unsteady computations by using different SOR factors in this algorithm to simulate each computational test. Different numbers of inner iterations are needed to converge to the same criterion for different SOR factors and optimal choice of SOR factor can improve the computational efficiency greatly.

Reconstruction of the Unsteady Supersonic Flow around a Spike Using the Colored Background Oriented Schlieren Technique

In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Background Oriented Schlieren technique (BOS) allows the measurement of the light deflection caused by density gradients in a compressible flow. For this purpose the distortion of the image of a background pattern observed through the flow is used. In order to increase the performance of the conventional Background Oriented Schlieren technique, the monochromatic background is replaced by a colored dot pattern. The different colors are treated separately using suitable correlation algorithms. Therefore, the precision and the spatial resolution can be highly increased. Furthermore a special arrangement of the different colored dot patterns in the background allows astigmatism in the region with high density gradients to be overcome. For the first time an algebraic reconstruction technique (ART) is then used to reconstruct the density field of unsteady flows around a spike-tipped model from CBOS measurements. The obtained images reveal the interaction between the free-stream flow and the high-pressure region in front of the model, which leads to large-scale instabilities in the flow.

Gridless Solution Method for Two-Dimensional Unsteady Flow

The main purpose of this paper is to develop a gridless method for unsteady flow simulation. A quadrantal point infilling strategy is developed to generate point and combine clouds of points automatically. A point-moving algorithm is introduced to ensure the clouds of points following the movements of bodyboundaries. A dual time method for solving the two-dimenslonal Euler equations in Arbitrary Lagrangian-Eulerian （ALE） formulation is presented. Dual time method allows the real-time step to be chosen on the basis of accuracy rather than stability. It also permits the acceleration techniques, which are commonly used to speed up steady flow calculations, to be used when marching the equations in pseudo time. The spatial derivatives, which are used to estimating the inviscid flux, are directly approximated by using local least-squares curve method. An explicit multistage Runge-Kutta algorithm is used to advance the flow equations in pseudo time. In order to accelerate the solution to convergence, local time stepping technique and residual averaging are employed. The results of NACA0012 airfoil in transonic steady flow are presented to verify the accuracy of the present spatial discretization method. Finally, two AGARD standard test cases in which NACA0012 airfoil and NACA64A010 airfoil oscillate in transonic flow are simulated. The computational results are compared with the experimental data to demonstrate the validity and practicality of the presented method.

A Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows

The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 （2004））. By analyzing the distribution of the finite-time Lyapunov exponent （FTLE） along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.

Momentum and heat transfer of a special case of the unsteady stagnation-point flow

This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy equation is solved with the closed form solutions for prescribed wall temperature and prescribed wall heat flux conditions.The wall temperature and heat flux have power dependence on both time and spatial distance.The solution domain,the velocity distribution,the flow field,and the temperature distribution in the fluids are studied for different controlling parameters.These parameters include the Prandtl number,the mass transfer parameter at the wall,the wall moving parameter,the time power index,and the spatial power index.It is found that two solution branches exist for certain combinations of the controlling parameters for the flow and heat transfer problems.The heat transfer solutions are given by the confluent hypergeometric function of the first kind,which can be simplified into the incomplete gamma functions for special conditions.The wall heat flux and temperature profiles show very complicated variation behaviors.The wall heat flux can have multiple poles under certain given controlling parameters,and the temperature can have significant oscillations with overshoot and negative values in the boundary layers.The relationship between the number of poles in the wall heat flux and the number of zero-crossing points is identified.The difference in the results of the prescribed wall temperature case and the prescribed wall heat flux case is analyzed.Results given in this paper provide a rare closed form analytical solution to the entire unsteady NS equations,which can be used as a benchmark problem for numerical code validation.

Computation of Unsteady Flow Past a Biomimetic Fin

The unsteady hydrodynamics of a biomimetic fin attached to a cylindrical body has been studied numerically using a computational fluid dynamic (CFD) simulator based on an in-house solver of the Navier-Stokes equations, combined with a recently developed multi-block, overset grid method. The fin-body CFD model is based on a mechanical pectoral fin device, which consists of a cylindrical body and an asymmetric fin and can mimic flapping, rowing and feathering motions of the pectoral fins in fishes. First the multi-block, overset grid method incorporated into the NS solver was verified through an extensive study of unsteady flows past a single fin undergoing rowing and feathering motion. Then unsteady flows past the biomimetic fin-body model undergoing the same motions were computed and compared with the measurements of forces of the mechanical pectoral fin, which shows good agreement in both time-varying and time-averaged hydrodynamic forces. The relationship between force generation and vortex dynamics points to the importance of the match in fin kinematics between power and recovery strokes and implies that an optimal selection of parameters of phase lags between and amplitudes of rowing and feathering motions can improve the performance of labriform propulsion in terms of either maximum force generation or minimum mechanical power.

Study of steady and unsteady wet steam condensing flows in a turbine stage

Objective To develop the numerical method for the steady and unsteady wet steam condensing flow in turbine stage. Methods An Eulerian/Eulerian numerical model is used to describe the spontaneous condensation flow in the steam turbine. For the steady condensing flow computations, the mixing plane model was used. For the unsteady condensing flow computations, the sliding mesh method was used to simulate the rotor-stator interactions. Results The numerical results showed the obvious differences between non-condensing and condensing flows. The results also showed the unsteadiness effect due to rotor-stator interactions had a deep influence on the formation and growth process of water droplets. Conclusion The numerical methods presented in this paper are valid for the condensing flow in the turbine stage.

Entropy analysis of unsteady magneto-nanofluid flow past accelerating stretching sheet with convective boundary condition

The unsteady laminar magnetohydrodynamics （MHD） boundary layer flow and heat transfer of nanofluids over an accelerating convectively heated stretching sheet are numerically studied in the presence of a transverse magnetic field with heat source/sink The unsteady governing equations are solved by a shooting method with the Runge-Kutta- Fehlberg scheme. Three different types of water based nanofluids, containing copper, aluminium oxide, and titanium dioxide, are taken into consideration. The effects of the pertinent parameters on the fluid velocity, the temperature, the entropy generation num- ber, the Bejan number, the shear stress, and the heat transfer rate at the sheet surface are graphically and quantitatively discussed in detail. A comparison of the entropy generation due to the heat transfer and the fluid friction is made with the help of the Bejan number. It is observed that the presence of the metallic nanoparticles creates more entropy in the nanofluid flow than in the regular fluid flow.

Unsteady Cavitation Analysis Using Phase Averaging and Conditional Approaches in a 2D Venturi Flow

The present study refers to a cavitating Venturi type section geometry characterized by a convergent angle of 18°?and a divergent angle of about 8°?where the sheet cavity presents typical self-oscillation behavior with quasi-periodic vapor clouds shedding. This work is an extension of previous works concerning void ratio measurements and velocity fields using double optical probe and constitutes a complete analysis of the two-phase structure of unsteady cavitating flow. This paper provides a new method based on conditional and phase averaging technique with wall pressure signal to treat experimental data in order to evaluate more precisely time-averaged and rms values of the void ratio and instantaneous velocity fields. Conditional analysis shows a different behavior of the two-phase flow dynamics leading to highlight high void ratio events linked to the break-off cycle. Unsteady phase averaging of the optical probe signal gives the evolution of the void ratio at each studied location in the venturi and shows that the fluctuations close to the wall (where the re-entrant jet is predominant) are in phase with the upper part of the cavity instead of the thickness of the cavity which is unchanged.