The double vertical-axis wind turbine (VAWT) system serves as a high-performance design solution. The Nautilus wind turbine investigated in this research imitated the structure of a Nautilus shell and is a vertical-axis drag-type wind turbine that exhibits relatively low efficiency. Therefore, the improvement of its wind energy efficiency is of paramount importance. This paper utilizes Computational Fluid Dynamics (CFD) software, based on the Reynolds-averaged Navier–Stokes equations and dynamic meshing techniques, to conduct numerical investigations on the aerodynamic performance of the Nautilus wind turbine array layout. The effects of wind direction, spacing ratio, and rotation direction are individually studied, and the interpretations and explanations are provided based on flow field characteristics. The results show that when the wind direction is 90°, i.e., a transverse layout, the closer the spacing between the transverse turbines, the higher the average power coefficient of the entire wind turbine system, with little effect from the three rotation directions. The maximum average power coefficient reached 28.9% and the power gain factor (TPGF) reached 11.1%. The enhancement effect primarily originates from the wake interaction among neighboring turbines. The experimental results showed a deviation of 8.1% compared to the CFD simulation results, thus validating the accuracy of the numerical CFD modeling. Ultimately, several array layouts are proposed, based on the prevalent wind direction and spacing ratio research. The enhancement of the wind turbine array’s situation could significantly increase the average efficiency of the entire wind turbine cluster. Consequently, this study provides a reference for the practical application of biomimetic vertical-axis drag-type wind turbine systems in actual engineering.