Development of a 3D Navier Stokes Solver for Application to all Types of Turbomachinery

Author(s):  
W. N. Dawes

This paper describes the current stage of development of a code aimed at solving the 3D Navier-Stokes equations in any type of turbomachinery geometry. The basic algorithm time marches the fully 3D unsteady equations of motion expressed in finite volume form with a two step explicit / one step implicit method. Full multigrid acceleration is used to reduce solution time and maintain code performance on fine meshes. Turbulence modelling is via mixing-length closure and the widely used Baldwin-Lomax model. The generality and robustness of the code is demonstrated by application to five different test cases, three axial and two radial configurations. Also included is a grid independence study which demonstrates near grid independent solutions for transonic compressor cascade flow (albeit with the actual result subject to transition modelling constraints). For two of the axial cases (transonic compressor in cascade, secondary flow in a high speed compressor) and one radial case (Eckardt high speed impellor) sufficient mesh is employed for the predictions to be essentially quantitative. The other two cases (radial inflow turbine with clearance and compressor stator with hub clearance) are really simulations rather than predictions, but are included as the flows are novel and provide much physical insight.

1990 ◽  
Author(s):  
P.-A. Chevrin ◽  
C. Vuillez

Accurate prediction of the flow in turbomachinery requires numerical solution of the Navier-Stokes equations. A two-dimensional Navier-Stokes solver developed at ONERA for the calculation of the flow in turbine and compressor cascades was adapted at SNECMA to run on different types of grid. The solver uses an explicit, time-marching, finite-volume technique, with a multigrid acceleration scheme. A multi-domain approach is used to handle difficulties due to the geometry of the flow. An H-C grid was used in the calculations. Two turbulence models, based on the mixing length approach, were used. The flow in a transonic compressor cascade, a subsonic and a transonic turbine cascade were computed. Comparison with experiments is presented.


Author(s):  
W. N. Dawes

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodolgy. The fully three dimensional, compressible, Reynolds averaged Navier-Stokes equations with k-ε turbulence modelling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.


1992 ◽  
Vol 114 (3) ◽  
pp. 528-537 ◽  
Author(s):  
W. N. Dawes

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodology. The fully three-dimensional, compressible, Reynolds-averaged Navier–Stokes equations with k–ε turbulence modeling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.


Author(s):  
Muheng Zhang ◽  
Yongsheng Lian

Coulter counters are analytical microfluidic instrument used to measure the size and concentration of biological cells or colloid particles suspended in electrolyte. The underlying working mechanism of Coulter counters is the Coulter principle which relies on the fact that when low-conductive cells pass through an electric field these cells cause disturbances in the measurement (current or voltage). Useful information about these cells can be obtained by analyzing these disturbances if an accurate correlation between the measured disturbances and cell characteristics. In this paper we use computational fluid dynamics method to investigate this correlation. The flow field is described by solving the Navier-Stokes equations, the electric field is represented by a Laplace’s equation in which the conductivity is calculated from the Navier-Stokes equations, and the cell motion is calculated by solving the equations of motion. The accuracy of the code is validated by comparing with analytical solutions. The study is based on a coplanar Coulter counter with three inlets that consist of two sheath flow inlet and one conductive flow inlet. The effects of diffusivity, cell size, sheath flow rate, and cell geometry are discussed in details. The impacts of electrode size, gap between electrodes and electrode location on the measured distribution are also studied.


Author(s):  
Wolfgang Höhn

During the design of the compressor and turbine stages of today’s aeroengines, aerodynamically induced vibrations become increasingly important since higher blade load and better efficiency are desired. In this paper the development of a method based on the unsteady, compressible Navier-Stokes equations in two dimensions is described in order to study the physics of flutter for unsteady viscous flow around cascaded vibrating blades at stall. The governing equations are solved by a finite difference technique in boundary fitted coordinates. The numerical scheme uses the Advection Upstream Splitting Method to discretize the convective terms and central differences discretizing the viscous terms of the fully non-linear Navier-Stokes equations on a moving H-type mesh. The unsteady governing equations are explicitly and implicitly marched in time in a time-accurate way using a four stage Runge-Kutta scheme on a parallel computer or an implicit scheme of the Beam-Warming type on a single processor. Turbulence is modelled using the Baldwin-Lomax turbulence model. The blade flutter phenomenon is simulated by imposing a harmonic motion on the blade, which consists of harmonic body translation in two directions and a rotation, allowing an interblade phase angle between neighboring blades. Non-reflecting boundary conditions are used for the unsteady analysis at inlet and outlet of the computational domain. The computations are performed on multiple blade passages in order to account for nonlinear effects. A subsonic massively stalled unsteady flow case in a compressor cascade is studied. The results, compared with experiments and the predictions of other researchers, show reasonable agreement for inviscid and viscous flow cases for the investigated flow situations with respect to the Steady and unsteady pressure distribution on the blade in separated flow areas as well as the aeroelastic damping. The results show the applicability of the scheme for stalled flow around cascaded blades. As expected the viscous and inviscid computations show different results in regions where viscous effects are important, i.e. in separated flow areas. In particular, different predictions for inviscid and viscous flow for the aerodynamic damping for the investigated flow cases are found.


Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 168 ◽  
Author(s):  
Agostino Lauria ◽  
Giancarlo Alfonsi ◽  
Ali Tafarojnoruz

Ski jump spillways are frequently implemented to dissipate energy from high-speed flows. The general feature of this structure is to transform the spillway flow into a free jet up to a location where the impact of the jet creates a plunge pool, representing an area for potential erosion phenomena. In the present investigation, several tests with different ski jump bucket angles are executed numerically by means of the OpenFOAM® digital library, taking advantage of the Reynolds-averaged Navier–Stokes equations (RANS) approach. The results are compared to those obtained experimentally by other authors as related to the jet length and shape, obtaining physical insights into the jet characteristics. Particular attention is given to the maximum pressure head at the tailwater. Simple equations are proposed to predict the maximum dynamic pressure head acting on the tailwater, as dependent upon the Froude number, and the maximum pressure head on the bucket. Results of this study provide useful suggestions for the design of ski jump spillways in dam construction.


2020 ◽  
Vol 2020 (6) ◽  
Author(s):  
Alexander Braginsky

Abstract In this paper, we study the vortex motion of a continuous medium, which is described by forces obtained from the principle of least action. It is shown that in a continuous medium the vortex force components are proportional to the velocity and pressure gradient components. This article gives a description of the 2D vortex motion of air in zones of high and low pressure. If the pressure decreases, the angular velocity of rotation of the continuous medium increases, whereas if the pressure increases, the angular velocity fades. The lifting force is obtained due to the vortex movement of air in the form of a funnel. It is shown that the vortex force contains a vortex term of the Euler hydrodynamic equations with a relative factor equal to the velocity of the continuous medium squared divided by the sound velocity squared. To describe the motion of a continuous medium correctly it is necessary to replace the forces obtained by Euler with the forces obtained from the minimum of action in the equations of motion. It is concluded that vortex motions and turbulence are described by the obtained equations of motion, and not by the Navier–Stokes equations. Most likely, this is related to the Problem of the Millennium description of turbulence announced at the International Congress of Mathematics in 2000.


Author(s):  
K M Guleren ◽  
A Pinarbasi

The main goal of the present work is to analyse the numerical simulation of a centrifugal pump by solving Navier-Stokes equations, coupled with the ‘standard k-∊’ turbulence model. The pump consists of an impeller having five curved blades with nine diffuser vanes. The shaft rotates at 890r/min. Flow characteristics are assumed to be stalled in the appropriate region of flowrate levels of 1.31-2.861/s. Numerical analysis techniques are performed on a commercial FLUENT package program assuming steady, incompressible flow conditions with decreasing flowrate. Under stall conditions the flow in the diffuser passage alternates between outward jetting when the low-pass-filtered pressure is high to a reverse flow when the filtered pressure is low. Being below design conditions, there is a consistent high-speed leakage flow in the gap between the impeller and the diffuser from the exit side of the diffuser to the beginning of the volute. Separation of this leakage flow from the diffuser vane causes the onset of stall. As the flowrate decreases both the magnitude of the leakage within the vaneless part of the pump and reverse flow within a stalled diffuser passage increase. As this occurs, the stall-cell size extends from one to two diffuser passages. Comparisons are made with experimental data and show good agreement.


Author(s):  
Y Wang ◽  
S Komori

A pressure-based finite volume procedure developed previously for incompressible flows is extended to predict the three-dimensional compressible flow within a centrifugal impeller. In this procedure, the general curvilinear coordinate system is used and the collocated grid arrangement is adopted. Mass-averaging is used to close the instantaneous Navier-Stokes equations. The covariant velocity components are used as the main variables for the momentum equations, making the pressure-velocity coupling easier. The procedure is successfully applied to predict various compressible flows from subsonic to supersonic. With the aid of the k-ɛ turbulence model, the flow details within a centrifugal impeller are obtained using the present procedure. Predicted distributions of the meridional velocity and the static pressure are reasonable. Calculated radial velocities and flow angles are favourably compared with the measurements at the exit of the impeller.


Sign in / Sign up

Export Citation Format

Share Document