Lecture -6.pdf

Modelling and Simulation in
aspect of Natural Gas
Dr Sourav Poddar
Department of Chemical Engineering
National Institute of Technology
Warangal, TS

Distillation column
% MATLAB Program for solving Distillation Program
A=[0.07 0.18 0.15 0.24
0.04 0.24 0.10 0.65
0.54 0.42 0.54 0.10
0.35 0.16 0.21 0.01]
f = [0.15*70 0.25*70 0.40*70 0.2*70]
disp('Solutions for D1 B1 D2 B2 is:')
X=Af'
D1=X(1)
B1=X(2)
D2=X(3)
B2=X(4)
disp {'Solve for column 2'}
D=D1+B1%43.75 MOL/MIN
XDx=(0.07*D1+0.18*B1)/D %0.114 mole fraction
XDs=(0.04*D1+0.24*B1)/D %0.120 mole fraction
XDt=(0.54*D1+0.42*B1)/D %0.492 mole fraction
XDb=(0.35*D1+0.16*B1)/D %0.274 mole fraction
%The mole fraction for column 3 are solved for directly
by evaluating
%Equation
D2=X(3)
B2=X(4)
disp('Solve for colum 3')
B=D2+B2 %26.25 mol/min
XBx=(0.15*D2+0.24*B2)/B %0.2100 mole fraction
XBs=(0.10*D2+0.65*B2)/B %0.4667 mole fraction
XBt=(0.54*D2+0.10*B2)/B %0.2467 mole fraction
XBb=(0.21*D2+0.01*B2)/B %0.0767 mole fraction

Result
dist
0.0700 0.1800 0.1500 0.2400
0.0400 0.2400 0.1000 0.6500
0.5400 0.4200 0.5400 0.1000
0.3500 0.1600 0.2100 0.0100
A =
f = 10.5000 17.5000 28.0000 14.0000
Solutions for D1 B1 D2 B2 is:
X = 26.2500, 17.5000, 8.7500, 17.5000
D1 = 26.2500; B1 = 17.5000; D2 = 8.7500; B2 = 17.5000
{'Solve for column 2'}
D = 43.7500; XDx = 0.1140; XDs = 0.1200; XDt = 0.4920; XDb = 0.2740; D2 = 8.7500
B2 = 17.5000
Solve for colum 3
B = 26.2500; XBx = 0.2100; XBs = 0.4667; XBt = 0.2467; XBb = 0.0767

HEAT EXCHANGER
% Figure(1) refers to Case 1
% Click "Run and Advance" to generate the Matlab 3D co-ordinates for all 3
% cases.
% Click "Run" wouldl only generate Figure(1).
% The utility source is taken to be in the middle of the plant.
% The values [0.5 0.5 0.5] and [0.9290 0.6940 0.1250] are RGB colours for
% grey (heat exchangers) and yellow-orange (other equipment) as Matlab does
% not have a pre-built value for them.
% The vectors are sometimes given a +0.2 or -0.2 to indicate the
% approximate suggested coordinates for the new heaters/heat exchangers to
% be installed (+0.2 or -0.2 from a nearby existing equipment in line with
% the stream).
%HE33E061
x1 = 2.9400;
y1 = 2.9620;
z1 = 0.1560;
%HE31E01
x2 = 4.2300;
y2 = 0.1500;
z2 = 1.1900;
%HE31E02
x3 = 5.3000;
y3 = 0.1600;
z3 = 0.4600;
%HE31E03 (G - heater)
x4 = 4.7800;
y4 = 1.4332;
z4 = 2.6200;
%HE31E04
x5 = 5.5300;
y5 = 0.3800;
z5 = 1.4450;

% HE31E05 (C - cooler)
x6 = 5.5690;
y6 = 0.4300;
z6 = 0.6500;
% HE32E01 (H - heater)
x7 = 4.4800;
y7 = 1.5250;
z7 = 1.0700;
% HE32E02 (I - 32E02)
x8 = 2.0143;
y8 = 1.6650;
z8 = 1.0620;
% HE32E03 (EN - heat exchanger)
x9 = 1.3868;
y9 = 1.6000;
z9 = 1.7600;
% HE32E04 (DN - heat exchanger)
x10 = 3.7444;
y10 = 1.2000;
z10 = 1.0400;
% HE32E05
x11 = 2.7800;
y11 = 2.0375;
z11 = 1.0900;
% HE32E06 (Q - heater)
x12 = 2.0275;
y12 = 1.1350;
z12 = 1.0620;
% HE32E07 (D,M - heat exchanger)
x13 = 1.3868;
y13 = 1.2000;
z13 = 1.5300;
% HE32E08
x14 = 6.6900;
y14 = 3.0200;
z14 = 0.0900;
% HE33E01
x15 = 6.2800;
y15 = 0.6000;
z15 = 1.3500;
% HE33E02 (FO - heat exchanger)
x16 = 6.2800;
y16 = 0.6000;
z16 = 4.6000;
% HE33E03
x17 = 5.7800;
y17 = 0.4500;
z17 = 4.6000;

% HE33E04 (J - heater)
x18 = 3.2800;
y18 = 0.2117;
z18 = 1.4700;
% HE33E05 (AP - heat exchanger)
x19 = 2.9800;
y19 = 0.2160;
z19 = 1.5900;
% HE33E06 (A - heater)
x20 = 2.6800;
y20 = 0.2336;
z20 = 1.4900;
% HE33E07 (BG - heat exchanger)
x21 = 2.9800;
y21 = 0.2250;
z21 = 0.0900;
% HE33E08 (B - cooler)
x22 = 3.2800;
y22 = 0.1860;
z22 = 0.1200;
% HE33E09 (K - heater)
x23 = 5.3755;
y23 = 0.2000;
z23 = 3.0100;
% HE33E10 (CL - heat exchanger)
x24 = 5.3700;
y24 = 1.6500;
z24 = 1.0450;
% HE33E11
x25 = 4.7550;
y25 = 1.2000;
z25 = 1.0450;
% HE33E12
x26 = 5.2800;
y26 = 1.6120;
z26 = 2.6500;
% HE33E13
x27 = 5.6400;
y27 = 1.6212;
z27 = 2.5950;
%HE33E14
x28 = 5.6400;
y28 = 1.6212;
z28 = 2.6950;
%HE33E15
x29 = 5.4300;
y29 = 1.6120;
z29 = 2.6000;
%Utility source
x30 = 4.2500;
y30 = 2.2500;
z30 = 2.7500;
sz = 230;
figure(1)

scatter3(x1, y1, z1, sz, [0.9290 0.6940 0.1250], 'filled',
'MarkerEdgeColor','k'); hold on;
scatter3(x2, y2, z2, sz, 'b', 'filled', 'MarkerEdgeColor','k')
'MarkerEdgeColor','k')
scatter3(x4, y4, z4, sz, 'b', 'filled', 'MarkerEdgeColor','k');
scatter3(x6, y6, z6, sz, 'r', 'filled', 'MarkerEdgeColor','k');
scatter3(x7, y7, z7, sz, 'r', 'filled', 'MarkerEdgeColor','k')
scatter3(x9, y9, z9, sz, [0.5,0.5,0.5], 'filled',
'MarkerEdgeColor','k');

hold off
grid on
xlabel('x-coordinates')
ylabel ('y-coordinates')
zlabel ('Floor Level')
figure(2)

scatter3(x13, y13, z13, sz, [0.5,0.5,0.5], 'filled', 'MarkerEdgeColor','k')
scatter3(x23, y23, z23, sz, 'r', 'filled',
scatter3(x25, y25, z25, sz, [0.9290 0.6940 0.1250],
'filled', 'MarkerEdgeColor','k')
scatter3(x26, y26, z26, sz, [0.9290 0.6940 0.1250],
scatter3(x27, y27, z27, sz, [0.9290 0.6940 0.1250],
scatter3(x28, y28, z28, sz, [0.9290 0.6940 0.1250],
scatter3(x29, y29, z29, sz, [0.9290 0.6940 0.1250],
scatter3(x30, y30, z30,[300],'r', 'd', 'filled',
%1. HE for K and A (HE33E05 & HE33E09)

v1_1 = [x23, y23-0.2, z23];
v1_2 = [x23, y23-0.2, z19];
v1_3 = [x23, y19-0.2, z19];
v2_1 = [x23, y23-0.2, z19];
v2_2 = [x23, y19-0.2, z19];
v2_3 = [x19, y19-0.2, z19];
v1=[v1_1;v1_2;v1_3];
v2=[v2_1;v2_2;v2_3];
plot3(v1(:,1),v1(:,2),v1(:,3),'g', 'Linewidth', 2)
%2. HE for H01 and K (HE33E09)
v3_1 = [x30, y30, z30];
v3_2 = [x30, y30, z23];
v3_3 = [x30, y23+0.2, z23];
v4_1 = [x30, y30, z23];
v4_2 = [x30, y23+0.2, z23];
v4_3 = [x23, y23+0.2, z23];
v3=[v3_1;v3_2;v3_3];
v4=[v4_1;v4_2;v4_3];
plot3(v3(:,1),v3(:,2),v3(:,3),'r', 'Linewidth', 2)
%3. HE for H02 and L (HE33E10)
v5_1 = [x30, y30, z30];
v5_2 = [x30, y30, z24];
v5_3 = [x30, y24+0.2, z24];
v6_1 = [x30, y30, z24];
v6_2 = [x30, y24+0.2, z24];
v6_3 = [x24, y24+0.2, z24];
v5=[v5_1;v5_2;v5_3];
v6=[v6_1;v6_2;v6_3];
%4. HE for H03 and L (HE33E05)
v7_1 = [x30, y30, z30];
v7_2 = [x30, y30, z24];
v7_3 = [x30, y24-0.2, z24];
v8_1 = [x30, y30, z24];
v8_2 = [x30, y24-0.2, z24];
v8_3 = [x24, y24-0.2, z24];
v7=[v7_1;v7_2;v7_3];
v8=[v8_1;v8_2;v8_3];

hold off
grid on
figure(3)

v1_1 = [x23, y23-0.2, z23];
v1_2 = [x23, y23-0.2, z19];
v1_3 = [x23, y19-0.2, z19];
v2_1 = [x23, y23-0.2, z19];
v2_2 = [x23, y19-0.2, z19];
v2_3 = [x19, y19-0.2, z19];
v1=[v1_1;v1_2;v1_3];
v2=[v2_1;v2_2;v2_3];
hold off
grid on

Results
Figure -1
Figure – 2

Hetran - (to simulate an heat exchanger), The B-JAC input file name that contains the heat exchanger specification., The Aerotran
block is the interface to the B-JAC Aerotran program for designing and simulating air-cooled heat exchangers

Lecture -6.pdf

More Related Content

Similar to Lecture -6.pdf (20)

More from Sourav Poddar (20)

Recently uploaded (20)

Lecture -6.pdf