Tommaso Benacchio

ESCAPE-2 project

Phone:+39 02 2399 4508
Fax: +39 02 2399
Office: 14 - II floor

  • Available MOX Reports
  • Theses
  • Thesis Proposals
  • MOX Projects

Simulation of wave propagation in unbounded domains

The project aims to devise and implement numerical discretizations of hyperbolic systems of geophysical fluid dynamics (shallow water equations, compressible Euler equations) on multidimensional semi-infinite strips, building on existing models for the one-dimensional case developed by T. Benacchio and L. Bonaventura. Continuous and discontinuous finite element discretizations will be used in the finite directions and Laguerre-based spectral methods will be used in the unbounded direction. Subject to good performance on standalone tests, the model will be coupled to a finite-element based discretization on a finite domain. The coupled setup will be used to perform wave propagation experiments, including tests in thermally stratified environments.
As one of several applications, the resulting code will supply efficient absorbing boundary conditions at model tops of state-of-the-art numerical weather prediction systems. The project is in collaboration with G. Tumolo (European Centre for Medium-Range Weather Forecasts, UK).


Benacchio, T., and L. Bonaventura, 2013: Absorbing boundary conditions: a spectral collocations approach. International Journal of Numerical Methods in Fluids, 72, 913-936, doi: 10.1002/fld.3768.

Benacchio, T. and L. Bonaventura, 2019: An extension of DG methods for hyperbolic problems to one-dimensional semi-infinite domains. Applied Mathematics and Computation, 350, 266-282, doi: 10.1016/j.amc.2018.12.057.

Tumolo, G. and L. Bonaventura, 2015: A semi-implicit, semi-Lagrangian discontinuous Galerkin framework for adaptive numerical weather prediction. Quarterly Journal of the Royal Meteorological Society, 141, 2582-2601, doi: 10.1002/qj.2544.

Efficient and scalable solvers for numerical weather prediction

Timely and trustworthy weather forecasts rely on accurate and efficient underlying numerical methods, and numerical weather prediction simulations are the ideal testbed for high-performance computing hardware towards the exascale. The project aims at exploring and extending numerical algorithms used in state-of-the-art models of atmospheric dynamics by:

1. On existing code bases, performing efficiency and scalability studies with linear solvers applied to the solution of simplified models of atmospheric flow;
2. Based on the development at point 1., applying the techniques to the implementation of solvers within semi-implicit discretizations of more complex equation sets and analyze their performance on benchmarks at different spatial and temporal scales.

The programming activity will be carried out using recently developed mini-apps ('dwarfs') produced by the EU project ESCAPE, and will be strongly connected to the current activities of the ESCAPE-2 EU project at MOX. Knowledge of, or willingness to learn, Fortran and shell scripting would be required for the project.

The project includes a partially funded research stay at the European Centre for Medium-Range Weather Forecasts in Reading, UK. This will give the opportunity to interact with scientists in the Numerical Methods team and to present the results of the work in a Research Seminar at the Centre.


ESCAPE-2 project:

Melvin, T., T. Benacchio, B. Shipway, N. Wood, J. Thuburn, and C. J. Cotter, 2018: A mixed finite-element, finite-volume, semi-implicit discretisation for atmospheric dynamics: Cartesian geometry. Q. J. Roy. Meteor. Soc., in press, doi:

Mueller, A., et al., 2019: The ESCAPE project: Energy-efficient Scalable Algorithms for Weather Prediction at Exascale, Geosci. Model Dev. Discuss.,, in review, 2019.

Tumolo, G. and L. Bonaventura, 2015: A semi-implicit, semi-Lagrangian discontinuous Galerkin framework for adaptive numerical weather prediction. Quarterly Journal of the Royal Meteorological Society, 141, 2582-2601, doi: 10.1002/qj.2544.

Escape-2 Energy-efficient Scalable Algorithms For Weather And Climate Prediction At Exascale,European Commission
ESCAPE-2 will develop world-class, extreme-scale computing capabilities for European operational numerical weather and climate prediction systems. It continues the pioneering work of the ESCAPE project. The project aims to attack all three sources of enhanced computational performance at once, namely (i) developing and testing bespoke numerical methods that optimally trade off accuracy, resilience and performance, (ii) developing generic programming approaches that ensure code portability and performance portability, (iii) testing performance on HPC platforms offering different processor technologies. ESCAPE-2 will prepare weather and climate domain benchmarks that will allow a much more realistic assessment of application specific performance on large HPC systems than current generic benchmarks such as HPL and HPCG. These benchmarks are specifically geared towards the pre-exascale and exascale HPC infrastructures that the European Commission and Member States will invest in through the European Commission Joint Undertaking. The role of MOX in ESCAPE-2 is the development of efficient semi-implicit solvers for Discontinuous Galerkin discretizations of the Euler equations. ESCAPE-2 is funded by the European Commission under the Future and Emerging Technologies - High-Performance Computing call for research and innovation actions, grant agreement 800897.