We consider continuous and discontinuous Galerkin time stepping methods of arbitrary order as applied to nonlinear initial value problems in real Hilbert spaces. Our only assumption is that the nonlinearities are continuous; in particular, we include the case of unbounded nonlinear operators. Specifically, we develop new techniques to prove general Peano-type existence results for discrete solutions. In particular, our results show that the existence of solutions is independent of the local approximation order, and only requires the local time steps to be sufficiently small (independent of the polynomial degree). The uniqueness of (local) solutions is addressed as well. In addition, our theory is applied to finite time blow-up problems with nonlinearities of algebraic growth. For such problems we develop a time step selection algorithm for the purpose of numerically computing the blow-up time, and provide a convergence result.This effect is commonly termed (finite-time) blow-up.Galerkin Time Stepping. Galerkin-type time stepping methods for initial-value problems are based on weak formulations. For both the cG and the dG time stepping schemes, the test spaces consist of polynomials that are discontinuous at the time nodes. In this way, the discrete Galerkin formulations decouple into local problems on each time step, and the discretizations can therefore be understood as implicit one-step schemes. Galerkin time stepping methods have been analyzed for ordinary differential equations (ODEs), e.g., in [3,[5][6][7][8]10].A key feature of Galerkin time stepping methods is their great flexibility with respect to the size of the time steps and the local approximation orders, thereby naturally leading to an hp-version Galerkin framework. The hp-versions of the cG and dG time stepping schemes were introduced and analyzed in the works [12,13,15,19]. In particular, in the articles [12,19], which focus on ordinary initial value problems with uniform Lipschitz nonlinearities, the use of the contraction mapping theorem made it possible to prove existence and uniqueness results for discrete Galerkin solutions, which are independent of the local approximation orders. We emphasize that the hp-approach is well-known for its ability to approximate smooth solutions with possible local singularities at high algebraic or even exponential rates of convergence; see, e.g., [4,13,14,18] for the numerical approximation of problems with start-up singularities.