Natural convection in cavity with built-in plates has been widely used in fluid mechanics and heat transfer engineering, such as heat exchangers, electronics, solar collector, growing crystals, etc. However, the natural convection flow physics and heat transfer characteristics in cavity with plate arrays remain to be investigated further. In this study, numerical analysis and computational fluid dynamics simulation are conducted for natural convection in cavity with plate arrays by finite volume method. Key influence factors impacting heat transfer are analyzed, including the line dip-angle of plate arrays θ, rotary angle of plates to the array φ, with relationships comparison between the Nusselt number and Rayleigh number in different built-in plate arrays. Moreover, for illustrative application example, the natural convection in solar greenhouse with PV arrays roof is investigated, with inside temperature distribution comparison under different design parameters. The preliminary results indicate that plates spacing could remarkably impacts the number of streamlines between ”late’arrays. The dip-angle of arrays θ leads to different shape flow in the cavity, and the rotary angle of plates φ notably affects the fluid flow direction in the vicinity areas of plates, especially when φ = 90°. Solar greenhouse case study indicates that increasing array group numbers can improve air temperature distribution uniformity in the cavity, which is favorable for built environment accurate control and operation. This work can provide modeling method support and reference for natural heat convection applications.