An analysis of over-relaxation in a kinetic approximation of systems of conservation laws

The over-relaxation approach is an alternative to the Jin–Xin relaxation method in order to apply the equilibrium source term in a more precise way. This is also a key ingredient of the lattice Boltzmann method for achieving second-order accuracy. In this work, we provide an analysis of the over-relaxation kinetic scheme. We compute its equivalent equation, which is particularly useful for devising stable boundary conditions for the hidden kinetic variables.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

Restrictions of Quadratic Forms to Lagrangian Planes,Quadratic Matrix Equations,and Gyroscopic Stabilization

We discuss the symplectic geometry of linear Hamiltonian systems with nondegenerate Hamiltonians. These systems can be reduced to linear second-order differential equations characteristic of linear oscillation theory. This reduction is related to the problem on the signatures of restrictions of quadratic forms to Lagrangian planes. We study vortex symplectic planes invariant with respect to linear Hamiltonian systems. These planes are determined by the solutions of quadratic matrix equations of a special form. New conditions for gyroscopic stabilization are found.     

A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

带有给定凸切线多边形的保形五次样条逼近

本讨论带有给定切线多边形的保形逼近问题．给出了一条与给定切线多边形相切的保形五次参数祥条曲线。     

Dirac constraint analysis and symplectic structure of anti-self-dual Yang-Mills equations

We present the explicit form of the symplectic structure of anti-self-dual Yang-Mills (ASDYM) equations in Yang’s J- and K-gauges in order to establish the bi-Hamiltonian structure of this completely integrable system. Dirac’s theory of constraints is applied to the degenerate Lagrangians that yield the ASDYM equations. The constraints are second class as in the case of all completely integrable systems which stands in sharp contrast to the situation in full Yang-Mills theory. We construct the Dirac brackets and the symplectic 2-forms for both J- and K-gauges. The covariant symplectic structure of ASDYM equations is obtained using the Witten-Zuckerman formalism. We show that the appropriate component of the Witten-Zuckerman closed and conserved 2-form vector density reduces to the symplectic 2-form obtained from Dirac’s theory. Finally, we present the Bäcklund transformation between the J- and K-gauges in order to apply Magri’s theorem to the respective two Hamiltonian structures.     

On Covariant Phase Space and the Variational Bicomplex

The notion of a phase space in classical mechanics is well known. The extension of this concept to field theory however, is a challenging endeavor, and over the years numerous proposals for such a generalization have appeared in the literature. In this paper We review a Hamiltonian formulation of Lagrangian field theory based on an extension to infinite dimensions of J.-M. Souriau's symplectic approach to mechanics. Following G. Zuckerman, we state our results in terms of the modern geometric theory of differential equations and the variational bicomplex. As an elementary example, we construct a phase space for the Monge–Ampere equation.     

Shape preservation regions for six-dimensional spaces

We analyze the critical length for design purposes of six-dimensional spaces invariant under translations and reflections containing the functions 1, cos t and sin t. These spaces also contain the first degree polynomials as well as trigonometric and/or hyperbolic functions. We identify the spaces whose critical length for design purposes is greater than 2π and find its maximum 4π. By a change of variables, two biparametric families of spaces arise. We call shape preservation region to the set of admissible parameters in order that the space has shape preserving representations for curves. We describe the shape preserving regions for both families. To our friend Mariano Gasca on occasion of his 60th birthday Research partially supported by the Spanish Research Grant MTM2006-03388, by Gobierno de Aragón and Fondo Social Europeo.