Related papers: Fully Implicit Multidimensional Hybrid Upwind Sche…
Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by…
An optimisation scheme is developed to accurately represent the sub-grid scale forcing of a high dimensional chaotic ocean system. Using a simple parameterisation scheme, the velocity components of a 30km resolution shallow water ocean…
In the context of flow in porous media, up-scaling is the coarsening of a geological model and it is at the core of water resources research and reservoir simulation. An ideal up-scaling procedure preserves heterogeneities at different…
A conservative finite-volume framework, based on a collocated variable arrangement, for the simulation of flows at all speeds, applicable to incompressible, ideal-gas and real-gas fluids is proposed in conjunction with a fully-coupled…
The reliability of machine learning in multiscale physical systems depends on how physical structure is embedded into the learning process. We investigate this in the context of turbulent multiphase flows, focusing on the prediction of…
The pressure and flow statistics of Darcy flow through a random permeable medium are expressed in a form suitable for evaluation by the method of simulated annealing. There are several attractive aspects to using simulated annealing: (i)…
Accurate representation of interfaces and flux exchange is vital for coupled multiphysics simulations across a broad range of applications. Currently, coupling approaches are limited by the underlying discretization or to specific physical…
While various phase-field models have recently appeared for two-phase fluids with different densities, only some are known to be thermodynamically consistent, and practical stable schemes for their numerical simulation are lacking. In this…
We present a hybrid particle/grid approach for simulating incompressible fluids on collocated velocity grids. We interchangeably use particle and grid representations of transported quantities to balance efficiency and accuracy. A novel…
In this paper we carry out numerical analysis for a family of simplified gas transport models with hydrate formation and dissociation in subsurface, in equilibrium and non-equilibrium conditions. These models are adequate for simulation of…
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 develop second-order path-conservative central-upwind (PCCU) schemes for the hyperbolic shallow water linearized moment equations (HSWLME), which are an extension of standard depth-averaged models for free-surface flows. The proposed…
Developing high-order numerical schemes for two-phase flow in porous media that preserve key physical properties remains a significant challenge in numerical analysis. In this article, we propose a general framework to construct fully…
In this paper the projection hybrid FV/FE method presented in Busto et al. 2014 is extended to account for species transport equations. Furthermore, turbulent regimes are also considered thanks to the $k-\varepsilon$ model. Regarding the…
This study introduces a hybrid fluid simulation approach that integrates generative diffusion models with physics-based simulations, aiming at reducing the computational costs of flow simulations while still honoring all the physical…
We propose Coadjoint Orbit FLIP (CO-FLIP), a high order accurate, structure preserving fluid simulation method in the hybrid Eulerian-Lagrangian framework. We start with a Hamiltonian formulation of the incompressible Euler Equations, and…
In this paper, we construct and analyze new first- and second-order implicit-explicit (IMEX) schemes for the unsteady Navier-Stokes-Darcy model to describe the coupled free flow-porous media system, which is based on the scalar auxiliary…
The unified gas kinetic scheme (UGKS) is a direct modeling method based on the gas dynamical model on the mesh size and time step scales. With the implementation of particle transport and collision in a time-dependent flux function, the…
The direct-forcing immersed boundary method (DF-IBM) algorithm previously developed by the authors is extended by coupling the Navier-Stokes equations with the Newton-Euler equations for rigid body dynamics within the DF-IBM framework. This…
A numerical algorithm for solving mantle convection problems with strongly variable viscosity is presented. Equations for conservation of mass and momentum for highly viscous and incompressible fluids are solved iteratively by a multigrid…