Lagrangian Formulation Research Articles

Schrödinger connections: from mathematical foundations towards Yano–Schrödinger cosmology

Abstract Schrödinger connections are a special class of affine connections, which despite being metric incompatible, preserve length of vectors under autoparallel transport. In the present paper, we introduce a novel coordinate-free formulation of Schrödinger connections. After recasting their basic properties in the language of differential geometry, we show that Schrödinger connections can be realized through torsion, non-metricity, or both. We then calculate the curvature tensors of Yano–Schrödinger geometry and present the first explicit example of a non-static Einstein manifold with torsion. We generalize the Raychaudhuri and Sachs equations to the Schrödinger geometry. The length-preserving property of these connections enables us to construct a Lagrangian formulation of the Sachs equation. We also obtain an equation for cosmological distances. After this geometric analysis, we build gravitational theories based on Yano–Schrödinger geometry, using both a metric and a metric-affine approach. For the latter, we introduce a novel cosmological hyperfluid that will source the Schrödinger geometry. Finally, we construct simple cosmological models within these theories and compare our results with observational data as well as the ΛCDM model.

Test Mass Capture Control for Drag-Free Satellite Based on State-Dependent Riccati Equation Method

The drag-free satellite plays an important role in the space-based gravitational wave observatory. The capture control of test mass after release is a crucial technology that can affect the success of the mission. The test mass must be released to the center of the electrostatic suspension cage accurately. This paper presents a nonlinear dynamic model of drag-free satellites in Lagrange formalism. A capture control scheme for test mass release phase is proposed based on the state-dependent Riccati equation (SDRE) strategy. To deal with the actuator saturation problem, a nonlinear saturation model is introduced to the dynamics of satellite, while the SDRE strategy is applied to the non-affine system. The effectiveness of the proposed methodology is verified by the numerical simulation for the drag-free satellite.

Theoretical and Numerical Modeling of Optical Switching of Epitaxial Nanostructures Based on Iron-Garnet Films

The paper presents a theoretical analysis of magnetization switching in a gadolinium ferrite garnet film due to the demagnetizing effect of a femtosecond laser pulse. Using the Lagrange formalism for a two-sublattice ferrimagnet, the effective Lagrangian, thermodynamic potential, and Rayleigh dissipative function are obtained. The phase diagram of the ferrite film is analyzed, and the main states of the system are identified. Magnetization switching diagrams and trajectories of the order parameter dynamics of the magnet are constructed. The ranges of magnetic fields, temperatures, and demagnetization values for the most efficient magnetization switching are analyzed.

Conformal Yang-Mills field in (A)dS space

Ordinary-derivative (second-derivative) Lagrangian formulation of classical conformal Yang-Mills field in the (A)dS space of six, eight, and ten dimensions is developed. For such conformal field, we develop two gauge invariant Lagrangian formulations which we refer to as generic formulation and decoupled formulation. In both formulations, the usual Yang-Mills field is accompanied by additional vector and scalar fields where the scalar fields are realized as Stueckelberg fields. In the generic formulation, the usual Yang-Mills field is realized as a primary field, while the additional vector fields are realized as auxiliary fields. In the decoupled formulation, the usual Yang-Mills field is realized as massless field, while the additional vector fields together with the Stueckelberg are realized as massive fields. Some massless/massive fields appear with the wrong sign of kinetic terms, hence demonstrating explicitly that the considered models are not unitary. The use of embedding space method allows us to treat the isometry symmetries of (A)dS space manifestly and obtain conformal transformations of fields in a relatively straightforward way. By accompanying each vector field by the respective gauge parameter, we introduce an extended gauge algebra. Levy-Maltsev decomposition of such algebra is noted. Use of the extended gauge algebra setup allows us to present concise form for the Lagrangian and gauge transformations of the conformal Yang-Mills field. Higher-derivative representation of the Lagrangian is also obtained.

Spur gear tooth root stress analysis by a 3D flexible multibody approach and a full-FE contact-based formulation

This paper proposes an original method to determine the gear tooth root stresses from a 3D finite element (FE) flexible multibody approach and a full-FE contact-based formulation. The contact problem is dealt with an augmented Lagrangian formulation whereas the analysis is performed by a preconditioned gradient solver (PCG). Tooth flank modifications are directly introduced within the 3D model. This one is thus able to take into account straightforwardly tooth bending and Hertzian-like deformations as well as the micro-geometry effect. Simulations are performed for several mesh periods, without making any assumptions about load distribution, tooth and gear blank flexibilities, and possible premature or delayed contacts between tooth pairs in quasi-static conditions. A precise distribution of tooth root stresses associated with instantaneous contacts conditions is then computed. For this study, a single stage spur gear with micro-geometry modifications corresponding to an arc-shaped profile crowning is modeled. Several output torques are considered. The obtained results are compared to those obtained using a 2D FE ISO-based model, where external forces are applied along the theoretical line of action.

A cortical field theory – dynamics and symmetries

We characterise cortical dynamics using partial differential equations (PDEs), analysing various connectivity patterns within the cortical sheet. This exploration yields diverse dynamics, encompassing wave equations and limit cycle activity. We presume balanced equations between excitatory and inhibitory neuronal units, reflecting the ubiquitous oscillatory patterns observed in electrophysiological measurements. Our derived dynamics comprise lowest-order wave equations (i.e., the Klein-Gordon model), limit cycle waves, higher-order PDE formulations, and transitions between limit cycles and near-zero states. Furthermore, we delve into the symmetries of the models using the Lagrangian formalism, distinguishing between continuous and discontinuous symmetries. These symmetries allow for mathematical expediency in the analysis of the model and could also be useful in studying the effect of symmetrical input from distributed cortical regions. Overall, our ability to derive multiple constraints on the fields — and predictions of the model — stems largely from the underlying assumption that the brain operates at a critical state. This assumption, in turn, drives the dynamics towards oscillatory or semi-conservative behaviour. Within this critical state, we can leverage results from the physics literature, which serve as analogues for neural fields, and implicit construct validity. Comparisons between our model predictions and electrophysiological findings from the literature — such as spectral power distribution across frequencies, wave propagation speed, epileptic seizure generation, and pattern formation over the cortical surface — demonstrate a close match. This study underscores the importance of utilizing symmetry preserving PDE formulations for further mechanistic insights into cortical activity.

A derivation of the Schrödinger equation from Feynman’s path-integral formulation of quantum mechanics

The equation of motion in the standard formulation of non-relativistic quantum mechanics, the Schrödinger equation, is based on the Hamiltonian. In contrast, in Feynman’s path-integral formulation of quantum mechanics, the equation of motion is the propagation equation, which is based on the Lagrangian. That these two different equations of motion are equivalent was shown by Feynman, who provided a derivation of the Schrödinger equation from the propagation equation. Surprisingly, however, while in classical mechanics there exists a simple relationship between the Hamiltonian and the Lagrangian, there is nothing in Feynman’s derivation that gives any indication of this relationship. Here I show that the equations of motion in the Hamiltonian and Lagrangian formulations of quantum mechanics are, in fact, simply related to each other.

Singularity formation of hydromagnetic waves in cold plasma

We study C1 blow-up of the compressible fluid model introduced by Gardner and Morikawa, which describes the dynamics of a magnetized cold plasma. We propose sufficient conditions that lead to C1 blow-up. In particular, we find that smooth solutions can break down in finite time even if the gradient of initial velocity is identically zero. The density and the gradient of the velocity become unbounded as time approaches the lifespan of the smooth solution. The Lagrangian formulation reduces the singularity formation problem to finding a zero of the associated second-order ODE.

Discrete time stability of augmented Lagrangian formalism based underactuated inverse dynamics control method

Stability problems of robotic systems arise sometimes suddenly, seemingly for no reason. The digital time sampling is often the main cause of these instabilities. Discrete models, which are capable of the prediction of stability, are available for low-degree-of-freedom template models of linear position and force control. However, for the inverse dynamics control of underactuated systems, the literature has a lack of generally applicable results related to the effect of time discretization. Several control approaches are available in the literature out of which a widely used one, the augmented Lagrangian formalism and its stability properties are analysed in this work. Theoretical stability properties are obtained for a generally usable, linear, underactuated, two-degree-of-freedom constrained template model. The actuator dynamics, the finite difference approximation of the feedback velocity and the filtering of the feedback data are considered in the model. These phenomena strongly affects the stability properties. The theoretically obtained stability maps are experimentally validated on an underactuated crane-like indoor robot. The position and orientation accuracy of the robot were assessed: the absolute position error was below 30 mm and the orientation error was below 3°.

The Dynamics Stability Judging Method for a Gantry-type Welding Robot

In order to improve the dynamic performance of gantry-type welding robot, this paper introduces the structural dynamics model and Lyapunov stability criterion of the gantry-type welding robot frame with nine degrees of freedom. The dynamic equations of the welding robot satisfy the Euler–Lagrange form, and the calculation index of welding robot Lyapunov stability is judged by matrix eigenvalues. Combining with the specific design parameters of an industrial welding robot, the motion stability of the gantry-type welding robot frame and the displacement response of the welding torch are analyzed. The influence of the change of the frame design structure on the stability of the welding torch for end effector of welding robot is discussed, furthermore, the improvement direction of vibration damping and anti-interference of welding robot system is pointed out. Finally, the dynamics simulation of the optimized gantry-type welding robot end effector is carried out by using ADAMS software, and the torque variation curve of the end effector is obtained, which provides a strong basis for the selection of motor models and the design of control system.

Experimental investigation of denoising electrocardiogram using lagrange form of hermite interpolating polynomial with chebyshev nodes

Experimental investigation of denoising electrocardiogram using lagrange form of hermite interpolating polynomial with chebyshev nodes

Exploring arteriolar atherosclerosis: laminar blood flow across stenosis with fluid-structure interaction and gravitational effects

Abstract In response to the unanswered relevant questions surrounding atherosclerosis, it becomes imperative to investigate arterioles using sophisticated mathematical modelling techniques to shed light on critical stress and strain patterns influenced by gravity. The primary objective of this study is to scrutinize flow characteristics and probe stress and strain distributions experienced by the intima layer of arterioles, encompassing coronary, renal, cerebral, mesenteric, and pulmonary arteries, under gravitational forces. This investigation employs a fluid-structure interaction methodology utilizing arbitrary Eulerian–Lagrangian formulation. The study delves into blood flow characteristics within coronary, renal, cerebral, mesenteric, and pulmonary arterioles using the fluid-structure interaction technique, employing an arbitrary Eulerian–Lagrangian formulation. It thoroughly examines various biomechanical parameters such as the Cauchy–Green stress tensor, Principal strain, Piola–Kirchoff stress tensor, deformation tensor, and volume strain along the intima layer under the gravitational influence, elucidating vulnerable regions prone to endothelial dysfunction. Higher values of δV are found at the left shoulder and in the intima’s post stenosis area due to the pressure gradient along the flow channel, whereas other intima regions show a null volume strain. A thorough understanding of stress distribution is essential to create focused therapies to lessen vascular health problems. The stress in the post-stenosis region seems to affect the endothelial layer to a significant extent.

Nonlinear analysis of spatial trusses with different strain measures and compressible solid

This paper investigates the nonlinear behavior of spatial truss elements under finite deformations, focusing on the impact of various strain measures in compressible materials. We examine both Total Lagrangian (using engineering and Green–Lagrange strains) and Eulerian formulations (using natural, Biot, and Almansi strains). The analysis assumes a linear spatial hyperelastic material where Cauchy stress is proportional to axial natural strain via Young’s modulus. For infinitesimal strains, Young’s modulus remains consistent across different stress/strain pairs. In the finite strain regime, we derive a nonlinear secant modulus based on Young’s modulus. Internal force vectors and tangent stiffness matrices are computed using the direction cosines of the truss element in its deformed state. The paper demonstrates that for infinitesimal deformations, adjusting the modulus of elasticity when using different stress/strain pairs is unnecessary. However, for finite deformations, it is essential to adjust the modulus of elasticity. Numerical simulations validate the performance of the proposed 3D truss element against established formulations. This research offers critical insights into the nonlinear response of spatial trusses, guiding the selection of appropriate strain measures for enhanced accuracy in engineering applications. These findings contribute to more reliable and efficient structural designs, especially in scenarios involving finite deformations and compressible materials.

A modular finite element approach to saturated poroelasticity dynamics: Fluid–solid coupling with Neo-Hookean material and incompressible flow

Several methods have been developed to model the dynamic behavior of saturated porous media. However, most of them are suitable only for small strain and small displacement problems and are built in a monolithic way, so that individual improvements in the solution of the solid or fluid phases can be difficult. This study shows a macroscopic approach through a partitioned fluid–solid coupling, in which the skeleton solid is considered to behave as a Neo-Hookean material and the interstitial flow is incompressible following the Stokes–Brinkman model. The porous solid is numerically modeled with a total Lagrangian position-based finite element formulation, while an Arbitrary Lagrangian-Eulerian stabilized finite element approach is employed for the porous medium flow dynamics. In both fields, an averaging procedure is applied to homogenize the problem, resulting in a macroscopic continuous phase. The solid and fluid homogenized domains are overlapped and strongly coupled, based on a block-iterative solution scheme. Two-dimensional simulations of wave propagation in saturated porous media are employed to validate the proposed formulation through a comprehensive comparison with analytical and numerical results from the literature. The analyses underscore the proposed formulation as a robust and precise modular approach for addressing dynamic problems in poroelasticity.

Interlaminar fracture toughness of metal/composite bonded joints under high-speed mode I loading considering the elastic vibration

In the present work, the mode I fracture toughness at the initiation of a crack in a metal-composite bonded joint under high-speed loading is calculated, considering the effect of the elastic vibration of the specimen. A new hybrid analytical-experimental method is developed, which integrates kinematics from Digital Image Correlation (DIC) analysis to an analytical model to determine the fracture toughness and its individual components. The analytical model is based on Timoshenko beam theory and incorporates a Lagrangian contact formulation. The conditions under which the hybrid method can be applied are analysed. The experiments are conducted on a novel test set-up, featuring a wedge accelerated by a projectile from an air gun that is able to load the Double Cantilever Beam (DCB) specimen with an opening displacement rate up to 2m/s. The results reveal a correlation between the fracture toughness, the displacement rate, and the vibrational characteristics of the specimen.

Asymptotic symmetries from the string worldsheet

In IIB string theory on AdS3 background with NS-NS fluxes, we show that Brown-Henneaux asymptotic Killing vectors can be derived by requiring both the worldsheet equations of motion and Virasoro constraints are preserved near the asymptotic boundary of the target spacetime. The charges on the worldsheet that generate the corresponding transformations can be written down in both the Lagrangian formalism and Hamiltonian formalism. This provides a method of studying asymptotic symmetry of the target spacetime directly from worldsheet string theories, without using results from the supergravity limit. As an example, we apply this method to flat spacetime in three dimensions and obtain the BMS3 generators on the worldsheet theory.

Existence of non-singular stellar solutions within the context of electromagnetic field: a comparison between minimal and non-minimal gravity models

In this paper, we explore the existence of various non-singular compact stellar solutions influenced by the Maxwell field within the matter-geometry coupling based modified gravity. We start this analysis by considering a static spherically symmetric spacetime which is associated with the isotropic matter distribution. We then determine the field equations corresponding to two specific functions of this modified theory. Along with these models, we also adopt different forms of the matter Lagrangian. We observe several unknowns in these equations such as the metric potentials, charge and fluid parameters. Thus, the embedding class-one condition and a particular realistic equation of state is used to construct their corresponding solutions. The former condition provides the metric components possessing three constants, and we calculate them through junction conditions. Further, four developed models are graphically analyzed under different parametric values. Finally, we find all our developed solutions well-agreeing with the physical requirements, offering valuable insights for future explorations of the stellar compositions in this theory.

No time to derive: unraveling total time derivatives in in-in perturbation theory

The in-in formalism provides a way to systematically organize the calculation of primordial correlation functions. Although its theoretical foundations are now firmly settled, the treatment of total time derivative interactions, incorrectly trivialized as “boundary terms”, has been the subject of intense discussions and conceptual mistakes. In this work, we demystify the use of total time derivatives — as well as terms proportional to the linear equations of motion — and show that they can lead to artificially large contributions cancelling at different orders of the in-in operator formalism. We discuss the treatment of total time derivative interactions in the Lagrangian path integral formulation of the in-in perturbation theory, and we showcase the importance of interaction terms proportional to linear equations of motion. We then provide a new route to the calculation of primordial correlation functions, which avoids the generation of total time derivatives, by working directly at the level of the full Hamiltonian in terms of phase-space variables. Instead of integrating by parts, we perform canonical transformations to simplify interactions. We explain how to retrieve correlation functions of the initial phase-space variables from the knowledge of the ones after canonical transformations. As an important first application, we find the explicit sizes of Hamiltonian cubic interactions in single-field inflation with canonical kinetic terms and for any background evolution, straight in terms of the primordial curvature perturbation and its canonical conjugate momentum, as well as the corresponding ones in the tensor sector, and the ones mixing scalars and tensors. We also briefly comment on quartic interactions. Our results are important for performing complete calculations of exchange diagrams in inflation, such as the (scalar and tensor) exchange trispectrum and the one-loop power spectrum. Being already written in a form amenable to characterize quantum properties of primordial fluctuations, they also promise to shed light on the non-linear dynamics of quantum states during inflation.

Lagrangian formulation for perfect fluid equations with the ℓ -conformal Galilei symmetry

Lagrangian formulation for perfect fluid equations which hold invariant under the ℓ-conformal Galilei group with half-integer ℓ is proposed. It is based on a Clebsch-type parametrization and reproduces Lagrangian description of the Euler fluid equations for ℓ=12. The transition from the Lagrangian formulation to the Hamiltonian one is analyzed in detail. Published by the American Physical Society 2024

A Discrete Hamilton–Jacobi Theory for Contact Hamiltonian Dynamics

In this paper, we propose a discrete Hamilton–Jacobi theory for (discrete) Hamiltonian dynamics defined on a (discrete) contact manifold. To this end, we first provide a novel geometric Hamilton–Jacobi theory for continuous contact Hamiltonian dynamics. Then, rooting on the discrete contact Lagrangian formulation, we obtain the discrete equations for Hamiltonian dynamics by the discrete Legendre transformation. Based on the discrete contact Hamilton equation, we construct a discrete Hamilton–Jacobi equation for contact Hamiltonian dynamics. We show how the discrete Hamilton–Jacobi equation is related to the continuous Hamilton–Jacobi theory presented in this work. Then, we propose geometric foundations of the discrete Hamilton–Jacobi equations on contact manifolds in terms of discrete contact flows. At the end of the paper, we provide a numerical example to test the theory.