models/convdiff/convdiff_descr.m

function descr=convdiff_descr(dummy)

  descr.name = 'convdiff_descr';
  descr.rb_problem_type = 'LinEvol';
  %% data that is definitely used outside of the detailed simulations
  descr.mu_names           = {'c_init','beta','k'};
  descr.mu_ranges          = {[0 1],[0 1],[0 5e-8]};

  %% data that might be of interest outside of the detailed simulations
  descr.T             = 0.5;
  descr.nt            = 200;
  %descr.dt;
  %descr.mu;

  %descr.N;
  %descr.Nmax;
  %descr.RB_numintervals;

  descr.xrange        = [0,1000e-6]; % grid x-coordinate range
  descr.yrange        = [0,200e-6];  % grid y-coordinate range
  descr.xnumintervals = 200;  % number of grid cellls in x-direction
  descr.ynumintervals = 40;   % number of grid cellls in y-direction

  descr.operators_ptr              = @fv_operators_implicit_explicit;
  descr.operators_diff_implicit    = @fv_operators_diff_implicit_gradient;
  descr.operators_conv_explicit    = @fv_operators_conv_explicit_lax_friedrichs;
  descr.operators_neumann_explicit = @fv_operators_neumann_explicit;

  %% specify initial data function
  %    name_init_values = 'homogeneous';
  %    c_init = 0.0;
  descr.init_values_ptr = @init_values_wave;
  descr.freq_x = 5*2*pi/1000e-6;
  descr.freq_y = 0;
  descr.c_init = 1.0;
  %    name_init_values = 'homogeneous';
  %    c_init = 0.3;

  %    name_init_values = 'leftright';
  %    c_init_left = 0.3;
  %    c_init_right = 0.7;
  %    init_middle = (params.xrange(2)+params.xrange(1))/2;

  %% convective flux specification
  %    conv_flux_ptr = @conv_flux_gdl2;
  descr.conv_flux_ptr = @conv_flux_velocity_matrixfile;
  descr.lxf_lambda              = 1.0194e+003;
  descr.divclean_mode           = false;
  descr.flux_quad_degree        = 1;
  descr.use_velocity_matrixfile = 1;

  %descr.velocity_matrixfile;
  descr.filecache_velocity_matrixfile_extract = 0;

  %% data that is relevant for detailed simulations only

  %% specify grid type
  % generate cartesian grid
  descr.gridtype = 'rectgrid';

  %% specify dirichlet boundary values
  descr.dirichlet_values_ptr = @dirichlet_values_weighted_boxes;
  descr.c_dir                 = 1.0; % maximum value
  descr.beta                  = 0.0;  % weighting parameter: 1=box1, 0 = box2

  %% specify neuman boundary values
  descr.neumann_values_ptr = @neumann_values_rightflow;
  %    name_neuman_values = 'zero';
  %    name_neuman_values = 'homogeneous';
  %    c_neu = 0;

  %% specify further data functions
  descr.flux_linear      = 1;
  descr.diffusivity_ptr  = @diffusivity_homogeneous;
  %    k                       = 0.000000;   % diffusion coefficient
  %    k                       = 0.000001;   % diffusion coefficient
  %    k                       = 1e-15;      % diffusion coefficient
  %    k                       = 1e-10;      % diffusion coefficient
  descr.k                     = 1e-8;       % diffusion coefficient

  % precomputed, divcleaned velocity field
  descr.name_flux               = 'gdl2';

  descr.lambda = 2.0e-7;   % v = - lambda * grad p
  %  descr.inner_product_matrix_algorithm = @fv_inner_product_matrix;
  descr.energy_norm_gamma    = 0.5;
  descr.alpha_over_k         = 1.9915e7;
  descr.coercivity_bound_ptr = @fv_coercivity_bound;


  % fields to ignore, such that filecaching works in basis generation
  descr.filecache_ignore_fields_in_descr = {'N','Nmax'};

  %bnd_rect_corner1;
  %bnd_rect_corner2;
  %bnd_rect_index   = [];

  descr.bnd_rect_corner1 = [ descr.xrange(1) ,         descr.yrange(2) ; ...
                             descr.xrange(1)+ 250e-6 , descr.yrange(2); ...
                             descr.xrange(1) ,         descr.yrange(1); ...
                             descr.xrange(1) ,         descr.yrange(1);
                             descr.xrange(2) ,         descr.yrange(1)
                           ];
  descr.bnd_rect_corner1 = descr.bnd_rect_corner1' - eps;
  descr.bnd_rect_corner2 = [ descr.xrange(2),        descr.yrange(2); ...
                             descr.xrange(2)-250e-6, descr.yrange(2); ...
                             descr.xrange(2),        descr.yrange(1);
                             descr.xrange(1),        descr.yrange(2);
                             descr.xrange(2),        descr.yrange(2)
                           ];
  descr.bnd_rect_corner2 = descr.bnd_rect_corner2' + eps;
  descr.bnd_rect_index   = [-1, -2, -2 , -1, -2 ...
                           ];

  descr.use_velocity_matrixfile = 1;

  %in case of later matrix generation: downscale original velocities
  %    descr.divclean_downscale = 1;
  %    descr.divclean_downscale_quota = 0.04;

  % set matrix-file name for optional generation or reading of file
  descr.velocity_matrixfile = [ 'vel_', descr.name_flux,'_',...
                                num2str( descr.xnumintervals ),'x',...
                                num2str( descr.ynumintervals ),...
                                '_l',num2str(descr.lambda),'.mat'];

  descr.dir_box_xrange = {[ descr.xrange(1)-eps, ...
                            0.5*( descr.xrange(1) + descr.xrange(2) )], ...
                          [ 0.5*( descr.xrange(1) + descr.xrange(2) )-eps, ...
                            descr.xrange(2)]...
                         };
  descr.dir_box_yrange = {[ descr.yrange(2)-eps,     descr.yrange(2)+eps], ...
                          [ descr.yrange(2)-eps,     descr.yrange(2)+eps] ...
                         };

  descr = model_default(descr);

  descr.verbose = 5;
  descr.debug   = 0;


  %% basis generation options
  descr.RB_train_num_intervals = 4;

  % In order to prevent name clashes the basis generation options can also be
  % prefixed with BG_
  descr.RB_val_size         = 10;

  descr.RB_stop_epsilon     = 1e-8;

  descr.RB_refinement_mode  = 'adaptive';

end