1 function descr = twophase_descr(params)
3 if nargin < 1 || ~isfield(params, 'size')
7 if ~isfield(params, 'mode')
8 params.mode = 'normal';
12 descr.rb_problem_type = 'TwoPhaseFlow';
14 descr.xnumintervals = 100 * params.size;
15 descr.ynumintervals = 100 * params.size;
18 descr.bnd_rect_corner1 = [-1,-1];
19 descr.bnd_rect_corner2 = [Inf, Inf];
20 descr.bnd_rect_index = -2;
22 descr.mu_names = {
'c_init'};
23 descr.mu_ranges = {[0.4, 0.6]};
26 descr.nt = 200 * params.size; % guess and adjust later
33 descr.name =
'single_test';
36 % % parameters
for data functions
37 % descr.c_init_min = 0.0;
38 % descr.c_init_max = 1.0;
39 % descr.c_init_phase_x = 0.0;
40 % descr.c_init_phase_y = 0.0;
41 % descr.c_init_freq_x = 2*pi;
42 % descr.c_init_freq_y = 2*pi;
44 descr.init_values_ptr = @init_values_homogeneous;
47 descr.mean_ptr = @(X)harmmean(max(X,0), 2);
50 descr.two_phase = TwoPhaseData.Michel;
52 descr.two_phase = TwoPhaseData.MichelSimple;
54 descr.two_phase = TwoPhaseData.MichelPrsLin;
57 descr.linear_velocity =
false;
59 descr = descr.two_phase.default_descr(descr);
62 % descr.two_phase_class = @TwoPhaseData.BrooksCorey;
64 % Neumann boundary values???