demo_rb_riemann_burgers.m

Go to the documentation of this file.

%demo_rb_riemann_burgers
%
% Script demonstrating the RB-model for the discontinuous 
% initial data for the nonlinear burgers equation. 
% Both moving shock front and rarefaction waves can be modelled.
% discretization is with explicit fv scheme and empirical
% interpolation of the space operator.
% The model demonstrates the space-dimension reduction: the 
% subgrid that is extracted in the offline phase is exactly the row
% of cells along the bottom line.
%
% To start: set u_left or u_right to some different value.
% Then observe either a moving shock or a rarefaction wave
% depending on the direction of the velocity.

% Bernard Haasdonk 10.9.2008

help demo_rb_riemann_burgers

detailedfname = 'riemann_burgers_detailed.mat';
load(detailedfname);

% change/add some settings, if datafile has been generated by earlier versions
%extend = 1;
extend = 0;
if extend
  detailed_data.grid = construct_grid(model);
  %  detailed_data.grid = extend_class_from_struct('triagrid',detailed_data.grid);
  model.newton_solver = 0;
  model.implicit_nonlinear = 0;
  model.local_stencil_size = 1;
  model.get_rb_size = @get_rb_size_default;
  model.enable_error_estimator = 0;
  model.get_mu = @get_mu_default;
  %  plot_params.axis_equal = 1;
  %  plot_params.axis_tight = 1;
%  save(fullfile(rbmatlabhome,'demos','datafiles',detailedfname),'detailed_data','plot_params','model');
end;

%[model,detailed_data,plot_params] = renew_model(model,detailed_data);

%model = riemann_burgers_model;

%model.N = model.get_rb_size(model,detailed_data);
%size(detailed_data.RB,2);
%model.M = size(detailed_data.QM{1},2);

rb_plot_interpolation_points(detailed_data,model);

%model.gen_reduced_data = @nonlin_evol_gen_reduced_data;
%model.rb_operators = @nonlin_evol_rb_operators;
%rb_plot_interpolation_points(detailed_data,model);
%model.set_time = @set_time_default;
%model.set_mu = @set_mu_default;
%model.reduced_data_subset = @nonlin_evol_reduced_data_subset;
%%model.rb_simulation = @nonlin_evol_rb_simulation;

% temporary: transform to "model"
% [model, detailed_data, plot_params] = renew_model(params, detailed_data);

demo_rb_gui(model,detailed_data); %,[],plot_params,'riemann burgers');