4 function kw = water_permeability(
this, glob, S, model)
6 kw = 0.01*ones(size(S));
9 function dkw = water_permeability_derivative(this, glob, S, model)
14 function ko = oil_permeability(this, glob, S, model)
19 function dko = oil_permeability_derivative(this, glob, S, model)
24 function pc = capillary_pressure(this, glob, S, model)
26 pc = model.tp_pb * sqrt((1-S)./S);
29 function dpc = capillary_pressure_derivative(this, glob, S, model)
31 dpc = model.tp_pb .* (-1/2)./sqrt(S.^3.-S.^4);
34 function ddpc = capillary_pressure_second_derivative(this, glob, S, model)
35 ddpc = model.tp_pb .* (1/4).*(S.^3-S.^4).^(-3/2).*(3.*S.^2-4.*S.^3);
36 % ddpc = model.tp_pb .* (1/2).*(4./S - 1./(S.^2-S.^3));
37 % ddpc = model.tp_pb .* (-1/2).*(-1.5.*S.^(0.5) + 2.*S)./(S.^(1.5)-S.^2).^2;
40 function c = injection_concentration(this, model)
42 c = model.tp_injection_c;
45 function s_under = lower_source(this, elemin, loc, grid)
47 glob = local2global(grid, elemin, loc, this);
52 D3 = ((X-0.8).^2 + (Y-0.5).^2) <= 0.01;
55 % s_under = zeros(size(X));
59 function s_above = upper_source(this, elemin, loc, grid)
61 glob = local2global(grid, elemin, loc, this);
66 D1 = ((X-0.5).^2 + (Y-0.8).^2) <= 0.01;
67 D2 = ((X-0.2).^2 + (Y-0.2).^2) <= 0.01;
68 s_above = 10 * D1 + 20 * D2;
70 % s_above = zeros(size(X));
74 function descr = default_descr(this, descr)
76 descr.tp_kw_factor = 0.15;
77 descr.tp_ko_factor = 1;
78 descr.tp_injection_c = 0.80;