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 CuuDuongThanCong com https //fb com/tailieudientucn[.]
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
CuuDuongThanCong.com https://fb.com/tailieudientucntt
Trang 2MUA 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.
Trang 3Ứ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ừ???
CuuDuongThanCong.com https://fb.com/tailieudientucntt
Trang 4Ứng dụng đơn giản: Vấn đề nội suy trong
kỹ thuật hóa học
Tính toán?
Trang 5Ứ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 k-1 x k x n
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!!!
CuuDuongThanCong.com https://fb.com/tailieudientucntt
Trang 6Ứ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 k-1 x k x n
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
k k
k
k k
k
y y
y y
x
x
x
x
)
1
1
k k
k k
k
x x
y
y y
y
Trang 7Ứ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
k k
k k
k
x x
y
y y
y
CuuDuongThanCong.com https://fb.com/tailieudientucntt
Trang 8Ứ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];
Trang 9Ứ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;
CuuDuongThanCong.com https://fb.com/tailieudientucntt
Trang 10Ứ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 = ?