simtech_fuelcell_gdl.m

Go to the documentation of this file.

% small script demonstrating RB-approach for the fuelcell-gdl with
% the simtech image initial data

% 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
detailedfname = 'simtech_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 = 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);
 otherwise
  error('step number unknown!');  
end;


return;


% TO BE ADJUSTED TO NEW SYNTAX
%| \docupdate