Achieving High Speed CFD simulations: Optimization, Parallelization, and FPGA Acceleration for the unstructured DLR TAU Code

Today, large scale parallel simulations are fundamental tools to handle complex problems. The number of processors in current computation platforms has been recently increased and therefore it is necessary to optimize the application performance and to enhance the scalability of massively-parallel systems. In addition, new heterogeneous architectures, combining conventional processors with specific hardware, like FPGAs, to accelerate the most time consuming functions are considered as a strong alternative to boost the performance. In this paper, the performance of the DLR TAU code is analyzed and optimized. The improvement of the code efficiency is addressed through three key activities: Optimization, parallelization and hardware acceleration. At first, a profiling analysis of the most time-consuming processes of the Reynolds Averaged Navier Stokes flow solver on a three-dimensional unstructured mesh is performed. Then, a study of the code scalability with new partitioning algorithms are tested to show the most suitable partitioning algorithms for the selected applications. Finally, a feasibility study on the application of FPGAs and GPUs for the hardware acceleration of CFD simulations is presented

The GNAT method for nonlinear model reduction: effective implementation and application to computational fluid dynamics and turbulent flows

The Gauss--Newton with approximated tensors (GNAT) method is a nonlinear model reduction method that operates on fully discretized computational models. It achieves dimension reduction by a Petrov--Galerkin projection associated with residual minimization; it delivers computational efficency by a hyper-reduction procedure based on the `gappy POD' technique. Originally presented in Ref. [1], where it was applied to implicit nonlinear structural-dynamics models, this method is further developed here and applied to the solution of a benchmark turbulent viscous flow problem. To begin, this paper develops global state-space error bounds that justify the method's design and highlight its advantages in terms of minimizing components of these error bounds. Next, the paper introduces a `sample mesh' concept that enables a distributed, computationally efficient implementation of the GNAT method in finite-volume-based computational-fluid-dynamics (CFD) codes. The suitability of GNAT for parameterized problems is highlighted with the solution of an academic problem featuring moving discontinuities. Finally, the capability of this method to reduce by orders of magnitude the core-hours required for large-scale CFD computations, while preserving accuracy, is demonstrated with the simulation of turbulent flow over the Ahmed body. For an instance of this benchmark problem with over 17 million degrees of freedom, GNAT outperforms several other nonlinear model-reduction methods, reduces the required computational resources by more than two orders of magnitude, and delivers a solution that differs by less than 1% from its high-dimensional counterpart

The prospect of using LES and DES in engineering design, and the research required to get there

In this paper we try to look into the future to divine how large eddy and detached eddy simulations (LES and DES, respectively) will be used in the engineering design process about 20-30 years from now. Some key challenges specific to the engineering design process are identified, and some of the critical outstanding problems and promising research directions are discussed.Comment: accepted for publication in the Royal Society Philosophical Transactions

Physics based GMRES preconditioner for compressible and incompressible Navier-Stokes equations

This paper presents the implementation of a local physics preconditioning mass matrix [8] for an unified approach of 3D compressible and incompressible Navier-Stokes equations using an SUPG finite element formulation and GMRES implicit solver. During the last years a lot of effort has been dedicated to finding a unified approach for compressible and incompressible flow in order to treat fluid dynamic problems with a very wide range of Mach and Reynolds numbers [10,26,37]. On the other hand, SUPG finite element formulation and GMRES implicit solver is one of the most robust combinations to solve state of the art CFD problems [1,6,9,22,29,30,31]. The selection of a good preconditioner and its performance on parallel architecture is another open problem in CFD community. The local feature of the preconditioner presented here means that no communication among processors is needed when working on parallel architectures. Due to these facts we consider that this research can make some contributions towards the development of a unified fluid dynamic model with high rates of convergence for any combination of Mach and Reynolds numbers, being very suitable for massively parallel computations. Finally, it is important to remark that while this kind of preconditioning produces stabilized results in nearly incompressible regimes the standard version exhibits some numerical drawbacks that lead to solutions without physical meaning

Numerical simulation of a supersonic ejector for vacuum generation with explicit and implicit solver in openfoam

Supersonic ejectors are used extensively in all kind of applications: compression of refrigerants in cooling systems, pumping of volatile fluids or in vacuum generation. In vacuum generation, also known as zero-secondary flow, the ejector has a transient behaviour. In this paper, a numerical and experimental research of a supersonic compressible air nozzle is performed in order to investigate and to simulate its behaviour. The CFD toolbox OpenFOAM 6 was used, with two density-based solvers: explicit solver rhoCentralFoam, which implements Kurganov Central-upwind schemes, and implicit solver HiSA, which implements the AUSM+up upwind scheme. The behaviour of the transient evacuation ranges between adiabatic polytropic exponent at the beginning of the process and isothermal at the end. A model for the computation of the transient polytropic exponent is proposed. During the evacuation, two regimes are encountered in the second nozzle. In the supercritic regime, the secondary is choked and sonic flow is reached. In the subcritic regime, the secondary flow is subsonic. The final agreement is good with the two different solvers, although simulation tends to slightly overestimate flow rate for large values region.Peer ReviewedPostprint (published version

CFD investigation of a complete floating offshore wind turbine

This chapter presents numerical computations for floating offshore wind turbines for a machine of 10-MW rated power. The rotors were computed using the Helicopter Multi-Block flow solver of the University of Glasgow that solves the Navier-Stokes equations in integral form using the arbitrary Lagrangian-Eulerian formulation for time-dependent domains with moving boundaries. Hydrodynamic loads on the support platform were computed using the Smoothed Particle Hydrodynamics method. This method is mesh-free, and represents the fluid by a set of discrete particles. The motion of the floating offshore wind turbine is computed using a Multi-Body Dynamic Model of rigid bodies and frictionless joints. Mooring cables are modelled as a set of springs and dampers. All solvers were validated separately before coupling, and the loosely coupled algorithm used is described in detail alongside the obtained results

Demonstration of a coupled floating offshore wind turbine analysis with high-fidelity methods

This paper presents results of numerical computations for floating off-shore wind turbines using, as an example, a machine of 10-MW rated power. The aerodynamic loads on the rotor are computed using the Helicopter Multi-Block flow solver developed at the University of Liverpool. The method solves the Navier–Stokes equations in integral form using the arbitrary Lagrangian–Eulerian formulation for time-dependent domains with moving boundaries. Hydrodynamic loads on the support platform are computed using the Smoothed Particle Hydrodynamics method, which is mesh-free and represents the water and floating structures by a set of discrete elements, referred to as particles. The motion of the floating offshore wind turbine is computed using a Multi-Body Dynamic Model of rigid bodies and frictionless joints. Mooring cables are modelled as a set of springs and dampers. All solvers were validated separately before coupling, and the results are presented in this paper. The importance of coupling is assessed and the loosely coupled algorithm used is described in detail alongside the obtained results