PHƢƠNG PHÁP SỐ TRONG CÔNG NGHỆ HÓA HỌC Mã học phần CH3454 TS Nguyễn Đặng Bình Thành BM Máy & TBCN Hóa chất Numerical Methods in Chemical Engineering Tuần 3 CuuDuongThanCong com https //fb com/tailieud[.]
Trang 1PHƯƠNG PHÁP SỐ TRONG CÔNG NGHỆ HÓA HỌC
Mã học phần: CH3454
TS Nguyễn Đặng Bình Thành
BM:Máy & TBCN Hóa chất
Numerical Methods in Chemical Engineering
Trang 2Ứng dụng nội suy trong tính toán số đĩa
lý thuyết và chiều cao tháp chƣng luyện
Chương trình
Program chungluyen_1;
uses crt;
type
mX = array [1 50] of real;
var
X,Y:mX;
F,P,W,xF,xP,xW,R:real;
NLT,NLTC,NLTL,n,i,k:integer;
{Các chương trình con}
Trang 3Program chungluyen_1;
…
{Các chương trình con}
Procedure NOSUY(xs:real;VAR ys:real;Y,X:mX); Begin
k:=0;
Repeat k:=k+1;
Until xs<x[k];
ys:=y[k-1]+ (y[k]-y[k-1])*(xs-x[k-1]) /(x[k]-x[k-1]);
End;
Trang 4Ứng dụng nội suy trong tính toán số đĩa
lý thuyết và chiều cao tháp chƣng luyện
Chương trình
Program chungluyen_1;
…
{Các chương trình con}
Function DLVL(xs:real):real;
Begin
DLVL:=R/(R+1)*xs+xP/(R+1);
End;
Function DLVC(xs:real):real;
Begin
DLVC:=(R+F/P)*xs/(R+1)+(1-F/P)/(R+1)*xW;
Trang 5Program chungluyen_1;
…
{Chương trình chính}
BEGIN
clrscr;
{số liệu đầu}
write (‘Nhập số điểm thực nghiệm n = ’);
readln(n);
{Nhập các giá trị thực nghiệm}
Trang 6Ứng dụng nội suy trong tính toán số đĩa
lý thuyết và chiều cao tháp chƣng luyện
Chương trình
Program chungluyen_1;
…
BEGIN
…
{Nhập các giá trị thực nghiệm}
For i:=1 to n do
Begin write (‘X[’,i,‘] = ’);readln(x[i]);
End;
Trang 7Program chungluyen_1;
…
BEGIN
…
{Nhập các giá trị thực nghiệm}
For i:=1 to n do
Begin write (‘Y[’,i,‘] = ’);readln(Y[i]);
End;
Trang 8Ứng dụng nội suy trong tính toán số đĩa
lý thuyết và chiều cao tháp chƣng luyện
Chương trình
Program chungluyen_1;
…
BEGIN
…
{Nhập các số liệu yêu cầu}
write (‘Nhập F = ’);readln(F);
write (‘Nhập xF = ’);readln(xF);
write (‘Nhập xP = ’);readln(xP);
write (‘Nhập xW = ’);readln(xF);
Trang 9Program chungluyen_1;
…
BEGIN
…
{Tính toán số đĩa lý thuyết đoạn luyện}
ys:=xP;
NLTL:=0;
Repeat
NOISUY(ys,xs,X,Y);
ys:=DLVL(xs);
NLTL:=NLTL+1;
Until xs <= xF;
tháp:
F = P + W F.x P:=F*(xP-xF)/(xP-xW);F = P.xP + W.xW W:=F-P;
Trang 10Ứng dụng nội suy trong tính toán số đĩa
lý thuyết và chiều cao tháp chƣng luyện
Chương trình
Program chungluyen_1;
…
BEGIN
…
{Tính toán số đĩa lý thuyết đoạn chưng}
yf:=DLVC(xf);
ys:=yf;
NLTC:=0;
Repeat
NOISUY(ys,xs,X,Y);
ys:=DLVC(xs);
NLTC:=NLTC+1;