Creation of the Center for Modeling and Simulation in Strasbourg

posted Feb 28, 2013, 10:16 PM by Christophe Prud'homme   [ updated Feb 28, 2013, 10:17 PM ]

Since January 2013, the Center for Modeling and Simulation in Strasbourg has started its activities. Several collaboration projects have already started:
    • Blood flow simulation and blood flow rheology
  • High field magnets
  • Tumor detection
      • Tumor ablation
  • Moving fluid inclusion 
  • Cartography
  • Cabin aerothermy
For more information you can contact me ( or Vincent Huber (

HP-Feel++ awarded 60 000 000 core hours by PRACE on SUPERMUC

posted Feb 28, 2013, 10:06 PM by Christophe Prud'homme

It is my pleasure to announce that HP-Feel++ (High Performance Feel++) has been awarded 60 000 000 core hours on SUPERMUC (GAUSS@LRZ,Germany)[1] by the PRACE 6th regular call [2]. Among 88 projects submitted, 57 have been selected. HP-Feel++ has received fully the requested core hours. This project will run from March 2013 and March 2014.
HP-Feel++ is a collaboration between U. of Strasbourg(France), U. Joseph Fourier(Grenoble, France), CNRS and U. Coimbra(Portugal).

The HP-Feel++ project aims at developing two research applications that require now access of TIER-0 computing resources: blood flow rheology and high field magnets.

Although these domains are quite different they have been thoroughly developed for the past few years within the Feel++ project ( They share the same mathematical kernel that encompasses a large range of numerical methods to solve partial differential equations such as (i) arbitrary order continuous and discontinuous Galerkin methods in 1D, 2D and 3D, (ii) domain decomposition methods, (iii) fictitious domain methods, (iv) level-set methods or (iv) certified reduced basis methods.  These methods are developed and used easily using a domain specific language embedded in C++ mimicking the mathematical language associated to Galerkin methods. This language allows physicists, engineers and mathematicians to focus on the numerical methods as well the physics whilst it hides the computer science details (e.g. parallelism) or algebraic solvers and enables the user to ramp up very quickly from rapid prototyping numerical methods to large scale computations. Within this context, blood flow rheology and high field magnets are the two domains driving Feel++ developments.

In blood flow rheology, we are interested in simulating suspensions of red blood cells (RBC) in arteries and veins and in studying the fluid properties (i.e. the fluid apparent viscosity) either in healthy contexts (our current focus) or pathological contexts (in the longer term). Not only the RBC are deformable entities, arteries and veins deform also during blood pulse; in both cases fluid structure interaction modeling and simulations are required. We have developed two main alternatives to tackle these problems: (i) fluid structure interaction within the so-called Arbitrary Lagrangian Eulerian framework coupled with a fictitious domain method to handle the RBC and (ii) fluid structure interaction using level-set methods. In both cases, the computational and storage costs for realistic simulations require using the TIER-0 infrastructures.

As to high field magnets (i.e. magnetic intensity greater than 24T), they are being developed by a large scale equipment laboratory (Laboratoire national des champs magnetiques intenses) and they are accessible to the international scientific community through project calls. Studies range from solid physics to applied supra-conductivity and magneto-science. The design and optimisation of these high field magnets require the solution of large scale multi-physics (and mildly multi-scale) non-linear partial differential equations. Moreover to ensure a robust design, we need to assess uncertainties through quantile estimations and sensitivity analysis. The latter is built on the former as it requires hundred or thousands evaluations of the former. We have developed the so-called certified reduced basis in this context to reduce the computational cost within the uncertainty quantification and optimisation processes from millions of degrees of freedom to a few tens or hundreds. This huge computational gain requires however the acceptance of an intensive offline stage allowing to get the independence with respect to the costly (typically finite element) underlying models and which demands now the access to TIER-0 infrastructures.

2 Years Engineer Position @ Strasbourg

posted Oct 17, 2012, 2:20 AM by Christophe Prud'homme   [ updated Oct 17, 2012, 2:32 AM ]

A 2 years engineering position is open at université de Strasbourg starting in December 2012 or January 2013. The engineer will work at setting up the center of mathematical modeling and simulation of Strasbourg
We look for a candidate having a Phd in applied mathematics/scientific computing or having a strong applied mathematics/scientific computing content.
More detailed information about the required profile is available in the Job Positions page.

ANR VivaBrain

posted Aug 27, 2012, 2:33 AM by Christophe Prud'homme

The project VivaBrain has been selected by ANR for a 4 years funding. The main part of the funding goes toParitioning and Pressure coloring of a part of the brain vascular system one Phd starting no later than September 2013 and a three year engineer contract. The images on the left are the first simulation results obtained by V. Chabannes  on a part of the vascular system of the brain using Feel++. Feel++ is the selected tool for numerical simulation for VivaBrain.

Professor in Strasbourg

posted Jul 8, 2012, 1:58 PM by Christophe Prud'homme

I am now a Professor at Université de Strasbourg in the Mathematics Institute. I will be definitely joining the Partial Differential Equantion and Control Theory in September after spending 1 year and a half on leave there.
Several projects are awaiting me like
  • building a modeling and simulation center and developing  interactions with entreprises as well as other disciplines 
  • working on the virtual vascular brain project (Vivabrain) which got selected this year at the call for project MN of ANR
  • wotking on multi-disciplinary projects (tumor detecting, tumor ablation, tumor growth monitoring...)

Newcomers in the research group

posted Mar 16, 2012, 2:54 PM by Christophe Prud'homme

4 newcomers are joining my research group in March and April 2012
  • Djamal Diallo, ATER from U. de Strasbourg, on tumor thermal ablation. He will work with C. Essert, Amir Hussein and myself.
  • Anthony Lombard, Magister student from U. de Strasbourg, on diffuse optimal tomography for tumor detection. He will work jointly with M. Torregrossa, Z. Belhachmi and myself.
  • Ranine Tarabay, Master student from U. de Strasbourg, on blood flows in the brain. She will wotk jointly with V. Chabannes, M. Szopos and myself.
  • Marc Sidibe, Master student from U. de Grenoble, on level set methods. He will work with M. Ismail, V. Doyeux and myself.

Feel++ User Day and 1st HAMM Workshop

posted Mar 16, 2012, 2:49 PM by Christophe Prud'homme

I have organized the first Feel++ User Day on March 12 and the first HAMM workshop on March 13.
The slides are available for both workshops:

LaBeX IRMIA accepted

posted Feb 19, 2012, 10:23 PM by Christophe Prud'homme

The LabeX IRMIA (Laboratoire IRMA, Université de Strasbourg) was selected among other projects(about 50 accepted over  about 200 submitted) but the Investments for the Future funding program.
This project has two components that I contributed to 
  • develop high performance computing @ Strasbourg with the computer science ICPS from LSIIT
  • create an agency for mathematics in interaction with enterprises and society that would become the local relay for AMIES

Feel++ goes parallel

posted Feb 19, 2012, 8:46 AM by Christophe Prud'homme   [ updated Feb 19, 2012, 9:05 AM ]

Feel++  is now parallel. It uses MPI as the underlying message passing interface and PETSc. The image displays Paraview with the partitioning (done by Gmsh) of a cube into 6 subdomains (top left). Then we solve for the laplacian on this cube (solution bottom left). The 2D counterpart is available on the right views: top right displays the domain partitioning and bottom right the solution of a similar equation to the 3D one.

1-9 of 9