Second-order constraints for equations of motion of constrained systems

1998 ◽  
Vol 3 (3) ◽  
pp. 240-248 ◽  
Author(s):  
Ye-Hwa Chen
2021 ◽  
Vol 502 (3) ◽  
pp. 3976-3992
Author(s):  
Mónica Hernández-Sánchez ◽  
Francisco-Shu Kitaura ◽  
Metin Ata ◽  
Claudio Dalla Vecchia

ABSTRACT We investigate higher order symplectic integration strategies within Bayesian cosmic density field reconstruction methods. In particular, we study the fourth-order discretization of Hamiltonian equations of motion (EoM). This is achieved by recursively applying the basic second-order leap-frog scheme (considering the single evaluation of the EoM) in a combination of even numbers of forward time integration steps with a single intermediate backward step. This largely reduces the number of evaluations and random gradient computations, as required in the usual second-order case for high-dimensional cases. We restrict this study to the lognormal-Poisson model, applied to a full volume halo catalogue in real space on a cubical mesh of 1250 h−1 Mpc side and 2563 cells. Hence, we neglect selection effects, redshift space distortions, and displacements. We note that those observational and cosmic evolution effects can be accounted for in subsequent Gibbs-sampling steps within the COSMIC BIRTH algorithm. We find that going from the usual second to fourth order in the leap-frog scheme shortens the burn-in phase by a factor of at least ∼30. This implies that 75–90 independent samples are obtained while the fastest second-order method converges. After convergence, the correlation lengths indicate an improvement factor of about 3.0 fewer gradient computations for meshes of 2563 cells. In the considered cosmological scenario, the traditional leap-frog scheme turns out to outperform higher order integration schemes only when considering lower dimensional problems, e.g. meshes with 643 cells. This gain in computational efficiency can help to go towards a full Bayesian analysis of the cosmological large-scale structure for upcoming galaxy surveys.


2005 ◽  
Vol 128 (3) ◽  
pp. 177-183 ◽  
Author(s):  
Sébastien Fouques ◽  
Harald E. Krogstad ◽  
Dag Myrhaug

Synthetic aperture radar (SAR) imaging of ocean waves involves both the geometry and the kinematics of the sea surface. However, the traditional linear wave theory fails to describe steep waves, which are likely to bring about specular reflection of the radar beam, and it may overestimate the surface fluid velocity that causes the so-called velocity bunching effect. Recently, the interest for a Lagrangian description of ocean gravity waves has increased. Such an approach considers the motion of individual labeled fluid particles and the free surface elevation is derived from the surface particles positions. The first order regular solution to the Lagrangian equations of motion for an inviscid and incompressible fluid is the so-called Gerstner wave. It shows realistic features such as sharper crests and broader troughs as the wave steepness increases. This paper proposes a second order irregular solution to these equations. The general features of the first and second order waves are described, and some statistical properties of various surface parameters such as the orbital velocity, slope, and mean curvature are studied.


Author(s):  
Simon A. Neild ◽  
Andrea Cammarano ◽  
David J. Wagg

In this paper we discuss a theoretical technique for decomposing multi-degree-of-freedom weakly nonlinear systems into a simpler form — an approach which has parallels with the well know method for linear modal analysis. The key outcome is that the system resonances, both linear and nonlinear are revealed by the transformation process. For each resonance, parameters can be obtained which characterise the backbone curves, and higher harmonic components of the response. The underlying mathematical technique is based on a near identity normal form transformation. This is an established technique for analysing weakly nonlinear vibrating systems, but in this approach we use a variation of the method for systems of equations written in second-order form. This is a much more natural approach for structural dynamics where the governing equations of motion are written in this form as standard practice. In fact the first step in the method is to carry out a linear modal transformation using linear modes as would typically done for a linear system. The near identity transform is then applied as a second step in the process and one which identifies the nonlinear resonances in the system being considered. For an example system with cubic nonlinearities, we show how the resulting transformed equations can be used to obtain a time independent representation of the system response. We will discuss how the analysis can be carried out with applied forcing, and how the approximations about response frequencies, made during the near-identity transformation, affect the accuracy of the technique. In fact we show that the second-order normal form approach can actually improve the predictions of sub- and super-harmonic responses. Finally we comment on how this theoretical technique could be used as part of a modal testing approach in future work.


2010 ◽  
Vol 25 (02) ◽  
pp. 135-141
Author(s):  
H. A. ELEGLA ◽  
N. I. FARAHAT

Motivated by the Hamilton–Jacobi approach of constrained systems, we analyze the classical structure of a four-dimensional superstring. The equations of motion for a singular system are obtained as total differential equations in many variables. The path integral quantization based on Hamilton–Jacobi approach is applied to quantize the system, and the integration is taken over the canonical phase space coordinates.


2018 ◽  
Vol 11 (1) ◽  
Author(s):  
Xi Kang ◽  
Xinsheng Zhang ◽  
Jian S. Dai

Reconfiguration identification of a mechanism is essential in design and analysis of reconfigurable mechanisms. However, reconfiguration identification of a multiloop reconfigurable mechanism is still a challenge. This paper establishes the first- and second-order kinematic model in the queer-square mechanism to obtain the constraint system by using the sequential operation of the Lie bracket in a bilinear form. Introducing a bilinear form to reduce the complexity of first- and second-order constraints, the constraint system with first- and second-order kinematics of the queer-square mechanism is attained in a simplified form. By obtaining the solutions of the constraint system, six motion branches of the queer-square mechanism are identified and their corresponding geometric conditions are presented. Moreover, the initial configuration space of the mechanism is obtained.


Entropy ◽  
2015 ◽  
Vol 17 (12) ◽  
pp. 6643-6662 ◽  
Author(s):  
Aimeric Colléaux ◽  
Sergio Zerbini

2005 ◽  
Vol 02 (03) ◽  
pp. 485-504 ◽  
Author(s):  
EQAB M. RABEI ◽  
EYAD H. HASSAN ◽  
HUMAM B. GHASSIB ◽  
S. MUSLIH

A general theory is given for quantizing both constrained and unconstrained systems with second-order Lagrangian, using the WKB approximation. In constrained systems, the constraints become conditions on the wave function to be satisfied in the semiclassical limit. This is illustrated with two examples.


Sign in / Sign up

Export Citation Format

Share Document