Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
4967751 | Journal of Computational Physics | 2017 | 38 Pages |
Abstract
We present two new cell-centered nonlinear finite-volume methods for the heterogeneous, anisotropic diffusion problem. The schemes split the interfacial flux into harmonic and transversal components. Specifically, linear combinations of the transversal vector and the co-normal are used that lead to significant improvements in terms of the mesh-locking effects. The harmonic component of the flux is represented using a conventional monotone two-point flux approximation; the component along the parameterized direction is treated nonlinearly to satisfy either positivity of the solution as in [29], or the discrete maximum principle as in [9]. In order to make the method purely cell-centered, we derive a homogenization function that allows for seamless interpolation in the presence of heterogeneity following a strategy similar to [46]. The performance of the new schemes is compared with existing multi-point flux approximation methods [3], [5]. The robustness of the scheme with respect to the mesh-locking problem is demonstrated using several challenging test cases.
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Kirill M. Terekhov, Bradley T. Mallison, Hamdi A. Tchelepi,