In this work, a stabilized formulation to solve the inductionless magnetohydrodynamic (MHD) problem using the finite element (FE) method is presented. The MHD problem couples the Navier-Stokes and a Darcy-type problem for the electric potential via Lorentz’s force in the momentum equation of the Navier-Stokes equations and the currents generated by the moving fluid in Ohm’s law. The key feature...

We consider a covolume method for a system of first order PDEs resulting from the mixed formulation of the variable-coefficient-matrix Poisson equation with the Neumann boundary condition. The system may be used to represent the Darcy law and the mass conservation law in anisotropic porous media flow. The velocity and pressure are approximated by the lowest order Raviart-Thomas space on rectang...

We use low order approximations, piecewise linear, continuous velocities and piecewise constant pressures to compute solutions to Brinkman’s equation of porous media flow, applying an edge stabilization term to avoid locking. In order to handle the limiting case of Darcy flow, when only the velocity component normal to the boundary can be prescribed, we impose the boundary conditions weakly usi...

The pressure drop in microchannels filled with porous media formed by square arrays of cylinders (micro-porous channels) is investigated. Combining the Brinkman equation and the existing models for permeability of regular arrays of cylinders, the pressure drop in the considered micro-porous channels is calculated theoretically. Soft lithography method is used to fabricate several Polydimethylsi...

The multigrid solution of coupled porous media and Stokes flow problems is considered. The Darcy equation as the saturated porous medium model is coupled to the Stokes equations by means of appropriate interface conditions. We focus on an efficient multigrid solution technique for the coupled problem, which is discretized by finite volumes on staggered grids, giving rise to a saddle point linea...

We introduce a new class of higher order mixed finite volume methods for elliptic problems. We start from the usual way of changing the given equation into a mixed system using the Darcy’s law, u = −K∇p. By integrating the system of equations with some judiciously chosen test spaces on each element, we define new mixed finite volume methods of higher order. We show that these new schemes are eq...

Homogenization techniques are used to derive a double porosity model for single phase flow in a reservoir with a preferred direction of fracture. The equations in the microscopic model are the usual ones derived from Darcy’s law in the fractures and matrix (rock). The permeability coefficients over the matrix domain are scaled, using a parameter ε, based on the fracture direction in the reservo...

An efficient nonlinear multigrid method for a mixed finite element method of the Darcy-Forchheimer model is constructed in this paper. A Peaceman-Rachford type iteration is used as a smoother to decouple the nonlinearity from the divergence constraint. The nonlinear equation can be solved element-wise with a closed formulae. The linear saddle point system for the constraint is reduced into a sy...

We present the development of a two-dimensional Mixed-Hybrid Finite Element (MHFE) model for the solution of the nonlinear equation of variably saturated ow in groundwater on unstructured triangular meshes. By this approach the Darcy velocity is approximated using lowest order Raviart-Thomas (RT 0) elements and is \exactly" mass-conserving. Hybridization is used to overcome the ill-conditioning...

In this paper we present the development of the macroscopic model describing the hydro-mechanical coupling of damaged porous media containing cracks or/and vugs, by using the asymptotic expansion method. The analysis starts at the mesoscopic scale at which we assume a generic microstructure and the validity of the Biot model in the micro-porous domain saturated by a fluid. In the crack/vug doma...