Numerical models are very useful for environmental applications. Several difficulties must be handled simultaneously, in a multidisciplinary context. For example, in geophysics, media are highly heterogeneous and only few data are available. Stochastic models are thus necessary. Some computational domains are complex 3D geometries, requiring adapted space discretization. Equations modeling flow and transport are transient, requiring also adapted time discretization. Moreover, these equations can be coupled together or with other equations in a global nonlinear system. These large-scale models are very time and memory consuming. High performance computing is thus required to run these types of scientific simulations. Supercomputers and clusters are quite powerful, provided that the numerical models are written with a parallel paradigm.

The team SAGE undertakes research on environmental applications and high performance computing and deals with two subjects:

numerical algorithms, involving parallel and grid computing,

numerical models applied to hydrogeology and physics.

These two subjects are highly interconnected: the first topic aims at designing numerical algorithms, which lead to high efficiency on parallel and grid architectures; these algorithms are applied to geophysical models.

Moreover, the team SAGE, in collaboration with other partners, develops the software platform H2OLab for groundwater numerical simulations in heterogeneous subsurface.

Linear algebra is at the kernel of most scientific applications, in particular in physical or chemical engineering.
For example, steady-state flow simulations in porous media are discretized in space and lead to a large sparse linear system.
The target size is

In geophysics, a main concern is to solve inverse problems in order to fit the measured data with the model. Generally, this amounts to solve a linear or nonlinear least-squares problem. Complex models are in general coupled multi-physics models. For example, reactive transport couples advection-diffusion with chemistry. Here, the mathematical model is a set of nonlinear Partial Differential Algebraic Equations. At each timestep of an implicit scheme, a large nonlinear system of equations arise. The challenge is to solve efficiently and accurately these large nonlinear systems.

Approximation in Krylov subspace is in the core of the team activity since it provides efficient iterative solvers for linear systems and eigenvalue problems as well. The later are encountered in many fields and they include the singular value problem which is especially useful when solving ill posed inverse problems.

The team Sage is strongly involved in numerical models for hydrogeology and physics. There are many scientific challenges in the area of groundwater simulations. This interdisciplinary research is very fruitful with cross-fertilizing subjects. For example, high performance simulations were very helpful for finding out the asymptotic behaviour of the plume of solute transported by advection-dispersion. Numerical models are necessary to understand flow transfer in fractured media.

The team develops stochastic models for groundware simulations. Numerical models must then include Uncertainty Quantification methods, spatial and time discretization. Then, the discrete problems must be solved with efficient algorithms. The team develops parallel algorithms for complex numerical simulations and conducts performance analysis. Another challenge is to run multiparametric simulations. They can be multiple samples of a non intrusive Uncertainty Quantification method, or multiple samples of a stochastic method for inverse problems, or multiple samples for studying the sensitivity to a given model parameter. Thus these simulations are more or less independent and are well-suited to grid computing but each simulation requires powerful CPU and memory resources.

A strong commitment of the team is to develop the scientific software platform H2OLab for numerical simulations in heterogeneous hydrogeology.

The team has chosen a particular domain of application, which is geophysics. In this domain, many problems require solving large scale systems of equations, arising from the discretization of coupled models. Emphasis is put on hydrogeology, but the team also investigates geodesy, heat and mass transfer in soil, and granular materials. One of the objectives is to use high performance computing in order to tackle 3D large scale computational domains with complex physical models.

This is joint work with Geosciences Rennes at OSUR, Pprime at University of Poitiers and CDCSP at University of Lyon. It is also done in the context of the group Momas and Andra grants.

Many environmental studies rely on modelling geo-chemical and hydrodynamic processes. Some issues concern water resources, aquifer contamination, underground waste disposal, clean-up of former waste deposits, acid mine drainage remediation. Other issues, also related to energy, concern geothermy, unconventional gas, enhanced oil recovery, underground storage of CO2, underground storage of nuclear waste.

Simulation of contaminant transport in groundwater is a highly complex problem, governed by coupled linear or nonlinear PDAEs. Moreover, due to the lack of experimental data, stochastic models are used for dealing with heterogeneity. The main objective of the team is to design and to implement efficient and robust numerical models, including Uncertainty Quantification methods.

Recent research showed that rock solid masses are in general fractured and that fluids can percolate through networks of inter-connected fractures. Fractured media are by nature very heterogeneous and multi-scale, so that homogenisation approaches are not relevant. The team develops a numerical model for fluid flow and contaminant transport in three-dimensional porous fractured media.

An important output is the parallel scientific platform H2OLab, running on clusters, grids and machines available in supercomputing centers.

The team SAGE ended in December 2015.

N. Nassif, J. Erhel and B. Philippe published a book entitled "introduction to computational linear algebra" . E. Gallopoulos, B. Philippe and A. Sameh published a book entitled "Parallelism in Matrix Computations".

Functional Description

GENFIELD allows the generation of log-normal correlated fields. It is based on a spectral method and uses the FFTW library. Parallelism is implemented using MPI communications. GENFIELD is used in hydrogeology to model natural fields, like hydraulic conductivity or porosity fields.

Authors: Grégoire Lecourt, Jocelyne Erhel, Jean-Raynald De Dreuzy and Géraldine Pichot

Contact: Jocelyne Erhel

Grains with Liquid Meniscus under Heating

Functional Description

The GLiMuH code is devoted to the understanding of how heat diffuses in an assembly of solid grains separated by air and water. In the pendular regime, the quantity of water is very small, leading to liquid bridges between the grains. In the current approximation, the grains are spherical in shape, and the numerical simulation is done in a 3D axisymmetric coordinate system. The shape of the liquid/gas interface is computed by integrating a differential algebraic system of equations, with a given quantity of water per unit volume of the porous medium, and under the constraint of a given contact angle between the liquid/gas interface and the solid boundaries. The numerical results allow us to estimate the effective thermal conductivity of a real wet granular medium, which is required to establish more realistic models for the HeMaTiS code.

Authors: Édouard Canot, Salwa Mansour and Renaud Delannay

Contact: Édouard Canot

Keywords: Advection - Simulation - Scientific calculation - Dispersion - Geochemistry

APP: version 2.0, April 2014, registered

Programming language: C

Current work: refactoring

Functional Description

Reactive transport modeling has become an essential tool for understanding complex environmental problems. It is an important issue for MoMaS and C2S@EXA partners, in particular Andra. We have developed a method coupling transport and chemistry, based on a method of lines such that spatial discretization leads to a semi-discrete system of algebraic differential equations (DAE system). The main advantage is to use a complex DAE solver, which controls simultaneously the timestep and the convergence of Newton algorithm. The approach SIA uses a fixed-point method to solve the nonlinear system at each timestep, whereas the approach SNIA uses an explicit scheme.

Authors: Yvan Crenner, Caroline De Dieuleveult, Jocelyne Erhel, Souhila Sabit, Nadir Soualem

Partner: ANDRA

Contact: Jocelyne Erhel

Heat and Mass Transfer in Soils

Functional Description

HeMaTiS is a set of Finite Volume programs (variants concern different geometrical configurations: 1D, 1D-radial, 2D, 3D-axisymmetric) for computing the transient heat diffusion in soils when there is a phase change of water. Currently, the soil is modelled by a heterogeneous porous medium having constant thermo-physical properties, and the porous medium is saturated with water. The phase change is treated by means of the Apparent Heat Capacity method. In the near future, we plan to use an unsaturated model (but limited to small water content), and an effective thermal conductivity which depends on the local humidity (this latter law may reveal hysteresis behaviour). The software is written in Fortran 95 and is based on the Muesli library. A Computer Algebra System (Maple or Maxima) is used to compute the Jacobian matrix.

Authors: Édouard Canot, Mohamad Muhieddine, Salwa Mansour and Renaud Delannay

Contact: Édouard Canot

Functional Description PALMTREE is an easy-to-use library for the parallelization of Lagrangian methods for partial differential equations and general Monte Carlo methods. The code aims at satisfying three properties:

Reduction of computation time by using parallel architecture,

Simplicity as the user just has to add the algorithm governing the behaviour of the particles,

Portability since one has the possibility to use the package with any compiler and OS,

Action-Replay which provides the ability of the package to replay a selected batch of particles.

The software also now possesses a beta version which allows to run in parallel hybrid solvers, that is solvers which use both deterministic methods (FEM, FDM, etc...) and probabilistic methods. Moreover, the engineering being this software was published in MCQMC 2014 .

Authors: Lionel Lenôtre, Géraldine Pichot.

Contact: Lionel Lenôtre.

Skew Brownian Motion

Functional Description

SBM is a code allowing exact or approximated simulations of the Skew Brownian Motion. This code is used for the simulation, with a Monte-Carlo approach, of a 1D diffusion process with a discontinuous diffusion coefficient. Several benchmark tests are also implemented.

Authors: Antoine Lejay and Géraldine Pichot

Contact: Géraldine Pichot

Thermal Properties by Inverse Problem

Functional Description

TPIP is a program which aims at estimating the thermo-physical of a saturated porous medium after a strong heating which leads to the phase change of the water contained in the pores, knowing the experimental heating curves history at few selected points. The least-square criterion is used, in which sensitivity coefficients are the solution of a huge, complex PDE system in order to take into account the phase change of water. These equations for the sensitivity coefficients are therefore obtained via a Computer Algebra System (Maple or Maxima). In many aspects, the forward problem is similar to the HeMaTiS code, and like it, is based on Muesli. Two different minimization algorithms may be used, Damped Gauss-Newton or Levenberg-Marquardt. A special procedure has been applied in order to obtain a robust convergence, by changing some parameters of the forward problem during the iterations.

Authors: Édouard Canot, Mohamad Muhieddine, Salwa Mansour and Renaud Delannay

Contact: Édouard Canot

Functional Description

Zohour is a node-based adaptive 2D mesh algorithm, written in Fortran 2003. A basic rectangular, regular set of nodes is recursively refined. Then the cells come from the Voronoi tessellation. While the domain is currently limited to a rectangular shape, its strength is three-fold:first, computing the flux via a Finite Element or Finite Volume method is both simple and accurate because each cell-side of cells is the bisection of two nodes,

second, the transition between zones of different levels of refinement is more progressive than other methods, leading to a smaller number of nodes for the whole mesh,

third, during successive refinements when dealing with a transient problem, interpolation is needed only by the new nodes, limiting the numerical errors.

Author: Édouard Canot

Contact: Édouard Canot

Abstract: The book "Introduction to Computational Linear Algebra" presents classroom-tested material on computational linear algebra and its application to numerical solutions of partial and ordinary differential equations. The book is designed for senior undergraduate students in mathematics and engineering as well as first-year graduate students in engineering and computational science.

The text first introduces BLAS operations of types 1, 2, and 3 adapted to a scientific computer environment, specifically MATLAB. It next covers the basic mathematical tools needed in numerical linear algebra and discusses classical material on Gauss decompositions as well as LU and Cholesky’s factorizations of matrices. The text then shows how to solve linear least squares problems, provides a detailed numerical treatment of the algebraic eigenvalue problem, and discusses (indirect) iterative methods to solve a system of linear equations. The final chapter illustrates how to solve discretized sparse systems of linear equations. Each chapter ends with exercises and computer projects.

Grants and projects: EXA2CT , EoCoE , C2S@EXA

Publications: in preparation.

Abstract: Sparse linear systems arise in computational science and engineering. The goal is to reduce the memory requirements and the computational cost, by means of high performance computing algorithms. Krylov methods combined with Domain Decomposition are very efficient for both fast convergence and fast computations.

Abstract: Shielding studies in neutron transport, with Monte Carlo codes, yield challenging problems of small-probability estimation. The particularity of these studies is that the small probability to estimate is formulated in terms of the distribution of a Markov chain, instead of that of a random vector in more classical cases. Thus, it is not straightforward to adapt classical statistical methods, for estimating small probabilities involving random vectors, to these neutron-transport problems. A recent interacting-particle method for small-probability estimation, relying on the Hastings-Metropolis algorithm, is presented. It is shown how to adapt the Hastings-Metropolis algorithm when dealing with Markov chains. A convergence result is also shown. Then, the practical implementation of the resulting method for small-probability estimation is treated in details, for a Monte Carlo shielding study. Finally, it is shown, for this study, that the proposed interacting-particle method considerably outperforms a simple Monte Carlo method, when the probability to estimate is small.

Abstract: We present some investigations on the parallelization of a stochastic Lagrangian simulation. For the self sufficiency of this work, we start by recalling the stochastic methods used to solve Parabolic Partial Differential Equations with a few physical remarks. Then, we exhibit different object-oriented ideas for such methods. In order to clearly illustrate these ideas, we give an overview of the library PALMTREE that we developed. After these considerations, we discuss the importance of the management of random numbers and argue for the choice of a particular strategy. To support our point, we show some numerical experiments of this approach, and display a speedup curve of PALMTREE. Then, we discuss the problem in managing the parallelization scheme. Finally, we analyze the parallelization of hybrid simulation for a system of Partial Differential Equations. We use some works done in hydrogeology to demonstrate the power of such a concept to avoid numerical diffusion in the solution of Fokker-Planck Equations and investigate the problem of parallelizing scheme under the constraint entailed by domain decomposition. We conclude with a presentation of the latest design that was created for PALMTREE and give a sketch of the possible work to get a powerful parallelized scheme.

Grants: HYDRINV , H2MN04

Publications: .

Abstract: Uncertainty quantification often requires the generation of large realizations of stationary Gaussian random field over a regular grid.

We compare the classical methods used to simulate the field defined by its covariance function, namely the Discrete Spectral method, the Circulant Embedding approach, and the Discrete Karhunen-Loève approximation. We design and implement a parallel algorithm related to the Discrete Spectral method.

Grants: ECOS Sud Chili (ARPHYMAT project)

Abstract: This study is devoted to the heat transfer between two spherical grains separated by a small gap; dry air is located around the grains and a liquid water meniscus is supposed to be present between them. This problem can be seen as a micro-scale cell of an assembly of solid grains, for which we are looking for the effective thermal conductivity. For a fixed contact angle and according to the volume of the liquid meniscus, two different shapes are possible for the meniscus, giving a “contacting” state (when the liquid makes a true bridge between the two spheres) and a “non-contacting” one (when the liquid is split in two different drops, separated by a thin air layer); the transition between these two states occurs at different times when increasing or decreasing the liquid volume, thus leading to a hysteresis behavior when computing the thermal flux across the domain. We consider also another process where humidity varies, for example during an evaporation or condensation process; in this situation, the shape of the menisci changes a lot, because some liquid bridges may break, and this can strongly affect the effective thermal conductivity. Then, the reorganization of the liquid menisci is predicted, especially their surface area variation; it is an important parameter for a global model of the evaporation phenomenon in wet porous media.

Software: TPIP ()

Publications: ,

Abstract: This study concerns the inverse problem which consists of the estimation of thermophysical properties of the soil knowing the temperature at few selected points of the domain. In order to solve this inverse problem, we used the least square criterion where we try to minimize the error function between real measures and simulated ones. The coupled system composed of the energy equation together with the three sensitivity boundary initial problems resulting from differentiating the basic energy equation with respect to the soil properties must be solved. To overcome the stiffness of our problem (due to the use of Apparent Heat Capacity method), the high nonlinearity of the coupled system and the problem of large residuals we used the Damped Gauss Newton and Levenberg-Marquardt methods. To take into account uncertainties of the position of the sensors, some constraints have been added to the least square problem. Results are good when the number of sensors is sufficiently large.

Grants: ECOS Sud Chili (ARPHYMAT project)

Software: HeMaTis ()

Publications: ,

Publications: submitted.

Abstract: We present several benchmark tests for Monte Carlo methods for simulating diffusion in one-dimensional discontinuous media, such as the ones arising the geophysics and many other domains. These benchmarks tests are developed according to their physical, statistical, analytic and numerical relevance. We then perform a systematic study on four numerical methods.

Abstract: The study of skew diffusion is of primary concern for their implication in the modeling and simulation of diffusion phenomenons in media with interfaces. First, we provide results on one-dimensional processes with discontinuous coefficients and their connections with the Feller theory of generators as well as the one of stochastic differential equations involving local time. Second, in view of developing new simulation techniques, we give a method to compute the density and the resolvent kernel of skew diffusions. Explicit closed-form are given for some particular cases.

Grants and projects: H2MNO4 .

Publications: .

Abstract: Two new numerical schemes are created for Skew Diffusions processes. Both algorithms rely on a more generic numerical scheme that can be used for any kind of Feller processes. The proof of convergence for this generic numerical scheme is performed.

Grants and projects: H2MNO4 .

Publications: .

Abstract: Some significant results on the distribution of the marginal processes of multidimensional Skew Diffusions are found together with new formula. In addition, totally analytical proofs of some results and algorithms given by A. Lejay are given.

Grants and projects: HYDRINV , H2MN04

Publications: .

Abstract: The stochastic collocation method based on the anisotropic sparse grid has become a significant tool to solve partial differential equations with stochastic inputs. The aim is to seek a vector of weights and a convenient level of interpolation for the method. The classical approach uses an a posteriori approach on the solution, which causes an additional prohibitive cost.

In this work, we discuss an adaptive approach of this method to calculate the statistics of the solution. It is based on an adaptive approximation of the *inverse* diffusion parameter. We construct an efficient error indicator which is an upper bound of the error on the solution. In the case of unbounded variables, we use an appropriate error estimation to compute suitable weights for the method. Numerical examples are presented to confirm the efficiency of the approach, and to show that the cost is considerably reduced without loss of accuracy.

Publications: .

Abstract: Reactive transport models are very useful for groundwater studies such as water quality, safety analysis of waste disposal, remediation, and so on. The MoMaS group defined a benchmark with several test cases. We present results obtained with a global method and show through these results the efficiency of our numerical model.

Abstract: Reactive transport models are very useful to study the fate of contaminants in grounwater. These models couple transport equations with geochemistry equations. In this talk, we focus on precipitation and dissolution chemical reactions, because they induce numerical difficulties.

We consider a set of solute species and minerals, with precipitation occuring when a saturation threshold is reached. A challenge is to detect which minerals are dissolved and which minerals are precipitated. This depends on the total quantities of chemical species. We propose an analytical approach to build a phase diagram, which provides the interfaces between the different possible cases. We illustrate our method with three examples arising from brine media and acid mine drainage.

Publications: .

Abstract: In a lot of geological environments, permeability is dominated by the existence of fractures and by their degree of interconnections. Flow properties depend mainly on the statistical properties of the fracture population (length, apertures, orientation), on the network topology, as well as on some detailed properties within fracture planes. Based on an extensive analysis of 2D and 3D DFNs as well as on reference connectivity structures, we investigate the relation between the local fracture structures and the effective permeability. Defined as the relative weight between the two extreme harmonic and arithmetic means, the power-law averaging exponent gives a compact way to compare fracture network hydraulics. It may further lead to some comprehensive upscaling rules.

Contract with ANDRA (National Agency for Nuclear Waste)

Duration: three years from November 2015.

Title: reactive transport in fractured porous media

Coordination: Jocelyne Erhel.

Partners: Geosciences Rennes.

Web page: http://

Abstract: Even in small numbers, fractures must be carefully considered for the geological disposal of radioactive wastes. They critically enhance diffusivity, speed up solute transport, extend mixing fronts and, in turn, modify the physicochemical conditions of reactivity around possible storage sites. Numerous studies in various fields have shown that fractures cannot be simply integrated within an equivalent po- rous medium with a simple enhancement of its petro-physical properties (porosity and permeability). We propose a combined numerical and experimental approach to determine the influence on reactivity of typical fracture patterns found in some radioactive waste applications.

Contract with ANR, program Modèles Numériques

Duration: four years from November 2012.

Title: Original Optimized Object Oriented Numerical Model for Heterogeneous Hydrogeology.

Coordination: Jocelyne Erhel and Géraldine Pichot, with Fabienne Cuyollaa.

Partners: Geosciences Rennes, University of Poitiers, University of Lyon 1, Andra, Itasca.

International collaborations: University of San Diego (USA), UPC, Barcelona (Spain)

Web page: http://

Abstract: The project H2MNO4 develops numerical models for reactive transport in heterogeneous media. It defines six mathematical and computational challenges and three applications for environmental problems with societal impact. We organized a project meeting in February.

Title: C2S@EXA - Computer and Computational Scienecs at Exascale

Duration: from January 2012.

Coordination: S. Lanteri, Nachos team.

Partners: Inria teams working on HPC; external partners: ANDRA and CEA.

Webpage: http://

Abstract: The C2S@Exa Inria Project Lab is concerned with the development of numerical modeling methodologies that fully exploit the processing capabilities of modern massively parallel architectures in the context of a number of selected applications related to important scientific and technological challenges for the quality and the security of life in our society. The team participated in several workshops.

Title: Numerical models for hydrogeology

Duration: 2015

Coordination: J. Erhel and G. Pichot.

Webpage: http://

Abstract: To run large scale simulations, we defined a project, based on the platform H2OLab and on a new GMRES solver. We obtained and used computing time on machines located at GENCI supercomputing centers.

Title: Workshops on multiphase flow and reactive transport

Duration: 2015

Coordination: J. Erhel

Partner: IFPEN

Webpage: http://

Abstract: The working group MOMAS includes many partners from CNRS, Inria, universities, CEA, ANDRA, EDF and BRGM. It covers many subjects related to mathematical modeling and numerical simulations for nuclear waste disposal problems. The team participated in a workshop on multiphase flow and organized an international workshop on reactive transport, in Paris, with IFPEN.

Title: EXascale Algorithms and Advanced Computational Techniques

Programm: FP7

Duration: September 2013 - August 2016

Coordinator: S. Ashby, IMEC, Belgium

Partners:

Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V (Germany)

Interuniversitair Micro-Electronica Centrum Vzw (Belgium)

Intel Corporations (France)

Numerical Algorithms Group Ltd (United Kingdom)

Systems Solutions for Research (Germany)

Universiteit Antwerpen (Belgium)

Universita della Svizzera italiana (Switzerland)

Universite de Versailles Saint-Quentin-En-Yvelines. (France)

Vysoka Skola Banska - Technicka Univerzita Ostrava (Czech Republic)

Inria contact: Luc Giraud

Abstract: Numerical simulation is a crucial part of science and industry in Europe. The advancement of simulation as a discipline relies on increasingly compute intensive models that require more computational resources to run. This is the driver for the evolution to exascale. Due to limits in the increase in single processor performance, exascale machines will rely on massive parallelism on and off chip, with a complex hierarchy of resources. The large number of components and the machine complexity introduce severe problems for reliability and programmability. The former of these will require novel fault-aware algorithms and support software. In addition, the scale of the numerical models exacerbates the difficulties by making the use of more complex simulation algorithms necessary, for numerical stability reasons. A key example of this is increased reliance on solvers. Such solvers require global communication, which impacts scalability, and are often used with preconditioners, increasing complexity again. Unless there is a major rethink of the design of solver algorithms, their components and software structure, a large class of important numerical simulations will not scale beyond petascale. This in turn will hold back the development of European science and industry which will fail to reap the benefits from exascale. The EXA2CT project brings together experts at the cutting edge of the development of solvers, related algorithmic techniques, and HPC software architects for programming models and communication. It will take a revolutionary approach to exascale solvers and programming models, rather than the incremental approach of other projects. We will produce modular open source proto-applications that demonstrate the algorithms and programming techniques developed in the project, to help boot-strap the creation of genuine exascale codes.

Program: EINFRA-5-2015

Project acronym: EoCoE

Project title: Energy oriented Center of Excellence for computer applications

Duration: 36 months

Coordinator: CEA

Other partners: organisme, labo (pays) : 12 other partners

Abstract: the EoCoE objectives aims at firstly, to design, test and spread new methodological and organisational paradigms (Objectives 1, 3, and 4) driven by the users communities and, secondly, to contribute to mathematical and computer sciences challenges on the whole HPC tool chain (Objective 2).

Program: Laboratoire International de Recherche en Informatique et Mathématiques Appliquées

Title: Problèmes Inverses et Contrôle

Inria principal investigator: Houssem Haddar, Defi team

International Partner (Institution - Laboratory - Researcher): ENIT, University of Tunis, Tunisia - LAMSIN - Amel ben Abda

Duration: 2011-2015

See also: http://

Abstract: The team deals with nonlinear and inverse problems.

Program: CONICYT

Title: Processus de formation et transformation de structures de combustion archéologique : un regard interdisciplinaire

Inria principal investigator: Édouard CANOT

International Partner (Institution - Laboratory - Researcher): Universidad de Tarapaca (Chili)

Duration: Jan 2014 - Dec 2016

Abstract: Multidisciplinary study of prehistoric fire traces in South America, by means of different approaches: taphonomy of the soil, physical processes involved during the heat transfer, modeling and numerical simulations.

Program: Euromediterranean 3+3

Title: Direct and inverse problems in subsurface flow and transport

Coordination: H. ben Ameur, ENIT, Tunisia and J. Jaffré, Inria, Paris

Inria-Rennes principal investigator: Jocelyne Erhel

Duration: Jan 2012 - Dec 2015

International Partners (Institution - Laboratory - Researcher):

Université Ibn Tofail - Faculté des Sciences de Kénitra (Morocco) - Laboratoire Interdisciplinaire en Ressources Naturelles et en Environnement - Zoubida Mghazli

Ecole Nationale d'Ingénieurs de Tunis (Tunisia) - Laboratoire de Modélisation en Hydraulique et Environnement - Rachida Bouhlila

Universidad de Sevilla (Spain) - Department Ecuaciones Diferenciales y Anålisis Numérico - Tomas Chacon Rebollo

Universitat Politècnica de Catalunya (Spain) - Department of Geotechnical Engineering and Geo-Sciences - Xavier Sànchez Vila

University Centre of KHEMIS MILIANA (Algeria) - Laboratoire de l'Energie et des Systèmes Intelligents - Mohammed Hachama

Ecole Mohammadia d'Ingénieurs (Morocco) - LERMA - Rajae Aboulaich

Ecole Nationale d'Ingénieurs de Tunis (Tunisia) - Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur - Hend Ben Ameur

The management of water resources is a problem of great importance in all countries, and is particularly acute around the Mediterranean sea. The goal is to find a reasonable balance between these resources and demand while preserving the quality of water. Towards this goal it is essential to understand and simulate flow and transport in the subsurface. The science corresponding to this topic is hydrogeology. Since models become more and more complicated and quantitative answers must be given, numerical modeling become more and more sophisticated and mathematicians must also be involved. This project brings together hydrogeologists and mathematicians from France, Spain, Algeria, Morocco and Tunisia in order to develop, analyze, and validate numerical methods for several problems arising from modeling flow and transport in the subsurface. The emphasis is put on direct nonlinear problems (air-water flow, density driven flow related to salinization, transport with chemistry) and on inverse problems.

Nabil Nassif, American University of Beirut, Lebanon, April, 1 week

Lamia Guellouz, ENIT, Tunisia, May, 1 week

Marwen ben Refifa, University of Tunis, 5 months, April-July 2015

Salwa Mansour, Lebanese University, 8 months, Feb-Sep 2015

Lionel Lenôtre visited Pr. Dr. Sylvie Roelly and her students at Potsdam University, Potsdam, Germany, one week, December 2015.

J. Erhel organized with T. Faney and A. Michel, from IFPEN, an international MoMaS workshop on reactive transport (Paris, France, November 2015).

J. Erhel is a member of the international advisory committee of the parallel CFD conferences (Trondheim, Norway, May 2014).

J. Erhel was a member of the scientific committee of the international conference MAMERN 2015 (Pau, France, June 2015).

J. Erhel was a member of the program committee of the workshop Visualization in Environmental Sciences 2015 (co-event of EuroVis 2015, Cagliari, Italy, May 2015)

J. Erhel is member of the editorial board of ETNA.

J. Erhel is member of the editorial board of ESAIM:Proceedings and Surveys.

J. Erhel is the scientific coordinator of the website Interstices (since June 2012).

J. Erhel was reviewer for the journals ADWR, CONHYD, GMDD, LAA.

L. Lenôtre was reviewer for the journal ETNA.

M. Oumouni was reviewer for the journal JOMP.

G. Pichot was reviewer for the journal Computers& Geosciences.

J. Erhel was reviewer for a proposal submitted to ICT15, Vienna Science and Technology fund.

É. Canot is member of the CLHSCT (Commité Local Hygiène Sécurité Conditions de Travail), of Inria-Rennes, from September 2010.

J. Erhel is the correspondant of Maison de la Simulation for Inria Rennes.

J. Erhel is the AMIES correspondant for Inria Rennes, from September 2015.

J. Erhel was a member of the Inria local committee Post-Doc 2015.

G. Pichot is responsible for the domain "environment" at IRISA, from June 2013 until August 2015.

G. Pichot is member of the Conseil de département MAM of Polytech Lyon.

G. Pichot was a member of the Inria local committee CORDI-S 2015.

J. Erhel: Master M2; title: Cours de modélisation et calcul scientifique; 12 hours; INSA, Rennes, France.

É. Canot: Master M2; title: TP de modélisation et calcul scientifique; 12 hours; INSA, Rennes, France.

PhD: L. Lenôtre, University of Rennes 1, November 2015, co-advisors A. Lejay (Inria Nancy) and G. Pichot, with J. Erhel.

PhD: S. Mansour, University of Rennes 1 with LIU and AUB (Beiruth, Lebanon), December 2015, co-advisors É. Canot, M. Muhieddine and N. Nassif.

PhD in progress: B. Delfino, University of Rennes 1, November 2015, co-advisors J.-R. de Dreuzy and J. Erhel.

PhD in progress: P.-M. Gibert, University of Lyon, October 2015, co-advisors D. Tromeur-Dervout and J. Erhel.

PhD in progress: M. ben Refifa, University of Tunis, October 2013, advisors R. Bouhlila and É. Canot.

PhD: A. Abeduwia, University of Littoral Côte d'Opale (Calais), Mathematics, December 2015. Reviewer J. Erhel.

PhD: C. Robinson, University of Rennes 1, Mathematics, December 2015. Committee member J. Erhel.

HdR: A. Beaudoin, University of Poitiers, Geosciences, December 2015. Committee member J. Erhel.

J. Erhel is co-author of two documents published in Intersitces, about modelling ice sheet evolution and estimating the contribution of Antartica and Greenland to sea-level rise. ,

J. Erhel was one of the coordinators of the Interstices publications devoted to climate change and COP21. She is co-author of the quiz on this subject.

J. Erhel gave a talk at the show "solutions COP21", held in Grand Palais, Paris, December 2015.

J. Erhel contributed to the new version of the spelling book for computer science (l'informatique de A à Z), published in November 2015. She was one of the coordinators of the project and she wrote several cards.

J. Erhel gave a talk in 180 seconds entitled "Enjeux pour les Eaux souterraines: Energie et Environnement", during the Inria scientific days, Nancy, June 2015.

J. Erhel presented the challenges for groundwater at the meeting between Inria and Sénat, January 2015.

L. Lenôtre gave a talk at "Lycée René Descartes de Rennes" on the diffusion phenomena and their underlying stochastic processes, February 2015.

L. Lenôtre participated in the "Fête de la science de Rennes" at Inria stand.