A transient boundary element method model of Schroeder diffuser scattering using well mouth impedance

Room acoustic diffusers can be used to treat critical listening environments to improve sound quality. One popular class is Schroeder diffusers, which comprise wells of varying depth separated by thin fins. This paper concerns a new approach to enable the modelling of these complex surfaces in the time domain. Mostly, diffuser scattering is predicted using steady-state, single frequency methods. A popular approach is to use a frequency domain Boundary Element Method (BEM) model of a box containing the diffuser, where the mouth of each well is replaced by a compliant surface with appropriate surface impedance. The best way of representing compliant surfaces in time domain prediction models, such as the transient BEM is, however, currently unresolved. A representation based on surface impedance yields convolution kernels which involve future sound, so is not compatible with the current generation of time-marching transient BEM solvers. Consequently, this paper proposes the use of a surface reflection kernel for modelling well behaviour and this is tested in a time domain BEM implementation. The new algorithm is verified on two surfaces including a Schroeder diffuser model and accurate results are obtained. It is hoped that this representation may be extended to arbitrary compliant locally reacting materials

A space-time mixed Galerkin marching-on-in-time scheme for the time-domain combined field integral equation

The time domain combined field integral equation (TD-CFIE), which is constructed from a weighted sum of the time domain electric and magnetic field integral equations (TD-EFIE and TD-MFIE) for analyzing transient scattering from closed perfect electrically conducting bodies, is free from spurious resonances. The standard marching-on-in-time technique for discretizing the TD-CFIE uses Galerkin and collocation schemes in space and time, respectively. Unfortunately, the standard scheme is theoretically not well understood: stability and convergence have been proven for only one class of space-time Galerkin discretizations. Moreover, existing discretization schemes are nonconforming, i.e., the TD-MFIE contribution is tested with divergence conforming functions instead of curl conforming functions. We therefore introduce a novel space-time mixed Galerkin discretization for the TD-CFIE. A family of temporal basis and testing functions with arbitrary order is introduced. It is explained how the corresponding interactions can be computed efficiently by existing collocation-in-time codes. The spatial mixed discretization is made fully conforming and consistent by leveraging both Rao-Wilton-Glisson and Buffa-Christiansen basis functions and by applying the appropriate bi-orthogonalization procedures. The combination of both techniques is essential when high accuracy over a broad frequency band is required

Stability Properties of the Time Domain Electric Field Integral Equation Using a Separable Approximation for the Convolution with the Retarded Potential

The state of art of time domain integral equation (TDIE) solvers has grown by leaps and bounds over the past decade. During this time, advances have been made in (i) the development of accelerators that can be retrofitted with these solvers and (ii) understanding the stability properties of the electric field integral equation. As is well known, time domain electric field integral equation solvers have been notoriously difficult to stabilize. Research into methods for understanding and prescribing remedies have been on the uptick. The most recent of these efforts are (i) Lubich quadrature and (ii) exact integration. In this paper, we re-examine the solution to this equation using (i) the undifferentiated form of the TD-EFIE and (ii) a separable approximation to the spatio-temporal convolution. The proposed scheme can be constructed such that the spatial integrand over the source and observer domains is smooth and integrable. As several numerical results will demonstrate, the proposed scheme yields stable results for long simulation times and a variety of targets, both of which have proven extremely challenging in the past.Comment: 9 pages, 13 figures. To be published in IEEE Transactions on Antennas and Propagatio

Parallel Implementation of the Discrete Green's Function Formulation of the FDTD Method on a Multicore Central Processing Unit

Parallel implementation of the discrete Green's function formulation of the finite-difference time-domain (DGF-FDTD) method was developed on a multicore central processing unit. DGF-FDTD avoids computations of the electromagnetic field in free-space cells and does not require domain termination by absorbing boundary conditions. Computed DGF-FDTD solutions are compatible with the FDTD grid enabling the perfect hybridization of FDTD with the use of time-domain integral equation methods. The developed implementation can be applied to simulations of antenna characteristics. For the sake of example, arrays of Yagi-Uda antennas were simulated with the use of parallel DGF-FDTD. The efficiency of parallel computations was investigated as a function of the number of current elements in the FDTD grid. Although the developed method does not apply the fast Fourier transform for convolution computations, advantages stemming from the application of DGF-FDTD instead of FDTD can be demonstrated for one-dimensional wire antennas when simulation results are post-processed by the near-to-far-field transformation

Distributed-memory parallelization of an explicit time-domain volume integral equation solver on Blue Gene/P

Two distributed-memory schemes for efficiently parallelizing the explicit marching-on in-time based solution of the time domain volume integral equation on the IBM Blue Gene/P platform are presented. In the first scheme, each processor stores the time history of all source fields and only the computationally dominant step of the tested field computations is distributed among processors. This scheme requires all-to-all global communications to update the time history of the source fields from the tested fields. In the second scheme, the source fields as well as all steps of the tested field computations are distributed among processors. This scheme requires sequential global communications to update the time history of the distributed source fields from the tested fields. Numerical results demonstrate that both schemes scale well on the IBM Blue Gene/P platform and the memory efficient second scheme allows for the characterization of transient wave interactions on composite structures discretized using three million spatial elements without an acceleration algorithm

Cosimulation of Electromagnetics-Circuit Systems Exploiting DGTD and MNA

published_or_final_versio

Fast time- and frequency-domain finite-element methods for electromagnetic analysis

Fast electromagnetic analysis in time and frequency domain is of critical importance to the design of integrated circuits (IC) and other advanced engineering products and systems. Many IC structures constitute a very large scale problem in modeling and simulation, the size of which also continuously grows with the advancement of the processing technology. This results in numerical problems beyond the reach of existing most powerful computational resources. Different from many other engineering problems, the structure of most ICs is special in the sense that its geometry is of Manhattan type and its dielectrics are layered. Hence, it is important to develop structure-aware algorithms that take advantage of the structure specialties to speed up the computation. In addition, among existing time-domain methods, explicit methods can avoid solving a matrix equation. However, their time step is traditionally restricted by the space step for ensuring the stability of a time-domain simulation. Therefore, making explicit time-domain methods unconditionally stable is important to accelerate the computation. In addition to time-domain methods, frequency-domain methods have suffered from an indefinite system that makes an iterative solution difficult to converge fast. The first contribution of this work is a fast time-domain finite-element algorithm for the analysis and design of very large-scale on-chip circuits. The structure specialty of on-chip circuits such as Manhattan geometry and layered permittivity is preserved in the proposed algorithm. As a result, the large-scale matrix solution encountered in the 3-D circuit analysis is turned into a simple scaling of the solution of a small 1-D matrix, which can be obtained in linear (optimal) complexity with negligible cost. Furthermore, the time step size is not sacrificed, and the total number of time steps to be simulated is also significantly reduced, thus achieving a total cost reduction in CPU time. The second contribution is a new method for making an explicit time-domain finite-element method (TDFEM) unconditionally stable for general electromagnetic analysis. In this method, for a given time step, we find the unstable modes that are the root cause of instability, and deduct them directly from the system matrix resulting from a TDFEM based analysis. As a result, an explicit TDFEM simulation is made stable for an arbitrarily large time step irrespective of the space step. The third contribution is a new method for full-wave applications from low to very high frequencies in a TDFEM based on matrix exponential. In this method, we directly deduct the eigenmodes having large eigenvalues from the system matrix, thus achieving a significantly increased time step in the matrix exponential based TDFEM. The fourth contribution is a new method for transforming the indefinite system matrix of a frequency-domain FEM to a symmetric positive definite one. We deduct non-positive definite component directly from the system matrix resulting from a frequency-domain FEM-based analysis. The resulting new representation of the finite-element operator ensures an iterative solution to converge in a small number of iterations. We then add back the non-positive definite component to synthesize the original solution with negligible cost

Explicit local time-stepping methods for time-dependent wave propagation

Semi-discrete Galerkin formulations of transient wave equations, either with conforming or discontinuous Galerkin finite element discretizations, typically lead to large systems of ordinary differential equations. When explicit time integration is used, the time-step is constrained by the smallest elements in the mesh for numerical stability, possibly a high price to pay. To overcome that overly restrictive stability constraint on the time-step, yet without resorting to implicit methods, explicit local time-stepping schemes (LTS) are presented here for transient wave equations either with or without damping. In the undamped case, leap-frog based LTS methods lead to high-order explicit LTS schemes, which conserve the energy. In the damped case, when energy is no longer conserved, Adams-Bashforth based LTS methods also lead to explicit LTS schemes of arbitrarily high accuracy. When combined with a finite element discretization in space with an essentially diagonal mass matrix, the resulting time-marching schemes are fully explicit and thus inherently parallel. Numerical experiments with continuous and discontinuous Galerkin finite element discretizations validate the theory and illustrate the usefulness of these local time-stepping methods.Comment: overview paper, typos added, references updated. arXiv admin note: substantial text overlap with arXiv:1109.448