Assessment of the Academic CFD Code Galatea-I With the DLR-F11 Model in High Lift Configuration

Author(s):  
Sotirios S. Sarakinos ◽  
Georgios N. Lygidakis ◽  
Ioannis K. Nikolos

In this study the development and assessment of an academic CFD (Computational Fluid Dynamics) code, named Galatea-I, is reported. The proposed solver employs the RANS (Reynolds-Averaged Navier-Stokes) approach, modified by the artificial compressibility method, along with the SST (Shear Stress Transport) turbulence model to predict steady or unsteady turbulent incompressible flow phenomena on three-dimensional unstructured hybrid grids, composed of prismatic, tetrahedral and pyramidal elements. Parallel processing and an agglomeration multigrid method have been included for the acceleration of the solver’s methodologies. Galatea-I is evaluated against a test case of the HiLiftPW-2 (Second High Lift Prediction Workshop). In particular, the low Mach number flow at 7° incidence angle over the DLR-F11 aircraft configuration of Case 1 of the aforementioned workshop was examined; it considers a three-element wing with a leading edge slat and a trailing edge flap attached on a body pod, without including though any of the support brackets used in the wind tunnel experiments. The obtained results are close to the available experimental data, as well as the numerical results of other reference solvers, indicating the proposed methodology’s potential to predict accurately such low Mach number flows over complex geometries.

Author(s):  
Ouafa Bouloumou ◽  
Eric Serre ◽  
Jochen Fröhlich

A three-dimensional spectral method based on Chebyshev-Chebyshev-Fourier discretizations is presented in the framework of the low Mach number approximation of Navier-Stokes equations. The working fluid is assumed to be a perfect gas with constant thermodynamic properties. The generalized Stokes problem, which arises from the time discretization by a second-order semi-implicit scheme, is solved by a preconditioned iterative Uzawa algorithm. Several validation results are presented in the case of steady and unsteady flows. This model is also evaluated for natural convection flows with large density variations in the case of a tall differentially heated cavity of aspect ratio 8. It is found that on contrary to convection at small temperature differences (Boussinesq), the 2D unsteady solution at Ra = 3.4 x 105 is unstable to 3D perturbations.


Author(s):  
Sotirios S. Sarakinos ◽  
Georgios N. Lygidakis ◽  
Ioannis K. Nikolos

In this study an academic CFD code, named Galatea-I, is presented, capable for simulating inviscid, viscous laminar and viscous turbulent incompressible fluid flows. It employs the RANS (Reynolds-Averaged Navier-Stokes) approach along with the SST (Shear Stress Transport) turbulence model to predict turbulent flow phenomena, such as recirculations and separations of flow, on three-dimensional unstructured hybrid grids, composed of prismatic, tetrahedral and pyramidal elements. Discretization of the governing equations is obtained with a node-centered finite-volume scheme. Parallel processing and agglomeration multigrid scheme are implemented for the acceleration of the numerical process. As the title of this paper reveals, the solver is validated against the test cases of the DARPA SUBOFF program; in particular, flows over the SUBOFF bare hull submarine geometry at two incident angles and the SUBOFF hull with fairwater configuration are examined. The obtained results, compared to available in open literature experimental data as well as results computed by reference solvers, indicate the proposed methodology’s potential to accurately simulate complex fluid flows.


1993 ◽  
Vol 22 (4-5) ◽  
pp. 441-451 ◽  
Author(s):  
Jörn Sesterhenn ◽  
Bernhard Müller ◽  
Hans Thomann

Author(s):  
David Maltese ◽  
Antonín Novotný

Abstract We investigate the error between any discrete solution of the implicit marker-and-cell (MAC) numerical scheme for compressible Navier–Stokes equations in the low Mach number regime and an exact strong solution of the incompressible Navier–Stokes equations. The main tool is the relative energy method suggested on the continuous level in Feireisl et al. (2012, Relative entropies, suitable weak solutions, and weak–strong uniqueness for the compressible Navier–Stokes system. J. Math. Fluid Mech., 14, 717–730). Our approach highlights the fact that numerical and mathematical analyses are not two separate fields of mathematics. The result is achieved essentially by exploiting in detail the synergy of analytical and numerical methods. We get an unconditional error estimate in terms of explicitly determined positive powers of the space–time discretization parameters and Mach number in the case of well-prepared initial data and in terms of the boundedness of the error if the initial data are ill prepared. The multiplicative constant in the error estimate depends on a suitable norm of the strong solution but it is independent of the numerical solution itself (and of course, on the discretization parameters and the Mach number). This is the first proof that the MAC scheme is unconditionally and uniformly asymptotically stable in the low Mach number regime.


Author(s):  
Wei Li ◽  
Hua Ouyang ◽  
Zhao-hui Du

To give insight into the clocking effect and its influence on the wake transportation and its interaction, the unsteady three-dimensional flow through a 1.5-stage axial low pressure turbine is simulated numerically using a density-correction based, Reynolds-Averaged Navier-Stokes equations commercial CFD code. The 2nd stator clocking is applied over ten equal tangential positions. The results show that the harmonic blade number ratio is an important factor affecting the clocking effect. The clocking effect has a very small influence on the turbine efficiency in this investigation. The efficiency difference between the maximum and minimum configuration is nearly 0.1%. The maximum efficiency can be achieved when the 1st stator wake enters the 2nd stator passage near blade suction surface and its adjacent wake passes through the 2nd stator passage close to blade pressure surface. The minimum efficiency appears if the 1st stator wake impinges upon the leading edge of the 2nd stator and its adjacent wake of the 1st stator passed through the mid-channel in the 2nd stator.


Author(s):  
Vijay K. Garg

A multi-block, three-dimensional Navier-Stokes code has been used to compute heat transfer coefficient on the blade, hub and shroud for a rotating high-pressure turbine blade with 172 film-cooling holes in eight rows. Film cooling effectiveness is also computed on the adiabatic blade. Wilcox’s k-ω model is used for modeling the turbulence. Of the eight rows of holes, three are staggered on the shower-head with compound-angled holes. With so many holes on the blade it was somewhat of a challenge to get a good quality grid on and around the blade and in the tip clearance region. The final multi-block grid consists of 4784 elementary blocks which were merged into 276 super blocks. The viscous grid has over 2.2 million cells. Each hole exit, in its true oval shape, has 80 cells within it so that coolant velocity, temperature, k and ω distributions can be specified at these hole exits. It is found that for the given parameters, heat transfer coefficient on the cooled, isothermal blade is highest in the leading edge region and in the tip region. Also, the effectiveness over the cooled, adiabatic blade is the lowest in these regions. Results for an uncooled blade are also shown, providing a direct comparison with those for the cooled blade. Also, the heat transfer coefficient is much higher on the shroud as compared to that on the hub for both the cooled and the uncooled cases.


Author(s):  
Weiliang Fu ◽  
Jie Gao ◽  
Chen Liang ◽  
Fukai Wang ◽  
Qun Zheng ◽  
...  

The flow in high endwall-angle turbine is complex, and it is different from the ordinary turbine flow in characteristics. In order to study the flow field characteristics of high endwall-angle turbines, the annular sector cascade experimental study of high endwall-angle turbines is carried out. The blade is studied experimentally in the form of annular sector cascade. The cascade includes 7 blades, and makes up 6 flow passages, in order to simulate full cascade flow. The experimental Mach number is adjusted by the way of changing inlet total pressure, and the Mach number influence (0.7, 0.8 and 0.9) on annular sector cascade flow is studied. Based on it, the inlet incidence angle (−15°, −7.5°, 0°, 7.5° and 15° )is changed with the way of changing sector straight pipes upstream of the cascade, and its influence on turbine flow fields is studied at the Mach number of 0.8. Here, five-hole probes are used to measure aerodynamic parameters distributions downstream of the cascade, and static pressure taps are positioned on the blade surface to measure surface static pressure distribution. The auto-traversing system and pressure sensors were operated by a self-compiled program based control program. The results indicate that there are two passage vortices inside the turbine cascade flow passage under the high Mach number condition, and the passage vortex near the high endwall-angle region is bigger. As Mach number increases, the passage vortices inside turbine cascade passage will become strong, and moves towards the blade mid-span. Besides, it is shown that the way of changing sector straight pipes can achieve the variation of inlet incidence angles. And, the blade profile with big leading-edge radius has good design and off-design performance. Detailed results and analyses are presented in the paper.


Sign in / Sign up

Export Citation Format

Share Document