Theory and Applications of Kalman Filtering
Kalman filtering applications in space navigation, underwater operations, surface navigation, and Gauss-Markov processes
C. T. Leondes↗
SEARCH · Engineering Papers
Results for “operator theory”
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.
Fast Particle-Wave Interactions and Alfvén Eigenmodes in JET Tokamak Plasmas
This document serves as the closeout report of DOE Grant Award No. DE-FG02-99ER54563 with project period 1 April 2015 through 31 March 2020. The project comprised the international collaboration between MIT and EU scientists on the JET facility to improve our understanding of the physics of energetic particle-wave interactions by measuring the damping rates of stable Alfvén Eigenmodes (AEs) and unstable energetic particle driven modes. More specifically, this project involved the continued participation of MIT, the Culham Centre for Fusion Energy, the Swiss Plasma Center, and theorists from various European laboratories and from UC Irvine. Past contributions from the University of São Paulo are gratefully acknowledged. The Alfvén Eigenmode Active Diagnostic (AEAD) on JET discharges was successfully upgraded in this period and has extracted physically useful information, in primarily deuterium plasma discharges. It is expected that such experiments would be continued during the DT campaign in CY 2021 to assess the damping rates of similar modes in the presence of alpha particles, a product of fusing burning plasma. Only JET would carry out such experiments in the near term in the world. It is important to note that this collaboration was continued under DOE Grant Award No. DE-SC0014264 from 1 April 2020 through the present. In this grant period, the AEAD was upgraded with a set of individual amplifiers for each of a set of six antennas in two toroidally opposite locations. These new amplifiers allowed targeted selection of the antennas’ toroidal spectrum for AEs of interest with toroidal mode number | n | ≤ 20. The resonant detection and measurement of the damping rates of AEs was obtained from magnetic probes that could be compared with theory and simulations. This is a key topic of investigation for ITER and all other next-step fusion experiments, where such modes will interact with energetic particles produced by fusion reactions (alpha particles), Neutral Beam Injection (NBI) and Ion Cyclotron Resonance Heating (ICRH). The majority of the observations during this period were of Toroidal Alfvén Eigenmodes (TAEs) as these are most commonly observed on JET; a number of measurements were made during dedicated TAE experiments. However, a new set of lower frequency band filters were procured with the goal of studying Geodesic Acoustic Modes (GAMs), Beta (Acoustic) AEs (BAE/BAAEs), and Reverse Shear AEs (RSAEs), also predicted by theory. Commissioning and optimizations were successfully completed in this period, and the diagnostic has been in operation during the more recent JET campaigns. Initial results have been obtained in dedicated TAE experiments and successfully compared to drift-kinetic theory. A wide range of theoretical studies have been undertaken in support of the upcoming JET campaigns experiments. To supplement the ongoing use of ideal MHD codes, such as MISHKA during the studies of TAEs, gyrokinetic simulations of low frequency AEs were performed in collaboration with UC Irvine. The Gyrokinetic Toroidal Code (GTC) was used to determine the structure, frequencies and stability of AEs in JET plasmas. Thus, a solid scientific foundation has been laid for future DT campaigns in JET in the 2021 operation period when the damping rates of relevant modes in the presence of alpha particles could be assessed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Fermionic mean-field theory as a tool for studying spin Hamiltonians
The Jordan–Wigner transformation permits one to convert spin 1/2 operators into spinless fermion ones, or vice versa. In some cases, it transforms an interacting spin Hamiltonian into a noninteracting fermionic one, which is exactly solved at the mean-field level. Even when the resulting fermionic Hamiltonian is interacting, its mean-field solution can provide surprisingly accurate energies and correlation functions. Furthermore, Jordan–Wigner is, however, only one possible means of interconverting spin and fermionic degrees of freedom. Here, we apply several such techniques to the XXZ and J 1 –J 2 Heisenberg models, as well as to the pairing or reduced Bardeen–Cooper–Schrieffer Hamiltonian, with the aim of discovering which of these mappings is most useful in applying fermionic mean-field theory to the study of spin Hamiltonians.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Data-driven learning of nonlocal physics from high-fidelity synthetic data
A key challenge to nonlocal models is the analytical complexity of deriving them from first principles, and frequently their use is justified a posteriori. Here, we extract nonlocal models from data, circumventing these challenges and providing data-driven justification for the resulting model form. Extracting data-driven surrogates is a major challenge for machine learning (ML) approaches, due to nonlinearities and lack of convexity — it is particularly challenging to extract surrogates which are provably well-posed and numerically stable. Our scheme not only yields a convex optimization problem, but also allows extraction of nonlocal models whose kernels may be partially negative while maintaining well-posedness even in small-data regimes. To achieve this, based on established nonlocal theory, we embed in our algorithm sufficient conditions on the non-positive part of the kernel that guarantee well-posedness of the learnt operator. These conditions are imposed as inequality constraints to meet the requisite conditions of the nonlocal theory. We demonstrate this workflow for a range of applications, including reproduction of manufactured nonlocal kernels; numerical homogenization of Darcy flow associated with a heterogeneous periodic microstructure; nonlocal approximation to high-order local transport phenomena; and approximation of globally supported fractional diffusion operators by truncated kernels.
42 ENGINEERING↗
The Physics Imposed on a Streaming Operator by Spherical Transport Problems
The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.
integral operator↗
Quantum computing for energy correlators
In recent years, energy correlators have emerged as powerful observables for probing the fragmentation dynamics of high-energy collisions. We introduce the first numerical strategy for calculating energy correlators using the Hamiltonian lattice approach, providing access to the intriguing nonperturbative dynamics of these observables. Furthermore, motivated by rapid advances in quantum computing hardware and algorithms, we propose a quantum algorithm for calculating energy correlators in quantum field theories. This algorithm includes ground state preparation, the application of source, sink, energy flux and real-time evolution operators, and the Hadamard test. We validate our approach by applying it to the SU(2) pure gauge theory in 2 + 1 dimensions on 3 × 3 and 5 × 5 honeycomb lattices with 𝑗 max = $\frac{1}{2}$ at various couplings, utilizing both classical methods and the quantum algorithm, the latter tested using the IBM emulator for specific configurations. The results are consistent with the expected behavior of the strong coupling regime and motivate a more comprehensive study to probe the confinement dynamics across the weak and strong coupling regimes.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Manifestly causal in-in perturbation theory about the interacting vacuum
In-In perturbation theory is a vital tool for cosmology and nonequilibrium physics. Here, we reconcile an apparent conflict between two of its important aspects with particular relevance to De Sitter/inflationary contexts: (i) the need to slightly deform unitary time evolution with an iϵ prescription that projects the free (“Bunch-Davies”) vacuum onto the interacting vacuum and renders vertex integrals well-defined, and (ii) Weinberg’s “nested commutator” reformulation of in-in perturbation theory which makes manifest the constraints of causality within expectation values of local operators, assuming exact unitarity. We show that a modified iϵ prescription maintains the exact unitarity on which the derivation of (ii) rests, while nontrivially agreeing with (i) to all orders of perturbation theory.
79 ASTRONOMY AND ASTROPHYSICS↗
Stability of a rigid rotor supported on oil-film journal bearings under dynamic load
Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film thickness from the predicted journal trajectories. These do not give any information about the subsynchronous whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, necessary to know whether the bearing system operation is stable or not under such an operating condition. The stability characteristics of the system are analyzed. A linearized perturbation theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient analysis of a rigid rotor supported on oil journal bearings under: (1) a unidirectional constant load, (2) a unidirectional periodic load, and (3) variable rotating load are performed. The hydrodynamic forces are calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, the journal trajectories are obtained for several different operating conditions.
Majumdar, B. C.↗
Role of Chiral Two-Body Currents in 6 Li Magnetic Properties in Light of a New Precision Measurement with the Relative Self-Absorption Technique
A direct measurement of the decay width of the excited $0^+_1$ state of 6 Li using the relative self-absorption technique is reported. Our value of $Γ_{γ,0^{+}_{1}} → 1^{+}_{1}$ = 8.17 (14) stat. (11) syst. eV provides sufficiently low experimental uncertainties to test modern theories of nuclear forces. The corresponding transition rate is compared to the results of ab initio calculations based on chiral effective field theory that take into account contributions to the magnetic dipole operator beyond leading order. This allows for a precision test of the impact of two-body currents that enter at next-to-leading order.
36 MATERIALS SCIENCE↗
Renormalization of the Standard Model Effective Field Theory from geometry
S-matrix elements are invariant under field redefinitions of the Lagrangian. They are determined by geometric quantities such as the curvature of the field-space manifold of scalar and gauge fields. We present a formalism where scalar and gauge fields are treated together, with a metric on the combined space of both types of fields. Scalar and gauge scattering amplitudes are given by the Riemann curvature Rijkl of this combined space, with indices i, j, k, l chosen to be scalar or gauge indices depending on the type of external particle. One-loop divergences can also be computed in terms of geometric invariants of the combined space, which greatly simplifies the computation of renormalization group equations. We apply our formalism to the Standard Model Effective Field Theory (SMEFT), and compute the renormalization group equations for even-parity bosonic operators to mass dimension eight.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Causality, unitarity, and the weak gravity conjecture
We consider the shift of charge-to-mass ratio for extremal black holes in the context of effective field theory, motivated by the Weak Gravity Conjecture. We constrain extremality corrections in different regimes subject to unitarity and causality constraints. In the asymptotic IR, we demonstrate that for any supersymmetric theory in flat space, and for all minimally coupled theories, logarithmic running at one loop pushes the Wilson coefficient of certain four-derivative operators to be larger at lower energies, guaranteeing the existence of sufficiently large black holes with Q > M. We identify two exceptional cases of nonsupersymmetric theories involving large numbers of light states and Planck-scale nonminimal couplings, in which the sign of the running is reversed, leading to black holes with negative corrections to Q/M in the deep IR, but argue that these do not rule out extremal black holes as the requisite charged states for the WGC. We separately show that causality and unitarity imply that the leading threshold corrections to the effective action from integrating out massive states, in any weakly coupled theory, can be written as a sum of squares and is manifestly positive for black hole backgrounds. Quite beautifully, the shift in the extremal Q/M ratio is directly proportional to the shift in the on-shell action, guaranteeing that these threshold corrections push Q > M in compliance with the WGC. Our results apply for black holes with or without dilatonic coupling and charged under any number of U(1)s.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions
The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.
Horiuchi, W. [Osaka Metropolitan University (Japan↗
Application of bifurcation theory to axial flow compressor instability
When a compression system goes unstable, the mode of response depends on the operating and system parameters, such as throttle setting and B parameter. Previous numerical work on the model developed by Moore and Greitzer (1986) has provided a limited picture of the parametric effects. Applying bifurcation theory to a single-harmonic version of the model has supplied much more complete information, defining the boundaries of each mode of response in the parameter space. Specifically, this is shown in a plot of B vs. throttle setting, which compares well with the corresponding map produced experimentally. The importance of the shape of the rotating stall characteristic is emphasized. The analysis shows a qualitative difference between classic surge and deep surge.
Mccaughan, F. E.↗
Applying Dynamical Systems Theory to Optimize Libration Point Orbit Stationkeeping Maneuvers for WIND
NASA's WIND mission has been operating in a large amplitude Lissajous orbit in the vicinity of the interior libration point of the Sun-Earth/Moon system since 2004. Regular stationkeeping maneuvers are required to maintain the orbit due to the instability around the collinear libration points. Historically these stationkeeping maneuvers have been performed by applying an incremental change in velocity, or (delta)v along the spacecraft-Sun vector as projected into the ecliptic plane. Previous studies have shown that the magnitude of libration point stationkeeping maneuvers can be minimized by applying the (delta)v in the direction of the local stable manifold found using dynamical systems theory. This paper presents the analysis of this new maneuver strategy which shows that the magnitude of stationkeeping maneuvers can be decreased by 5 to 25 percent, depending on the location in the orbit where the maneuver is performed. The implementation of the optimized maneuver method into operations is discussed and results are presented for the first two optimized stationkeeping maneuvers executed by WIND.
WIND↗
Strategies for simulating the time evolution of Hamiltonian lattice field theories
Simulating the time evolution of quantum field theories given some Hamiltonian H requires developing algorithms for implementing the unitary operator e -iHt . A variety of techniques exist that accomplish this task, with the most common technique used so far being Trotterization, which is a special case of the application of a product formula. However, other techniques exist that promise better asymptotic scaling in certain parameters of the theory being simulated, the most efficient of which are based on the concept of block encoding. In this work we study the performance of such algorithms in simulating lattice field theories. We derive and compare the asymptotic gate complexities of several commonly used simulation techniques in application to Hamiltonian lattice field theories. Using the scalar $\hat{φ}$ 4 theory as a test, we also perform numerical studies and compare the gate costs required by product formulas and signal-processing-based techniques to simulate time evolution. For the latter, we use the linear combination of unitaries (LCU) construction augmented with the quantum Fourier transform circuit to switch between the field and momentum eigenbases, which leads to immediate order-of-magnitude improvement in the cost of preparing the block encoding. Further, this paper also includes a pedagogical review of the techniques used, in particular product formulas, LCU, qubitization, quantum signal processing, as well as the technique for simulating geometrically-local Hamiltonians developed by Haah, Hastings, Kothari, and Low.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
An analysis for high speed propeller-nacelle aerodynamic performance prediction. Volume 2: User's manual
A user's manual for the computer program developed for the prediction of propeller-nacelle aerodynamic performance reported in, An Analysis for High Speed Propeller-Nacelle Aerodynamic Performance Prediction: Volume 1 -- Theory and Application, is presented. The manual describes the computer program mode of operation requirements, input structure, input data requirements and the program output. In addition, it provides the user with documentation of the internal program structure and the software used in the computer program as it relates to the theory presented in Volume 1. Sample input data setups are provided along with selected printout of the program output for one of the sample setups.
Egolf, T. Alan↗
Adaptive control of robotic manipulators
The author presents a novel approach to adaptive control of manipulators to achieve trajectory tracking by the joint angles. The central concept in this approach is the utilization of the manipulator inverse as a feedforward controller. The desired trajectory is applied as an input to the feedforward controller which behaves as the inverse of the manipulator at any operating point; the controller output is used as the driving torque for the manipulator. The controller gains are then updated by an adaptation algorithm derived from MRAC (model reference adaptive control) theory to cope with variations in the manipulator inverse due to changes of the operating point. An adaptive feedback controller and an auxiliary signal are also used to enhance closed-loop stability and to achieve faster adaptation. The proposed control scheme is computationally fast and does not require a priori knowledge of the complex dynamic model or the parameter values of the manipulator or the payload.
Seraji, H.↗
Universality of Rényi Entropy in Conformal Field Theory
We use the thermal effective theory to prove that, for the vacuum state in any conformal field theory in 𝑑 dimensions, the 𝑛th Rényi entropy 𝑆$^{(𝑛)}_{𝐴}$ behaves as 𝑆$^{(𝑛)}_{𝐴}$ = [𝑓/(2𝜋𝑛) 𝑑−1 ][Area(∂𝐴)/(𝑑−2)𝜀 𝑑−2 ](1+𝑂(𝑛)) in the 𝑛 → 0 limit when the boundary of the entanglement domain 𝐴 is spherical with the UV cutoff 𝜀. The theory dependence is encapsulated in the cosmological constant 𝑓 in the thermal effective action. Using this result, we estimate the density of states for large eigenvalues of the modular Hamiltonian for the domain 𝐴. In two dimensions, we can use the hot spot idea, which describes the effective action in the high-temperature limit when the temperature is position-dependent, to derive more powerful formulas valid for arbitrary positive 𝑛. We discuss the difference between two and higher dimensions and clarify the applicability of the hot spot idea. We also use the thermal effective theory to derive an analog of the Cardy formula for boundary operators in higher dimensions.
Conformal field theory↗