The purpose of this paper is to adapt the methodology of dynamical systems to the study of a wide class of mathematical models currently used to solve problems of steady 1-D two-phase flows. A detailed study of the geometrical features of the ensemble of solutions is then used for two purposes. First, it makes it possible to understand the physical characteristics of such flows without the need to produce complete solutions. Secondly, the methodology gives valuable indications as to how to supplement computer codes which become inadequate in the neighborhood of singular points. This leads to avoidable numerical difficulties and incorrect interpretations. This methodology is particularly useful in the study of the phenomena of cooking. The important contribution of this general analysis is to show that, regardless of the number of equations in the model, the generic solutions involve only three types of singular points: saddle points, spirals or nodes. By way of example, the method is applied to the mathematically simplest case, the homogeneous flow model with adiabatic boundary conditions. The channel consists of a convergent portion with transition to a divergent portion through a smooth throat. The divergent portion possesses an inflection point after which the rate of area divergence decreases. The fluid is a mixture of water and steam in thermodynamic equilibrium. The expected saddle point, located just downstream from the throat, is followed by an unexpected spiral. The phase space consisting of pressure P, enthalpy h and spatial coordinate z divides itself into four distinct areas. In area A of figure 13 all solutions are single-valued, possess a pressure minimum and correspond to physically acceptable subcritical flows. In area B the solutions of the differential equations (all possessing a turning point) are irrelevant for the physical problem at hand. Through the saddle point there passes a trajectory which is unique in the subcritical portion and continues with two branches, one subcritical ending at pressure P′, the other supercritical, ending at pressure P″. Area C between those branches contains states which are not described by the postulated model. It is conjectured that flows against a back-pressure, P″ < P a < P′, evidently observed in nature, must presumably be described by an extension to the model involving enhanced entropy production. The remaining area D contains supercritical state points and is of no present interest.