rbmatlab  1.16.09
 All Classes Namespaces Files Functions Variables Modules Pages
fv_coercivity_bound.m
1 function alpha = fv_coercivity_bound(model)
2 %function alpha = fv_coercivity_bound(model)
3 %
4 % function determining the coercivity-alpha value of the implicit
5 % operator of the current simulation and current parameters.
6 %
7 % currently only for implicit diffusivity with linear dependence on
8 % k. In more general case, extension must be performed.
9 %
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
14 
15 % Bernard Haasdonk 29.8.2006
16 
17 % only if diffusion is discretized implicitly
18 alpha = model.alpha_over_k * model.k;
19 
20 
21 
22