Optimization of unsteady PDE systems using a multi