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 MUA BÀI THÍ NGHIỆM Lớp cử một người (Lớp trưởng) lên Bộ môn Quá trình thiết bị Tại C4-109 Thứ năm ngày 16 tháng 2 năm 2012 Sáng 9h-11h30 Chiều 13h30-15h30 Gặp cô Hoa mua tài liệu. 25.000/5 bài. Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Tại sao phải sử dụng nội suy trong tính toán các quá trình CN Hóa học??? Các đường cong này được xây dựng từ??? Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Tính toán? Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Các thuật toán nội suy: Tuyến tính, Lagrance, Newton, … x s y s y 1 x 1 x n x k-1 x k y k-1 y k y n 0 y x Nhưng… Không có số liệu thực nghiệm!!! Giả thiết đường cong nối giữa hai điểm là đường thẳng Đó là??? Nội suy tuyến tính!!! Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Nội suy tuyến tính x s y s y 1 x 1 x n x k-1 x k y k-1 y k y n 0 y x Phương trình đường thẳng đi qua hai điểm (x k-1 ,y k-1 ) và (x k ,y k ): 1 1 1 1 − − − − − − = − − kk k kk k yy yy xx xx )( 1 1 1 1 − − − − − − − += ks kk kk ks xx xx yy yy Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Nội suy tuyến tính Thuật toán: 1. Chỉ ra khoảng (x k-1 ,x k ) chứa giá trị x s  giá trị của k 2. Đưa giá trị của k tìm được vào biểu thức nội suy tuyến tính )( 1 1 1 1 − − − − − − − += ks kk kk ks xx xx yy yy Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Nội suy tuyến tính Procedure NOSUY(xs:real;VAR ys:real;Y,X:mX); Begin k:=0; {so sánh xs với các giá trị x1, …, xn} {Số vòng lặp sẽ là không xác định!!!?} {Sử dụng cấu trúc:} {Repeat … Until hoặc While … End} Repeat k:=k+1; Until xs < x[k]; Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Nội suy tuyến tính Procedure NOSUY(xs:real;VAR ys:real;Y,X:mX); Begin k:=0; Repeat k:=k+1; Until xs < x[k]; {Ra khỏi vòng lặp trên đã tìm được giá trị k} ys:=y[k-1]+(y[k]-y[k-1])*(xs-x[k-1]) /(x[k]-x[k-1]); End; Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Ví dụ 1: Cho hỗn hợp lỏng Bezne – Toluen, biết hàm lượng Benzen trong pha lỏng x = 0,4 (phần mol). Hãy xác định hàm lượng Benze trong pha hơi ở trạng thái cân bằng. x = 0,4 y CB = ? [...]... chính} Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Ví dụ 2: {Chương trình chính} BEGIN clrscr; writeln (‘Nhập số điểm thực nghiệm n = ’); readln (n); {Nhập các số liệu của pha lỏng x[i]} For i:=1 to n do Begin writeln (‘x[‘,i,’] =‘);readln (x[i]); End; Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Ví dụ 2: {Chương trình chính} BEGIN … {Nhập các số liệu của pha hơi ở TTCB y[i]} For... đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Ví dụ 3: Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Ví dụ 3: {Chương trình chính} BEGIN … {Tìm áp suất hơi bão hòa tại ts=119.6} writeln (‘Nhập giá trị ts =’);readln(ts); NOISUY (ts,ps,P,T); {Hiển thị kết quả} writeln (‘Áp suất hơi bão hòa, pS =’,ps); readln; END Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Ví dụ 4: Cho hệ hơi... giản: Vấn đề nội suy trong kỹ thuật hóa học Ví dụ 1: Procedure NOSUY(xs:real;VAR ys:real;Y,X:mX); Begin k:=0; Repeat k:=k+1; Until xs . 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 MUA BÀI THÍ NGHIỆM Lớp. Hoa mua tài liệu. 25.000/5 bài. Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Tại sao phải sử dụng nội suy trong tính toán các quá trình CN Hóa học? ?? Các đường cong này được. xây dựng từ??? Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Tính toán? Ứng dụng đơn giản: Vấn đề nội suy trong kỹ thuật hóa học Các thuật toán nội suy: Tuyến tính, Lagrance,