Python Project 1 : Flow over Bicycle

Procedure to calculate Drag Force in Python :-

For the calculation of Drag force all Input variables are predefined. So for that first we will assign the values of the Input variable . After the assignment of input variable values we will define the drag force according to the formula. This will give the value of drag force. Now we have to plot the graph between velocity and drag force, for this plotting we will define only single value to drag coefficient and multiple values to velocity variable. For this multiple velocity, Pytho will compute multiple values for Drag force. Then by using plotting commands we will visualize the relation between velocity and drag force.

For create a plot between drag coefficient and drag force we will define a single value to velocity and multiple values to drag force. This multiple values of drag coefficient gives multiple values of Drag force. For the multiple values of Drag coefficient we will visualize the relation between drag coefficient and drag force by using plotting command.

Python program :-

Google Colab is used for this programming

"""Project Title:- Flow over bicycle"""
# Importing matplot function for plotting
import matplotlib.pyplot as plt

# Input Data
# Density of Air (rho) in kg/m^3
rho = 1.225

# Frontal Area(A) in m^2
A= 0.1

# For plot between velocity and Drag force(Drag coefficient constant and
# velocity is varying.)

# Drag coefficient(c_d contain drag coefficient 0.8 for bicycle)
c_d = 0.8 

#Velocity range (velocity taken from 1 to 35 m/s)
velocities = [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35]

# Calculation of Drag force of varying velocity
# Using formula of Drag force
drag_force_velocity1 = []

for velocity in velocities:
     drag_force_velocity1.append((0.5*rho*A*c_d*velocity*velocity))

# plot of velocity vs drag force
plt.plot(velocities,drag_force_velocity1)
plt.xlabel('Velocity in m/s')
plt.ylabel('Drag force in Newton')
plt.title('Velocity v/s Drag force')
plt.show()

# Input for plot Dag coefficient v/s drag force where velocity is constant
# and drag coefficient has multiple values
# Drag coefficient range
c_d = [0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7]

# Velocity in m/s
v = 25 

# Drag force for varying drag coefficient
drag_force_coefficientrange=[]

for c_d_range in c_d:
    drag_force_coefficientrange.append(0.5*rho*A*v*v*c_d_range)

# plot of Drag coefficient vs drag force.
plt.plot(c_d,drag_force_coefficientrange)
plt.xlabel('Drag coefficient')
plt.ylabel('Drag force in Newton')
plt.title('Drag coefficient v/s Drag force')
plt.show()
#end of the program​