CFD Online Discussion Forums

CFD Online Discussion Forums (
-   System Analysis (
-   -   1D solution of shock tube (

immortality May 30, 2013 10:08

1D solution of shock tube
shock tube is a basic CFD problem for testing systems with discontinuities.
its valuable that provide a analytic program for 1D inviscid shock tube that gives flow properties at each sections whenever the user wants,in any form like a Matlab program or what form that anyone interested thinks its good is worthy.

kaushik.kiit June 18, 2013 06:24

reducing CFL no.
1 Attachment(s)
I have run the shock tube problem for various schemes and I found that when the CFL number is reduced, the error increases...can you explain why??

BTW the analytical solution of shock tube is attached here

immortality June 18, 2013 06:36

thats a valuable work.could you add some explanations in different parts?
you compare error with what?

kaushik.kiit June 18, 2013 06:45

I have found the error in pressure and density by evaluating the diiference between the numerical solution and the exact solution and I have carried out the process for two values of CFL no, 0.5 and 0.1. The errors associated with CFL 0.1 are larger than CFL 0.5 and I would like to know why....

immortality June 18, 2013 06:57

why this error occurs?

ehsan@Ehsan-com:~/Desktop$ gcc analytical_WR.c -o WR
analytical_WR.c:4:18: fatal error: conio.h: No such file or directory
compilation terminated.

immortality June 27, 2013 09:01

ShockTube program execution error
1 Attachment(s)
what does errors mean?and how to resolve?
can anyone help?

chegdan June 27, 2013 09:47

1 Attachment(s)
Added the stdlib.h header and compiled it correctly with


gcc -o output ShockTube.CPP -lm
Run code with



immortality June 27, 2013 11:25

1 Attachment(s)
Thanks and whats the problem with this file?could you have a look into it?

chegdan June 27, 2013 12:25

Did you try to Google the error?


gcc -o output shockTube2.CPP -lm
shockTube2.CPP:4:18: fatal error: conio.h: No such file or directory
compilation terminated.

here is a link to what i found.

All times are GMT -4. The time now is 10:13.