# Euler Equations, Sod shock tube & conservation

 Register Blogs Members List Search Today's Posts Mark Forums Read

 June 27, 2014, 23:15 Euler Equations, Sod shock tube & conservation #1 New Member   Join Date: Jun 2014 Posts: 1 Rep Power: 0 Is momemtum 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 where pressure is related to the conserved quantities by The Sod shock problem splits the domain into two regions separated by a density and pressure discontinuity with initial velocity zero. That is, 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!

 July 17, 2014, 13:30 #2 Member   Join Date: Jul 2013 Posts: 49 Rep Power: 6 In the way I see Sod shock tube, there is a seal between the two region (left and right). Since the pressure on both side are not equal, there is a net force applied by the fluid on the seal. Using Newton's 3rd law, the seal also applies a net force on the fluid. Momentum is only conserved if no external forces are applied on the system, which is not the case here. As soon as the seal is removed from the system (or broken) there is a discontinuity in pressure and density. The very first fluid element, adjacent to this pressure discontinuity, moves at a very high velocity, with all the momentum of the domain.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Vinoo_P FLUENT 18 July 29, 2016 18:40 jojo81 Main CFD Forum 0 October 15, 2012 12:05 Astaria OpenFOAM Running, Solving & CFD 5 March 4, 2012 04:07 technophobe Main CFD Forum 5 April 28, 2009 15:11 Jan Ramboer Main CFD Forum 2 August 19, 1999 01:58

All times are GMT -4. The time now is 21:01.