Context Navigation

sample_column.mso @ 580

Last change on this file since 580 was 580, checked in by gerson bicca, 15 years ago
still more updates
Property svn:eol-style set to `native` Property svn:keywords set to `Id`
File size: 21.6 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 580 2008-07-26 19:34:58Z bicca $
	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.P = 185 * 'kPa';
	77	reb.T = 328.12 * 'K';
[1]	78	reb.z = [0.001848, 0.9982];
	79
	80	sec.trays.Emv = 1;
[398]	81	sec.trays(1).OutletV.F = 150 * 'kmol/h';
[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
[530]	96	sec.trays.OutletL.T = 270 *'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
[398]	159	# reb.F = 153 * 'kmol/h';
[187]	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;
[398]	165	sec.trays(1).OutletV.F = 150 * 'kmol/h';
[1]	166
	167	SET
	168	sec.NTrays = 8;
	169	#COLUMN
[187]	170	sec.trays.V = 4 * 'ft^3';
	171	sec.trays.Ah = 0.394 * 'ft^2';
	172	sec.trays.lw = 20.94 * 'in';
	173	sec.trays.hw = 0.125 * 'ft';
	174	sec.trays.Q = 0 * 'kW';
[1]	175	sec.trays.beta = 0.6;
	176	sec.trays.alfa = 4;
[187]	177	sec.trays.Ap = 3.94 * 'ft^2';
[1]	178
	179	INITIAL
[530]	180	sec.trays.OutletL.T = [290:(300-290)/(sec.NTrays-1):300] *'K';
[1]	181	sec.trays.Level = 0.3 * sec.trays.hw;
	182	sec.trays.OutletL.z(1) = 0.5;
	183
	184	OPTIONS
[125]	185	TimeStep = 1;
	186	TimeEnd = 100;
[208]	187	#GuessFile="SectionColumn_Test_with8tray.rlt";
	188	#Dynamic = false;
[1]	189	end
	190
	191
	192	FlowSheet Distillation_kettle_cond_Test
	193	PARAMETERS
[213]	194	PP as Plugin(Brief="Physical Properties",
	195	Type="PP",
[218]	196	Components = [ "isobutane", "n-pentane", "propylene",
	197	"benzene", "isobutene" ],
[213]	198	LiquidModel = "PR",
	199	VapourModel = "PR"
	200	);
[1]	201	NComp as Integer;
	202
	203	VARIABLES
[310]	204	Qc as energy_source (Brief="Heat rate removed from condenser");
	205	Qr as energy_source (Brief="Heat rate supplied to reboiler");
[1]	206
	207	SET
	208	NComp = PP.NumberOfComponents;
	209
	210	DEVICES
	211	col as Distillation_kettle_cond;
[125]	212	feed as source;
[1]	213	zero as stream;
	214
	215	CONNECTIONS
[187]	216	feed.Outlet to col.trays(5).Inlet;
[1]	217	zero to col.reb.Inlet;
	218	zero to col.trays([1:4]).Inlet;
	219	zero to col.trays([6:col.NTrays]).Inlet;
[313]	220	Qc.OutletQ to col.cond.InletQ;
	221	Qr.OutletQ to col.reb.InletQ;
[1]	222
	223	SPECIFY
[580]	224	feed.F = 113.4 * 'kmol/h';
	225	feed.T = 291 * 'K';
	226	feed.P = 168.3 * 'kPa';
	227	feed.Composition = 1/NComp;
[1]	228
[187]	229	zero.F = 0 * 'kmol/h';
	230	zero.T = 300 * 'K';
	231	zero.P = 1 * 'atm';
[1]	232	zero.z = 1/NComp;
	233	zero.v = 0;
[187]	234	zero.h = 0 * 'J/mol';
[1]	235
[187]	236	col.sptop.Outlet2.F = 85 * 'kmol/h';
	237	col.reb.OutletL.F = 28.4 * 'kmol/h';
[19]	238	col.sptop.frac = 0.444445;
[187]	239	col.cond.OutletV.F = 0 * 'kmol/h';
[555]	240	Qr.OutletQ = 3.7743e6 * 'kJ/h';
	241	Qc.OutletQ = -3.71e6 * 'kJ/h';
[187]	242	col.pump1.dP = 16 * 'kPa';
[1]	243	col.trays.Emv = 1;
	244
[398]	245	col.alfaTopo = 2;
	246
[1]	247	SET
	248	col.NTrays = 8;
[187]	249	col.cond.V = 2 * 'm^3';
	250	col.cond.Across = 1 * 'm^2';
	251	col.trays.V = 4 * 'ft^3';
	252	col.trays.Ah = 0.394 * 'ft^2';
	253	col.trays.lw = 20.94 * 'in';
	254	col.trays.hw = 0.125 * 'ft';
	255	col.trays.Q = 0 * 'kW';
[1]	256	col.trays.beta = 0.6;
	257	col.trays.alfa = 4;
[187]	258	col.trays.Ap = 3.94 * 'ft^2';
	259	col.reb.V = 2 * 'm^3';
	260	col.reb.Across = 1 * 'm^2';
[1]	261
	262	INITIAL
	263	# condenser
[187]	264	col.cond.OutletL.T = 260 *'K';
	265	col.cond.Level = 1 * 'm';
[1]	266	col.cond.OutletL.z([1:4]) = [0.65, 0.05, 0.01, 0.01];
	267
	268	# reboiler
[530]	269	col.reb.OutletL.T = 350 *'K';
[187]	270	col.reb.Level = 1 * 'm';
[1]	271	col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
	272
	273	# column trays
[187]	274	col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
[19]	275	col.trays.Level = 1.2 * col.trays.hw;
[1]	276	col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];
	277
	278	OPTIONS
[187]	279	TimeStep = 0.1;
[125]	280	TimeEnd = 50;
[208]	281	#GuessFile="Distillation_kettle_cond_Test.rlt";
	282	#Dynamic = false;
[1]	283	end
[111]	284
[270]	285
[111]	286	FlowSheet Column_ctrl
	287	PARAMETERS
[213]	288	PP as Plugin(Brief="Physical Properties",
	289	Type="PP",
	290	Components = [ "isobutane", "n-pentane", "propylene",
	291	"benzene", "isobutene" ],
	292	LiquidModel = "PR",
	293	VapourModel = "PR"
	294	);
[111]	295	NComp as Integer;
	296
	297	Qcmin as heat_rate (Brief="Minimum Condenser Heat supplied");
	298	Qcmax as heat_rate (Brief="Maximum Condenser Heat supplied");
	299	Qrmin as heat_rate (Brief="Minimum Reboiler Heat supplied");
	300	Qrmax as heat_rate (Brief="Maximum Reboiler Heat supplied");
	301	Frmin as flow_mol (Brief="Minimum bottom flow rate");
	302	Frmax as flow_mol (Brief="Maximum bottom flow rate");
	303	Fcmin as flow_mol (Brief="Minimum reflux flow rate");
	304	Fcmax as flow_mol (Brief="Maximum reflux flow rate");
	305	Hmint as length (Brief="Minimum liquid level in top tank");
	306	Hmaxt as length (Brief="Maximum liquid level in top tank");
	307	Hminb as length (Brief="Minimum liquid level in reboiler");
	308	Hmaxb as length (Brief="Maximum liquid level in reboiler");
	309	Pmax as pressure (Brief="Maximum column pressure");
	310	Pmin as pressure (Brief="Minimum column pressure");
	311	Tmax as temperature (Brief="Maximum column temperature");
	312	Tmin as temperature (Brief="Minimum column temperature");
	313
	314	VARIABLES
[310]	315	Qc as energy_source (Brief="Heat rate removed from condenser");
	316	Qr as energy_source (Brief="Heat rate supplied to reboiler");
[111]	317	Had_top as Real (Brief="Dimensionless condenser level");
	318	Had_bot as Real (Brief="Dimensionless reboiler level");
	319	Pad as Real (Brief="Dimensionless pressure");
	320	Tad as Real (Brief="Dimensionless temperature");
[112]	321	RR as positive (Brief="Reflux ratio");
	322
[111]	323	SET
	324	NComp = PP.NumberOfComponents;
	325
	326	DEVICES
	327	col as Distillation_kettle_cond;
[125]	328	feed as source;
[111]	329	zero as stream;
[297]	330	TCcond as PIDIncr;
	331	LCtop as PIDIncr;
	332	LCbot as PIDIncr;
	333	PC as PIDIncr;
[111]	334
	335	CONNECTIONS
[134]	336	feed.Outlet to col.trays(5).Inlet;
[111]	337	zero to col.reb.Inlet;
	338	zero to col.trays([1:4]).Inlet;
	339	zero to col.trays([6:col.NTrays]).Inlet;
[313]	340	Qc.OutletQ to col.cond.InletQ;
	341	Qr.OutletQ to col.reb.InletQ;
[111]	342
[574]	343	SET
	344	TCcond.tau = 0*'s';
	345	TCcond.tauSet = 0*'s';
	346	TCcond.alpha = 0.3;
	347	TCcond.bias = 0.5;
	348	TCcond.gamma = 1;
	349	TCcond.beta = 1;
[555]	350	TCcond.Action = "Direct";
	351	TCcond.Clip = "Clipped";
	352	TCcond.Mode = "Automatic";
[574]	353	TCcond.intTime = 60*'s';
	354	TCcond.gain = 0.6;
	355	TCcond.derivTime = 1*'s';
	356
	357	PC.tau = 0*'s';
	358	PC.tauSet = 0*'s';
	359	PC.alpha = 0.3;
	360	PC.bias = 0;
	361	PC.gamma = 1;
	362	PC.beta = 1;
	363	PC.Action = "Reverse";
	364	PC.Clip = "Clipped";
	365	PC.Mode = "Automatic";
	366	PC.intTime = 50*'s';
	367	PC.gain = 0.5;
	368	PC.derivTime = 1*'s';
	369
	370	LCtop.tau = 0*'s';
	371	LCtop.tauSet = 0*'s';
	372	LCtop.alpha = 0.3;
	373	LCtop.bias = 0.5;
	374	LCtop.gamma = 1;
	375	LCtop.beta = 1;
	376	LCtop.Action = "Reverse";
	377	LCtop.Clip = "Clipped";
	378	LCtop.Mode = "Automatic";
	379	LCtop.intTime = 10*'s';
	380	LCtop.gain = 1;
	381	LCtop.derivTime = 0*'s';
	382
	383	LCbot.tau = 0*'s';
	384	LCbot.tauSet = 0*'s';
	385	LCbot.alpha = 0.3;
	386	LCbot.bias = 0.5;
	387	LCbot.gamma = 1;
	388	LCbot.beta = 1;
	389	LCbot.Action = "Reverse";
	390	LCbot.Clip = "Clipped";
	391	LCbot.Mode = "Automatic";
	392	LCbot.intTime = 100*'s';
	393	LCbot.gain = 1;
	394	LCbot.derivTime = 0*'s';
	395
	396	EQUATIONS
[177]	397	TCcond.Ports.setPoint = ((15+273.15) * 'K' - Tmin)/(Tmax-Tmin);
[111]	398	TCcond.Ports.input = Tad;
	399	Tad = (col.cond.OutletL.T-Tmin)/(Tmax-Tmin);
[555]	400	Qc.OutletQ = Qcmin+(Qcmax-Qcmin)*TCcond.Ports.output;
[111]	401
[346]	402	PC.Ports.setPoint = (4.0*'bar'-Pmin)/(Pmax-Pmin);
[111]	403	PC.Ports.input = Pad;
	404	Pad = (col.cond.OutletV.P-Pmin)/(Pmax-Pmin);
	405	col.cond.OutletV.F = (Fcmin+(Fcmax-Fcmin)*PC.Ports.output);
	406
[177]	407	LCtop.Ports.setPoint = (1.0 * 'm' - Hmint)/(Hmaxt-Hmint);
[111]	408	LCtop.Ports.input = Had_top;
	409	Had_top = (col.cond.Level-Hmint)/(Hmaxt-Hmint);
	410	col.sptop.Outlet1.F = Fcmin + (Fcmax-Fcmin) * LCtop.Ports.output;
	411
[177]	412	LCbot.Ports.setPoint = (1.0 * 'm' - Hminb)/(Hmaxb-Hminb);
[111]	413	LCbot.Ports.input = Had_bot;
	414	Had_bot = (col.reb.Level-Hminb)/(Hmaxb-Hminb);
	415	col.reb.OutletL.F = Frmin + (Frmax-Frmin) * LCbot.Ports.output;
	416
[112]	417	RR * (col.cond.OutletV.F + col.sptop.Outlet1.F) = col.sptop.Outlet2.F;
	418
[177]	419	if time < 1 * 'h' then
[237]	420	col.sptop.Outlet2.F = 75 * 'kmol/h'; # reflux
[112]	421	else
[177]	422	col.sptop.Outlet2.F = 85 * 'kmol/h'; # reflux
[112]	423	end
	424
[111]	425	SPECIFY
[580]	426	feed.F = 113.4 * 'kmol/h';
	427	feed.T = 291 * 'K';
	428	feed.P = 5 * 'bar';
	429	feed.Composition = 1/NComp;
[111]	430
[177]	431	zero.F = 0 * 'kmol/h';
	432	zero.T = 300 * 'K';
	433	zero.P = 1 * 'atm';
[111]	434	zero.z = 1/NComp;
	435	zero.v = 0;
[177]	436	zero.h = 0 * 'J/mol';
[111]	437
[555]	438	Qr.OutletQ = 4e6 * 'kJ/h';
[177]	439	col.pump1.dP = 16 * 'kPa';
[111]	440	col.trays.Emv = 1;
[398]	441	col.alfaTopo = 2;
[111]	442
	443	SET
	444	col.NTrays = 8;
[177]	445	col.cond.V = 2 * 'm^3';
	446	col.cond.Across = 1 * 'm^2';
	447	col.trays.V = 4 * 'ft^3';
	448	col.trays.Ah = 0.394 * 'ft^2';
	449	col.trays.lw = 20.94 * 'in';
	450	col.trays.hw = 0.125 * 'ft';
	451	col.trays.Q = 0 * 'kW';
[111]	452	col.trays.beta = 0.6;
	453	col.trays.alfa = 4;
[177]	454	col.trays.Ap = 3.94 * 'ft^2';
	455	col.reb.V = 2 * 'm^3';
	456	col.reb.Across = 1 * 'm^2';
[111]	457
[297]	458	# Controllers type
	459	TCcond.PID_Select = "Ideal_AW";
	460	PC.PID_Select = "Ideal_AW";
	461	LCtop.PID_Select = "Ideal_AW";
	462	LCbot.PID_Select = "Ideal_AW";
	463
[177]	464	Qrmax = 5e6 * 'kJ/h';
	465	Qrmin = 1e6 * 'kJ/h';
	466	Frmin = 0 * 'kmol/h';
	467	Frmax = 60 * 'kmol/h';
	468	Fcmin = 0 * 'kmol/h';
	469	Fcmax = 120 * 'kmol/h';
	470	Hmint = 0 * 'm';
	471	Hmaxt = 2 * 'm';
	472	Hminb = 0 * 'm';
	473	Hmaxb = 2 * 'm';
	474	Pmin = 0.5 * 'bar';
[346]	475	Pmax = 6 * 'bar';
[177]	476	Qcmax = -5e5 * 'kJ/h';
	477	Qcmin = -5e6 * 'kJ/h';
	478	Tmax = (30+273.15) * 'K';
	479	Tmin = (-20+273.15) * 'K';
[111]	480
	481	INITIAL
	482	# condenser
[177]	483	col.cond.OutletL.T = 260 *'K';
	484	col.cond.Level = 1 * 'm';
[111]	485	col.cond.OutletL.z([1:4]) = [0.2, 0.2, 0.4, 0.05];
	486
	487	# reboiler
[177]	488	col.reb.OutletL.T = 350 *'K';
	489	col.reb.Level = 1 * 'm';
[111]	490	col.reb.OutletL.z([1:4]) = [0.1, 0.4, 0.1, 0.3];
[218]	491
[111]	492	# column trays
[177]	493	col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
[111]	494	col.trays.Level = 1.2 * col.trays.hw;
	495	col.trays.OutletL.z([1:4]) = [0.15, 0.3, 0.25, 0.2];
	496
	497	OPTIONS
[218]	498	TimeStep = 0.1;
	499	TimeEnd = 5;
[177]	500	TimeUnit = 'h';
[218]	501	InitialFile = "Column_ctrl.rlt";
[530]	502	DAESolver(File="dassl");
[208]	503	#GuessFile = "Column_ctrl.rlt";
	504	#Dynamic = false;
[111]	505	end
	506
[498]	507	# packed column section with 8 trays
	508	FlowSheet PackedSectionColumn
	509	PARAMETERS
	510	PP as Plugin(Brief="Physical Properties",
	511	Type="PP",
	512	Components = [ "isobutane", "n-pentane", "propylene",
	513	"benzene", "isobutene" ],
	514	LiquidModel = "PR",
	515	VapourModel = "PR"
	516	);
	517	NComp as Integer;
[270]	518
[498]	519	SET
	520	NComp = PP.NumberOfComponents;
	521
	522	DEVICES
	523	sec as Packed_Section_Column;
	524	feed as liquid_stream;
	525	reb as vapour_stream;
	526	cond as liquid_stream;
	527	zero as stream;
	528
	529	CONNECTIONS
	530	feed to sec.stage(5).Inlet;
	531	zero to sec.stage([1:4]).Inlet;
	532	zero to sec.stage([6:sec.NStages]).Inlet;
	533	reb to sec.stage(8).InletV;
	534	cond to sec.stage(1).InletL;
	535
	536	SPECIFY
	537	feed.F = 113.4 * 'kmol/h';
	538	feed.T = 291 * 'K';
	539	feed.P = 1.66 * 'atm';
	540	feed.z = [0.2, 0.2, 0.2, 0.2, 0.2];
	541
	542	zero.F = 0 * 'kmol/h';
	543	zero.T = 300 * 'K';
	544	zero.P = 1 * 'atm';
	545	zero.z = [0.2, 0.2, 0.2, 0.2, 0.2];
	546	zero.v = 0;
	547	zero.h = 0 * 'J/mol';
	548
	549	cond.F = 85 * 'kmol/h';
	550	cond.P = 2.33 * 'atm';
	551	cond.T = 283.5 * 'K';
	552	cond.z = [0.599, 0.044, 0.035, 0.007, 0.315];
	553
	554	reb.F = 137.57 * 'kmol/h';
	555	reb.P = 2.46 * 'atm';
	556	reb.T = 325 * 'K';
	557	reb.z = [0.16, 0.542, 0.013, 0.008, 0.277];
	558
[515]	559	sec.dP = 0.008 * 'atm'; #queda de pressao na secao!
[498]	560
	561	SET
[515]	562	sec.H = 4 * 'm';#altura do recheio
	563	sec.NStages = 8; #numero de estagios teoricos do recheio
[498]	564	sec.stage.Q = 0 * 'kW';
	565	sec.stage.d = 1.009 * 'm';
	566	sec.stage.Cpo = 0.763;
	567	sec.stage.e = 0.951;
	568	sec.stage.a = 112.6 * 'm^2/m^3';
[515]	569	sec.stage.Qsil = 0.1;#coeficiente de resistencia do recheio
[498]	570
	571	INITIAL
	572	sec.stage.OutletL.T =[283:(325-283)/(sec.NStages-1):325] *'K';
	573	sec.stage.ML = 0.01 * 'kmol';
	574	sec.stage.OutletL.z([1:4]) = [0.3, 0.2, 0.002, 0.1];
	575
	576	OPTIONS
	577	DAESolver(File="dassl");
	578	NLASolver(File="nlasolver");
	579	TimeStep = 10;
	580	TimeEnd = 100;
	581	end
	582
	583	FlowSheet Packed_kettle_cond_Test
	584	PARAMETERS
	585	PP as Plugin(Brief="Physical Properties",
	586	Type="PP",
	587	Components = [ "isobutane", "n-pentane", "propylene",
	588	"benzene", "isobutene" ],
	589	LiquidModel = "PR",
	590	VapourModel = "PR"
	591	);
	592	NComp as Integer;
	593
	594	VARIABLES
	595	Qc as energy_source (Brief="Heat rate removed from condenser");
	596	Qr as energy_source (Brief="Heat rate supplied to reboiler");
	597
	598	SET
	599	NComp = PP.NumberOfComponents;
	600
	601	DEVICES
	602	col as PackedDistillation_kettle_cond;
	603	feed as source;
	604	zero as stream;
	605
	606	CONNECTIONS
	607	feed.Outlet to col.stage(5).Inlet;
	608	zero to col.reb.Inlet;
	609	zero to col.stage([1:4]).Inlet;
	610	zero to col.stage([6:col.NStages]).Inlet;
	611	Qc.OutletQ to col.cond.InletQ;
	612	Qr.OutletQ to col.reb.InletQ;
	613
	614	SPECIFY
[580]	615	feed.F = 113.4 * 'kmol/h';
	616	feed.T = 291 * 'K';
	617	feed.P = 168.3 * 'kPa';
	618	feed.Composition = 1/NComp;
[498]	619
	620	zero.F = 0 * 'kmol/h';
	621	zero.T = 300 * 'K';
	622	zero.P = 1 * 'atm';
	623	zero.z = 1/NComp;
	624	zero.v = 0;
	625	zero.h = 0 * 'J/mol';
	626
	627	col.sptop.Outlet2.F = 85 * 'kmol/h';
	628	col.reb.OutletL.F = 28.4 * 'kmol/h';
	629	col.sptop.frac = 0.444445;
	630	col.cond.OutletV.F = 0 * 'kmol/h';
[555]	631	Qr.OutletQ = 3.7743e6 * 'kJ/h';
	632	Qc.OutletQ = -3.71e6 * 'kJ/h';
[498]	633	col.pump1.dP = 16 * 'kPa';
	634
[515]	635	col.dP = 0.024 * 'atm'; #queda de pressao entre o topo e o fundo da coluna.
[498]	636
	637	SET
[515]	638	col.H = 3.5 * 'm'; # altura de recheio
	639	col.NStages = 8; # numero de estagios teoricos do recheio
[498]	640	col.stage.Q = 0 * 'kW';
	641	col.stage.d = 2.24 * 'ft';
	642	col.stage.Cpo = 0.763;
	643	col.stage.e = 0.951;
	644	col.stage.a = 112.6 * 'm^2/m^3';
	645
	646	col.cond.V = 2 * 'm^3';
	647	col.cond.Across = 1 * 'm^2';
	648	col.reb.V = 2 * 'm^3';
	649	col.reb.Across = 1 * 'm^2';
[515]	650	col.stage.Qsil = 2; # coeficiente de resistencia do recheio
[498]	651
	652	INITIAL
	653	# condenser
	654	col.cond.OutletL.T = 260 *'K';
	655	col.cond.Level = 1 * 'm';
	656	col.cond.OutletL.z([1:4]) = [0.65, 0.05, 0.01, 0.01];
	657
	658	# reboiler
	659	col.reb.OutletL.T = 330 *'K';
	660	col.reb.Level = 1 * 'm';
	661	col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
	662
	663	# column trays
	664	col.stage.OutletL.T = [290:(330-290)/(col.NStages-1):330] * 'K';
	665	col.stage.ML = 0.1 * 'kmol';
	666	col.stage.OutletL.z([1:4]) = [0.15, 0.5, 0.001, 0.1];
	667
	668	OPTIONS
	669	TimeStep = 10;
	670	TimeEnd = 500;
	671	end
	672
	673	FlowSheet PackedColumn_ctrl
	674	PARAMETERS
	675	PP as Plugin(Brief="Physical Properties",
	676	Type="PP",
	677	Components = [ "isobutane", "n-pentane", "propylene",
	678	"benzene", "isobutene" ],
	679	LiquidModel = "PR",
	680	VapourModel = "PR"
	681	);
	682	NComp as Integer;
	683
	684	Qcmin as heat_rate (Brief="Minimum Condenser Heat supplied");
	685	Qcmax as heat_rate (Brief="Maximum Condenser Heat supplied");
	686	Qrmin as heat_rate (Brief="Minimum Reboiler Heat supplied");
	687	Qrmax as heat_rate (Brief="Maximum Reboiler Heat supplied");
	688	Frmin as flow_mol (Brief="Minimum bottom flow rate");
	689	Frmax as flow_mol (Brief="Maximum bottom flow rate");
	690	Fcmin as flow_mol (Brief="Minimum reflux flow rate");
	691	Fcmax as flow_mol (Brief="Maximum reflux flow rate");
	692	Hmint as length (Brief="Minimum liquid level in top tank");
	693	Hmaxt as length (Brief="Maximum liquid level in top tank");
	694	Hminb as length (Brief="Minimum liquid level in reboiler");
	695	Hmaxb as length (Brief="Maximum liquid level in reboiler");
	696	Pmax as pressure (Brief="Maximum column pressure");
	697	Pmin as pressure (Brief="Minimum column pressure");
	698	Tmax as temperature (Brief="Maximum column temperature");
	699	Tmin as temperature (Brief="Minimum column temperature");
	700
	701	VARIABLES
	702	Qc as energy_source (Brief="Heat rate removed from condenser");
	703	Qr as energy_source (Brief="Heat rate supplied to reboiler");
	704	Had_top as Real (Brief="Dimensionless condenser level");
	705	Had_bot as Real (Brief="Dimensionless reboiler level");
	706	Pad as Real (Brief="Dimensionless pressure");
	707	Tad as Real (Brief="Dimensionless temperature");
	708	RR as positive (Brief="Reflux ratio");
	709
	710	SET
	711	NComp = PP.NumberOfComponents;
	712
	713	DEVICES
	714	col as PackedDistillation_kettle_cond;
	715	feed as source;
	716	zero as stream;
	717	TCcond as PIDIncr;
	718	LCtop as PIDIncr;
	719	LCbot as PIDIncr;
	720	PC as PIDIncr;
	721
	722	CONNECTIONS
	723	feed.Outlet to col.stage(5).Inlet;
	724	zero to col.reb.Inlet;
	725	zero to col.stage([1:4]).Inlet;
	726	zero to col.stage([6:col.NStages]).Inlet;
	727	Qc.OutletQ to col.cond.InletQ;
	728	Qr.OutletQ to col.reb.InletQ;
[555]	729
	730	SET
	731	TCcond.Action = "Direct";
	732	TCcond.Clip = "Clipped";
	733	TCcond.Mode = "Automatic";
[574]	734	TCcond.tau = 0*'s';
	735	TCcond.tauSet = 0*'s';
	736	TCcond.alpha = 0.3;
	737	TCcond.bias = 0.5;
	738	TCcond.gamma = 1;
	739	TCcond.beta = 1;
	740	TCcond.intTime = 60*'s';
	741	TCcond.gain = 0.6;
	742	TCcond.derivTime = 1*'s';
[498]	743
[555]	744	PC.Action = "Reverse";
	745	PC.Clip = "Clipped";
	746	PC.Mode = "Automatic";
[574]	747	PC.tau = 0*'s';
	748	PC.tauSet = 0*'s';
	749	PC.alpha = 0.3;
	750	PC.bias = 0;
	751	PC.gamma = 1;
	752	PC.beta = 1;
	753	PC.intTime = 50*'s';
	754	PC.gain = 0.5;
	755	PC.derivTime = 1*'s';
[555]	756
	757	LCtop.Action = "Reverse";
	758	LCtop.Clip = "Clipped";
	759	LCtop.Mode = "Automatic";
[574]	760	LCtop.tau = 0*'s';
	761	LCtop.tauSet = 0*'s';
	762	LCtop.alpha = 0.3;
	763	LCtop.bias = 0.5;
	764	LCtop.gamma = 1;
	765	LCtop.beta = 1;
	766	LCtop.intTime = 10*'s';
	767	LCtop.gain = 1;
	768	LCtop.derivTime = 0*'s';
	769
[555]	770	LCbot.Action = "Reverse";
	771	LCbot.Clip = "Clipped";
	772	LCbot.Mode = "Automatic";
[574]	773	LCbot.tau = 0*'s';
	774	LCbot.tauSet = 0*'s';
	775	LCbot.alpha = 0.3;
	776	LCbot.bias = 0.5;
	777	LCbot.gamma = 1;
	778	LCbot.beta = 1;
	779	LCbot.intTime = 100*'s';
	780	LCbot.gain = 1;
	781	LCbot.derivTime = 0*'s';
[555]	782
[498]	783	EQUATIONS
	784	TCcond.Ports.setPoint = ((15+273.15) * 'K' - Tmin)/(Tmax-Tmin);
	785	TCcond.Ports.input = Tad;
	786	Tad = (col.cond.OutletL.T-Tmin)/(Tmax-Tmin);
[555]	787	Qc.OutletQ = Qcmin+(Qcmax-Qcmin)*TCcond.Ports.output;
[498]	788
	789	PC.Ports.setPoint = (4.0*'bar'-Pmin)/(Pmax-Pmin);
	790	PC.Ports.input = Pad;
	791	Pad = (col.cond.OutletV.P-Pmin)/(Pmax-Pmin);
	792	col.cond.OutletV.F = (Fcmin+(Fcmax-Fcmin)*PC.Ports.output);
	793
	794	LCtop.Ports.setPoint = (1.0 * 'm' - Hmint)/(Hmaxt-Hmint);
	795	LCtop.Ports.input = Had_top;
	796	Had_top = (col.cond.Level-Hmint)/(Hmaxt-Hmint);
	797	col.sptop.Outlet1.F = Fcmin + (Fcmax-Fcmin) * LCtop.Ports.output;
	798
	799	LCbot.Ports.setPoint = (1.0 * 'm' - Hminb)/(Hmaxb-Hminb);
	800	LCbot.Ports.input = Had_bot;
	801	Had_bot = (col.reb.Level-Hminb)/(Hmaxb-Hminb);
	802	col.reb.OutletL.F = Frmin + (Frmax-Frmin) * LCbot.Ports.output;
	803
	804	RR * (col.cond.OutletV.F + col.sptop.Outlet1.F) = col.sptop.Outlet2.F;
	805
	806	if time < 1 * 'h' then
	807	col.sptop.Outlet2.F = 75 * 'kmol/h'; # reflux
	808	else
	809	col.sptop.Outlet2.F = 85 * 'kmol/h'; # reflux
	810	end
	811
	812	SPECIFY
[580]	813	feed.F = 113.4 * 'kmol/h';
	814	feed.T = 291 * 'K';
	815	feed.P = 5 * 'bar';
	816	feed.Composition = 1/NComp;
[498]	817
	818	zero.F = 0 * 'kmol/h';
	819	zero.T = 300 * 'K';
	820	zero.P = 1 * 'atm';
	821	zero.z = 1/NComp;
	822	zero.v = 0;
	823	zero.h = 0 * 'J/mol';
	824
[555]	825	Qr.OutletQ = 4e6 * 'kJ/h';
[498]	826	col.pump1.dP = 16 * 'kPa';
	827
[515]	828	col.dP = 0.024 * 'atm'; #queda de pressao entre o topo e o fundo da coluna
[498]	829
	830	SET
[515]	831	col.H = 3.5 * 'm'; #altura do recheio
	832	col.NStages = 8; #numero de estagios teoricos do recheio
[498]	833	col.stage.Q = 0 * 'kW';
	834	col.stage.d = 2.24 * 'ft';
	835	col.stage.Cpo = 0.763;
	836	col.stage.e = 0.951;
	837	col.stage.a = 112.6 * 'm^2/m^3';
	838
	839	col.cond.V = 2 * 'm^3';
	840	col.cond.Across = 1 * 'm^2';
	841	col.reb.V = 2 * 'm^3';
	842	col.reb.Across = 1 * 'm^2';
	843
[515]	844	col.stage.Qsil = 2; #coeficiente de resistencia do recheio
[498]	845
	846	# Controllers type
	847	TCcond.PID_Select = "Ideal_AW";
	848	PC.PID_Select = "Ideal_AW";
	849	LCtop.PID_Select = "Ideal_AW";
	850	LCbot.PID_Select = "Ideal_AW";
	851
	852	Qrmax = 5e6 * 'kJ/h';
	853	Qrmin = 1e6 * 'kJ/h';
	854	Frmin = 0 * 'kmol/h';
	855	Frmax = 60 * 'kmol/h';
	856	Fcmin = 0 * 'kmol/h';
	857	Fcmax = 120 * 'kmol/h';
	858	Hmint = 0 * 'm';
	859	Hmaxt = 2 * 'm';
	860	Hminb = 0 * 'm';
	861	Hmaxb = 2 * 'm';
	862	Pmin = 0.5 * 'bar';
	863	Pmax = 6 * 'bar';
	864	Qcmax = -5e5 * 'kJ/h';
	865	Qcmin = -5e6 * 'kJ/h';
	866	Tmax = (30+273.15) * 'K';
	867	Tmin = (-20+273.15) * 'K';
	868
	869	INITIAL
	870	# condenser
	871	col.cond.OutletL.T = 260 *'K';
	872	col.cond.Level = 1 * 'm';
	873	col.cond.OutletL.z([1:4]) = [0.2, 0.2, 0.4, 0.05];
	874
	875	# reboiler
	876	col.reb.OutletL.T = 350 *'K';
	877	col.reb.Level = 1 * 'm';
	878	col.reb.OutletL.z([1:4]) = [0.1, 0.4, 0.1, 0.3];
	879
	880	# column trays
	881	col.stage.OutletL.T = [290:(330-290)/(col.NStages-1):330] * 'K';
	882	col.stage.ML = 0.1 * 'kmol';
	883	col.stage.OutletL.z([1:4]) = [0.15, 0.3, 0.25, 0.2];
	884
	885	OPTIONS
	886	TimeStep = 0.1;
	887	TimeEnd = 5;
	888	TimeUnit = 'h';
[499]	889	InitialFile = "Packed_Column_ctrl.rlt";
	890	#GuessFile = "Packed_Column_ctrl.rlt";
	891	#Dynamic = false;
[498]	892	end

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

Download in other formats:

Original Format