MoRePaS09_fuelcell_gdl.m

Go to the documentation of this file.

% small script demonstrating RB-approach for the fuelcell-gdl with
% the MoRePaS09 initial data, i.e. Logo for the Workshop
% Model Reduction of Parametrized Systems

% Bernard Haasdonk 24.5.2008

%step = 1; % single simulation
%step = 2; % dummy reduced basis & comparison of reduced to detailed simulation
%step = 3; % real reduced basis generation and store as detailed_data
%step = 4; % time comparison
%step = 5; % demo_rb_gui with modified fuel cell
step = 6; % export some pictures
%detailedfname = 'simtech_detailed.mat';
detailedfname = 'MoRePaS09_detailed.mat';

load demo_lin_evol_params;
params.RB_generation_mode: 'uniform_fixed'
params.RB_numintervals = [10 4 4];
params.mu_names ={'c_init'  'beta'  'k'};
params.mu_ranges ={[0 1]  [0 1]  [0 5.0000e-008]}

% display params
%params.name_convective_num_flux = 'enquist_osher';

params.name_init_values = 'grey_image';
%params.c_init_filename = 'simtech_binary.bmp';
params.c_init_filename = 'MoRePaS09_binary.bmp';
%params.c_init = 0.0;
params.c_init = 1.0;
params.c_init_xmin = params.xrange(1);
params.c_init_xmax = params.xrange(2);
params.c_init_ymin = params.yrange(1);
params.c_init_ymax = params.yrange(2);
params.c_init_xdivisions = (0.1:0.1:1.0) *...
    (params.xrange(2)-params.xrange(1));

params.k = 0;

params.name_convective_num_flux = 'engquist-osher';

grid = construct_grid(params);
%params.lxf_lambda = 1.0194e+003;
%fv_search_max_lxf_lambda([],params);
params.nt = params.nt;

% output settings
params.verbose = 10;
params.axis_equal = 1;
params.clim = [0,1];

detailed_data.grid = grid;

switch step
 case 1
  disp('performing single detailed simulation')
  U = detailed_simulation(detailed_data.grid, params);
  plot_element_data_sequence(detailed_data.grid,U,params);
 case 2 % construct dummy reduced basis by single trajectory and simulate
  U =  detailed_simulation(detailed_data.grid,params);
  A = feval(params.inner_product_matrix_algorithm,detailed_data.grid,params);
  UON = orthonormalize(U,A);
  detailed_data.RB = UON;
  disp('reduced simulation:')
  offline_data = rb_offline_prep(detailed_data,params);
  params.N = size(detailed_data.RB,2);
  reduced_data = rb_online_prep(offline_data, params);
  simulation_data = rb_simulation(reduced_data,params);
  Uappr = rb_reconstruction(detailed_data,simulation_data);

  params.title = 'reduced simulation result';   
  plot_element_data_sequence(detailed_data.grid,Uappr,params);
  params.title = 'detailed simulation result';  
  plot_element_data_sequence(detailed_data.grid,U,params);
 case 3 % generate good reduced basis
  disp('constructing reduced basis')
  detailed_data = rb_basis_generation(detailed_data, ...
                                      params);
  save(fullfile(rbmatlabresult,detailedfname),...
       'detailed_data','params');
  plot(detailed_data.RB_info.max_err_sequence);
  set(gca,'Yscale','log');
  title('RB-generation error convergence');
  
 case 4
  load(fullfile(rbmatlabresult,detailedfname));
  disp('reduced simulation:')
  offline_data = rb_offline_prep(detailed_data,params);
  params.N = size(detailed_data.RB,2);
  params.k = 0.0001;
  reduced_data = rb_online_prep(offline_data, params);
  tic;
  simulation_data = rb_simulation(reduced_data,params); 
  t = toc;
  disp(['time for online phase: t = ',num2str(t)]);
  
  disp('full simulation:')
  tic;
  U = detailed_simulation(detailed_data.grid, params);
  t = toc;
  disp(['time for detailed simulation: t = ',num2str(t)]);

 case 5
  load(fullfile(rbmatlabresult,detailedfname));
  demo_rb_gui(detailed_data,[],params);
 case 6
  %save image series
  disp('performing single detailed simulation')
  U = detailed_simulation(detailed_data.grid, params);
  % slice 26 is good  
  h = hot(256);
  p = pink(256);
    slices = [1, 26, 51, 76, 101, 126,151,176,201];
    colormaps = {h(end:-1:1,:)};
%  slices = [26];
%  colormaps = {h(end:-1:1,:),h,p,p(end:-1:1,:),jet(256)};
  for s=1:length(slices)
    for cm = 1:length(colormaps)
      params.show_colorbar=0;
      plot_element_data(grid,U(:,slices(s)),params);
      colormap(colormaps{cm});
      %      a = reshape(U(:,slices(s)),200,40)';
      %      a = uint8(a(end:-1:1,:)*255);
      %      imshow(a,colormaps{cm});
      set(gca,'XTickLabel',{})
      set(gca,'YTickLabel',{})
      saveas(gca,['MoRePaS09_logo_t',num2str(slices(s)),'_col',num2str(cm)],...
              'epsc2')
      close(gcf);
%      keyboard;
    end;
  end;
 otherwise
  error('step number unknown!');  
end;


return;


% TO BE ADJUSTED TO NEW SYNTAX
%| \docupdate