------------------------------
CHANGES TO .mmp FILE
------------------------------
$CAPILLARY_PRESSURE
------------------------------
--> OLD FORMAT
; pc_model -- pb (or alpha for van_genuchten)
6 5000.0
!Problems with this:
1) Pc cannot be independent of relative permeability values. Actually, they should be totally independent.
--> NEW FORMAT
; pc_model -- pb -- slr -- slm -- max_pc -- alpha_switch (if alpha_switch>0, Pb is taken as "alpha", where Pb = rho*g/alpha)
6 5000.0 0.0 1.0 1.0e12 0
!Notes:
Pc is no longer linked to relative permeability. It carries it's own values.
$PERMEABILITY_SATURATION
------------------------------
--> OLD FORMAT
; model -- min_sat -- max_sat -- exponent -- min_kr
6 0.0 1.0 2.0
44 0.0 1.0 0.5 1.0e-9
!Problems with this (these aren't big problems, but since we're reconfiguring Pc and k_rel anyway, may as well make the optimization):
1) Each phase should have k_rel entries that are independent of other phases.
--> NEW FORMAT
; model -- slr -- sgr -- exponent -- min_kr
6 0.25 0.025 2.0 1.0e-9
4 0.25 0.025 2.0 1.0e-9
! Notes.
1) Phases are totally independent. Each has it's own slr, sgr, exponent, minimum_k.
2) Exponents are the same for van Genuchten or Brooks-Corey (i.e., always greater than 1.0. thus for van Genuchten, m = 1 - 1/exponent)
3) No need to distinguish phase specific models.