Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational modes”

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 325 records · Page 18

Data-Centric Approach to Capture Non-Polynomial Nonlinear Dynamics

We propose an analytical construction of observable functions in the extended dynamic mode decomposition (EDMD) algorithm. EDMD is a numerical method for approximating the spectral properties of the Koopman operator. The choice of observable functions is fundamental for applying EDMD to nonlinear problems arising in systems and control. Existing methods either start from a set of dictionary functions and look for the subset that best fits the underlying nonlinear dynamics or rely on machine learning algorithms to “learn” observable functions. Conversely, in this paper, we start from the dynamical system model and lift it through the Lie derivatives, rendering it into a polynomial form. This proposed transformation into a polynomial form is exact and provides an adequate set of observable functions. The strength of the proposed approach is its applicability to a broader class of nonlinear dynamical systems, particularly those with nonpolynomial functions and compositions thereof. Moreover, it retains the physical interpretability of the underlying dynamical system and can be readily integrated into existing numerical libraries. We demonstrate the proposed approach with an application to electric power systems. The modeled system consists of a single generator connected to an infinite bus, where nonlinear terms include sine and cosine functions. The results demonstrate the effectiveness of the proposed procedure in off-attractor nonlinear dynamics for estimation and prediction; the observable functions obtained from the proposed construction outperform methods that use dictionary functions comprising monomials or radial basis functions.

extended dynamic mode decomposition↗

Synergy of coupled kink and tearing modes in fast ion transport

One aspect of the interaction between fast ions and magnetohydrodynamic (MHD) instabilities is the fast ion transport. Coupled kink and tearing MHD instabilities have also been reported to cause fast ion transport. Recently, the 'kick' model has been developed to compute the evolution of the fast ion distribution from neutral beam injection using instabilities as phase-space resonance sources. Overall, the goal of this paper is to utilize the kick model to understand the physics of fast ion transport caused by the coupled kink and tearing modes. Soft x-ray diagnostics are used to identify the mode parameters in the National Spherical Torus Experiment. The comparison of neutron rates measured and computed from time-dependent TRANSP simulation with the kick model shows that the coupling of kink and tearing mode is important in determination of the fast ion transport. The numerical scan of the mode parameters shows that the relative phase of the kink and tearing modes and the overlapping of kink and tearing mode resonances in the phase space can affect the fast ion transport, suggesting that the synergy of the coupled modes may be causing the fast ion transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Provably accurate simulation of gauge theories and bosonic systems

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ − 1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ = polylog ( ϵ − 1 ) .

Tong, Yu↗

A new efficient grain growth model using a random Gaussian-sampled mode filter

This paper presents the use of a Gaussian neighborhood mode filter for predicting grain growth in a manner similar to the solutions obtained by a Monte Carlo Potts model. This flexible grain growth model can quickly utilize modern, computationally optimized data science strategies on graphics processing units to simulate grain growth up to 100 times faster than the state-of-the-art, publicly available Monte Carlo Potts model. We show that, given the correct neighborhood, the mode filter can replicate normal grain growth in two or three dimensions. In addition, the paper briefly demonstrates the ability to model limited anisotropic in grain boundary energy and mobility. Anisotropic grain boundary energy is modeled by defining a weighted mode filter operation. Anisotropic grain boundary mobility is modeled by scaling and orienting the Gaussian neighborhood in a particular direction.

Anisotropy↗

Mechanistic studies of small molecule ligands selective to RNA single G bulges

Abstract Small-molecule RNA binders have emerged as an important pharmacological modality. A profound understanding of the ligand selectivity, binding mode, and influential factors governing ligand engagement with RNA targets is the foundation for rational ligand design. Here, we report a novel class of coumarin derivatives exhibiting selective binding affinity towards single G RNA bulges. Harnessing the computational power of all-atom Gaussian accelerated molecular dynamics simulations, we unveiled a rare minor groove binding mode of the ligand with a key interaction between the coumarin moiety and the G bulge. This predicted binding mode is consistent with results obtained from structure-activity relationship studies and transverse relaxation measurements by nuclear magnetic resonance spectroscopy. We further generated 444 molecular descriptors from 69 coumarin derivatives and identified key contributors to the binding events, such as charge state and planarity, by lasso (least absolute shrinkage and selection operator) regression. Our work deepened the understanding of RNA-small molecule interactions and integrated a new framework for the rational design of selective small-molecule RNA binders.

Biochemistry & Molecular Biology↗

Cosmological collider signatures of Higgs-R 2 inflation

We study the cosmological collider signatures in the Higgs-R 2 inflation model. We consider two distinct types of signals: one originating from the inflaton coupling to Standard Model fermions and gauge bosons, and another arising from the isocurvature mode interaction with the inflaton. In the former case, we determine that the signal magnitude is likely too small for detection by upcoming probes, primarily due to suppression by both the Planck scale and slow-roll parameters. However, we provide a detailed computation of the signal which could be potentially applicable to various Higgs inflation variants. For the isocurvature mode signals, we observe that the associated couplings remain unsuppressed when the isocurvature mode is relatively light or comparable to the inflationary scale. In this case, we study the Higgs-R 2 inflation parameter space that corresponds to the quasi-single-field inflation regime and find that the signal strength could be as large as |f NL | > 1, making Higgs-R 2 inflation a viable candidate for observation by future 21-cm surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗

Dynamic Frequency Regulation Improvement in Hydropower-Hybrid System using Variational Mode Decomposition

Providing frequency response, especially fast frequency response such as Red D in the Pennsylvania-New Jersey-Maryland Interconnection (PJM) market, is challenging for many generation plants to deliver on their own. If they have adequate flexibility, hydropower plants are typically able to provide slower regulation support (e.g., Reg A in PJM), but do not respond fast enough to provide Reg D. The ability to provide Reg D would improve their revenue because this service is typically more valuable. This work presents a control approach to use hydropower, battery, and ultracapacitor systems to provide fast regulation in a way that uses the response contribution of each asset. In particular, the proposed control architecture applies a variational mode decomposition (VMD) technique on the incoming Reg D signal to extract multiple dynamic-regulation components with non-overlapping frequencies. With the response-speed dependent alignment, these regulation components are fed to the hydrogenator and hybrid energy storage system (HESS). The paper evaluates the proposed approach on a direct-current (DC)-coupled active system by computing several performance measures and analyzing sensitivity based on HESS component proportional capacities. The results reveal that the proposed VMD-based signal conditioning performs well and that optimizing the sizing of the battery and ultracapacitor components further enhances performance.

13 HYDRO ENERGY↗

Nonadiabatic excited-state intramolecular proton transfer in 3-hydroxyflavone: S2 state involvement via multi-mode effect

The excited-state intramolecular proton transfer process in 3-hydroxyflavone is investigated based on the computed structural parameters and energetics of stationary points of vibronically coupled S 1 -S 2 potential energy surfaces. A conical intersection close to the Franck-Condon point on S1 is identified. The minimum energy of the conical intersection is found to be near-degenerate with the equilibrium minimum of S 2 . Quantum nuclear wavepacket simulations revealed a small population transfer from the “bright” S 1 to “dark” S 2 on a time scale shorter than the O—H stretching vibrational period. Such a nonadiabatic transition opens up the possibility of new photophysical and photochemical pathways, including the proton transfer via S 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fiber optic computing using distributed feedback

Abstract The widespread adoption of machine learning and other matrix intensive computing algorithms has renewed interest in analog optical computing, which has the potential to perform large-scale matrix multiplications with superior energy scaling and lower latency than digital electronics. However, most optical techniques rely on spatial multiplexing, requiring a large number of modulators and detectors, and are typically restricted to performing a single kernel convolution operation per layer. Here, we introduce a fiber-optic computing architecture based on temporal multiplexing and distributed feedback that performs multiple convolutions on the input data in a single layer. Using Rayleigh backscattering in standard single mode fiber, we show that this technique can efficiently apply a series of random nonlinear projections to the input data, facilitating a variety of computing tasks. The approach enables efficient energy scaling with orders of magnitude lower power consumption than GPUs, while maintaining low latency and high data-throughput.

97 MATHEMATICS AND COMPUTING↗

Performance Potential of Mixed Data Management Modes for Heterogeneous Memory Systems

Many high-performance systems now include different types of memory devices within the same compute platform to meet strict performance and cost constraints. Such heterogeneous memory systems often include an upper-level tier with better performance, but limited capacity, and lower-level tiers with higher capacity, but less bandwidth and longer latencies for reads and writes. To utilize the different memory layers efficiently, current systems rely on hardware-directed, memory -side caching or they provide facilities in the operating system (OS) that allow applications to make their own data-tier assignments. Since these data management options each come with their own set of trade-offs, many systems also include mixed data management configurations that allow applications to employ hardware- and software-directed management simultaneously, but for different portions of their address space. Despite the opportunity to address limitations of stand-alone data management options, such mixed management modes are under-utilized in practice, and have not been evaluated in prior studies of complex memory hardware. In this work, we develop custom program profiling, configurations, and policies to study the potential of mixed data management modes to outperform hardware- or software-based management schemes alone. Our experiments, conducted on an Intel ® Knights Landing platform with high-bandwidth memory, demonstrate that the mixed data management mode achieves the same or better performance than the best stand-alone option for five memory intensive benchmark applications (run separately and in isolation), resulting in an average speedup compared to the best stand-alone policy of over 10 %, on average.

Effler, Chad↗

Computational Spectroscopy of the Cr–Cr Bond in Coordination Complexes

In this work, we report the accurate computational vibrational analysis of the Cr–Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr–Cr vibrational stretching modes, ν(Cr 2 ), have suggested weaker bonding, even for so-called ultrashort Cr–Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr 2 ) with CASPT2 for three well-characterized complexes, Cr 2 (O 2 CCH 3 ) 4 (H 2 O) 2 , Cr 2 (mhp) 4 , and Cr 2 (dmp) 4 . We obtain CASPT2 harmonic ν(Cr 2 ) modes in good agreement with experiment at 282 cm –1 for Cr 2 (mhp) 4 and 353 cm –1 for Cr 2 (dmp) 4 , compute 50 Cr and 54 Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr–Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr 2 ), leading to an anharmonic mode of 272 cm –1 for Cr 2 (mhp) 4 and 333 cm –1 for Cr 2 (dmp) 4 .

36 MATERIALS SCIENCE↗

Editorial: Edge computation and digital distribution networks

The development of digital technologies is penetrating all areas of energy revolution. Based on the in-depth integration of advanced digital technologies, distribution networks are gradually transforming into digital distribution networks (DDNs) with tremendous changes from the structure to the operation mode. DDNs is the digitalized appearance of the physical distribution network, in which ubiquitous connections and massive data are the basic characteristics (Huo et al., 2022). It is an important task to utilize the massive data and propose novel operation modes to construct more efficient and intelligent distribution networks (Jian et al., 2022). Among the advanced digital technologies in DDNs, edge computing has received wide attention (Zhao et al., 2022). It has superior performance in local sensing and intelligent computation, which can effectively relieve huge communication pressure. However, the limited computing resources and the complex computing tasks at the edge side significantly challenge the collaboration of distribution network regulation and advanced digital technologies (Hu et al., 2022). It is necessary to find out proper methods to utilize advanced digital technology to construct DDNs. This Research Topic is organized to introduce the recent progress in the construction, operation and advanced computational methods for DDNs. Finally, seven papers have been accepted, which can be sorted into the following three categories: 1) Evolution and technical features of DDNs, 2) Intelligent operation control of DDNs, 3) Advanced simulation for large-scale DDNs. The three sections below respectively introduce the major research and contributions of the papers covered in each category.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Sliding Mode Network Perimeter Control

Urban traffic congestion is a chronic problem faced by many cities in the US and worldwide. It results in inefficient infrastructure use as well as increased vehicle fuel consumption and emission levels. Congestion is intertwined with delay, as road users waste precious hours on the road, which in turn reduces productivity. Researchers have developed, and continue to design, tools and systems to alleviate this problem. Network perimeter control is one such tool that has been studied extensively. It attempts to control the flow of vehicles entering a protected area to ensure that the congested regime predetermined by the Network Fundamental Diagram (NFD) is not reached. In this paper, a method derived from sliding mode control theory is presented. Its main advantages over proportional-integral controllers include (1) minimal tuning, (2) no linearization of the governing equations, (3) no assumptions with regard to the shape of the NFD, and (4) ability to handle various demand profiles without the need to retune the controller. A sliding mode controller was implemented and tested on a congested grid network. The results show that the proposed controller produces network-wide delay savings and disperses congestion effectively.

42 ENGINEERING↗

Automatic Differentiation of C++ Codes on Emerging Manycore Architectures with Sacado

Automatic differentiation (AD) is a well-known technique for evaluating analytic derivatives of calculations implemented on a computer, with numerous software tools available for incorporating AD technology into complex applications. However, a growing challenge for AD is the efficient differentiation of parallel computations implemented on emerging manycore computing architectures such as multicore CPUs, GPUs, and accelerators as these devices become more pervasive. In this work, we explore forward mode, operator overloading-based differentiation of C++ codes on these architectures using the widely available Sacado AD software package. In particular, we leverage Kokkos, a C++ tool providing APIs for implementing parallel computations that is portable to a wide variety of emerging architectures. Here we describe the challenges that arise when differentiating code for these architectures using Kokkos, and two approaches for overcoming them that ensure optimal memory access patterns as well as expose additional dimensions of fine-grained parallelism in the derivative calculation. We describe the results of several computational experiments that demonstrate the performance of the approach on a few contemporary CPU and GPU architectures. We then conclude with applications of these techniques to the simulation of discretized systems of partial differential equations.

97 MATHEMATICS AND COMPUTING↗

Energetics of the nucleation and glide of disconnection modes in symmetric tilt grain boundaries

Grain boundaries (GBs) evolve by the nucleation and glide of disconnections, which are dislocations with a step character. In this work, motivated by recent success in predicting GB properties such as the shear coupling factor and mobility from the intrinsic properties of disconnections, we develop a systematic method to calculate the energy barriers for the nucleation and glide of individual disconnection modes under arbitrary driving forces and a quasi-2D setting. This method combines tools from bicrystallography to enumerate disconnection modes and the Nudged elastic band (NEB) method to calculate their energetics, yielding minimum energy paths and atomistic mechanisms for the nucleation and glide of each disconnection mode. We apply the method to accurately predict shear coupling factors of $[001]$ symmetric tilt grain boundaries in Cu. Particular attention is paid to the boundaries where the dislocation-based disconnection nucleation model produces incorrect nucleation barriers. We demonstrate that the method can accurately compute energy barriers and predict shear-coupling factors in the low-temperature regime. For certain disconnection modes in which the assumptions underlying our method do not hold, we report upper bounds on the energy barriers for disconnection nucleation and glide. In addition, the NEB trajectories reveal interesting phenomena such as the dissociation of a higher energy mode into lower energy modes, and in some cases, shear coupling being mediated by partial disconnections, wherein the GB structure temporarily changes to a metastable state before reverting back to its original structure. Graphical abstract

36 MATERIALS SCIENCE↗

Inverse Design of Plasma Metamaterial Devices for Optical Computing

We apply inverse-design methods to produce two-dimensional plasma metamaterial (PMM) devices. Backpropagated finite-difference frequency-domain (FDFD) simulations are used to design waveguides and demultiplexers operating under both TE and TM modes. Demultiplexing and waveguiding are demonstrated for devices composed of plasma elements with reasonable plasma densities of approximately 7 GHz, allowing for future in situ training and experimental realization of these designs. We also explore the possible applicability of PMMs to nonlinear boolean operations for use in optical computing. Functionally complete logical connectives (or and and) are achieved in the TM mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

One-step sputtering of MoSSe metastable phase as thin film and predicted thermodynamic stability by computational methods

Abstract We present the fabrication of a MoS 2−x Se x thin film from a co-sputtering process using MoS 2 and MoSe 2 commercial targets with 99.9% purity. The sputtering of the MoS 2 and MoSe 2 was carried out using a straight and low-cost magnetron radio frequency sputtering recipe to achieve a MoS 2−x Se x phase with x = 1 and sharp interface formation as confirmed by Raman spectroscopy, time-of-flight secondary ion mass spectroscopy, and cross-sectional scanning electron microscopy. The sulfur and selenium atoms prefer to distribute randomly at the octahedral geometry of molybdenum inside the MoS 2−x Se x thin film, indicated by a blue shift in the A 1g and E 1 g vibrational modes at 355 cm −1 and 255 cm −1 , respectively. This work is complemented by computing the thermodynamic stability of a MoS 2−x Se x phase whereby density functional theory up to a maximum selenium concentration of 33.33 at.% in both a Janus-like and random distribution. Although the Janus-like and the random structures are in the same metastable state, the Janus-like structure is hindered by an energy barrier below selenium concentrations of 8 at.%. This research highlights the potential of transition metal dichalcogenides in mixed phases and the need for further exploration employing low-energy, large-scale methods to improve the materials’ fabrication and target latent applications of such structures.

36 MATERIALS SCIENCE↗