Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian mechanics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Electromagnetic and gravitational interactions from Lagrangian mechanics

Highlights: • Background Maxwell and Einstein fields are universal features of Lagrangian systems. • The Lagrangian function is promoted to a dynamic field. • The Einstein–Maxwell theory is derived from Lagrangian mechanics. Background fields of electromagnetic and gravitational type emerge in the low kinetic energy limit of any regular Lagrangian system and, in particular, in the corresponding limit of any spacetime theory in which the free motion of test particles is described by an unspecified regular Lagrangian. Electromagnetic and gravitational type interactions are therefore a universal feature of low kinetic energy Lagrangian systems. These background interactions can be consistently turned into dynamic Einstein–Maxwell fields by promoting the Lagrangian function to a dynamic scalar field on the tangent bundle of the configuration space. Accordingly, Einstein–Maxwell theory can be deduced from the assumption that the motion of elementary test particles in spacetime is described by Lagrangian mechanics. For higher kinetic energy-type values, identified with the square of the invariant mass of the particle, the Lagrangian induces higher rank interactions that seem however too weak to have been detected in spacetime physics, but which might prove relevant at the Planck scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

One-Dimensional Multi-Velocity Capabilities for Arbitrary Lagrangian-Eulerian Normal Contact Mechanics

Lagrangian and Arbitrary Lagrangian Eulerian (ALE) hydrodynamics codes such as FLAG form the backbone of many mission-critical multi physics simulations at Los Alamos National Laboratory. Critical to pre forming high fidelity simulations with these codes are Lagrangian and ALE contact algorithms, which allow materials to collide, slide, and sep arate throughout a simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Ergodic Lagrangian dynamics in a superhero universe

We present a fictional scenario that, while undeniably whimsical, provides the foundation for a unique exercise in extended problem solving, physics analysis, and quantitative model development. Starting with the foundational premise of the Wild Cards shared-world superhero universe, we demonstrate how a variety of concepts appropriate to the advanced undergraduate level—ergodicity, functional analysis, Lagrangian mechanics, and the ever-important simplifying approximation—can be combined into a rich, coherent mathematical model. The goal of this case study is to develop a useful pedagogical exercise in exploring an open-ended research question that presents, at first glance, no clear path forward. Being both eclectic and lengthy, this exercise offers a unique way for students to apply their core physics and mathematics education. It is perhaps best used within a senior honors seminar or within a brief (e.g., January term) elective class.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Light-Dependent Impedance Spectra and Transient Photoconductivity in a Ruddlesden–Popper 2D Lead–Halide Perovskite Revealed by Electrical Scanned Probe Microscopy and Accompanying Theory

Electric force microscopy was used to record the light-dependent impedance spectrum and the probe transient photoconductivity of a film of butylammonium lead iodide, BA 2 PbI 4 , a 2D Ruddlesden–Popper perovskite semiconductor. The impedance spectrum of BA 2 PbI 4 showed modest changes as the illumination intensity was varied up to 1400 mW cm –2 , in contrast with the comparatively dramatic changes seen for 3D lead–halide perovskites under similar conditions. BA 2 PbI 4 ’s light-induced conductivity had a rise time and decay time of ~100 μs, 10 4 slower than expected from direct electron–hole recombination yet 10 5 faster than the conductivity-recovery times recently observed in 3D lead–halide perovskites and attributed to the relaxation of photogenerated vacancies. What sample properties are probed by electric force microscope measurements remains an open question. A Lagrangian-mechanics treatment of the electric force microscope experiment was recently introduced by Dwyer, Harrell, and Marohn which enabled the calculation of steady-state electric force microscope signals in terms of a complex sample impedance. Here this impedance treatment of the tip–sample interaction is extended, through the introduction of a time-dependent transfer function, to include time-resolved electrical scanned probe measurements. Here, it is shown that the signal in a phase-kick electric force microscope experiment, and therefore also the signal in a time-resolved electrostatic force microscope experiment, can be written explicitly in terms of the sample’s time-dependent resistance (i.e., conductivity).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Phase transitions in Zr at sub-nanosecond time scales

Here solid-solid phase transitions are investigated in Zr thin films using shock compression induced by a short laser pulse (<1 ns). Shock wave profiles are measured at free surfaces for films of different thicknesses (a few micrometers) using chirped-pulse line velocimetry with 10 ps time resolution. Experiments are performed at pressures up to ~50 GPa, which is sufficient to reach the ω and β phases under equilibrium conditions. The shock wave structures are analyzed using a general Lagrangian analysis method, which allows for the calculation of stress–strain paths and assessments of phase transition behavior. In agreement with recent short laser pulse experiments using ultra-fast x-ray diffraction, we do not find any clear evidence of the α–ω transition, though this would be expected from the phase diagram. Instead, we infer a direct transformation to a metastable β-phase at lower shock pressures (<20 GPa) and equilibrium β at higher pressures. Through the velocimetry analysis, we find α–β transformation onset times of less than ~100 ps and completion times of less than ~200 ps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A Vibrational Energy Harvesting Sensor Based on Linear and Rotational Electromechanical Effects

In this investigation, a magnetically coupled double-spring design is presented for harvesting low-level non-stationary random vibrational energy. The sensor relies on multimodal coupling between the translation and rotation of a two-spring magnet and coil system to widen the harvesting bandwidth. Energy methods are used to develop a model to characterize the electromechanical response of the system, the solution of which is obtained using stochastic techniques based on a particle swarm algorithm. This approach provides an efficient method to estimate system parameters that otherwise are difficult or impossible to determine with independent measurements. The experimental results demonstrate agreement with the theoretical predictions over a limited bandwidth. The sensor can effectively harvest non-stationary vibration energy down to 10 -4 g within a limited bandwidth of 130–150 Hz. The sensor prototype has an operational volume of 2.6 cm 3 with a calculated power density of 0.2 W/cm 3 . The sensor’s small size results in a coupling efficiency of approximately 6% across the tested bandwidth.

42 ENGINEERING↗
DOE OSTI2025

What can solve the strong CP problem?

Three possible strategies have been advocated to solve the strong CP problem. The first is the axion, a dynamical mechanism that relaxes any initial value of the CP violating angle $\overline{θ}$ to zero. The second is the imposition of new symmetries that are believed to set $\overline{θ}$ to zero in the UV. The third is the acceptance of the fine tuning of parameters. We argue that the latter two solutions do not solve the strong CP problem. The θ term of QCD is not a parameter — it does not exist in the Hamiltonian. Rather, it is a property of the quantum state that our universe finds itself in, arising from the fact that there are CP violating states of a CP preserving Hamiltonian. It is not eliminated by imposing parity as a symmetry since the underlying theory is already parity symmetric and that does not preclude the existence of CP violating states. Moreover, since the value of θ realized in our universe is a consequence of measurement, it is inherently random and cannot be fine tuned by choice of parameters. Rather any fine tuning would require a tuning between parameters in the theory and the random outcome of measurement. Our results considerably strengthen the case for the existence of the axion and axion dark matter. The confusion around θ arises from the fact that unlike classical mechanics, the Hamiltonian and Lagrangian are not equivalent in quantum mechanics. The Hamiltonian defines the differential time evolution, whereas the Lagrangian is a solution to this evolution. Consequently, initial conditions could in principle appear in the Lagrangian but not in the Hamiltonian. This results in aspects of the initial condition such as θ misleadingly appearing in the Lagrangian as parameters. We comment on the similarity between the θ vacua and the violations of the constraint equations of classical gauge theories in quantum mechanics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
DOE OSTI2020

An augmented Lagrangian filter method

Here, we introduce a filter mechanism to enforce convergence for augmented Lagrangian methods for nonlinear programming. In contrast to traditional augmented Lagrangian methods, our approach does not require the use of forcing sequences that drive the first-order error to zero. Instead, we employ a filter to drive the optimality measures to zero. Our algorithm is flexible in the sense that it allows for equality-constrained quadratic programming steps to accelerate local convergence. We also include a feasibility restoration phase that allows fast detection of infeasible problems. We provide a convergence proof that shows that our algorithm converges to first-order stationary points. We provide preliminary numerical results that demonstrate the effectiveness of our proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Insights on sources and formation mechanisms of liquid-bearing clouds over MOSAiC examined from a Lagrangian framework

Understanding Arctic stratiform liquid-bearing cloud life cycles and properly representing these life cycles in models is crucial for evaluations of cloud feedbacks as well as the faithfulness of climate projections for this rapidly warming region. Examination of cloud life cycles typically requires analyses of cloud evolution and origins on short time scales, on the order of hours to several days. Measurements from the recent Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition provide a unique view of the current state of the central Arctic over an annual cycle. Here, we use the MOSAiC radiosonde measurements to detect liquid-bearing cloud layers over full atmospheric columns and to examine the cloud-generating air masses’ properties. We perform 5-day (120 h) back-trajectory calculations for every detected cloud and cluster them using a unique set of variables extracted from these trajectories informed by ERA5 reanalysis data. This clustering method enables us to separate between the air mass source regions such as ice-covered Arctic and midlatitude open water. We find that moisture intrusions into the central Arctic typically result in multilayer liquid-bearing cloud structures and that more than half of multilayer profiles include overlying liquid-bearing clouds originating in different types of air masses. Finally, we conclude that Arctic cloud formation via prolonged radiative cooling of elevated stable subsaturated air masses circulating over the Arctic can occur frequently (up to 20% of detected clouds in the sounding data set) and may lead to a significant impact of ensuing clouds on the surface energy budget, including net surface warming in some cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Magnetic reconnection and thermal equilibration

When a magnetic field is forced to evolve on a time scale sev, as by footpoint motions driving the solar corona or non-axisymmetric instabilities in tokamaks, the magnetic field lines undergo large-scale changes in topology on a time scale approximately an order of magnitude longer than $\tau$ ev . But the physics that allows such changes operates on a time scale eight or more orders of magnitude slower. An analogous phenomenon occurs in air. Temperature equilibration occurs on a time scale approximately an order of magnitude longer than it takes air to cross a room, $\tau$ ev , although the physical mechanism that allows temperature equilibration is approximately four orders of magnitude slower than $\tau$ ev . The use of Lagrangian coordinates allows the fundamental equations to be solved and both phenomena explained. The theories of thermal equilibration and magnetic reconnection are developed in parallel to help readers obtain an understanding of the importance and implications of analyses using Lagrangian coordinates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method

Reverse Osmosis (RO) is a promising technology to address the impending water-crisis in the upcoming decades. RO at high salinities and high pressures is challenging due to membrane compaction that changes its porosity and permeability. In this study, we present a simulation methodology for membrane structural mechanics to understand pore size distribution and permeability variations under high pressure. We use the material point method (MPM), that solves the solid mechanics equations in a Lagrangian framework. MPM provides many features that make it well-suited for simulating mesoporous membranes. The Lagrangian framework allows for large deformations, easy integration of constitutive models and direct import of complex geometries as particles. The spatial discretization in our MPM solver is achieved using linear or cubic-spline shape functions while the time integration is carried out using the explicit Euler method. A series of images containing detailed pore structures obtained from X-ray tomography is first converted to a collection of material points. Compressive loads are applied to the top layer of the membrane to simulate the application of pressure. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compact in detail.

compaction↗
DOE OSTI2021

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The Dynamics of Megafire Smoke Plumes in Climate Models: Why a Converged Solution Matters for Physical Interpretations

As the climate system warms, megafires have become more frequent with devastating effects. A byproduct of these events is the creation of smoke plumes that can rise into the stratosphere and spread across the globe where they reside for many months. To gain a deeper understanding of the plume dynamics, global climate simulations of a megafire were performed at a wide range of grid spacings from 2.0° down to 7 km, including a 7 km nonhydrostatic experiment. The analysis focuses on how the resolved dynamics affects the specification of the plume characteristics such as injection height and black carbon (BC) mass. Prior studies initialize the smoke plume at one or a few grid points and this is shown here to produce severely dissipative dynamics. In order to validate such simulations with observations, enhancements of the plume characteristics to offset the dissipation is necessary. Using a numerically converged simulation, sensitivity tests show that to approximate the observed stratospheric lifetime, a reduction in BC fraction by 50% is necessary for external mixtures. The vorticity dynamics of the plume is also analyzed with a Lagrangian budget to understand the mechanisms responsible for the evolution of a collocated anticyclonic vortex. The results can be distilled down into a simple conceptual model. As the plume rises, the air diverges at the top of the updraft where the largest concentrations of smoke are found. This divergence induces a dilution of the background cyclonic absolute vorticity producing an anticyclonic vortex. Vortex decay occurs from opposite arguments.

54 ENVIRONMENTAL SCIENCES↗