185 research outputs found
R-adaptive multisymplectic and variational integrators
Moving mesh methods (also called r-adaptive methods) are space-adaptive
strategies used for the numerical simulation of time-dependent partial
differential equations. These methods keep the total number of mesh points
fixed during the simulation, but redistribute them over time to follow the
areas where a higher mesh point density is required. There are a very limited
number of moving mesh methods designed for solving field-theoretic partial
differential equations, and the numerical analysis of the resulting schemes is
challenging. In this paper we present two ways to construct r-adaptive
variational and multisymplectic integrators for (1+1)-dimensional Lagrangian
field theories. The first method uses a variational discretization of the
physical equations and the mesh equations are then coupled in a way typical of
the existing r-adaptive schemes. The second method treats the mesh points as
pseudo-particles and incorporates their dynamics directly into the variational
principle. A user-specified adaptation strategy is then enforced through
Lagrange multipliers as a constraint on the dynamics of both the physical field
and the mesh points. We discuss the advantages and limitations of our methods.
Numerical results for the Sine-Gordon equation are also presented.Comment: 65 pages, 13 figure
Computational Engineering
This Workshop treated a variety of finite element methods and applications in computational engineering and expanded their mathematical foundation in engineering analysis. Among the 53 participants were mathematicians and engineers with focus on mixed and nonstandard finite element schemes and their applications
Entropies from coarse-graining: convex polytopes vs. ellipsoids
We examine the Boltzmann/Gibbs/Shannon and the
non-additive Havrda-Charv\'{a}t / Dar\'{o}czy/Cressie-Read/Tsallis \
\ and the Kaniadakis -entropy \ \
from the viewpoint of coarse-graining, symplectic capacities and convexity. We
argue that the functional form of such entropies can be ascribed to a
discordance in phase-space coarse-graining between two generally different
approaches: the Euclidean/Riemannian metric one that reflects independence and
picks cubes as the fundamental cells and the symplectic/canonical one that
picks spheres/ellipsoids for this role. Our discussion is motivated by and
confined to the behaviour of Hamiltonian systems of many degrees of freedom. We
see that Dvoretzky's theorem provides asymptotic estimates for the minimal
dimension beyond which these two approaches are close to each other. We state
and speculate about the role that dualities may play in this viewpoint.Comment: 63 pages. No figures. Standard LaTe
The Magnus expansion and some of its applications
Approximate resolution of linear systems of differential equations with
varying coefficients is a recurrent problem shared by a number of scientific
and engineering areas, ranging from Quantum Mechanics to Control Theory. When
formulated in operator or matrix form, the Magnus expansion furnishes an
elegant setting to built up approximate exponential representations of the
solution of the system. It provides a power series expansion for the
corresponding exponent and is sometimes referred to as Time-Dependent
Exponential Perturbation Theory. Every Magnus approximant corresponds in
Perturbation Theory to a partial re-summation of infinite terms with the
important additional property of preserving at any order certain symmetries of
the exact solution. The goal of this review is threefold. First, to collect a
number of developments scattered through half a century of scientific
literature on Magnus expansion. They concern the methods for the generation of
terms in the expansion, estimates of the radius of convergence of the series,
generalizations and related non-perturbative expansions. Second, to provide a
bridge with its implementation as generator of especial purpose numerical
integration methods, a field of intense activity during the last decade. Third,
to illustrate with examples the kind of results one can expect from Magnus
expansion in comparison with those from both perturbative schemes and standard
numerical integrators. We buttress this issue with a revision of the wide range
of physical applications found by Magnus expansion in the literature.Comment: Report on the Magnus expansion for differential equations and its
applications to several physical problem
Lectures on Mechanics
Publisher's description: The use of geometric methods in classical mechanics has proven fruitful, with wide applications in physics and engineering. In this book, Professor Marsden concentrates on these geometric aspects, especially on symmetry techniques. The main points he covers are: the stability of relative equilibria, which is analyzed using the block diagonalization technique; geometric phases, studied using the reduction and reconstruction technique; and bifurcation of relative equilibria and chaos in mechanical systems. A unifying theme for these points is provided by reduction theory, the associated mechanical connection and techniques from dynamical systems. These methods can be applied to many control and stabilization situations, and this is illustrated using rigid bodies with internal rotors, and the use of geometric phases in mechanical systems. To illustrate the above ideas and the power of geometric arguments, the author studies a variety of specific systems, including the double spherical pendulum and the classical rotating water molecule
Novel Discretization Schemes for the Numerical Simulation of Membrane Dynamics
Motivated by the demands of simulating flapping wings of Micro Air Vehicles, novel numerical methods were developed and evaluated for the dynamic simulation of membranes. For linear membranes, a mixed-form time-continuous Galerkin method was employed using trilinear space-time elements, and the entire space-time domain was discretized and solved simultaneously. For geometrically nonlinear membranes, the model incorporated two new schemes that were independently developed and evaluated. Time marching was performed using quintic Hermite polynomials uniquely determined by end-point jerk constraints. The single-step, implicit scheme was significantly more accurate than the most common Newmark schemes. For a simple harmonic oscillator, the scheme was found to be symplectic, frequency-preserving, and conditionally stable. Time step size was limited by accuracy requirements rather than stability. The spatial discretization scheme employed a staggered grid, grouping of nonlinear terms, and polygon shape functions in a strong-form point collocation formulation. Validation against existing experimental data showed the method to be accurate until hyperelastic effects dominate
- …