models/twophase_descr.m

function descr = twophase_descr(params)

if nargin < 1 || ~isfield(params, 'size')
  params.size = 1;
end

if ~isfield(params, 'mode')
  params.mode = 'normal';
end

descr = TwoPhaseFlow.descr_default;
descr.rb_problem_type = 'TwoPhaseFlow';
descr.gridtype = 'rectgrid';
descr.xnumintervals = 100 * params.size;
descr.ynumintervals = 100 * params.size;
descr.xrange = [0,1];
descr.yrange = [0,1];
descr.bnd_rect_corner1 = [-1,-1];
descr.bnd_rect_corner2 = [Inf, Inf];
descr.bnd_rect_index   = -2;

descr.mu_names  = {'c_init'};
descr.mu_ranges = {[0.4, 0.6]};

descr.T  = 0.5;
descr.nt = 200 * params.size; % guess and adjust later
%descr.T  = 0.1;
%descr.nt = 10;

descr.verbose = 1;
descr.debug = 0;

descr.name = 'single_test';

% descr.init_values_ptr = @init_values_waveproduct;
% % parameters for data functions
% descr.c_init_min = 0.0;
% descr.c_init_max = 1.0;
% descr.c_init_phase_x = 0.0;
% descr.c_init_phase_y = 0.0;
% descr.c_init_freq_x = 2*pi;
% descr.c_init_freq_y = 2*pi;

descr.init_values_ptr = @init_values_homogeneous;
descr.c_init = 0.50;

descr.mean_ptr  = @(X)harmmean(max(X,0), 2);
switch params.mode
  case 'normal'
    descr.two_phase = TwoPhaseData.Michel;
  case 'test'
    descr.two_phase = TwoPhaseData.MichelSimple;
  case 'linear'
    descr.two_phase = TwoPhaseData.MichelPrsLin;
end

descr.linear_velocity = false;

descr = descr.two_phase.default_descr(descr);

descr = model_default(descr);
% descr.two_phase_class = @TwoPhaseData.BrooksCorey;

% Neumann boundary values???