ĐẠI HỌC QUOC GIA ĐẠI HỌC BÁCH KHOA TP.HCM KHOA CƠ KHÍ BỘ MÔN CÔNG NGHỆ NHIỆT LẠNH a KHONG MINH TRUONG PHUN BANG PHAN MEM ANSYS FLUENT Chuyén Nganh: CONG NGHE NHIET Mã Số : 605280 LUẬN VĂN THẠC SĨ Tp Hồ Chí Minh, tháng 12 năm 2016 CÔNG TRÌNH ĐƯỢC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BACH KHOA - DAI HOC QUOC GIA TP.HCM Cán bộ hướng dẫn khoa học: TS. HÀ ANH TUNG (Ghi rõ họ, tên, học hàm, học vi và chữ ký) Cán bộ cham nhận xét 1: TS. NGUYEN VĂN TUYẾN (Ghi rõ họ, tên, học hàm, học vi và chữ ký) Cán bộ chấm nhận xét 2: PGS.TS DANG THANH TRUNG (Ghi rõ họ, tên, học hàm, học vi va chữ ky) Luận văn thạc sĩ được bảo vệ tại Trường Dai Học Bách Khoa — ĐHQG Tp.HCM, ngày 05 tháng 01 năm 2017. Thanh phan hội đồng đánh giá luận văn thạc sĩ bao gồm: 1.
Ủy viên — Phản biện 1: 4. Ủy viên — Phản biện 2: 5. Ủy viên: Xác nhận của Chủ tịch Hội đồng đánh giá luận văn và Trưởng Khoa quản lý chuyên ngành sau khi luận văn đã được sửa chữa (nếu có). CHỦ TỊCH HỘI ĐÔNG TRƯỞNG KHOA CƠ KHÍ (Họ tên và chữ ký) (Họ tên và chữ ký) ĐẠI HỌC QUỐC GIA TP.HCM CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM ĐẠI HỌC BÁCH KHOA Độc Lập — Tự Do — Hạnh Phúc Tp.
HCM, ngày 05 tháng 12 năm 2016 NHIEM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: KHÔNG MINH TRƯỚNG MSHV: 12054829 Ngày, tháng, năm sinh: 27/07/1986 Nơi sinh: Hải Phòng Chuyên nghành: CÔNG NGHỆ NHIỆT Mã số: 605280 I. TÊN ĐÈ TÀI TÍNH TOÁN LÝ THUYÉT, THỰC NGHIỆM VÀ MÔ PHỎNG QUÁ TRÌNH SAY DUNG DICH SỮA TƯƠI TRONG THIET BỊ SAY PHUN BANG PHAN MEM ANSYS FLUENT II. NHIEM VU VÀ NOI DUNG - Nghiên cứu tong quan về công nghệ say phun, phương pháp tính toán lý thuyết quá trình say và thiết kế thiết bi say phun băng phương pháp truyền thống. - Thực nghiệm sấy mẫu trên thiết bi sấy phun thực tế để xác định các thông số cần thiết cho quá trình tính kiểm tra dựa theo phương pháp tính toán truyền thống.
- _ Nghiên cứu tổng quan phương số trong thủy động lực học và các mô hính toán cho quá trình say phun. - Ung dụng phần mềm ANSYS FLUENT để mô phỏng quá trình sấy. So sánh, đánh giá kết quả thu được từ mô hình mô phỏng và thực nghiệm. - - Kết luận và kiến nghị.
NGÀY GIAO NHIỆM VU: 15/08/2016 IV. NGÀY HOÀN THÀNH NHIỆM VU: 05/12/2016 V. HO VA TÊN CÁN BO HUONG DAN: TS. HÀ ANH TUNG CAN BỘ HƯỚNG DAN CHỦ NHIỆM BỘ MÔN (Họ tên và chữ ký) ĐÀO TẠO (Họ tên và chữ ký) TS.
HÀ ANH TÙNG GS. LE CHÍ HIỆP TRƯỞNG KHOA CƠ KHÍ (Họ tên và chữ ký) LOI CAM ON Hoc viên xin bay to long biét on sâu sac dén: > Cán bộ hướng dẫn: TS HA ANH TÙNG đã tận tinh hướng dan, chi bảo và giúp đỡ học viên hoàn thành các nội dung nghiên cứu trong luận văn này. Xin chân thành cảm ơn quý Thay Cô trong bộ môn Công nghệ nhiệt lạnh — Truong Dai học Bách Khoa Tp.HCM đã tận tâm truyền đạt kiến thức cũng như kinh nghiệm trong suôt quá trình học viên học tập và nghiên cứu tại trường. Cam ơn gia đình đã tạo điều kiện thuận lợi về thời gian, vật chat, chia sẻ, động viên hoc viên trong suôt thời gian qua.
Các bạn bè thân thiết, đồng nghiệp đã tạo điều kiện giúp đỡ dé học viên hoàn thành các đo đạc thực nghiệm trong luận văn. Học viên KHONG MINH TRUONG TOM TAT NỘI DUNG LUẬN VĂN Luận văn “Tính toán lý thuyết, thực nghiệm và mô phỏng quá trình say dung dich sữa tươi trong thiết bị say phun bằng phan mềm Ansys Fluent” sẽ bao gồm những nội dung như Sau: > Tổng quan về công nghệ say phun nói chung và công nghệ say phun sữa, các phương pháp tính toán và mô hình hóa quá trình say phun. Xác định mục tiêu và đối tượng nghiên cứu của đê tài. > Tổng Hợp phương pháp tính toán lý thuyết cho quá trình sấy phun theo phương pháp truyền thống.
> Tiến hành thực nghiệm sấy phun dung dịch sữa tươi nguyên chất để xác định các thông sô đặc trưng thực tê của quá trình say. > Tính kiểm tra thiết bị sấy phun thực tế bằng phương pháp tính toán lý thuyết, để xác định mức độ chính xác của phương pháp tính. > Giới thiệu cơ sở cho việc thiết lập các phương trình chủ đạo nhằm mô hình hóa dong chat lỏng nói chung và cho quá trình say phun theo mô hình Particle-Source-In Cell. > Thực hiện tính toán và mô phỏng quá trình sấy phun trên phần mềm Ansys Fluent.
Sau đó so sánh với kêt quả thực nghiệm đề mức độ chính xác của quá trình mô phỏng. > Kết luận và kiến nghị. ABSTRACT The thesis “Theoretical calculations, experimental of whole milk drying in spray dryer and simulation by using Ansys Fluent software” includes the following main contents: > Overview of spray drying technology and dairy spray drying technology, the methods of calculation and modeling the spray drying process. > Theoretical methods (black box method) for the spray drying process.
> Experimental whole milk spray drying to determine the actual drying propertises. > Overview of the computational modelling of spray drying process. > Simulation of the spray drying process by using ANSYS FLUENT software. Comparison between measurements and predictions by simulation.
LỜI CAM ĐOAN Học viên xin cam đoan răng những nội dung kiên thức và các sô liệu thực nghiệm được trình bay trong luận văn này là do học viên tìm hiệu, nghiên cứu và đo đạc được trên thiết bị thực tế. Trong quá trình làm luận văn, học viên có sử dụng một sô nguôn tài liệu của các tác giả khác ở Việt Nam cũng như trên thê giới, nhưng có trích dân nguôn goc rõ ràng. Tuy nhiên, cũng không thê tránh khỏi những thiêu xót, rat mong được sự góp ý dé tác giả có thê điều chỉnh và đính chính lại. Học viên Không Minh Trưởng DANH MỤC HÌNH ẢNH Hình 1.1(a-b): Sơ dé hệ thống thiết bị say phun hở Hình 1.2: Sơ đồ hệ thống say phun kín sử dụng tác nhân say là nitơ (N2) Hình 2.1: Mô hình cân băng vật chat cho buông say Hình 2.2: Mô hình cân bang năng lượng cho buông say Hình 2.3: Biểu diễn quá trình say Hình 2.4: Biểu diễn quá trình trao đối nhiệt giữa TNS và VLS Hình 2.5: Mối quan hệ giữa thông số cơ bản của vòi phun áp lực - hãng BETE Hình 2.6: Quá trình tốn thất nhiệt ra môi trường.7: Thiết bị say phun SD-05 của hãng LabPlant Hình 2.8: So đồ nguyên lý cấu tạo thiết bi say phun SD-05 Hình 2.9: Kích thước chi tiết của buồng say phun Hình 2.9: (a) Sơ đồ lắp đặt cảm biến; (b) Chi tiết các dòng chất lỏng vào ra thiết bị Hình 3.1: Ma trận Tensor ứng suất Hình 3.2: Mô tả tác động bởi ứng suất pháp tuyến và tiếp tuyến trên phần tử Hình 3.3: Biéu diễn tat cả các ứng suất tác động lên phan tử chất lỏng Hình 3.4: Biểu diễn các dòng nhiệt và công của phân tử theo phương ox Hình 3.5: Mô tả phương pháp kết hợp hai hướng Hình 3.6: Mô tả các thành phần của vận tốc tức thời Hình 3.7: Mô phỏng theo mô hình Stochastic Tracking Hình 3.8: Mô phỏng theo mô hình Particle Cloud Hình 3.9: Biêu diễn các giai đoạn của quá trình sây phun Hình 3.10: Thể hiện cầu trúc của các giọt lỏng và hạt âm Hình 3.11: Quá trình trao đổi nhiệt, chất và động lượng giữa hai pha Hình 4.1: Kích thước chi tiết của buồng say phun Hình 4.2: Các bộ giải có san trên Fluent Hình 4.3: Sơ d6 quy trình tính toán theo bộ giải Pressure Based solver Hình 4.4: Mô hình 3D buông sấy dược xây dựng trên Ansys Design-Modeler, và vị trị gốc hệ trục tọa độ theo mồ hình 3D.5: Mô hình sau khi được chia lưới phan tử Hình 4.6: Các sơ đồ rời rạc hóa Hình 4.7: Việc thực hiện tính toán trên Ansys Workbench Hình 4.8: Bước hội tụ cua pha khí cho trường hợp 6 Sc20 - T = 186,45°C Hình 4.9: Phân bố vận tốc và cường độ rối của pha khí của trường hợp Sc20 - T = 186,45°C Hình 4.10: Phân bố đường kính giọt lỏng với 2 trường hợp lưu lượng khác nhau Hình 4.1 1: Phân bố nhiệt độ với trường hợp lưu lượng phun ở thang đo Sc15(1.12: Phân bố thành phan khối lượng hơi nước với trường hợp lưu lượng phun ở thang đo Sc15(1.13: Phân bố vận tốc pha phân tán cho các trường hợp 4.14: Phân bố đường kính trung bình, vận tốc pha phân tán và cường độ rồi cho trường hợp I( Sc15-T144) Hình 4.15: Phân bố vận tốc hạt theo chiều cao buông say Hình 4.16: Phân bồ vận tốc hạt theo đường kính buông sây có tâm tại tọa độ (0;0,1;0) Hình 4.17: Phân bồ vận tốc hạt theo đường kính buông sây có tâm tại tọa độ (0;0.18: Phân bố vận tốc hạt theo đường kính buông sây có tâm tại tọa độ (0;0,3;0) Hình 4.19: Phân bố nhiệt độ với trường hợp lưu lượng phun ở thang đo Se20(1.20: Phân bố thành phần hơi nước với trường hợp lưu lượng phun ở thang đo Sc20(1 7E-4) Hình 4.21: Phân bó vận tốc pha phân tán cho các trường hợp 4,5,6 - Sc20(1.22: Nhiệt độ phân bồ trên đoạn thang đi qua tâm miệng ra của buồng sấy cho trường hợp 6 (Sc20-T=186,45°C) với tọa độ hai đầu đoạn thăng này là (0,155;0.23: Thành phan hơi nước phân bố trên đoạn thang đi qua tâm miệng ra của buồng say cho trường hợp 6 (Sc20-T=186,45°C) với tọa độ hai đầu đoạn thăng này là (0,155;0.1: Sơ đồ bước giải theo hai cách tinh khác nhau của bộ giải Pressure Based Hình P.2: Nhiệt độ phân bố trên đoạn thắng đi qua tâm miệng ra của buông sấy cho trường hợp 1 (Scl5-T=144°C) với tọa độ hai đầu đoạn thăng này là (0,155;0.3: Nhiệt độ phân bồ trên đoạn thang đi qua tâm miệng ra của buông sấy cho trường hợp 2 (Scl5-T=ló0°C) với tọa độ hai đầu đoạn thăng này là (0,155;0.4: Nhiệt độ phân bồ trên đoạn thang di qua tâm miệng ra của buông sấy cho trường hợp 3 (Scl5-T=I80°C) với tọa độ hai đầu đoạn thang này là (0.5: Nhiệt độ phân bố trên đoạn thắng đi qua tâm miệng ra của buông sấy cho trường hợp 4 (Sc20-T=143°C) với tọa độ hai đầu đoạn thăng này là (0,155;0.6: Nhiệt độ phân bố trên đoạn thắng đi qua tâm miệng ra của buông sấy cho trường hợp 5 (Scl5-T=l66°C) với tọa độ hai đầu đoạn thang này là (0.7: Thanh phan hơi nước phân bố trên đoạn thang đi qua tâm miệng ra của buông say cho trường hop | (Sc15-T=144°C) với tọa độ hai đầu đoạn thăng này là (0,155;0.