rbmatlab  1.16.09
 All Classes Namespaces Files Functions Variables Modules Pages
twophase_descr.m
1 function descr = twophase_descr(params)
2 
3 if nargin < 1 || ~isfield(params, 'size')
4  params.size = 1;
5 end
6 
7 if ~isfield(params, 'mode')
8  params.mode = 'normal';
9 end
10 
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;
16 descr.xrange = [0,1];
17 descr.yrange = [0,1];
18 descr.bnd_rect_corner1 = [-1,-1];
19 descr.bnd_rect_corner2 = [Inf, Inf];
20 descr.bnd_rect_index = -2;
21 
22 descr.mu_names = {'c_init'};
23 descr.mu_ranges = {[0.4, 0.6]};
24 
25 descr.T = 0.5;
26 descr.nt = 200 * params.size; % guess and adjust later
27 %descr.T = 0.1;
28 %descr.nt = 10;
29 
30 descr.verbose = 1;
31 descr.debug = 0;
32 
33 descr.name = 'single_test';
34 
35 % descr.init_values_ptr = @init_values_waveproduct;
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;
43 
44 descr.init_values_ptr = @init_values_homogeneous;
45 descr.c_init = 0.50;
46 
47 descr.mean_ptr = @(X)harmmean(max(X,0), 2);
48 switch params.mode
49  case 'normal'
50  descr.two_phase = TwoPhaseData.Michel;
51  case 'test'
52  descr.two_phase = TwoPhaseData.MichelSimple;
53  case 'linear'
54  descr.two_phase = TwoPhaseData.MichelPrsLin;
55 end
56 
57 descr.linear_velocity = false;
58 
59 descr = descr.two_phase.default_descr(descr);
60 
61 descr = model_default(descr);
62 % descr.two_phase_class = @TwoPhaseData.BrooksCorey;
63 
64 % Neumann boundary values???
65 
function model = model_default(model, T, nt)
model = model_default(model)
Definition: model_default.m:17
function U0 = init_values_waveproduct(glob, params)
product of two sinus waves in both coordinate directions