29 research outputs found
Higher resolution total velocity Vt and Va finite-volume formulations on cell-centred structured and unstructured grids
Novel cell-centred finite-volume formulations are presented for incompressible and immiscible two-phase flow with both gravity and capillary pressure effects on structured and unstructured grids. The Darcy-flux is approximated by a control-volume distributed multipoint flux approximation (CVD-MPFA) coupled with a higher resolution approximation for convective transport. The CVD-MPFA method is used for Darcy-flux approximation involving pressure, gravity, and capillary pressure flux operators. Two IMPES formulations for coupling the pressure equation with fluid transport are presented. The first is based on the classical total velocity Vt fractional flow (Buckley Leverett) formulation, and the second is based on a more recent Va formulation. The CVD-MPFA method is employed for both Vt and Va formulations. The advantages of both coupled formulations are contrasted. The methods are tested on a range of structured and unstructured quadrilateral and triangular grids. The tests show that the resulting methods are found to be comparable for a number of classical cases, including channel flow problems. However, when gravity is present, flow regimes are identified where the Va formulation becomes locally unstable, in contrast to the total velocity formulation. The test cases also show the advantages of the higher resolution method compared to standard first-order single-point upstream weighting
On the upstream mobility scheme for two-phase flow in porous media
When neglecting capillarity, two-phase incompressible flow in porous media is
modelled as a scalar nonlinear hyperbolic conservation law. A change in the
rock type results in a change of the flux function. Discretizing in
one-dimensional with a finite volume method, we investigate two numerical
fluxes, an extension of the Godunov flux and the upstream mobility flux, the
latter being widely used in hydrogeology and petroleum engineering. Then, in
the case of a changing rock type, one can give examples when the upstream
mobility flux does not give the right answer.Comment: A preprint to be published in Computational Geoscience
A theory of -dissipative solvers for scalar conservation laws with discontinuous flux
We propose a general framework for the study of contractive semigroups
of solutions to conservation laws with discontinuous flux. Developing the ideas
of a number of preceding works we claim that the whole admissibility issue is
reduced to the selection of a family of "elementary solutions", which are
certain piecewise constant stationary weak solutions. We refer to such a family
as a "germ". It is well known that (CL) admits many different contractive
semigroups, some of which reflects different physical applications. We revisit
a number of the existing admissibility (or entropy) conditions and identify the
germs that underly these conditions. We devote specific attention to the
anishing viscosity" germ, which is a way to express the "-condition" of
Diehl. For any given germ, we formulate "germ-based" admissibility conditions
in the form of a trace condition on the flux discontinuity line (in the
spirit of Vol'pert) and in the form of a family of global entropy inequalities
(following Kruzhkov and Carrillo). We characterize those germs that lead to the
-contraction property for the associated admissible solutions. Our
approach offers a streamlined and unifying perspective on many of the known
entropy conditions, making it possible to recover earlier uniqueness results
under weaker conditions than before, and to provide new results for other less
studied problems. Several strategies for proving the existence of admissible
solutions are discussed, and existence results are given for fluxes satisfying
some additional conditions. These are based on convergence results either for
the vanishing viscosity method (with standard viscosity or with specific
viscosities "adapted" to the choice of a germ), or for specific germ-adapted
finite volume schemes
Notation
Abstract. The swelling and shrinkage of biological tissues are modelled by a four-component mixture theory. This theory results in a coupled system of nonlinear parabolic differential equations together with an algebraic constraint for electroneutrality. In this model, it is desirable to obtain accurate approximations of the fluid flow and ions flow. Such accurate approximations can be determined by the mixed finite element method. The solid displacement, fluid and ions flow and electro-chemical potentials are taken as degrees of freedom. In this article the lowest-order mixed method is discussed. This results into a first-order nonlinear algebraic equation with an indefinite coefficient matrix. The hybridization technique is then used to reduce the list of degrees of freedom and to speed up the numerical computation. The mixed hybrid finite element method is then validated for small deformations using the analytical solutions for one-dimensional confined consolidation and swelling. Two-dimensional results are shown in a swelling cylindrical hydrogel sample. 1991 Mathematics Subject Classification. 35K55, 35M10, 65F10, 65M60. The dates will be set by the publisher
Numerical simulation of deformations and electrical potentials in a cartilage substitute
Cartilage exhibits a swelling and shrinking behaviour that influences the function of the cells inside the tissue. This behaviour is caused by mechanical, chemical and electrical loads. It is described by the electrochemomechanical mixture theory, in which the tissue is represented by four components: a charged porous solid, a fluid, cations and anions. By distinguishing between the cations and anions, electrical phenomena can be modelled. This mixture theory is verified by fitting the deformations and the electrical potentials in a uniaxial confined swelling and compression experiment to a mixed finite element simulation. The fitted stiffness, permeability, diffusion coefficients, and osmotic coefficients are in the same range as reported in literature