A Hardware-Oriented Finite-Difference Time-Domain Algorithm for Sound Field Rendering

Sound field renderings are data-intensive and computation-intensive applications. An alternative solution is to directly implement sound field rendering algorithms by using hardware. In this paper, a hardware-oriented finite-difference time-domain (FDTD) algorithm named HO-FDTD is proposed for sound field rendering, which has no complex operations involved, and consumes small hardware resources. In a sound space with 32,768 elements surrounded by rigid walls, the hardware simulation results are in good agreement with the software simulation results except for the one-cycle delay. In the software simulation, when the element scale is 32×32×32 and the time steps are 20,000, the HO-FDTD speeds up computations by 19% against the updated digital Huygen's model (DHM) and Yee-FDTD, and by 132% against the original DHM. Compared with the software simulation, the hardware systems with the parallel architecture and the time-sharing architecture enhance their calculation performance significantly in the case of different element scales, and provide a higher data throughput.

[1]  T. Tsuchiya Three-Dimensional Sound Field Rendering with Digital Boundary Condition using Graphics Processing Unit , 2010 .

[2]  Dick Botteldooren,et al.  ACOUSTICAL FINITE-DIFFERENCE TIME-DOMAIN SIMULATION IN A QUASI-CARTESIAN GRID , 1994 .

[3]  A. Krokstad,et al.  Calculating the acoustical room response by the use of a ray tracing technique , 1968 .

[4]  M. Aliabadi,et al.  Boundary‐Element Method , 2009 .

[5]  Murray Hodgson,et al.  Investigation of the validity of radiosity for sound-field prediction in cubic rooms. , 2004, The Journal of the Acoustical Society of America.

[6]  Kan Okubo,et al.  Examination on Fast Algorithms of Compact Finite Difference Calculation for Finite Difference Time Domain Acoustic Wave Simulation , 2011 .

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

[8]  Takao Tsuchiya,et al.  Finite Element Analysis of the Standing Wave Field in a Tube with Lossy Wall for Thermoacoustic System , 2011 .

[9]  Takao Tsuchiya,et al.  Numerical Simulation of Acoustic Imaging Using a Combination of Finite Difference Time Domain and Boundary Integral Equation Methods , 2009 .

[10]  T. Lewers,et al.  A combined beam tracing and radiatn exchange computer model of room acoustics , 1993 .

[11]  Kan Okubo,et al.  An Improvement of Absorbing Boundary Condition Based on Algorithm of Constrained Interpolation Profile Method for Finite-Difference Time-Domain Acoustic Simulation , 2010 .

[12]  Takao Tsuchiya,et al.  Digital equivalent circuits for acoustic field based on discrete huygens' modeling , 2005 .

[13]  Takao Tsuchiya Numerical Simulation of Sound Wave Propagation with Sound Absorption Using Digital Huygens' Model , 2007 .

[14]  Takao Tsuchiya,et al.  Perfectly Matched Layer for the Wave Equation Finite Difference Time Domain Method , 2012 .

[15]  Thomas Funkhouser,et al.  A beam tracing method for interactive architectural acoustics. , 2004, The Journal of the Acoustical Society of America.

[16]  Naoki Matsuda,et al.  A Finite-Difference Time-Domain Technique for Nonlinear Elastic Media and Its Application to Nonlinear Lamb Wave Propagation , 2012 .

[17]  Yukio Iwaya,et al.  A real-time sound field renderer based on digital Huygens' model , 2011 .

[18]  Jont B. Allen,et al.  Image method for efficiently simulating small‐room acoustics , 1976 .

[19]  Erik Lindholm,et al.  NVIDIA Tesla: A Unified Graphics and Computing Architecture , 2008, IEEE Micro.