A hypothesized diaphragm between the two cylinders is ruptured to begin the flow. Time steps are either accepted or rejected based on whether or not elemental error indicators exceed a prescribed tolerance. Rejected time steps are repeated, subsequent to adaptive h-refinement and rebalancing. Coarsening is essential to keep mesh sizes manageable, as areas of interest move through the domain. Upon h-refinement, the solution is interpolated to the new mesh and the integration continues with, possibly, a new time step.
Accepted steps continue with, possibly, a new global time step.
Flow features compare favorably with experimental and numerical results of Nagamatsu, Carofano, et al. The partitioning is on 16 processors and is obtained using Parallel Sort Inertial Recursive Bisection. Mach number with velocity vectors after the 125th. We have used these tools to examine the Euler flow in intersecting cylinders, simulating a perforated muzzle brake. The tools and techniques referenced here are generic and are being used to solve problems in many areas including materials processing, crystal growth, and biomechanics.