demo_rb_stokes_fem.m

Go to the documentation of this file.

function demo_rb_stokes_fem(mesh_number)
%function demo_rb_stokes_fem(mesh_number)
% This function demonstrates the reduced basis approximation of the
% stationary Navier-Stokes equations (2-d CFD benchmark problem).

if (nargin < 1)
    
    mesh_number = 1;
end


params = [];
params.mesh_number = mesh_number;

% create model
model = laminar_flow_model(params);
model.verbose = 2;

% generate model data
model_data = gen_model_data(model);

% generate reduced basis
fprintf('\n');
disp('-> generate radial basis interpolant for inf-sup constant approximation');
model_data.infsup_constant = stokes_infsup_rbf_interpolant(model, model_data);
fprintf('\n');
disp('-> approximate continuity constant (sobolev embedding)');
model_data.rho_sqr = stokes_sobolev_embedding_constant(model, model_data);
fprintf('\n');
disp('-> reduced basis greedy algorithm');
detailed_data = gen_detailed_data(model,model_data);
reduced_data = gen_reduced_data(model,detailed_data);

% plot solution
plot_params = [];
% specify functions to be plotted
%plot_params.modes = { ...
%    'pressure', ...
%    'velocity_xcomp', ...
%    'velocity_ycomp', ...
%    'velocity_abs', ...
%    'velocity_vec'};
plot_params.modes = {'pressure', 'velocity_abs'};
plot_params.subsampling_level = 0;
plot_params.vector_plot_distance = 0.02;
plot_params.no_lines=1;
plot_params.rbgui=1;
plot_params.hide_N_ruler=1;

%GUI demo
fprintf('\n');
disp('Please call demo_rb_stokes_fem(3) for a better resolution');

model.descr.is_stationary = 1;
model.enable_error_estimator = 0;
model.descr.gridtype ='triagrid';
model.N  = get_rb_size(model,detailed_data);

demo_rb_gui(model,detailed_data,reduced_data,plot_params);

end