Reactor Modeling Tools - STANJAN

GBH Enterprises, Ltd.

Process Engineering Guide:
GBHE-PEG-RXT-820

Tools for Reactor Modeling:
“THE ELEMENT POTENTIAL METHOD FOR CHEMICAL
EQUILIBRIUM ANALYSIS: STANJAN”

Process Disclaimer
Information contained in this publication or as otherwise supplied to Users is
believed to be accurate and correct at time of going to press, and is given in
good faith, but it is for the User to satisfy itself of the suitability of the information
for its own particular purpose. GBHE gives no warranty as to the fitness of this
information for any particular purpose and any implied warranty or condition
(statutory or otherwise) is excluded except to the extent that exclusion is
prevented by law. GBHE accepts no liability resulting from reliance on this
information. Freedom under Patent, Copyright and Designs cannot be assumed.

Refinery Process Stream Purification Refinery Process Catalysts Troubleshooting Refinery Process Catalyst Start-Up / Shutdown
Activation Reduction In-situ Ex-situ Sulfiding Specializing in Refinery Process Catalyst Performance Evaluation Heat & Mass
Balance Analysis Catalyst Remaining Life Determination Catalyst Deactivation Assessment Catalyst Performance
Characterization Refining & Gas Processing & Petrochemical Industries Catalysts / Process Technology - Hydrogen Catalysts /
Process Technology – Ammonia Catalyst Process Technology - Methanol Catalysts / process Technology – Petrochemicals
Specializing in the Development & Commercialization of New Technology in the Refining & Petrochemical Industries
Web Site: www.GBHEnterprises.com

Process Engineering Guide:

Reactor Modeling Tools
STANJAN

CONTENTS

1

SCOPE

2

SUMMARY

3

INTRODUCTION

4

EXAMPLES

4.1 CARBON-RICH C-0 SYSTEM
4.2 EXAMPLE WITH TWO COMPLEX PHASES
4.3 GAS TURBINE ENGINE EXAMPLE
4.4 OTHER APPLICATIONS

APPENDIX

FIGURES
5.1

EXAMPLE RUN LOG FOR CARBON-RICH C-O SYSTEM

5.2

OUTPUT FOR EXAMPLE WITH TWO COMPLEX PHASES

5.3

FIRST STEP IN THE TURBINE EXAMPLE: CALCULATION OF THE
ENTHALPY OF THE REACTANTS

5.4

SECOND STEP IN THE TURBINE EXAMPLE: CALCULATION OF THE
ADIABATIC FLAME TEMPERATURE


5.5

THIRD STEP IN THE TURBINE EXAMPLE: CALCULATION OF THE
NOZZLE EXIT STATE

AVAILABILITY AND IMPLEMENTATION OF STANJAN

REFERENCES


THE ELEMENT POTENTIAL METHOD FOR CHEMICAL EQUILIBRIUM
ANALYSIS: IMPLEMENTATION IN THE INTERACTIVE PROGRAM STANJAN

1

SCOPE

This Process Engineering Guide summarizes STANJAN, an interactive program
for chemical equilibrium analysis.
The program is an interactive program designed for use on either a desktop or a
mainframe computer. Its name is derived from Stanford University, where it was
developed, and the JANNAF thermochemical data tables, from which the basic
data are taken. It relates mol fractions of each species to element potentials for
each independent atom in the system.
The program assumes that the gas phase is a mixture of ideal gases and that the
condensed phases are ideal solutions. Its origin lies in combustion systems.
Section 4 gives a more detailed description and examples

2

SUMMARY

The element-potential method for chemical equilibrium analysis is a powerful
technique that is virtually unknown in the thermodynamics community. It provides
a superior means for solution of complicated problems, especially those involving
several phases. The concept of element potentials is so useful that it should be
included, if not preferred, in any advanced instruction on chemical equilibrium.
This Process Engineering Guide describes the basic method, a new algorithm for
its effective computer implementation, and a general-purpose interactive
freeware program that can be used to solve chemical equilibrium problems in
single- or multiphase systems. The method is extremely flexible and robust, and
the programs have been tested by many users in difficult problems in sooting
combustion, optical materials fabrication, and other multi-phase chemical
problems.


3

INTRODUCTION

The solution of chemical equilibrium problems has posed a tough challenge for
numerical computation. The problem may be formulated in several ways. If one
uses the concept of equilibrium constants, then it is necessary to identify the set
of reactions that take place and to determine the associated equilibrium
constants. One then has to solve a set of nonlinear algebraic equations for the
mol numbers of each species, a difficult task if the system is large. Other
methods based on the minimization of the Gibbs function adjust the mols of each
species, consistent with atomic constraints, until the minimum Gibbs function
state is found. Again, there are many variables involved, and great care must be
taken to be sure that all mols are non-negative. When there are important rare
species in the system, this can be a very difficult task.
In search of a better way to solve combustion-equilibrium problems, the author
reinvented what he later discovered to be a "lost" method, the so-called method
of element potentials. Early development of the method was done by Powell [1].
The RAND method for equilibrium calculation described by Clasen [2] is
essentially an early implementation of the method. White [3] pointed out some
computational advantages of the method. Bigelow [4] extended the nonlinear
programming theory of the method. The author's contribution is the development
of the dual problem and a powerful numerical implementation.
For pre-STANJAN history of the method, see Van Zeggeren and Storey [5]. The
purpose of this paper is to make the method of element potentials known to the
combustion community and to outline an interactive computer program based on
the method that is available for solving chemical equilibrium problems.
The method of element potentials uses theory to relate the mol fractions of each
species to quantities called element potentials. There is one element potential for
each independent atom in the system, and these element potentials, plus the
total number of mols in each phase, are the only variables that must be adjusted
for the solution. In large problems this is a much smaller number than the number
of species, and hence far fewer variables need be adjusted. There are many
other advantages to the element-potential method that quickly become obvious
when one begins to use it. We believe that element potentials should be part
of modern instruction in thermodynamics, but the concept is not widely known or
taught.


The analysis, and the program, assumes that the gas phase is a mixture of ideal
gases and that condensed phases are ideal solutions. These are good
approximations for many practical problems of interest. The concept of element
potentials is not limited to these models, and we believe that it may be very
helpful in dealing with non-ideal systems.
The program, called STANJAN because of its roots at Stanford and its
connection with the JANNAF thermochemical data tables, is an interactive
program designed for use with either desktop or mainframe computers. The
basic data are taken from the JANNAF tables, and data for a selection of species
accompany the program. A companion program, JANFILE, can be used to
prepare data for other species from the JANNAF table data. Both are very robust,
User-friendly interactive programs.
With STANJAN, the user selects the species to be included in each phase of the
system, sets the atomic populations and state parameters, and then STANJAN
solves for the equilibrium state using the method of element potentials.
This is extremely rapid, and, with an i7 processor on a PC, solutions for typical
combustion problems are returned almost immediately. The results include the
composition of each phase (mols and mol fractions), and the thermodynamic
properties of the system, including (if desired) the speed of sound.

Thermodynamic cycle analysis Is easily executed with STANJAN, because the
user may specify the state parameters in a variety of ways, including:
(i)

temperature and pressure,

(ii)

pressure and entropy,

(iii)

enthalpy and pressure same as last run,

(iv)

volume, entropy same as last run.

The equilibrium composition can be calculated, or a frozen composition can be
specified.


STANJAN can be used to compute adiabatic flame temperatures for reactions at
constant pressure (or volume). The approach is first to make a run with the
reactants at the Inlet (or initial) state, which calculates the enthalpy (and energy).
This is followed by a run in which the products of combustion are considered,
with the state specified as having the same pressure and enthalpy (or volume
and energy) as the previous run.
A series of calculations can be made over a matrix of T and P values. These
tabulated results can be stored in a file for later processing; special provisions
are made for creating output tables readable by spreadsheet programs,
particularly Excel.
In summary, STANJAN is a powerful and easy-to-use program for analysis of
chemical equilibrium in single- or multiple-phase systems. The executable PC
program disks are freeware and may be freely copied by any institution for its
use. The FORTRAN source programs are also available at reasonable cost for
users who wish to recompile for other machines or uses.

4

EXAMPLES

4.1

CARBON-RICH C-0 SYSTEM

The first example involves a carbon-rich mixture CO of CO2, O2 and O, with the
possibility of a solid carbon phase.
The user begins by calling STANJAN. After an opening opportunity to see a brief
description of the element potential method, the dialog begins. The user's
responses are underlined in the log in Fig. 5.1. The margin numbers refer to the
notes that follow. This example is similar to the first two examples in Section 3,
and the reader may find it useful to make the comparison.
1)

Note that the first step is to get the data for the species to be used. In this
case the user selects the general-purpose combustion file COMB.DAT,
from which he selects the desired species. When running a lot of problems
with the same species, it is faster to first use the program JANNFILE to set
up a data set containing only the desired species, and then use the singlekeystroke defaults * and # to make the phase selection from this species
set.


2)

In this example the user had made a previous run with the same atoms,
and so is given a choice of using the same populations or making a setup
change. In the first run STANJAN offered two ways to specify the atom
populations, either by giving them directly or by naming some reactants in
the data file and then specifying the number of mols of each put into the
system.

3) The isentropic sound speed, defined as:

is calculated. Equilibrium states will be used to calculate c unless the frozen
composition option is selected, in which case it becomes the sound speed with
frozen chemistry.
4)

STANJAN includes an optional instructive monitor, which can be used to
learn about the method. It is included here to display the features of the
initializer and dual problem solution described in previous sections.

5)

This is what the user would have seen had he chosen to check the atomic
composition.

6)

These are the properties as computed from the JANNAF table data.

7)

The initializer first finds that CO will be the dominant species. The
balancing requires the presence of CO2, in the amount estimated. The
initializer then redistributes the atoms to the phases to provide the target
mol fraction of the balancing species CO.

8)

Note how close the solution is after the initialization. Only two iterations
are required to reach a solution accurate in mol fractions to one part in
108.

9)

These are the data taken from the JANNAF tables.

10)

This is the primary output.


4.2

EXAMPLE WITH TWO COMPLEX PHASES

The next example is taken from the field of fiber optics manufacture. It involves a
complex set of species in the gas phase, and a condensed phase, treated here
as an ideal solution. The JANNAF data and operating conditions were provided
to the author by the AT&T Bell Laboratories. For brevity only the final output is
given in Fig. 5.2. Note that only five iterations were required.

4.3

GAS TURBINE ENGINE EXAMPLE

This example shows how STANJAN can be used to calculate the adiabatic flame
temperature in a gas turbine engine combustor, and then the composition after
isentropic expansion in the turbine nozzle.
The first step is to get the enthalpy of the reactants by a run at the combustor
inlet state. Here we took T = 400°K, P = 6 atm, and assumed that the reactants
were
CH4 + 202 + 7.52 N2
The results are shown in Fig. 5.3.
The next step is to get the adiabatic flame temperature by finding the state of the
products at the same enthalpy and pressure as the last run. The reactants and a
set of products are allowed species. The results are shown in Fig. 5.4.
The final step is to get the temperature following isentropic expansion in the
turbines by finding the equilibrium state for the same species at the case entropy
and a specified pressure (here 1 atm.). The results are given in Fig. 5.5.

4.4

OTHER APPLICATIONS

An early version of STANJAN was used for fiber-optics fabrication analysis by
McAffee et al [7]. This version, which did not use the full dual problem, had
difficulties in converging with multiple phases, and had to be extended by them to
cover ideal solutions. The present version runs their problems very quickly, with
no difficulty.
More recently, McAffee et al [8] studied other such problems using the element
potential method. At that time, the full dual problem had not yet been developed,
so they developed their own algorithm for making the steepest descent.

APPENDIX
FIGURES
5.1
EXAMPLE RUN LOG FOR CARBON-RICH C-O SYSTEM


5.2

OUTPUT FOR EXAMPLE WITH TWO COMPLEX PHASES


5.3

FIRST STEP IN THE TURBINE EXAMPLE: CALCULATION OF
THE ENTHALPY OF THE REACTANTS

5.4

SECOND STEP IN THE TURBINE EXAMPLE: CALCULATION OF THE
ADIABATIC FLAME TEMPERATURE


5.5

THIRD STEP IN THE TURBINE EXAMPLE: CALCULATION OF THE
NOZZLE EXIT STATE

AVAILABILITY AND IMPLEMENTATION OF STANJAN
Special Note:
STANJAN: A general purpose multiphase tool which appears to be no
longer freely available. The author Bill Reynolds died in January 2004. A
version was rewritten to use the CHEMKIN subroutines - under the name
EQUIL or or DEQUIL. The latest version of CHEMKIN is commercial,
although version II used to be free. It is very difficult to find now.


However, STANJAN is available from Colorado State University:
Chemical Equilibrium Calculator
BIOANALYTICAL MICROFLUIDICS PROGRAM
CHEMICAL AND BIOLOGICAL ENGINEERING
CHEMICAL EQUILIBRIUM CALCULATION
This spreadsheet will calculate the chemical equilibrium state of an ideal gas
mixture, subject to necessary constraints on two intrinsic variables. Depending
on the constraint chosen, the calculation invokes STANJAN to minimize the
appropriate derived property—Gibbs energy, Helmholz energy, internal energy,
enthalpy—or maximize entropy for the user-supplied gas mixture.
At this time all species must be gaseous; the calculator will not consider multiple
phases. Thus, while solid and liquid species are present in the database they
cannot be included in the equilibrium calculation.
navier.engr.colostate.edu/~dandy/code/code-4/
Depending on the constraint chosen, the calculation invokes STANJAN to
minimize the appropriate derived property—Gibbs energy, Helmholz energy,
internal ...


REFERENCES
1

Powell, N R, and Sarner, S F, "The Use of Element Potentials in Analysis
of Chemical Equilibrium," Vol. 1, General Electric Co., Report R59/FPD
(1959).

2

Clasen, R J, "The Numerical Solution of the Chemical Equilibrium
Problem," RAND Memorandum RM-4345-PR, AD-609-904 (1965).

3

White, W S, "Numerical Determination of Chemical Equilibrium and the
Partitioning of Free Energy," J. Chem. Phys. Vol. 46 (1967), p. 4171.

4

Bigelow, J R, "Computing Equilibrium Compositions of Ideal Chemical
Systems," TR 70-4, Operations Research House, March, 1970. See also
TR70-1, TR-70-2, TR-70-3 (1970).

5

Van Zeggeren, F, and Storey, S H, The Computation of Chemical
Equilibria, Cambridge Univ. press (1970).

6

Veinott, A F, Mathematical Studies in Management Sciences, McMillan,
New York, pp. 5-8 (1965).

7

McAffee, K 5, Laudise, R A, and Hozack, R S, "Equilibria Concentrations
in the Oxidation of SiCl4 and GeCl4 for Optical Fibers," J Lightwave
Technology, VLT-1, No. 4 (1983), pp. 555-561.

8

McAffee, K B, Walker, J R, Laudise, R A, and Hozack, R S, "Dependence
of Equilibria in the Modified Chemical Vapor Deposition Process on SiCl4
and GeCl4 and O2” J. Am. Ceramic Soc. Vol. 67, No. 6 (1984), pp. 420444.


Reactor Modeling Tools - STANJAN

More Related Content

What's hot (20)

Viewers also liked (20)

Similar to Reactor Modeling Tools - STANJAN (20)

More from Gerard B. Hawkins (20)

Recently uploaded (20)

Reactor Modeling Tools - STANJAN