CFD simulation of blast effects in a confined environment. The domain is a cylinder of 1 meter in diameter and 1 meter in height, with rigid walls, modeled as 2D axisymmetric; the initial bulk conditions are 1 bar and 300 K; the blast is generated by patching a high pressure (1000 bar) and high temperature (2300 K) spherical region, at the center of the domain. The simulation was performed using a density-based, unsteady explicit solver, with 3rd order spatial discretization and 4th order temporal discretization schemes, on a perfectly uniform mesh of 1 mm resolution and a time step size of 2e-08 s. Material properties correspond to dry air, density was calculated using the ideal gas law, and all other properties are temperature-dependent.
The first animation represents the flow field colored by temperature, and the second one, flow field colored by Mach number.
It is very interesting to see the interaction between the incident and reflected shock waves that generates strong turbulence, starting from small initial disturbances in the flow.