scholarly journals A GPU-Accelerated Compressible RANS Solver for Fluid-Structure Interaction Simulations in Turbomachinery

Author(s):  
L. Mangani ◽  
E. Casartelli ◽  
G. Romanelli ◽  
Magnus Fischer ◽  
A. Gadda ◽  
...  

Computational Fluid Dynamics (CFD) is a fundamental tool for the aerodynamic development in industrial applications. In the usual approach structural deformation due to aerodynamic and thermal loads is often neglected. However, in some cases, where power efficiency is the ultimate goal, an accurate prediction of the structure-flow interaction is essential. This is particularly true for trim and flutter analysis of aircrafts, helicopter and turbomachinery blades. Particularly, turbomachinery trim and flutter predictions still represent a challenge due to phenomena like rotor-stator interaction, separations and shock waves. The usual time-linearised, frequency-domain strategies can be inadequate when this kind of strong non-linear phenomena occur in the flow, making necessary full non-linear time-domain simulations or the harmonic balance technique. Beside flutter, another important aspect, not yet adequately investigated, is the trim analysis, which is fundamental for an accurate steady simulation that aims to consider static blade elasticity for the performance evaluation of turbomachines. Moreover, alongside the obvious contribution given by centrifugal loads to the blade deformation, a not less important source of blade displacement is the thermal effect due to the heat exchanged between the solid and the fluid domains. In particular, for some geometries and operating conditions, thermal effects can be more important than centrifugal effects for the blade deformations. Considering multiple sources of blade deformation (elastic, centrifugal and thermal) in a what is often called “multiphysics” approach is nowadays more and more important, if the goal of the analysis is geometry optimization. To achieve this, next to result’s accuracy also computational efficiency is required, when hundreds of aeroelastic simulations have to be performed in a typical optimization loop. Modern GPUs can be exploited to pursue this goal thanks to their high peak computational power available at relatively low costs and low power consumption with respect to the usual CPUs. In this paper a pioneer work describing the impact of static deformation due to blade elasticity, thermal and centrifugal effects on the performances and power efficiency will be provided. Alongside with accurate results, computational efficiency is taken into account. The purpose of this article is to show the architecture of a GPU-accelerated Fluid-Structure Interaction (FSI) solver for compressible viscous flows. The proposed approach is validated with a typical industrial case, i.e. a turbocharger transonic centrifugal-compressor provided by ABB. The effects of trimmed solutions on the most important integral quantities (i.e. mass flow, characteristic curves, mass-averaged outflow profiles) are investigated and a comparison with pure aerodynamic results is provided. Due to the high blade stiffness and thus the very small displacements obtained with the trim solutions, for the particular case presented in the paper the aeroelastic solutions basically provide nearly the same results as the pure aerodynamic solutions.

Author(s):  
Bhuiyan Shameem Mahmood Ebna Hai ◽  
Markus Bause

Advanced composite materials such as Carbon Fiber Reinforced Plastics (CFRP) are being applied to many aircraft structures in order to improve performance and reduce weight. Most composites have strong, stiff fibers in a matrix which is weaker and less stiff. However, aircraft wings can break due to Fluid-Structure Interaction (FSI) oscillations or material fatigue. This paper focuses on the analysis of a non-linear fluid-structure interaction problem and its solution in the finite element software package DOpElib: the deal.II based optimization library. The principal aim of this research is to explore and understand the behaviour of the fluid-structure interaction during the impact of a deformable material (e.g. an aircraft wing) on air. Here we briefly describe the analysis of incompressible Navier-Stokes and Elastodynamic equations in the arbitrary Lagrangian-Eulerian (ALE) frameworks in order to numerically simulate the FSI effect on a double wedge airfoil. Since analytical solutions are only available in special cases, the equation needs to be solved by numerical methods. This coupled problem is defined in a monolithic framework and fractional-step-θ time stepping scheme are implemented. Spatial discretization is based on a Galerkin finite element scheme. The non-linear system is solved by a Newton method. The implementation using the software library package DOpElib and deal.II serves for the computation of different fluid-structure configurations.


Author(s):  
Prabu Ganesh Ravindren ◽  
Kirti Ghia ◽  
Urmila Ghia

Recent studies of the joined-wing configuration of the High Altitude Long Endurance (HALE) aircraft have been performed by analyzing the aerodynamic and structural behaviors separately. In the present work, a fluid-structure interaction (FSI) analysis is performed, where the fluid pressure on the wing, and the corresponding non-linear structural deformation, are analyzed simultaneously using a finite-element matrix which couples both fluid and structural solution vectors. An unsteady, viscous flow past the high-aspect ratio wing causes it to undergo large deflections, thus changing the domain shape at each time step. The finite element software ANSYS 11.0 is used for the structural analysis and CFX 11.0 is used for the fluid analysis. The structural mesh of the semi-monocoque joined-wing consists of finite elements to model the skin panel, ribs and spars. Appropriate mass and stress distributions are applied across the joined-wing structure [Kaloyanova et al. (2005)], which has been optimized in order to reduce global and local buckling. The fluid region is meshed with very high mesh density at the fluid-structure interface and where flow separation is predicted across the joint of the wing. The FSI module uses a sequentially-coupled finite element equation, where the main coupling matrix utilizes the direction of the normal vector defined for each pair of coincident fluid and structural element faces at the interface [ANSYS 11.0 Documentation]. The k-omega turbulence model captures the fine-scale turbulence effects in the flow. An angle of attack of 12°, at a Mach number of 0.6 [Rangarajan et al. (2003)], is used in the simulation. A 1-way FSI analysis has been performed to verify the proper transfer of loads across the fluid-structure interface. The CFX pressure results on the wing were transferred across the comparatively coarser mesh on the structural surface. A maximum deflection of 16 ft is found at the wing tip with a calculated lift coefficient of 1.35. The results have been compared with the previous study and have proven to be highly accurate. This will be taken as the first step for the 2-way simulation. The effect of a coupled 2-way FSI analysis on the HALE aircraft joined wing configuration will be shown. The structural deformation history will be presented, showing the displacement of the joined-wing, along the wing span over a period of aerodynamic loading. The fluid-structure interface meshing and the convergence at each time step, based on the quantities transferred across the interface will also be discussed.


Author(s):  
Mina Shahi ◽  
Jim B. W. Kok ◽  
P. R. Alemela

The liner of a gas turbine combustor is a very flexible structure that is exposed to the pressure oscillations that occur in the combustor. These pressure oscillations can be of very high amplitude due to thermoacoustic instability, when the fluctuations of the rate of heat release and the acoustic pressure waves amplify each other. The liner structure is a dynamic mechanical system that vibrates at its eigenfrequencies and at the frequencies by which it is forced by the pressure oscillations to which it is exposed. On the other hand the liner vibrations force a displacement of the flue gas near the wall in the combustor. The displacement is very small but this acts like a distributed acoustic source which is proportional to the liner wall acceleration. Hence liner and combustor are a coupled elasto-acoustic system. When this is exposed to a limit cycle oscillation the liner may fail due to fatigue. In this paper the method and the results will be presented of the partitioned simulation of the coupled acousto-elastic system composed of the liner and the flue gas domain in the combustor. The partitioned simulation uses separate solvers for the flow domain and the structural domain, that operate in a coupled way. In this work 2-way fluid structure interaction is studied for the case of a model combustor for the operating conditions 40–60 kW with equivalence ratio of 0.625. This is done in the framework of the LIMOUSINE project. Computational fluid dynamics analysis is performed to obtain the thermal loading of the combustor liner and finite element analysis renders the temperature, stress distribution and deformation in the liner. The software used is ANSYS workbench V13.0 software, in which the information (pressure and displacement) is also exchanged between fluid and structural domain transiently.


Author(s):  
Shuo Yang ◽  
Raymond K. Yee

As a common phenomenon in liquid motions, sloshing usually happens in a partially filled liquid tank of moving vehicle or structure. The objectives of this paper are to study sloshing behavior in rigid tank and deformable tank, and to develop a better performance baffle design in the tank under seismic excitations. The tank is surged with a sinusoidal oscillation about horizontal x-direction. The hydro-elasticity effect of sloshing pressure on the tank wall was taken into consideration due to the fluid-structure interaction between impact pressures and tank structures. ABAQUS finite element program using Coupled Eulerian-Lagrangian (CEL) technique was employed to simulate fluid sloshing. The sloshing phenomenon was studied in rigid tank and deformable tank models with three different water levels, and the effect of wall thickness of the deformable tank on sloshing behavior was discussed. One way to minimize the effect of sloshing in a tank, baffles are used and installed in the middle of the tank, and then various heights and material types of baffle were evaluated. The simulation results show that higher water level case creates greater pressure impact on the tank wall than lower water level case, and the elasticity of the tank structure would reduce the impact pressure of the wall. For the simulation tank model with size of 1m (H) × 1m (W) × 0.2m (D), better performance baffle was found to be the one with the height of 0.35m and was made of acrylic material. Moreover, the conclusion of this study can be extrapolated to other dimensions of the model based on similarity theory. This paper also can serve as an aid in further studying sloshing phenomenon. The findings of this study can be applied to restrain or minimize sloshing motions inside a tank.


Author(s):  
Joost Sterenborg ◽  
Nicola Grasso ◽  
Rogier Schouten ◽  
Arjen Tjallema

Abstract One of the aims of The Ocean Cleanup is to develop technologies to extract plastic pollution from the world’s oceans. Several concepts of passive floating systems were considered that are supposed to confine plastics to ease their collection. Such concepts consist of a floating member and a submerged flexible skirt and have in common that their span is generally more than 500 meters. Consequently, fluid-structure interaction plays an important role in the response of such a floating system. To support numerical simulations, MARIN carried out extensive model tests on a 120 meter system section of the final concept, with focus on the fluid-structure interaction (FSI) of the submerged skirt in operating conditions and in towing configuration. The ability to capture plastics was not investigated in these model tests. Novel for wave-basin tests were non-intrusive measurements using underwater Digital Image Correlation (DIC) to obtain the displacements and deformations of the flexible skirt. DIC proved to be a capable measurement technique for this type of structure in combination with a wave basin. Detailed quantitative data on skirt motions and deformations were delivered and the last concept of the cleanup system was tested in the towing configuration and operational configuration.


Sign in / Sign up

Export Citation Format

Share Document