data for Cd and terminal velocity are done

The Rocketry Forum

Help Support The Rocketry Forum:

This site may earn a commission from merchant affiliate links, including eBay, Amazon, and others.

Lt72884

Well-Known Member
TRF Supporter
Joined
Sep 11, 2022
Messages
506
Reaction score
174
Location
Utah
Ok, so i finally finished my personal project, BUT it will be used for my senior paper since its interesting results.
Here is the background; I took sensor data that was given to me from various members here on the forums. I then used excel in 2 different ways to find the Drag Coefficient (Cd) of the rocket.

A little about Cd for beginners, this number usually has been calculated from HUNDREDS and THOUSANDS of tests in a windtunnel. IE, they take the rocket, put it in the wind tunnel and find the FORCE the wind pushes on the rocket at low, medium and HIGH speeds. Using the Drag Force equation

1677711941386.png
You then massage the above equation a little bit to solve for Cd. Since you have the force of drag(Fd), area, velocity, and air density. you can calculate the Cd quickly with wind tunnel data.
But what if you didnt have a wind tunnel? This is where i come in, i kind of made a "virtual" wind tunnel in a way, but not really. I used excel to calculate theoretical (modeled) velocity values based off of total force the rocket is feeling and acceleration as in the next set of equations: Frr is for friction of a car, not a rocket. its not needed for rockets
{aid:null,type:lalign*,code:\\begin{lalign*}\n&{Fd\\,=\\,\\frac{\\left(Cd\\right)\\left(A\\right)\\left(v^{2}\\right)\\left(\\rho\\right)}{2}}\\\\\n&{}\\\\\n&{Frr\\,=\\,Crr\\left(m\\right)\\left(g\\right)}\\\\\n&{}\\\\\n&{Ft=\\,Fd+Frr}\\\\\n&{}\\\\\n&{F=ma\\,\\,}\\\\\n&{}\\\\\n&{Ft=ma\\,\\to\\,a=\\frac{Ft}{m}\\to\\,\\frac{Fd+Frr}{m}}\\\\\n&{}\\\\\n&{\\Delta v\\,=\\,a\\left(\\Delta t\\right)\\,\\to\\,a=\\frac{\\Delta v}{\\Delta t}}\t\n\\end{lalign*},id:1,font:{family:Arial,color:#000000,size:14},backgroundColor:#ffffff,ts:1674327143399,cs:IQ7hpZqro4vxgePE4AAYKQ==,size:{width:334,height:382}}


This gave me a list of MODELED aka calculated velocities. I then had excel compute the difference between my calculated velocities and the ACTUAL velocities then squared that value for my error. I gave excel an initial guess for Cd, because you ALWAYS have to have an initial guess to work with, even if its WAY off. Then using the excel add-in called solver, i let it adjust the Cd until my error was close to 0. Once it found this Cd value, i graphed the modeled velocities to the actual and they were right on top of each other, meaning, that i found the Cd of that rocket for that particular set of data.

To verify this even more, i calculated the Cd from a first order diff eq for EACH velocity, so about 389 different entries, using the following equations:
{aid:null,code:\\begin{lalign*}\n&{F=ma\\to m\\frac{dv}{dt}\\,where\\,dv/dt\\,is\\,aceleration}\\\\\n&{Ft=\\,-Fd-Fg}\\\\\n&{m\\frac{dv}{dt}\\,=\\,-Fd-Fg}\\\\\n&{Fd\\,=\\,Cd\\left(A\\right)\\left(v^{2}\\right)\\left(\\rho\\right)\\left(\\frac{1}{2}\\right)}\\\\\n&{let\\,B\\,=\\,\\left(Cd\\right)\\left(A\\right)\\left(\\rho\\right)\\left(\\frac{1}{2}\\right)}\\\\\n&{then\\,Fd\\,=\\,Bv^{2}}\\\\\n&{m\\frac{dv}{dt}=-Bv^{2}-mg}\\\\\n\\end{lalign*},backgroundColor:#ffffff,font:{size:11,family:Arial,color:#000000},id:2,type:lalign*,ts:1677696150897,cs:UfwOT7W80JXQLcyquAxGhQ==,size:{width:324,height:256}}


I then solved for B
1677712655931.png
Once i got a value back, such as 0.00647 for B, i then back substituted it back into the B=(Cd)(A)(rho)(0.5) and solved for Cd yielding 0.67 and open rocket says its 0.63.. so i am pretty close. Here is the graph of beta (y axis) and velocity (x axis) you can see how around 200m/s, this is where the Cd doesnt change much
1677712802151.png
1677712812659.png
SUPER grateful for the sensor data and ALL the help i received on this from the physics forums. I hope you enjoy. This was an excellent learning project for me.
 
Last edited:
Oh, you're posting on here too. It's funny - I saw your stuff on the physics forums and was referencing some of my old data from here, and noticed this. Glad you seem to have gotten it mostly figured out and your data looks reasonable (if a bit noisy, but that's kinda inevitable with speeds this low and a high sectional density rocket). Got any more details about the rocket itself or the flight?
 
Interesting. I will need to copy this out and study/play with the equations.
Which rocket data did you use for above graphs?
Did you do this with other rockets and motors?
In the data I gave you there are different motors on the same rocket. This should give a variation of maximum velocity.
 
Oh, you're posting on here too. It's funny - I saw your stuff on the physics forums and was referencing some of my old data from here, and noticed this. Glad you seem to have gotten it mostly figured out and your data looks reasonable (if a bit noisy, but that's kinda inevitable with speeds this low and a high sectional density rocket). Got any more details about the rocket itself or the flight?
LOL, i was wondering if it was the same person haha. I posted a pic of the rocket on the physics forums. its a snapshot from OR.
Thanks for all your help on it. its been a blast to see and do. Physics forums is a great place to be and i have loved my stay there. Same with this place.
 
Last edited:
Interesting. I will need to copy this out and study/play with the equations.
Which rocket data did you use for above graphs?
Did you do this with other rockets and motors?
In the data I gave you there are different motors on the same rocket. This should give a variation of maximum velocity.
so far it has only been with one of the rockets. Now that this is polished out, i can use the other rocket data and compare. Here are the files i created. This will let you see what i did. Enjoy.
 

Attachments

  • cd for rocket from raw data.xlsx
    85.1 KB · Views: 0
  • time and velocity from burnout to apogee.xlsx
    191.8 KB · Views: 0
This maybe a dumb question so sorry in advance, are you calculating Cd for the entire rockets itself and not just the parachute? Also do you think your calculations would work for calculating Cd of a parachute? Thank you!
 
Back
Top