Related papers: A new formulation for coupled magma/mantle dynamic…
We present a new model for two phase Darcy flows in fractured media, in which fractures are modelled as submanifolds of codimension one with respect to the surrounding domain (matrix). Fractures can act as drains or as barriers, since…
Matched layers are commonly used in numerical simulations of wave propagation to model (semi-)infinite domains. Attenuation functions describe the damping in layers, and provide a matching of the wave impedance at the interface between the…
The energy balance of a partially molten rocky planet can be expressed as a non-linear diffusion equation using mixing length theory to quantify heat transport by both convection and mixing of the melt and solid phases. In this formulation…
The spherical Jeans equation is widely used to estimate the mass content of a stellar systems with apparent spherical symmetry. However, this method suffers from a degeneracy between the assumed mass density and the kinematic anisotropy…
A model of the three-dimensional rotating compressible Euler equations on the cubed sphere is presented. The model uses a mixed mimetic spectral element discretization which allows for the exact exchanges of kinetic, internal and potential…
We describe the implementation of a second-order accurate volume-of-fluid interface tracking algorithm in the open source finite element code ASPECT, which is designed to model convection in the Earth's mantle. This involves the solution of…
We will present qualitative and numerical results on a partial differential equation (PDE) system which models a certain fluid-structure dynamics. The wellposedness of this PDE model is established by means of constructing for it a…
A variational model of pressure-dependent plasticity employing a time-incremental setting is introduced. A novel formulation of the dissipation potential allows one to construct the condensed energy in a variationally consistent manner. For…
In this paper, we investigate a simple model where two, initially unconnected, flux systems are forced to interact in response to the imposed boundary driving by solving the non-ideal 3D MHD equations numerically. The reconnection rate of…
The molecular simulations solve the equation of motion of molecular systems, making 3D shapes of molecules four-dimensional by adding the time coordinate. These methods have a great potential in drug discovery because they can realistically…
A novel phase-flip model is proposed for thermodynamically consistent and computationally efficient description of spallation and cavitation in pure liquids within the framework of ideal hydrodynamics. Aiming at ultra-fast dynamic loads,…
This work proposes a novel model and numerical formulation for lubricated contact problems describing the mutual interaction between two deformable 3D solid bodies and an interposed fluid film. The solid bodies are consistently described…
Accurate numerical simulation of coupled fracture/fault deformation and fluid flow is crucial to the performance and safety assessment of many subsurface systems. In this work, we consider the discretization and enforcement of contact…
The solution of partial differential equations (PDEs) on complex domains often presents a significant computational challenge by requiring the generation of fitted meshes. The Diffuse Domain Method (DDM) is an alternative which reformulates…
Seismic imaging of the mantle has revealed large and small scale heterogeneities in the lower mantle; specifically structures known as large low shear velocity provinces (LLSVP) below Africa and the South Pacific. Most interpretations…
We introduce pointwise map smoothness via the Dirichlet energy into the functional map pipeline, and propose an algorithm for optimizing it efficiently, which leads to high-quality results in challenging settings. Specifically, we first…
We introduce a new micro/macro decomposition of collisional kinetic equations which naturally incorporates the exact space boundary conditions. The idea is to write the distribution fonction $f$ in all its domain as the sum of a Maxwellian…
Primal-dual methods for solving convex optimization problems with functional constraints often exhibit a distinct two-stage behavior. Initially, they converge towards a solution at a sublinear rate. Then, after a certain point, the method…
The immersed boundary method is a mathematical formulation and numerical method for solving fluid-structure interaction problems. For many biological problems, such as models that include the cell membrane, the immersed structure is a…
An efficient numerical approach to equilibrium properties of strongly coupled systems which include a subsystem of fermionic quantum particles and a subsystem of classical particles is presented. It uses an improved path integral…