CFD Online Discussion Forums

CFD Online Discussion Forums (
-   Main CFD Forum (
-   -   Two solutions of Sod's shock tube problem (

orxan.shibli January 23, 2013 18:52

Two solutions of Sod's shock tube problem

When I was googling I found two solutions of the Sod's shock tube problem with the same initial conditions. In the given references, especially density profiles are different. My solution is the same as the second reference. Do you have any idea how there could be two solutions with the same initial conditions? The references are:


praveen January 24, 2013 00:42

This is strange. I always get the first solution for this problem. This is a standard test case and you can see sample results in most books like Toro which show the first solution. Toro's book also has exact Riemann solution which you use for this problem.

orxan.shibli January 25, 2013 15:26

Interestingly I get the first result with my 1D code but I get second one with my 3D code. I cannot figure out why this happens.

AeroMike January 30, 2013 09:27

This is very strange. I concur that the first solution is the same as the one in Toro and is the result I obtain with my code using 1D or 3D tests.

If the results differs between 1 and 3D solvers is it to do with the storage and/or order of flux splitting?

The second page is very interesting, though I'm somewhat concerned for anyone trying to replicate those results, for a 1D shock tube okay we should be able to do it, but for example the Rayleigh-Taylor case, this result is very dependent on the order of accuracy of the code (and even which implementation e.g MUSCL or WENO etc) so I think while you would reproduce the basic bubble/spike image to get the exact same detail is most unlikely. Liska and Wendroff have a good paper showing it "Comparison of several difference schemes on 1D and 2D test problems for the Euler equations".

Sorry, a bit off topic but I'm keen to hear a solution to the OP's question and reignite some interest in methods :-)


orxan.shibli January 30, 2013 16:32

I tried 2nd order accuracy in space and 4th order in time however that only makes the solution to look sharper and "better". Still my answer is same as the second resource.

AeroMike January 30, 2013 16:35

Sorry Orxan, my comments weren't really directed at the solution of the Sod problem, rather the other solutions on that page. I was just worried that if the latter data was a bit misleading the former could be too.

Certainly I wouldn't expect the order of accuracy to effect the Sod solution very much (definitely not in the manner here). This second solution is a mystery to me and I am eager to find the answer!


praveen January 31, 2013 00:51

Are the 3-d simulations in cartesian or spherical coordinates ? Perhaps the 3d simulations are obtaining radial solutions which would be different from plane shock tube solutions.

immortality February 2, 2013 07:45

are there any experimental results?

immortality February 2, 2013 09:51

the diagrams is told that are in t=.1s but it seems too much time and at this time equilibrium is almost reached.what time should be set accurately?

bob91 July 7, 2013 06:40

ı want to make difference of pressure to be shock but ı don how to do this, so can you help me ?

immortality February 6, 2014 11:21

is there any data about after reflection of shock from the wall? (pressure,density,temperature)

gaurav11 November 12, 2014 04:48

Instability in solution
Hi all,

I am trying to simulate sod's shock tube with Mccormack scheme (second order). I am getting infinite values, and solution becomes unstable very soon. Anyone has any idea what is going wrong.

All times are GMT -4. The time now is 16:26.