Superconducting radio frequency cavities used for accelerating charged particle beams are commonly used in accelerator facilities around the world. The design and optimization of modern superconducting RF cavities requires intensive numerical simulations. Vast number of operational parameters must be calculated to ensure appropriate functioning of the accelerating structures. In this study, we primarily focus on estimation and behavior of higher order modes in superconducting RF cavities connected in chains. To calculate large RF models the state-space concatenation scheme, an efficient hybrid method, is employed.