Nilpotent normal form for divergence-free vector fields and volume-preserving maps

We study the normal forms for incompressible flows and maps in the neighborhood of an equilibrium or fixed point with a triple eigenvalue. We prove that when a divergence free vector field in $\mathbb{R}^3$ has nilpotent linearization with maximal Jordan block then, to arbitrary degree, coordinates can be chosen so that the nonlinear terms occur as a single function of two variables in the third component. The analogue for volume-preserving diffeomorphisms gives an optimal normal form in which the truncation of the normal form at any degree gives an exactly volume-preserving map whose inverse is also polynomial inverse with the same degree.Comment: laTeX, 20 pages, 1 figur

Nilpotent normal form for divergence-free vector fields and volume-preserving maps

We study the normal forms for incompressible flows and maps in the neighborhood of an equilibrium or fixed point with a triple eigenvalue. We prove that when a divergence free vector field in $\mathbb{R}^3$ has nilpotent linearization with maximal Jordan block then, to arbitrary degree, coordinates can be chosen so that the nonlinear terms occur as a single function of two variables in the third component. The analogue for volume-preserving diffeomorphisms gives an optimal normal form in which the truncation of the normal form at any degree gives an exactly volume-preserving map whose inverse is also polynomial inverse with the same degree.Comment: laTeX, 20 pages, 1 figur

Palindromic 3-stage splitting integrators, a roadmap

The implementation of multi-stage splitting integrators is essentially the same as the implementation of the familiar Strang/Verlet method. Therefore multi-stage formulas may be easily incorporated into software that now uses the Strang/Verlet integrator. We study in detail the two-parameter family of palindromic, three-stage splitting formulas and identify choices of parameters that may outperform the Strang/Verlet method. One of these choices leads to a method of effective order four suitable to integrate in time some partial differential equations. Other choices may be seen as perturbations of the Strang method that increase efficiency in molecular dynamics simulations and in Hybrid Monte Carlo sampling.Comment: 20 pages, 8 figures, 2 table