MichelSimple.m

classdef MichelSimple < TwoPhaseData.Interface

  methods
    function kw = water_permeability(this, glob, S, model)

      kw = model.tp_kw_factor * ones(size(S));
    end

    function dkw = water_permeability_derivative(this, glob, S, model)

      dkw = 0;
    end

    function ko = oil_permeability(this, glob, S, model)

      ko = model.tp_ko_factor * ones(size(S));
    end

    function dko = oil_permeability_derivative(this, glob, S, model)

      dko = 0;
    end

    function pc = capillary_pressure(this, glob, S, model)

      pc = model.tp_pb * S;
    end

    function dpc = capillary_pressure_derivative(this, glob, S, model)

      dpc = model.tp_pb * ones(size(S));
    end

    function ddpc = capillary_pressure_second_derivative(this, glob, S, model)

      ddpc = zeros(size(S));
    end

    function c = injection_concentration(this, model)

      c = model.tp_injection_c;
    end

    function s_under = lower_source(this, glob, model)

      X = glob(:,1);
      Y = glob(:,2);

      D3 = ((X-0.8).^2 + (Y-0.5).^2) <= 0.01;
      s_under = 30 * D3;
    end

    function s_above = upper_source(this, glob, model)

      X = glob(:,1);
      Y = glob(:,2);

      D1 = ((X-0.5).^2 + (Y-0.8).^2) <= 0.01;
      D2 = ((X-0.2).^2 + (Y-0.2).^2) <= 0.01;
      s_above = 10 * D1 + 20 * D2;
    end

    function descr = default_descr(this, descr)
      descr.tp_pb          = -0.5;
      descr.tp_kw_factor   = 1;
      descr.tp_ko_factor   = 1;
      descr.tp_injection_c = 0.8;
    end
  end
end