Xem mẫu

  1. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) XÂY DỰNG CÔNG CỤ MÔ PHỎNG VỊ TRÍ MẶT TRỜI VÀ CƯỜNG ĐỘ BỨC XẠ TỚI TRÁI ĐẤT THEO THỜI GIAN VÀ VỊ TRÍ ĐỊA LÝ DEVELOPE A PROGRAM TO SIMULATE THE SUN MOTION AND RADIATION RELATED TO DAY'S TIME AND GEOGRAPHIC LOCATION ON THE EARTH Phạm Anh Tuân Trường Đại học Điện lực Ngày nhận bài: 04/01/2019, Ngày chấp nhận đăng: 28/03/2019, Phản biện: TS. Phạm Mạnh Hải Tóm tắt: Các ứng dụng năng lượng mặt trời ngày càng chiếm một vị trí quan trọng đối với cuộc sống, vì đây là nguồn năng lượng sạch, bao phủ khắp Trái Đất và vô tận. Tính đến năm 2050 nhu cầu sử dụng năng lượng từ điện mặt trời sẽ chiếm khoảng 20 % tổng nhu cầu điện năng trên toàn thế giới [1,2]. Năng lượng mặt trời không đồng nhất trên Trái Đất; nó phụ thuộc vào tọa độ, thời gian trong ngày, ngày trong năm, dữ liệu khí tượng... Trong bài toán dự báo, bức xạ mặt trời có thể được tính trên cơ sở là các phương trình mô tả chuyển động tương đối của Trái Đất - Mặt Trời, kết hợp với các thông tin khác như vị trí địa lý, ngày trong năm, sự hấp thụ và phản xạ của các tầng mây... Trong bài báo này, chúng tôi trình bày các phương trình thiên văn học mô tả chuyển động của Trái Đất - Mặt Trời, qua đó lập trình để tính toán và mô phỏng nhằm đưa ra các đồ thị mô tả: toạ độ và góc phương vị của Mặt Trời theo thời gian, năng lượng bức xạ chiếu đến Trái Đất trên nền matlab... Kết quả tính toán và mô phỏng được ứng dụng để tính và đưa ra một số đồ thị kỹ thuật nhằm mô tả chuyển động của Mặt Trời và lượng bức xạ tại các toạ độ khác nhau tại vị trí toạ độ bất kỳ trên Trái Đất và là cơ sở cho việc chọn hướng lắp đặt dàn pin mặt trời hoặc các thiết bị thu năng lượng mặt trời và dự báo tổng năng lượng sẽ thu được. Từ khóa: Bức xạ mặt trời, mô phỏng năng lượng mặt trời, điện mặt trời, Matlab/GUI. Abstract: Solar energy applications increasingly occupy an important position for life, as this is a clean energy, renewable resource and covering all over the earth. By 2050, the demand for solar energy will account for about 20% of the total electricity demand worldwide [1,2]. Solar energy is not uniform on earth; it depends on coordinates, time of day, day of year, meteorological data, etc. For the energy prediction, solar radiation can be calculated on the basis of equations describing relative motion of the earth-sun, combined with other information such as geographic location, day of the year, absorption and reflection of cloud layers, etc. In this paper, we present astronomical equations describing the motion of the earth-sun, thereby programming to calculate and simulate to produce descriptive graphs: coordinates and angles the position of the sun over time, radiant energy to the earth under Matlab/GUI framework, etc. The simulation results are applied to calculate and give a number of technical graphs to describe the sun's movement and the amount of radiation at different coordinates at any position on earth and is the basis for choosing the direction of installing solar panels or solar collectors and forecasting the total energy to be collected. 48 Số 19
  2. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) Keywords: Solar radiation, solar energy simulation, solar electric system, Matlab/GUI. 1. GIỚI THIỆU CHUNG 8766 giờ. Quĩ đạo của Trái Đất quanh Năng lượng mặt trời có thể chia thành 2 Mặt Trời tại các thời điểm trong năm dạng gồm năng lượng bức xạ nhiệt và được mô tả như hình 1 [4-6]. năng lượng ánh sáng. Các ứng dụng của năng lượng mặt trời bao gồm cả nhiệt năng và điện năng ngày càng trở nên quan trọng với cuộc sống. Điểm mạnh năng lượng mặt trời so với nguồn năng lượng khác là miễn phí, sạch và phong phú. Hơn nữa, năng lượng mặt trời có thể khai thác ở hầu hết các nơi điểm trên bề mặt trái đất. Nó càng quan trọng hơn khi chi phí Hình 1. Vị trí tương đối của Trái Đất và Mặt Trời nhiên liệu hóa thạch cao và những tác động làm ô nhiễm môi trường từ việc Đối với một điểm quan sát cụ thể trên đốt chúng. Năng lượng mặt trời được Trái Đất, vị trí mặt trời có thể được xác truyền đến trái đất thông qua bức xạ, định bởi hai góc, bao gồm: góc cao độ (α) bao gồm hai phần: phần bức xạ nằm trên và góc phương vị (θs); các góc này được bầu khí quyển và phần bức xạ dưới bầu mô tả như trong hình 2, trong đó: khí quyển. Phần năng lượng bức xạ dưới bầu khí quyển có thể được đo bằng Góc cao độ α (hay còn gọi là góc nhập xạ các thiết bị đo năng lượng mặt trời như: - Altitude) là góc tạo bởi chùm tia sáng nhật xạ kế (pyranometer, solarimeter), mặt trời so với mặt phẳng xích đạo có thể trực xạ kế (pyrheliometer), nhật xạ ký được xác định bởi phương trình 1: (actinography) [3]. Ngoài ra, nó cũng có sin sin . sin cos . cos . cos thể được tính toán dựa trên các nguyên lý thiên văn học và hình học, thông qua đó (1) sẽ mô hình hóa nguồn mặt trời chiếu tới Với: L là vĩ độ tại vị trí quan sát, δ là góc trái đất. Trong nghiên cứu này chúng tôi nghiêng tại vị trí quan sát (được tính theo sẽ xây dựng công cụ để tính toán mô phương trình 2), và ω là góc giờ tại vị trí phỏng vị trí và cường độ bức xạ năng quan sát (được tính theo phương trình 8). lượng mặt trời theo thời gian và vị trí địa lý trên nền Matlab/GUI.  23,45°. sin (2) 2. CƠ SỞ LÝ THUYẾT VÀ MÔ HÌNH ở đây N là ngày qui đổi trong năm, N được TOÁN HỌC tính qui ước từ ngày 1 tháng 1 (N=1); các Trái Đất xoay quanh Mặt Trời trong một ngày tiếp theo sẽ được tính tròn theo ngày, quỹ đạo hình elip, thời gian của mỗi vòng ví dụ ngày 2 tháng 1 thì N=2, ngày 3 tháng Trái Đất quay quanh Mặt Trời là khoảng 1 thì N=3… Số 19 49
  3. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) Trái Đất nghiêng so với mặt phẳng quĩ 15°. (7) đạo của nó quanh Mặt Trời; phương trình Trong phương trình 7, hệ số 15° chính là thực nghiệm để điều chỉnh độ lệch tâm góc quay được của Trái Đất mỗi giờ (một (EoT) do độ nghiêng trục này được tính như sau (theo đơn vị phút): ngày Trái Đất quay được 360° trong vòng 24 giờ); LMST được hiểu là thời gian qui đổi theo phút của múi giờ (giờ GMT). 9,87. sin 2 7,53. cos 1,5. sin Trên cơ sở đó có thể tính được góc giờ tại (3) vị trí quan sát: ớ 81 (4) 15°. 12 (8) Một số nhận xét: vào ngày 21/3 hoặc 4°. 23/9, Mặt Trời sẽ lên thiên đỉnh tại xích (9) đạo và vào ngày 22/6 hoặc ngày 22/12, Mặt Trời lên thiên đỉnh tại chí tuyến Bắc Theo qui ước này, vào buổi sáng góc giờ và chí tuyến Nam. sẽ mang giá trị qui ước là âm và vào buổi chiều giờ góc sẽ mang giá trị qui ước là Góc phương vị (Azimuth) là góc giữa dương. vector bắc và mặt phẳng chiếu vuông góc của Mặt Trời xuống mặt phẳng xích đạo; Năng lượng ánh sáng và cường độ bức xạ góc phương vị có thể được tính theo công phát ra từ Mặt Trời sẽ tỉ lệ với bình thức: phương khoảng cách tới điểm nhận bức   xạ. Cường độ bức xạ chiếu tới Trái Đất sẽ sin  (5)  tỉ lệ với bình phương khoảng cách thực Hệ số hiệu chỉnh thời gian (tính bằng giữa Trái Đất và Mặt Trời (Gex) được tính phút) tại vị trí khảo sát trong múi giờ (giờ bởi phương trình sau [5]: GMT) của vị trí đó khi xét đến ảnh hưởng của độ lệch tâm EoT được tính theo . (10) phương trình 6 (với LOD là kinh độ của trong ®ã Rav lμ kho¶ng c¸ch trung b×nh vị trí khảo sát): gi÷a MÆt Trêi vμ Tr¸i §Êt vμ R lμ kho¶ng LMT = 4.(LOD – LSMT) + EoT (6) c¸ch thùc gi÷a MÆt Trêi vμ Tr¸i §Êt øng víi thêi ®iÓm cÇn kh¶o s¸t. Theo c¸c sè liÖu tÝnh to¸n thèng kª th× tØ lÖ kho¶ng c¸ch thùc so víi kho¶ng c¸ch trung b×nh sÏ phô thuéc vμo ngμy qui ®æi trong n¨m vμ cã gi¸ trÞ kho¶ng: 1 0,0333. cos (11) Thay thế phương trình 11 vào phương trình 10 ta có: Hình 2. Mô tả góc cao độ và góc phương vị của Mặt Trời . 1 0,0333. cos (12) 50 Số 19
  4. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) Gex chÝnh lμ cường ®é bøc x¹ n¨ng lượng ViÖc tÝnh to¸n phÇn n¨ng lượng t¸n x¹ mÆt trêi chiÕu xuèng Tr¸i §Êt. thùc tÕ rÊt khã so víi tÝnh to¸n bøc x¹ mÆt trêi trùc tiÕp. PhÇn t¸n x¹ phô thuéc vμo Cường ®é bøc x¹ cña MÆt Trêi trªn mÆt líp khÝ quyÓn vμ tÇng m©y. Mét phÇn ®Êt (GT) được hiÓu lμ bøc x¹ cña MÆt Trêi n¨ng lượng chiÕu xuèng Tr¸i §Êt sÏ bÞ xuyªn qua bÇu khÝ quyÓn xuèng tíi ®é cao ph¶n x¹ tõ bÒ mÆt trë l¹i bÇu trêi vμ mét mùc nước biÓn. Bøc x¹ mÆt trêi trªn mÆt phÇn t¸n x¹ xuèng mÆt ®Êt. Theo m« h×nh ®Êt bao gåm hai thμnh phÇn: phÇn chiÕu ®¼ng hướng (coi lượng t¸n x¹ ®Õn mét ®Þa xuèng trùc tiÕp (GB) vμ phÇn t¸n x¹ bëi ®iÓm cã cường ®é b»ng nhau tõ mäi c¸c ®¸m m©y vμ tÇng khÝ quyÓn (GD). hướng), n¨ng lượng t¸n x¹ cã thÓ xÊp xØ (GD).GT=GB+GD (13) b»ng:  360  N 100 GB (18) Khi chïm bøc x¹ mÆt trêi ngoμi Tr¸i §Êt GD  0,095  0,04sin  ®i qua bÇu khÝ quyÓn, mét phÇn cña chïm  365  tia nμy bÞ hÊp thô bëi bÇu khÝ quyÓn. Th«ng thường ®èi víi mét ngμy trêi 3. XÂY DỰNG CÔNG CỤ MÔ PHỎNG VỊ quang, 70% bøc x¹ mÆt trêi xuèng Tr¸i TRÍ MẶT TRỜI VÀ CƯỜNG ĐỘ SÁNG §Êt lμ bøc x¹ mÆt trêi trùc tiÕp. Theo m« TỚI TRÁI ĐẤT THEO THỜI GIAN VÀ VỊ TRÍ ĐỊA LÝ h×nh ASHRAE, bøc x¹ mÆt trêi trùc tiÕp tíi bÒ mÆt tr¸i ®Êt (GB, nom) cã thÓ được Trên cơ sở lý thuyết đã trình bày ở trên, tÝnh như sau: chúng tôi xây dựng thuật toán và lập trình K xây dựng công cụ mô phỏng cường độ GB ,nom  Ae sin  (14) bức xạ năng lượng mặt trời theo thời gian và vị trí địa lý trên Matlab/GUI. Công cụ Trong đó A là một thông lượng biểu kiến này cho phép tính toán và mô phỏng vị trí chiếu tới bề mặt trái đất và K là độ sâu tương đối của Trái Đất và Mặt Trời, quang học (đại lượng này không thứ cường độ bức xạ của mặt trời trên bề mặt nguyên). A và K phụ thuộc vào ngày qui trái đất... Sơ đồ khối của thuật toán này đổi trong năm và có thể được tính theo được mô tả trong hình 3. công thức thực nghiệm như sau: Hình 4 là giao diện chính của công cụ mô  360 A  1160  75sin   N  275   (15) phỏng vị trí và cường độ bức xạ của Mặt  365  Trời tới Trái Đất.  360 K  0,174  0, 035sin   N  100   (16) D÷ liÖu ®Çu vμo gåm: täa ®é ®Þa lý, ngμy  365  th¸ng, mói giê, sè ®iÓm tÝnh trªn ®å thÞ... Ngoμi ra ®Ó thuËn tiÖn trong sö dông, Theo đó, cường độ bức xạ mặt trời trực c«ng cô nμy cho phÐp người dïng chän tiếp trên độ cao mực nước biển GB có thể nhËp vÞ trÝ cÇn kh¶o s¸t theo mét sè ®Þa được tính theo công thức 17. danh được t¹o s½n trªn nÒn file excel, tiÕp GB=GB,nomsin (17) ®ã d÷ liÖu sÏ tù ®éng cËp nhËp c¸c th«ng Số 19 51
  5. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) tin vÒ kinh ®é vμ vÜ ®é, mói giê cña ®iÓm tÝnh được cμi ®Æt, sau ®ã khi tÝnh to¸n cÇn kh¶o s¸t. c«ng cô cho phÐp xuÊt ®å thÞ m« pháng (cöa sæ bªn tr¸i giao diÖn chÝnh) hoÆc xuÊt D÷ liÖu ®Çu ra gåm: gãc cao ®é, gãc d÷ liÖu dưới d¹ng sè. phương vÞ, cường ®é bøc x¹ mÆt trêi t¹i c¸c ®iÓm trªn bÒ mÆt khÝ quyÓn, t¹i c¸c Người dùng có thể khảo sát các dữ liệu ®iÓm ®o trªn Tr¸i §Êt. D÷ liÖu nμy được đầu ra hoặc dựng đồ thị trên cùng một tÝnh theo tõng ngμy víi kho¶ng thêi gian trục trong nhiều ngày. Hình 3. Sơ đồ khối thuật toán mô phỏng vị trí của Mặt Trời và bức xạ tới Trái Đất theo thời gian Hình 4. Giao diện của phần mềm mô phỏng vị trí của Mặt Trời và bức xạ tới Trái Đất theo thời gian 52 Số 19
  6. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) 4. KẾT QUẢ Quan s¸t trªn ®å thÞ thÊy gãc cao ®é lín nhÊt (®iÓm Al-2) gÇn trïng víi gãc phương vÞ b»ng 0 (®iÓm Az-2, MÆt Trêi trïng víi phương nh×n hướng vÒ phÝa b¾c) vμo thêi ®iÓm kho¶ng 11 giê 45 phót trong ngμy (tøc lμ vμo kho¶ng gi¸ trÞ 705 phót qui ®æi øng víi trôc tung cña ®å thÞ). Ngoμi ra đường phương vÞ cũng cho thÊy thêi ®iÓm r¹ng s¸ng (®iÓm Al-, kho¶ng thêi gian 380 phót quy ®æi - gãc cao ®é vượt qua 0) còng lμ thêi ®iÓm gãc phương vÞ øng víi gi¸ trÞ trªn trôc tung vμo kho¶ng -65o (®iÓm Az-1); ®iÒu nμy cho biÕt MÆt Trêi kh«ng ë hướng chÝnh ®«ng Hình 5A. mô phỏng góc cao độ (Alitude angle) và góc phương vị (Azimuth angle) tại thành mμ ë gãc kho¶ng 65o ®«ng nam. Gãc cao phố Hà Nội ngày 1 tháng 2 ®é còng cho thÊy MÆt Trêi lÆn (®iÓm A1-3, gãc cao ®é tiÕn tíi 0) vμo kho¶ng Mét kÕt qu¶ m« pháng gãc cao ®é vμ gãc 1035 thêi ®iÓm phót quy ®æi, tøc kho¶ng phương vÞ cña thμnh phè Hμ Néi (vÜ ®é 17 giê 15 phót. H×nh 5B minh häa vÞ trÝ 21o 02’ 00” b¾c, kinh ®é 105o 51’ 00” cña MÆt Trêi nh×n tõ ®iÓm kh¶o s¸t. ®«ng) ngμy 1 th¸ng 12 như trong h×nh 5A; trong ®ã trôc hoμnh lμ trôc thêi gian trong Một kết quả mô phỏng để so sánh về ngμy qui ®æi theo phót (0 giê 00 phót sÏ cường độ bức xạ của Mặt Trời tại mặt đất tương øng víi gi¸ trÞ phót 0), trôc tung lμ cho khu vực thành phố Hà Nội và khu vực gi¸ trÞ gãc, tÝnh theo ®é (tÝnh chung cho c¶ thành phố Hồ Chí Minh (vĩ độ 10° 50’ gãc cao ®é vμ gãc phương vÞ). 00” bắc, kinh độ 106° 37’ 58” đông) vào ngày 1 tháng 5 như hình 6. Kết quả này được thực hiện mô phỏng trong điều kiện bỏ qua sự suy giảm bức xạ do hấp thụ và phản xạ của mây. Các kết quả trên hình thể hiện mặt trời xuất hiện vào khoảng 6 giờ sáng (ứng với điểm qui đổi khoảng 360 phút), và lặn vào khoảng 18 giờ (ứng với điểm qui đổi khoảng 1080 phút) và cường độ bức xạ lớn nhất tại thời điểm khoảng 12 h, trưa khi đó cường độ bức xạ vào khoảng 1050 W/m2 đối với Hà Nội và Hình 5B. Minh họa vị trí của Mặt Trời 1100 W/m2 đối với thành phố Hồ Chí tại thành phố Hà Nội ngày 1 tháng 12 Minh). Số 19 53
  7. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) Một kết quả khác tính toán mô phỏng bức xạ mặt trời trong ngày mồng 1 tại các tháng khác nhau từ tháng 1 đến tháng 6 tại thành phố Hà Nội như hình 8 (các đường được đánh số từ 1 đến 6 ứng với các tháng). Kết quả cho thấy cường độ bức xạ thay đổi rõ rệt từ mức cường độ bức xạ khoảng 740 W/m2 trong tháng 1 đã tăng lên đến khoảng 1130 W/m2 trong tháng 6. Ngoài ra thông qua các đường Hình 6. Mô phỏng so sánh cường độ sáng bức xạ này ta cũng có thể dễ dàng quan tại Hà Nội và Thành phố Hồ Chí Minh sát thấy thời gian xuất hiện cường độ bức xạ (ban ngày) của các tháng mùa hè tăng lên nhiều so với các tháng mùa đông. Hình 7. Mô phỏng bức xạ mặt trời ngày 18 trên tại bề mặt khí quyển và Trái Đất tại Hà Nội Hình 8. Mô phỏng bức xạ mặt trời KÕt qu¶ tÝnh to¸n m« pháng kÕt qu¶ so trong các ngày mồng 1 của các tháng 1 s¸nh cường ®é bøc x¹ t¹i bÒ mÆt khÝ đến tháng 6 tại Hà Nội quyÓn vμ bÒ mÆt tr¸i ®Êt t¹i Hμ Néi ngμy 1 5. KẾT LUẬN th¸ng 8 như trong h×nh 7. M« pháng nμy cho thÊy t¹i thêi ®iÓm cường ®é bøc x¹ Các thiết bị đo bức xạ mặt trời như lín nhÊt vμo kho¶ng 12 giê. Đường phÝa nhật xạ kế (pyranometer, solarimeter), trªn (đường xanh) lμ cường ®é bøc x¹ t¹i trực xạ kế (pyrheliometer), nhật xạ ký bÒ mÆt khÝ quyÓn; gi¸ trÞ lín nhÊt vμo (actinography) thường được lắp đặt tại kho¶ng 1330 W/m2. Đường phÝa dưới các địa điểm được lựa chọn ở các vùng cụ (đường mμu vμng) lμ cường ®é bøc x¹ ë thể. Trên thực tế, để thực hiện đo và tổng bÒ mÆt tr¸i ®Êt, gi¸ trÞ lín nhÊt chØ cßn hợp dữ liệu đo này có chi phí khá lớn nên kho¶ng 1030 W/m2 (do bÞ suy gi¶m bëi không thể đặt chúng tại mọi điểm trên khÝ quyÓn). Trái Đất. Các tính toán này cũng được 54 Số 19
  8. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) dùng để dự báo năng lượng mặt trời bức nhiều ngày; mô phỏng này giúp người xạ ở mọi địa điểm trên Trái Đất. Nó có ý dùng dễ dàng so sánh, phân tích các kết nghĩa đặc biệt trong việc xác định dữ liệu quả. Kết quả này là cơ sở cho việc tính tại những vị trí không có hoặc chưa thể chọn hướng lắp đặt các thiết bị khai thác lắp đặt được thiết bị đo năng lượng mặt năng lượng mặt trời (dàn pin mặt trời, trời. thiết bị thu nhiệt…) ở một vị trí bất kỳ để đảm bảo năng lượng bức xạ chiếu tới nó Các giá trị mà công cụ mô phỏng này tính là lớn nhất. Ngoài ra, các kết quả mô được mặc dù không mới, tuy nhiên kết phỏng cũng là cơ sở cho việc tích phân để quả mô phỏng giúp người dùng nhanh dự báo tổng năng lượng chiếu đến một chóng có được đồ thị mô phỏng trong điểm bất kỳ. TÀI LIỆU THAM KHẢO [1] https://www.equinor.com/en/how-and-why/sustainability/energy-perspectives.html, 2018. [2] IRENA, Global Energy transformation: A roadmap to 2050 (2019 edition). International Renewable Energy Agency, 2019. [3] S. Radiation and D. Models, Solar Radiation and Daylight Models. Elsevier, 2004. [4] V. Badescu, Modeling Solar Radiation at the Earth’s Surface. Spinger, 2008. [5] T. Khatib, Modeling of photovoltaic systems using MATLAB. Wiley, 2016. [6] “https://www.pveducation.org/pvcdrom/terrestrial-solar-radiation,” 2018. Giới thiệu tác giả: Tác giả Phạm Anh Tuân nhận bằng Thạc sĩ ngành kỹ thuật điện năm 2006, nhận bằng Tiến sĩ ngành khoa học vật liệu năm 2017 tại Trường Đại học Bách khoa Hà Nội. Tác giả hiện là giảng viên Khoa Kỹ thuật điện, Trường Đại học Điện lực. Lĩnh vực nghiên cứu: vật liệu pin mặt trời và mô phỏng hệ thống điện mặt trời. Số 19 55
  9. TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ NĂNG LƯỢNG - TRƯỜNG ĐẠI HỌC ĐIỆN LỰC (ISSN: 1859 - 4557) 46 Số 19
nguon tai.lieu . vn