An efficient GPU-based time domain solver for the acoustic wave equation

An efficient algorithm for time-domain solution of the acoustic wave equation for the purpose of room acoustics is presented. It is based on adaptive rectangular decomposition of the scene and uses analytical solutions within the partitions that rely on spatially invariant speed of sound. This technique is suitable for auralizations and sound field visualizations, even on coarse meshes approaching the Nyquist limit. It is demonstrated that by carefully mapping all components of the algorithm to match the parallel processing capabilities of graphics processors (GPUs), significant improvement in performance is gained compared to the corresponding CPU-based solver, while maintaining the numerical accuracy. Substantial performance gain over a high-order finite-difference time-domain method is observed. Using this technique, a 1 s long simulation can be performed on scenes of air volume 7500 m³ till 1650 Hz within 18 min compared to the corresponding CPU-based solver that takes around 5 h and a high-order finite-difference time-domain solver that could take up to three weeks on a desktop computer. To the best of the authors’ knowledge, this is the fastest time-domain solver for modeling the room acoustics of large, complex-shaped 3D scenes that generates accurate results for both auralization and visualization.     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

Sound propagation above an inhomogeneous plane: Boundary integral equation methods

A boundary integral equation method is used to compute the sound pressure emitted by a harmonic source above an inhomogeneous plane. First, the theoretical aspects of the problem (behaviour of the pressure around the discontinuities,…) are studied. Then, a comparison between theoretical levels and experimental levels obtained in an anechoic room is presented. It shows that the boundary integral equation (BIE) method is quite convenient for solving this kind of problem. Two interesting results are pointed out: (i) if only a prediction of maximum sound levels is needed, the attenuation is the same for a cylindrical source, a spherical source and N spherical sources, and so it is possible to transform some three-dimensional problems into two-dimensional ones; (ii) a numerical method of computation of the sound field above an inhomogeneous plane does not provide a correct prediction if each part of the plane is not accurately described by the boundary condition chosen.     

Waveguide mode solver based on Neumann-to-Dirichlet operators and boundary integral equations

For optical waveguides with high index-contrast and sharp corners, existing full-vectorial mode solvers including those based on boundary integral equations typically have only second or third order of accuracy. In this paper, a new full-vectorial waveguide mode solver is developed based on a new formulation of boundary integral equations and the so-called Neumann-to-Dirichlet operators for sub-domains of constant refractive index. The method uses the normal derivatives of the two transverse magnetic field components as the basic unknown functions, and it offers higher order of accuracy where the order depends on a parameter used in a graded mesh for handling the corners. The method relies on a standard Nyström method for discretizing integral operators and it does not require analytic properties of the electromagnetic field (which are singular) at the corners.     

Thermodynamic functions of dense equilibrium plasmas using the allnatt integral equation

Solving the nonlinear Allnatt integral equation in the HNC-approximation the radial distribution functions (r df) and thermodynamic functions for the subsystem of free electrons and ions in a dense plasma were calculated.     

Free vibration of annular sector plates by an integral equation technique

A numerical method is presented for the free vibration analysis of polar orthotropic clamped annular sector plates. The results are compared with analytical and experimental values of other investigators. A parametric study has been done by varying the sector angle and radii ratio. The frequencies for isotropic and orthotropic cases are presented in the form of graphs.     

Vibration analysis of composite annular plates by an integral equation technique

In layered plates, coupling between bending and stretching makes the analysis complicated. A numerical method, namely an integral equation technique which has been used for solving static and dynamic problems, is extended to vibration analysis of layered annular plates which are fixed at the inner and outer radii. The results obtained are compared with those of other investigators and a parametric study is made of the effect of layer thickness, the layers being made up of two different materials.     

用时域积分方程法计算不规则形状目标的瞬态声散射

瞬态声散射信号包含有较丰富的目标信息。本文研究了克服时域积分方程法计算瞬态声散射中存在的长时间解发散问题的方法,给出了一种实用的算法,使得计算结果在时间较长时都稳定。将时域积分方程法推广到阻抗表面的情况,计算了带半椭球帽的阻抗圆柱体的瞬态反向声散射。计算结果表明,目标的瞬态回波具有亮点特征.     

Comparison of an integral equation on energy and the ray-tracing technique in room acoustics

This paper deals with a comparison of two room acoustic models. The first one is an integral formulation stemming from power balance and the second is the ray-tracing technique with a perfectly diffuse reflection law. The common assumptions to both models are the uncorrelated wave hypothesis and the perfectly diffuse reflection law. The latter allows the use of these methods for nondiffuse fields beyond the validity domain of Sabine's formula. Comparisons of numerical simulations performed with the softwares RAYON and CeReS point out that these results are close to each other and finally, a formal proof is proposed showing that both methods are actually equivalent.     

Fast analysis of transient acoustic wave scattering from rigid bodies using the multilevel plane wave time domain algorithm

The analysis of transient wave scattering from rigid bodies using integral equation-based techniques is computationally intensive: if carried out using classical schemes, the evaluation of the velocity potential on the surface of a three-dimensional scatterer, represented in terms of Ns spatial basis functions for Nt time steps, requires O(NtNs2) operations. The recently developed plane wave time domain (PWTD) algorithm permits the rapid evaluation of transient fields that are generated by bandlimited source distributions. It has been shown that incorporation of the PWTD algorithm into integral equation-based solvers in a two-level setting reduces the computational complexity of a transient analysis to O(NtNs1.5 logNs). In this paper, it is shown that casting the PWTD scheme into a multilevel framework permits the analysis of transient acoustic surface scattering phenomena in O(NtNslog2Ns) operations using O(NtNs) memory. Numerical examples that demonstrate the efficacy of the multilevel implementation are also presented.     

Phase integral theory,coupled wave equations,and mode conversion

Phase integral or WKB theory is applied to multicomponent wave equations, i.e., wave equations in which the wave field is a vector, spinor, or tensor of some kind. Specific examples of physical interest often have special features that simplify their analysis, when compared with the general theory. The case of coupled channel equations in atomic or molecular scattering theory in the Born-Oppenheimer approximation is examined in this context. The problem of mode conversion, also called surface jumping or Landau-Zener-Stuckelberg transitions, is examined in the multidimensional case, and cast into normal form. The group theoretical principles of the normal form transformation are laid out, and shown to involve both the Lorentz group and the symplectic group.     

High-frequency plane wave diffraction by an ideal strip at oblique incidence: Parabolic equation approach

The problem of diffraction of a high-frequency plane wave by a strip with ideal boundary conditions is considered for the case of oblique incidence. The study is based on the parabolic approximation, which is used to construct an expression for the directional pattern in terms of single quadratures. A similar result is obtained using the embedding formula. It is shown that the derived expression approximates the classical Michaeli result. A proof of the optical theorem for the parabolic problem is presented.     

Path integral solution of the time-domain progressive wave equation and wave propagation in random media

A time-domain progressive wave equation is derived from the usual linear acoustic wave equation, and it is shown that the solution to this new equation can be expressed as a Feynman path integral. This path integral representation is used to derive the time-dependent statistics of acoustic fields propagating through random media.     

Path integral treatment for relativistic particles in external fields and quantized plane wave

The dynamics of relativistic particles of spin 0 and 1/2, interacting with an external electromagnetic field and a quantized plane wave, is studied using the path integral framework. We take advantage of the existing properties of the interaction to introduce a delta functional in order to calculate Green's functions. This simply reduces the problem to two coupled oscillators. The energy spectrum and wave functions are calculated for the spin 0 case and the analogy with Jaynes‐Cummings model is made to derive the energy spectrum for the spin 1/2 case.     

Transmission loss prediction of reactive silencers using 3-D time-domain CFD approach and plane wave decomposition technique

A predictive method is proposed to determine the transmission loss of reactive silencers using the three-dimensional (3-D) time-domain computational fluid dynamics (CFD) approach and the plane wave decomposition technique. Firstly, a steady flow computation is performed with a mass-flow-inlet boundary condition, which provides an initial condition for the following two unsteady flow computations. The first unsteady flow computation is conducted by imposing an impulse (acoustic excitation) superimposed on the constant mass flow at the inlet of the model and then adding the non-reflecting boundary condition (NRBC) when the impulse completely propagates into the silencer. The second unsteady flow computation is conducted for the case without acoustic excitation at the inlet. The time histories of pressure and velocity at the upstream monitoring point as well as history of pressure at the downstream monitoring point are recorded during the two transient computations. The differences between the two unsteady flow computational results are the corresponding acoustic quantities. Therefore, the incident sound pressure signal is obtained by using plane wave decomposition at upstream, while the transmitted sound pressure signal is just the sound pressure at downstream. Finally, those two sound pressure signals in the time-domain are transformed into the frequency-domain by Fast Fourier Transform (FFT) and then the transmission loss (TL) of silencer is determined. For the straight-through perforated tube silencers with and without flow, the numerical results agree well with the published measurements.     

On-line detection of fatigue cracks using an automatic mode tracking technique

Experimental fatigue tests usually require large testing times. In addition to the resulting increased time-to-market, the large fatigue test time also implies that any structural health monitoring technique that is used should be automatic. When using the modal parameters as damage indicators, an important amount of user interaction is still needed to separate physical poles from computational ones. In this paper, an experimental framework will be developed to automatically track the health of the structure on-line with the performance of fatigue tests. The modal parameters are tracked using a combination of the maximum likelihood estimator and an auto-regressive model. Since confidence levels on the modal parameter are available it is possible to detect if damage is present. In addition, the quasi-static stiffness with computed confidence levels is also used as a damage indicator. The proposed techniques are demonstrated on a steel beam with a propagating fatigue crack.     

Momentum space distribution of electrons in an atom using hydrogenic wave functions

A formulation is developed to derive exact analytic expressions for electron-electron correlation and density of electrons in momentum space using hydrogenic wave functions. It is shown that for large atoms the expression for density of electrons has a simple form.     

Path integral for spinning particle in the plane wave field: Global and local projections

The Green function related to the problem of a Dirac particle interacting with a plane wave is calculated via the path integral formalism proposed recently by Alexandrou et al. according to the two so-called global and local projections. With the help of the incorporation of two simple identities, it is shown that the contribution to the calculation of the integrals comes essentially from classical solutions projected along the direction of wave propagation.     

超声平面波同相正交信号频域波束形成算法

针对频域平面波波束形成算法尚不能直接应用于采用正交解调的超声成像系统中的实际问题,提出一种平面波同相正交信号频域波束形成方法。以平面波射频信号频域波束形成算法为基础,根据射频信号与同相正交信号波数空间变量之间的关系,推导了平面波同相正交信号与对应图像的频谱映射表达式,从而实现对正交解调后信号的频域超声图像重建。采用2016年Plane Wave Imaging Challenge in Medical Ultrasound (PICMUS)发布的仿真、仿体及在体原始回波数据对这一方法进行验证,并评估图像的分辨率和对比噪声比。仿真和仿体实验结果表明,相对于射频信号的频域波束形成算法,该算法可获得与之相同的图像质量,并将成像速度提高了约4倍。在体实验数据的成像结果也验证了该算法在实际应用中的有效性。