Search NASA⌕ Search

SEARCH · Search NASA

Results for “rotational transform”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Vibrational Cooper minima used to verify the dissociation pathway induced by a short intense laser pulse

Vibrational Cooper minima due to weak dipole coupling were previously observed in strong-field dissociation of H 2⁡ + as minima in the kinetic-energy-release (KER) spectrum [McKenna et al., Phys. Rev. Lett. 103, 103006 (2009)]. We demonstrate in two independent experiments that these vibrational Cooper minima occur more generally and can be used as “fingerprints” to identify the final electronic state in the dissociation pathway. In the first experiment, we show the wavelength dependence of the vibrational Cooper minima observed in the KER spectrum of O 2⁢ + dissociation induced by 396- and 264-nm photons. These minima appear at the locations predicted by our first-order perturbation-theory calculations. In the second experiment, a delay-dependent KER spectrum is obtained from O 2 using an ionizing pump and a dissociating probe, then Fourier-transformed to generate a spin-rotation quantum beat spectrum of O 2 ⁢+ . As a result, this spectrum shows vibrational Cooper minima at the expected locations, confirming their presence in another observable. In this case, these minima are used to identify the dissociation pathway.

Molecular dissociation↗
DOE OSTI2020

X Fe 4 Ge 2 ( X = Y , Lu ) and Mn 3 Pt : Filling-enforced magnetic topological metals

Magnetism, coupled with nontrivial band topology, can bring about many interesting and exotic phenomena, so that magnetic topological materials have attracted persistent research interest. However, compared with nonmagnetic topological materials (TMs), the magnetic TMs are less studied, since their magnetic structures and topological phase transitions are usually complex and the first-principles predictions are usually sensitive on the effect of Coulomb interaction. Here we present a comprehensive investigation of XFe 4 Ge 2 (X = Y, Lu) and Mn 3 Pt, and find these materials to be filling-enforced magnetic topological metals. Our first-principles calculations show that XFe 4 Ge 2 (X = Y, Lu) host Dirac points near the Fermi level at high symmetry point S. These Dirac points are protected by $\text{P}\mathcal{T}$ symmetry ($\text{P}$ and $\mathcal{T}$ are inversion and time-reversal transformations, respectively) and a twofold screw rotation symmetry. Moreover, through breaking $\text{P}\mathcal{T}$ symmetry, the Dirac points would split into Weyl nodes. Mn 3 Pt is found to host fourfold degenerate band crossings in the whole high symmetry path of A – Z . We also utilize the GGA + U scheme to take into account the effect of Coulomb repulsion and find that the filling-enforced topological properties are naturally insensitive on U .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

FLIT: A Generic Fortran Library based on Interfaces and Templates

This Fortran code consists of multiple modules with a focus on simplifying array operations, image processing, and numerical computation especially for computational geophysics applications. We intend to use this code to demonstrate the application and usefulness of Fortran interface and templates for generic programming, especially for computational geophysics and seismology applications. The code has several notable features. Firstly, it is based on a modularized structure, where each module contains multiple functions but with a focus of functionality. Secondly, it heavily uses interfaces and templates for improving the genericness and convenience of the resulting code, where a same function interface can enclose a group of functions that perform the same functionality but with inputs/output variables of different data types. Thirdly, it includes a variety of generic functions with an emphasis on array operations, such as rotation, flipping, cropping, padding, fast Fourier transform, Gaussian blurring, interpolation, and so on. We name this package FLIP – a generic Fortran Library based on Interfaces and Templates.

Gao, Kai↗
DOE OSTI2023

The C 2 H 4 O isomers in the oxidation of ethylene

A combined rotational spectroscopy, thermochemistry, and kinetic modeling study explores the mechanism of acetaldehyde (ethanal), vinyl alcohol (ethenol), and ethylene oxide (oxirane) formation in the oxidation of ethylene (ethene). Multiplexed quantitative detection of oxidation intermediates and products exiting a SiO 2 / SiC microreactor at 1700 K is demonstrated with the BrightSpec W-band chirped-pulse Fourier transform millimeter-wave spectrometer. The broadband rotational spectrum contains transitions of formaldehyde (CH 2 O), methoxy (CH 3 O), methanol (CH 3 OH), ketene (CH 2 CO), acetaldehyde (CH 3 CHO), syn-vinyl alcohol (syn- CH 2 CHOH), anti-vinyl alcohol (anti-CH 2 CHOH), oxirane (c-CH 2 OCH 2 ), propyne (CH 3 CCH), and syn-propanal (syn-CH 3 CH 2 CHO). We focus on the three C 2 H 4 O species and deduce their concentration ratio [CH 3 CHO]:[CH 2 CHOH]:[c-CH 2 OCH 2 ] = 1:0.7(2):0.06(2) by comparing the observed line intensities to those simulated with the PGOPHER software. Detailed thermochemistry of the C 2 H 4 O isomers and conformers is provided by Active Thermochemical Tables. The observed excess concentrations of vinyl alcohol and oxirane relative to the more stable acetaldehyde compared to the equilibrium ratio at 1700 K (1:0.087:0.000024), point to direct chemical pathways to these higher energy isomers. The mechanism for the formation of the three C 2 H 4 O isomers is analyzed using a 0-D homogenous reactor kinetics simulation for ethylene oxidation. The ratios of the C 2 H 4 O isomers concentrations predicted by the kinetic model are compared to the experimental values.

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

Getting to the Core of PARAFAC2, A Nonnegative Approach

In this paper, the authors present a novel method of performing PARAFAC2 factorization of three-way data using a compact representation of that data. In the standard PARAFAC2 algorithm, two modes of the data are recovered directly during the decomposition while the third mode is returned as a transformation matrix, which is then used to rotate sets of orthogonal third-mode basis factors into interpretable factors. In our new method, the data are first decomposed into a core matrix and orthogonal factor loading matrices in the first two modes as well as sets of orthogonal factors in the third mode (as in standard PARAFAC2). The core matrix is then decomposed using a the standard PARAFAC2 strategy to produce transformation matrices in all three modes. The algorithm is particularly useful for very large data sets and essentially permits imposition of nonnegativity in all three modes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Effective Model for Fractional Topological Corner Modes in Quasicrystals

High-order topological insulators (HOTIs), as generalized from topological crystalline insulators, are characterized with lower-dimensional metallic boundary states protected by spatial symmetries of a crystal, whose theoretical framework based on band inversion at special k points cannot be readily extended to quasicrystals because quasicrystals contain rotational symmetries that are not compatible with crystals, and momentum is no longer a good quantum number. In this work, we develop a low-energy effective model underlying HOTI states in 2D quasicrystals for all possible rotational symmetries. By implementing a novel Fourier transform developed recently for quasicrystals and approximating the long-wavelength behavior by their large-scale average, we construct an effective k · p Hamiltonian to capture the band inversion at the center of a pseudo-Brillouin zone. We show that an in-plane Zeeman field can induce mass kinks at the intersection of adjacent edges of a 2D quasicrystal topological insulators and generate corner modes (CMs) with fractional charge, protected by rotational symmetries. Our model predictions are confirmed by numerical tight-binding calculations. Furthermore, when the quasicrystal is proximitized by an s-wave superconductor, Majorana CMs can also be created by tuning the field strength and chemical potential. Our work affords a generic approach to studying the low-energy physics of quasicrystals, in association with topological excitations and fractional statistics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mechanical systems and kinematics

Kinematics is the study of how motion is achieved without reference to the forces that create the motion. This includes the basic structures of all robotic systems, such as links and joints, as well as their mathematical representations. This chapter will provide a brief overview of the types of motion used in additive manufacturing (AM) systems, as well as the mathematical equations that govern them. Rotations are explained using both Euler rotations and quaternions. Combined rotation and translations are shown using homogeneous transform matrices. The robot Jacobian, which relates joint velocities to the end effector translational velocities, is explained with an example calculation. This section will give the reader the basic mathematical foundation in kinematics that is necessary to understand the mechanical underpinnings of path planning for AM systems.

Wang, Peter↗
DOE OSTI2020

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. We present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method.

Zhang, Xin↗
DOE OSTI2022

Phenomenological R -Matrix parameterization of direct, doorway, and compound nuclear reactions [Abstract]

Although formal expressions for scattering matrix accounting for direct, doorway, and compound nuclear (CN) resonant reactions have been derived several decades ago in both the transition ( T -)matrix formalism and the reactance ( K -)matrix formalism, the absence of corresponding expressions in phenomenological R -matrix formalism has limited the application of the latter to CN resonant reactions only. We remove this limitation by parameterizing direct, doorway, and CN resonant reactions in a phenomenological R -matrix scattering matrix, and provide a parameterization for a corresponding Reich-Moore approximation of eliminated capture channels. Direct reactions induce (previously neglected) mixing among the incoming or outgoing R -matrix channel wave functions, parameterized by real and orthonormal channel-rotation matrix, M , whereby the original scattering matrix U is transformed into M T UM . Any real and orthonormal matrix, M , can be equivalently expressed as e η , where η is a real and skew-symmetric 2 rotation-generating matrix that subsequently yields a more intuitive parameterization of eliminated direct capture reactions in Reich-Moore approximation. A phenomenological R -matrix parameterization of doorway reactions is inferred by equating the expression for reactance ( K -)matrix, given in terms of Brune’s alternative R -matrix parameterization, to a corresponding expression derived using Feshbach’s projection operator formalism. Assuming that all doorway states, just like CN states, are confined within spheres defined by R -matrix channel radii, a new R -matrix-like term induced by doorway states is gleaned, wherein each doorway state is parameterized by its energy, width, and the strength of its coupling to each CN state. Since a Reich-Moore approximation for retained-channel scattering matrix ought to approximate the effect of eliminated capture channels taking place via direct, doorway, or CN reactions, each of the three kinds of reactions contributing to the capture entails a corresponding Reich-Moore parameterization in a first-order approximation: direct contribution is parameterized by introducing finite diagonal elements of a retained-channel rotation-generating matrix, doorway contribution is parameterized by doorway capture widths, while CN contribution is parameterized by conventional Reich-Moore capture widths. We will present evidence of direct and doorway reactions observed in recent measurements of resolved resonance cross sections at the Gaerttner LINAC Center at Rensselaer Polytechnic Institute, and will outline a path for implementing this new R -matrix parameterization into the SAMMY nuclear data evaluation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Generalized Quantum Signal Processing

Quantum signal processing (QSP) and quantum singular value transformation (QSVT) currently stand as the most efficient techniques for implementing functions of block-encoded matrices, a central task that lies at the heart of most prominent quantum algorithms. However, current QSP approaches face several challenges, such as the restrictions imposed on the family of achievable polynomials and the difficulty of calculating the required phase angles for specific transformations. In this paper, we present a generalized quantum signal processing (GQSP) approach, employing general SU(2) rotations as our signal-processing operators, rather than relying solely on rotations in a single basis. Our approach lifts all practical restrictions on the family of achievable transformations, with the sole remaining condition being that | P | ≤ 1 , a restriction necessary due to the unitary nature of quantum computation. Furthermore, GQSP provides a straightforward recursive formula for determining the rotation angles needed to construct the polynomials in cases where P and Q are known. In cases where only P is known, we provide an efficient optimization algorithm capable of identifying in under a minute of GPU time, a corresponding Q for polynomials of degree on the order of 10 7 . We further illustrate GQSP simplifies QSP-based strategies for Hamiltonian simulation, offer an optimal solution to the ϵ -approximate fractional query problem that requires O ( ( 1 / δ ) + log ( 1 / ϵ ) ) queries to perform where O ( 1 / δ ) is a proved lower bound, and introduces novel approaches for implementing bosonic operators. Moreover, we propose a novel framework for the implementation of normal matrices, demonstrating its applicability through synthesis of diagonal matrices, as well as the development of a new algorithm for convolution through synthesis of circulant matrices using only O ( d log N + log 2 N ) 1 and 2-qubit gates for a filter of lengths d . Published by the American Physical Society 2024

Motlagh, Danial↗
DOE OSTI2024

System Modeling and Robust Stability Region Analysis for Multi-inverters Based on VSG

As the number of inverter-interfaced distributed generators (IIDGs) increases, the decrease in the inertia of the microgrid affects the stability of microgrid voltage and current. The virtual synchronous generator (VSG) control method simulates the output characteristics of the synchronous machine to improve the inertia and damping of the microgrid. However, low-frequency oscillations can occur in systems with multi-VSGs connected in parallel, primarily when operating in island mode. Therefore, this article develops the small-signal state-space model for multi-inverters based on VSG and analyses the robust stability region (RSR) of the control parameters. Firstly, an extensible small-signal state-space model is established and simplified in the dq axis, which is a necessary preprocessing for the subsequent stability analysis. Secondly, the guardian map method based on bialternate product operation is constructed, and the RSR is established through translation mapping and rotation mapping. Therefore, the solution of RSR is transformed into the problem of identifying the Hurwitz matrix. Finally, the simulations and experiments verify that the calculated region can ensure stable operation during startup and meet the requirement to suppress low-frequency oscillations. In conclusion, the method of solving RSR proposed in this paper provides a basis for selecting controller parameters in practical engineering.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

ACDC (Automated Campbell Diagram Code) [SWR-26-042]

This application provides a web-based graphical user interface to generating Campbell Diagrams and visualizing mode shapes for OpenFAST turbine models. Determining the aeroelastic stability and dynamic characteristics of wind turbines is a critical step in turbine design and analysis. Historically, extracting natural frequencies and mode shapes from OpenFAST—the industry-standard whole-turbine simulation code—has been a fragmented and tedious process. It required manual model configuration, command-line linearization execution, and complex post-processing via proprietary scripts to handle rotating-frame dynamics. To address these workflow bottlenecks, we present the Automated Campbell Diagram Code (ACDC), an open-source graphical software tool developed by the National Laboratory of the Rockies (NLR) under the DOE-funded Distributed Wind Aeroelastic Modeling (dWAM) project. ACDC streamlines the end-to-end linearization and stability analysis workflow into a single, intuitive cross-platform application. The software guides users through OpenFAST model configuration, definition of operating points, and the automated execution of steady-state trim and linearization simulations. Under the hood, ACDC automates the complex mathematical post-processing steps required for rotating systems, including Multi-Blade Coordinate (MBC) transformations, eigenanalysis, and advanced modal tracking utilizing the Modal Assurance Criterion (MAC) and spectral clustering. Finally, ACDC processes these results to automatically generate Campbell diagrams and features a robust 3D visualization engine to animate full-system mode shapes. By eliminating the reliance on external post-processing environments and manual data manipulation, ACDC significantly accelerates dynamic analysis and lowers the barrier to entry for wind energy researchers and engineers.

Summerville, Brent [National Laboratory of the Roc↗
DOE OSTI2021

Understanding the phase transformation mechanisms that affect the dynamic response of Fe-based microstructures at the atomic scales

Large-scale molecular dynamics (MD) simulations were carried out to investigate the shock-induced evolution of microstructure in Fe-based systems comprising single-crystal and layered Cu/Fe alloys with a distribution of interfaces. The shock compression of pure single-crystal Fe oriented along [110] above a threshold pressure results in a BCC (α)→HCP (ε) phase transformation behavior that generates a distribution of ε phase variants in the phase transformed region of the microstructure behind the shock front. The propagation of the release wave through a phase transformed ε phase causes a reverse ε→α phase transformation and renders a distribution of twins for the [110] oriented Fe that serve as void nucleation sites during spall failure. The simulations reveal that the α→ε→α transformation-induced twinning for shock loading along the [110] direction is due to a dominant ε phase variant formed during compression that rotates on the arrival of the release wave followed by a reverse phase transformation to twins in the α phase. The modifications in the evolution of the ε phase variants and twins in Fe behavior are also studied for Cu–Fe layered microstructures due to the shock wave interactions with the Cu/Fe interfaces using a newly constructed Cu–Fe alloy potential. Here, the MD simulations suggest that interfaces affect the observed variants during shock compression and, hence, distributions of twins during shock release that affects the void nucleation stresses in the Fe phase of Cu/Fe microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kepler and TESS Observations of PG 1159-035

Abstract PG 1159-035 is the prototype of the PG 1159 hot (pre-)white dwarf pulsators. This important object was observed during the Kepler satellite K2 mission for 69 days in 59 s cadence mode and by the TESS satellite for 25 days in 20 s cadence mode. We present a detailed asteroseismic analysis of those data. We identify a total of 107 frequencies representing 32 ℓ = 1 modes, 27 frequencies representing 12 ℓ = 2 modes, and eight combination frequencies. The combination frequencies and the modes with very high k values represent new detections. The multiplet structure reveals an average splitting of 4.0 ± 0.4 μ Hz for ℓ = 1 and 6.8 ± 0.2 μ Hz for ℓ = 2, indicating a rotation period of 1.4 ± 0.1 days in the region of period formation. In the Fourier transform of the light curve, we find a significant peak at 8.904 ± 0.003 μ Hz suggesting a surface rotation period of 1.299 ± 0.002 days. We also present evidence that the observed periods change on timescales shorter than those predicted by current evolutionary models. Our asteroseismic analysis finds an average period spacing for ℓ = 1 of 21.28 ± 0.02 s. The ℓ = 2 modes have a mean spacing of 12.97 ± 0.4 s. We performed a detailed asteroseismic fit by comparing the observed periods with those of evolutionary models. The best-fit model has T eff = 129, 600 ± 11 100 K, M * = 0.565 ± 0.024 M ⊙ , and log g = 7.41 − 0.54 + 0.38 , within the uncertainties of the spectroscopic determinations. We argue for future improvements in the current models, e.g., on the overshooting in the He-burning stage, as the best-fit model does not predict excitation for all of the pulsations detected in PG 1159-035.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Lorentz Symmetry Group, Retardation, Intergalactic Mass Depletion and Mechanisms Leading to Galactic Rotation Curves

The general theory of relativity (GR) is symmetric under smooth coordinate transformations, also known as diffeomorphisms. The general coordinate transformation group has a linear subgroup denoted as the Lorentz group of symmetry, which is also maintained in the weak field approximation to GR. The dominant operator in the weak field equation of GR is thus the d’Alembert (wave) operator, which has a retarded potential solution. Galaxies are huge physical systems with dimensions of many tens of thousands of light years. Thus, any change at the galactic center will be noticed at the rim only tens of thousands of years later. Those retardation effects are neglected in the present day galactic modelling used to calculate rotational velocities of matter in the rims of the galaxy and surrounding gas. The significant differences between the predictions of Newtonian instantaneous action at a distance and observed velocities are usually explained by either assuming dark matter or by modifying the laws of gravity (MOND). In this paper, we will show that, by taking general relativity seriously without neglecting retardation effects, one can explain the radial velocities of galactic matter in the M33 galaxy without postulating dark matter. It should be stressed that the current approach does not require that velocities v are high; in fact, the vast majority of galactic bodies (stars, gas) are substantially subluminal—in other words, the ratio of $\frac{v}{c} \ll 1$. Typical velocities in galaxies are 100 km/s, which makes this ratio 0.001 or smaller. However, one should consider the fact that every gravitational system, even if it is made of subluminal bodies, has a retardation distance, beyond which the retardation effect cannot be neglected. Every natural system, such as stars and galaxies and even galactic clusters, exchanges mass with its environment, for example, the sun loses mass through solar wind and galaxies accrete gas from the intergalactic medium. This means that all natural gravitational systems have a finite retardation distance. The question is thus quantitative: how large is the retardation distance? For the M33 galaxy, the velocity curve indicates that the retardation effects cannot be neglected beyond a certain distance, which was calculated to be roughly 14,000 light years; similar analysis for other galaxies of different types has shown similar results. We demonstrate, using a detailed model, that this does not require a high velocity of gas or stars in or out of the galaxy and is perfectly consistent with the current observational knowledge of galactic and extra galactic material content and dynamics

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Redox Chemistry of the Subphases of α-CsPbI2Br and β-CsPbI2Br: Theory Reveals New Potential for Photostability

The logic in the design of a halide-mixed APb(I 1-x Br x ) 3 perovskite is quite straightforward: to combine the superior photovoltaic qualities of iodine-based perovskites with the increased stability of bromine-based perovskites. However, even small amounts of Br doped into the iodine-based materials leads to some instability. In the present report, using first-principles computations, we analyzed a wide variety of α-CsPbI 2 Br and β-CsPbI 2 Br phases, compared their mixing enthalpies, explored their oxidative properties, and calculated their hole-coupled and hole-free charged Frenkel defect (CFD) formations by considering all possible channels of oxidation. Nanoinclusions of bromine-rich phases in α-CsPbI 2 Br were shown to destabilize the material by inducing lattice strain, making it more susceptible to oxidation. The uniformly mixed phase of α-CsPbI 2 Br was shown to be highly susceptible towards a phase transformation into β-CsPbI 2 Br when halide interstitial or halide vacancy defects were introduced into the lattice. The rotation of PbI 4 Br 2 octahedra in α-CsPbI2Br allows it either to transform into a highly unstable apical β-CsPbI 2 Br, which may phase-segregate and is susceptible to CFD, or to phase-transform into equatorial β-CsPbI 2 Br, which is resilient against the deleterious effects of hole oxidation (energies of oxidation >0 eV) and demixing (energy of mixing <0 eV). Thus, the selective preparation of equatorial β-CsPbI 2 Br offers an opportunity to obtain a mixed perovskite material with enhanced photostability and an intermediate bandgap between its constituent perovskites.

36 MATERIALS SCIENCE↗