Flow over Bicycle

%Incorporate a function to plot Drag force vs Velocity:

clear all
close all
clc

%inputs

c_d=0.8; % drag coefficient
v=linspace(5,20,30); % velocity in m/sec
A=0.1; % frontal area in m^2
rho=1.2; % density in kg/m^3
% drag force calculation
drag_force=(c_d*rho*A*v.^2*0.5);
% plotting
plot(v,drag_force,'color','g')
xlabel('Velocity')
ylabel('Drag force') grid on

%drag coefficient for various shapes and try to plot drag force as a function of drag coeffiecient:

clear all
close all
clc

%inputs

v=100; % velocity in m/sec
A=0.1; % frontal area in m^2
rho=1.2; % density in kg/m^3

% flat plate has Cd = 1.28, %a wedge shaped prism with the wedge facing downstream has Cd = 1.14, %a sphere has a Cd that varies from .07 to .5, a bullet Cd = .295, %and a typical airfoil Cd= .045.

c_d=[1.28,1.14,.07,.295] % for various shapes drag_force=(rho*A*0.5*c_d.*v^2)
plot(c_d,drag_force,'color','b')
xlabel('c_d')
ylabel('drag force') grid on