Publica
Hier finden Sie wissenschaftliche Publikationen aus den FraunhoferInstituten. On convergence of certain finite volume difference discretizations for 1D poroelasticity interface problems
 Numerical methods for partial differential equations 23 (2007), No.3, pp.652671 ISSN: 0749159X 

 English 
 Journal Article 
 Fraunhofer ITWM () 
Abstract
In the article two finite difference schemes for the ID poroelasticity equations (Biot model) with discontinuous coefficients are derived, analyzed, and numerically tested. A recent discretization [Gaspar et al., Appl Numer Math 44 (2003), 487506] of these equations with constant coefficients on a staggered grid is used as a basis. Special attention is given to the interfaces and as a result a scheme with harmonic averaging of the coefficients is derived. Convergence rate of O(h (3)/(2)) in a discrete H1 norm for both the pressure and the displacement is established in the case of an arbitrary position of the interface. Further, rate of 0 (h 2) is proven for the case when the interface coincides with a grid node. Following an approach applied to secondorder elliptic equations in [Ewing et al., SIAM J Sci Comp 23(4) (2001), 13341350] we derive a modified and more accurate discretization that gives secondorder convergence of the fluid velocity and the stress of the solid. Finally, numerical experiments of model problems that confirm the theoretical considerations are presented.