Glacier lake outburst floods (GLOFs) initiate with the rapid outburst of a glacier lake, endangering downstream populations, land, and infrastructure. The flow initiates as a mud flow; however, with the entrainment of additional solid material, the flood will often transform into a debris flow. As the run-out slope flattens, the coarse solid material deposits and the flow de-waters. The flow transforms back into a muddy, hyperconcentrated flow of fine sediments in suspension. These flow transitions change the flow composition dramatically and influence both the overall mass balance and flow rheology of the event. In this paper, we apply a two-phase/layer model to simulate flow transitions, solid–fluid phase separations, entrainment, and run-out distances of glacier lake outburst floods. A key feature of the model is the calculation of dilatant actions in the solid–fluid mixture which control flow transitions and phase separations. Given their high initial amount of fluid within the flow, GLOFs are sensitive to slope changes inducing flow transitions, which also implies changes in the flow rheology. The changes in the rheology are computed as a function of the flow composition and do not need any adaptation by ad-hoc selection of friction coefficients. This procedure allows the application of constant rheological input parameters from initiation to run-out. Our goal is to increase the prediction reliability of debris flow modeling. We highlight the problems associated with initial and boundary (entrainment) conditions. We test the new model against the well-known Lake 513 (Peru, 2010), Lake Palcacocha (Peru, 1941), and Lake Uchitel in the Aksay Valley (Kyrgyzstan) GLOF events. We show that flow transition modeling is essential when studying areas that have significant variations in slope.