×

openPSTD: the open source pseudospectral time-domain method for acoustic propagation. (English) Zbl 1380.65472

Summary: An open source implementation of the Fourier pseudospectral time-domain (PSTD) method for computing the propagation of sound is presented, which is geared towards applications in the built environment. Being a wave-based method, PSTD captures phenomena like diffraction, but maintains efficiency in processing time and memory usage as it allows to spatially sample close to the Nyquist criterion, thus keeping both the required spatial and temporal resolution coarse. In the implementation it has been opted to model the physical geometry as a composition of rectangular two-dimensional subdomains, hence initially restricting the implementation to orthogonal and two-dimensional situations. The strategy of using subdomains divides the problem domain into local subsets, which enables the simulation software to be built according to Object-Oriented Programming best practices and allows room for further computational parallelization. The software is built using the open source components, Blender, Numpy and Python, and has been published under an open source license itself as well. For accelerating the software, an option has been included to accelerate the calculations by a partial implementation of the code on the Graphical Processing Unit (GPU), which increases the throughput by up to fifteen times. The details of the implementation are reported, as well as the accuracy of the code.

MSC:

65Y15 Packaged methods for numerical algorithms
65M70 Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs
35Q51 Soliton equations
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Välimäki, V.; Parker, J. D.; Savioja, L.; Smith, J. O.; Abel, J. S., IEEE Trans. Audio Speech Lang. Process., 20, 1421-1448 (2012)
[2] Vorländer, M., J. Acoust. Soc. Am., 133, 1203-1213 (2013)
[3] Hornikx, M.; Hak, C.; Wenmaekers, R., J. Build. Perform. Simul., 8, 1, 26-38 (2015)
[4] Mehra, R.; Raghuvanshi, N.; Savioja, L.; Lin, M.; Manocha, D., Appl. Acoust., 73, 2, 83-94 (2012)
[6] Kowalczyk, K.; van Walstijn, M., IEEE Trans. Audio Speech Lang. Process., 19, 1, 34-46 (2011)
[8] Hornikx, M.; Waxler, R.; Forssén, J., J. Acoust. Soc. Am., 128, 1632-1646 (2010)
[9] Hornikx, M.; De Roeck, W.; Desmet, W., J. Comput. Phys., 231, 14, 4759-4774 (2012)
[10] Hornikx, M.; Dragna, D., J. Acoust. Soc. Am., 138, 1, 425-435 (2015)
[11] Hornikx, M.; De Roeck, W.; Toulorge, T.; Desmet, W., Comput. Fluids, 116, 176-191 (2015)
[12] Van Renterghem, T.; Hornikx, M.; Forssén, J.; Botteldooren, D., Build. Environ., 61, 34-44 (2013)
[13] Bashir, I.; Hill, T.; Taherzadeh, S.; Attenborough, K.; Hornikx, M., Appl. Acoust., 83, 1-15 (2014)
[17] Ostashev, V. E.; Wilson, D. K.; Liu, L.; Aldridge, D. F., J. Acoust. Soc. Am., 117, 503-517 (2005)
[20] Oliphant, T. E., Comput. Sci. Eng., 9, 10-20 (2007)
[24] Bogey, C.; Bailly, C., J. Comput. Phys., 194, 194-214 (2004)
[25] Escolano, J.; Spa, C.; Garriga, A.; Mateos, T., Appl. Acoust., 74, 818-822 (2013)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.