Three-Dimensional Sound Field Analysis Using Compact Explicit-Finite Difference Time Domain Method with Graphics Processing Unit Cluster System

In this study, the compact explicit-finite difference time domain (CE-FDTD) method is applied to the three-dimensional sound field analysis to reduce computer resources. There are various derivative schemes in the CE-FDTD method. They are first examined theoretically to evaluate the numerical accuracy. As a theoretical result, it is found that the interpolated wide band (IWB) scheme has the widest bandwidth in which the cut-off frequency is in agreement with the Nyquist frequency. The calculation performance is theoretically estimated, then experimentally evaluated with the graphics processing unit cluster system. As a result, it is found that the memory usage of the IWB scheme is less than one-third of that of the standard leapfrog (SLF) scheme to achieve the same cut-off frequency. It is also found that the calculation time of the IWB scheme with the shared memory is about 19% compared with that of the SLF scheme with the graphics processing unit (GPU) cluster system. The impulse response is calculated for a large room with a volume capacity of about 4500 m3 in which the sampling rate was 40 kHz. It is confirmed that the three-dimensional sound field with the natural reverberation can be calculated by the IWB scheme.

[1]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[2]  Evaluation of Acoustic Simulation Using Wave Equation Finite Difference Time Domain Method with Compact Finite Differences , 2012 .

[3]  K. Mori,et al.  Tidal Effect in Small-Scale Sound Propagation Experiment , 2012 .

[4]  Tetsuya Sakuma,et al.  Approximate theory of reverberation in rectangular rooms with specular and diffuse reflections. , 2012, The Journal of the Acoustical Society of America.

[5]  Basic Study of Properties of Planate Acoustic Lens Constructed with Phononic Crystal Structure , 2012 .

[6]  Kan Okubo,et al.  Examination of Sub-Grid Technique for Simulation of Sound Wave Propagation Using Constrained Interpolation Profile Method with Method of Characteristics , 2011 .

[7]  Hanako Ogasawara,et al.  Design and Convergence Performance Analysis of Aspherical Acoustic Lens Applied to Ambient Noise Imaging in Actual Ocean Experiment , 2011 .

[8]  Maarten van Walstijn,et al.  Wideband and Isotropic Room Acoustics Simulation Using 2-D Interpolated FDTD Schemes , 2010, IEEE Transactions on Audio, Speech, and Language Processing.

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

[10]  Two-Dimensional Numerical Analysis of Nonlinear Sound Wave Propagation Using Constrained Interpolation Profile Method Including Nonlinear Effect in Advection Equation , 2011 .

[11]  T. Tsuchiya,et al.  Numerical Analysis of Long Range Sound Wave Propagation in Ocean by Wave Equation Finite Difference Time Domain Method with Graphics Processing Unit , 2012 .

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

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

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

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

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

[17]  Staggered Grid with Collocated Grid Points of Velocities for Modeling Propagation of Elastic Waves in Anisotropic Solids by Finite-Difference Time Domain Method , 2012 .

[18]  Performance Comparison of Subgrid Techniques in Acoustic Simulation Using the Type-M and Type-C Constrained Interpolation Profile Methods , 2012 .

[19]  M. Matsukawa,et al.  Effect of Boundary Condition on the Two-Wave Propagation in Cancellous Bone , 2011 .