This paper presents a new approach to verify accuracy of computational simulations. We develop mathematical theorems which can serve as robust a posteriori error estimation techniques to identify numerical pollution, check the performance of adaptive meshes, and verify numerical solutions. We demonstrate performance of this methodology on problems from flow thorough porous media. However, one c...

In this paper we consider a modification to Darcy equation by taking into account the dependence of viscosity on the pressure. We present a stabilized mixed formulation for the resulting governing equations. Equal-order interpolation for the velocity and pressure is considered, and shown to be stable (which is not the case under the classical mixed formulation). The proposed mixed formulation i...

This paper concerns the discretization of multiphase Darcy flows, in the case of heterogeneous anisotropic porous media and general 3D meshes used in practice to represent reservoir and basin geometries. An unconditionally coercive and symmetric vertex centred approach is introduced in this paper. This scheme extends the Vertex Approximate Gradient scheme (VAG), already introduced for single ph...

We present a method to efficiently simulate coronary perfusion in subject-specific models of the heart within clinically relevant time frames. Perfusion is modelled as a Darcy porous-media flow, where the permeability tensor is derived from homogenization of an explicit anatomical representation of the vasculature. To account for the disparity in length scales present in the vascular network, i...

We propose and analyze a third-order in time numerical algorithm for the Stokes-Darcy system for flows in fluid saturated karst aquifers. A novel third-order Adams-Moulton scheme is used for the discretization of the dissipative term while a third-order explicit Adams-Bashforth scheme is utilized for the time discretization of the interface term. The scheme is efficient in the sense that one ne...

Dissolution mechanisms in porous media may lead to unstable dissolution fronts (“wormholing”). It has been shown in the literature that Darcy-scale models may reproduce all the characteristics of such dissolution patterns. This paper considers the core-scale averaged behavior of these Darcy-scale dissolution models. The form of core-scale equations is discussed based on the volume averaging of ...

We propose a novel second order in time, decoupled and unconditionally stable numerical scheme for solving the Cahn-Hilliard-Darcy (CHD) system which models two-phase flow in porous medium or in a Hele-Shaw cell. The scheme is based on the ideas of second order convex-splitting for the Cahn-Hilliard equation and pressure-correction for the Darcy equation. We show that the scheme is uniquely sol...

Unconventional hydrocarbon resources from low-permeability formation, i.e., tight sands and shales, are currently received great attention because of their potential to supply the entire world with sufficient energy for the decades to come. In the past few years, as a result of industry-wide R&D effort, progresses are being made towards commercial development of gas and oil from such unconventi...

Linear stability of double diffusive convection of Darcy-Maxwell fluid with Soret and Dufour effects is investigated. The effects of the Soret and Dufour numbers, Lewis number, relaxation time and solutal Darcy Rayleigh number on the stationary and oscillatory convection are presented graphically. The Dufour number enhances the stability of Darcy-Maxwell fluid for stationary convection while it...

In this paper we investigate the accuracy of a Darcy-scale model for competitive adsorption in a porous medium through comparison with a smoothed particle hydrodynamics (SPH) pore-scale model. SPH is a Lagrangian, particle based modeling method which uses the particles as interpolation points to directly discretize the governing equations of the system. The models consider a binary system of co...