Related papers: Boundary elements method for microfluidic two-phas…
In this paper, we will provide the the finite element method for the electro-osmotic flow in micro-channels, in which a convection-diffusion type equation is given for the charge density $\rho^e$. A time-discrete method based on the…
We introduce an adaptive superconvergent finite element method for a class of mixed formulations to solve partial differential equations involving a diffusion term. It combines a superconvergent postprocessing technique for the primal…
The Boundary Element Method (BEM) is implemented using piecewise linear elements to solve the two-dimensional Dirichlet problem for Laplace's equation posed on a disk. A benefit of the BEM as opposed to many other numerical solution…
A low-Mach-number flow, in the laminar regime, has intrinsically two characteristic spatial scales for a given time scale, or two characteristic temporal scales for a given spatial scale, and these dual scales are very different due to the…
A two-phase, low-Mach-number flow solver is created and verified for variable-density liquid and gas with phase change. The interface is sharply captured using a split Volume-of-Fluid method generalized for a non-divergence-free liquid…
Microfluidic devices offer unique opportunities to directly observe multiphase flow in porous media. However, as a direct representation of flow in geological pore networks, conventional microfluidics face several challenges. One is that…
Double emulsion formation in a hierarchical flow-focusing channel is systematically investigated using a free energy ternary lattice Boltzmann model. A three dimensional formation regime diagram is constructed based on the capillary numbers…
We develop a spectral element lattice Boltzmann method (SELBM) with the flux bounce-back (FBB) scheme, to enable accurate simulations of single-phase fluid dynamics in unstructured mesh. We adopt an Eulerian description of the streaming…
Despite its numerical challenges, finite element method is used to compute viscous fluid flow. A consensus on the cause of numerical problems has been reached; however, general algorithms---allowing a robust and accurate simulation for any…
This paper proposes a new method, in the frequency domain, to define absorbing boundary conditions for general two-dimensional problems. The main feature of the method is that it can obtain boundary conditions from the discretized equations…
The implementation of boundary conditions is among the most challenging parts of modeling fluid flow through channels and complex media. Here, we show that the existing methods to deal with liquid-wall interactions using multicomponent…
We present a hybrid spectral element-Fourier spectral method for solving the coupled system of Navier-Stokes and Cahn-Hilliard equations to simulate wall-bounded two-phase flows in a three-dimensional domain which is homogeneous in at least…
Trajectories of a buoyant spherical solid particle in a linear shear flow were investigated at low Reynolds numbers. A two-dimensional CFD analysis was performed to simulate the solid-fluid flows. Our numerical model, the discrete phase…
In this paper we study a system of advection-diffusion equations in a bulk domain coupled to an advection-diffusion equation on an embedded surface. Such systems of coupled partial differential equations arise in, for example, the modeling…
We present a novel approach for solving the shallow water equations using a discontinuous Galerkin spectral element method. The method we propose has three main features. First, it enjoys a discrete well-balanced property, in a spirit…
The Richards equation, a nonlinear elliptic parabolic equation, is widely used to model infiltration in porous media. We develop a finite element method for solving the Richards equation by introducing a new bounded auxiliary variable to…
A coupled boundary spectral element method (BSEM) and spectral element method (SEM) formulation for the propagation of small-amplitude water waves over variable bathymetries is presented in this work. The wave model is based on the…
Central finite difference schemes have long been avoided in the context of two-phase flows for the advection of the phase indicator function due to numerical overshoots and undershoots associated with their dispersion errors. We will show…
An algorithm is proposed for generalized mean curvature flow of closed two-dimensional surfaces, which include inverse mean curvature flow, powers of mean and inverse mean curvature flow, etc. Error estimates are proven for semi- and full…
This paper presents a numerical study of immiscible, compressible two-phase flows in porous media, that takes into account heterogeneity, gravity, anisotropy, and injection/production wells. We formulate a fully implicit stable…