Related papers: A Complex-boundary Treatment Method for Finite Vol…
Numerical simulation of compressible fluid flows is performed using the Euler equations. They include the scalar advection equation for the density, the vector advection equation for the velocity and a given pressure dependence on the…
The paper develops a method for the numerical simulation of a free-surface flow of incompressible viscous fluid around a streamlined body. The body is a rigid stationary construction partially submerged in the fluid. The application we are…
This paper presents a new and efficient nodal scheme for cell-centered compressible flows in Lagrangian formulation. A single pressure and the velocity at cell vertex are computed by the nodal solver in a least-squares sense, and both…
A general, compact way of achieving second-order in finite-volume numerical methods is to perform a MUSCL-like, piecewise linear reconstruction of flow properties at each cell interface. To avoid the surge of spurious oscillations in the…
A high order finite difference method is proposed for unstructured meshes to simulate compressible inviscid/viscous flows with/without discontinuities. In this method, based on the strong form equation, the divergence of the flux on each…
Numerical solutions of 2-D steady incompressible flow inside a triangular cavity are presented. For the purpose of comparing our results with several different triangular cavity studies with different triangle geometries, a general triangle…
We propose a mathematical model for fluids in multiphase flows in order to establish a solid theoretical foundation for the study of their complex topology, large geometric deformations, and topological changes such as merging. Our modeling…
In this work we are interested in dealing with single-phase flows in fractured porous media for underground processes. We focus our attention on domains where the presence of faults, with thickness several orders of magnitude smaller than…
We consider the numerical approximation of single phase flow in porous media by a mixed finite element method with mass lumping. Our work extends previous results of Wheeler and Yotov, who showed that mass lumping together with an…
A fourth-order finite volume embedded boundary (EB) method is presented for the unsteady Stokes equations. The algorithm represents complex geometries on a Cartesian grid using EB, employing a technique to mitigate the "small cut-cell"…
We present a novel method for calculating interface curvature on 3D unstructured meshes from piecewise-linear interface reconstructions typically generated in the volume of fluid method. Interface curvature is a necessary quantity to…
Defined mathematically as critical points of surface area subject to a volume constraint, constant mean curvatures (CMC) surfaces are idealizations of interfaces occurring between two immiscible fluids. Their behavior elucidates phenomena…
We present a computational framework for modeling large-scale particle-laden flows in complex domains with the goal of enabling simulations in medical-image derived patient specific geometries. The framework is based on a volume-filtered…
Buoyancy-induced (Rayleigh-Benard) convection of a fluid between two horizontal plates is a central paradigm for studying the transition to complex spatiotemporal dynamics in sustained nonequilibrium systems. To improve the analysis of…
A finite difference scheme is used to develop a numerical method to solve the flow of an unbounded viscoelastic fluid with zero to moderate inertia around a prolate spheroidal particle. The equations are written in prolate spheroidal…
To numerically solve a generic elliptic equation on two-dimensional domains with rectangular Cartesian grids, we propose a cut-cell geometric multigrid method that features (1) general algorithmic steps that apply to two-dimensional…
The treatment of complex geometries in Computational Fluid Dynamics applications is a challenging endeavor, which immersed boundary and cut-cell techniques can significantly simplify by alleviating the meshing process required by…
Numerous infinite dimensional dynamical systems arising in different fields have been shown to exhibit a gradient flow structure in the Wasserstein space. We construct Two Point Flux Approximation Finite Volume schemes discretizing such…
In a companion study \cite{patterson2020computing2D}, we present a numerical method for simulating 2D viscous flow through an open compliant closed channel, drive by pressure gradient. We consider the highly viscous regime, where fluid…
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…