ulm_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

disp('THIS SCRIPT STILL USES VERY OLD MECHANISM')
disp('MUST BE UPDATED!');

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

%load demo_lin_evol_params;
model = convdiff_model;
model.detailedfname = 'ulm_detailed.mat';
model.RB_generation_mode = 'greedy_uniform_fixed';
model.RB_numintervals = [10 4 4];
model.mu_names ={'c_init'  'beta'  'k'};
model.mu_ranges ={[0 1]  [0 1]  [0 5.0000e-008]};
model.RB_stop_Nmax = 100;

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

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

model.k = 0;

model.operators_conv_explicit = @fv_operators_conv_explicit_engquist_osher;

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

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

model = model;
model_data = gen_model_data(model);

%detailed_data.grid = grid;

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

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

 case 5
  load(fullfile(rbmatlabresult,model.detailedfname));
  plot_params = [];
  plot_params.axis_tight = 1;
  demo_rb_gui(model,detailed_data,[],plot_params);
 otherwise
  error('step number unknown!');  
end;


return;


% TO BE ADJUSTED TO NEW SYNTAX
%| \docupdate