Related papers: A Moving Boundary Flux Stabilization Method for Ca…
We develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for…
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must confront the challenge of controlling errors in the discrete divergence of the magnetic field. One approach that has been…
In this paper, we present an efficient numerical method to address a thermodynamically consistent gas flow model in porous media involving compressible gas and deformable rock. The accurate modeling of gas flow in porous media often poses…
We develop a linearly-scaling variant of the Force Coupling Method [K. Yeo and M. R. Maxey, J. Fluid Mech. 649, 205-231 (2010)] for computing hydrodynamic interactions among particles confined to a doubly-periodic geometry with either a…
In this paper we establish a stability barrier of a class of high-order Hermite-type discretization of 1D advection equations underlying the hybrid-variable (HV) and active flux (AF) methods. These methods seek numerical approximations to…
The constrained transport (CT) method reflects the state of the art numerical technique for preserving the divergence-free condition of magnetic field to machine accuracy in multi-dimensional MHD simulations performed with Godunov-type, or…
State redistribution is an algorithm that stabilizes cut cells for embedded boundary grid methods. This work extends the earlier algorithm in several important ways. First, state redistribution is extended to three spatial dimensions.…
A general derivation is proposed for several boundary conditions arisen in the lattice Boltzmann simulations of various physical problems. Pair-wise moment conservations are proposed to enforce the boundary conditions with given macroscopic…
We employ the summation-by-parts (SBP) framework to extend the recent domain-of-dependence (DoD) stabilization for cut cells to linear kinetic models in diffusion scaling. Numerical methods for these models are challenged by increased…
Distant spreading of primary lesions is modulated by the vascular dynamics of circulating tumor cells (CTCs) and their ability to establish metastatic niches. While the mechanisms regulating CTC homing in specific tissues are yet to be…
We develop arbitrarily high-order, stationarity-preserving stabilized finite element methods for multidimensional nonlinear hyperbolic balance laws on Cartesian grids. We aim at approximating all the steady states of the problem at hand,…
The Active Flux scheme is a Finite Volume scheme with additional degrees of freedom. It makes use of a continuous reconstruction and does not require a Riemann solver. An evolution operator is used for the additional degrees of freedom on…
In this paper, we develop a new free-stream preserving (FP) method for high-order upwind conservative finite-difference (FD) schemes on the curvilinear grids. This FP method is constrcuted by subtracting a reference cell-face flow state…
We present a new discretization for advection-diffusion problems with Robin boundary conditions on complex time-dependent domains. The method is based on second order cut cell finite volume methods introduced by Bochkov et al. to discretize…
The conventional no-slip boundary condition leads to a non-integrable stress singularity at a moving contact line. This makes numerical simulations challenging, especially when capillary effects are essential for the dynamics of the flow.…
The numerical simulation of multiphase flows involving dispersed components with large scale disparities, such as the collisions between millimeter-sized bubbles and micron-sized mineral particles in flotation, poses a significant…
A series of third- and fifth-order hybrid compact least-squares central weighted essentially non-oscillatory schemes are proposed and applied to curvilinear structured grids for the finite volume method. In smooth regions, compact…
Numerical methods for hyperbolic PDEs require stabilization. For linear acoustics, divergence-free vector fields should remain stationary, but classical Finite Difference methods add incompatible diffusion that dramatically restricts the…
High-order discontinuous Galerkin (DG) methods equipped with subcell finite-volume (FV) limiters provide an efficient framework for the simulation of nonlinear hyperbolic balance laws featuring shocks and complex flow structures. However,…
Two methods for solid body representation in flow simulations available in the Pencil Code are the immersed boundary method and overset grids. These methods are quite different in terms of computational cost, flexibility and numerical…