(Luận văn thạc sĩ) mô phỏng hệ thống phát điện dùng từ thuỷ động lực trong nhà máy nhiệt điện

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ ĐỒN TRUNG TẮNG MƠ PHỎNG HỆ THỐNG PHÁT ĐIỆN DÙNG TỪ THỦY ĐỘNG LỰC TRONG NHÀ MÁY NHIỆT ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 8520201 SKC006584 Tp Hồ Chí Minh, tháng 05/2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ ĐỒN TRUNG TẮNG MƠ PHỎNG HỆ THỐNG PHÁT ĐIỆN DÙNG TỪ THỦY ĐỘNG LỰC TRONG NHÀ MÁY NHIỆT ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 8520201 Tp Hồ Chí Minh, tháng 05/2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ ĐỒN TRUNG TẮNG MÔ PHỎNG HỆ THỐNG PHÁT ĐIỆN DÙNG TỪ THỦY ĐỘNG LỰC TRONG NHÀ MÁY NHIỆT ĐIỆN NGÀNH: KỸ THUẬT ĐIỆN - 8520201 Hướng dẫn khoa học: PGS.TS LÊ CHÍ KIÊN Tp Hồ Chí Minh, tháng 05/2020 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng QUYẾT ĐỊNH GIAO ĐỀ TÀI Trang i Báo cáo chuyên đề HVTH: Đoàn Trung Tắng BIÊN BẢN HỘI ĐỒNG Trang ii Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Trang iii Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Trang iv Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Trang v Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Trang vi Báo cáo chuyên đề HVTH: Đoàn Trung Tắng LÝ LỊCH KHOA HỌC I LÝ LỊCH SƠ LƯỢC: Họ & tên: Đoàn Trung Tắng Giới tính: Nam Ngày, tháng, năm sinh: 01/05/1986 Nơi sinh: BR-VT Quê quán: Đất Đỏ, Bà Rịa - Vũng Tàu Dân tộc: Kinh Chỗ riêng địa liên lạc: KP.Thanh Bình, Đất Đỏ, BR-VT Điện thoại quan: Fax: II QUÁ TRÌNH ĐÀO TẠO: Trung học chuyên nghiệp: Hệ đào tạo: Nơi học (trường, thành phố): Ngành học: Đại học: Hệ đào tạo: Chính quy Nơi học (trường, thành phố): Đại học Sư Phạm Kỹ Thuật Tp.HCM Ngành học: Điện khí hóa - Cung cấp điện Tên đồ án, luận án môn thi tốt nghiệp: Ngày & nơi bảo vệ đồ án, luận án thi tốt nghiệp: Người hướng dẫn: Thạc sĩ: Hệ đào tạo: Chính Quy Thời gian đào tạo từ : 04/2019 đến 10/2020 Nơi học (trường, thành phố): Trường Đại Học Sư Phạm Kỹ Thuật Thành Phố Hồ Chí Minh Ngành học: Kỹ Thuật Điện Tên luận văn: Mô hệ thống phát điện dùng từ thủy động lực nhà máy nhiệt điện Ngày & nơi bảo vệ luận văn: Ngày 24/05/2020 Trường ĐH Sư Phạm KỹThuật TPHCM Người hướng dẫn: PGS.TS Lê Chí Kiên Trang vii Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Pref=1.03e+005; deltaQMHD=5.00e-003; p41=p31*(1-nEE/nMHD)^(y/(y-1)); T41=T31*(1-nEE-deltaQMHD); %phan tich thiet bi trao doi nhiet toc=1+(1/ns)*(pic^((y-1)/(y*N))-1); T11=T61*toc; T21=(T41-deltaTTDN-T11)*nTDN+T11; T51=T41-((T21-T11)/(1-deltaQTDN)); T71=T11; T81=T61; T91=T11; T101=T61; p51=p41*(1-deltaQTDN*nTDN); %phan tich thiet bi lam lanh p61=p51*(1-deltaQLL*nLL); %phan tich may nen p11=pic*p61; piS=pic^(1/N); tos=1+(1/ns)*(piS^((y-1)/y)-1); p71=p61*piS; p81=p71; p91=p81*piS; p101=p91; %phan tich tuabin %deltaTTDN=50; p21=p11*(1-deltaQTDN*nTDN); pit1=p31/p21; %ty so ap suat - vao tuabin tot=pit1^(et*(y-1)/y); % ty so nhiet vao tuabin T2A1=T21*tot; nt=(1-tot)/(1-pit1^((y-1)/y)); %hieu suat cua tuabin %phan tich nhiet luong G=Qin/(Cp*(T31-T2A1));%luu luong cua chat qua may phat MHD %nhiet luong tai tung nut chu trinh Q11=G*Cp*T11; Q21=G*Cp*T21; Q2A1=G*Cp*T2A1; Q31=G*Cp*T31; Q41=G*Cp*T41; Q51=G*Cp*T51; Q61=G*Cp*T61; PMHD=Q31*nEE; %Dien nang thoat khoi may phat MHD Pc1=N*(T11-T61)*G*Cp; %nang luong ma may nen yeu cau Trang 70 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Wion=G/(MassNo*Seed*60220*1620*12.13); Pion=Wion/nion W11=PMHD-Pc1-Pion W21=Q21-Q2A1 n1=(W11+W21)/Qin %tinh entropy p2A1=p31; S1=(Cp*log(T31/Tref)-Cp*((y-1)/y)*log(p31/Pref))*G; S31=S1+(Cp*log(T31/T2A1)-Cp*((y-1)/y)*log(p31/p2A1))*G; S41=S31+(Cp*log(T41/T31)-Cp*((y-1)/y)*log(p41/p31))*G; S51=S41+(Cp*log(T51/T41)-Cp*((y-1)/y)*log(p51/p41))*G; S61=S51+(Cp*log(T61/T51)-Cp*((y-1)/y)*log(p61/p51))*G; S71=S61+(Cp*log(T71/T61)-Cp*((y-1)/y)*log(p71/p61))*G; S81=S71+(Cp*log(T81/T71)-Cp*((y-1)/y)*log(p81/p71))*G; S91=S81+(Cp*log(T91/T81)-Cp*((y-1)/y)*log(p91/p81))*G; S101=S91+(Cp*log(T101/T91)-Cp*((y-1)/y)*log(p101/p91))*G; S11=S101+(Cp*log(T11/T101)-Cp*((y-1)/y)*log(p11/p101))*G; S21=S11+(Cp*log(T21/T11)-Cp*((y-1)/y)*log(p21/p11))*G; S2A1=S21+(Cp*log(T2A1/T21)-Cp*((y-1)/y)*log(p2A1/p21))*G; x1=[S31 S41 S51 S61 S71 S81 S91 S101 S11 S21 S2A1 S31]; y1=[T31 T41 T51 T61 T71 T81 T91 T101 T11 T21 T2A1 T31]; %thu nghiem T62 p32=4.00e+005; T32=2400; T62=350; p42=p32*(1-nEE/nMHD)^(y/(y-1)); T42=T32*(1-nEE-deltaQMHD); %phan tich thiet bi trao doi nhiet %toc=1+(1/ns)*(pic^((y-1)/(y*N))-1); T12=T62*toc; T22=(T42-deltaTTDN-T12)*nTDN+T12; T52=T42-((T22-T12)/(1-deltaQTDN)); T72=T12; T82=T62; T92=T12; T102=T62; p52=p42*(1-deltaQTDN*nTDN); %phan tich thiet bi lam lanh p62=p52*(1-deltaQLL*nLL); %phan tich may nen p12=pic*p62; %piS=pic^(1/N); %tos=1+(1/ns)*(piS^((y-1)/y)-1); Trang 71 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng p72=p62*piS; p82=p72; p92=p82*piS; p102=p92; %phan tich tuabin %deltaTTDN=50; p22=p12*(1-deltaQTDN*nTDN); pit2=p32/p22; %ty so ap suat - vao tuabin tot=pit2^(et*(y-1)/y); % ty so nhiet vao tuabin T2A2=T22*tot; nt=(1-tot)/(1-pit2^((y-1)/y)); %hieu suat cua tuabin %phan tich nhiet luong G=Qin/(Cp*(T32-T2A2));%luu luong cua chat qua may phat MHD %nhiet luong tai tung nut chu trinh Q12=G*Cp*T12; Q22=G*Cp*T22; Q2A2=G*Cp*T2A2; Q32=G*Cp*T32; Q42=G*Cp*T42; Q52=G*Cp*T52; Q62=G*Cp*T62; PMHD=Q32*nEE; %Dien nang thoat khoi may phat MHD Pc2=N*(T12-T62)*G*Cp; %nang luong ma may nen yeu cau W12=PMHD-Pc2-Pion W22=Q22-Q2A2 n2=(W12+W22)/Qin %tinh entropy p2A2=p32; S2=(Cp*log(T32/Tref)-Cp*((y-1)/y)*log(p32/Pref))*G; S32=S2+(Cp*log(T32/T2A2)-Cp*((y-1)/y)*log(p32/p2A2))*G; S42=S32+(Cp*log(T42/T32)-Cp*((y-1)/y)*log(p42/p32))*G; S52=S42+(Cp*log(T52/T42)-Cp*((y-1)/y)*log(p52/p42))*G; S62=S52+(Cp*log(T62/T52)-Cp*((y-1)/y)*log(p62/p52))*G; S72=S62+(Cp*log(T72/T62)-Cp*((y-1)/y)*log(p72/p62))*G; S82=S72+(Cp*log(T82/T72)-Cp*((y-1)/y)*log(p82/p72))*G; S92=S82+(Cp*log(T92/T82)-Cp*((y-1)/y)*log(p92/p82))*G; S102=S92+(Cp*log(T102/T92)-Cp*((y-1)/y)*log(p102/p92))*G; S12=S102+(Cp*log(T12/T102)-Cp*((y-1)/y)*log(p12/p102))*G; S22=S12+(Cp*log(T22/T12)-Cp*((y-1)/y)*log(p22/p12))*G; S2A2=S22+(Cp*log(T2A2/T22)-Cp*((y-1)/y)*log(p2A2/p22))*G; x2=[S32 S42 S52 S62 S72 S82 S92 S102 S12 S22 S2A2 S32]; y2=[T32 T42 T52 T62 T72 T82 T92 T102 T12 T22 T2A2 T32]; %thu nghiem T63 Trang 72 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng p33=4.00e+005; T33=2400; T63=400; p43=p33*(1-nEE/nMHD)^(y/(y-1)); T43=T33*(1-nEE-deltaQMHD); %phan tich thiet bi trao doi nhiet %toc=1+(1/ns)*(pic^((y-1)/(y*N))-1); T13=T63*toc; T23=(T43-deltaTTDN-T13)*nTDN+T13; T53=T43-((T23-T13)/(1-deltaQTDN)); T73=T13; T83=T63; T93=T13; T103=T63; p53=p43*(1-deltaQTDN*nTDN); %phan tich thiet bi lam lanh p63=p53*(1-deltaQLL*nLL); %phan tich may nen p13=pic*p63; piS=pic^(1/N); %tos=1+(1/ns)*(piS^((y-1)/y)-1); p73=p63*piS; p83=p73; p93=p83*piS; p103=p93; %phan tich tuabin %deltaTTDN=50; p23=p13*(1-deltaQTDN*nTDN); pit3=p33/p23; %ty so ap suat - vao tuabin tot=pit3^(et*(y-1)/y); % ty so nhiet vao tuabin T2A3=T23*tot; nt=(1-tot)/(1-pit3^((y-1)/y)); %hieu suat cua tuabin %phan tich nhiet luong G=Qin/(Cp*(T33-T2A3));%luu luong cua chat qua may phat MHD %nhiet luong tai tung nut chu trinh Q13=G*Cp*T13; Q23=G*Cp*T23; Q2A3=G*Cp*T2A3; Q33=G*Cp*T33; Q43=G*Cp*T43; Q53=G*Cp*T53; Q63=G*Cp*T63; PMHD=Q33*nEE; %Dien nang thoat khoi may phat MHD Trang 73 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Pc3=N*(T13-T63)*G*Cp; %nang luong ma may nen yeu cau W13=PMHD-Pc3-Pion W23=Q23-Q2A3 n3=(W13+W23)/Qin %tinh entropy p2A3=p33; S3=(Cp*log(T33/Tref)-Cp*((y-1)/y)*log(p33/Pref))*G; S33=S2+(Cp*log(T33/T2A3)-Cp*((y-1)/y)*log(p33/p2A3))*G; S43=S32+(Cp*log(T43/T33)-Cp*((y-1)/y)*log(p43/p33))*G; S53=S42+(Cp*log(T53/T43)-Cp*((y-1)/y)*log(p53/p43))*G; S63=S52+(Cp*log(T63/T53)-Cp*((y-1)/y)*log(p63/p53))*G; S73=S62+(Cp*log(T73/T63)-Cp*((y-1)/y)*log(p73/p63))*G; S83=S72+(Cp*log(T83/T73)-Cp*((y-1)/y)*log(p83/p73))*G; S93=S82+(Cp*log(T93/T83)-Cp*((y-1)/y)*log(p93/p83))*G; S103=S92+(Cp*log(T103/T93)-Cp*((y-1)/y)*log(p103/p93))*G; S13=S102+(Cp*log(T13/T103)-Cp*((y-1)/y)*log(p13/p103))*G; S23=S12+(Cp*log(T23/T13)-Cp*((y-1)/y)*log(p23/p13))*G; S2A3=S22+(Cp*log(T2A3/T23)-Cp*((y-1)/y)*log(p2A3/p23))*G; x3=[S33 S43 S53 S63 S73 S83 S93 S103 S13 S23 S2A3 S33]; y3=[T33 T43 T53 T63 T73 T83 T93 T103 T13 T23 T2A3 T33]; %h2=plot(x2,y2,'b-*') h1=plot(x1,y1,'m-*',x2,y2,'bp ',x3,y3,'r+-'); grid xlabel('Truc Entropy') ylabel('Truc Nhiet Do (K)') bar([T61,T62,T63],[n1,n2,n3],0.2) …………………………….End…………………………… Trang 74 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng NỘI DUNG BÀI BÁO MÔ PHỎNG HỆ THỐNG PHÁT ĐIỆN DÙNG TỪ THỦY ĐỘNG LỰC TRONG NHÀ MÁY NHIỆT ĐIỆN SIMULATION OF POWER GENERATION SYSTEMS USING HYDRODYNAMICS IN THERMAL POWER PLANTS Lê Chí Kiên , Đồn Trung Tắng Trường đại học Sư phạm Kỹ thuật Tp.HCM Học viên cao Trường đại học Sư phạm Kỹ thuật Tp.HCM TÓM TẮT Hiện nay, ngành điện nước ta đối mặt với nhiệm vụ khó khăn cung ứng nhu cầu ngày cao khách hàng Bên cạnh vấn đề nhiễm mơi trường, khủng hoảng cung cầu lượng diễn gay gắt Việt Nam mà xuất nhiều quốc gia giới Trong đó, cơng nghệ phát điện máy phát MHD có bước phát triển Nó ứng dụng rộng rãi chu trình phát điện hỗn hợp để nâng cao hiệu suất nhà máy điện Đề tài: "Mô hệ thống phát điện dùng từ thủy động lực nhà máy nhiệt điện" nghiên cứu với mục tiêu kết hợp máy phát điện MHD vào nhà máy nhiệt điện truyền thống nhằm nâng cao hiệu suất phát nhà máy nhiệt điện Các kết mơ tính tốn luận văn cho thấy hiệu suất có máy phát MHD tham gia hiệu suất cải thiện tốt so với mơ hình nhiệt điện truyền thống Từ khóa: máy phát điện MHD, nhiệt điện MHD ABSTRACT Currently, power industry in our country is facing a very difficult task that is to supply increasing demands of customers In addition, the problem of environmental pollution, the crisis of energy supply and demand are taking place very seriously not only in Vietnam but also in many countries around the world In particular, the power generation technology of MHD generators has got new developments It has been widely used in mixed generation cycles to improve the efficiency of power plants Topic: "Simulation of power generation systems using hydrodynamics in thermal power plants" has been studied with the goal of combining MHD generators into traditional thermal power plants in order to improve generation efficiency of thermal power plants The simulation results and calculations in the thesis show that the performance when the MHD generator is involved, the performance is improved better than traditional thermal models Keywords: MHD generators, thermal power plants and MHD Trang 75 Báo cáo chuyên đề đặc biệt trọng đến vấn đề tìm nguồn lượng tái tạo mới, có phương pháp phát điện phát minh nghiên cứu để hoàn thiện nhiều nước đặc biệt ý, phương pháp phát điện từ thủy động lực (MHD) Khác với loại máy phát điện theo nguyên lý cũ, máy phát điện (MHD) vận hành nhiệt độ cao khơng có phận chuyển động MHD phát triển mạnh mẽ nhiệt từ máy phát điện MHD gia nhiệt cho lò nhà máy nhiệt điện Hiện nay, nhà máy nhiệt điện nước ta phát triển số lượng chất lượng, đề tài thực nhằm nghiên cứu kết hợp hệ thống phát điện từ thủy động lực nhiệt điện để tăng cơng suất phát điện cho nhà máy nhiệt điện CƠ SỞ LÝ HVTH: Đồn Trung Tắng Hình 2.2: Đĩa phát MHD [2] Chu trình Brayton Ngày chu trình Brayton[3] nhắc đến động tuốc bin khí Động có ba phần:  Buồng nén khí THUYẾT Nhiệt điện  Buồng đốt tuabin khí  Buồng giãn nở làm quay tuốc bin Turbin khí đại tích hợp công nghệ cao gồm nén, đốt, turbin máy phát (hình 2.3) Hình 2.3: Chu trình Brayton lý tưởng [4] Hiệu suất chu trình Brayton là: Hình 2.1: Sơ đồ nguyên lý nhiệt điện dùng turbin khí chu trình đơn [1] Máy phát dùng đĩa Máy dùng đĩa có hiệu suất cao, máy thiết kế nhỏ gọn Năm 1994, Viện Kỹ thuật Tokyo chế tạo máy dùng đĩa có hiệu suất 22% 1 Hiệu suất chu trình tuabin khí Chu trình tuabin khí thường dùng chu trình Brayton, chu trình bao gồm hai trình đẳng áp hai trình đoạn nhiệt biểu diễn đồ thị: Trang 76 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Áp suất P Vậy ta có η= - Nhiệt độ T P = số ( ) (2.21) PHÂN TÍCH CHU TRÌNH Chu trình phát điện kết hợp máy phát từ thủy động lực (MHD) tuabin khí chu trình kín hoạt động dựa chu trình Brayton 3 Hình 2.4: Quá trình biến đổi trạng thái chất khí chu trình Brayton - 4: Q trình nén đoạn nhiệt - : Quá trình nung nóng đẳng áp – 2: Q trình giãn nở đoạn nhiệt – 3: Quá trình làm lạnh đẳng áp Nhiệt lượng cung cấp cho hệ thống CP.(T1 – T4).Nhiệt lượng thu từ hệ thống C P (T2 – T3) Trong Cp nhiệt dung riêng đẳng áp Hình 3.1: Chu trình kết hợp MHD – tuabin khí Hiệu suất chu trình Brayton xác định thông qua nhiệt dung riêng, nhiệt độ lưu chất Áp suất vào máy phát P 3, áp suất khỏi máy phát P4 ta có tỉ lệ áp suất – vào MHD: = (1 − ) (3.1) Và nhiệt độ chất khí khỏi máy phát MHD xác định bởi: T4 = T3.( - △Q ) (3.2) Phân tích thiết bị trao đổi nhiệt Thiết bị trao đổi nhiệt làm nhiệm vụ thực trao nhiệt chất khí chu trình Áp suất khỏi thiết bị trao đổi nhiệt là: Phân tích máy phát MHD Máy dùng đĩa có hiệu suất ηEE = 35% [6] W ( MHD (2 mà suy T1 = T2( P5 = P4( 1- △QTDN L ) Và nhiệt độ chất khí khỏi thiết bị trao đổi nhiệt xác định bởi: Và ta có suy T4 = T3 ( T5=T4- ) Thay T1 T4 vào phương trình (2.20) ta có η=1– – ( =1− ( ) ) T2– T3 ( Phân tích thiết bị làm lạnh ) Áp suất khỏi thiết bị làm lạnh xác định bởi: (T2 – T3) P6 = P5( 1- △QLL ) (3.5) Trang 77 Báo cáo chuyên đề Phân tích máy nén Áp suất vào máy nén P6 máy nén P1 ta có tỉ số nén máy nén là: Πc = Nếu máy nén có N tầng nén ta có mối quan hệ tỉ số nén máy nén tỉ số nén tầng nén là: Πc = ΠS N Nhiệt độ qua tầng nén xác định bởi: Tỉ số nhiệt độ – vào máy nén là: T1 = T6 ‫ح‬c Phân tích tuabin khí ‫ح‬c = = + ( Hiệu suất η máy nén tỉ số biến đổi enthalpy lý tưởng enthalpy thực tế Tương tự phân tích ta có áp suất khỏi thiết bị trao đổi nhiệt để vào tuabin : P = P ( 1- △Q ) (3.12) Nhiệt độ chất khí vào tuabin: T = (T - △T – T ).η +T (3.13) Tỉ số áp suất – vào tuabin khí là: Thay ℎ ℎ −ℎ = CP(T1i – T6) ℎ = CP(T1 – T6) ta cơng thức tính hiệu = ‫ح‬ Theo công thức (3.9) nhiệt độ chất khí sau nén là: ‫ح‬sj = + ( ηc = HVTH: Đoàn Trung Tắng ηc = − suất ηc = TDN TDN TDN Y Πt = Tỉ số nhiệt độ – vào tuabin tính sau: Từ hiệu suất et = = có = Cuối hiệu suất máy nén xác định công thức sau: Suy et = ηc = (3.10) ‫ح‬ ‫ح‬t = Nhiệt độ chất khí khỏi tuabin là: T2A = T2 ‫ح‬t Hiệu suất tuabin tỉ số biến đổi enthalpy thực tế enthalpy lý tưởng Nếu tính theo ec = = Tti = P T2 ta Thực Lý tưởng P2A T2A Từ ta có ‫ح‬c = cơng thức (3.10) viết lại sau T 2Ai Hình 3.2: Đồ thị T – s tuabin khí Trang 78 Báo cáo chuyên đề ηt = ‫ح‬ ηt = Thay công thức (3.15) vào (3.17) ta có hiệu suất tuabin xác định sau: ηt = Phân tích nhiệt lượng Lưu lượng chất khí qua máy phát MHD G= HVTH: Đoàn Trung Tắng (3.18) f ( )  Nhiệt lượng nút thứ i xác định  Tại nút 4:  Tại nút 5:  Tại nút 6:  Tại tầng nén : S5 = S4 + [CP ln – CP.( Qi = G.CP.Ti Phần điện khỏi máy phát ) ln ].G (3.27) S6 = S5 +[CP ln – CP.( ) ln ].G (3.28) ) ln ].G (3.29) ) ln ].G (3.30) (3.31) PMHD = Q3 Nhiệt lượngcung cấp cho chu trình tuabin : Q4 = Q3 – Q3 (3.21) Năng lượng cung cấp cho máy nén: S Tại nút 3: S4 = S3 + [CP ln – CP.( sau:  là:  S3 = S + [CP ln – CP.( S7 = S6 (3.32) S8 = S7 (3.33) S9 = S8 – Q3 △QMHD PC = N.(T1 – T6).G.CP Năng lượng cần thiết để ion hóa chất khí S10 = S9 + [CP ln là: Pion = (3.34) f f  Điện đưa lên lưới sau chu trình MHD: W1 = PMHD – PC – P ion Điện lên lưới sau chu trình tuabin khí:  Tại nút 1:  Tại nút 2:  Tại nút 2A: S1 = S10 + [CP ln – CP.( S2 = S1 + [CP ln – CP.( S2A = S2 + [CP ln – CP.( W2 = Q2 – Q2A 3.3.7 Tính entropy Tính entropy mẫu: ) ln ].G (3.35) ) ln ].G (3.36) ) ln ].G (3.37) (3.25) ].G (3.26) Trong Tref Pref nhiệt độ áp suất mẫu, thông thường áp suất khí nhiệt độ mơi trường Tính entropy nút: S = [CP ln – CP.( ) ln s Trang 79 Báo cáo chuyên đề TÍNH TỐN VÀ MƠ PHỎNG Bảng 4.2: Kết tính tốn thơng số với T3 = 1800 K 4.1 Bài tốn Dữ liệu tính tốn [2],[5] [6] Bảng 4.1: Dữ liệu tính tốn tốn Đồ thị T –S với T3 = 1800 K Kết tính tốn HVTH: Đồn Trung Tắng Hình 4.1: Đồ thị T –S với T3 = 1800 K Hiệu suất chu trình η= So sánh kết Nghiên cứu phân tích báo cáo Đại học Nagaoka Niigata Nhật Bản điều kiện nhiệt độ áp suất η = 55,2% Hình 4.2: Kết phân tích chu trình MHD dạng đơn [7] 4.2 Bài tốn Giữ thơng số đầu vào toán Cho T3 = 2400oK, thay đổi T6 = 350oK 400oK Dữ liệu tính tốn tốn Bảng 4.3: Dữ liệu tính tốn tốn Kết thông số nút Bảng 4.4: Kết tính tốn thơng số với T = 350 K Trang 80 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng Bảng 4.5: Kết tính tốn thơng số Truc Nhiet Do (K) với T6 = 400 K Hình 4.3: Đồ thị T –S với T6 0 thay đổi 300 K, 350 K 400 K Hiệu suất chu trình T6 = 350 K η= Hiệu suất chu trình T6 = 400 K K.K η= = 60,76% ' 0.7 0.6 0.5 0.4 0.3 0.2 0.1 300 350 400 Hình 4.4: Biểu đồ biểu thị mối quan hệ hiệu suất chu trình nhiệt độ vào máy nén (T6) KẾT LUẬN Từ kết mơ tính toán hai toán thấy hiệu suất chu trình kết hợp cao chu trình đơn Hiệu suất phụ thuộc vào nhiệt độ mơi chất vào máy phát MHD mà cịn phụ thuộc vào nhiệt độ môi chất đưa vào máy nén Do tính tốn thiết kế thiết bị chu trình người thiết kế cần lưu ý đến tính chất vật lý mơi chất, quan hệ nhiệt độ entropy môi chất nút chu trình nhằm tính tốn thiết bị xác để hiệu suất chu trình nâng cao góp phần giải tốn lượng tương lai Bên cạnh hiệu suất chu trình cịn phụ thuộc lớn vào nhiệt độ môi chất vào máy nén khí sau làm lạnh Nếu mơi chất làm lạnh sâu hiệu suất chu trình cao nhiên việc làm lạnh mơi chất có giới hạn khơng phải vơ tận cịn phụ thuộc vào nhiệt độ mơi trường Ngồi yếu tố nhiệt độ tỷ số nén số tầng nén máy nén tác động đến hiệu suất chu trình Trang 81 Báo cáo chuyên đề HVTH: Đoàn Trung Tắng TÀI LIỆU THAM KHẢO [1] Nguyễn Công Hân, Nhà máy nhiệt điện, NXB Khoa học kỹ thuật, 2002, 200 trang [2] Richard J Rosa “ Magnetohydrodynamic Energy Conversion” copyright 1987 by hemisphere publishing coporation, Printed in USA, 234 pages [3] https://vi.wikipedia.org/ [4] Motoo Ishikawa, Susumu Takebe, Fukasi Kumura Yoshitaka Inui and Juro Umoto“ Application of MHD – Brayton cycle to fusion reactors”vol4cap36 [5] Ron J Litchford and Nobuhiro Harada, 2011“ Multi-MW Closed Cycle MHD Nuclear Space Power Via Nonequilibrium He/Xe Working Plasma”, Proceedings of Nuclear and Emerging Technologies for Space 2011 [6] Jack D.Mattingly“ Elements of Gas Turbine Propulsion” copyright 1996 by McGraw – Hill, Inc, Printed in Singapore, 960 pages [7] Nob Harada, Le Chi Kien, and M Hishikawa “Basic Studies on Closed Cycle MHD Power Generation System for Space Application” AIAA 2004-2365 Tác giả chịu trách nhiệm viết: Họ tên: Lê Chí Kiên Đơn vị: Khoa Điện – Điện tử, ĐHSPKT TP.HCM Điện thoại: 0987 67 30 30 Email: kienlc@hcmute.edu.vn Trang ... hệ thống phát điện từ thủy động lực nhiệt điện để tăng cơng suất phát điện cho nhà máy nhiệt điện 1.2 Giới thiệu máy phát điện MHD nghiên cứu MHD Máy phát điện MHD hệ thống chuyển nhiệt hay động. .. hình hệ thống phát điện dùng từ thủy động lực kết hợp nhiệt điện Chương 4: Tính tốn mơ Tính tốn mơ hình mơ hình hệ thống phát điện dùng từ thủy động lực kết hợp nhiệt điện thông số cụ thể, mô thông... điện máy phát MHD có bước phát triển Nó ứng dụng rộng rãi chu trình phát điện hỗn hợp để nâng cao hiệu suất nhà máy điện Đề tài: "Mô hệ thống phát điện dùng từ thủy động lực nhà máy nhiệt điện"