[dnm2017]

Hi Everyone,

I'm attempting to setup a 1D, spherically symmetric shock expansion in OpenFoam. To do this I have setup a wedge type mesh with one cell in the y and y directions and say 100 in the x direction, so basically a pyramid shape with the origin of the blast at the apex. I have set the boundary condition at the origin as empty, as suggested in other 2D axi-symmetric examples. The other boundaries are set as wedge types apart from the end which is transmissive. the angle is around 5 degrees as suggested. The solver I will be using is rhoCentralFoam

When I set initial pressures in the first few cells I see a shock propagating but it appears to be 'leaking' out of the empty boundary at the apex.

Can anyone give any advice how best to set up this kind of problem, is it appropriate to collapse all points at the origin into one, I have seen in other (autodyn) tutorials the point is clipped, but I don't know what effect this has.

Thanks in advance for your advice