Đoạn chương trình xác định mô hình tháp C-02 với đầu ra là nhiệt độ tháp
NT=30; NF=10; M0(1)=86.36; i=2:NT-1; M0(i)=8.1*ones(1,NT-2); M0(NT)=138.23; taul=0.063; F0=16.2; qF0 = 1; L0=19.575; L0b=L0 + qF0*F0; lambda=0; V0=32.625;V0t=V0+(1-qF0)*F0 x=X(1:NT)'; M=X(NT+1:2*NT)';
LT = U(1); VB = U(2); D = U(3); B = U(4); F = U(5); zF = U(6); qF = U(7);
i=1:NT-1;
V(i)=VB*ones(1,NT-1); i=NF:NT-1;
V(i)=V(i) + (1-qF)*F; i=2:NF;
L(i) = L0b + (M(i)-M0(i))./taul + lambda.*(V(i-1)-V0); i=NF+1:NT-1;
L(i) = L0 + (M(i)-M0(i))./taul + lambda.*(V(i-1)-V0t); i=2:NF-1;
dMdt(i) = L(i+1)- L(i)+ V(i-1)- V(i);
dMxdt(i)= (V(i-1).*Sv(i-1).*x(i-1) - V(i).*Sv(i).*x(i) + L(i+1).*Sl(i+1).*x(i+1) - L(i).*Sl(i).*x(i))./Sl(i); i=NF+1:NT-1;
dMdt(i) = L(i+1)- L(i) + V(i-1) - V(i);
dMxdt(i)= (V(i-1).*Sv(i-1).*x(i-1) - V(i).*Sv(i).*x(i) + L(i+1).*Sl(i+1).*x(i+1) - L(i).*Sl(i).*x(i))./Sl(i); dMdt(NF) = dMdt(NF) + F; dMdt(1) = L(2)- V(1) - B; dMxdt(1)= (V(1).*(Sl(1)-Sv(1)).*x(1) + L(2).*(Sl(2).*x(2) - Sl(1).*x(1)))./Sl(1); dMdt(NT) = V(NT-1) - LT - D; dMxdt(NT)= ((Sl(NT-1).*x(NT-1) - Sl(NT).*x(NT)).*V(NT-1) + L(NT).*Sl(NT).*x(NT))./Sl(NT);
dxdt(i) = (dMxdt(i) - x(i).*dMdt(i) )./M(i); xprime=[dxdt';dMdt'];
2.3.4. Khảo sát động học của tháp C02
Các chế độ làm việc của tháp chưng cất, động học dòng lỏng, ảnh hưởng của dòng trong và dòng ngoài cũng như ảnh hưởng của nhiễu nguồn cấp tới nồng độ sản phẩm đỉnh, đáy tháp, nhiệt độ tháp được khảo sát và phân tích một cách sâu sắc dựa trên các mô hình phi tuyến đã xây dựng từ phần trên.
Mô phỏng tháp với các số liệu vận hành
Sơ đồ mô phỏng tháp C-02 được cho trên hình 2.15 với các số liệu vận hành của lưu lượng hồi lưu, lưu lượng hơi cấp nhiệt, lưu lượng sản phẩm đỉnh và đáy. Nồng độ sản phẩm đỉnh và đáy tháp được xác định khi mô phỏng tháp với mô hình phi tuyến. Đáp ứng của sản phẩm đỉnh tháp yD với các số liệu vận hành được cho trên hình 2.16.
Mô hình tháp chưng với cấu hình LV, kết quả mô phỏng tháp với các số liệu
vận hành cho thấy: sản phẩm đỉnh đạt giá trị xác lập là khoảng 0.993 sau thời gian xấp xỉ 150 phút Hình 2.15. Sơ đồ mô phỏng tháp C02 Thời gian [phút] y D [ p h ần mo l]