742 research outputs found
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
operations over 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
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
On-surface radiation condition for multiple scattering of waves
The formulation of the on-surface radiation condition (OSRC) is extended to
handle wave scattering problems in the presence of multiple obstacles. The new
multiple-OSRC simultaneously accounts for the outgoing behavior of the wave
fields, as well as, the multiple wave reflections between the obstacles. Like
boundary integral equations (BIE), this method leads to a reduction in
dimensionality (from volume to surface) of the discretization region. However,
as opposed to BIE, the proposed technique leads to boundary integral equations
with smooth kernels. Hence, these Fredholm integral equations can be handled
accurately and robustly with standard numerical approaches without the need to
remove singularities. Moreover, under weak scattering conditions, this approach
renders a convergent iterative method which bypasses the need to solve single
scattering problems at each iteration.
Inherited from the original OSRC, the proposed multiple-OSRC is generally a
crude approximate method. If accuracy is not satisfactory, this approach may
serve as a good initial guess or as an inexpensive pre-conditioner for Krylov
iterative solutions of BIE
Unsteady-flow-field predictions for oscillating cascades
The unsteady flow field around an oscillating cascade of flat plates with zero stagger was studied by using a time marching Euler code. This case had an exact solution based on linear theory and served as a model problem for studying pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step size was shown for a moderate reduced frequency. Results show that an approximate nonreflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer to the airfoils than when reflective boundaries are used. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates captures the unsteady flow field better than when uniform grids are used as long as the 'Courant Friedrichs Levy' (CFL) number is less than 1 for a sufficient portion of the grid. Finally, a solution based on an optimization of grid, CFL number, and boundary conditions shows good agreement with linear theory
XTRAN2L: A program for solving the general-frequency unsteady transonic small disturbance equation
A program, XTRAN2L, for solving the general-frequency unsteady transonic small disturbance potential equation was developed. It is a modification of the LTRAN2-NLR code. The alternating-direction-implicit (ADI) method of Rizzetta and Chin is used to advance solutions of the potential equation in time Engquist-Osher monotone spatial differencing is used in the ADI solution algorithm. As a result, the XTRAN2L code is more robust and more efficient than similar codes that use Murman-Cole type-dependent spatial differencing. Nonreflecting boundary conditions that are consistent with the general-frequency equation have been developed and implemented at the far-field boundaries. Use of those conditions allow the computational boundaries to be moved closer to the airfoil with no loss of accuracy. This makes the XTRAN2L code more economical to use
Analytic structure of radiation boundary kernels for blackhole perturbations
Exact outer boundary conditions for gravitational perturbations of the
Schwarzschild metric feature integral convolution between a time-domain
boundary kernel and each radiative mode of the perturbation. For both axial
(Regge-Wheeler) and polar (Zerilli) perturbations, we study the Laplace
transform of such kernels as an analytic function of (dimensionless) Laplace
frequency. We present numerical evidence indicating that each such
frequency-domain boundary kernel admits a "sum-of-poles" representation. Our
work has been inspired by Alpert, Greengard, and Hagstrom's analysis of
nonreflecting boundary conditions for the ordinary scalar wave equation.Comment: revtex4, 14 pages, 12 figures, 3 table
- …