Masters Theses

Date of Award


Degree Type


Degree Name

Master of Science


Nuclear Engineering

Major Professor

J.C. Robinson

Committee Members

T.W. Kerlin, Hale C. Roland


A numerical algorithm is formulated to solve the first order, nonlinear differential equations that describe a multistage flash evaporator. The nonlinearities appearing in the formulation of the algorithm are products of up to three terms with each term being a dependent variable raised to some power.

To develop the algorithm, the first order differential equations are written in integral form. The dependent variables are then assumed to have a purely exponential dependence over a finite time step thereby allowing for the explicit integration of all terms. The solution of the differential equations is then reduced to the determination of the exponential dependences. The exponential dependences are determined by an iterative method.

A computer code based upon the aforementioned algorithm was written. Before the algorithm was used to obtain solutions to a flash evaporator system, it was applied to several differential equations with known solutions. The algorithm was then used to obtain solutions for two perturbations in the twenty-third order system that describes a three stage flash evaporator. These solutions are compared with solutions obtained by other methods.

Files over 3MB may be slow to open. For best results, right-click and select "save as..."