All Projects

Earthquakes can be devastating, both in terms of human and material damage. Although their existence has been known since the dawn of time, the physics of earthquakes is still poorly understood. Natural faults and earthquake characteristics are known to follow scaling power-law. The origin of this phenomenon has strong implications on the physical mechanisms driving slip events. However, it is not yet clear. The emergence of complexity can be related to the disorder of the system. Understanding if the observed complexity comes from the inherent complexity of the frictional motion or the system’s complexity is essential to better understand – and one day eventually predict – earthquakes. It has been shown numerically with a simple system without any disorder that resulting slip events follow a power-law distribution for the small events – like natural slip events – and a log-normal distribution for the larger ones. This project aims to study how adding disorder in this simple system will influence the transition between the power-law and the log-normal distribution of slip events. To do so, the student will use a finite element software developed in the lab (Akantu).

Supervisors:
Ferry Roxane Mathilde Suzanne, Jean-François Molinari

Extreme loads on solids lead to the formation of a multitude of cracks that propagate, branch and coalesce to form fragments. This process is called dynamic fragmentation. This process is of importance in many domains of engineering, where it is fundamental to predict the outcome of high velocity impacts or explosions. Often, one would like to extract statistics such as fragment size distribution. This project will feature experiments on object breaking into pieces to extract experimental statistics on fragments. The student will then use a finite element software (Akantu) to simulate crack propagation using different methods such as phase-field modelling of fracture or cohesive elements. The statistics obtained numerically will be compared to the experimental ones to highlight  the advantage and limitations of the different simulation methods.

Supervisors:
Thibault Ghesquière-Diérickx,  Shad  Durussel, Jean-François Molinari

This project seeks to combine the traditional finite element method with new machine-learning tools. The tentative goals are, first, to analyze new methods recently proposed in the literature  (e.g., the “Deep Ritz Method”) and, second, to adapt them to solve large-deformation problems in simple geometries.
 
Supervisors:
Joaquin Garcia Suarez, Jean-François Molinari
 
Reference: Yu, Bing. “The deep Ritz method: a deep learning-based numerical algorithm for solving variational problems.” Communications in Mathematics and Statistics 6.1 (2018): 1-12.
Wave propagation is a well-understood phenomenon, but gripping complexity can arise from it when it takes place within a layered medium, due to the continuous reflection, transmission and superposition of what originally may have been a coherent wavefront. On this topic, we propose two types of projects: (1) mostly theoretical ones (intended for students with a taste for applied mathematics) aimed at exploring how abstract algebra can help to unravel this complexity, (2) mostly numerical ones (for students interested in coding and simulations) aimed at implementing and testing new methods to simulate wave propagation in structured materials.
 
Supervisors:
Joaquin Garcia Suarez
 
References: email or come talk to Dr. Garcia-Suarez ([email protected], GC A2 515)
The disruptive work of Karniadakis and colleagues has created a new way to find approximate solution of ODEs and PDEs with relevance to physics and engineering. In this project, we would seek to survey the most efficient ways to train PINNs, and apply these findings to solve problems featuring stick-slip: . If the student is interested, a mathematical twist is also possible by focusing on analyzing the convergence and error estimates of PINNs.  
 
Supervisors:
Joaquin Garcia Suarez, Jean-François Molinari
 
Reference: Rucker, Cody, and Brittany A. Erickson. “Physics-informed deep learning of rate-and-state fault friction.” Computer Methods in Applied Mechanics and Engineering 430 (2024): 117211.
We are developing a new numerical method to solve non-linear problems in continuum mechanics and structural engineering. The new method has a strong geometrical component and comes with remarkable advantages: easy parallelization and adaptability. This project is ideal for students who want to put their coding (either Python or C++) and computational mechanics abilities (FEM) to work.
 
Supervisors:
Joaquin Garcia Suarez
 
Reference: Garcia-Suarez, J. “Phase-space iterative solvers.” arXiv preprint arXiv:2309.14031 (2023).

The impact of  a drop on a solid surface is a canonical problem in fluid mechanics of fundamental significance in numerous natural and industrial processes, such as ink-jet printing, aircraft icing and spray cooling. Recently we found out soft solids display a similar behavior when colliding with a rigid surface. Namely, the contact is not made on the tip, but on an annular radius, with air trapped in between. This project  will explore the scenario of highly viscous droplets and soft solids impacting on each other. The student will use the finite element software (Comsol Multiphysics) to simulate the dynamics using knowledge of both fluid and solid mechanics. Depending on the interest of the student the project will focus either on full 3D simulations to capture symmetry breaking or on axisymmetric ones to investigate the feasibility of using level-set or phase field simulation for droplet-air interface. 

Supervisors:
Jacopo Bilotto, Jean-François Molinari

Granular materials play a critical role in civil, chemical, and mechanical engineering, influencing sectors such as food processing, pharmaceuticals, soil mechanics, and mixing technologies. The Discrete Element Method (DEM) is a powerful tool for simulating the behavior of particle systems by applying Newton’s laws of motion and friction. However, one significant limitation of DEM is its ability to accurately model materials that undergo large deformations.
To address this challenge, a recent reduced-order model, developed by Professor Mollon, offers an innovative approach to extend traditional DEM codes into the deformable regime. This project aims to evaluate the feasibility and stability of this novel method.
The student will implement the model using Python, leveraging the Taichi and NumPy libraries to efficiently code the method and test various interaction laws. Depending on the student’s interests, the project can be tailored to focus on either high-performance GPU computing or the mechanics of specific test cases.

Supervisors:
Jacopo Bilotto, Jean-François Molinari