preloading

    Flexi at a glance

    Multiscale Numerics

    The Flexi framework is built to solve the unsteady compressible Navier-Stokes equations (and more!) in a highly efficient way. Target applications are DNS and LES of internal and external compressible flows.

    Open and user-friendly

    Flexi and its tools are open-source under GPL V3.0, so everybody can enjoy it free of charge. Even better, the clear algorithmic design, the active community and the extensive documentation make it very easy to adapt or extend Flexi to your needs.

    Not just a solver

    The Flexi framework offers a complete CFD solution, and consists of the high-order mesh generator and preprocessor HOPR, the solver Flexi and a converter to the Paraview format for visualization and postprocessing.

    Fully unstructured

    Flexi supports fully unstructured curved hexahedral meshes with non-conforming element interfaces, allowing for complex geometries and meshing flexibility. Its preprocessor HOPR can generate curved meshes from a variety of linear meshes from open-source or commercial grid generators.

    High-order accurate

    Flexi is based on the highly efficient Discontinuous Galerkin Spectral Element Method, allowing you to exploit spectral convergence with very high orders of accuracy. Simulations utilizing O(20) and more are possible.

    Massively parallel

    The Flexi framework is designed to fully exploit modern parallel architectures up to maximum scale, with actual production runs over 100.000 cores. The solver itself and its I/O scale efficiently to multiple hundred thousand cores, the tool chain allows for parallel post-processing and visualization of huge amounts of data.

hp-Adaptive discretization

A key challenge when simulating complex flow phenomena is a temporal and spacial variation in the require resolution. In turbulent flow problems, the appearance of boundary layers, shear layers or complex wake flows and the associated decay of flow structures along a wide range of scales impose a very high resolution requirement. Further, in the transsonic and supersonic regime, as well as in the presence of phase interfaces, discontinuous solution features like shocks or material interfaces need to be captured and localized.
To tackle these challenges, we developed a hybrid hp-adaptive Discontinuous Galerkin/Finite Volume operator. It exploits the high-order convergence and scale resolving capabilities of the efficient Discontinuous Galerkin Spectral Element Method (DGSEM) with the robust Finite Volume (FV) method, well suited for shock and interface capturing. Based on a error and smoothness indicator, we employ a p-adaptive DGSEM operator in smooth areas, while discontinuous solution features are treated with a second order TVD FV operator on a h-refined sub-cell grid.
While this hybrid, adaptive discretization strategy allows for an efficient localization of the degrees of freedom, it entails significant load variations among processor the units in parallel computations. Using a dynamic load balancing scheme that repartitions the domain along a space filling curve, we can maintain favorable scaling properties and thus profit from the vast computational resources of todays massively parallel high performance architectures.

Literature:

  • An Efficient hp-Adaptive Strategy for a Level-Set Ghost-Fluid Method, Journal of Scientific Computing Journal of Scientific Computing, 2023
  • A p-Adaptive Discontinuous Galerkin Method with hp-Shock Capturing, Journal of Scientific Computing Journal of Scientific Computing, 2022

Multiphase flow simulation

Multi-phase flows are omnipresent in nature and relevant for many technical applications like the fuel injection systems of rocket or airplane engines. With the multiphase extension of FLEXI, we are at the forefront of studying complex multiphase problems like phase transition under extreme conditions, high speed multi-component jet-flows and droplet grouping. Key challenges in the field of multi-phase flow dynamics are the thermodynamic and numerical modelling of the interface and the work with complex real material equations of state (EOS).
The multi-phase FLEXI provides implementations of both a diffuse interface approach to simulate complex mixing and jet flows and a sharp-interface method to study the dynamics of individual droplets. In the sharp-interface context, interface modeling is achieved with the solution of an interfacial Riemann problem, that allows the incorporation of local thermodynamic phase transition models. A key building block of the multiphase code framework is an EOS-module, that supports a wide range of real material EOS from simple cubic material laws to multi-parameter eos. Affordable EOS-evaluation is facilitated with an efficient tabulation approach.

Literature:

  • On the use of tabulated equations of state for multi-phase simulations in the homogeneous equilibrium limit, Shock Waves, 2019
  • Recent Advances and Complex Applications of the Compressible Ghost-Fluid Method, Springer International Publishing, 2021
  • Riemann solvers for phase transition in a compressible sharp-interface method, Applied Mathematics and Computation, 2023

AI-Integrated CFD

Aritificial Intelligence (AI) and especially Machine learning (ML) have emerged as powerful tools in Computational Fluid Dynamics (CFD) that promise to accelerate the pace of innovation and discovery in fluid dynamics. ML models can be utilized in CFD as surrogates for physical models, for uncertainty optimization or to detect and extract flow patterns within highly complex turbulent flow fields. Hence, the application of ML models allows to reduce the computational cost while maintaining or even improving the accuracy of complex high-fidelity simulations.
As computational requirements and capabilities continue to advance in both fields, integrating AI and CFD efficiently is challenging due to their differences in the applied software and hardware. To address this, FLEXI incorporates ML models efficiently using for instance our TensorFlow-Fortran-Bridge that allows to deploy trained AI models in large-scale simulations with FLEXI. Moreover, the Relexi Framework can be used to train AI models with Reinforcement Learning using CFD simulations with FLEXI as the training environment to derive novel surrogate models or for applications in active low control. These approaches have been applied successfully for a wide variety of different applications ranging from shock localization and capturing methods to advanced to advanced data-driven turbulence models. With this, FLEXI yields the foundation for a novel class of intelligent, optimizable flow solvers.

Literature:

  • Deep neural networks for data-driven LES closure models, Journal of Computational Physics, 2019
  • A neural network based shock detection and localization approach for discontinuous Galerkin methods, Journal of Computational Physics, 2019
  • A data-driven high order sub-cell artificial viscosity for the discontinuous Galerkin spectral element method, Journal of Computational Physics, 2021
  • Relexi — A scalable open source reinforcement learning framework for high-performance computing, Software Impacts, 2922
  • Deep reinforcement learning for turbulence modeling in large eddy simulations, International Journal of Heat and Fluid Flow, 2023
  • Toward discretization-consistent closure schemes for large eddy simulation using reinforcement learning, Physics of Fluids, 2023
  • Towards Exascale CFD Simulations Using the Discontinuous Galerkin Solver FLEXI, Preprint, 2023

Complex Simulations of Turbulent Flows

Most flows in engineering and nature are turbulent. A good unsterstanding of turbulence is thus essential for improving the performance and efficiency of engineering systems as well as predicting the evolution of earth's weather and climate. However, turbulent flow is characterized by a strong multiscale behavior where the active small features span multiple orders of magnitude in size. Resolving such flow phenomena in numerical simulations is challenging, due to the massive spatial resolution that is required to capture the nonlinear, chaotic dynamics of the flow.
The high-order accuracy and excellent scaling properties of FLEXI renders it the perfect candidate for computationally intensive, high-fidelity simulations of complex turbulent flow. For this, FLEXI provides a wide variety of different features to allow for the accurate simulation of such flows. For instance, FLEXI allows for moving domains using a sliding mesh approach or a Chimera overset mesh approach in combination with an arbitrary Lagrangian-Eulerian (ALE) formulation, which makes it suitable to capture fluid-structure interactions. Moreover, different shock indicators and shock-capturing methods based on a finite-volume subcell method are implemented in FLEXI to yield accurate and stable representations of shocks in compressible flow. Lastly, sponge zones, various turbulence models and a wide variety of different boundary conditions including turbulent inflow methods complete the rich set of features available in FLEXI.

Literature:

  • FLEXI: A high order discontinuous Galerkin framework for hyperbolic–parabolic conservation laws, Computers & Mathematics with Applications, 2021
  • An efficient sliding mesh interface method for high-order discontinuous Galerkin schemes, Computers & Fluids, 2021
  • Boundary-Layer Dynamics in Wall-Resolved LES Across Multiple Turbine Stages, AIAA Journal, 2021
  • Increasing the Flexibility of the High Order Discontinuous Galerkin Framework FLEXI Towards Large Scale Industrial Applications, High Performance Computing in Science and Engineering '20, 2022
  • Development of turbulent inflow methods for the high order HPC framework FLEXI, High Performance Computing in Science and Engineering '21, 2023
  • Towards Exascale CFD Simulations Using the Discontinuous Galerkin Solver FLEXI, Preprint, 2023