Nghiên cứu quá trình cháy dầu fo trong thiết bị vòi phun dòng rối xoáy

BỘ GIÁO DỤC VÀ ĐÀO TẠO ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN THANH HÀO NGHIÊN CỨU QUÁ TRÌNH CHÁY DẦU FO TRONG THIẾT BỊ VÒI PHUN DÒNG RỐI XOÁY CHUYÊN NGÀNH: CHẾ TẠO MÁY MÃ SỐ NGÀNH: 2.01.00 LUẬN VĂN THẠC SĨ THÀNH PHỐ HỒ CHÍ MINH - 2004 CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH CÁN BỘ HƯỚNG DẪN KHOA HỌC: PGS.TS NGUYỄN THANH NAM CÁN BỘ CHẤM NHẬN XÉT 1: TS NGUYỄN VĂN TUYÊN CÁN BỘ CHẤM NHẬN XÉT 2: TS BÙI NGỌC HÙNG Luận văn thạc só bảo vệ HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ – TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày tháng năm 2004 Có thể tìm hiểu luận văn Thư Viện Cao Học - Trường Đại Học Bách Khoa Đại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM PHÒNG ĐÀO TẠO SAU ĐẠI HỌC ĐỘC LẬP – TỰ DO – HẠNH PHÚC TP.HCM, ngày tháng năm 2004 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên : NGUYỄN THANH HÀO Phái: Nam Ngày tháng năm sinh: 06 / 11 / 1973 Nơi sinh: Đắc Lắc Chuyên ngành : CƠ KHÍ CHẾ TẠO MÁY Mã số: 2.01.00 I- TÊN ĐỀ TÀI: NGHIÊN CỨU QUÁ TRÌNH CHÁY DẦU FO TRONG THIẾT BỊ VÒI PHUN DÒNG RỐI XOÁY II- NHIỆM VỤ VÀ NỘI DUNG: Thiết lập mô hình số trình cháy nhiên liệu dầu FO thiết bị vòi phun dòng rối xoáy Xác định tối ưu hoá thông số buồng đốt đảm bảo cho trình cháy nâng cao tuổi thọ cho thiết bị III- NGÀY GIAO NHIỆM VỤ: 14 / / 2004 IV- NGÀY HOÀN THÀNH NHIỆM VỤ: 10 / / 2004 V- HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: PGS.TS NGUYỄN THANH NAM CÁN BỘ HƯỚNG DẪN PGS.TS.NGUYỄN THANH NAM CHỦ NHIỆM NGÀNH BỘ MÔN QUẢN LÝ CHUYÊN NGÀNH TS.PHẠM NGỌC TUẤN PGS.TS.TRẦN DOÃN SƠN Nội dung đề cương luận văn thạc só Hội Đồng Chuyên Ngành thông qua Ngày tháng năm 2004 PHÒNG ĐÀO TẠO SAU ĐẠI HỌC KHOA QUẢN LÝ NGÀNH Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM LỜI CẢM ƠN Xin chân thành cảm ơn hướng dẫn tận tình, chu đáo PGS TS Nguyễn Thanh Nam suốt trình thực luận văn Xin chân thành cảm ơn Thầy Cô dành thời gian quý báu để duyệt đóng góp ý kiến giúp cho luận văn hoàn thiện Xin chân thành cảm ơn đến Thầy Cô giảng dạy chuyên ngành Cao học Chế tạo máy Trường đại học Bách khoa Tp Hồ Chí Minh niên khoá 2002 ÷ 2004 nhiệt tình truyền thụ kiến thức kinh nghiệm suốt thời gian học tập nghiên cứu trường Xin chân thành cảm ơn giúp đỡ Anh Chị đồng nghiệp trình thực luận văn Xin chân thành cảm ơn công ty AUREOLE CSD INC tạo điều kiện thuận lợi suốt trình học tập thực luận văn Xin gửi lời cảm ơn đến gia đình người thân giúp đỡ, động viên, hỗ trợ khích lệ tinh thần suốt thời gian học tập Sau cùng, xin chúc quý Thầy Cô, Anh Chị đồng nghiệp sức khỏe thành đạt Tp Hồ Chí Minh, ngày 10 tháng năm 2004 Người thực Nguyễn Thanh Hào Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM LỜI NÓI ĐẦU nh hưởng xoáy lên dòng lưu chất biết đến từ nhiều năm trước Các nghiên cứu dòng xoáy ứng dụng nhiều ngành kỹ thuật cyclone, thiết bị phun dùng nông nghiệp, thiết bị trao đổi nhiệt… Tuy nhiên ứng dụng cho hiệu không cao mong muốn trước khả nghiên cứu sở để thiết kế thiết bị máy móc thực nghiệm, điều đòi hỏi tiêu tốn nhiều công sức phương tiện thực mô hình Ứng dụng quan trọng dòng xoáy công nghiệp có lẽ thiết bị buồng đốt turbines, lò hơi, động đốt, lò nung công nghiệp số chủng loại thiết bị khác Thông thường dòng xoáy có tác dụng làm tăng mức độ rối buồng đốt, điều đồng nghóa với việc làm tăng hiệu hoà trộn không khí nhiên liệu Bên cạnh đó, nhiệt trị dầu FO cao, giá thành rẻ ngày việc sử dụng lượng nhiệt từ việc đốt dầu FO phổ biến Tuy nhiên, để đốt dầu FO cần phải dùng thiết bị sấy phun để phun dầu thành hạt bụi nhỏ vào buồng đốt, đường kính hạt dầu bé trình cháy tốt Do việc phun dầu thành hạt bụi quan trọng Có thể phun dầu bơm áp lực cao, nước khí nén có áp suất cao… kết hợp với việc tăng nhanh trình khuếch tán hỗn lưu đối lưu dầu với không khí cách tạo dòng xoáy đảm bảo cho dòng dầu sau miệng phun dễ dàng bắt lửa Số lượng khối lượng nghiên cứu thực nghiệm trình cháy dầu FO thiết bị vòi phun rối xoáy giảm thiểu đáng kể việc thiết kế thiết bị buồng đốt dầu FO đơn giản nhờ vào tính toán kiểm chứng Mô hình số trình cháy dầu FO thiết bị vòi phun rối xoáy hoàn toàn thực phương pháp tính toán đại phương pháp sai phân hữu hạn phương pháp khối hữu hạn với hỗ trợ máy tính Các mô hình tính toán cho phép xác định thông số động học lửa, dự báo giới hạn thay đổi thông số, đảm bảo tối ưu hiệu trình cháy Xuất phát từ lý trên, tác giả thực đề tài "Nghiên cứu trình cháy dầu FO thiết bị vòi phun dòng rối xoáy" Khoa khí Trường đại học Bách khoa Tp Hồ Chí Minh phân công PGS TS Nguyễn Thanh Nam trực tiếp hướng dẫn Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM SUMMARY "Study the process of burning FO fuel in swirling turbulent combustion chamber" Burning FO fuel in swirling turbulent combustion chamber is encountered in many applications such as turbines, boilers, diesel engines, industrial furnaces and liquid-fuelled rocket engines In spite of its importance, our understanding of the many aspects of spray combustion phenomena, such as particle dispersion, turbulence modification by particles, evaporation/mixing and combustion are relatively limited This is because the experimental investigation of such complex particle-laden flows is often very time consuming, making it impossible to conduct a detailed interrogation of such flows On the other hand, many of the numerical studies performed in the past have employed approximate models which are not close to reality, thus limiting the scope of their application Analytical studies are often restricted to simple flows, making it difficult to extend the results of these studies to any practical flows This study introduced the numerical model of the burning FO fuel in swirling turbulent combustion chamber with the real initial velocity distribution The results are compared and agreed with the experimental data confirming the possibility of the used numerical model in practice The thesis includes chapers that presented following details: Chapter 1: Fundamentals of combustion processes Chapter 2: Mathemetical model of non-isothermal swirling turbulent flows in combustion chamber Chapter 3: Numerical model of non-isothermal swirling turbulent flows in combustion chamber Chapter 4: Analytical results and definition optimal parameters Chapter 5: Estimation the possibility of the used in practice Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM MỤC LỤC Trang LỜI NÓI ĐẦU CHƯƠNG TỔNG QUAN VỀ QUÁ TRÌNH CHÁY DÒNG PHUN DẦU FO RỐI XOÁY 1.1 Cơ sở lý thuyết cháy 1.1.1 Cơ sở nhiệt động hoá học 12 12 a) Tốc độ phản ứng hoá học 12 b) nh hưởng nồng độ đến tốc độ phản ứng hoá học 13 c) nh hưởng nhiệt độ đến tốc độ phản ứng hoá học 14 d) nh hưởng áp suất đến tốc độ phản ứng hoá học 15 e) Cân hoá học 16 1.1.2 Lý thuyết bắt cháy nhiên liệu 17 a) Giai đoạn bắt lửa 17 b) Giai đoạn cháy kiệt 18 1.2 Cơ sở lý thuyết dòng phun dầu FO rối xoáy 19 1.2.1 Động học dòng phun dầu FO rối xoáy 19 1.2.2 Phương thức cháy dầu FO rối xoáy 26 1.3 Khảo sát số công trình nghiên cứu trình cháy 28 1.3.1 Mô hình FSD (Flamelet Surface Density Models) 28 a) Mô hình đại số Bray-Libby-Moss (1990) 29 b) Mô hình sử dụng phương trình chuyển đổi để xác định giá trị Σ 29 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM 1.3.2 Mô hình PDF (Probability Density Function Models) 30 1.3.3 Mô hình EBU (Eddy Break Up Models) 30 1.3.4 Mô hình EDM (Eddy Dissipation Combustion Models) 31 1.3.5 Mô hình EDC (Eddy Dissipation Concept Model) 31 1.3.6 Mô hình EDK (Eddy Dissipation Kinetic Model) 32 1.3.7 Mô hình RANS (Reynolds Averaged Navier-Stokes Model) 32 1.3.8 Mô hình số DNS (Direct Numerical Simulation Model) 33 1.3.9 Mô hình số LES (Large Eddy Simulation Model) 34 1.4 Nhiệm vụ luận văn tốt nghiệp 34 CHƯƠNG MÔ HÌNH TOÁN CỦA DÒNG PHUN RỐI XOÁY KHÔNG ĐẲNG NHIỆT TRONG BUỒNG ĐỐT 2.1 Phương trình liên tục 37 2.2 Phương trình cân moment động lượng - phương trình Reynolds 38 2.3 Phương bảo toàn lượng 38 2.4 Phương trình hệ số hỗn hợp cháy 39 2.5 Điều kiện biên 41 CHƯƠNG MÔ PHỎNG SỐ QUÁ TRÌNH CHÁY CỦA DÒNG PHUN RỐI XOÁY KHÔNG ĐẲNG NHIỆT TRONG BUỒNG ĐỐT 3.1 Phương pháp giải hệ phương trình dòng phun rối xoáy hai pha không đẳng nhiệt 44 3.2 Thuật toán giải hệ phương trình dòng phun rối xoáy hai pha không đẳng nhiệt 50 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM 3.3 Trình tự tính toán 52 3.3.1 Giải phương trình tìm biến f 53 3.3.2 Giải phương trình tìm biến w 57 3.3.3 Giải phương trình tìm biến Ω 61 3.3.4 Giải phương trình tìm biến ψ 65 3.3.5 Giải phương trình tìm biến u 69 CHƯƠNG PHÂN TÍCH KẾT QUẢ TÍNH TOÁN XÁC ĐỊNH CÁC THÔNG SỐ TỐI ƯU 4.1 Các kết tính toán 71 4.2 Phân tích kết tính toán 74 4.2.1 Kết đo profile vận tốc từ thực nghiệm Luc Vervisch 74 4.2.2 Kết profile vận tốc từ phương pháp tính toán mô hình số mô dòng phun dầu FO rối xoáy hai pha không đẳng nhiệt 75 4.2.3 So sánh kết phương pháp đo từ thực nghiệm phương pháp mô số dòng phun dầu FO rối xoáy hai pha không đẳng nhiệt 76 4.3 Xác định thông số tối ưu phục vụ thiết kế chế tạo thiết bị buồng đốt đảm bảo trình cháy 77 4.3.1 Xác định thông số hình học tối ưu hệ số cường độ xoáy S = 78 4.3.2 Xác định thông số hình học tối ưu hệ số cường độ xoáy S = 80 4.3.3 Xác định thông số hình học tối ưu hệ số cường độ xoáy thay đổi 81 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM CHƯƠNG ĐÁNH GIÁ KHẢ NĂNG ỨNG DỤNG TRONG THỰC TẾ 5.1 Mô dòng phun dầu FO rối xoáy buồng đốt 84 5.2 Tính toán kích thước buồng đốt 89 KẾT LUẬN 92 PHỤ LỤC 95 TÀI LIỆU THAM KHẢO 127 10 Luận Văn Cao Học end Trường Đại Học Bách Khoa Tp.HCM X41=inv(B41)*C41; C42=(CW42.*X41)+(CE42.*X43)+DP42; B42=spdiags(A42,[-1,0,1],4,4); full(B42); X42=inv(B42)*C42; C43=(CW43.*X42)+(CE43.*X44)+DP43; B43=spdiags(A43,[-1,0,1],4,4); full(B43); X43=inv(B43)*C43; C44=(CW44.*X43)+DP44; B44=spdiags(A44,[-1,0,1],4,4); full(B44); X44=inv(B44)*C44; Mt04=[X40 X41 X42 X43 X44]; 113 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM function Mt05=matran05() clear all home format short % % Giai phuong trinh van toc theo phuong doc truc % dong phun roi xoay dau FO % matran04; mtbb5=[1 1 1]; mtbn5=[1 1 1]; DPb5=[0.4286 0.2126 0.2520 0.2211 0]; DPn5=[0.3333 0.2098 0.2519 0.1841 0]; source05=[mtbn5;matran04;mtbb5]; n=5; % So nut theo phuong Tay Dong m=4; % So nut theo phuong Nam Bac rP=1.5; for j=1:n for i=1:m DP511(i,j)=(source05(i,j)-source05(i+2,j))/2; end end for j=1:n for i=1:m DP522(i,j)=DP511(i,j)/rP; end rP=rP+1; end DP522(2,:)=[]; DP522(1,:)=[]; DP533=[DPb5;DP522;DPn5]; Mt05=DP533; 114 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM % -% Do thi he so hon hop chay f % dong phun roi xoay dau FO % -function dothi01(x,y,z) matran01; x=0:75:300; y=0:50/3:50; [X,Y]=meshgrid(x,y); Z= matran01; mesh(X,Y,Z) title('He so hon hop chay dau FO') xlabel('z') ylabel('r') zlabel('f') % -% Do thi van toc theo phuong tiep tuyen % dong phun roi xoay dau FO % -function dothi02(x,y,z) matran02; x=0:75:300; y=0:50/3:50; [X,Y]=meshgrid(x,y); Z= matran02; mesh(X,Y,Z) title('Van toc theo phuong tiep tuyen') xlabel('z') ylabel('r') zlabel('w') % % Do thi cuong xoay Omega % dong phun roi xoay dau FO % function dothi03(x,y,z) matran03; x=0:75:300; y=0:50/3:50; [X,Y]=meshgrid(x,y); 115 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM Z=matran03; mesh(X,Y,Z) title('Cuong xoay hon hop dau FO roi xoay') xlabel('z') ylabel('r') zlabel('Omega') % % Do thi ham dong Nyu % dong phun roi xoay dau FO % function dothi04(x,y,z) matran04; x=0:75:300; y=0:50/3:50; [X,Y]=meshgrid(x,y); Z=matran04; mesh(X,Y,Z) title('Ham dong hon hop dau FO roi xoay') xlabel('z') ylabel('r') zlabel('Nyu') % % Do thi van toc theo phuong doc truc % dong phun roi xoay dau FO % function dothi05(x,y,z) matran05; x=0:75:300; y=0:50/3:50; [X,Y]=meshgrid(x,y); Z=matran05; mesh(X,Y,Z) title('Van toc theo phuong doc truc') xlabel('z') ylabel('r') zlabel('u') 116 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM % -% Tao giao dien % Mo phong dau FO roi xoay % -function varargout = MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY(varargin) % MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY M-file for MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY.fig % MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY, by itself, creates a new MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY or raises the existing % singleton* % % H = MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY returns the handle to a new MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY or the handle to % the existing singleton* % % MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY('CALLBACK',hOb ject,eventData,handles, ) calls the local % function named CALLBACK in MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY.M with the given input arguments % % MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY('Property','Value', ) creates a new MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY or raises the % existing singleton* Starting from the left, property value pairs are % applied to the GUI before MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY_OpeningFunction gets called An % unrecognized property name or invalid value makes property application % stop All inputs are passed to MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY_OpeningFcn via varargin % % *See GUI Options on GUIDE's Tools menu Choose "GUI allows only one % instance to run (singleton)" % 117 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM % See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to modify the response to help MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY % Last Modified by GUIDE v2.5 04-May-2004 10:17:21 % Begin initialization code - DO NOT EDIT gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, 'gui_Singleton', gui_Singleton, 'gui_OpeningFcn', @MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY_OpeningFcn, 'gui_OutputFcn', @MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY_OutputFcn, 'gui_LayoutFcn', [] , 'gui_Callback', []); if nargin & isstr(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); end % End initialization code - DO NOT EDIT % % - Executes just before MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY is made visible function MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY_OpeningFcn(hObj ect, eventdata, handles, varargin) % This function has no output args, see OutputFcn % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % varargin command line arguments to MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY (see VARARGIN) 118 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM % Choose default command line output for MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY handles.output = hObject; % Update handles structure guidata(hObject, handles); % This sets up the initial plot - only when we are invisible % so window can get raised using MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY if strcmp(get(hObject,'Visible'),'off') matran01; x=(0:75:300); y=0:50/3:50; mesh(X,Y,matran01); title('HE SO HON HOP CHAY DAU FO') xlabel('z') ylabel('r') end % UIWAIT makes MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY wait for user response (see UIRESUME) % uiwait(handles.figure1); % % - Outputs from this function are returned to the command line function varargout = MOPHONGSODONGPHUNHONHOPDAUFOROIXOAY_OutputFcn(hObjec t, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT); % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Get default command line output from handles structure varargout{1} = handles.output; % function FileMenu_Callback(hObject, eventdata, handles) % hObject handle to FileMenu (see GCBO) 119 Luận Văn Cao Học Trường Đại Học Baùch Khoa Tp.HCM % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % function OpenMenuItem_Callback(hObject, eventdata, handles) % hObject handle to OpenMenuItem (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) file = uigetfile('*.fig'); if ~isequal(file, 0) open(file); end % function PrintMenuItem_Callback(hObject, eventdata, handles) % hObject handle to PrintMenuItem (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) printdlg(handles.figure1) % function CloseMenuItem_Callback(hObject, eventdata, handles) % hObject handle to CloseMenuItem (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) selection = questdlg(['Close ' get(handles.figure1,'Name') '?'], ['Close ' get(handles.figure1,'Name') ' '], 'Yes','No','Yes'); if strcmp(selection,'No') return; end delete(handles.figure1) % % - Executes during object creation, after setting all properties function vantocu_CreateFcn(hObject, eventdata, handles) % hObject handle to vantocu (see GCBO) 120 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER usewhitebg = 1; if usewhitebg set(hObject,'BackgroundColor','white'); else set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); end % function vantocu_Callback(hObject, eventdata, handles) % hObject handle to vantocu (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of vantocu as text % str2double(get(hObject,'String')) returns contents of vantocu as a double vantocu = str2double(get(hObject, 'String')); if isnan(vantocu) set(hObject, 'String', 0); errordlg('Input must be a number','Error'); end data = getappdata(gcbf, 'metricdata'); data.vantocu = vantocu; setappdata(gcbf, 'metricdata', data); % % - Executes during object creation, after setting all properties function vantocw_CreateFcn(hObject, eventdata, handles) % hObject handle to vantocw (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER 121 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM usewhitebg = 1; if usewhitebg set(hObject,'BackgroundColor','white'); else set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); end % function vantocw_Callback(hObject, eventdata, handles) % hObject handle to vantocw (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of vantocw as text % str2double(get(hObject,'String')) returns contents of vantocw as a double vantocw = str2double(get(hObject, 'String')); if isnan(vantocw) set(hObject, 'String', 0); errordlg('Input must be a number','Error'); end data = getappdata(gcbf, 'metricdata'); data.vantocw = vantocw; setappdata(gcbf, 'metricdata', data); % % - Executes during object creation, after setting all properties function vantocv_CreateFcn(hObject, eventdata, handles) % hObject handle to vantocv (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER usewhitebg = 1; if usewhitebg set(hObject,'BackgroundColor','white'); else set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 122 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM end % function vantocv_Callback(hObject, eventdata, handles) % hObject handle to vantocv (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of vantocv as text % str2double(get(hObject,'String')) returns contents of vantocv as a double vantocv = str2double(get(hObject, 'String')); if isnan(vantocv) set(hObject, 'String', 0); errordlg('Input must be a number','Error'); end data = getappdata(gcbf, 'metricdata'); data.vantocv = vantocv; setappdata(gcbf, 'metricdata', data); % % - Executes on button press in pushbutton1 function pushbutton1_Callback(hObject, eventdata, handles) % hObject handle to pushbutton1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) axes(handles.axes1); cla; popup_sel_index = get(handles.popupmenu1, 'Value'); switch popup_sel_index case matran01; x=(0:75:300); y=0:50/3:50; mesh(X,Y,matran01); title('HE SO HON HOP CHAY DAU FO') xlabel('z') ylabel('r') case matran02; 123 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM data = getappdata(gcbf, 'metricdata'); matran02=data.vantocw*matran02; x=(0:75:300); y=0:50/3:50; mesh(X,Y,matran02); title('VAN TOC PHUONG TIEP TUYEN CUA HON HOP DAU FO') xlabel('z') ylabel('r') case matran03; x=(0:75:300); y=0:50/3:50; mesh(X,Y,matran03); title('CUONG DO XOAY CUA HON HOP DAU FO') xlabel('z') ylabel('r') case matran04; x=(0:75:300); y=0:50/3:50; mesh(X,Y,matran04); title('DO THI HAM DONG CUA HON HOP DAU FO') xlabel('z') ylabel('r') case matran05; data = getappdata(gcbf, 'metricdata'); matran05=data.vantocw*matran05; x=(0:75:300); y=0:50/3:50; mesh(X,Y,matran05); title('VAN TOC PHUONG DOC TRUC CUA HON HOP DAU FO') xlabel('z') ylabel('r') end % % - Executes during object creation, after setting all properties function popupmenu1_CreateFcn(hObject, eventdata, handles) % hObject handle to popupmenu3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB 124 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM % handles empty - handles not created until after all CreateFcns called % Hint: popupmenu controls usually have a white background on Windows % See ISPC and COMPUTER if ispc set(hObject,'BackgroundColor','white'); else set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); end set(hObject, 'String', {'mesh(X,Y,matran01)', 'mesh(X,Y,matran02)', 'mesh(X,Y,matran03)', 'mesh(X,Y,matran04)', 'mesh(X,Y,matran05)'}); % % - Executes on selection change in popupmenu3 function popupmenu1_Callback(hObject, eventdata, handles) % hObject handle to popupmenu3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: contents = get(hObject,'String') returns popupmenu3 contents as cell array % contents{get(hObject,'Value')} returns selected item from popupmenu3 % % - Executes on button press in pushbutton2 function pushbutton2_Callback(hObject, eventdata, handles) % hObject handle to pushbutton2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) initialize_gui(gcbf, handles); % function initialize_gui(fig_handle, handles) data.vantocu = 0; data.vantocw = 0; data.vantocv = 0; setappdata(fig_handle, 'metricdata', data); 125 Luận Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM set(handles.vantocu, 'String', data.vantocu); set(handles.vantocw, 'String', data.vantocw); set(handles.vantocv, 'String', data.vantocv); set(handles.text22, 'String', 'm/s'); set(handles.text23, 'String', 'm/s'); set(handles.text21, 'String', 'm/s'); 126 Luaän Văn Cao Học Trường Đại Học Bách Khoa Tp.HCM LÝ LỊCH TRÍCH NGANG Họ tên: Nguyễn Thanh Hào Năm sinh: 06-11-1973 Nơi sinh: Đắc Lắc Địa liên lạc: Phòng số 06, lầu 19, 115 Nguyễn Huệ, Quận 1, Tp Hồ Chí Minh Điện thoại: 0913645203 08.8219383 Quá trình đào tạo: 1992 - 1997: Sinh viên Khoa Cơ Khí Trường Đại Học Bách Khoa Tp Hồ Chí Minh 2002 - 2004: Học viên cao học chuyên ngành Chế Tạo Máy Trường Đại Học Bách Khoa Tp Hồ Chí Minh Nơi công tác nay: Bộ phận thiết kế, Công ty Aureole CSD INC 131 ... ĐỀ TÀI: NGHIÊN CỨU QUÁ TRÌNH CHÁY DẦU FO TRONG THIẾT BỊ VÒI PHUN DÒNG RỐI XOÁY II- NHIỆM VỤ VÀ NỘI DUNG: Thiết lập mô hình số trình cháy nhiên liệu dầu FO thiết bị vòi phun dòng rối xoáy Xác... lượng nghiên cứu thực nghiệm trình cháy dầu FO thiết bị vòi phun rối xoáy giảm thiểu đáng kể việc thiết kế thiết bị buồng đốt dầu FO đơn giản nhờ vào tính toán kiểm chứng Mô hình số trình cháy dầu. .. không đẳng nhiệt không gian ba chiều Thiết lập mô hình số trình cháy nhiên liệu dầu FO thiết bị vòi phun dòng rối xoáy dựa sở nghiên cứu toán dòng phun rối xoáy hai pha không đẳng nhiệt không