1 function alpha = fv_coercivity_bound(model)
2 %
function alpha = fv_coercivity_bound(model)
4 %
function determining the coercivity-alpha value of the implicit
5 %
operator of the current simulation and current parameters.
7 % currently only
for implicit diffusivity with linear dependence on
8 % k. In more general
case, extension must be performed.
10 % required fields of model:
11 % alpha_over_k : alpha is assumed to depend linearly on diffusivity k
12 % value can be determined by fv_estimate_coercivity_alpha.m
13 % and dividing it by k
15 % Bernard Haasdonk 29.8.2006
17 % only
if diffusion is discretized implicitly
18 alpha = model.alpha_over_k * model.k;