Skip to content
Snippets Groups Projects
Commit f968cc90 authored by David's avatar David
Browse files

test and fix vanGenuchten implementation

parent b95c0f32
No related branches found
No related tags found
No related merge requests found
......@@ -131,13 +131,12 @@ def vanGenuchten_S_sym(pc, n_index, alpha):
# derivative of S-pc relationship with respect to pc. This is needed for the
# construction of a analytic solution.
def vanGenuchten_S_prime(pc, n_index, alpha):
def vanGenuchten_S_prime_sym(pc, n_index, alpha):
# inverse capillary pressure-saturation-relationship
enumerator = -1.0*(alpha*(n_index - 1)*(alpha*pc)**(n_index - 1))
denominator = ((1 + (alpha*pc)**n_index)**((2*n_index - 1)/n_index))
return enumerator/denominator
def generate_Spc_dicts(
Spc_on_subdomains: tp.Dict[int, tp.Dict[str, tp.Dict[str, float]]]
)-> tp.Dict[str, tp.Dict[int, tp.Callable]]:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment