A tightly coupled fluid-structure interaction (FSI) solution technique incorporating fluid and solid mechanics, phase change and chemical reactions is presented. The continuum equations are solved with a cell-centered, multi-material ICE solution method. This formulation is integrated with a Lagrangian, particle based, solid mechanics technique, known as the Material Point Method, as described by Kashiwa et al. [1] and Guilkey et al. [2]. The combined method can handle large deformations and phase change within a single grid, without the need of separate domains for fluids and solids, or the passing of boundary conditions. This paper discusses algorithmic issues involved in accounting for chemical reactions and phase transition among material phases (e.g., solid --> gas). Validation is presented as are simulations showing large deformation with phase change. These simulations were performed within a computational framework that contains tools for parallelization, performance analysis, data management, algorithm integration, and data visualization. Features of this framework are described.