Journées Scientifiques du GNR MoMaS
CIRM Marseille
du 2 au 4 novembre 2011



Mercredi 2 novembre


Session 1 (chaire A. Ern)

10h30 - 11h15   Conférence invitée: M. Ohlberger, Universitaet Muenster.
                           Error control and adaptivity for numerical multiscale methods.

11h15 - 11h40   M. Vohralik
                           A framework for robust a posteriori error control and adaptivity for unsteady nonlinear convection-diffusion problems.

11h40 - 12h05   D. Moreau
                           Hexahedral adaptive remeshing.

12h05 - 12h30   Présentation 3, thème C: C. Chainais
                           Simulation numérique de la corrosion d'acier dans les conditions de stockage souterrain.


12h30 - 13h45 déjeuner

Session 2 (chaire A. Bourgeat)

14h00 - 14h45   Conférence invitée: R. Helmig, University of Stuttgart.
                           Efficient Modeling of Flow and Transport in Porous Media Using Multiphysics and Multiscale Approaches.

14h45 - 15h10   J. Erhel et M. Kern
                           MoMaS: 10 years of reactive transport.

15h10 - 15h35   R. Herbin et F. Hubert
                           3D benchmark on anisotropic and heterogeneous diffusion.


15h35 - 16h00 pause café

Session 3 (chaire O. Le Maître)

16h00 - 16h25   Présentation partenaire: Andra: L. Loth
                           Contribution de l'Andra à la prospective dans le domaine des mathématiques appliquées.

16h25 - 16h50   Présentation partenaire IRSN: M. Dymitrowska and H. Amor
                           Numerical tools for monophasic and biphasic flows and solute transport in porous media developed by IRSN.

16h50 - 17h30   PACEN, MoMaS et Discussion

(17h45 - 19h30  Réunion du bureau du GNR)



19h30 - 20h30 diner

Session Posters (de 20h30 à 22h00)



Jeudi 3 novembre


Session 4 (chaire Ph. Montarnal)

09h00 - 09h45   Conférence invitée: F. Nobile, EPF-Lausanne.
                           Stochastic Polynomial approximation of PDEs with random coefficients.

09h45 - 10h10   N. Fajraoui, A. Younes et Th. Mara
                           Parameter identification and global sensitivity indices using sparse polynomial chaos expansions for reactive transport problems.

10h10 - 10h35   O. Le Maître and Ch. Soize
                           Non-parametric approach for modeling uncertainties in elliptic problems.


10h35 - 11h00 pause café

Session 5 (chaire G. Allaire)

11h00 - 11h45   Conférence invitée: T. Lelièvre, ENPC.
                           Numerical methods in molecular dynamics.

11h45 - 12h10   Ph. Ackerer
                           Multiscale approach for parameter identification.


12h30 - 14h00 déjeuner

Session 6 (chaire A. Mikelic)

14h00 - 14h45   Conférence invitée: A. Prohl, Universität Tübingen.
                           Structure-preserving discretizations of elektrokinetic flows.

14h45 - 15h10   J.-F. Dufrêche (avec GNR Paris)
                           To be announced.

15h10 - 15h35   G. Allaire
                           Ion transport in porous media: derivation of the macroscopic equations using upscaling and properties of the effective coefficients .

15h35 - 16h00   Présentation partenaire EdF: Fernandes, Granet, Mondoloni, Santamaria, Rupp
                           Study of the Thermo-Hydro-mechanical couplings around a nuclear waste tunnel with Syrthes and Code_Aster.


16h00 - 16h25 pause café

Session 7 (chaire S. Granet)

16h25 - 16h50   J.-J. Marigo
                           De l'initiation à la rupture: les propriétés générales et quelques exemples d'utilisation des modèles à gradient d'endommagement.

16h50 - 17h15   Présentation Partenaire CEA: P. Omnès
                           A posteriori error estimation and domain decomposition: towards efficient computations for waste repositories.

17h15 - 17h40   D. Stemmelen
                           Hydrodynamique et dispersion en milieu poreux étudiées par IRM.

17h40 - 18h05   M.-C. Néel
                           Non-Fickian effects and fractional models for diffusion in porous media.


19h30 - 21h00 diner



Vendredi 4 novembre


Session 8 (chaire A. Burnol)

09h00 - 09h45   Conférence invitée: B. Noetinger, IFPEN.
                           Modélisation avancée des réservoirs d'hydrocarbures souterrains, quelques problèmes ouverts.

09h45 - 10h10   M. Saad, F. Caro, K. Kahlil et B. Saad
                           Modèles diphasique bicomposant : analyse mathématique et numérique.

10h10 - 10h35   D. Hilhorst, K. Brenner et C. Cancès
                           Méthodes numériques pour les écoulements diphasiques en milieu poreux hétérogène.


10h35 - 11h00 pause café

Session 9 (chaire R. Herbin)

11h00 - 11h25   A. Pazdniakou et P.M. Adler
                           Modélisation des mélanges à deux phases à l'échelle du pore par méthode de Boltzmann.

11h25 - 11h45   I. Panfilova
                           Analyse hydrodynamique des régimes d'étalement et diffusion multi-composante des gaz en stockage souterrains de déchets radioactifs.

11h45 - 12h00   Présentation partenaire BRGM: A. Burnol et M. Parmentier
                           Problèmes de transfert de masse et d'énergie dans les milieux géologiques au BRGM: bilan et perspectives.

12h00 - 12h20   S. Granet
                           Benchmark multiphasique.


12h30 - 14h00 déjeuner

Session 10 (chaire A. Ern)

14h00 - 14h25   B. Amaziane et J. Jaffré
                           Modélisation mathématique et simulation numérique de la migration de Gaz dans le système d'un stockage souterrain de déchets radioactifs.

14h25 - 14h50   J. Bodin, M. Ghesmoune, A. Mikelic and M. Panfilov
                           Combination between the models of negative saturations and vector capillarity for multi-component fluids in tight porous media.

14h50 - 15h00   Clôture des journées: A. Ern

Abstracts of talks by invited speakers

  • R. Helmig
    Flow and transport processes in porous media including multiple fluid phases are the governing processes in a large variety of geological and technical systems. In general, these systems include processes of varying complexity occurring in different parts of the domain of interest. The different processes mostly also take place on different spatial and temporal scales. It is extremely challenging to model such systems in an adequate way accounting for the spatially varying and scale-dependent character of these processes. In this presentation, we give a brief overview of existing upscaling, multiscale, and multiphysics methods, and we present mathematical models and model formulations for multiphase flow in porous media including compositional and non-isothermal flow. Finally, we show simulation results for two-phase flow using a multiphysics method and a multiscale multiphysics algorithm.

  • T. Lelièvre
    Molecular dynamics is now a very widely used tool to study the matter at the molecular level. It is used in various fields, such as biology, chemistry or materials science. The aim is in particular to understand the relationships between the macroscopic properties of a molecular system and its atomistic features. For example, one would like to compute the constitutive relations for materials from molecular models, or predict the most likely conformations of a protein in a solvent from its amino acid sequence. One of the difficulty to reach this aim is related to timescales: the typical timescale of a molecular dynamics simulation is much smaller than the typical timescale at which the crucial events, from a macroscopic viewpoint, occur. This is related to the metastability of a molecular dynamics trajectory: the system stays for a very long time in some metastable state, before hopping to another one, and it is difficult to observe and simulate such rare events. An associated feature is the multimodality of the statistical ensemble (a probability measure) sampled by the molecular dynamics trajectories.
    Many methods have been proposed in the molecular dynamics community to deal with these difficulties, and we will focus on two prototypical ones for which a mathematical analysis gives useful insights. We will first present adaptive importance sampling techniques, which have been proposed to sample efficiently statistical ensembles. Then, we will propose a mathematical analysis of the parallel replica algorithm which has been introduced by A.F. Voter to generate efficiently metastable dynamics.

  • F. Nobile
    We consider the problem of numerically approximating statistical moments of the solution of a partial differential equation (PDE), whose input data (coefficients, forcing terms, boundary conditions, geometry, etc.) are uncertain and described by a finite or countably infinite number of random variables. This situation includes the case of infinite dimensional random fields suitably expanded in e.g Karhunen-Loève or Fourier expansions.
    We focus on polynomial chaos approximations of the solution with respect to the underlying random variables by either Galerkin projection; or collocation on sparse grids of Gauss points.
    We discuss in particular the proper choice of the polynomial space for linear elliptic PDEs with random diffusion coefficient. Numerical results showing the effectiveness and limitations of the approaches will be presented.
    We will present also some recent results on discrete projection on polynomial spaces from random evaluations.

  • B. Noetinger
    les coûts et les risques sans cesse croissants de l'exploitation des hydrocarbures souterrains justifient le recours à la modélisation numérique des écoulements. Ces modéles servent de réceptacle à la connaissance du géologue, du géophysicien et du spécialiste des écoulements. Ils permettent donc de prendre les décisions de forage ou d'exploitation les plus éclairées possibles. Des outils sont développés pour gérer le manque de connaissance exhaustive du réservoir, ce qui amène à la résolution de problèmes inverses et à gerer des incertitudes. Des problèmes fondamentaux de changement d'échelle sont aussi posés, en particulier dans le cas spécifique de roches fracturées. La présentation passera en revue ces problématiques, et quelques approches, tendances et résultats récents.

  • M. Ohlberger
    In this talk we address a posteriori error estimation and adaptivity for the heterogeneous multiscale method applied to elliptic and parabolic multi-scale problems. Moreover, inspired by the reduced basis approach, we present a new framework for an efficient treatment of heterogeneous multiscale problems. The approach is based on the idea of considering heterogeneous multiscale problems as parametrized partial differential equations where the parameters are smooth functions. We then construct in an offline phase, a suitable localized reduced basis that is used in an online phase to efficiently compute approximations of the multiscale problem by means of a Discontinuous Galerkin method on a coarse grid. We present our new approach for elliptic multiscale problems and discuss an a posteriori error estimate that is used in the construction process of the localized reduced basis. Numerical experiments are given to demonstrate the efficiency of the new approach.

  • A. Prohl
    Elektrokinetic flows of fluids dispersed with electric charges may be modelled by the Navier-Stokes equations coupled with classical semi-conductor equations. Corresponding solutions are non-negative, bounded charge concentrations, as well as an incompressible velocity field, which satisfy an energy- and an entropy principle.
    In the talk we propose discretizations which inherit some of these properties, and report on difficulties to construct-preserving discretizations. - This is joint work with M. Schmuck (Imperial College London).