Related papers: Smooth Implicit Hybrid Upwinding for Compositional…
Many numerical schemes for hyperbolic systems require a piecewise polynomial reconstruction of the cell averaged values, and to simulate perturbed steady states accurately we require a so called 'well balanced' reconstruction scheme. For…
We study three different time integration methods for a dynamic pore network model for immiscible two-phase flow in porous media. Considered are two explicit methods, the forward Euler and midpoint methods, and a new semi-implicit method…
This paper proposes a computationally efficient algorithm to address the Full-Waveform Inversion (FWI) problem with a Total Variation (TV) constraint, designed to accurately reconstruct subsurface properties from seismic data. FWI, as an…
This paper presents the numerical solution of immiscible two-phase flows in porous media, obtained by a first-order finite element method equipped with mass-lumping and flux up-winding. The unknowns are the physical phase pressure and phase…
Coupled nonlinear system of reaction-diffusion equations describing multi-component (species) interactions with heterogeneous coefficients is considered. Finite volume method based approximation for the space is used to construct…
The goal of this paper is to study approaches to bridge the gap between first-order and second-order type methods for composite convex programs. Our key observations are: i) Many well-known operator splitting methods, such as…
When modeling astrophysical fluid flows, it is often appropriate to discard the canonical magnetohydrodynamic approximation thereby freeing the magnetic field to diffuse with respect to the bulk velocity field. As a consequence, however,…
We develop a new finite volume method using unstructured mesh-vertex grids for coupled systems modeling shallow water flows and solute transport over complex bottom topography. Novel well-balanced positivity preserving discretization…
The complexity of binary droplet collisions increases for the collision of immiscible liquids with the occurrence of triple lines and thin encapsulating films. The Volume of Fluid (VOF) method is extended with an efficient interface…
We propose a variational framework for the resolution of a non-hydrostatic Saint-Venant type model with bottom topography. This model is a shallow water type approximation of the incompressible Euler system with free surface and slightly…
We propose the complex-plane generalization of a powerful algebraic sequence acceleration algorithm, the Method of Weighted Averages (MWA), to guarantee exponential-cum-algebraic convergence of Fourier and Fourier-Hankel (F-H) integral…
We study a one dimensional model for two-phase flows in heterogeneous media, in which the capillary pressure functions can be discontinuous with respect to space. We first give a model, leading to a system of degenerated non-linear…
In this work we present a mass conservative numerical scheme for two-phase flow in porous media. The model for flow consists on two fully coupled, non-linear equations: a degenerate parabolic equation and an elliptic equation. The proposed…
An open-sourced multiphase Darcy-Brinkman approach is proposed to simulate two-phase flow in hybrid systems containing both solid-free regions and porous matrices. This micro-continuum model is rooted in elementary physics and volume…
We develop a new second-order unstaggered path-conservative central-upwind (PCCU) scheme for ideal and shallow water magnetohydrodynamics (MHD) equations. The new scheme possesses several important properties: it locally preserves the…
We propose and analyze a combined finite volume--nonconforming finite element scheme on general meshes to simulate the two compressible phase flow in porous media. The diffusion term, which can be anisotropic and heterogeneous, is…
In this work, we design and analyze a Hybrid High-Order (HHO) discretization method for incompressible flows of non-Newtonian fluids with power-like convective behaviour. We work under general assumptions on the viscosity and convection…
An issue for molecular dynamics simulations is that events of interest often involve timescales that are much longer than the simulation time step, which is set by the fastest timescales of the model. Because of this timescale separation,…
The majority of available numerical algorithms for interfacial two-phase flows either treat both fluid phases as incompressible (constant density) or treat both phases as compressible (variable density). This presents a limitation for the…
In this work we introduce a novel semi-implicit structure-preserving finite-volume/finite-difference scheme for the viscous and resistive equations of magnetohydrodynamics (MHD) based on an appropriate 3-split of the governing PDE system,…