Methods and Systems for Generating Acoustic Impulse Responses
A method for generating an impulse response for a listening point in a room includes: receiving a 3D model of the room, the position of at least one sound source in the 3D model, and acoustic properties of at least one boundary in the 3D model; using a wave based solver for determining a wave based impulse response of a wave based propagation of an impulse emitted at the at least one sound source in the 3D model and received at the listening point within a first acoustic frequency range; using a geometrical acoustics based solver for determining a geometrical impulse response of the ray based propagation of an impulse emitted at the at least one sound source in the 3D model and received at the listening point within a second acoustic frequency range; and generating the impulse response by merging the wave impulse response and the geometrical impulse response.
1 . A computer implemented method for generating an impulse response for a listening point in a room, the method comprising:
receiving a 3D model of the room, a position of at least one sound source in the 3D model of the room, and acoustic properties of at least one boundary in the 3D model of the room;
processing the 3D model of the room to generate a 3D mesh model of the room, such that an inner volume of the 3D model of the room is represented as a volumetric 3D mesh model;
partitioning the volumetric 3D mesh model into at least one sub-domain using an acoustic wave equation module;
applying a wave-based solver to each of the at least one sub-domain, such that each of the at least one sub-domain is processed independently on a central processing unit (CPU) and/or a graphical processing unit (GPU), and wherein the CPU and/or the GPU solves at least one partial differential equation associated to the at least one sub-domain; and
determining, using the wave-based solver, a wave-based impulse response of a wave-based propagation of an impulse emitted at the at least one sound source in the 3D model of the room and received at the listening point within a first acoustic frequency range.
2 . The computer implemented method according to claim 1 , wherein the CPU is at least a first CPU central processing unit.
3 . The computer implemented method according to claim 1 , further comprising performing a qualitative evaluation on the volumetric 3D mesh model evaluating computational complexity, wherein the computational complexity is correlated to critical areas of sound propagation.
4 . The computer implemented method according to claim 2 , wherein a message passing interface (MPI) is configured to handle communication between CPUs.
5 . The computer implemented method according to claim 4 , wherein the communication handled by MPI is performed following a halo exchange, wherein the halo exchange is performed across neighbouring subdomains.
6 . The computer implemented method according to claim 1 , wherein the wave-based solver comprises a time marching sub-module comprising at least a first time marching method.
7 . The computer implemented method according to claim 6 , wherein the at least first time marching method is an explicit time stepping method, a low-storage explicit Runge-Kutta algorithm (LSERK), an implicit-explicit time marching algorithm, or any combinations thereof.
8 . The computer implemented method according to claim 1 , wherein the at least one partial differential equation is further processed into discrete algebraic equations to be solved on the central processing unit.
9 . The computer implemented method according to claim 1 , further comprising:
determining, using a geometrical acoustics-based solver, a geometrical impulse response of a ray-based propagation of an impulse emitted at the at least one sound source in the 3D model of the room and received at the listening point within a second acoustic frequency range; and
generating the impulse response by merging the wave-based impulse response and the geometrical impulse response.
10 . The computer implemented method according to claim 1 , wherein the 3D model of the room comprises at least one directive sound source configured for emitting sound in a defined direction.
11 . The computer implemented method according to claim 1 , wherein the 3D mesh model is a 3D curvilinear mesh model.
12 . The computer implemented method according to claim 1 , wherein the wave-based solver applies a discontinuous Galerkin finite element method (DGFEM) a time-domain discontinuous Galerkin finite element method, a finite element method (FEM), or a spectral element method (SEM).
13 . The computer implemented method according to claim 9 , further comprising adjusting a power level of the at least one sound source such that a sound level received at a predetermined distance from the at least one sound source is the same in the wave-based solver and in the geometrical acoustic solver.
14 . The computer implemented method according to claim 9 , wherein using the wave-based solver or the geometrical acoustic solver comprises extracting at least one wave impulse response or at least one geometrical impulse response based on a simulation of a propagation of sound, wherein the propagation of sound is the wave-based propagation or the ray-based propagation, and wherein the at least one wave-based impulse response or the at least one geometrical impulse response are at least one spatial impulse response.
15 . The computer implemented method according to claim 14 , wherein the at least one spatial impulse response comprises a plurality of single channel impulse responses, wherein each one of the plurality of single channel impulse responses corresponds to an auxiliary wave-based impulse response weighted with a specific angular distribution from the same listening point.
16 . The computer implemented method according to claim 1 , wherein a spherical receiver array is arranged around the listening point, and wherein the spherical receiver array comprises a plurality of receivers.
17 . The computer implemented method according to claim 16 , wherein the spherical receiver array is an open spherical array of cardioid receivers.
18 . The computer implemented method according to claim 16 , wherein a number of the plurality of receivers is determined based on a maximum ambisonics truncation order N, such that the number of the plurality of receivers is higher or equal to (N+1) 2 .
19 . The computer implemented method according to claim 1 , further comprising rendering a base audio signal by convolving a base audio signal with the impulse response or the wave-based impulse response, thereby creating a rendered audio signal.
20 . A system for generating an impulse response for a listening point in a room, the system comprising:
a computer system having processing circuitry coupled to a memory, the processing circuitry configured to:
receive a 3D model of the room, a position of at least one sound source in the 3D model of the room, and acoustic properties of at least one boundary in the 3D model of the room;
process the 3D model of the room to generate a 3D mesh model of the room, such that an inner volume of the 3D model of the room is represented as a volumetric 3D mesh model;
partition the volumetric 3D mesh model into at least one sub-domain using an acoustic wave equation module;
apply a wave-based solver to each of the at least one sub-domain, such that each of the at least one sub-domain is processed independently on a central processing unit (CPU) and/or a graphical processing unit (GPU), and wherein the CPU and/or the GPU solves at least one partial differential equation associated to the at least one sub-domain; and
determine, using a wave-based solver, a wave-based impulse response of a wave-based propagation of an impulse emitted at the at least one sound source in the 3D model of the room and received at the listening point within a first acoustic frequency range.