595 research outputs found
Uniformly high order accurate essentially non-oscillatory schemes 3
In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear
Study of space environment fabrication and repair techniques Final technical summary report
Joining techniques for fabrication and repair of structures in extraterrestrial environmen
Numerical Methods for Multilattices
Among the efficient numerical methods based on atomistic models, the
quasicontinuum (QC) method has attracted growing interest in recent years. The
QC method was first developed for crystalline materials with Bravais lattice
and was later extended to multilattices (Tadmor et al, 1999). Another existing
numerical approach to modeling multilattices is homogenization. In the present
paper we review the existing numerical methods for multilattices and propose
another concurrent macro-to-micro method in the numerical homogenization
framework. We give a unified mathematical formulation of the new and the
existing methods and show their equivalence. We then consider extensions of the
proposed method to time-dependent problems and to random materials.Comment: 31 page
Field theory simulation of Abelian-Higgs cosmic string cusps
We have performed a lattice field theory simulation of cusps in Abelian-Higgs
cosmic strings. The results are in accord with the theory that the portion of
the strings which overlaps near the cusp is released as radiation. The radius
of the string cores which must touch to produce the evaporation is
approximately in natural units. In general, the modifications to the
string shape due to the cusp may produce many cusps later in the evolution of a
string loop, but these later cusps will be much smaller in magnitude and more
closely resemble kinks.Comment: 9 pages, RevTeX, 13 figures with eps
Numerical Modeling of Transient Wave Propagation for High Frequency NDT
Electromagnetic nondestructive testing (NDT) methods use frequencies ranging from low (dc) to high (microwave) frequencies [1]. Applications of numerical methods to model two- and three-dimensional low-frequency (dc or eddy current) nondestructive testing phenomena, where displacement currents can be omitted, were extensively examined, [2,3]. These are all interior boundary value problems. Finite element study of ultrasonic wave propagation and scattering in metals, which is again an interior boundary value problem, was recently reported in [4]. However, modeling of wave propagation for high-frequency NDT problems have not yet been attempted. Recently, finite difference methods in time domain have been successfully applied to solve transient electromagnetic wave propagation problems over the atmosphere and the ground [5], and time-dependent eddy current problems [6]. The method used here is a generalization of this work and is designed for numerical modeling of high-frequency electromagnetic wave propagation arising from nondestructive testing applications. The physical situation includes examination of the scattering effects by cracks inside a piece of material (especially dielectrics) or due to surface variations when the material is illuminated by a TM plane wave. This leads to an interface type problem
Absorbing boundary conditions for the Westervelt equation
The focus of this work is on the construction of a family of nonlinear
absorbing boundary conditions for the Westervelt equation in one and two space
dimensions. The principal ingredient used in the design of such conditions is
pseudo-differential calculus. This approach enables to develop high order
boundary conditions in a consistent way which are typically more accurate than
their low order analogs. Under the hypothesis of small initial data, we
establish local well-posedness for the Westervelt equation with the absorbing
boundary conditions. The performed numerical experiments illustrate the
efficiency of the proposed boundary conditions for different regimes of wave
propagation
Matching Conditions in Atomistic-Continuum Modeling of Materials
A new class of matching condition between the atomistic and continuum regions
is presented for the multi-scale modeling of crystals. They ensure the accurate
passage of large scale information between the atomistic and continuum regions
and at the same time minimize the reflection of phonons at the interface. These
matching conditions can be made adaptive if we choose appropriate weight
functions. Applications to dislocation dynamics and friction between
two-dimensional atomically flat crystal surfaces are described.Comment: 6 pages, 4 figure
Recommended from our members
Fast Hybrid Algorithms For High Frequency Scattering
This paper deals with numerical methods for high frequency wave scattering. It introduces a new hybrid technique that couples a directional fast multipole method for a subsection of it scattering surface to an asymptotic formulation over the rest of the scattering domain. The directional fast multipole method is new and highly efficient for the solution of the boundary integral formulation of a general scattering problem but it requires at least a few unknowns per wavelength on the boundary. The asymptotic method that was introduced by Bruno and collaborators requires much fewer unknowns. Oil the other hand the scattered field must have a simple structure. Hybridization of these two methods retains their best properties for the solution of the full problem, Numerical examples are given for the solution of the Helmholtz equation in two space dimensions.Mathematic
On the well posedness of the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's field equations
We give a well posed initial value formulation of the
Baumgarte-Shapiro-Shibata-Nakamura form of Einstein's equations with gauge
conditions given by a Bona-Masso like slicing condition for the lapse and a
frozen shift. This is achieved by introducing extra variables and recasting the
evolution equations into a first order symmetric hyperbolic system. We also
consider the presence of artificial boundaries and derive a set of boundary
conditions that guarantee that the resulting initial-boundary value problem is
well posed, though not necessarily compatible with the constraints. In the case
of dynamical gauge conditions for the lapse and shift we obtain a class of
evolution equations which are strongly hyperbolic and so yield well posed
initial value formulations
A Novel Long Range Spin Chain and Planar N=4 Super Yang-Mills
We probe the long-range spin chain approach to planar N=4 gauge theory at
high loop order. A recently employed hyperbolic spin chain invented by
Inozemtsev is suitable for the SU(2) subsector of the state space up to three
loops, but ceases to exhibit the conjectured thermodynamic scaling properties
at higher orders. We indicate how this may be bypassed while nevertheless
preserving integrability, and suggest the corresponding all-loop asymptotic
Bethe ansatz. We also propose the local part of the all-loop gauge transfer
matrix, leading to conjectures for the asymptotically exact formulae for all
local commuting charges. The ansatz is finally shown to be related to a
standard inhomogeneous spin chain. A comparison of our ansatz to semi-classical
string theory uncovers a detailed, non-perturbative agreement between the
corresponding expressions for the infinite tower of local charge densities.
However, the respective Bethe equations differ slightly, and we end by refining
and elaborating a previously proposed possible explanation for this
disagreement.Comment: 48 pages, 1 figure. v2, further results added: discussion of the
relationship to an inhomogeneous spin chain, normalization in sec 3 unified,
v3: minor mistakes corrected, published versio
- …