The following energy industry HPC codes are being used and optimised by ENERXICO:

WRF (Weather Research and Forecasting Model) 

WRF is a mesoscale numerical weather prediction system designed for both atmospheric research and operational forecasting applications. It features two dynamical cores, a data assimilation system, and a software architecture supporting parallel computation and system extensibility. The model serves a wide range of meteorological applications across scales from tens of meters to thousands of kilometres.

Link: https://github.com/wrf-model/WRF/releases

BSIT (Barcelona Subsurface Imaging Tools) 

BSIT is a production geophysical imaging application with the capacity to run on extremely large systems implementing Full Waveform Inversion. FWI is a cutting-edge technique that aims to acquire the physical properties of the subsoil from a set of seismic measurements. Starting from a guess (initial model) of the variables being inverted (e.g., sound transmission velocity), the stimulus introduced and the recorded signals, Full Waveform Inversion performs several phases of iterative computations to reach the real value of the set of variables being inverted with an acceptable error threshold.

Link: https://github.com/albertfc/FWI_enerxico


Alya is a code developed by the Barcelona Supercomputing Center (BSC). It is a parallel multi-physics CFD code of the PRACE Benchmark Suite for HPC. The numerical discretization is based on a second-order spatial low-dissipation finite element scheme with an explicit temporal third-order Runge-Kutta method for momentum and scalar transport. Alya has been designed to run on leading HPC systems and includes both the MPI and OpenMP models to take advantage of the distributed and the shared memory paradigms. Accelerators like GPUs are also exploited to further enhance the performance of the code. 

The combustion problem is solved using a flamelet approach for which transport equations for user-defined controlling variables (i.e. mixture fraction and progress variable) and their corresponding variances are solved using a presumed-shape probability density function that accounts for turbulence-chemistry interactions in the subgrid scale. This model is coupled to a Lagrangian solver to account for the liquid phase of the spray with a two-way coupling approach. A Newmark/Newton–Raphson scheme is used to solve the kinematic equations which is coupled to heat and mass transfer models describing droplet heating and evaporation. Various combinations of heat and mass transfer models are available, with increasing complexity considering the interaction of the two phenomena and non-equilibrium effects. 

Link: https://gitlab.bsc.es/alya/alya/-/wikis/home

ExaHyPE (An Exascale Hyperbolic PDE Engine) 

ExaHyPE is an engine for solving systems of first-order hyperbolic partial differential equations. Due to the robustness and shock capturing abilities of ExaHyPE’s numerical methods, both linear and non-linear hyperbolic PDEs can be simulated with very high accuracy. 

Link: https://gitlab.lrz.de/exahype/ExaHyPE-Engine

DualSPHysics code 

DualSPHysics is used to develop an extension that will be the first SPH based code for the numerical simulation of oil reservoirs and which has also important benefits versus commercial codes based on other numerical techniques. 

This BH code (for Black Hole) is a large-scale massively parallel reservoir simulator capable of performing simulations with billions of “particles” or fluid elements that represents the system under study. It contains improved multi-physics modules that automatically combine the effects of interrelated physical and chemical phenomena to accurately simulate in-situ recovery processes. This leads to the development of a graphical user interface multiple-platform application for code execution and visualization, and for carrying out simulations with data provided by industrial partners and performing comparisons with available commercial packages.

Furthermore, a large effort is being made to simplify the process of setting up the input for reservoir simulations from exploration data by means of a workflow fully integrated in our industrial partners’ software environment.

Link: https://github.com/DualSPHysics/DualSPHysics/tree/develop


SEM46 is a seismic full waveform modelling and inversion code developed at Universite Grenoble Alpes. It aims to provide high resolution 3D models of the subsurface physical properties, from active seismic recordings. Such models can then be interpreted as geological formations, rock types, reservoirs and/or fluid content, in order to evaluate resources quantities and monitor production. 

Contact: romain.brossier@univ-grenoble-alpes.fr


SeiSol is a software package for simulating wave propagation and dynamic rupture based on the arbitrary high-order accurate derivative discontinuous Galerkin method (ADER-DG). Characteristics of the SeisSol simulation software are:
    • use of tetrahedral meshes to approximate complex 3D model geometries (faults & topography) and rapid model generation
    • use of elastic, viscoelastic and viscoplastic material to approximate realistic geological subsurface properties.

Link: https://github.com/SeisSol/SeisSol