105 research outputs found
Existence and stability of viscoelastic shock profiles
We investigate existence and stability of viscoelastic shock profiles for a
class of planar models including the incompressible shear case studied by
Antman and Malek-Madani. We establish that the resulting equations fall into
the class of symmetrizable hyperbolic--parabolic systems, hence spectral
stability implies linearized and nonlinear stability with sharp rates of decay.
The new contributions are treatment of the compressible case, formulation of a
rigorous nonlinear stability theory, including verification of stability of
small-amplitude Lax shocks, and the systematic incorporation in our
investigations of numerical Evans function computations determining stability
of large-amplitude and or nonclassical type shock profiles.Comment: 43 pages, 12 figure
A note on the convergence of parametrised non-resonant invariant manifolds
Truncated Taylor series representations of invariant manifolds are abundant
in numerical computations. We present an aposteriori method to compute the
convergence radii and error estimates of analytic parametrisations of
non-resonant local invariant manifolds of a saddle of an analytic vector field,
from such a truncated series. This enables us to obtain local enclosures, as
well as existence results, for the invariant manifolds
A familial case of myoclonus-epilepsy of Unferricht-Lundborg
The aim of this work was to highlight the clinical features of a rare familial variant of Unferricht-Lundborg myoclonus-epilepsy in comparison with the manifestations of single cases described in the literature.Цель работы — выделение клинических особенностей редкого семейного варианта миоклонус-эпилепсии Унферрихта-Лундборга в сравнении с проявлениями одиночных случаев, описанными в литературе
An integral method for solving nonlinear eigenvalue problems
We propose a numerical method for computing all eigenvalues (and the
corresponding eigenvectors) of a nonlinear holomorphic eigenvalue problem that
lie within a given contour in the complex plane. The method uses complex
integrals of the resolvent operator, applied to at least column vectors,
where is the number of eigenvalues inside the contour. The theorem of
Keldysh is employed to show that the original nonlinear eigenvalue problem
reduces to a linear eigenvalue problem of dimension .
No initial approximations of eigenvalues and eigenvectors are needed. The
method is particularly suitable for moderately large eigenvalue problems where
is much smaller than the matrix dimension. We also give an extension of the
method to the case where is larger than the matrix dimension. The
quadrature errors caused by the trapezoid sum are discussed for the case of
analytic closed contours. Using well known techniques it is shown that the
error decays exponentially with an exponent given by the product of the number
of quadrature points and the minimal distance of the eigenvalues to the
contour
Analysis of the shearing instability in nonlinear convection and magnetoconvection
Numerical experiments on two-dimensional convection with or without a vertical magnetic field reveal a bewildering variety of periodic and aperiodic oscillations. Steady rolls can develop a shearing instability, in which rolls turning over in one direction grow at the expense of rolls turning over in the other, resulting in a net shear across the layer. As the temperature difference across the fluid is increased, two-dimensional pulsating waves occur, in which the direction of shear alternates. We analyse the nonlinear dynamics of this behaviour by first constructing appropriate low-order sets of ordinary differential equations, which show the same behaviour, and then analysing the global bifurcations that lead to these oscillations by constructing one-dimensional return maps. We compare the behaviour of the partial differential equations, the models and the maps in systematic two-parameter studies of both the magnetic and the non-magnetic cases, emphasising how the symmetries of periodic solutions change as a result of global bifurcations. Much of the interesting behaviour is associated with a discontinuous change in the leading direction of a fixed point at a global bifurcation; this change occurs when the magnetic field is introduced
Lin's method for heteroclinic chains involving periodic orbits
We present an extension of the theory known as Lin's method to heteroclinic
chains that connect hyperbolic equilibria and hyperbolic periodic orbits. Based
on the construction of a so-called Lin orbit, that is, a sequence of continuous
partial orbits that only have jumps in a certain prescribed linear subspace,
estimates for these jumps are derived. We use the jump estimates to discuss
bifurcation equations for homoclinic orbits near heteroclinic cycles between an
equilibrium and a periodic orbit (EtoP cycles)
Equation-Free Analysis of Macroscopic Behavior in Traffic and Pedestrian Flow
Equation-free methods make possible an analysis of the evolution of a few
coarse-grained or macroscopic quantities for a detailed and realistic model
with a large number of fine-grained or microscopic variables, even though no
equations are explicitly given on the macroscopic level. This will facilitate a
study of how the model behavior depends on parameter values including an
understanding of transitions between different types of qualitative behavior.
These methods are introduced and explained for traffic jam formation and
emergence of oscillatory pedestrian counter flow in a corridor with a narrow
door
Three-dimensional coherent X-ray diffraction imaging of a ceramic nanofoam: determination of structural deformation mechanisms
Ultra-low density polymers, metals, and ceramic nanofoams are valued for
their high strength-to-weight ratio, high surface area and insulating
properties ascribed to their structural geometry. We obtain the labrynthine
internal structure of a tantalum oxide nanofoam by X-ray diffractive imaging.
Finite element analysis from the structure reveals mechanical properties
consistent with bulk samples and with a diffusion limited cluster aggregation
model, while excess mass on the nodes discounts the dangling fragments
hypothesis of percolation theory.Comment: 8 pages, 5 figures, 30 reference
Computing the common zeros of two bivariate functions via Bézout resultants
The common zeros of two bivariate functions can be computed by finding the common zeros of their polynomial interpolants expressed in a tensor Chebyshev basis. From here we develop a bivariate rootfinding algorithm based on the hidden variable resultant method and Bézout matrices with polynomial entries. Using techniques including domain subdivision, Bézoutian regularization, and local refinement we are able to reliably and accurately compute the simple common zeros of two smooth functions with polynomial interpolants of very high degree (≥ 1000). We analyze the resultant method and its conditioning by noting that the Bézout matrices are matrix polynomials. Two implementations are available: one on the Matlab Central File Exchange and another in the roots command in Chebfun2 that is adapted to suit Chebfun’s methodology
- …