# Two solutions of Sod's shock tube problem

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

 January 23, 2013, 18:52 Two solutions of Sod's shock tube problem #1 New Member   Orxan Shibliyev Join Date: Aug 2011 Posts: 18 Rep Power: 8 Hi 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: 1. http://www.csun.edu/~jb715473/exampl...1d.htm#density 2. http://irfu.cea.fr/Projets/Site_hera...est_suite.html

 January 24, 2013, 00:42 #2 Super Moderator     Praveen. C Join Date: Mar 2009 Location: Bangalore Posts: 273 Blog Entries: 6 Rep Power: 11 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. __________________ http://cpraveen.github.io http://twitter.com/cfdlab

 January 25, 2013, 15:26 #3 New Member   Orxan Shibliyev Join Date: Aug 2011 Posts: 18 Rep Power: 8 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.

 January 30, 2013, 09:27 #4 New Member   Join Date: Aug 2012 Posts: 22 Rep Power: 7 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 :-) Mike

 January 30, 2013, 16:32 #5 New Member   Orxan Shibliyev Join Date: Aug 2011 Posts: 18 Rep Power: 8 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.

 January 30, 2013, 16:35 #6 New Member   Join Date: Aug 2012 Posts: 22 Rep Power: 7 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! Mike

 January 31, 2013, 00:51 #7 Super Moderator     Praveen. C Join Date: Mar 2009 Location: Bangalore Posts: 273 Blog Entries: 6 Rep Power: 11 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. __________________ http://cpraveen.github.io http://twitter.com/cfdlab

 February 2, 2013, 07:45 #8 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 20 are there any experimental results?

 February 2, 2013, 09:51 #9 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 20 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?

 July 7, 2013, 06:40 shock_tube #10 New Member   Join Date: Mar 2013 Location: turkey Posts: 8 Rep Power: 6 hi, ı want to make difference of pressure to be shock but ı don how to do this, so can you help me ?

 February 6, 2014, 11:21 #11 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,210 Rep Power: 20 is there any data about after reflection of shock from the wall? (pressure,density,temperature) __________________ Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked.

 November 12, 2014, 04:48 Instability in solution #12 New Member   gaurav kumar yadav Join Date: Jul 2013 Posts: 19 Rep Power: 6 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.

May 20, 2016, 20:51
#13
New Member

Matt
Join Date: May 2016
Posts: 2
Rep Power: 0
Quote:
 Originally Posted by gaurav11 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.
Hey did you ever find the solution to this?

If so, can I see your code?
I am having the same issue and it is driving me crazy.

May 20, 2016, 21:32
#14
Member

LUQILIN
Join Date: May 2016
Posts: 34
Rep Power: 3
Quote:
 Originally Posted by gaurav11 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.
Please give more details. Based on my previous experiences, maybe:
1. CFL is too large.
2. Incorrect way to compute flux.
3. inappropriate boundary condition.
4. simple mistake: for example , take "i" as “j".

If you code blows very soon, there must be some fatal errors. I suggest you deal with the simple case first. Give an artificial initial conditions, say let density be a linear function in the whole computational domain. Debug step by step, see what you can dig.

May 20, 2016, 21:41
#15
Member

LUQILIN
Join Date: May 2016
Posts: 34
Rep Power: 3
Quote:
 Originally Posted by AeroMike 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 :-) Mike
Hello, Mike
I would like to ask if you have any experiences dealing discontinuities with minmod limiters. I got some problems with my code for 1D shock nozzle. I can get a convergent solution without limiters. But Lax-Wendroff oscillations occur at the discontinuity. Unfortunately, my code blows after adding a minmod limiter.
By the way, I am trying a high-order method named Correction Procedure via Reconstruction (CPR). Finite Volume method coincides well with the exact solution.
Good luck.

 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 manukamin Main CFD Forum 3 March 22, 2016 02:02 Accelerator Main CFD Forum 0 August 26, 2012 15:52 Astaria OpenFOAM Running, Solving & CFD 5 March 4, 2012 04:07 AB Siemens 0 November 21, 2004 19:43 Ravi FLUENT 2 October 27, 2004 12:05

All times are GMT -4. The time now is 03:40.