Improved Finite Difference Schemes for a 3-D Viscothermal Wave Equation on a GPU

Viscothermal effects in air lead to a damping of high frequencies over time. Such effects cannot be neglected in large-scale room acoustics simulations for the full audible bandwidth. In this study, full-bandwidth room acoustics is modelled using a variant of the three-dimensional wave equation including viscothermal losses in air following from a simplification of the Navier-Stokes equations suitable for room acoustics applications. The model equation is numerically solved using time domain finite difference methods. A three-step parameterised finite difference scheme is proposed to model T60 decay times as a function of frequency more accurately than two-step schemes. Timing results from parallelised implementations on a graphics processing unit (GPU) device are presented.

[1]  B. B. Bauer,et al.  Fundamentals of acoustics , 1963 .

[2]  Allan D. Pierce,et al.  Acoustics , 1989 .

[3]  D. Botteldooren.  Finite‐difference time‐domain simulation of low‐frequency room acoustic problems , 1995 .

[4]  C. Groot‐Hedlin.  Finite-difference time-domain synthesis of infrasound propagation through an absorbing atmosphere. , 2008 .

[5]  J. Strikwerda.  Finite Difference Schemes and Partial Differential Equations , 1989 .

[6]  Stefan Bilbao,et al.  Computing room acoustics with CUDA - 3D FDTD schemes with boundary losses and viscosity , 2011, 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[7]  Brian Hamilton,et al.  ROOM ACOUSTICS MODELLING USING GPU-ACCELERATED FINITE DIFFERENCE AND FINITE VOLUME METHODS ON A FACE-CENTERED CUBIC GRID , 2013 .

[8]  Dingwen Deng,et al.  Accuracy improvement of a multistep splitting method for nonlinear viscous wave equations , 2014, Int. J. Comput. Math..

[9]  Stefan Bilbao,et al.  Optimized FDTD Schemes for 3-D Acoustic Wave Propagation , 2012, IEEE Transactions on Audio, Speech, and Language Processing.

[10]  Maarten van Walstijn,et al.  Room Acoustics Simulation Using 3-D Compact Explicit FDTD Schemes , 2011, IEEE Transactions on Audio, Speech, and Language Processing.

[11]  Craig J. Webb.  Parallel computation techniques for virtual acoustics and physical modelling synthesis , 2014 .

[12]  Zhiyue Zhang.  The multistep finite difference fractional steps method for a class of viscous wave equations , 2010 .

[13]  D. Borup,et al.  Simulation of acoustic wave propagation in dispersive media with relaxation losses by using FDTD method with PML absorbing boundary condition , 1999, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[14]  C de Groot-Hedlin,et al.  Finite-difference time-domain synthesis of infrasound propagation through an absorbing atmosphere. , 2008, The Journal of the Acoustical Society of America.

[15]  R. M. Alford,et al.  ACCURACY OF FINITE‐DIFFERENCE MODELING OF THE ACOUSTIC WAVE EQUATION , 1974 .

[16]  K. Yoshihisa.  Calculation of the Absorption of Sound by the Atmosphere , 1997 .

[17]  P. Brandimarte.  Finite Difference Methods for Partial Differential Equations , 2006 .

[18]  Stefan Bilbao,et al.  Wave and scattering methods for the numerical integration of partial differential equations , 2001 .

[19]  Stefan Bilbao,et al.  Modeling of Complex Geometries and Boundary Conditions in Finite Difference/Finite Volume Time Domain Room Acoustics Simulation , 2013, IEEE Transactions on Audio, Speech, and Language Processing.

[20]  Damian T. Murphy,et al.  Room Impulse Response Synthesis and Validation Using a Hybrid Acoustic Model , 2013, IEEE Transactions on Audio, Speech, and Language Processing.

[21]  Stefan Bilbao,et al.  ON FINITE DIFFERENCE SCHEMES FOR THE 3-D WAVE EQUATION USING NON-CARTESIAN GRIDS , 2013 .

[22]  Tapio Lokki,et al.  The perceptual effects of dispersion error on room acoustic model auralization , 2011 .

[23]  Michel Bruneau,et al.  General formulation of the dispersion equation in bounded visco-thermal fluid, and application to some simple geometries , 1989 .

[24]  Mark Staniland Wochner.  Numerical simulation of multi-dimensional acoustic propagation in air including the effects of molecular relaxation , 2006 .

[25]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[26]  Craig J. Webb,et al.  Large-scale Virtual Acoustics Simulation at Audio Rates Using Three Dimensional Finite Difference Time Domain and Multiple GPUs , 2013 .

[27]  Lauri Savioja,et al.  Integrating finite difference schemes for scalar and vector wave equations , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

[28]  M. van Walstijn,et al.  On the numerical solution of the 2d wave equation with compact fdtd schemes , 2008 .

[29]  Michael J Buckingham.  Causality, Stokes' wave equation, and acoustic pulse propagation in a viscous fluid. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .

[31]  M. Cramer.  Numerical estimates for the bulk viscosity of ideal gases , 2012 .

[32]  Kan Okubo,et al.  Three-Dimensional Sound Field Analysis Using Compact Explicit-Finite Difference Time Domain Method with Graphics Processing Unit Cluster System , 2013 .

[33]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[34]  Allen Taflove,et al.  Application of the Finite-Difference Time-Domain Method to Sinusoidal Steady-State Electromagnetic-Penetration Problems , 1980, IEEE Transactions on Electromagnetic Compatibility.

[35]  J. Crank.  Difference Methods for Initial-Value Problems 2nd edn , 1968 .