Comparison of Different Acceleration Techniques and Methods for Periodic Boundary Treatment in Unsteady Turbine Stage Flow Simulations

1999 ◽  
Vol 122 (2) ◽  
pp. 234-246 ◽  
Author(s):  
Martin von Hoyningen-Huene ◽  
Alexander R. Jung

This paper studies different acceleration techniques for unsteady flow calculations. The results are compared with a nonaccelerated, fully explicit solution in terms of time-averaged pressure distributions, the unsteady pressure and entropy in the frequency domain, and the skin friction factor. The numerical method solves the unsteady three-dimensional Navier–Stokes equations via an explicit time-stepping procedure. The flow in the first stage of a modern industrial gas turbine is chosen as a test case. After a description of the numerical method used for the simulation, the test case is introduced. The purpose of the comparison of the different numerical algorithms for explicit schemes is to facilitate the decision as to which acceleration technique should be used for calculations with regard to accuracy and computational time. The convergence acceleration methods under consideration are explicit time-stepping with implicit residual averaging, explicit time-consistent multigrid, and implicit dual time stepping. The investigation and comparison of the different acceleration techniques apply to all explicit unsteady flow solvers. This paper also examines the influence of the stage blade count ratio on the flowfield. For this purpose, a simulation with a stage pitch ratio of unity is compared with a calculation using the real ratio of 78:80, which requires a more sophisticated method for periodic boundary condition treatment. This paper should help to decide whether it is crucial from the turbine designer’s point of view to model the real pitch ratio in unsteady flow simulations in turbine stages. [S0889-504X(00)00702-9]

Author(s):  
Martin von Hoyningen-Huene ◽  
Alexander R. Jung

This paper studies different acceleration techniques for unsteady flow calculations. The results are compared with a non-accelerated, fully-explicit solution in terms of time-averaged pressure distributions, the unsteady pressure and entropy in the frequency domain and the skin friction factor. The numerical method solves the unsteady three-dimensional Navier-Stokes equations via an explicit time-stepping procedure. The flow in the first stage of a modern industrial gas turbine is chosen as a test case. After a description of the numerical method used for the simulation, the test case is introduced. The comparison of the different numerical algorithms for explicit schemes is intended to ease the decision about which acceleration technique to use for calculations as far as accuracy and computational time are concerned. The convergence acceleration methods under consideration are, respectively, explicit time-stepping with implicit residual averaging, explicit time-consistent multigrid and implicit dual time stepping. The investigation and comparison of the different acceleration techniques are applicable to all explicit unsteady flow solvers. As another point of interest, the influence of the stage blade count ratio on the flow field is investigated. For this purpose, a simulation with a stage pitch ratio of unity is compared with a calculation using the real ratio of 78:80, which requires a more sophisticated method for periodic boundary condition treatment. This paper should help to decide whether it is vital from the turbine designer’s point of view to model the real pitch ratio in unsteady flow simulations in turbine stages.


Author(s):  
L. Eça ◽  
G. Vaz ◽  
M. Hoekstra ◽  
S. Pal ◽  
E. Muller ◽  
...  

Abstract Two workshops were held at the ASME V&V Symposiums of 2017 and 2018 dedicated to Iterative Errors in Unsteady Flow Simulations. The focus was on the effect of iterative errors on numerical simulations performed with implicit time integration, which require the solution of a nonlinear set of equations at each time-step. The main goal of these workshops was to create awareness to the problem and to confirm that different flow solvers exhibited the same trends. The test case was a simple two-dimensional, laminar flow of a single-phase, incompressible, Newtonian fluid around a circular cylinder at the Reynolds number of 100. A set of geometrically similar multiblock structured grids was available and boundary conditions to perform the simulations were proposed to the participants. Results from seven flow solvers were submitted, but not all of them followed exactly the proposed conditions. One set of results was obtained with adaptive grid and time refinement using triangular elements (CADYF) and another used a compressible flow solver with a dual time stepping technique and a Mach number of 0.2 (DLR-Tau). The remaining five submissions were obtained with five different incompressible flow solvers (ansyscfx 14.5, pimplefoam, refresco, saturne, starccm+ v12.06.010-R8) using implicit time integration in the proposed grids. The results obtained in this simple test case showed that iterative errors may have a significant impact on the numerical accuracy of unsteady flow simulations performed with implicit time integration. Iterative errors can be significantly larger (one to two orders of magnitude) than the residuals and/or solution changes used as convergence criteria at each time-step. The Courant number affected the magnitude of the iterative errors obtained in the proposed exercise. For the same iterative convergence criteria at each time-step, increasing the Courant number tends to increase the iterative error.


2003 ◽  
Vol 9 (5) ◽  
pp. 353-361 ◽  
Author(s):  
Eun Seok Lee ◽  
George S. Dulikravich ◽  
Brian H. Dennis

An axial turbine rotor cascade-shape optimization with unsteady passing wakes was performed to obtain an improved aerodynamic performance using an unsteady flow, Reynolds-averaged Navier-Stokes equations solver that was based on explicit, finite difference; Runge-Kutta multistage time marching; and the diagonalized alternating direction implicit scheme. The code utilized Baldwin-Lomax algebraic andk-εturbulence modeling. The full approximation storage multigrid method and preconditioning were implemented as iterative convergence-acceleration techniques. An implicit dual-time stepping method was incorporated in order to simulate the unsteady flow fields. The objective function was defined as minimization of total pressure loss and maximization of lift, while the mass flow rate was fixed during the optimization. The design variables were several geometric parameters characterizing airfoil leading edge, camber, stagger angle, and inter-row spacing. The genetic algorithm was used as an optimizer, and the penalty method was introduced for combining the constraints with the objective function. Each individual's objective function was computed simultaneously by using a 32-processor distributedmemory computer. The optimization results indicated that only minor improvements are possible in unsteady rotor/stator aerodynamics by varying these geometric parameters.


1994 ◽  
Vol 29 (1-2) ◽  
pp. 53-61
Author(s):  
Ben Chie Yen

Urban drainage models utilize hydraulics of different levels. Developing or selecting a model appropriate to a particular project is not an easy task. Not knowing the hydraulic principles and numerical techniques used in an existing model, users often misuse and abuse the model. Hydraulically, the use of the Saint-Venant equations is not always necessary. In many cases the kinematic wave equation is inadequate because of the backwater effect, whereas in designing sewers, often Manning's formula is adequate. The flow travel time provides a guide in selecting the computational time step At, which in turn, together with flow unsteadiness, helps in the selection of steady or unsteady flow routing. Often the noninertia model is the appropriate model for unsteady flow routing, whereas delivery curves are very useful for stepwise steady nonuniform flow routing and for determination of channel capacity.


Author(s):  
J. C. Jaeger

The object of this note is to indicate a numerical method for finding periodic solutions of a number of important problems in conduction of heat in which the boundary conditions are periodic in the time and may be linear or non-linear. One example is that of a circular cylinder which is heated by friction along the generators through a rotating arc of its circumference, the remainder of the surface being kept at constant temperature; here the boundary conditions are linear but mixed. Another example, which will be discussed in detail below, is that of the variation of the surface temperature of the moon during a lunation; in this case the boundary condition is non-linear. In all cases the thermal properties of the solid will be assumed to be independent of temperature. Only the semi-infinite solid will be considered here, though the method applies equally well to other cases.


Sign in / Sign up

Export Citation Format

Share Document