太陽能發電為目前國際上積極發展之重要綠色能源之一，各地設置於建築物屋頂或地面的太陽光電系統日漸增多。目前常見的太陽能光電系統大致分成五類，分別為斜屋頂平貼型、屋頂距置型、地面距置型、地面單斜式棚架型與平屋頂單斜式棚架型。依據國內相關規範，當地面單斜式棚架型之高寬比(棚架的屋頂平均高度除以屋頂弦長水平投影)介於0.25到1之間時，風壓係數為定值。而在此參數範圍外，規範未規定其風壓係數。因此，本研究針對地面單斜式太陽能棚架，使用計算流體力學方法模擬其平均受力特性；紊流模型採用穩態SST k-ω 模型，來流之上游假設為開闊地況；將依據風洞實驗結果來驗證模擬的準確性。再分別模擬棚架屋頂傾角10度、30度，以及不同高寬比之棚架，比較其平均風壓係數的差異性，作為地面單斜式棚架型耐風設計上之參考。
本研究發現網格加密後，屋頂大部分平均風壓係數差異不大，但在屋頂邊緣會產生較平滑且合理的趨勢。當屋頂傾角較大時，CFD結果與風洞實驗結果較為接近，判斷是由於此時屋頂後方產生較大尺度雙渦流，SST k-ω 模型會有較準確之結果，同時風洞實驗屋頂下柱子的影響不大。當高寬比較小時，高寬比的變化會造成屋頂最大正風壓係數及最大負風壓係數較明顯的變化。

Solar power generation is one of the important green energy sources currently actively developed internationally. Solar photovoltaic systems installed on the ground or building roofs are increasing. Ac-cording to the wind resistant design code, when the height-to-width ra-tio (average height of the roof of the canopy divided by the horizontal projection of the chord length of the roof) of the ground monoslope canopy is between 0.25 and 1, the wind pressure coefficient is given as a constant; no values are specified outside this parameter range. There-fore, this study aims at exploring the characteristics of mean wind pres-sures of ground monoslope canopies using CFD; steady-state SST k-ω turbulence model is adopted and the upstream of the incoming flow is assumed to be open terrain; the associated wind tunnel test results will be used to verify the accuracy of the CFD simulation. It is found that finer grid arrangements give smoother and more reasonable prediction for the wind pressures, especially at the roof edge. When the roof in-clined angle gets larger, larger-scale eddies behind the roof make the SST k-ω model more accurate, and the pillars under the roof in wind tunnel tests have little effect. Moreover, the variation of the height-to-width ratio has a larger impact on the extremes of mean wind pressures when the height-to-width ratio is smaller.

[1] ANSYS,Inc. “ANSYS FLUENT 14.0,User’s Guide”, 2011.
[2] ANSYS,Inc. “ANSYS FLUENT 12.0 Theory Guide”, 2009.
[3] ASCE7. “Minimum Design Loads And Associated Criteria For Buildings And Other Structures”, 2016.
[4] Ayodeji Abiola-Ogedengbe. “Experimental investigation of wind effect on solar panels”, 2013.
[5] A. M. Aly., G. Bitsuamlak., “Aerodynamic loads on ground-mounted solar panels: multi-scale computational and experimental investigations”, 2014.
[6] Blocken, B., Stathopoulo, T., Carmeliet, J. “CFD simulation of the atmospheric boundary layer: wall function problems”,Atmos. Environ, 2007, pp. 238-252.
[7] Chowdhury Mohammad Jubayer, Horia Hangan. “Numerical simulation of wind effects on a stand-alone ground mounted pho-tovoltaic (PV) system”, 2014.
[8] Duthinh. D., J. A. Main, B. M. Phillips. “Analysis of Wind Pressure Data on Components and Cladding of Low-Rise Buildings” ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems, 2017, Part A: Civil Engineering, ASCE, ISSN 2376-7642.
[9] Delaunay, B.“Sur La sphère Vide.” Izv. Akad. Nauk SSSR, Otdelenie Matematicheskii I Estestvennyka Nauk, 1934, 7(793-800):1-2.
[10] Delaunay, B. “On the empty sphere. ” In Memory of Georges Voronoi. Bulletin of the USSR Academy of Sciences, 1934, Section: Mathematics and Natural Sciences 6: 793-800.
[11] Emil Simiu, Robert H. Scanlen, “Wind Effect On Building: Fundamentals and Applications to Design”, 1996.
[12] Emil Simiu, Dong-Hau Yeo. “Wind Effect On Building: Modern Structural Design for Wind”, 2019.
[13] Fage, A., Johansen, F.C. On the flow of air behind an inclined flat plate of infinite span (Containing Papers of a Mathematical and Physical Character). Proc. R. Soc. Lond. , 1927, Ser. A 116, 170–197.
[14] G.T. Bitsuamlak, A.K. Dagnew, J.Erwin. “Evaluation of wind loads on solar panel modules using CFD”, in Proceedings of the Fifth International Symposium on Computational Wind Engineer-ing, Chapel Hill, North Carolina, USA, 2010.
[15] JIS, “JIS-C8955” .Loading Design Guide on Structures For Pho-tovoltaic Array, Japanese Standard Association”, 2017.
[16] MATLAB. MATLAB documentation 2018b. The Mathworks Inc., 2018.
[17] Menter, F.R. “Two-Equation Eddy-Viscosity Turbulence Models For Engineering Applications”, AIAA J., 1994, Vol. 32, pp. 1598-1605.
[18] SEAOC.Solar Photovoltaic Systems Committee, “Wind Design For Low-Profile Solar Photovoltaic Arrays On Flat Roofs”, 2012.
[19] Shademan, M., Barron, R.M., Balachandar, R., Hangan, H. “Nu-merical simulation of wind loading on ground-mounted solar panels at different flow configurations. ” Can. J. Civ. Eng. 41, 2014, 728–738.
[20] Voronoi, G. “New Applications of Continuous Parameter in the Theory of Quadratic Forms.” Journal of Pure and Applied Mathematics 133, 1908, 97-178.
[21] 內政部營建署，建築物耐風設計規範及解說，2015年。
[22] 方富民、黎益肇、賴冠廷、李唯佳，應用高精度數值地形模型進行CFD風場模擬，內政部建築研究所委託研究報告，2016年。
[23] 中華民國風工程學會，風工程理論與應用，科技圖書，2016年。
[24] 朱佳仁，風工程概論，科技圖書，2006年。
[25] 林逸崧，棚架型太陽光電系統極值風壓係數之估計，國立台灣科技大學營建工程學系碩士學位論文，2019年。
[26] 陳瑞華、傅仲麟、蔡宇勛、沈朝斌、林逸崧，太陽光電系統之耐風設計規範研擬，內政部建築研究所委託研究報告，2020年。