varyingspeeds.m

close all;
clear;
clc;
format short;

source_speed = 0.05; % strength of source in hood, ~0.05 m3/s
show_graphs = false; %show intermediate graphs or not
threshold = 0.0002; % amount of gas (m3) at entrance before gas has escaped


% Compare low height with lower speeds
speeds_lo = 0:0.01:0.20; % range of vent speeds to check
volumes = zeros(1, length(speeds_lo));
height = 0.1;

for i = 1:length(speeds_lo)
    vent_speed = speeds_lo(i);
    volumes(i) = stationary(height, source_speed, vent_speed, show_graphs);
    disp([num2str(vent_speed), ',', num2str(volumes(i))]);
end
figure
plot(speeds_lo, volumes);

% Compare high height with higher speeds
speeds_hi = 0.20:0.01:0.50; % range of vent speeds to check
volumes = zeros(1, length(speeds_hi));
height = 0.65;

for i = 1:length(speeds_hi)
    vent_speed = speeds_hi(i);
    volumes(i) = stationary(height, source_speed, vent_speed, show_graphs);
    disp([num2str(vent_speed), ',', num2str(volumes(i))]);
end
figure
plot(speeds_hi, volumes);