Nghiên Cứu Hệ Thống Phát Điện Từ Thủy Động Lực (Mhd) Sử Dụng Năng Lượng Địa Nhiệt Với Chu Trình Kết Hợp.pdf

Untitled 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Ĩ LÂM HỮU DŨNG NGHIÊN CỨU HỆ THỐNG PHÁT ÐIỆN TỪ THỦY ÐỘNG LỰC (MHD) SỬ DỤNG NĂNG LƯỢNG ÐỊA NHIỆT VỚ[.]

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Ĩ LÂM HỮU DŨNG NGHIÊN CỨU HỆ THỐNG PHÁT ÐIỆN TỪ THỦY ÐỘNG LỰC (MHD) SỬ DỤNG NĂNG LƯỢNG ÐỊA NHIỆT VỚI CHU TRÌNH KẾT HỢP NGÀNH: KỸ THUẬT ĐIỆN - 60520202 S K C0 Tp Hồ Chí Minh, tháng 04/2018 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Ĩ LÂM HỮU DŨNG NGHIÊN CỨU HỆ THỐNG PHÁT ĐIỆN TỪ THỦY ĐỘNG LỰC (MHD) SỬ DỤNG NĂNG LƯỢNG ĐỊA NHIỆT VỚI CHU TRÌNH KẾT HỢP NGÀNH: KỸ THUẬT ĐIỆN - 60520202 Tp Hồ Chí Minh, tháng 04/2018 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Ĩ LÂM HỮU DŨNG NGHIÊN CỨU HỆ THỐNG PHÁT ĐIỆN TỪ THỦY ĐỘNG LỰC (MHD) SỬ DỤNG NĂNG LƯỢNG ĐỊA NHIỆT VỚI CHU TRÌNH KẾT HỢP NGÀNH: KỸ THUẬT ĐIỆN - 60520202 Hướng dẫn khoa học: PGS.TS LÊ CHÍ KIÊN Tp Hồ Chí Minh, tháng 04/2018 Luận Văn LÝ LỊCH KHOA HỌC I LÝ LỊCH SƠ LƯỢC: Họ tên: Lâm Hữu Dũng Giới tính: Nam Ngày, tháng, năm sinh: 20 – 07 – 1979 Nơi sinh: Quảng Nam Quê quán: Quảng Nam Dân tộc: Kinh Chổ riêng địa liên lạc: 1500/15B, Đơng Thành, Tân Đơng Hiệp, Dĩ An, Bình Dương Điện thoại quan: Điện thoại riêng: 0909196793 Fax: E-mail: lamdung.ptc44@gmail.com II QUÁ TRÌNH ĐÀO TẠO: Đại học: Hệ đào tạo: VLVH Thời gian đào tạo từ 10/2006 đến 05/2010 Nơi học: Đại học sư phạm kỹ thuật thành phố Hồ Chí Minh Ngành học: Điện cơng nghiệp Năm tốt nghiệp: 2010 Cao học: Hệ đào tạo: Chính quy Thời gian đào tạo từ 10/2016 đến 4/2018 Nơi học (trường, thành phố): Đại Học Sư Phạm Kỹ Thuật Tp Hồ Chí Minh Ngành học: Kỹ Thuật Điện Tên Luận văn tốt nghiệp: Nghiên Cứu Hệ Thống phát điện Từ Thủy Động Lực (MHD) Sử Dụng Năng Lượng địa nhiệt với chu trình kết hợp Thời gian & nơi bảo vệ Luận văn: Tháng năm 2018, Đại Học Sư Phạm Kỹ Thuật Tp Hồ Chí Minh Người hướng dẫn: PGS.TS Lê Chí Kiên HVTH: Lâm Hữu Dũng i GVHD: PGS.TS Lê Chí Kiên Luận Văn III Q TRÌNH CƠNG TÁC CHUN MƠN KỂ TỪ KHI TỐT NGHIỆP ĐẠI HỌC: Thời gian Nơi công tác Công việc đảm nhiệm 09/2010 Công Ty Lưới Điện Cao Thế Điều hành viên trạm biến đến HVTH: Lâm Hữu Dũng TPHCM ii áp 220KV GVHD: PGS.TS Lê Chí Kiên Luận Văn LỜI CAM ĐOAN Tôi cam đoan công trình nghiên cứu tơi Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình khác Tp Hồ Chí Minh, ngày 04 tháng 02 năm 2018 (Ký tên ghi rõ họ tên) Lâm Hữu Dũng HVTH: Lâm Hữu Dũng iii GVHD: PGS.TS Lê Chí Kiên Luận Văn LỜI CẢM ƠN Sau thời gian học tập nghiên cứu hướng dẫn tận tình thầy PGS.TS Lê Chí Kiên, học viên hoàn thành luận văn thời hạn mục tiêu đề ban đầu Học viên xin trân trọng bày tỏ lòng biết ơn sâu sắc đến thầy PGS.TS Lê Chí Kiên kiến thức q báu phương pháp nghiên cứu mà thầy truyền đạt Học viên xin gửi lịng tri ân đến q thầy trường Đại Học Sư Phạm Kỹ Thuật Tp Hồ Chí Minh trang bị cho học viên lượng kiến thức bổ ích, đặc biệt xin chân thành cảm ơn q thầy Khoa Điện – Điện Tử tạo điều kiện thuận lợi hỗ trợ cho học viên nhiều trình học tập thời gian làm luận văn Cuối học viên xin gửi lời cám ơn đến đồng nghiệp, gia đình, bạn bè giúp đỡ học viên hỗ trợ học viên mặt tinh thần để học viên hồn thành luận văn Tp Hồ Chí Minh, ngày 04 tháng 02 năm 2018 Học viên thực Lâm Hữu Dũng HVTH: Lâm Hữu Dũng iv GVHD: PGS.TS Lê Chí Kiên Luận Văn TĨM TẮT Điện sản xuất từ nhà máy điện sử dụng nguồn nhiên liệu hóa thạch (than, dầu, khí đốt…) lượng mới: lượng tái tạo (gió, mặt trời, địa nhiệt, sinh khối,…), lượng nguyên tử Các nhà máy điện tuabin khí, tuabin sử dụng nhiên liệu hóa thạch loại sử dụng chủ yếu giới với nhược điểm: ô nhiễm môi trường, giá nhiên liệu cao dần cạn kiệt Bởi vậy, nguồn lượng tái tạo xu hướng thiết thực gần vô tận xem giải pháp tốt thời điểm Trong địa nhiệt nguồn lượng vơ tận chịu ảnh hưởng thời tiết người sử dụng từ lâu chủ yếu sưởi ấm Với cơng nghệ người sử dụng nguồn nhiệt từ giếng địa nhiệt để tạo điện cách sử dụng nóng làm quay tuabin máy phát Nhưng hiệu suất máy phát điện khơng cao bị tổn thất thành phần khí chuyển động Cơng nghệ phát điện từ thủy động hình thành từ lâu, nghiên cứu, thử nghiệm đưa vào sử dụng quốc gia phát triển giới Nó có ưu điểm phát điện trực tiếp từ lượng nhiệt plasma (khí ion hóa) kim loại lỏng có tính dẫn điện cao Nak,…khi xun qua từ trường mạnh mà không cần đến thành phần khí chuyển động Đây phương pháp có hiệu cao độc phát điệ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 “Nghiên cứu hệ thống phát điện từ thủy động lực sử dụng lượng địa nhiệt với chu trình kết hợp” đề xuất mơ hình hệ thống phát điện LMMHD kết hợp lượng địa nhiệt tạo hiệu suất phát điện cao mơ hình địa nhiệt điện truyền thống sử dụng tuabin Kết mơ hình đề xuất có hiệu suất phát điện 34,81% cao hệ thống địa nhiệt truyền thống 9,01% hệ thống địa nhiệt tuabin cấp 6,79% HVTH: Lâm Hữu Dũng v GVHD: PGS.TS Lê Chí Kiên Luận Văn ABSTRACT Electrical power is produced from power plants which uses fossil fuels (coal, oil, gas etc.) or new energies: renewable energy (wind, solar, geothermal, biomass,…), nuclear energy etc Gas turbine, stream turbine power plants which use fossil fuels, are very popular and they have some disadvantages: CO2 emissions, high fuel cost and exhaustion etc Therefore, renewable energy sources are practical trend and are inexhaustible energy source is being seen as the best solution at this point In which geothermal energy is inexhaustible less affected by the weather and it was long used mainly heating Now people can use heat from geothermal wells to generate electricity by using steam that turns a turbine's generator But the efficiency of the turbine generator is not higher due to losses by the moving mechanical components Magnetohydrodynamic technology has formed a long time ago, has been studied, tested and put into use in developing countries around the world The electricity is directly extracted from thermal energy of plasma (ionized gas) or Highly conductive liquid metals such as Nak which is passing through the strong magnetic field without resorting to the mechanical motion components This is a highly effective method and unique for power generation Topic: “Research magnetohydrodynamic (MHD) generator system using geothermal energy with combined cycle” will create higher power generation efficiency models traditional geothermal plants use turbines The result is the performance of the proposed model is 34,81% higher than traditional geothermal systems is 9,01% and geothermal turbine level system is 6.79% HVTH: Lâm Hữu Dũng vi GVHD: PGS.TS Lê Chí Kiên Luận Văn Ta có mơ hình địa nhiệt – MHD kết hợp tuabin nước cơng bố [4] sau: Hình 4.5 Mơ hình địa nhiệt – MHD kết hợp tubin nước công bố Hiệu suất chu trình hình 4.5 cơng bố là: ηMHD-TUABIN = 0,1436 hay 14,36% Ta thấy hiệu suất chu trình địa nhiệt điện – LMMHD kết hợp tua bin nước luận văn (26,18%) cao hiệu suất chu trình Hình 4.5 (14,36%) 11,82% HVTH: Lâm Hữu Dũng 66 GVHD: PGS.TS Lê Chí Kiên Luận Văn Chương KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN CỦA ĐỀ TÀI 5.1 Kết luận Sự phát triển kinh tế, dân số tăng nhanh tốc độ đô thị hóa,… yếu tố ảnh hưởng mạnh mẽ đến nhu cầu lượng toàn giới Việt nam không ngoại lệ, bối cảnh nhu cầu lượng ngày gia tăng, khả cung cấp nguồn lượng nội địa hạn chế Trong đó, tiềm nguồn lượng tái tạo Việt Nam lớn, kèm theo nhu cầu sử dụng điện nhiệt cho sản xuất cao việc xem xét khai thác nguồn lượng tái tạo sẳn có cho sản xuất điện khả công nghệ lẫn hiệu kinh tế môi trường Với mơ hình kết hợp hệ thống phát điện từ thủy động sử dụng nguồn lượng địa nhiệt tạo hiệu suất phát điện cao mô hình truyền thống sử dụng tuabin Đây mơ hình phát điện vừa hiệu vừa không ô nhiễm môi trường phù hợp với xu hướng phát triển tương lai Luận văn “Nghiên cứu hệ thống phát điện từ thủy động lực (MHD) sử dụng lượng địa nhiệt với chu trình kết hợp” dựa vào định luật, phương trình nhiệt động lực học để tính tốn mô thông số hệ thống phát điện MHD kết hợp với lượng địa nhiệt dựa phần mềm Matlab Ta có kết quả: - Tính tốn hiệu suất phát điện chu trình địa nhiệt truyền thống 25,8% - Tính tốn hiệu suất phát điện chu trình địa nhiệt tuabin cấp 28,02% - Tính tốn hiệu suất phát điện chu trình kết hợp địa nhiệt – LMMHD 34,81% - Tính tốn hiệu suất phát điện chu trình kết hợp địa nhiệt – LMMHD kết hợp tuabin nước 26,18% HVTH: Lâm Hữu Dũng 67 GVHD: PGS.TS Lê Chí Kiên Luận Văn - So sánh hiệu suất phát điện chu trình kết hợp địa nhiệt – LMMHD kết hợp tuabin nước (26,18%) với hiệu suất chu trình kết hợp địa nhiệt – MHD kết hợp tuabin nước (14,36%) công bố Với kết nghiên cứu ta thấy rằng: - Hiệu suất hệ thống kết hợp địa nhiệt – LMMHD cao hệ thống địa nhiệt truyền thống 9,01%, hệ thống địa nhiệt tuabin cấp 6,79% Hiệu suất nâng cao có tham gia hệ thống phát điện MHD - Hiệu suất hệ thống địa nhiệt – LMMHD kết hợp tua bin nước luận văn cao hệ thống công bố [4] 11,82%, yếu tố:  Chọn thông số áp suất, nhiệt độ vào tối ưu tính tốn chu trình tua bin nước nên hiệu suất chu trình địa nhiệt truyền thống nâng cao  Chọn tỷ số nhiệt vào máy nén hợp lý làm tăng lưu lượng lưu chất chạy qua máy phát MHD làm nâng cao hiệu suất hệ thống phát điện LMMHD Sự kết hợp yếu tố làm cho hiệu suất chung hệ thống nâng cao hiệu suất hệ thống công bố [4] - Công việc xây dựng mơ hình tốn học tính tốn khơng q phức tạp Tuy nhiên luận văn số khuyết điểm sau: - Việc nghiên cứu dựa cơng thức tốn học tính tốn mơ chưa có điều kiện thực nghiệm - Việc nghiên cứu hướng tới hiệu suất hệ thống chưa sâu vào vấn đề hiệu kinh tế - Hiệu suất chu trình phụ thuộc vào nhiệt lượng cung cấp từ nguồn địa nhiệt tự nhiên nên áp suất nhiệt độ không ổn định 5.2 Hướng phát triển đề tài - Phân tích thêm thành phần ảnh hưởng đến hiệu suất chu trình làm việc - Nghiên cứu thêm phương pháp để tối ưu hóa thơng số nhằm nâng cao hiệu suất phát điện hệ thống MHD kết hợp với địa nhiệt - Tính tốn mơ thơng số chu trình kết hợp với nguồn lượng khác HVTH: Lâm Hữu Dũng 68 GVHD: PGS.TS Lê Chí Kiên Luận Văn TÀI LIỆU THAM KHẢO [1] Đỗ Văn Chương, “Chuyên đề lượng – VNGG”, chương 9, https://sites.google.com/site/vnggenergy/dianhiet, – – 2016 [2] Trần Huyên, “Năng lượng lượng địa nhiệt Việt Nam”, Hội dầu khí Việt Nam, http://petrotimes.vn/nang-luong-va-nang-luong-dia-nhiet-o-viet-namky-1-256111.html, – – 2016 [3] Đoàn Văn Tiến, Đinh Văn Toàn, Trịnh Việt Bắc, “Nghiên cứu nguồn địa nhiệt cho phát triển lượng Việt Nam”, Tạp chí khoa học trái đất, 2008 [4] Nguyễn Bá Sang, luận văn thạc sĩ “Nghiên cứu hệ thống phát điện MHD kết hợp với địa nhiệt điện”, 2013 [5] A Aali a, N Pourmahmoud b, V Zare “Exergoeconomic analysis and multi-objective optimization of a novel combined flash-binary cycle for Sabalan geothermal power plant in Iran”, Faculty of Mechanical Engineering, Urmia University of Technology, Urmia, Iran, 2017, Mechanical Engineering Department, Faculty of Engineering, Urmia University, Urmia, Iran, 2017 [6] Published by the World Energy Council 2017 “World Energy Issues Monitor 2017 -Full-report” [7] IRENA (2017), “Planning for the Renewable Future: Long-term modelling and tools to expand variable renewable power in emerging economies”, International Renewable Energy Agency, Abu Dhabi [8] World Energy Council, World Energy Reso Urces 2016 “World Energy Resources Geothermal 2016”, [9] Geothermal Energy Association, “2016 Annual U.S & Global Geothermal Power Production Report” [10] Sudhir Patel and Gangadharaiah Y H, “Review Note on Magnetohydrodynimics Power Generator”, Department of Mathematics, New Horizon College of Engineering, Bangalore, India, International Journal of Trend in Research and Development, Volume 3(1), ISSN: 2394-9333 www.ijtrd.com, 2016 HVTH: Lâm Hữu Dũng 69 GVHD: PGS.TS Lê Chí Kiên Luận Văn [11] Sufi a Khalili, Ali Jafarian Dehkordi*, Mohammad Hossein Giahi, “Investigating the effect of channel angle of a subsonic MHD (Magneto-HydroDynamic) generator on optimum efficiency of a triple combined cycle”, School of Mechanical Engineering, Tarbiat Modares University, Jalal Ale Ahmad Highway, P.O BOX: 14115-143, Tehran, Iran-2015 [12] Ajith Krishnan R, Jinshah B S, Magnetohydrodynamic Power Generation, International Journal of Scientific and Research Publications, Volume 3, Issue 6, June 2013, ISSN 2250-3153 [13] Hyungsul Moon and Sadiq J Zarrouk, “Efficiency of geothermal power plants: a worldwide review”, Deparment of Engineering Science, University of Auckland, New Zealand, 2013 [14] Ronald DiPippo, Ph D, Chancellor Professor “Geothermal Power Plants: Principles, Applications, Case Studies and Environmental Impact”, Emeritus University of Massachusetts Dartmouth North Dartmouth, Massachusetts [15] John W Lund, Geothermal energy, http://www.britannica.com/science/ geothermal-energy, – – 2016 [16] https://nrginitiative.wordpress.com/2014/04/21/ormat-combined-cycleunits-gccu-geothermal-power-plants/ [17] http://www.learnthermo.com/examples/example-problem.php?ch=9& lesson=B&problem=1 [18] Claus Borgnakke, Richard E Sonntag, University of Michigan, Fundamentals of thermodynamics, Don Fowley and Dan Sayre, 2009 [19] Bilal Masood, Malik Husnain Riaz and M Yasir, Integration of Magnetohydrodynamics (MHD) Power Generating Technology with Thermal Power Plants for Efficiency Improvement, World Applied Sciences Journal 32 (7): 1356-1363, 2014, ISSN 1818-4952, © IDOSI Publications, 2014, DOI: 10.5829/idosi.wasj.2014.32.07.503 [20] Reshmi Banerjee, Importance of Magneto Hydro Dynamic Generation, International Journal of Advanced Research in Electrical, Electronics and Instrumentation Engineering, Vol 4, Issue 7, July 2015 HVTH: Lâm Hữu Dũng 70 GVHD: PGS.TS Lê Chí Kiên Luận Văn PHỤ LỤC % CHƯƠNG TRÌNH TÍNH TỐN HIỆU SUẤT CHU TRÌNH ĐỊA NHIỆT TRUYỀN THỐNG clc clear all; %thong so cai dat ban dau P1=0.02;%MPa P2=0.7;%MPa P3=P2; h3=3375.1;%KJ/Kg s3=7.7884;%KJ/Kg.K h1=251.4;%KJ/Kg v1=0.0010171;%m^3/Kg Wp=v1*(P2-P1) h2=h1+Wp Qh=h3-h2 gamaSteam=0.8320; deltaTuabin=7.0752; x4=(s3-gamaSteam)/deltaTuabin; gamaCondenser=2357.5; h4=h1+x4*gamaCondenser %nang luong sinh WT=h3-h4 Wnet=WT-Wp %hieu suat he thong nuyHT=Wnet/Qh HVTH: Lâm Hữu Dũng 71 GVHD: PGS.TS Lê Chí Kiên Luận Văn %nhiet cac nut T3=723; toTuabin=0.68; T4=T3*toTuabin Nbom=1.07; T1=333.21 T2=T1*Nbom %do thi Ts x=[h3 h4 h1 h2 h3] y=[T3 T4 T1 T2 T3] h=plot(x,y,'b-*') grid xlabel('Truc Entropy') ylabel('Truc nhiet (K)') % CHƯƠNG TRÌNH TÍNH TỐN HIỆU SUẤT CHU TRÌNH ĐỊA NHIỆT SỬ DỤNG TUABIN CẤP clc clear all %thong so cai dat ban dau P1=0.02;%MPa P6=P1 P2=0.7;%MPa P3=P2; P4=0.3;%MPa P5=P4 h3=3375.1;%KJ/Kg s3=7.7884;%KJ/Kg.K h1=251.4;%KJ/Kg HVTH: Lâm Hữu Dũng 72 GVHD: PGS.TS Lê Chí Kiên Luận Văn v1=0.0010171;%m^3/Kg h4s=561.46; gamaSteam1=1.9921; deltaTuabin1=4.7149; gamaCondenser1=2065.56; h5=3380.2;%KJ/Kg s5=8.1848;%KJ/Kg.K x4=(s3-gamaSteam1)/deltaTuabin1 h4=h4s+x4*gamaCondenser1 Whp=h3-h4 gamaSteam2=0.832; deltaTuabin2=7.0752; x6=(s5-gamaSteam2)/deltaTuabin2 gamaCondenser2=2357.2; h6=h1+x6*gamaCondenser2 Wlp=h5-h6 Wt=(h3-h4)+(h5-h6) Wp=v1*(P2-P1) h2=h1+Wp Qh=(h3-h2)+(h5-h4) %nang luong sinh Wnet=Wt-Wp %hieu suat he thong nuyHT=Wnet/Qh %nhiet cac nut T3=723; toTuabin=0.68; T4=T3*toTuabin T5=T3 HVTH: Lâm Hữu Dũng 73 GVHD: PGS.TS Lê Chí Kiên Luận Văn Nbom=1.07; toCondenser=0.56; T1=333.21; T6=T5*toCondenser T2=T1*Nbom %do thi Ts x=[h3 h4 h5 h6 h1 h2 h3] y=[T3 T4 T5 T6 T1 T2 T3] h=plot(x,y,'b-*') grid xlabel('Truc Entropy') ylabel('Truc nhiet (K)') % CHƯƠNG TRÌNH TÍNH TỐN HIỆU SUẤT CHU TRÌNH KẾT HỢP ĐỊA NHIỆT – LMMHD clc clear all y=1.6667;%he so dan nhiet cua kim loai long Qin=3.4e+003;%nhiet luong he thong dia nhiet cung cap P2=4.0e+005; Ttraodoinhiet=723;%nhiet bo trao doi nhiet T2=Ttraodoinhiet Ptron=0.2e+004; ntach=0.9; nEE=0.4; nMHD=0.8; deltaQMHD=0.005; deltaQngung=0.01; nuyTngung=1; pinen=2;%ty so nen HVTH: Lâm Hữu Dũng 74 GVHD: PGS.TS Lê Chí Kiên Luận Văn tonen=1.05;%ty so nhiet vao may nen Cp=0.27;%nhiet dung rieng NaK deltaTtron=5; deltaTtach=5; %nut T3=T2-deltaTtron P3=Ptron %nut T4=T3*(1-nEE-deltaQMHD) P4=P3*((1-(nEE/nMHD))^(y/(y-1))) %nut P5=P4/2 P5a=P5 T5=T4-deltaTtach %nut Tout=700; Tin=660; P1=P5*(1-deltaQngung*nuyTngung) T1=T5-((Tout-Tin)/(1-deltaQngung)) %nut P2=P1*pinen T2=T1*tonen %luu luong kim loai long G=Qin/(Cp*(T3-T2)) %nhiet luong tai tung nut cua chu trinh Q1=G*Cp*T1 Q2=G*Cp*T2 Q3=G*Cp*T3 Q4=G*Cp*T4 HVTH: Lâm Hữu Dũng 75 GVHD: PGS.TS Lê Chí Kiên Luận Văn Q5=G*Cp*T5 %cong suat may phat MHD WMHD=Q3*nEE %cac cong suat tieu thu chu trinh hesobom=1.25;%thong so cai dat chu trinh Wtron=1.0e+003;%thong so cai dat chu trinh Wtach=0.5e+003;%thong so cai dat chu trinh Wbom=hesobom*(T2-T1)*G*Cp Wbomdt=140;%thong so tra cuu %hieu suat he thong MHD Wnet=WMHD-Wtron-Wtach-Wbom-Wbomdt nuytong=Wnet/Qin %tinh Entropy Tref=298; Pref=1.03e+005; %chu trinh MHD s=(Cp*log(T3/Tref)-Cp*((y-1)/y)*log(P3/Pref))*G s3=s+(Cp*log(T3/T2)-Cp*((y-1)/y)*log(P3/P2))*G s4=s3+(Cp*log(T4/T3)-Cp*((y-1)/y)*log(P4/P3))*G s5=s4+(Cp*log(T5/T4)-Cp*((y-1)/y)*log(P5/P4))*G s1=s5+(Cp*log(T1/T5)-Cp*((y-1)/y)*log(P1/P5))*G s2=s1+(Cp*log(T2/T1)-Cp*((y-1)/y)*log(P2/P1))*G x=[s3 s4 s5 s1 s2 s3]; y=[T3 T4 T5 T1 T2 T3]; h=plot(x,y,'b-*'); grid xlabel('Truc Entropy') ylabel('Truc nhiet (K)') HVTH: Lâm Hữu Dũng 76 GVHD: PGS.TS Lê Chí Kiên Luận Văn % CHƯƠNG TRÌNH TÍNH TỐN HIỆU SUẤT CHU TRÌNH KẾT HỢP ĐỊA NHIỆT – LMMHD VỚI TUABIN HƠI NƯỚC THEO CÁC THÔNG SỐ THAM KHẢO TRONG [4] clc clear all Wnet1=367210.5; Qh=1423527; y=1.6667;%he so dan nhiet cua kim loai long Qin=1.0e+004;%nhiet luong he thong dia nhiet cung cap P2=5.0e+004; Ttraodoinhiet=340;%nhiet bo trao doi nhiet T2=Ttraodoinhiet Ptron=0.2e+004; ntach=0.9; nEE=0.4; nMHD=0.8; deltaQMHD=0.005; deltaQngung=0.01; nuyTngung=1; pinen=2;%ty so nen tonen=1.05;%ty so nhiet vao may nen Cp=0.27;%nhiet dung rieng NaK deltaTtron=5; deltaTtach=5; %nut T3=T2-deltaTtron P3=Ptron %nut T4=T3*(1-nEE-deltaQMHD) HVTH: Lâm Hữu Dũng 77 GVHD: PGS.TS Lê Chí Kiên Luận Văn P4=P3*((1-(nEE/nMHD))^(y/(y-1))) %nut P5=P4/2 P5a=P5 T5=T4-deltaTtach %nut Tout=450; Tin=410; P1=P5*(1-deltaQngung*nuyTngung) T1=T5-((Tout-Tin)/(1-deltaQngung)) %nut P2=P1*pinen T2=T1*tonen %luu luong kim loai long G=Qin/(Cp*(T3-T2)) %nhiet luong tai tung nut cua chu trinh Q1=G*Cp*T1 Q2=G*Cp*T2 Q3=G*Cp*T3 Q4=G*Cp*T4 Q5=G*Cp*T5 %cong suat may phat MHD WMHD=Q3*nEE %cac cong suat tieu thu chu trinh hesobom=1.25;%thong so cai dat chu trinh Wtron=1.0e+003;%thong so cai dat chu trinh Wtach=0.5e+003;%thong so cai dat chu trinh Wbom=hesobom*(T2-T1)*G*Cp Wbomdt=140;%;thong so tra cuu HVTH: Lâm Hữu Dũng 78 GVHD: PGS.TS Lê Chí Kiên Luận Văn %hieu suat he thong MHD Wnet2=WMHD-Wtron-Wtach-Wbom-Wbomdt WdianhietMHD=Wnet1+Wnet2 nuytong=(Wnet1+Wnet2)/Qh %tinh Entropy Tref=298; Pref=1.03e+005; %chu trinh MHD s=(Cp*log(T3/Tref)-Cp*((y-1)/y)*log(P3/Pref))*G s3=s+(Cp*log(T3/T2)-Cp*((y-1)/y)*log(P3/P2))*G s4=s3+(Cp*log(T4/T3)-Cp*((y-1)/y)*log(P4/P3))*G s5=s4+(Cp*log(T5/T4)-Cp*((y-1)/y)*log(P5/P4))*G s1=s5+(Cp*log(T1/T5)-Cp*((y-1)/y)*log(P1/P5))*G s2=s1+(Cp*log(T2/T1)-Cp*((y-1)/y)*log(P2/P1))*G x=[s3 s4 s5 s1 s2 s3]; y=[T3 T4 T5 T1 T2 T3]; h=plot(x,y,'b-*'); grid xlabel('Truc Entropy') ylabel('Truc nhiet (K)') HVTH: Lâm Hữu Dũng 79 GVHD: PGS.TS Lê Chí Kiên S K L 0