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.

On the analytical solution of transient friction in channel flows

The transient friction in channel mean flows is the sum of two contributions,i.e.,the underlying laminar flow(ULF)and the purely turbulent component(PTC),and the contributions are analyzed separately by theoretical experiments.It is found that,the transient friction may be higher or remarkably lower than that in equal-Reynolds number steady-state flows.The universal time constant for plane-parallel laminar flows is reported,and the role of the time constant in a turbulent mean flow is examined.It is shown that the time constant is related to the turbulence's frozen time.Finally,a study of the logarithmic layer during the transient flow is accomplished,which shows that the logarithmic layer is destroyed.

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.

Computational Analysis of Centrifugal Pump Delivering Solid-liquid Two-phase Flow during Startup Period

The transient behavior of centrifugal pumps during transient operating periods, such as startup and stopping, has drawn more and more attention recently because of urgent needs in engineering. Up to now, almost all the existing studies on this behavior are limited to using water as working fluid. The study on the transient behavior related to solid-liquid two-phase flow has not been seen yet. In order to explore the transient characteristics of a high specific-speed centrifugal pump during startup period delivering the pure water and solid-liquid two-phase flow, the transient flows inside the pump are numerically simulated using the dynamic mesh method. The variable rotational speed and flow rate with time obtained from experiment are best fitted as the function of time, and are written into computational fluid dynamics （CFD） code-FLUENT by using a user defined function. The predicted heads are compared with experimental results when pumping pure water. The results show that the difference in the transient performance during startup period is very obvious between water and solid-liquid two-phase flow during the later stage of startup process. Moreover, the time for the solid-liquid two-phase flow to achieve a stable condition is longer than that for water. The solid-liquid two-phase flow results in a higher impeller shaft power, a larger dynamic reaction force, a more violent fluctuation in pressure and a reduced stable pressure rise comparing with water. The research may be useful to tmderstanding on the transient behavior of a centrifugal pump under a solid-liquid two-phase flow during startup period.

Numerical Study of Solid-Liquid Two-Phase Flow in StirredTanks with Rushton Impeller(Ⅱ) Prediction of Critical Impeller Speed

Three-dimensional solid-liquid flow is mathematically formulated by means of the 'two-fluid' approach and the two-phase k-ε-AP turbulence model. The turbulent fluctuation correlations appearing in the Reynolds time averaged governing equations are fully incorporated. The solid-liquid flow field and solid concentration distribution in baffled stirred tanks with a standard Rushton impeller are numerically simulated using an improved 'inner-outer'iterative procedure. The flow pattern is identified via the velocity vector plots and a recirculation loop with higher solid concentration is observed in the central vicinity beneath the impeller. Comparison of the simulation with experimental data on the mean velocities and the turbulence quantities of the solid phase is made and quite reasonable agreement is obtained except for the impeller swept volume. The counterpart of liquid phase is presented as well.The predicted solid concentration distribution for three experimental cases with the average solid concentration up to 20% is also found to agree reasonably with the experimental results published in the literature.

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.

Numerical Study of Solid-Liquid Two-Phase Flow in Stirred Tanks with Rushton Impeller (Ⅰ) Formulation and Simulation of Flow Field

The critical impeller speed, N_(JS), for complete suspension of solidparticles in the agitated solid-liquid two-phase system in baffled stirred tanks with a standardRushton impeller is predicted using the computational procedure proposed in Part Ⅰ. Three differentnumerical criteria are tested for determining the critical solid suspension. The predicted N_(JS)is compared with those obtained from several empirical correlations. It is suggested the mostreasonable criterion for determining the complete suspension of solid particles is the positive signof simulated axial velocity of solid phase at the location where the solid particles are mostdifficult to be suspended.

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.