Ethanol Rocket: Difference between revisions
(23 intermediate revisions by the same user not shown) | |||
Line 5: | Line 5: | ||
=== Ethanol 3d model === | === Ethanol 3d model === | ||
<gallery> | |||
Ethanol.png|The Ethanol | |||
Ethanol_0244.png|C2H6O | |||
Ethanol_0270.png|The Ethanol | |||
</gallery> | |||
{| class="wikitable" | {| class="wikitable" | ||
Line 36: | Line 41: | ||
|+ Van der Waals radius | |+ Van der Waals radius | ||
|- | |- | ||
! Atom !! Radius | ! Atom !! Radius !! Relative radius | ||
|- | |- | ||
| C || 170 pm || 1.42 | | C || 170 pm || 1.42 | ||
Line 58: | Line 63: | ||
The volume-% of oxygen in air is 20.9%. The volume of the bottle is <math>V = 500 </math> ml which gives the amount of oxygen to be <math>0.209 V = 0.209 \times 500 = 104.5</math> ml <math>=0.1045\times1.314 = 0.137 </math> g of oxygen, which gives the amount of ethanol <math>0.137 / 2.1 = 0.0659 </math> g <math>=0.0659/0.789= 0.0834</math> ml. '''OR''' directly using air--ethanol ratio we have <math>\rho V/9 = 1.2041 \times 0.5 /9 = 0.066</math> gram. That amount equals to <math>V/\rho_\textrm{Ethanol} = 0.066g/(789g/l) = 0.084 </math> ml <math>=84 </math>mm<sup>3</sup>. Almost the same result using the methods. See the attached spreadsheet for detailed calculations. | The volume-% of oxygen in air is 20.9%. The volume of the bottle is <math>V = 500 </math> ml which gives the amount of oxygen to be <math>0.209 V = 0.209 \times 500 = 104.5</math> ml <math>=0.1045\times1.314 = 0.137 </math> g of oxygen, which gives the amount of ethanol <math>0.137 / 2.1 = 0.0659 </math> g <math>=0.0659/0.789= 0.0834</math> ml. '''OR''' directly using air--ethanol ratio we have <math>\rho V/9 = 1.2041 \times 0.5 /9 = 0.066</math> gram. That amount equals to <math>V/\rho_\textrm{Ethanol} = 0.066g/(789g/l) = 0.084 </math> ml <math>=84 </math>mm<sup>3</sup>. Almost the same result using the methods. See the attached spreadsheet for detailed calculations. | ||
The energy released by burning ethanol is 17.9 kJ/ml. Thus, <math>0.084</math> ml of ethanol releases <math>1,50</math> kJ of energy. This energy is converted into heat, sound and projectile motion (plus others). | |||
{| class="wikitable" | {| class="wikitable" | ||
|+ | |+ Densities | ||
|- | |- | ||
! !! Density <math>\rho</math>!! At | ! !! Density <math>\rho</math>!! At | ||
Line 76: | Line 84: | ||
The inner diameter of the rocket bottle is 25 mm. The height of the ethanol in the cap need to be <math> V = \pi r^2 h \iff h = \frac{V}{\pi r^2} = \frac{84 mm^3 }{ \pi 12.5^2 mm^2} = 0.17 </math> mm. | The inner diameter of the rocket bottle is 25 mm. The height of the ethanol in the cap need to be <math> V = \pi r^2 h \iff h = \frac{V}{\pi r^2} = \frac{84 mm^3 }{ \pi 12.5^2 mm^2} = 0.17 </math> mm. | ||
== Exothermic Reaction and Energy release == | |||
=== Specific Heat === | |||
The energy is transferred into pressure, sound, etc. The isochoric specific heat <math>C_v</math> of air is <math>C_v = 0.7171</math> kJ/(kgK) at 18 centigrade. At 180 degrees Celsius <math>C_v = 0.7352</math> kJ/(kgK). Thus, the energy released heats | |||
<math> | |||
\begin{align} | |||
\Delta E &= C_v \Delta T m \\ | |||
\Delta T &= \frac{\Delta E}{C_v m} \\ | |||
&=\frac{1.5 \textrm{kJ}}{0.7 \textrm{kJ/(kgK)} \times 0.5 \times 10^{-3}} \textrm{kg} \\ | |||
&= 2100 K | |||
\end{align} | |||
</math> | |||
=== Ideal gas law === | |||
The simplest idea is to use | |||
ideal gas law for isochoric process <math>p/T = </math> constant gives <math>p_2 = \frac{T_2}{T_1}p_1</math> which gives force <math>F = pA</math> where <math>A</math> is the diameter of the nozzle. | |||
Thus, the propulsive force is | |||
<math> | |||
\begin{align} | |||
F &= pA \\ | |||
&= \frac{T_2}{T_1}p_1 A \\ | |||
&= \frac{T_1 + \Delta T}{T_1}p_1 A \\ | |||
&= (1 + \frac{\Delta T}{T_1}) p_1 A \\ | |||
&= (1 + \frac{\Delta E/ (C_v m)}{T_1}) p_1 A \\ | |||
\end{align} | |||
</math> | |||
=== Detonation velocity === | |||
Burn rate, detonation velocity. | |||
https://en.wikipedia.org/wiki/Table_of_explosive_detonation_velocities | |||
http://www.explosionsolutions.co.uk/110411016.pdf | |||
https://link.springer.com/article/10.1007/s00193-015-0554-7 <math>v=1500</math> m/s??? | |||
Horizontal(?) accelaration due to rapidly expanding air. | |||
Benjamin Robins: | |||
<math> | |||
F(x) = \frac{RPAc}{x} | |||
</math> | |||
where <math>x</math> distance in barrel, <math>R</math> is the initial ratio of hot gas pressure to atmospheric pressure, <math>P</math> is the atmospheric pressure, <math>A</math> is the cross-sectional area of the ball or bore, <math>c</math> is the length of the barrel occupied by the powder charge before ignition. | |||
<math> | |||
E_k = \int_0^L F(x) dx \iff v_0^2 = \frac{2RP}{m} \frac{\pi d^2 c}{4}\ln(L/c) | |||
</math> | |||
where <math>d</math> is the barrel diameter (the bore) | |||
The powder change <math>p</math> is given by | |||
<math>p= \frac{\pi d^2 c}{4}\eta </math> | |||
where <math>\eta</math> is the density of gunpowder. | |||
-- | |||
https://www.arc.id.au/RobinsOnBallistics.html | |||
The pressure falls as <math>1/x</math>. | |||
-- | |||
Bernoulli? | |||
-- | |||
<math>\sum F = 0 </math> | |||
<math>F = ma = pA \iff F = m \frac{d^2x}{dt^2} </math> | |||
<math> F = m \frac{dv}{dt} = mv \frac{dv}{dx} = pA \iff mv dv = pA dx </math> | |||
<math>\int mvdv = \int_0^L pA dx </math> | |||
where | |||
<math> | |||
\frac1L \int_0^L P dx = average pressure = \vec p | |||
</math> | |||
If <math>A</math> is constant | |||
<math> | |||
\frac12 mv^2 = A \int_0^L p dx = A \vec p L | |||
</math> | |||
then | |||
<math> | |||
v = \sqrt{\frac{2\vec p A L}{m c_f}} | |||
</math> | |||
and | |||
<math> | |||
\vec p = \frac{c_f mv^2}{2AL} | |||
</math> | |||
Projectile friction, rotational energy, heat transfer: correction factor <math>c_f</math>. | |||
The average pressure is 25% of peak pressure | |||
-- | |||
http://closefocusresearch.com/calculating-barrel-pressure-and-projectile-velocity-gun-systems | |||
https://www.arc.id.au/CannonBallistics.html | |||
== Maximum Flying Distance == | |||
By the conservation of energy, all explosive energy is transferred into the kinetic energy. | |||
<math>E_\textrm{explosion}=E_\text{kinetic}</math> | |||
The distance covered by a projectile with initial velocity $v_0$ is | |||
<math>s = \frac{v_0^2}{g}\sin(2\alpha) = \frac{}{}</math> | |||
The drag coefficient <math>c_d</math> needs to be found. For the circular disc (a coin) the drag coefficient is almost constant for all velocities (Reynold numbers). The coefficient of drag for a cylinder in this orientation is about 0.81 so long as the length to diameter ratio is greater than 2, see http://www.aerospaceweb.org/question/aerodynamics/q0231.shtml. The cone in either end gives some complications. | |||
The drag equation for the drag force is <math>F_D = \frac12 \rho v^2 c_d A</math>. | |||
=== Piezo Crystal === | === Piezo Crystal === | ||
To ignite the air--ethanol mixture, we use piezo crystal. | To ignite the air--ethanol mixture, we use piezo crystal. When tjhe piezo crystal is compressed, it will generate an electric charge which creates a spark. | ||
== References == | == References == | ||
Line 86: | Line 224: | ||
https://www.peacesoftware.de/einigewerte/o2_e.html | https://www.peacesoftware.de/einigewerte/o2_e.html | ||
http://www.users.miamioh.edu/sommerad/NSF%20Files/drag_coefficient_calculation.pdf | |||
DRAG COEFFICIENTS FOR FLAT PLATES , SPHERES, AND CYLINDERS MOVING AT LOW REYNOLDS | |||
NUMBERS IN A VISCOUS FLUID | |||
by | |||
ALVA MERLE JONES | |||
http://www.aerospaceweb.org/question/aerodynamics/q0231.shtml | |||
Ethanol rockets | |||
https://www.youtube.com/watch?v=zTwz6FGobCA Ethanol Rocket - Cool Science Experiment | |||
https://www.youtube.com/watch?v=4s-SZypWxeg Ethanol Explosion - Cool Science Experiment |
Latest revision as of 22:16, 31 May 2021
Introduction
Theory
Ethanol 3d model
-
The Ethanol
-
C2H6O
-
The Ethanol
-1.1712 | 0.2997 | 0.0000 | O |
-0.0463 | -0.5665 | 0.0000 | C |
1.2175 | 0.2668 | 0.0000 | C |
-0.0958 | -1.2120 | 0.8819 | H |
-0.0952 | -1.1938 | -0.8946 | H |
2.1050 | -0.3720 | -0.0177 | H |
1.2426 | 0.9307 | -0.8704 | H |
1.2616 | 0.9052 | 0.8886 | H |
-1.1291 | 0.8364 | 0.8099 | H |
HCH bonds are assumed to be undistorded tetraherdal angle 109.5 degrees. Actually the electrons repeal each other. The HOC bond is 104.5 deg because. . .
Valence shell electron-pair repulsion theory (VSEPR theory). . .
Atom | Radius | Relative radius |
---|---|---|
C | 170 pm | 1.42 |
H | 120 pm | 1.00 |
O | 152 pm | 1.27 |
Use CPK coloring convention, white (hydrogen), black (carbon) and red (oxygen).
Ethanol reaction with Oxygen and Air
Ethanol reaction with oxygen
The molecular weight of ethanol is g/mol, and the molar weight of oxygen is 32 g/mol. The oxygen--ethanol fuel ratio is . We need kg of oxygen to kg of ethanol. The air consists of 23.2 mass-% of oxygen, thus the air--ethanol ratio is .
The volume-% of oxygen in air is 20.9%. The volume of the bottle is ml which gives the amount of oxygen to be ml g of oxygen, which gives the amount of ethanol g ml. OR directly using air--ethanol ratio we have gram. That amount equals to ml mm3. Almost the same result using the methods. See the attached spreadsheet for detailed calculations.
The energy released by burning ethanol is 17.9 kJ/ml. Thus, ml of ethanol releases kJ of energy. This energy is converted into heat, sound and projectile motion (plus others).
Density | At | |
---|---|---|
Oxygen (g) | 1.429 g/l | STP |
Oxygen (g) | 1.314 g/l | 20 °C |
Ethanol (l) | 789.45 g/l | 20 °C |
Air (l) | 1.2041 g/l | 20 °C |
File:Ethanol oxygen combustion.ods
The inner diameter of the rocket bottle is 25 mm. The height of the ethanol in the cap need to be mm.
Exothermic Reaction and Energy release
Specific Heat
The energy is transferred into pressure, sound, etc. The isochoric specific heat of air is kJ/(kgK) at 18 centigrade. At 180 degrees Celsius kJ/(kgK). Thus, the energy released heats
Ideal gas law
The simplest idea is to use ideal gas law for isochoric process constant gives Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p_2 = \frac{T_2}{T_1}p_1} which gives force Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle F = pA} where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A} is the diameter of the nozzle.
Thus, the propulsive force is Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} F &= pA \\ &= \frac{T_2}{T_1}p_1 A \\ &= \frac{T_1 + \Delta T}{T_1}p_1 A \\ &= (1 + \frac{\Delta T}{T_1}) p_1 A \\ &= (1 + \frac{\Delta E/ (C_v m)}{T_1}) p_1 A \\ \end{align} }
Detonation velocity
Burn rate, detonation velocity. https://en.wikipedia.org/wiki/Table_of_explosive_detonation_velocities
http://www.explosionsolutions.co.uk/110411016.pdf
https://link.springer.com/article/10.1007/s00193-015-0554-7 Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v=1500} m/s???
Horizontal(?) accelaration due to rapidly expanding air.
Benjamin Robins:
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle F(x) = \frac{RPAc}{x} }
where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x} distance in barrel, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle R} is the initial ratio of hot gas pressure to atmospheric pressure, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle P} is the atmospheric pressure, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A} is the cross-sectional area of the ball or bore, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle c} is the length of the barrel occupied by the powder charge before ignition.
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_k = \int_0^L F(x) dx \iff v_0^2 = \frac{2RP}{m} \frac{\pi d^2 c}{4}\ln(L/c) }
where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle d} is the barrel diameter (the bore)
The powder change Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p} is given by Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p= \frac{\pi d^2 c}{4}\eta } where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \eta} is the density of gunpowder.
--
https://www.arc.id.au/RobinsOnBallistics.html
The pressure falls as Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 1/x} .
--
Bernoulli?
--
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sum F = 0 }
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle F = ma = pA \iff F = m \frac{d^2x}{dt^2} }
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle F = m \frac{dv}{dt} = mv \frac{dv}{dx} = pA \iff mv dv = pA dx }
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int mvdv = \int_0^L pA dx }
where
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac1L \int_0^L P dx = average pressure = \vec p }
If Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A} is constant
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac12 mv^2 = A \int_0^L p dx = A \vec p L }
then
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v = \sqrt{\frac{2\vec p A L}{m c_f}} } and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec p = \frac{c_f mv^2}{2AL} }
Projectile friction, rotational energy, heat transfer: correction factor Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle c_f} .
The average pressure is 25% of peak pressure
--
http://closefocusresearch.com/calculating-barrel-pressure-and-projectile-velocity-gun-systems
https://www.arc.id.au/CannonBallistics.html
Maximum Flying Distance
By the conservation of energy, all explosive energy is transferred into the kinetic energy.
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_\textrm{explosion}=E_\text{kinetic}}
The distance covered by a projectile with initial velocity $v_0$ is
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle s = \frac{v_0^2}{g}\sin(2\alpha) = \frac{}{}}
The drag coefficient Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle c_d} needs to be found. For the circular disc (a coin) the drag coefficient is almost constant for all velocities (Reynold numbers). The coefficient of drag for a cylinder in this orientation is about 0.81 so long as the length to diameter ratio is greater than 2, see http://www.aerospaceweb.org/question/aerodynamics/q0231.shtml. The cone in either end gives some complications.
The drag equation for the drag force is Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle F_D = \frac12 \rho v^2 c_d A} .
Piezo Crystal
To ignite the air--ethanol mixture, we use piezo crystal. When tjhe piezo crystal is compressed, it will generate an electric charge which creates a spark.
References
https://nptel.ac.in/content/storage2/courses/123106002/MODULE%20-%20I/Lecture%201.pdf
https://www.peacesoftware.de/einigewerte/o2_e.html
http://www.users.miamioh.edu/sommerad/NSF%20Files/drag_coefficient_calculation.pdf
DRAG COEFFICIENTS FOR FLAT PLATES , SPHERES, AND CYLINDERS MOVING AT LOW REYNOLDS NUMBERS IN A VISCOUS FLUID by ALVA MERLE JONES
http://www.aerospaceweb.org/question/aerodynamics/q0231.shtml
Ethanol rockets
https://www.youtube.com/watch?v=zTwz6FGobCA Ethanol Rocket - Cool Science Experiment
https://www.youtube.com/watch?v=4s-SZypWxeg Ethanol Explosion - Cool Science Experiment