771 research outputs found
The Erpenbeck high frequency instability theorem for ZND detonations
The rigorous study of spectral stability for strong detonations was begun by
J.J. Erpenbeck in [Er1]. Working with the Zeldovitch-von Neumann-D\"oring (ZND)
model, which assumes a finite reaction rate but ignores effects like viscosity
corresponding to second order derivatives, he used a normal mode analysis to
define a stability function V(\tau,\eps) whose zeros in
correspond to multidimensional perturbations of a steady detonation profile
that grow exponentially in time. Later in a remarkable paper [Er3] he provided
strong evidence, by a combination of formal and rigorous arguments, that for
certain classes of steady ZND profiles, unstable zeros of exist for
perturbations of sufficiently large transverse wavenumber \eps, even when the
von Neumann shock, regarded as a gas dynamical shock, is uniformly stable in
the sense defined (nearly twenty years later) by Majda. In spite of a great
deal of later numerical work devoted to computing the zeros of V(\tau,\eps),
the paper \cite{Er3} remains the only work we know of that presents a detailed
and convincing theoretical argument for detecting them.
The analysis in [Er3] points the way toward, but does not constitute, a
mathematical proof that such unstable zeros exist. In this paper we identify
the mathematical issues left unresolved in [Er3] and provide proofs, together
with certain simplifications and extensions, of the main conclusions about
stability and instability of detonations contained in that paper.
The main mathematical problem, and our principal focus here, is to determine
the precise asymptotic behavior as \eps\to \infty of solutions to a linear
system of ODEs in , depending on \eps and a complex frequency as
parameters, with turning points on the half-line
Reliable Computation of the Zeros of Solutions of Second Order Linear ODEs Using a Fourth Order Method
A fourth order fixed point method to compute the zeros of solutions of second order
homogeneous linear ODEs is obtained from the approximate integration of the Riccati equation
associated with the ODE. The method requires the evaluation of the logarithmic derivative of the
function and also uses the coefficients of the ODE. An algorithm to compute with certainty all the
zeros in an interval is given which provides a fast, reliable, and accurate method of computation.
The method is illustrated by the computation of the zeros of Gauss hypergeometric functions (including
Jacobi polynomials) and confluent hypergeometric functions (Laguerre polynomials, Hermite
polynomials, and Bessel functions included) among others. The examples show that typically 4 or 5
iterations per root are enough to provide more than 100 digits of accuracy, without requiring a priori
estimations of the roots
Multidomain Spectral Method for the Helically Reduced Wave Equation
We consider the 2+1 and 3+1 scalar wave equations reduced via a helical
Killing field, respectively referred to as the 2-dimensional and 3-dimensional
helically reduced wave equation (HRWE). The HRWE serves as the fundamental
model for the mixed-type PDE arising in the periodic standing wave (PSW)
approximation to binary inspiral. We present a method for solving the equation
based on domain decomposition and spectral approximation. Beyond describing
such a numerical method for solving strictly linear HRWE, we also present
results for a nonlinear scalar model of binary inspiral. The PSW approximation
has already been theoretically and numerically studied in the context of the
post-Minkowskian gravitational field, with numerical simulations carried out
via the "eigenspectral method." Despite its name, the eigenspectral technique
does feature a finite-difference component, and is lower-order accurate. We
intend to apply the numerical method described here to the theoretically
well-developed post-Minkowski PSW formalism with the twin goals of spectral
accuracy and the coordinate flexibility afforded by global spectral
interpolation.Comment: 57 pages, 11 figures, uses elsart.cls. Final version includes
revisions based on referee reports and has two extra figure
Basic Methods for Computing Special Functions
This paper gives an overview of methods for the numerical evaluation of special functions, that is, the functions that arise in many problems from mathematical physics, engineering, probability theory, and other applied sciences. We consider in detail a selection of basic methods which are
frequently used in the numerical evaluation of special functions: converging and asymptotic series, including Chebyshev expansions, linear recurrence relations, and numerical quadrature. Several other methods are available and some of these will be discussed in less detail. We give examples of recent software for special functions where these methods are used. We mention a list of new publications on computational aspects of special functions available on our website
Computational Methods for Nonlinear Systems Analysis With Applications in Mathematics and Engineering
An investigation into current methods and new approaches for solving systems of nonlinear equations was performed. Nontraditional methods for implementing arc-length type solvers were developed in search of a more robust capability for solving general systems of nonlinear algebraic equations. Processes for construction of parameterized curves representing the many possible solutions to systems of equations versus finding single or point solutions were established. A procedure based on these methods was then developed to identify static equilibrium states for solutions to multi-body-dynamic systems. This methodology provided for a pictorial of the overall solution to a given system, which demonstrated the possibility of multiple candidate equilibrium states for which a procedure for selection of the proper state was proposed. Arc-length solvers were found to identify and more readily trace solution curves as compared to other solvers making such an approach practical. Comparison of proposed methods was made to existing methods found in the literature and commercial software with favorable results. Finally, means for parallel processing of the Jacobian matrix inherent to the arc-length and other nonlinear solvers were investigated, and an efficient approach for implementation was identified. Several case studies were performed to substantiate results. Commercial software was also used in some instances for additional results verification
- …