Research

My research is focused on problems with moving interfaces in geophysics, mechanical and biomedical engineering. Most of the challenges related to the numerical solution of moving boundary problems can be bypassed by using the phase-field method, which reformulates the moving boundary problem as a set of PDEs posed on a known and fixed computational domain. The phase-field method introduces a new variable, the phase field, which identifies each phase (or domain) and smoothly transitions on the interface. I develop new accurate, robust, and efficient computational methods to study a wide range of problems of interest in geophysics, mechanics and biomechanics.

Cell migration

Both individual and collective cell migration regulate key biological functions such as tissue formation, wound healing, or cancer metastasis. A motile cell is an outstanding example of a moving boundary problem. Cell motion results from complex intra- and extra-cellular mechanochemical interactions. The phase-field method allows me to define the cell shape, localize the evolution equations of different compounds to moving volumes (cell's interior and exterior) and surfaces (cell's membrane), and control the cellular dynamics of multicellular systems through an energy functional. The pictures show individual cell chemotaxis in a fibrous environment (left) and cell co-attraction in non-confluent multicellular systems (right).

Snow hydrology

Snow is a heterogeneous porous material composed of ice, liquid water, and air. Complex processes such as preferential meltwater infiltration or ice piping affect the hydraulic, thermal, and mechanical properties of the snow. I study these processes by using the phase-field method. I model the phase transitions of the different phases (ice, liquid water, and air) at the microscopic scale to reproduce snow metamorphism. I also derived a model for meltwater infiltration at the macroscopic scale (~meter), which may lead to better prediction of meltwater infiltration and improved management of water resources. In the picture, water vapor concentration, ice, liquid water, and temperature during a snow melting event.

Soil carbon

Soil organic matter represents the second largest reservoir of carbon in Earth. Microbial respiration, which is mediated by soil moisture, temperature, and soil structure, dictates CH4 and CO2 fluxes into the atmosphere. I develop a Darcy-scale model that couples soil biogeochemical processes under a reactive-transport framework to capture soil decomposition and gas generation. The model accounts for soil heterogeneity, organic matter heterogeneity, and water infiltration. The figure in the left shows the water saturation, CO2 concentration, and O2 concentration caused by aerobic respiration during a drainage event. Red circles represent soil aggregates (hotspots for microbial activity). In the right, CO2 flux and accumulated CO2 released to the atmosphere.

Cancer growth

The phase-field method may be used to model different processes involved in cancer growth. In particular, I have studied angiogenesis and glioma growth. Angiogenesis is the growth of new capillaries induced by tumors. The new capillaries provide an extra supply of oxygen to the tumor, which results essential to continue the tumor growth. Capillaries are localized through the phase-field method. Glioma is a highly malignant tumor that forms in the glial cells of the brain. I used the phase-field method to localize the dynamics of tumor cells, vasculature, and angiogenic factors to the brain geometry. The picture shows the time evolution of the concentration of normoxic tumor cells.

Multiphase flow

In phase-transition problems, interfaces are intrinsic to the problem. I have developed the phase-field formulation and the numerical implementation for solidification of binary mixtures, phase separation of ternary mixtures, and methane hydrate formation and dissociation. The Ouzo effect refers to the spontaneous nucleation of oil microdroplets in an ethanol-oil-water mixture, in which oil-water miscibility depends on the ethanol concentration. Methane hydrate formation and dissociation constitutes a two-component (methane and water) three-phase (solid, liquid, and gas) flow problem. Hydrate natural deposits could represent an important energy source, which explains why the computational study of the extraction process is crucial. The pictures show the nucleation of oil microdroplets in an ethanol-oil-water mixture (left) and the initial times of methane hydrate formation (right).

Fracture dynamics

I have studied brittle fracture in elastic solids. The sharp crack discontinuity is modeled through the phase field, which smoothly transitions between the fully broken phase and the intact material phase. The momentum balance equation for the solid and the crack evolution equation are derived from a free energy functional. These equations can be coupled to a mass continuity equation for the fluid flow, which permits to model fluid-driven fracture in elastic media. In the picture, single-edge notched shear test when crack propagation is allowed under tension and compression loading states (isotropic formulation).