- #1
Antigravity324
- 1
- 0
Is momentum conserved?
I am considering the Euler equations in conservative form and solving the Sod shock tube problem I have written a Godunov finite volume type solver. It solves for density ρ, momentum ρu, and total energy E; therefore, I would expect all of these quantities to be conserved wrt time. Density and total energy are indeed conserved, however, momentum is not.
The Euler equations are given by
<br /> \begin{bmatrix}<br /> \rho \\ \rho u \\ E<br /> \end{bmatrix}_t<br /> +<br /> \begin{bmatrix}<br /> \rho u \\ \rho u^2 + p \\ u(E +p)<br /> \end{bmatrix}_x<br /> = 0<br />
where pressure is related to the conserved quantities by
<br /> p = (\gamma - 1)(E - 0.5 \rho u^2)<br />
The Sod shock problem splits the domain into two regions separated by a density and pressure discontinuity with initial velocity zero. That is, \rho_L = 1, \rho_R = .125; p_L=1, p_R =.1; u_L=u_R =0 These initial conditions imply that ρu=0, therefore, momentum should be zero throughout the simulation.
The solution profiles are well known and can be found here or here. We see that ρ >0 and u > 0, therefore, there is no way for momentum ρu = 0 (which it should be from the intial conditions). As a result, I do not even see why it is reasonable to expect that momentum would be conserved. A paper by Sod himself which surveys some methods for solution, on page 20, list a table which shows momentum is not conserved, but no reason why. I generally do not work in this area, so maybe I am missing something basic. Can anyone shed some light on this? Thanks!
I am considering the Euler equations in conservative form and solving the Sod shock tube problem I have written a Godunov finite volume type solver. It solves for density ρ, momentum ρu, and total energy E; therefore, I would expect all of these quantities to be conserved wrt time. Density and total energy are indeed conserved, however, momentum is not.
The Euler equations are given by
<br /> \begin{bmatrix}<br /> \rho \\ \rho u \\ E<br /> \end{bmatrix}_t<br /> +<br /> \begin{bmatrix}<br /> \rho u \\ \rho u^2 + p \\ u(E +p)<br /> \end{bmatrix}_x<br /> = 0<br />
where pressure is related to the conserved quantities by
<br /> p = (\gamma - 1)(E - 0.5 \rho u^2)<br />
The Sod shock problem splits the domain into two regions separated by a density and pressure discontinuity with initial velocity zero. That is, \rho_L = 1, \rho_R = .125; p_L=1, p_R =.1; u_L=u_R =0 These initial conditions imply that ρu=0, therefore, momentum should be zero throughout the simulation.
The solution profiles are well known and can be found here or here. We see that ρ >0 and u > 0, therefore, there is no way for momentum ρu = 0 (which it should be from the intial conditions). As a result, I do not even see why it is reasonable to expect that momentum would be conserved. A paper by Sod himself which surveys some methods for solution, on page 20, list a table which shows momentum is not conserved, but no reason why. I generally do not work in this area, so maybe I am missing something basic. Can anyone shed some light on this? Thanks!
