Context Navigation

sample_column.mso @ 247

Last change on this file since 247 was 247, checked in by Rafael de Pelegrini Soares, 16 years ago
Updated solver options for the new language
Property svn:eol-style set to `native` Property svn:keywords set to `Id`
File size: 12.2 KB

Rev	Line
[86]	1	#*-------------------------------------------------------------------
	2	* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
	3	*
	4	* This LIBRARY is free software; you can distribute it and/or modify
	5	* it under the therms of the ALSOC FREE LICENSE as available at
	6	* http://www.enq.ufrgs.br/alsoc.
	7	*
	8	* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
	9	* from http://www.rps.eng.br Copyright (C) 2002-2004.
	10	* All rights reserved.
	11	*
	12	* EMSO is distributed under the therms of the ALSOC LICENSE as
	13	* available at http://www.enq.ufrgs.br/alsoc.
	14	*
	15	*--------------------------------------------------------------------
	16	* Sample file for column model
	17	*--------------------------------------------------------------------
	18	*
[88]	19	* This sample file needs VRTherm DEMO(www.vrtech.com.br) to run
	20	* SectionColumn_Test and needs VRTherm full to run the distillation
	21	* column flowsheet.
[86]	22	*
	23	*----------------------------------------------------------------------
	24	* Author: Paula B. Staudt
	25	* $Id: sample_column.mso 247 2007-04-24 13:44:18Z rafael $
	26	--------------------------------------------------------------------#
	27
[1]	28	using "stage_separators/column";
[111]	29	using "controllers/PIDIncr";
[1]	30
	31	# column section with 2 trays
	32	FlowSheet SectionColumn_Test_with2tray
	33	PARAMETERS
[213]	34	PP as Plugin(Brief="Physical Properties",
	35	Type="PP",
	36	Components = [ "isobutane", "benzene"],
	37	LiquidModel = "PR",
	38	VapourModel = "PR"
	39	);
[1]	40	NComp as Integer;
	41
	42	SET
	43	NComp = PP.NumberOfComponents;
	44
	45	DEVICES
	46	sec as Section_Column;
[125]	47	feed as liquid_stream;
	48	reb as vapour_stream;
	49	cond as liquid_stream;
[1]	50	zero as stream;
	51
	52	CONNECTIONS
	53	feed to sec.trays(2).Inlet;
	54	zero to sec.trays(1).Inlet;
	55	reb to sec.trays(2).InletV;
	56	cond to sec.trays(1).InletL;
	57
	58	SPECIFY
[187]	59	feed.F = 113.4 * 'kmol/h';
	60	feed.T = 291 * 'K';
	61	feed.P = 168.3 * 'kPa';
[1]	62	feed.z = [0.5, 0.5];
	63
[187]	64	zero.F = 0 * 'kmol/h';
	65	zero.T = 300 * 'K';
	66	zero.P = 1 * 'atm';
[1]	67	zero.z = [0.5, 0.5];
	68	zero.v = 0;
[187]	69	zero.h = 0 * 'J/mol';
[1]	70
[187]	71	cond.F = 68 * 'kmol/h';
	72	cond.P = 150 * 'kPa';
	73	cond.T = 281.75 * 'K';
[1]	74	cond.z = [0.6664, 0.3336];
	75
[187]	76	reb.F = 153 * 'kmol/h';
	77	reb.P = 185 * 'kPa';
	78	reb.T = 328.12 * 'K';
[1]	79	reb.z = [0.001848, 0.9982];
	80
	81	sec.trays.Emv = 1;
	82
	83	SET
	84	sec.NTrays = 2;
	85	#COLUMN
[187]	86	sec.trays.V = 4 * 'ft^3';
	87	sec.trays.Ah = 0.394 * 'ft^2';
	88	sec.trays.lw = 20.94 * 'in';
	89	sec.trays.hw = 0.125 * 'ft';
	90	sec.trays.Q = 0 * 'kW';
[1]	91	sec.trays.beta = 0.6;
	92	sec.trays.alfa = 4;
[187]	93	sec.trays.Ap = 3.94 * 'ft^2';
[1]	94
	95	INITIAL
[187]	96	sec.trays.OutletL.T = 290 *'K';
[1]	97	sec.trays.Level = 0.9 * sec.trays.hw;
	98	sec.trays.OutletL.z(1) = 0.5;
	99
	100	OPTIONS
[125]	101	TimeStep = 10;
	102	TimeEnd = 1000;
[111]	103
	104	# After running few seconds of transient the steady-state
	105	# can be obtained by using the results from that transient:
	106
[208]	107	#GuessFile="SectionColumn_Test_with2tray.rlt";
	108	#Dynamic = false;
[1]	109	end
	110
	111	# column section with 8 trays
	112	FlowSheet SectionColumn_Test_with8tray
	113	PARAMETERS
[213]	114	PP as Plugin(Brief="Physical Properties",
	115	Type="PP",
	116	Components = [ "isobutane", "benzene"],
	117	LiquidModel = "PR",
	118	VapourModel = "PR"
	119	);
[1]	120	NComp as Integer;
	121
	122	SET
	123	NComp = PP.NumberOfComponents;
	124
	125	DEVICES
	126	sec as Section_Column;
[125]	127	feed as liquid_stream;
	128	reb as vapour_stream;
	129	cond as liquid_stream;
[1]	130	zero as stream;
	131
	132	CONNECTIONS
	133	feed to sec.trays(5).Inlet;
	134
	135	zero to sec.trays([1:4]).Inlet;
	136	zero to sec.trays([6:8]).Inlet;
	137
	138	reb to sec.trays(8).InletV;
	139	cond to sec.trays(1).InletL;
	140
	141	SPECIFY
[187]	142	feed.F = 113.4 * 'kmol/h';
	143	feed.T = 291 * 'K';
	144	feed.P = 168.3 * 'kPa';
[1]	145	feed.z = [0.5, 0.5];
	146
[187]	147	zero.F = 0 * 'kmol/h';
	148	zero.T = 300 * 'K';
	149	zero.P = 1 * 'atm';
[1]	150	zero.z = [0.5, 0.5];
	151	zero.v = 0;
[187]	152	zero.h = 0 * 'J/mol';
[1]	153
[187]	154	cond.F = 68 * 'kmol/h';
	155	cond.P = 150 * 'kPa';
	156	cond.T = 281.75 * 'K';
[1]	157	cond.z = [0.6664, 0.3336];
	158
[187]	159	reb.F = 153 * 'kmol/h';
	160	reb.P = 185 * 'kPa';
	161	reb.T = 328.12 * 'K';
[1]	162	reb.z = [0.001848, 0.9982];
	163
	164	sec.trays.Emv = 1;
	165
	166	SET
	167	sec.NTrays = 8;
	168	#COLUMN
[187]	169	sec.trays.V = 4 * 'ft^3';
	170	sec.trays.Ah = 0.394 * 'ft^2';
	171	sec.trays.lw = 20.94 * 'in';
	172	sec.trays.hw = 0.125 * 'ft';
	173	sec.trays.Q = 0 * 'kW';
[1]	174	sec.trays.beta = 0.6;
	175	sec.trays.alfa = 4;
[187]	176	sec.trays.Ap = 3.94 * 'ft^2';
[1]	177
	178	INITIAL
[187]	179	sec.trays.OutletL.T = [290:(330-290)/(sec.NTrays-1):330] *'K';
[1]	180	sec.trays.Level = 0.3 * sec.trays.hw;
	181	sec.trays.OutletL.z(1) = 0.5;
	182
	183	OPTIONS
[125]	184	TimeStep = 1;
	185	TimeEnd = 100;
[208]	186	#GuessFile="SectionColumn_Test_with8tray.rlt";
	187	#Dynamic = false;
[1]	188	end
	189
	190
	191	FlowSheet Distillation_kettle_cond_Test
	192	PARAMETERS
[213]	193	PP as Plugin(Brief="Physical Properties",
	194	Type="PP",
[218]	195	Components = [ "isobutane", "n-pentane", "propylene",
	196	"benzene", "isobutene" ],
[213]	197	LiquidModel = "PR",
	198	VapourModel = "PR"
	199	);
[1]	200	NComp as Integer;
	201
	202	VARIABLES
	203	Qc as heat_rate (Brief="Heat rate removed from condenser");
	204	Qr as heat_rate (Brief="Heat rate supplied to reboiler");
	205
	206	SET
	207	NComp = PP.NumberOfComponents;
	208
	209	DEVICES
	210	col as Distillation_kettle_cond;
[125]	211	feed as source;
[1]	212	zero as stream;
	213
	214	CONNECTIONS
[187]	215	feed.Outlet to col.trays(5).Inlet;
[1]	216	zero to col.reb.Inlet;
	217	zero to col.trays([1:4]).Inlet;
	218	zero to col.trays([6:col.NTrays]).Inlet;
	219	Qc to col.cond.Q;
	220	Qr to col.reb.Q;
	221
	222	SPECIFY
[187]	223	feed.Outlet.F = 113.4 * 'kmol/h';
	224	feed.Outlet.T = 291 * 'K';
	225	feed.Outlet.P = 168.3 * 'kPa';
[125]	226	feed.Outlet.z = 1/NComp;
[1]	227
[187]	228	zero.F = 0 * 'kmol/h';
	229	zero.T = 300 * 'K';
	230	zero.P = 1 * 'atm';
[1]	231	zero.z = 1/NComp;
	232	zero.v = 0;
[187]	233	zero.h = 0 * 'J/mol';
[1]	234
[187]	235	col.sptop.Outlet2.F = 85 * 'kmol/h';
	236	col.reb.OutletL.F = 28.4 * 'kmol/h';
[19]	237	col.sptop.frac = 0.444445;
[187]	238	col.cond.OutletV.F = 0 * 'kmol/h';
	239	Qr = 3.7743e6 * 'kJ/h';
	240	Qc = -3.71e6 * 'kJ/h';
	241	col.pump1.dP = 16 * 'kPa';
[1]	242	col.trays.Emv = 1;
	243
	244	SET
	245	col.NTrays = 8;
[187]	246	col.cond.V = 2 * 'm^3';
	247	col.cond.Across = 1 * 'm^2';
	248	col.trays.V = 4 * 'ft^3';
	249	col.trays.Ah = 0.394 * 'ft^2';
	250	col.trays.lw = 20.94 * 'in';
	251	col.trays.hw = 0.125 * 'ft';
	252	col.trays.Q = 0 * 'kW';
[1]	253	col.trays.beta = 0.6;
	254	col.trays.alfa = 4;
[187]	255	col.trays.Ap = 3.94 * 'ft^2';
	256	col.reb.V = 2 * 'm^3';
	257	col.reb.Across = 1 * 'm^2';
[1]	258
	259	INITIAL
	260	# condenser
[187]	261	col.cond.OutletL.T = 260 *'K';
	262	col.cond.Level = 1 * 'm';
[1]	263	col.cond.OutletL.z([1:4]) = [0.65, 0.05, 0.01, 0.01];
	264
	265	# reboiler
[187]	266	col.reb.OutletL.T = 330 *'K';
	267	col.reb.Level = 1 * 'm';
[1]	268	col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
	269
	270	# column trays
[187]	271	col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
[19]	272	col.trays.Level = 1.2 * col.trays.hw;
[1]	273	col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];
	274
	275	OPTIONS
[187]	276	TimeStep = 0.1;
[125]	277	TimeEnd = 50;
	278	#time = [0:0.01:1, 2:50];
[208]	279	#GuessFile="Distillation_kettle_cond_Test.rlt";
	280	#Dynamic = false;
[1]	281	end
[111]	282
	283	FlowSheet Column_ctrl
	284	PARAMETERS
[213]	285	PP as Plugin(Brief="Physical Properties",
	286	Type="PP",
	287	Components = [ "isobutane", "n-pentane", "propylene",
	288	"benzene", "isobutene" ],
	289	LiquidModel = "PR",
	290	VapourModel = "PR"
	291	);
[111]	292	NComp as Integer;
	293
	294	Qcmin as heat_rate (Brief="Minimum Condenser Heat supplied");
	295	Qcmax as heat_rate (Brief="Maximum Condenser Heat supplied");
	296	Qrmin as heat_rate (Brief="Minimum Reboiler Heat supplied");
	297	Qrmax as heat_rate (Brief="Maximum Reboiler Heat supplied");
	298	Frmin as flow_mol (Brief="Minimum bottom flow rate");
	299	Frmax as flow_mol (Brief="Maximum bottom flow rate");
	300	Fcmin as flow_mol (Brief="Minimum reflux flow rate");
	301	Fcmax as flow_mol (Brief="Maximum reflux flow rate");
	302	Hmint as length (Brief="Minimum liquid level in top tank");
	303	Hmaxt as length (Brief="Maximum liquid level in top tank");
	304	Hminb as length (Brief="Minimum liquid level in reboiler");
	305	Hmaxb as length (Brief="Maximum liquid level in reboiler");
	306	Pmax as pressure (Brief="Maximum column pressure");
	307	Pmin as pressure (Brief="Minimum column pressure");
	308	Tmax as temperature (Brief="Maximum column temperature");
	309	Tmin as temperature (Brief="Minimum column temperature");
	310
	311	VARIABLES
	312	Qc as heat_rate (Brief="Heat rate removed from condenser");
	313	Qr as heat_rate (Brief="Heat rate supplied to reboiler");
	314	Had_top as Real (Brief="Dimensionless condenser level");
	315	Had_bot as Real (Brief="Dimensionless reboiler level");
	316	Pad as Real (Brief="Dimensionless pressure");
	317	Tad as Real (Brief="Dimensionless temperature");
[112]	318	RR as positive (Brief="Reflux ratio");
	319
[111]	320	SET
	321	NComp = PP.NumberOfComponents;
	322
	323	DEVICES
	324	col as Distillation_kettle_cond;
[125]	325	feed as source;
[111]	326	zero as stream;
	327	TCcond as PIDIncr_Ideal_AW;
	328	LCtop as PIDIncr_Ideal_AW;
	329	LCbot as PIDIncr_Ideal_AW;
	330	PC as PIDIncr_Ideal_AW;
	331
	332	CONNECTIONS
[134]	333	feed.Outlet to col.trays(5).Inlet;
[111]	334	zero to col.reb.Inlet;
	335	zero to col.trays([1:4]).Inlet;
	336	zero to col.trays([6:col.NTrays]).Inlet;
	337	Qc to col.cond.Q;
	338	Qr to col.reb.Q;
	339
	340	EQUATIONS
	341	"Temperature Controller"
[177]	342	TCcond.Parameters.tau = 0*'s';
	343	TCcond.Parameters.tauSet = 0*'s';
[111]	344	TCcond.Parameters.alpha = 0.3;
	345	TCcond.Parameters.bias = 0.5;
	346	TCcond.Parameters.gamma = 1;
	347	TCcond.Parameters.beta = 1;
	348	TCcond.Options.action = 1;
	349	TCcond.Options.clip = 1;
	350	TCcond.Options.autoMan = 0;
[177]	351	TCcond.Parameters.intTime = 60*'s';
[111]	352	TCcond.Parameters.gain = 0.6;
[177]	353	TCcond.Parameters.derivTime = 1*'s';
	354	TCcond.Ports.setPoint = ((15+273.15) * 'K' - Tmin)/(Tmax-Tmin);
[111]	355	TCcond.Ports.input = Tad;
	356	Tad = (col.cond.OutletL.T-Tmin)/(Tmax-Tmin);
	357	Qc = Qcmin+(Qcmax-Qcmin)*TCcond.Ports.output;
	358
	359	"Pressure Controller"
[177]	360	PC.Parameters.tau = 0*'s';
	361	PC.Parameters.tauSet = 0*'s';
[111]	362	PC.Parameters.alpha = 0.3;
	363	PC.Parameters.bias = 0;
	364	PC.Parameters.gamma = 1;
	365	PC.Parameters.beta = 1;
	366	PC.Options.action = -1;
	367	PC.Options.clip = 1;
	368	PC.Options.autoMan = 0;
[177]	369	PC.Parameters.intTime = 50*'s';
[111]	370	PC.Parameters.gain = 0.5;
[177]	371	PC.Parameters.derivTime = 1*'s';
	372	PC.Ports.setPoint = (2.0*'bar'-Pmin)/(Pmax-Pmin);
[111]	373	PC.Ports.input = Pad;
	374	Pad = (col.cond.OutletV.P-Pmin)/(Pmax-Pmin);
	375	col.cond.OutletV.F = (Fcmin+(Fcmax-Fcmin)*PC.Ports.output);
	376
	377	"Ttop Level Controller"
[177]	378	LCtop.Parameters.tau = 0*'s';
	379	LCtop.Parameters.tauSet = 0*'s';
[111]	380	LCtop.Parameters.alpha = 0.3;
	381	LCtop.Parameters.bias = 0.5;
	382	LCtop.Parameters.gamma = 1;
	383	LCtop.Parameters.beta = 1;
	384	LCtop.Options.action = -1;
	385	LCtop.Options.clip = 1;
	386	LCtop.Options.autoMan = 0;
[177]	387	LCtop.Parameters.intTime = 10*'s';
[111]	388	LCtop.Parameters.gain = 1;
[218]	389	LCtop.Parameters.derivTime = 0*'s';
[177]	390	LCtop.Ports.setPoint = (1.0 * 'm' - Hmint)/(Hmaxt-Hmint);
[111]	391	LCtop.Ports.input = Had_top;
	392	Had_top = (col.cond.Level-Hmint)/(Hmaxt-Hmint);
	393	col.sptop.Outlet1.F = Fcmin + (Fcmax-Fcmin) * LCtop.Ports.output;
	394
	395	"Tbottom Level Controller"
[177]	396	LCbot.Parameters.tau = 0*'s';
	397	LCbot.Parameters.tauSet = 0*'s';
[111]	398	LCbot.Parameters.alpha = 0.3;
	399	LCbot.Parameters.bias = 0.5;
	400	LCbot.Parameters.gamma = 1;
	401	LCbot.Parameters.beta = 1;
	402	LCbot.Options.action = -1;
	403	LCbot.Options.clip = 1;
	404	LCbot.Options.autoMan = 0;
[177]	405	LCbot.Parameters.intTime = 100*'s';
[111]	406	LCbot.Parameters.gain = 1;
[218]	407	LCbot.Parameters.derivTime = 0*'s';
[177]	408	LCbot.Ports.setPoint = (1.0 * 'm' - Hminb)/(Hmaxb-Hminb);
[111]	409	LCbot.Ports.input = Had_bot;
	410	Had_bot = (col.reb.Level-Hminb)/(Hmaxb-Hminb);
	411	col.reb.OutletL.F = Frmin + (Frmax-Frmin) * LCbot.Ports.output;
	412
[112]	413	RR * (col.cond.OutletV.F + col.sptop.Outlet1.F) = col.sptop.Outlet2.F;
	414
[177]	415	if time < 1 * 'h' then
[237]	416	col.sptop.Outlet2.F = 75 * 'kmol/h'; # reflux
[112]	417	else
[177]	418	col.sptop.Outlet2.F = 85 * 'kmol/h'; # reflux
[112]	419	end
	420
[111]	421	SPECIFY
[177]	422	feed.Outlet.F = 113.4 * 'kmol/h';
	423	feed.Outlet.T = 291 * 'K';
	424	feed.Outlet.P = 168.3 * 'kPa';
[125]	425	feed.Outlet.z = 1/NComp;
[111]	426
[177]	427	zero.F = 0 * 'kmol/h';
	428	zero.T = 300 * 'K';
	429	zero.P = 1 * 'atm';
[111]	430	zero.z = 1/NComp;
	431	zero.v = 0;
[177]	432	zero.h = 0 * 'J/mol';
[111]	433
[177]	434	Qr = 3e6 * 'kJ/h';
	435	col.pump1.dP = 16 * 'kPa';
[111]	436	col.trays.Emv = 1;
	437
	438	SET
	439	col.NTrays = 8;
[177]	440	col.cond.V = 2 * 'm^3';
	441	col.cond.Across = 1 * 'm^2';
	442	col.trays.V = 4 * 'ft^3';
	443	col.trays.Ah = 0.394 * 'ft^2';
	444	col.trays.lw = 20.94 * 'in';
	445	col.trays.hw = 0.125 * 'ft';
	446	col.trays.Q = 0 * 'kW';
[111]	447	col.trays.beta = 0.6;
	448	col.trays.alfa = 4;
[177]	449	col.trays.Ap = 3.94 * 'ft^2';
	450	col.reb.V = 2 * 'm^3';
	451	col.reb.Across = 1 * 'm^2';
[111]	452
[177]	453	Qrmax = 5e6 * 'kJ/h';
	454	Qrmin = 1e6 * 'kJ/h';
	455	Frmin = 0 * 'kmol/h';
	456	Frmax = 60 * 'kmol/h';
	457	Fcmin = 0 * 'kmol/h';
	458	Fcmax = 120 * 'kmol/h';
	459	Hmint = 0 * 'm';
	460	Hmaxt = 2 * 'm';
	461	Hminb = 0 * 'm';
	462	Hmaxb = 2 * 'm';
	463	Pmin = 0.5 * 'bar';
	464	Pmax = 4 * 'bar';
	465	Qcmax = -5e5 * 'kJ/h';
	466	Qcmin = -5e6 * 'kJ/h';
	467	Tmax = (30+273.15) * 'K';
	468	Tmin = (-20+273.15) * 'K';
[111]	469
	470	INITIAL
	471	# condenser
[177]	472	col.cond.OutletL.T = 260 *'K';
	473	col.cond.Level = 1 * 'm';
[111]	474	col.cond.OutletL.z([1:4]) = [0.2, 0.2, 0.4, 0.05];
	475
	476	# reboiler
[177]	477	col.reb.OutletL.T = 350 *'K';
	478	col.reb.Level = 1 * 'm';
[111]	479	col.reb.OutletL.z([1:4]) = [0.1, 0.4, 0.1, 0.3];
[218]	480
[111]	481	# column trays
[177]	482	col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
[111]	483	col.trays.Level = 1.2 * col.trays.hw;
	484	col.trays.OutletL.z([1:4]) = [0.15, 0.3, 0.25, 0.2];
	485
	486	OPTIONS
[218]	487	TimeStep = 0.1;
	488	TimeEnd = 5;
[177]	489	TimeUnit = 'h';
[218]	490	InitialFile = "Column_ctrl.rlt";
[208]	491	#GuessFile = "Column_ctrl.rlt";
	492	#Dynamic = false;
[111]	493	end
	494

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format