Sedov-Taylor blast: ambient rho=1,p=1e-5; total energy E=1 deposited as high pressure in cells within r0 = 3.5*dx of the domain centre (cylindrical).
Energy normalisation uses the DISCRETE cell count (not pi*r0^2) so that exactly E=1 is deposited on any grid, avoiding the O(15%) shortfall that arises from rounding when the disk boundary crosses cell centres.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(solver_state_2d_t), | intent(inout) | :: | state | |||
| type(config_2d_t), | intent(in) | :: | cfg |