II. MƠ PHỎNG QUÁ TRÌNH CHUYỂN HĨA BIODIESEL TRÊN CƠ SỞ NGƠN NGỮ
2.Mơ hình 1 bình phản ứng chuyển hĩa gián đoạn
Thuật tốn mơ phỏng
Hướng giải của thuật tốn là sau khi giải ra được giá trị cho thời điểm t(n) bất kỳ, sẽ lấy giá trị đĩ làm điều kiện đầu để giải hệ phương trình vi phân tìm ra giá trị cho thời điểm tiếp theo t(n+1) . Quá trình giải cứ tiếp tục cho đến khi thỏa điều kiện về thời gian chuyển hĩa tndungđã được đề cập.
Các bước lập trình tính tốn tạo thành đoạn code trọng yếu khi mơ phỏng cho một thiết bịđược gọi là đoạn code TB1. Nếu cĩ 2 hay 3 thiết bị nối tiếp nhau sau code TB1 sẽ bổ sung các đoạn code tương tự cĩ kí hiệu TB2, TB3, ... Cách này được dùng để lập trình cho hệ thống gồm 2 hay 3 thiết bị phản ứng dạng khuấy lí tưởng hoạt động liên tục.
Trang 53 Sau khi thỏa điều kiện tính tốn đặt ra (vịng lặp kết thúc), các kết quả cần thiết như thời gian chuyển hĩa In, khối lượng ester tạo thành,... được truy xuất dưới dạng bảng dữ liệu và các đồ thị.
Trong kết quả tính tốn mơ phỏng, sử dụng giả thiết sau đây:
- Pha ester gồm: methyl ester và tồn bộ TG, (2/3 DG+1/3MG) chưa chuyển hĩa.
- Hàm lượng methyl ester trong pha ester = (khối lượng methyl ester / khối
lượng của cả pha ester) x 100.
- Chỉ số hiệu suất sản phẩm = (khối lượng E tạo thành / khối lượng mỡ cá) x 100. - Pha glyceine gồm tất cả các chất cịn lại (G, Me, (1/3 DG+2/3MG), KOH, FFA,
W,S).
- Phần trăm glycerine tổng trong pha ester = (khối lượng G được qui ra từ khối
lượng (TG+DG+MG) / khối lượng của cả pha ester) x 100.
- Phần trăm glycerine tổng trong pha glycerine = [(khối lượng G tự do+khối lượng G qui ra từ khối lượng MG+DG) / khối lượng của cả pha] x 100.
Trang 54
Hình 10 - Sơ đồ thuật tốn mơ phỏng cho quá trình chuyển hĩa gián đoạn
Đ Bắt đầu tn=tstep xC, xv X(1,:)=xC; x0=xC u2(1); Fr(1); Ue(1) i=length([0:tstep:tn]) tn <= tndung
Ue(i); u2(i); Fr(i)
tspan=[0 tstep] [T1,V1]=ode45(@heptvp,tspan,x0) X(i,:)=V1(end,:) x0=V1(end,:) tn=tn+tstep i=length([0:tstep:tn]) Kết quả:X, T, u2, mE,%E .... In=tn-0.5 Kết thúc S TB1
Trang 55