Finite-Difference and Pseudospectral Time-Domain Methods Applied to Backwards-Wave Metamaterials

Backwards-wave (BW) materials that have simultaneously negative real parts of their electric permittivity and magnetic permeability can support waves where phase and power propagation occur in opposite directions. These materials were predicted to have many unusual electromagnetic properties, among them amplification of the near-field of a point source, which could lead to the perfect reconstruction of the source field in an image [J. Pendry, Phys. Rev. Lett. \textbf{85}, 3966 (2000)]. Often systems containing BW materials are simulated using the finite-difference time-domain technique. We show that this technique suffers from a numerical artifact due to its staggered grid that makes its use in simulations involving BW materials problematic. The pseudospectral time-domain technique, on the other hand, uses a collocated grid and is free of this artifact. It is also shown that when modeling the dispersive BW material, the linear frequency approximation method introduces error that affects the frequency of vanishing reflection, while the auxiliary differential equation, the Z transform, and the bilinear frequency approximation method produce vanishing reflection at the correct frequency. The case of vanishing reflection is of particular interest for field reconstruction in imaging applications.Comment: 9 pages, 8 figures, accepted by IEEE Transactions on Antennas and Propagatio

Numerical Boundary Conditions for Schemes with Centered and Biased Differences in Subsonic Gas Dynamics

For finite-difference schemes of EBR class in multi-dimensional inviscid gas dynamics, the nonreflecting boundary conditions are analyzed and developed. The wave-reflection properties of discrete models differ significantly from each other and from the continuous Euler equations. For certain schemes there exist local boundary conditions which lead to small reflections of waves with arbitrary incidence angle. Numerical examples are shown both linear and nonlinear

Fast and Accurate Computation of Time-Domain Acoustic Scattering Problems with Exact Nonreflecting Boundary Conditions

This paper is concerned with fast and accurate computation of exterior wave equations truncated via exact circular or spherical nonreflecting boundary conditions (NRBCs, which are known to be nonlocal in both time and space). We first derive analytic expressions for the underlying convolution kernels, which allow for a rapid and accurate evaluation of the convolution with $O(N_t)$ operations over $N_t$ successive time steps. To handle the onlocality in space, we introduce the notion of boundary perturbation, which enables us to handle general bounded scatters by solving a sequence of wave equations in a regular domain. We propose an efficient spectral-Galerkin solver with Newmark's time integration for the truncated wave equation in the regular domain. We also provide ample numerical results to show high-order accuracy of NRBCs and efficiency of the proposed scheme.Comment: 22 pages with 9 figure

Boundary conditions for coupled quasilinear wave equations with application to isolated systems

We consider the initial-boundary value problem for systems of quasilinear wave equations on domains of the form $[0,T] \times \Sigma$, where $\Sigma$ is a compact manifold with smooth boundaries $\partial\Sigma$. By using an appropriate reduction to a first order symmetric hyperbolic system with maximal dissipative boundary conditions, well posedness of such problems is established for a large class of boundary conditions on $\partial\Sigma$. We show that our class of boundary conditions is sufficiently general to allow for a well posed formulation for different wave problems in the presence of constraints and artificial, nonreflecting boundaries, including Maxwell's equations in the Lorentz gauge and Einstein's gravitational equations in harmonic coordinates. Our results should also be useful for obtaining stable finite-difference discretizations for such problems.Comment: 22 pages, no figure

Acoustic Saturation in Bubbly Cavitating Flow Adjacent to an Oscillating Wall

Bubbly cavitating flow generated by the normal oscillation of a wall bounding a semi-infinite domain of fluid is computed using a continuum two-phase flow model. Bubble dynamics are computed, on the microscale, using the Rayleigh-Plesset equation. A Lagrangian finite volume scheme and implicit adaptive time marching are employed to accurately resolve bubbly shock waves and other steep gradients in the flow. The one-dimensional, unsteady computations show that when the wall oscillation frequency is much smaller than the bubble natural frequency, the power radiated away from the wall is limited by an acoustic saturation effect (the radiated power becomes independent of the amplitude of vibration), which is similar to that found in a pure gas. That is, for large enough vibration amplitude, nonlinear steepening of the generated waves leads to shocking of the wave train, and the dissipation associated with the jump conditions across each shock limits the radiated power. In the model, damping of the bubble volume oscillations is restricted to a simple "effective" viscosity. For wall oscillation frequency less than the bubble natural frequency, the saturation amplitude of the radiated field is nearly independent of any specific damping mechanism. Finally, implications for noise radiation from cavitating flows are discussed