Related papers: A fast mesh-free boundary integral method for two-…
An existing phase-field model of two immiscible fluids with a single soluble surfactant present is discussed in detail. We analyze the well-posedness of the model and provide strong evidence that it is mathematically ill-posed for a large…
We study the pore-scale transport of a conservative scalar forming an advancing mixing front, which can be re-interpreted to predict instantaneous mixing-limited bimolecular reactions. We investigate this using a set of two-dimensional,…
The stability of most surface-tension-driven interfacial flow simulations is governed by the capillary time-step constraint. This concerns particularly small-scale flows and, more generally, highly-resolved liquid-gas simulations with…
We propose a numerical method for fluid deformable surfaces governed by surface Stokes flow and Helfrich bending energy under active growth, aiming to model shape evolution of the epithelium in developmental processes. To prevent…
In this work, we develop a new algorithm to solve large-scale incompressible time-dependent fluid--structure interaction (FSI) problems using a matrix-free finite element method in arbitrary Lagrangian--Eulerian (ALE) frame of reference. We…
Pure advection of a conservative scalar is relevant to several applications including two-phase flow. Successful numerical schemes must capture the sharp interface between the phases while maintaining a smooth (wrinkle-free) interfacial…
We present a novel framework to deal with static and moving immersed boundaries (IB). In this strategy, called Volume-Filtering Immersed Boundary (VFIB) method, transport equations are derived by filtering the Navier-Stokes equations and…
Peskin's Immersed Boundary (IB) model and method are among one of the most important modeling tools and numerical methods. The IB method has been known to be first order accurate in the velocity. However, almost no rigorous theoretical…
Surfactants are critical in natural processes and engineering, but measuring their concentrations in non-equilibrium conditions and in the presence of flow is difficult. Therefore, computational methods are a key tool for improving our…
We consider primal-dual mixed finite element methods for the advection--diffusion equation. For the primal variable we use standard continuous finite element space and for the flux we use the Raviart-Thomas space. We prove optimal a priori…
We present a robust immersed boundary (IB) method for high density ratio multiphase flows that is capable of modeling complex wave-structure interaction (WSI) problems arising in marine and coastal engineering applications. The IB/WSI…
We present a parallel-scalable method for simulating non-dilute suspensions of deformable particles immersed in Stokesian fluid in three dimensions. A critical component in these simulations is robust and accurate collision handling. This…
This paper is concerned with numerical solution of transport problems in heterogeneous porous media. A semi-discrete continuous-in-time formulation of the linear advection-diffusion equation is obtained by using a mixed hybrid finite…
In the present work, we propose a novel hybrid explicit jump immersed interface approach in conjunction with a higher order compact (HOC) scheme for simulating transient complex flows governed by the streamfunction-vorticity…
Simulation of unsteady creeping flows in complex geometries has traditionally required the use of a time-stepping procedure, which is typically costly and unscalable. To reduce the cost and allow for computations at much larger scales, we…
We present a novel platform for the large-scale simulation of fibrous structures immersed in a Stokesian fluid and evolving under confinement or in free-space. One of the main motivations for this work is to study the dynamics of fiber…
We present the first convergence proof for an iso-parametric finite element discretization of two-phase Stokes flow in $\Omega \subset \mathbb{R}^d$, $d=2,3$, with interface dynamics governed by mean curvature. The proof relies on a crucial…
The incompressible Stokes equations can classically be recast in a boundary integral (BI) representation, which provides a general method to solve low-Reynolds number problems analytically and computationally. Alternatively, one can solve…
We study the combined effects of natural convection and rotation on the dissolution of a solute in a solvent-filled circular cylinder. The density of the fluid increases with the increasing concentration of the dissolved solute, and we…
We present a numerical method for interface-resolved simulations of evaporating two-fluid flows based on the volume-of-fluid (VoF) method. The method has been implemented in an efficient FFT-based two-fluid Navier-Stokes solver, using an…