Assessment of a 3-D boundary layer code to predict heat transfer and flow losses in a turbine

The prediction of the complete flow field in a turbine passage is an extremely difficult task due to the complex three dimensional pattern which contains separation and attachment lines, a saddle point and horseshoe vortex. Whereas, in principle such a problem can be solved using full Navier-Stokes equations, in reality methods based on a Navier-Stokes solution procedure encounter difficulty in accurately predicting surface quantities (e.g., heat transfer) due to grid limitations imposed by the speed and size of the existing computers. On the other hand the overall problem is strongly three dimensional and too complex to be analyzed by the current design methods based on inviscid and/or viscous strip theories. Thus there is a strong need for enhancing the current prediction techniques through inclusion of 3-D viscous effects. A potentially simple and cost effective way to achieve this is to use a prediction method based on three dimensional boundary layer (3-DBL) theory. The major objective of this program is to assess the applicability of such a 3-DBL approach for the prediction of heat loads, boundary layer growth, pressure losses and streamline skewing in critical areas of a turbine passage. A brief discussion of the physical problem addressed here along with the overall approach is presented

Analysis of airfoil leading edge separation bubbles

A local inviscid-viscous interaction technique was developed for the analysis of low speed airfoil leading edge transitional separation bubbles. In this analysis an inverse boundary layer finite difference analysis is solved iteratively with a Cauchy integral representation of the inviscid flow which is assumed to be a linear perturbation to a known global viscous airfoil analysis. Favorable comparisons with data indicate the overall validity of the present localized interaction approach. In addition numerical tests were performed to test the sensitivity of the computed results to the mesh size, limits on the Cauchy integral, and the location of the transition region

Development of a three-dimensional Navier-Stokes code on CDC star-100 computer

A three-dimensional code in body-fitted coordinates was developed using MacCormack's algorithm. The code is structured to be compatible with any general configuration, provided that the metric coefficients for the transformation are available. The governing equations are developed in primitive variables in order to facilitate the incorporation of physical boundary conditions and turbulence-closure models. MacCormack's two-step, unsplit, time-marching algorithm is used to solve the unsteady Navier-Stokes equations until steady-state solution is achieved. Cases discussed include (1) flat plate in supersonic free stream; (2) supersonic flow along an axial corner; (3) subsonic flow in an axial corner at M infinity = 0.95; and (4) supersonic flow in an axial corner at M infinity 1.5

ALESEP: A computer program for the analysis of airfoil leading edge separation bubbles

The ALESEP program for the analysis of the inviscid/viscous interaction which occurs due to the presence of a closed laminar transitional separation bubble on an airflow is presented. The ALESEP code provides a iterative solution of the boundary layer equations expressed in an inverse formulation coupled to a Cauchy integral representation of the inviscid flow. This interaction analysis is treated as a local perturbation to a known solution obtained from a global airfoil analysis. Part of the required input to the ALESEP code are the reference displacement thickness and tangential velocity distributions. Special windward differencing may be used in the reversed flow regions of the separation bubble to accurately account for the flow direction in the discretization of the streamwise convection of momentum. The ALESEP code contains a forced transition model based on a streamwise intermittency function and a natural transition model based on a solution of the integral form of the turbulent kinetic energy equation. Instructions for the input/output, and program usage are presented

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic

The use of Levy-Lees variables in three-dimensional boundary-layer flows

A method for solving a general class of three-dimensional boundary layer flows is developed. In the development, Levy-Lees variables are extended to three dimensions and equations are placed in these similarity variables. An implicit finite difference scheme which is stable for negative transverse velocities is used to solve these equations. The method developed is applied to obtain solutions for sharp and spherically blunted circular cones at angle of attack. Longitudinal and transverse distributions are presented for these cases. Good agreement is found with the results obtained by other numerical schemes and the experimental data of Tracy, for sharp circular cones at angle of attack. For spherically blunted cones at angle of attack, the results are in good agreement with axisymmetric sphere results up to the region where spherical symmetry holds

A multistage time-stepping scheme for the thin-layer Navier-Stokes equations

A finite-volume scheme for numerical integration of the Euler equations was extended to allow solution of the thin-layer Navier-Stokes equations in two and three dimensions. The extended algorithm, which is based on a class of four-stage Runge-Kutta time-stepping schemes, was made numerically efficient through the following convergence acceleration technique: (1) local time stepping, (2) enthalpy damping, and (3) residual smoothing. Also, the high degree of vectorization possible with the algorithm has yielded an efficient program for vector processors. The scheme was evaluated by solving laminar and turbulent flows. Numerical results have compared well with either theoretical or other numerical solutions and/or experimental data

Computational Analysis of Dual Radius Circulation Control Airfoils

The goal of the work is to use multiple codes and multiple configurations to provide an assessment of the capability of RANS solvers to predict circulation control dual radius airfoil performance and also to identify key issues associated with the computational predictions of these configurations that can result in discrepancies in the predicted solutions. Solutions were obtained for the Georgia Tech Research Institute (GTRI) dual radius circulation control airfoil and the General Aviation Circulation Control (GACC) dual radius airfoil. For the GTRI-DR airfoil, two-dimensional structured and unstructured grid computations predicted the experimental trend in sectional lift variation with blowing coefficient very well. Good code to code comparisons between the chordwise surface pressure coefficients and the solution streamtraces also indicated that the detailed flow characteristics were matched between the computations. For the GACC-DR airfoil, two-dimensional structured and unstructured grid computations predicted the sectional lift and chordwise pressure distributions accurately at the no blowing condition. However at a moderate blowing coefficient, although the code to code variation was small, the differences between the computations and experiment were significant. Computations were made to investigate the sensitivity of the sectional lift and pressure distributions to some of the experimental and computational parameters, but none of these could entirely account for the differences in the experimental and computational results. Thus, CFD may indeed be adequate as a prediction tool for dual radius CC flows, but limited and difficult to obtain two-dimensional experimental data prevents a confident assessment at this time

Summary of the 2004 CFD Validation Workshop on Synthetic Jets and Turbulent Separation Control

A computational fluid dynamics (CFD) validation workshop for synthetic jets and turbulent separation control (CFDVAL2004) was held in Williamsburg, Virginia in March 2004. Three cases were investigated: synthetic jet into quiescent air, synthetic jet into a turbulent boundary layer crossflow, and flow over a hump model with no-flow-control, steady suction, and oscillatory control. This paper is a summary of the CFD results from the workshop. Although some detailed results are shown, mostly a broad viewpoint is taken, and the CFD state-of-the-art for predicting these types of flows is evaluated from a general point of view. Overall, for synthetic jets, CFD can only qualitatively predict the flow physics, but there is some uncertainty regarding how to best model the unsteady boundary conditions from the experiment consistently. As a result, there is wide variation among CFD results. For the hump flow, CFD as a whole is capable of predicting many of the particulars of this flow provided that tunnel blockage is accounted for, but the length of the separated region compared to experimental results is consistently overpredicted

Stratégie de grille conforme octrée intersectées pour les Applications aux calculs Aéroacoustiques de LAGOON, Modèle de train d'Atterrissage, utilisant le Flow Solver CEDRE non structuré

International audienceAircraft noise is a societal concern and landing gears contribute significantly to the generated noise in approach and landing configurations. Landing gears are characterized by their complex geometry and numerous works have been carried out to develop and validate aeroacoustics simulations to predict the associated noise. Most of them associate a time resolved flow solution, to capture the acoustic sources, to an acoustic computation, to estimate the resulting far field noise. Due to the geometric complexity, unstructured grids are required and may necessitate meticulous work to optimize. In this context, Lattice Boltzmann Methods (LBM) have become popular as they propose to combine automatic grid generation and high CPU efficiency and produced remarked results. The automatic grid generation is facilitated by the use of advanced wall models that do not require resolution of complex details of boundary layer flow, ranging from attached to detached regimes, that are produced by the complex geometries and flow environment of landing gears. Navier-Stokes (NS) solvers on the contrary rely on precise boundary layer solution that require complex grids, even in the unstructured approach, to handle the attached boundary layer regimes, that require strong grid anisotropy, as well as detached regimes and their trailing flow, that require grid isotropy. The grid construction work can therefore become a complex process. The simplification of this process is then an important challenge for industrial applications. The present work details a multi-year effort at ONERA in that direction