This paper presents our initial work in performing large-eddy simulations of tidal turbine array flows. First, a horizontally-periodic precursor simulation is performed to create turbulent flow data. Then that data is used as inflow into a tidal turbine array two rows deep and infinitely wide. The turbines are modelled using rotating actuator lines, and the finite-volume method is used to solve the governing equations. In studying the wakes created by the turbines, we observed that the vertical shear of the inflow combined with wake rotation causes lateral wake asymmetry. Also, various turbine configurations are simulated, and the total power production relative to isolated turbines is examined. Staggering consecutive rows of turbines in the simulated configurations allows the greatest efficiency using the least downstream row spacing. Counter-rotating consecutive downstream turbines in a non-staggered array shows a small benefit. This work has identified areas for improvement, such as the use of a larger precursor domain to better capture elongated turbulent structures, the inclusion of salinity and temperature equations to account for density stratification and its effect on turbulence, improved wall shear stress modelling, and the examination of more array configurations.