15,225 research outputs found
Two-component {CH} system: Inverse Scattering, Peakons and Geometry
An inverse scattering transform method corresponding to a Riemann-Hilbert
problem is formulated for CH2, the two-component generalization of the
Camassa-Holm (CH) equation. As an illustration of the method, the multi -
soliton solutions corresponding to the reflectionless potentials are
constructed in terms of the scattering data for CH2.Comment: 22 pages, 3 figures, draft, please send comment
Continuous and discrete Clebsch variational principles
The Clebsch method provides a unifying approach for deriving variational
principles for continuous and discrete dynamical systems where elements of a
vector space are used to control dynamics on the cotangent bundle of a Lie
group \emph{via} a velocity map. This paper proves a reduction theorem which
states that the canonical variables on the Lie group can be eliminated, if and
only if the velocity map is a Lie algebra action, thereby producing the
Euler-Poincar\'e (EP) equation for the vector space variables. In this case,
the map from the canonical variables on the Lie group to the vector space is
the standard momentum map defined using the diamond operator. We apply the
Clebsch method in examples of the rotating rigid body and the incompressible
Euler equations. Along the way, we explain how singular solutions of the EP
equation for the diffeomorphism group (EPDiff) arise as momentum maps in the
Clebsch approach. In the case of finite dimensional Lie groups, the Clebsch
variational principle is discretised to produce a variational integrator for
the dynamical system. We obtain a discrete map from which the variables on the
cotangent bundle of a Lie group may be eliminated to produce a discrete EP
equation for elements of the vector space. We give an integrator for the
rotating rigid body as an example. We also briefly discuss how to discretise
infinite-dimensional Clebsch systems, so as to produce conservative numerical
methods for fluid dynamics
An Optimal Control Formulation for Inviscid Incompressible Ideal Fluid Flow
In this paper we consider the Hamiltonian formulation of the equations of
incompressible ideal fluid flow from the point of view of optimal control
theory. The equations are compared to the finite symmetric rigid body equations
analyzed earlier by the authors. We discuss various aspects of the Hamiltonian
structure of the Euler equations and show in particular that the optimal
control approach leads to a standard formulation of the Euler equations -- the
so-called impulse equations in their Lagrangian form. We discuss various other
aspects of the Euler equations from a pedagogical point of view. We show that
the Hamiltonian in the maximum principle is given by the pairing of the
Eulerian impulse density with the velocity. We provide a comparative discussion
of the flow equations in their Eulerian and Lagrangian form and describe how
these forms occur naturally in the context of optimal control. We demonstrate
that the extremal equations corresponding to the optimal control problem for
the flow have a natural canonical symplectic structure.Comment: 6 pages, no figures. To appear in Proceedings of the 39th IEEEE
Conference on Decision and Contro
The Hamiltonian structure and Euler-Poincar\'{e} formulation of the Vlasov-Maxwell and gyrokinetic systems
We present a new variational principle for the gyrokinetic system, similar to
the Maxwell-Vlasov action presented in Ref. 1. The variational principle is in
the Eulerian frame and based on constrained variations of the phase space fluid
velocity and particle distribution function. Using a Legendre transform, we
explicitly derive the field theoretic Hamiltonian structure of the system. This
is carried out with a modified Dirac theory of constraints, which is used to
construct meaningful brackets from those obtained directly from
Euler-Poincar\'{e} theory. Possible applications of these formulations include
continuum geometric integration techniques, large-eddy simulation models and
Casimir type stability methods.
[1] H. Cendra et. al., Journal of Mathematical Physics 39, 3138 (1998)Comment: 36 pages, 1 figur
Kinetic and ion pairing contributions in the dielectric spectra of electrolyte aqueous solutions
Understanding dielectric spectra can reveal important information about the
dynamics of solvents and solutes from the dipolar relaxation times down to
electronic ones. In the late 1970s, Hubbard and Onsager predicted that adding
salt ions to a polar solution would result in a reduced dielectric permittivity
that arises from the unexpected tendency of solvent dipoles to align opposite
to the applied field. So far, this effect has escaped an experimental
verification, mainly because of the concomitant appearance of dielectric
saturation from which the Hubbard-Onsager decrement cannot be easily separated.
Here we develop a novel non-equilibrium molecular dynamics simulation approach
to determine this decrement accurately for the first time. Using a
thermodynamic consistent all-atom force field we show that for an aqueous
solution containing sodium chloride around 4.8 Mol/l, this effect accounts for
12\% of the total dielectric permittivity. The dielectric decrement can be
strikingly different if a less accurate force field for the ions is used. Using
the widespread GROMOS parameters, we observe in fact an {\it increment} of the
dielectric permittivity rather than a decrement. We can show that this
increment is caused by ion pairing, introduced by a too low dispersion force,
and clarify the microscopic connection between long-living ion pairs and the
appearance of specific features in the dielectric spectrum of the solution
An Integrable Shallow Water Equation with Linear and Nonlinear Dispersion
We study a class of 1+1 quadratically nonlinear water wave equations that
combines the linear dispersion of the Korteweg-deVries (KdV) equation with the
nonlinear/nonlocal dispersion of the Camassa-Holm (CH) equation, yet still
preserves integrability via the inverse scattering transform (IST) method.
This IST-integrable class of equations contains both the KdV equation and the
CH equation as limiting cases. It arises as the compatibility condition for a
second order isospectral eigenvalue problem and a first order equation for the
evolution of its eigenfunctions. This integrable equation is shown to be a
shallow water wave equation derived by asymptotic expansion at one order higher
approximation than KdV. We compare its traveling wave solutions to KdV
solitons.Comment: 4 pages, no figure
The optimal P3M algorithm for computing electrostatic energies in periodic systems
We optimize Hockney and Eastwood's Particle-Particle Particle-Mesh (P3M)
algorithm to achieve maximal accuracy in the electrostatic energies (instead of
forces) in 3D periodic charged systems. To this end we construct an optimal
influence function that minimizes the RMS errors in the energies. As a
by-product we derive a new real-space cut-off correction term, give a
transparent derivation of the systematic errors in terms of Madelung energies,
and provide an accurate analytical estimate for the RMS error of the energies.
This error estimate is a useful indicator of the accuracy of the computed
energies, and allows an easy and precise determination of the optimal values of
the various parameters in the algorithm (Ewald splitting parameter, mesh size
and charge assignment order).Comment: 31 pages, 3 figure
Recommended from our members
Solving linear programs without breaking abstractions
We show that the ellipsoid method for solving linear programs can be implemented in a way that respects the symmetry of the program being solved. That is to say, there is an algorithmic implementation of the method that does not distinguish, or make choices, between variables or constraints in the program unless they are distinguished by properties definable from the program. In particular, we demonstrate that the solvability of linear programs can be expressed in fixed-point logic with counting (FPC) as long as the program is given by a separation oracle that is itself definable in FPC. We use this to show that the size of a maximum matching in a graph is definable in FPC. This settles an open problem first posed by Blass, Gurevich and Shelah [Blass et al. 1999]. On the way to defining a suitable separation oracle for the maximum matching program, we provide FPC formulas defining canonical maximum flows and minimum cuts in undirected capacitated graphs.Research supported by EPSRC grant EP/H026835.This is the author accepted manuscript. The final version is available from ACM via http://dx.doi.org/10.1145/282289
NASA technology utilization survey on composite materials
NASA and NASA-funded contractor contributions to the field of composite materials are surveyed. Existing and potential non-aerospace applications of the newer composite materials are emphasized. Economic factors for selection of a composite for a particular application are weight savings, performance (high strength, high elastic modulus, low coefficient of expansion, heat resistance, corrosion resistance,), longer service life, and reduced maintenance. Applications for composites in agriculture, chemical and petrochemical industries, construction, consumer goods, machinery, power generation and distribution, transportation, biomedicine, and safety are presented. With the continuing trend toward further cost reductions, composites warrant consideration in a wide range of non-aerospace applications. Composite materials discussed include filamentary reinforced materials, laminates, multiphase alloys, solid multiphase lubricants, and multiphase ceramics. New processes developed to aid in fabrication of composites are given
- …