Drag force variation at different velocities & Drag Co-efficient using python program

AIM: To find the variation of Drag foce with respect to different velocities & Drag Co-efficients using python programming.

Problem Statement:

        1.Write a program to find the variation of Drag force vs Velocity.

        2.Write a program to find the variation of Drag force vs Drag Co-efficients.

Theory:

        A Drag force is a resistance force which is caused by the motion of a body through a medium such as fluid-(Air,Water).The Drag force tends to act in the direction opposite to the motion of the body.It acts an opposing force to the object in motion.Due to this phenomenon,There is a relative motion existing between the Drag force and the velocity of the body.

The equation of drag force is stated as:

`Drag=frac{1}{2}cdotrhocdotAcdotv^2cdotCd`

Solution Procedure:

 Program to find the variation of Drag force acting on a bi-cyclist at different Velocities

""" Program to find the variation of Drag force acting on a bi-cyclist at different Velocities"""
#Importing required Libraries
import math
import matplotlib.pyplot as plt 
#Inputs

#Drag Co-efficient Unit(No Unit)
c_d=0.8

#Frontal Area Unit(m^2)
Area=0.1

#Density of Resisting Unit(kg/m^3)
rho=1.3

#List of Velocity Unit(m/s)
Velocity=[5,6,7,8,9,10,11,12]

#Empty List to insert Drag force values  Unit(kgm/s^2)
Drag_Force=[]

#Using for loop to append values of Velocity into Drag_Forces
for v in Velocity:
     Drag=(0.5*rho*Area*v*v*c_d)
     Drag_Force.append(Drag) #Calculating Drag Force
print(Drag_Force)	

#Plotting the result
plt.plot(Velocity,Drag_Force)
plt.xlabel('Velocity')
plt.ylabel('Drag Force')
plt.savefig('Drag Force vs Velocity')
plt.show()

Output Graph:

 Program to find the variation of Drag force acting on a bi-cyclist at different Drag Co-efficients:

""" Program to find the variation of Drag force acting on a bi-cyclist at different Drag Co-efficients"""
#Importing required Libraries
import math
import matplotlib.pyplot as plt 
#Inputs

#Velocity Unit(m/s)
Velocity=3

#Frontal Area Unit(m^2)
Area=1

#Density of Resisting Unit(kg/m^3)
rho=1.3

#List of Drag Co-efficients Unit(No Unit)
c_d=[0.2,0.38,0.48,0.67,0.87]

#Empty List to insert Drag force values  Unit(kgm/s^2)
Drag_Force=[]

#Using for loop to append values of Velocity into Drag_Forces
for cd in c_d:
        Drag_Force.append(0.5*rho*Area*Velocity*Velocity*cd) #Calculating Drag Force

print(Drag_Force)	

#Plotting the result
plt.plot(Drag_Force,c_d)
plt.xlabel('Drag Force')
plt.ylabel('Drag Co-efficients')
plt.savefig('Drag Force vs Drag Co-efficients')
plt.show()

Output Graph:

Results:

             1.From the first graph obtained through the programs,We can say that we obtain a parabolic  curve since Drag force is proportional to the square of the velocity of the cyclist.

             2.From the second graph obtained through the program,We can say that we obtain a straight  and linear curve since the Drag force is directly proportional the Drag Co-efficient.

Conclusion: 

             From this report we can conclude that the  Drag force is proportional to the square of the velocity of the cyclist.From the second graph obtained we can conclude that the

             Drag Co-efficient is directly proportional to the Drag force.