@article {Kantartzis:2007:1610-1928:870, title = "Compact Higher-Order Time-Domain Schemes for the Unconditionally Stable Modeling of Generalized 3-D Acoustic Propagation Problems", journal = "Acta Acustica united with Acustica", parent_itemid = "infobike://dav/aaua", publishercode ="dav", year = "2007", volume = "93", number = "6", publication date ="2007-11-01T00:00:00", pages = "870-887", itemtype = "ARTICLE", issn = "1610-1928", url = "https://www.ingentaconnect.com/content/dav/aaua/2007/00000093/00000006/art00003", author = "Kantartzis, Nikolaos V. and Prokopidis, Konstantinos P. and Antonopoulos, Christos S. and Tsiboukis, Theodoros D.", abstract = "A higher-order curvilinear finite-difference time-domain (FDTD) technique, combined with the alternating-direction implicit (ADI) concept, is presented in this paper for the unconditionally stable characterization and optimal design of modern 3-D acoustic propagation arrangements. Departing from a compact non-standard discretization discipline, applied directly to the general form of the linearized Euler's equations, the novel method develops a family of accurate spatial/temporal operators and introduces a dual-mesh framework for the correct association of velocity components and pressure values. Moreover, the presence of an arbitrary mean flow vector in lossy regions is efficiently treated via a stencil control process that selects the appropriate candidate according to the problem, under study. On the other hand, different media interfaces, not necessarily aligned with grid's cells, are elaborately analyzed in terms of a convergent algorithm which satisfies the physical continuity conditions and exhibits enhanced directional sensitivity. To truncate open-boundary regions and ensure satisfactory annihilation rates for the intricate vorticity and entropy waves, the higher-order technique introduces a parameterized perfectly matched layer (PML) absorber. Finally, the frequent case of bent ducts with flexible or rigid walls and changeable cross-section is rigorously manipulated by a multimodal decomposition process that decouples non-separable modes. Hence, the detrimental dispersion errors of customary approaches are drastically subdued and the prolonged duration of time-domain simulations is safely circumvented through time increments far above the Courant criterion. Numerical results, concerning various modern acoustic structures as well as extensive comparisons with reference data, confirm the preceding merits and yield optimized configurations.", }