## 1 Introduction
Furthermore, these tools have been employed for post-operative assessment after glottis-widening surgery to improve breathing in patients with bilateral vocal fold paralysis (BVFP) [43] and also in post-endoscopic sinus surgery of nasal cavity airflow [40]. In the field of voice disorders, only a few studies have been focused on the use of computational tools to directly support clinical work. A major problem is the incorporation of different scales in the flow and tissue dynamics and acoustics and their interactions in the phonation process. [PERSON] et al. [44] employed a high-resolution model based on the immersed boundary method for a fully coupled phonation process as a pre-surgery tool to identify the optimal placement position of an implant for the medialization of the paralyzed vocal fold. Additionally, they evaluated the sound sources with a hybrid approach. Unfortunately, they did not mention the computational costs of their simulations to evaluate the potential application in clinics. [PERSON] et al. [28; 45] later employed this model to investigate phonatory disorders based on vocal fold left-right asymmetry and tension imbalance and reported that a huge computational time and power were required for each of their simulations. [PERSON] and [PERSON] [27] also investigated numerically another pathology, subglottic stenosis (SGS), without providing the computational costs. In a recent study, [PERSON] [46] used a reduced-order model of the vocal fold dynamics based on the eigenmodes of oscillations to reduce the computational effort in order to find a real-time simulation model of voice production. However, for the flow excitation, they applied a quasi-one-dimensional model. Although they claimed that the effects of flow three-dimensionality have less importance in the phonation process, many previous studies proved the significance of these effects such as glottal jet instabilities, its axis switching and the transition to turbulence on laryngeal aerodynamics and acoustics [29; 47; 48; 49]. In a previous study [50], we showed that pure flow simulation in the numerical larynx model is able to capture all the characteristic features of laryngeal aerodynamics if realistic oscillation patterns of vocal folds are applied as the moving boundaries on the model. Furthermore, such models with forced motion of vocal fold oscillation have been reported in the literature as powerful tools for the analysis of sound generation within the larynx [29; 30; 31; 51] with the advantage of saving the computational time consumed by a fully coupled flow-structure interaction. In our study, it was demonstrated that the total wall time (elapsed real-time) to simulate laryngeal aerodynamics with high accuracy is in the range of weeks even when employing high-performance computing (HPC) resources. This was found to be in agreement with other studies. Therefore, CFD analysis may not yet be applicable in clinical diagnostics and treatment. The question arises of whether the accuracy of the computational models can be reduced without losing the ability to provide the relevant information for supporting physicians' planning and decision-making processes during the treatment of patients with voice disorders. To answer this question, we reduced systematically the spatial and temporal resolution of the model and evaluated the results considering the computational costs. By making optimal use of a high-performance parallel computer cluster, we performed Large Eddy Simulations (LES) of laryngeal aerodynamics with externally imposed vocal fold dynamics for different mesh resolutions. The validity of the simulations was further evaluated by comparing the results with experimental data obtained from a synthetic larynx model with silicone vocal folds.