HyPar has the capability to run ensemble simulations on multiple domains that of varying sizes and initial solution using the same discretization methods, time step sizes, final time, physics parameters, and boundary conditions. The following examples demonstrate this functionality:

1D Linear Advection - Sine Wave - Ensemble Simulation
2D Euler Equations - Isentropic Vortex Convection - Ensemble Simulation

Other examples can also be run with this feature; the modifications needed are as follows:

Need to have an input file simulation.inp that specifies the number of simulations.
Modify the code to generate the intial solution (and exact solution, if available) to write out initial solution files for each simulation with the filenames properly indexed, i.e., initial_<n>.inp, where n = 0, 1, ..., number of simulations-1.
Modify size and iproc in solver.inp to specify the grid dimensions and MPI rank distribution for each simulation.
Note that multiple output files will be written, one set for each simulation.