1 function descr = twophase_descr(params)
3 if nargin < 1 || ~isfield(params, 'size')
7 if ~isfield(params, 'mode')
8 params.mode = 'normal';
11 descr = TwoPhaseFlow.descr_default;
12 descr.rb_problem_type = 'TwoPhaseFlow';
13 descr.gridtype = 'rectgrid';
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);
57 descr.linear_velocity =
false;
59 descr = descr.two_phase.default_descr(descr);
64 % Neumann boundary values???
function model = model_default(model, T, nt)
model = model_default(model)
function U0 = init_values_waveproduct(glob, params)
product of two sinus waves in both coordinate directions