search

TDycores-Project / TDycore

BSD 2-Clause "Simplified" License
4 stars 0 forks source link

Debug mismatch between MPFA-O TPF and PFLOTRAN #130

Closed bishtgautam closed 3 years ago

bishtgautam commented 3 years ago

MPFA-O can be reduced to two-point flux approximation via:

-tdy_mpfao_gmatrix_method          MPFAO_GMATRIX_TPF \
-tdy_mpfao_boundary_condition_type MPFAO_DIRICHLET_BC

Debug why MPFA-O TPF does not match PFLTORAN.

bishtgautam commented 3 years ago

The contribution of the gravity term is incorrectly accounted for in the current implementation. The MPFA-O should only account for pressure and the gravity term should be accounted for as two-point flux (similar to PFLOTRAN's implementation). Here is a relevant reference https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2019WR025384?af=R&ai=1gvoi&mi=3ricys