burgers_fem.m

Go to the documentation of this file.

% hauptskript fuer burgers-fem-rb

params = [];

% Data settings (k is diffusivity):
params.name_dirichlet_values = 'gauss_convcomb';
params.udir_height = 3;
params.k = 1;   %linear test 1
params.cx = 0;
params.cy = 0;
params.udir_amplitude = 1;
params.udir_xscale = 1;
params.name_diffusivity = 'homogeneous'; 
params.name_velocity = 'linear'; 
params.T = 0.2;
params.nt = 2;
params.newton_eps = 1e-5;
params.max_newton_steps = 10;
params.fixed_newton_steps = 1; % turn on or off

% RB settings
params.rb_problem_type = 'burgers_fem';
params.mu_names = {'k','cx','cy','udir_height','udir_amplitude','udir_xscale'};
params.mu_ranges = {[1,2],[-5,5],[-5,5],[1,6],[1,2],[1,1000]};
params.offline_path = ...
    '/usr/people/naddel/fuelcell/probs/Burger/BURGER_OFFLINE_DATA';

params.verbose = 5;

% select here, whether subset-selection of reduced basis is to be
% performed or not:
select_offline_subset = 0;
%select_offline_subset = 1;

% RB Simulation:
detailed_data = rb_detailed_prep(params);
offline_data = rb_offline_prep(detailed_data, params);

params.N = 4;

if select_offline_subset
  % perhaps play around with accuracy threshold  1e-10 in the
  % following line
  indices = regular_submatrix_indices(offline_data.M,1e-10);
  % for test:
  %indices = [1,2];
  offline_data = rb_burgers_fem_offline_subset(offline_data,indices);
  old_N = params.N; 
  params.N = length(indices);
end;
  
reduced_data = rb_online_prep(offline_data, params);
params = set_mu([0.1, 0,0, 3,1,1],params);
%linear test 0.1,0,0,1,1,1

simulation_data = rb_simulation(reduced_data,params);

if select_offline_subset
  % extension of simulation output by null-entries for not used
  % basisvectors
  params.N = old_N;
  a = zeros(params.N,size(simulation_data.a,2));
  a(indices,:) = simulation_data.a;
  simulation_data.a = a;
end;

disp(simulation_data.a);


cd /usr/people/naddel/fuelcell/probs/Burger/BURGER_OFFLINE_DATA/   

fid = open_writefile('matrixB');  
store_dune_matrix(fid,simulation_data.b);
fclose(fid);

fid = open_writefile('matrixA');  
store_dune_matrix(fid,simulation_data.a); 
fclose(fid);  

% TO BE ADJUSTED TO NEW SYNTAX
%| \docupdate