Combustion Jet

The equation of motion for rockets is m r v ˙ r = u e m ˙ r , m_\textrm{r}\,\dot{v}_\textrm{r}= u_\textrm{e}\, \dot{m}_\textrm{r}, where u e u_\textrm{e} is the speed of the exhaust shooting out the back. The final velocity of the rocket is a few multiples of the exhaust speed, so it's quite important, but in physics textbooks, u e u_\textrm{e} is usually treated like a freely adjustable parameter whose origin and precise value is a detail left for engineers to worry about. With a little bit of approximation, we can do better than that.

Estimate u e u_\textrm{e} for the combustion of methane ( C H X 4 ) , (\ce{CH4}), the propellant used by SpaceX for its Raptor rocket engine program.

Details and Assumptions:

  • The specific enthalpy of combustion for methane is Δ c H 890 kJ / mol . \Delta_c H \approx \SI[per-mode=symbol]{890}{\kilo\joule\per\mole}.
  • Assume that all the energy released in combustion goes into the kinetic energy of the gas.
  • Beyond its performance, methane is nice because it could potentially be produced on other planets via the Sabatier reaction C O X 2 + 4 H X 2 C H X 4 + 2 H X 2 O . \ce{CO2 + 4H2 -> CH4 + 2H2O}.
50 m / s \approx\SI[per-mode=symbol]{50}{\meter\per\second} 500 m / s \approx\SI[per-mode=symbol]{500}{\meter\per\second} 5 000 m / s \approx\SI[per-mode=symbol]{5000}{\meter\per\second} 50 000 m / s \approx\SI[per-mode=symbol]{50000}{\meter\per\second} 500 000 m / s \approx\SI[per-mode=symbol]{500000}{\meter\per\second}

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

3 solutions

Josh Silverman Staff
Apr 22, 2018

The combustion of methane proceeds according to C H X 4 + 2 O X 2 C O X 2 + 2 H X 2 O , \ce{CH4 + 2O2 -> CO2 + 2H2O,} forming one molecule of carbon dioxide and two molecules of water for each molecule of methane combusted. If we assume that 100 % 100\% of the energy that's released goes into the kinetic energy of the products, then we have E total = 1 2 m C O X 2 u C O X 2 2 + 2 × 1 2 m H X 2 O u H X 2 O 2 . E_\textrm{total} = \frac12 m_{\ce{CO2}} u_{\ce{CO2}}^2 + 2\times \frac12 m_{\ce{H2O}}u_{\ce{H2O}}^2. Now, the atomic mass of carbon dioxide and water are a bit different ( 44 44 and 18 , 18, respectively) but by multiplying and dividing the right hand side by the sum of the masses, we can find a characteristic value for the speed of a given exhaust molecule. E total = 1 2 m C O X 2 u C O X 2 2 + 2 × m H X 2 O u H X 2 O 2 m C O X 2 + 2 m H X 2 O ( m C O X 2 + 2 m H X 2 O ) 2 E total = u e 2 ( m C O X 2 + 2 m H X 2 O ) \begin{aligned} E_\textrm{total} &= \frac12 \frac{m_{\ce{CO2}} u_{\ce{CO2}}^2 + 2\times m_{\ce{H2O}}u_{\ce{H2O}}^2}{m_{\ce{CO2}} + 2m_{\ce{H2O}}}\left(m_{\ce{CO2}} + 2m_{\ce{H2O}}\right) \\ 2E_\textrm{total} &= \langle u_\textrm{e}^2\rangle\left(m_{\ce{CO2}} + 2m_{\ce{H2O}}\right) \end{aligned} where we've defined u e 2 \langle u_\textrm{e}^2\rangle as the mass-weighted mean square velocity of the exhaust molecules. Taking the root, we have an estimate for the maximum exhaust speed for the combustion of methalox (methane + liquid oxygen):

u e 2 E total m C O X 2 + 2 m H X 2 O 2 × 890 000 J mol 1 / 0.08 kg mol 1 < 20 × 5 4 × 1 000 000 J k g 1 = 25 × 1 000 000 J k g 1 = 5 000 m s 1 . \begin{aligned} u_\textrm{e} &\approx \sqrt{\frac{2E_\textrm{total}}{m_{\ce{CO2}} + 2m_{\ce{H2O}}}} \\ &\approx \sqrt{2\times\SI{890000}{\joule\per\mole}/\SI{0.08}{\kilo\gram\per\mole}} \\ &< \sqrt{20\times\frac{5}{4}\times\SI{1000000}{\joule\per\kilo\gram}} \\ &= \sqrt{25\times\SI{1000000}{\joule\per\kilo\gram}} \\ &= \SI{5000}{\meter\per\second}. \end{aligned}

In practice, u e u_e is something less than 5 km s 1 , \SI{5}{\kilo\meter\per\second}, topping out near 3.6 km s 1 \SI{3.6}{\kilo\meter\per\second} in the vacuum of space, and a bit less on Earth. But this gets us in the ballpark.

Is there a systematic way of determining what fraction of the combustion energy ends up as kinetic energy of the reaction products? Or is it just an empirical test? And when I say "kinetic energy", I'm referring to ordered kinetic energy, rather than heat.

Steven Chase - 3 years, 1 month ago

Log in to reply

Steve, see my solution. I argue that the average velocity is the same for both molecules (this is the "ordered kinetic energy" in your terms), and there is also a thermal energy (heat), carried by disordered motion.

Laszlo Mihaly - 3 years, 1 month ago

Log in to reply

When I was reading some combustion engine lit to check if this estimate was anywhere close, I saw calculations that basically do what Laszlo is saying. If I come across it again I'll link.

Josh Silverman Staff - 3 years, 1 month ago

How can this be valid, would not the motion depend on the mass of the gas being expelled and the mass of the rocket?

Dennis Gelinas - 3 years, 1 month ago

Log in to reply

If the question is how fast will the rocket go, the masses are important. However, here the question was the speed of the rocket exhaust - for that one does not need masses.

Laszlo Mihaly - 3 years, 1 month ago

This is intermediate, Are you kidding me? I would have to look up such a big MASS of knowledge to be able to solve this...

Christian Herrmann - 3 years, 1 month ago

How did u get the mass-weighted mean square velocity of the exhaust molecules???Like I mean how did u figure out the formula????

erica phillips - 3 years, 1 month ago

Level 1? I did not even understand the question.

A Former Brilliant Member - 2 years, 4 months ago
Laszlo Mihaly
Apr 23, 2018

I solved it by assuming that u C O 2 = u H 2 O u_{CO_2}=u_{H_2O} , to get 4717m/s.

The exhaust gas is an ionized plasma, and there is a strong interaction between the particles. (We know that from the fact that the exhaust glows.) If there is a strong interaction between the particles, they will come to thermal equilibrium. The actual velocity of a given particle is u = u t h e r m a l + u a v e r a g e u=u_{thermal}+u_{average} , where the average is the same for all particles and the thermal velocity depends on the mass, calculated from the equipartition theorem, u t h e r m a l = 3 k B T m u_{thermal}=\sqrt{\frac{3k_B T}{m}} , where k B k_B is the Boltzmann constant, T T is the temperature and m m is the mass of the molecule. The number quoted above is valid only if we neglect the energy carried away by thermal motion. Taking into account the disordered motion in the hot gas (and the energy carried by the excited state of the ions) will make the actual exhaust velocity smaller than the estimated value, bringing it closer to the 3.6km/s value quoted by Josh.

Pau Cantos
May 4, 2018

I didn't know many about this type of reactions so what I did was:

Atomic mas for Oxigen: 16.

Atomic mas for Carbon: 12.

Atomic mas for Hidrogen: 1.

Atomic mas for all the reactives or products: 16 2+12 1+8*1=52.

Then with the converssion 52g/mol:

890kJ/mol (1mol/52g) (1000g/1kg)*(1000J/1kJ) to get the units J/kg=(m/s)^2. Then I suposed the result would be u^2 so I took the square root and got the correct answer with u=4137.07...m/s~5000m/s.

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...