In this thesis, a model has been developed to analyze the behavior of wave energy converter arrays in the frequency domain. This development serves the purpose of establishing an in-house model tool for array modelling for CorPower Ocean. A literature review has been conducted to compare various modelling approaches and prior studies on different configurations of wave energy converter arrays have been examined. The literature review is followed by theory that includes linear wave theory and the equation of motion. Following this, the model is introduced, accompanied by behavior validation in comparison to similar models. The selected case studies for this thesis are divided into two parts. Firstly, the investigation of different layout configurations with 16 Wave Energy Converters: 1) Circular configuration, 2) Double Circular configuration, 3) Double Row configuration (2x8), 4) Square Grid configuration (4x4), and 5) Randomized layout. The second chosen study case involves investigating layout 3) Double Row configuration (2x8) with four different distances between each Wave Energy Converter, 50, 100, 150, and 200 meters. The results demonstrate the presence of interactions within the wave energy converter arrays. The Double row layout exhibits the largest power outputs but power drops at certain waveperiods were also observed. Additionally, the results of different distances between the Wave Energy Converters demonstrate that larger spacings result in reduced interactions.