Delivered-To: skoge@chembio.ntnu.no X-Sender: engelien@pop.chembio.ntnu.no X-Mailer: QUALCOMM Windows Eudora Version 4.3.2 Date: Mon, 17 Sep 2001 15:40:23 +0200 To: skoge@chembio.ntnu.no From: engelien@pop.chembio.ntnu.no Subject: Matlab filer ----------------------------------------------------------------- Matlab Models for Multi-Effect Distillation Columns (paper presented at SIMS 2001 ----------------------------------------------------------------- Sigurd, Her er Matlab filene mine. Filene for de dynamiske simuleringene ligger i zip fila Dynamic_simulations.zip og optimaliseringsfilene i Self_optimizing.zip Følgende er 'summary' av filene: The following models are available for simulation with a multi-effect distillation column: Dynamic simulation run_multiefffect a program that integrates the derivatives from multieffect.m, and plots the column profiles and responses multieffect a program that calls and connects the two column models colamod1.m and colamod2.m and returns the derivatives xprime colamod1 the dynamic model for the HP column colamod2 the dynamic model for the LP column Optimisation run_optimisation runs the program optimisation.m, finds the optimum for the nominal (d0) operating point and the two disturbances d1 and d2 optimisation the program with the optimisation routine fmincon it calls the programs cola_lv_fun.m and cola_lv_nonlcon.m cola_lv_fun the objective function, calls colamod1 and colamod2 cola_lv_nonlcon the equality and inequality constraints, calls colmod1 and colamod2 colamod1 the dynamic model for the HP column colamod2 the dynamic model for the LP column Calculate Loss run_loss calculates the loss L for the disturbances d0, d1 and d2 for each of the controlled variables, calls the program optimisation.m optimisation the program with the optimisation routine fmincon it calls the programs cola_lv_fun.m and cola_lv_nonlcon.m cola_lv_fun the objective function, calls colamod1 and colamod2 cola_lv_nonlcon the equality and inequality constraints, calls colmod1 and colamod2 colamod1 the dynamic model for the HP column colamod2 the dynamic model for the LP column Vary Heatload vary_heatload program for varying heatload to HP column secant_cont calls multieffect_heatload and finds solution to f(x)= 0 while varying heatload, VB multieffect_heatload calls colamod1 and colamod2 colamod1 the dynamic model for the HP column colamod2 the dynamic model for the LP colum Self_optimizing.zip Dynamic_simulations.zip *************************************************** Hilde K. Engelien Dr.ing student Institutt for Kjemisk Prosessteknologi NTNU Tel : + 47 73 59 58 67 e-mail : engelien@chembio.ntnu.no ***************************************************