Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equations of motion”

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 145 records · Page 8

Soft de Sitter Effective Theory

Calculating the quantum evolution of a de Sitter universe on superhorizon scales is notoriously difficult. To address this challenge, we introduce the Soft de Sitter Effective Theory (SdSET). This framework holds for superhorizon modes whose comoving momentum is far below the UV scale, which is set by the inverse comoving horizon. The SdSET is formulated using the same approach that yields the Heavy Quark Effective Theory. The degrees of freedom that capture the long wavelength dynamics are identified with the growing and decaying solutions to the equations of motion. The operator expansion is organized using a power counting scheme, and loops can be regulated while respecting the low energy symmetries. For massive quantum fields in a fixed de Sitter background, power counting implies that all interactions beyond the horizon are irrelevant. Alternatively, if the fields are very light, the leading interactions are at most marginal, and resumming the associated logarithms using (dynamical) renormalization group techniques yields the evolution equation for canonical stochastic inflation. The SdSET is also applicable to models where gravity is dynamical, including inflation. In this case, diffeomorphism invariance ensures that all interactions are irrelevant, trivially implying the all-orders conservation of adiabatic density fluctuations and gravitational waves. We briefly touch on the application to slow-roll eternal inflation by identifying novel relevant operators. This work serves to demystify many aspects of perturbation theory outside the horizon, and has a variety of applications to problems of cosmological interest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Chapter Six - Derivation, cross-validation, and comparison of analytic formulas for electrostatic deflector aberrations

Here, we derived first and second order analytic aberration formulas in the horizontal transverse plane for an electrostatic deflector specified by the reference orbit radius, the central angle spanning the deflector, and its inhomogeneity coefficients. The derivation was performed using an iterative order-by-order perturbation method in a Frenet–Serret beamline coordinate system. We produced a C program edabrt for calculation of the first and second order aberrations using the formulas that we derived. The electrostatic deflector aberration formulas from Wollnik (1965) disagree with COSY INFINITY and the aberration formulas that we derived here, and they also deviate from the first and second order symplecticity conditions. The reason for this discrepancy is that the aberration formulas from Wollnik (1965) consider only motion in the main field of the deflector and do not properly account for fringe field effects, in particular the necessarily occurring change of kinetic energy due to the change in potential. The code GIOS uses electrostatic aberration formulas from Wollnik (1965) and exhibits the same discrepancy. We performed a comparison of first and second order electrostatic deflector aberrations for (1) the analytic formulas derived in this work, (2) differential-algebraic (DA) numerical integration of the equations of motion using the code COSY INFINITY, (3) the aberration formulas from Wollnik (1965), (4) aberrations computed using the code GIOS, and (5) the aberration formulas from Wollnik (1965) adjusted to account for the occurring change in potential. An electrostatic spherical deflector and an electrostatic cylindrical deflector were used as test cases for this comparison. There is excellent agreement between methods (1), (2), and (5), and these three methods satisfy the first and second order symplecticity conditions.

43 PARTICLE ACCELERATORS↗

Energy-Casimir, dynamically accessible, and Lagrangian stability of extended magnetohydrodynamic equilibria

The formal stability analysis of Eulerian extended magnetohydrodynamics (XMHD) equilibria is considered within the noncanonical Hamiltonian framework by means of the energy-Casimir variational principle and the dynamically accessible stability method. Specifically, we find explicit sufficient stability conditions for axisymmetric XMHD and Hall MHD (HMHD) equilibria with toroidal flow and for equilibria with arbitrary flow under constrained perturbations. The dynamically accessible, second-order variation of the Hamiltonian, which can potentially provide explicit stability criteria for generic equilibria, is also obtained. Moreover, we examine the Lagrangian stability of the general quasineutral two-fluid model written in terms of MHD-like variables, by finding the action and the Hamiltonian functionals of the linearized dynamics, working within a mixed Lagrangian-Eulerian framework. Upon neglecting electron mass, we derive a HMHD energy principle, and in addition, the perturbed induction equation arises from Hamilton’s equations of motion in view of a consistency condition for the relation between the perturbed magnetic potential and the canonical variables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Critical and near-critical relaxation of holographic superfluids

We investigate the relaxation of holographic superfluids after quenches, when the end state is either tuned to be exactly at the critical point, or very close to it. By solving the bulk equations of motion numerically, we demonstrate that in the former case the system exhibits a power law falloff, as well as an emergent discrete scale invariance. The latter case is in the regime dominated by critical slowing down, and we show that there is an intermediate time range before the onset of late-time exponential falloff, where the system behaves similarly to the critical point with its power law falloff. We further postulate a phenomenological Gross-Pitaevskii-like equation (corresponding to model F of Hohenberg and Halperin) that is able to make quantitative predictions for the behavior of the holographic superfluid after near-critical quenches into the superfluid and normal phase. Intriguingly, all parameters of our phenomenological equation, which describes the nonlinear time evolution, may be fixed with information from the static equilibrium solutions and linear response theory. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Real-time time-dependent density functional theory in Quantum Espresso and related codes

A set of subroutines that interact with an existing computer program, Quantum Espresso (QE), that takes a set of electronic orbitals from QE and propagates them forward in time in response to a number of perturbations. These perturbations include the motion of ions and a wide spectrum of applied electromagnetic fields (e.g., static magnetic fields, slowly-varying electric fields, x-rays, etc.). From the time evolution of the electronic orbitals comes time-varying values of physical observables like the electronic charge and current densities and atomic forces. These are further post-processed to yield properties of interest like stopping powers, conductivities, and the dynamic structure factor. The electronic orbitals are described in a plane wave basis and the equations of motion are carried out using a unitary or approximately unitary time propagation algorithm. Related work will include subroutines for interfacing this capability, as well as others that already exist within QE, with a fork of QE that is more extensible and modular.SAND2019-12611 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Baczewski, Andrew↗

Theoretical prediction of magnetic exchange coupling constants from broken-symmetry coupled cluster calculations

Exchange coupling constants ( J ) are fundamental to the understanding of spin spectra of magnetic systems. In this study, we investigate the broken-symmetry (BS) approaches of Noodleman and Yamaguchi in conjunction with coupled cluster (CC) methods to obtain exchange couplings. J values calculated from CC in this fashion converge smoothly toward the full configuration interaction result with increasing level of CC excitation. We compare this BS-CC scheme to the complementary equation-of-motion CC approach on a selection of bridged molecular cases and give results from a few other methodologies for context.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dispersion relation for gauge-free electromagnetic drift kinetics

Recently, a new approach to gyrokinetics, invariant under electromagnetic gauge transformations, was developed. The gyrocenter equations of motion are now expressed in terms of the perturbed fields instead of the potentials, in a form suitable for numerical simulations and analytic studies. In this paper, we verify that the long-wavelength limit, i.e., the drift-kinetic limit of the new gyrokinetic theory, is in line with existing work, providing a solid foundation for simulations. Here, we compute the dispersion relation of the new drift-kinetic theory in slab geometry and find agreement with a long-wavelength limit of the full Vlasov–Maxwell model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. Here, the non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH 3 demonstrate the accuracy of the non-hierarchical CDVR approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗

Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation

Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo

Abstract We present a discussion of recent progress in excited‐state‐specific quantum chemistry and quantum Monte Carlo alongside a demonstration of how a combination of methods from these two fields can offer reliably accurate excited state predictions across singly excited, doubly excited, and charge transfer states. Both of these fields have seen important advances supporting excited state simulation in recent years, including the introduction of more effective excited‐state‐specific optimization methods, improved handling of complicated wave function forms, and ways of explicitly balancing the quality of wave functions for ground and excited states. To emphasize the promise that exists at this intersection, we provide demonstrations using a combination of excited‐state‐specific complete active space self‐consistent field theory, selected configuration interaction, and state‐specific variance minimization. These demonstrations show that combining excited‐state‐specific quantum chemistry and variational Monte Carlo can be more reliably accurate than either equation of motion coupled cluster theory or multi‐reference perturbation theory, and that it can offer new clarity in cases where existing high‐level methods do not agree. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry

Otis, Leon↗

Non-standard axion electrodynamics and the dual Witten effect

Standard axion electrodynamics has two closely related features. First, the coupling of a massless axion field to photons is quantized, in units proportional to the electric gauge coupling squared. Second, the equations of motion tell us that a time-dependent axion field in a background magnetic field sources an effective electric current, but a time-dependent axion field in a background electric field has no effect. These properties, which manifestly violate electric-magnetic duality, play a crucial role in experimental searches for axions. Recently, electric-magnetic duality has been used to motivate the possible existence of non-standard axion couplings, which can both violate the usual quantization rule and exchange the roles of electric and magnetic fields in axion electrodynamics. We show that these non-standard couplings can be derived from SL(2,Z) duality, but that they come at a substantial cost: in non-standard axion electrodynamics, all electrically charged particles become dyons when the axion traverses its field range, in a dual form of the standard Witten effect monodromy. This implies that there are dyons near the weak scale, leads to a large axion mass induced by Standard Model fermion loops, and dramatically alters Higgs physics. We conclude that non-standard axion electrodynamics, although interesting to consider in abstract quantum field theory, is not phenomenologically viable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rounding out the story of higher derivative consistent truncations

At the two-derivative order, the group manifold reduction of heterotic supergravity on S 3 results in a half-maximal 7D gauged supergravity coupled to three vector multiplets, and a further truncation can be taken to remove the vector multiplets. We demonstrate that this truncation remains consistent at the four-derivative level; we do so both by analysis of the equations of motion and the supersymmetry variations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The no boundary density matrix

We discuss a no-boundary proposal for a subregion of the universe. In the classical approximation, this density matrix involves finding a specific classical solution of the equations of motion with no boundary. Beyond the usual no boundary condition at early times, we also have another no boundary condition in the region we trace out. We can find the prescription by starting from the usual Hartle-Hawking proposal for the wavefunction on a full slice and tracing out the unobserved region in the classical approximation. We discuss some specific subregions and compute the corresponding solutions. These geometries lead to phenomenologically unacceptable probabilities, as expected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Black hole perturbations of massive and partially massless spin-2 fields in (anti) de Sitter spacetime

We provide a systematic and comprehensive derivation of the linearized dynamics of massive and partially massless spin-2 particles in a Schwarzschild (anti) de Sitter black hole background, in four and higher spacetime dimensions. In particular, we show how to obtain the quadratic actions for the propagating modes and recast the resulting equations of motion in a Schrödinger-like form. In the case of partially massless fields in Schwarzschild de Sitter spacetime, we study the isospectrality between modes of different parity. In particular, we prove isospectrality analytically for modes with multipole number L = 1 in four spacetime dimensions, providing the explicit form of the underlying symmetry. We show that isospectrality between partially massless modes of different parity is broken in higher-dimensional Schwarzschild de Sitter spacetimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interaction of cosmological domain walls with large classical objects, like planets and satellites, and the flyby anomaly

Cosmological domain walls can be formed as a result of symmetry breaking at any epoch during the evolution of our universe. We study their interaction with a classical macroscopic object, like Earth or a satellite in Earth’s orbit. We set up an action that includes the interaction term between the massive classical object and the scalar field that the domain wall is made of. We use numerical calculations to solve the coupled equations of motion which describe the crossing between the domain wall and the classical object. Depending on the strength of the interaction, relative velocity and size, the object can be either stopped by the wall, or it can pass through it inducing deformations in the wall that cost energy. At the same time, the coupling to the scalar filed might change the object’s mass during the crossover. The fact that satellites in Earth’s orbit (or planets in Sun’s orbit) can change their mass and/or lose energy interacting with walls can be used as a new domain wall detection probe. For example, a typical velocity precision of a satellite is about 0.5 mm/s, which directly puts an upper limit on its mass change to ΔM/M ≲ 5 × 10 –17 . Alternatively, a known satellite flyby anomaly can easily be explained as an interaction with a closed domain wall. We also show that the presence of matter modifies the scalar filed potential and can locally create a bubble of the true vacuum, and thus trigger the decay of the false vacuum. For a critical bubble which is able to expand, such an interaction with the domain wall must be strong enough.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

New superstrata from three-dimensional supergravity

We find a two-parameter family of generalized superstrata that emerge as smooth, supersymmetric solutions in three-dimensional gauged supergravity coupled to additional scalar fields. This new family of generalized superstrata are smooth microstate geometries and may be thought of as supersymmetric Coulomb-branch extensions of the original superstrata in which the underlying supertube undergoes an elliptical deformation. These solutions had already been obtained numerically, and as series solutions, to the equations of motion, and some of them were conjectured to be supersymmetric. Here we prove the supersymmetry of an entire two-parameter family and we obtain a highly non-trivial analytic and smooth solution for a one-parameter limit in which the global symmetry of the metric is enhanced to SO(3). We also confirm that the other known families of microstrata are not supersymmetric. We conclude with a cursory analysis of some of the singular brane distributions that can be accessed from three-dimensional gauged supergravity while preserving the same supersymmetries as the superstratum, and therefore of the three-charge black hole.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗