scholarly journals Adaptive Unstructured 2D Navier-Stokes Solutions on Mixed Quadrilateral/Triangular Meshes

Author(s):  
Stuart D. Connell ◽  
D. Graham Holmes ◽  
Mark E. Braaten

This paper presents a solution adaptive scheme for solving the Navier-Stokes equations on an unstructured mixed grid of triangles and quadrilaterals. The solution procedure uses an explicit Runge-Kutta finite volume time marching scheme with an adaptive blend of second and fourth order smoothing. The governing equations are solved in a 2D, axisymmetric or quasi-3D form. In viscous regions quadrilateral elements are used to facilitate the one dimensional refinement required for the efficient resolution of boundary layers and wakes. The effect of turbulence is incorporated through using either a Baldwin-Lomax or k-ε turbulence model. Solutions are presented for several examples that illustrate the capability of the algorithm to predict viscous phenomena accurately. The examples are a transonic turbine, a nozzle and a combustor diffuser.

1988 ◽  
Vol 110 (3) ◽  
pp. 339-346 ◽  
Author(s):  
O. K. Kwon

A robust, time-marching Navier–Stokes solution procedure based on the explicit hopscotch method is presented for solution of steady, two-dimensional, transonic turbine cascade flows. The method is applied to the strong conservation form of the unsteady Navier–Stokes equations written in arbitrary curvilinear coordinates. Cascade flow solutions are obtained on an orthogonal, body-conforming “O” grid with the standard k–ε turbulence model. Computed results are presented and compared with experimental data.


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.


1987 ◽  
Author(s):  
Chen Nai-Xing ◽  
Dai Li-Hong

In this paper, a time marching method using a hopscotch algorithm and a stream-surface-fitted co-ordinate system for calculating steady viscous flow on the S2 stream surface is presented. It is convenient to express the Navier-Stokes equations by the non-orthogonal curvilinear co-ordinate system directly on the S2 stream surface because the blade force term which exists in the governing equations written on the meridional plane, disappears. Numerical results for the CAS single rotor research compressor of Institute of Engineering Thermophysics are compared with invisid calculation and experiment. It is shown that the computational results are agreed with experiment well.


1997 ◽  
Vol 52 (4) ◽  
pp. 358-368 ◽  
Author(s):  
Michio Nishida ◽  
Masashi Matsumotob

Abstract • This paper describes a computational study of the thermal and chemical nonequilibrium occuring in a rapidly expanding flow of high-temperature air transported as a free jet from an orifice into low-density stationary air. Translational, rotational, vibrational and electron temperatures are treated separately, and in particular the vibrational temperatures are individually treated; a multi-vibrational temperature model is adopted. The governing equations are axisymmetric Navier-Stokes equations coupled with species vibrational energy, electron energy and species mass conservation equations. These equations are numerically solved, using the second order upwind TVD scheme of the Harten-Yee type. The calculations were carried out for two different orifice temperatures and also two different orifice diameters to investigate the effects of such parameters on the structure of a nonequilibrium free jet.


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.


Author(s):  
T. Tanuma ◽  
N. Shibukawa ◽  
S. Yamamoto

An implicit time-marching higher-order accurate finite-difference method for solving the two-dimensional compressible Navier-Stokes equations was applied to the numerical analyses of steady and unsteady, subsonic and transonic viscous flows through gas turbine cascades with trailing edge coolant ejection. Annular cascade tests were carried out to verify the accuracy of the present analysis. The unsteady aerodynamic mechanisms associated with the interaction between the trailing edge vortices and shock waves and the effect of coolant ejection were evaluated with the present analysis.


1994 ◽  
Vol 271 ◽  
pp. 1-16 ◽  
Author(s):  
Peter Y. Huang ◽  
Jimmy Feng ◽  
Daniel D. Joseph

We do a direct two-dimensional finite-elment simulation of the Navier–Stokes equations and compute the forces which turn an ellipse settling in a vertical channel of viscous fluid in a regime in which the ellipse oscillates under the action of vortex shedding. Turning this way and that is induced by large and unequal values of negative pressure at the rear separation points which are here identified with the two points on the back face where the shear stress vanishes. The main restoring mechanism which turns the broadside of the ellipse perpendicular to the fall is the high pressure at the ‘stagnation point’ on the front face, as in potential flow, which is here identified with the one point on the front face where the shear stress vanishes.


It is shown that the boundary layer approximation to the flow of a viscous fluid past a flat plate of length l , generally valid near the plate when the Reynolds number Re is large, fails within a distance O( lRe -3/4 ) of the trailing edge. The appropriate governing equations in this neighbourhood are the full Navier- Stokes equations. On the basis of Imai (1966) these equations are linearized with respect to a uniform shear and are then completely solved by means of a Wiener-Hopf integral equation. The solution so obtained joins smoothly on to that of the boundary layer for a flat plate upstream of the trailing edge and for a wake downstream of the trailing edge. The contribution to the drag coefficient is found to be O ( Re -3/4 ) and the multiplicative constant is explicitly worked out for the linearized equations.


Sign in / Sign up

Export Citation Format

Share Document