Programming with ISIM

Numerous Dynamic Simulation Examples, Hands on Experience with ISIM at the PC s. Advance programming with a Workstation, Optimization with ACSL and SimuSolvc... 

Full details of the ISlM digital simulation programming language can be found in the appendix, and by reference to the ISIM programs associated with the simulation examples of Chapter 5. 

Non-linear parameter estimation is far from a trivial task, even though it is greatly simplified by the availability of user-friendly program packages such as a) SIMUSOLV (Steiner et al., 1986), b) ESL, c) a set of BASIC programs (supplied with the book of Nash and Walker-Smith, 1987) or d) by mathematical software (MATLAB). ISIM itself does not supply these advanced features, but ISIM programs can easily be translated into other more powerful languages. 

Each simulation example is identified by a file name and title, and each comprises the qualitative physical description with drawing, the model equation development, the nomenclature, the ISIM program, suggested exercises, sample graphical results and literature references. The diskette in the pocket at the back of the book contains the programs and the ISIM software. 

Program THERM solves the dynamic model equations. The initial values of concentration and temperature in the reactor can be changed after each run using the ISIM interactive commands. The plot statement causes a composite phase-plane graph of concentration versus temperature to be drawn. Note that for comparison both programs should be used with the same parameter values. 

The ISIM program for this example with typical output is shown. Example BENZHYD... 

The ISIM program RUN with KV = 0 represents operation of the reactor under zero venting (V = 0) conditions, i.e., the bursting disc remains closed. With KV > 0 ( 100 - 1000) the program simulates the performance under emergency venting conditions. The program runs stably until the user interacts with HOLD and causes the coolant failure to occur with HOLD, VAL FC = 0 and GO. 

In addition, we would like to thank the following people at the Saarland University Susan Lochow for help with the word processing and Patrick Cernko and Stefan Kiefer for converting most of the older ISIM programs to Madonna. 

Since r is contained in two nonlinear expressions, a global equation for the kinetics can not be derived. The equations can, however, be solved simultaneously by dynamic modeling with the program ISIM [18]. The corresponding ISIM program is shown in Figme 13-20. 

---