6,061 research outputs found
Hopf algebras in dynamical systems theory
The theory of exact and of approximate solutions for non-autonomous linear
differential equations forms a wide field with strong ties to physics and
applied problems. This paper is meant as a stepping stone for an exploration of
this long-established theme, through the tinted glasses of a (Hopf and
Rota-Baxter) algebraic point of view. By reviewing, reformulating and
strengthening known results, we give evidence for the claim that the use of
Hopf algebra allows for a refined analysis of differential equations. We
revisit the renowned Campbell-Baker-Hausdorff-Dynkin formula by the modern
approach involving Lie idempotents. Approximate solutions to differential
equations involve, on the one hand, series of iterated integrals solving the
corresponding integral equations; on the other hand, exponential solutions.
Equating those solutions yields identities among products of iterated Riemann
integrals. Now, the Riemann integral satisfies the integration-by-parts rule
with the Leibniz rule for derivations as its partner; and skewderivations
generalize derivations. Thus we seek an algebraic theory of integration, with
the Rota-Baxter relation replacing the classical rule. The methods to deal with
noncommutativity are especially highlighted. We find new identities, allowing
for an extensive embedding of Dyson-Chen series of time- or path-ordered
products (of generalized integration operators); of the corresponding Magnus
expansion; and of their relations, into the unified algebraic setting of
Rota-Baxter maps and their inverse skewderivations. This picture clarifies the
approximate solutions to generalized integral equations corresponding to
non-autonomous linear (skew)differential equations.Comment: International Journal of Geometric Methods in Modern Physics, in
pres
Algebraic Structures and Stochastic Differential Equations driven by Levy processes
We construct an efficient integrator for stochastic differential systems
driven by Levy processes. An efficient integrator is a strong approximation
that is more accurate than the corresponding stochastic Taylor approximation,
to all orders and independent of the governing vector fields. This holds
provided the driving processes possess moments of all orders and the vector
fields are sufficiently smooth. Moreover the efficient integrator in question
is optimal within a broad class of perturbations for half-integer global root
mean-square orders of convergence. We obtain these results using the
quasi-shuffle algebra of multiple iterated integrals of independent Levy
processes.Comment: 41 pages, 11 figure
The algebra of the box spline
In this paper we want to revisit results of Dahmen and Micchelli on
box-splines which we reinterpret and make more precise. We compare these ideas
with the work of Brion, Szenes, Vergne and others on polytopes and partition
functions.Comment: 69 page
Hydrodynamic Flows on Curved Surfaces: Spectral Numerical Methods for Radial Manifold Shapes
We formulate hydrodynamic equations and spectrally accurate numerical methods
for investigating the role of geometry in flows within two-dimensional fluid
interfaces. To achieve numerical approximations having high precision and level
of symmetry for radial manifold shapes, we develop spectral Galerkin methods
based on hyperinterpolation with Lebedev quadratures for -projection to
spherical harmonics. We demonstrate our methods by investigating hydrodynamic
responses as the surface geometry is varied. Relative to the case of a sphere,
we find significant changes can occur in the observed hydrodynamic flow
responses as exhibited by quantitative and topological transitions in the
structure of the flow. We present numerical results based on the
Rayleigh-Dissipation principle to gain further insights into these flow
responses. We investigate the roles played by the geometry especially
concerning the positive and negative Gaussian curvature of the interface. We
provide general approaches for taking geometric effects into account for
investigations of hydrodynamic phenomena within curved fluid interfaces.Comment: 14 figure
Einstein equations in the null quasi-spherical gauge III: numerical algorithms
We describe numerical techniques used in the construction of our 4th order
evolution for the full Einstein equations, and assess the accuracy of
representative solutions. The code is based on a null gauge with a
quasi-spherical radial coordinate, and simulates the interaction of a single
black hole with gravitational radiation. Techniques used include spherical
harmonic representations, convolution spline interpolation and filtering, and
an RK4 "method of lines" evolution. For sample initial data of "intermediate"
size (gravitational field with 19% of the black hole mass), the code is
accurate to 1 part in 10^5, until null time z=55 when the coordinate condition
breaks down.Comment: Latex, 38 pages, 29 figures (360Kb compressed
- …