Numerical Sedov-Taylor (blast-wave) solution

[jdwood983]

Homework Statement

I was asked to model the blast-wave of an explosion using Sedov's solution in MATLAB, but I'm not really sure where to begin. Most of the programming I have done has been

x=linspace(1,1000,1000);
for i=1:1000
     y(i)=x(i)^2;
end

type of code, and this one seems a lot more complex (I found a version of this code in C [sorta familiar with C, to some extent all object oriented languages are the same] which is long and has many functions in it which leads me to the 'complex' conclusion here).

Any suggestions on where to begin? If not, any suggestions on a better book than Sedov's Similarity and Dimensional Methods in Mechanics to get me started? Or maybe both?

Homework Equations

<br /> \frac{r}{r_2}=\left[\frac{(\nu+2)(\gamma+1)}{4}V\right]^{-2/(2+\nu)}\left[\frac{\gamma+1}{\gamma-1}\left(\frac{(\nu+2)\gamma}{2}V-1\right)\right]^{-\alpha_2}\left[\frac{(\nu+2)(\gamma+1)}{(\nu+2)(\gamma+1)-2[2+\nu(\gamma-1)]}\left(1-\frac{2+\nu(\gamma-1)}{2}V\right)\right]^{-\alpha_1}<br />
and a couple others

The Attempt at a Solution

Not sure where to begin, so not much of a solution attempt can exist. I did try using the equation above and changing V from its lower bound of 2/([\nu+2]\gamma) to its upper bound of 4/((\nu+2)(\gamma+1)) but that gave me bad data.

Any help would be appreciated.

 

[CFDFEAGURU]

Blandford and Thorne's free notes from CalTech discuss the blast wave. It doesn't contain any code but it might help.

See the link below

http://www.pma.caltech.edu/Courses/ph136/yr2008/0616.1.K.pdf

Thanks
Matt