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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Predictive Synthesis

Current solid state synthesis intrinsically involves a multidimensional space which is challenging to parametrize and predict. The diversity of extended structures comes from the diversity of basic properties of elements of the Periodic system which may exhibit a variety of bonding modes. The fundamental challenges of prediction of the preparative outcome are further complicated by practical synthetic issues. Current development in computational and experimental methods calls for collaborative efforts to make solid state synthesis more predictable. Overall, this Perspective discusses several steps in this direction, including integration of predictions of synthetic conditions with new structure predictions, widespread in situ studies to obtain a panoramic view of the reaction mechanism, and the creation of a synthetic database to properly document all synthetic efforts, including the unsuccessful ones.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A new Majorana platform in an Fe-As bilayer superconductor

Iron-chalcogenide superconductors have emerged as a promising Majorana platform for topological quantum computation. By combining topological band and superconductivity in a single material, they provide significant advantage to realize isolated Majorana zero modes. However, iron-chalcogenide superconductors, especially Fe(Te,Se), suffer from strong inhomogeneity which may hamper their practical application. In addition, some iron-pnictide superconductors have been demonstrated to have topological surface states, yet no Majorana zero mode has been observed inside their vortices, raising a question of universality about this new Majorana platform. In this work, through angle-resolved photoemission spectroscopy and scanning tunneling microscopy/spectroscopy measurement, we identify Dirac surface states and Majorana zero modes, respectively, for the first time in an iron-pnictide superconductor, CaKFe 4 As 4 . More strikingly, the multiple vortex bound states with integer-quantization sequences can be accurately reproduced by our model calculation, firmly establishing Majorana nature of the zero mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electromagnetic total- f algorithm for gyrokinetic particle-in-cell simulations of boundary plasma in XGC

We report the simplified δf mixed-variable/pullback electromagnetic simulation algorithm implemented in XGC for core plasma simulations by Cole et al. [Phys. Plasmas 28, 034501 (2021)] has been generalized to a total- f electromagnetic algorithm that can include, for the first time, the boundary plasma in diverted magnetic geometry with neutral particle recycling, turbulence, and neoclassical physics. The δf mixed-variable/pullback electromagnetic implementation is based on the pioneering work by Kleiber and Mischenko et al. [Kleiber et al., Phys. Plasmas 23, 032501 (2016); Mishchenko et al., Comput. Phys. Commun. 238, 194 (2019)]. An electromagnetic demonstration simulation is performed in a DIII-D-like, H-mode boundary plasma, including a corresponding comparative electrostatic simulation, which confirms that the electromagnetic simulation is necessary for a higher fidelity understanding of the electron particle and heat transport even at the low- β pedestal foot in the vicinity of the magnetic separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Phase transitions in HfO2 probed by first-principles computations

Ever since ferroelectricity was discovered in HfO2, the question of its origin remains controversial. Here, we probe this question using a combination of Landau theory of phase transitions and first-principles computations. In such an approach, the energy landscape associated with the phase transition between cubic and different experimentally demonstrated phases of HfO2 (tetragonal, monoclinic, orthorhombic Pbca, orthorhombic Pnma, and orthorhombic Pca21) is explored using density functional theory calculations. Computations revealed that stabilization of all but orthorhombic Pbca phase is driven by a single unstable zone-boundary antipolar mode X2−. When coupled with zone-center modes (Γ1+ and Γ3+), it stabilizes the tetragonal phase. Coupling with four additional modes (Γ5+, X3−, X5−, X5+) results in the monoclinic phase, which is the ground state of the material. If, however, Γ5+ mode is replaced with Γ4− mode, orthorhombic polar phase Pca21 is stabilized. The application of this framework to examine the effect of electric field on the ferroelectric phase of hafnia reveals that the field of 5 MV/cm is capable of stabilizing ferroelectric phase over the monoclinic one at 0 K.

Physics↗
DOE OSTI2023

Parity-odd and even trispectrum from axion inflation

The four-point correlation function of primordial scalar perturbations has parity-even and parity-odd contributions and the parity-odd signal in cosmological observations is opening a novel window to look for new physics in the inflationary epoch. We study the distinct parity-odd and even prediction from the axion inflation model, in which the inflaton couples to a vector field via a Chern-Simons interaction, and the vector field is considered to be either approximately massless (m A << Hubble scale H) or very massive (m A ~ H). The parity-odd signal arises due to one transverse mode of the vector field being predominantly produced during inflation. We adopt the in-in formalism to evaluate the correlation functions. Considering the vector field mode function to be dominated by its real part up to a constant phase, we simplify the formulas for numerical computations. The numerical studies show that the massive and massless vector fields give significant parity-even signals, while the parity-odd contribution is about one to two orders of magnitude smaller.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Index Theorems, Generalized Hall Currents, and Topology for Gapless Defect Fermions

In this work, we show how the index of the fermion operator from the Euclidean action can be used to uncover the existence of gapless modes living on defects (such as edges and vortices) in topological insulators and superconductors. The 1-loop Feynman diagram that computes the index reveals an analog of the quantum Hall current flowing on and off the defect—even in systems without conserved currents or chiral anomalies—and makes explicit the interplay between topology in momentum and coordinate space. We provide several explicit examples.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Instabilities driven by the fundamental crabbing mode

The Electron Ion Collider (EIC) requires crab cavities to correct for crossing angle effects. These cavities produce a time varying horizontal kick that makes the collision head on in the center of mass frame of the bunches. This results in a luminosity that is very close to the the luminosity obtained without a crossing angle. Producing the kick requires a transverse RF field. This field requires a superconducting resonator and reasonable power levels require loaded quality factors of order a million. This results is a very large transverse impedance for the crabbing mode. The horizontal beta function is large and this, coupled with the large transverse impedance of the resonator can lead to strong instabilities. In the first section of the paper a simple model is used to describe the situation without feedback. There are large ranges of tune leading to strong instabilities in both the electron and hadron storage rings of the EIC. These estimates are confirmed with simulations. In the second section, RF feedback is used to reduce the apparent impedance of the crabbing mode. Simulating the actual feedback would require both a significant rewrite of the tracking code and a significant reduction in computational speed. Instead, the apparent impedance is fitted using one pole filters, which allows use of the existing simulation code. In the third section a transverse damping system is outlined.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Towards Universal Non-Dimensional Characterization of the Oscillatory Dynamics of Wind Turbine Rotors of Multiple Sizes

One concern in the field of Horizontal Axis Wind Turbines (HAWTs) is what control strategies are needed to handle gust pulses in the wind to prevent extreme oscillations of the blades to reduce fatigue stress, prevent blade rupture, and extend the turbine’s operational life. In order to design innovative control strategies to modify the blade’s oscillatory response, it is crucial to establish the fundamental vibrational behavior of the blades when excited by gust pulses of different frequencies and amplitudes present in the fluctuating wind inflow. In a series of previous works, the authors presented a novel Reduced-Order Characterization (ROC) technique that provided an energy-based characterization of the fundamental modes of oscillation of wind turbine rotors when excited by combinations of wind gust pulses of different frequencies and amplitudes. The main focus of the present work is to extend these original notions of energy-based ROC to a universal technique expressed in terms of non-dimensional quantities that could be applied to turbines of any size, operating in any set of wind conditions, as long as they share geometrical and material similarity. The ROC technique provides a simple formula that is capable of predicting the dominant vibrational modes of a blade with sufficient precision to be useful in the determination of a control decision that can be computed in real time, an aspect of fundamental importance in dealing with rapid fluctuations in wind conditions.

Yates, North (ORCID:000900052204261X)↗
DOE OSTI2026

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗
DOE OSTI2025

Combining Theory and Experiment to Map the Atomic-Level Structure–Energy Pathways of Adsorbate-Mediated Phase Changes in a Cooperatively Flexible Metal–Organic Framework

An important subclass of metal–organic frameworks (MOFs) exhibits cooperative flexibility, wherein individual crystallites undergo global structural phase changes in response to external stimuli. Where cooperative flexibility results in reversible changes between crystalline states of distinct accessible porosity, these frameworks can exhibit rare yet desirable behaviors that cannot be explained by local dynamics alone. Yet, the chemical and structural origins of cooperative flexibility and how frameworks undergo these reversible phase changes at the atomic level remain poorly understood. Deliberate design for specific applications is therefore exceedingly difficult, and there is great impetus to develop a fundamental understanding of this phenomenon. Here, an effective and widely accessible computational approach is developed, which is designed to provide microscopic resolution via direct comparison to experimental data along the desorption-guided pathway. The strategy is applied to explain the desorption-induced phase change in an experimentally well-characterized framework, CdIF-13 (sod-Cd(benzimidazolate)2), where experiment alone was unable to resolve the atomistically detailed phase change landscape. Our findings reveal that the cooperative phase change pathways are adsorbate dependent with thermodynamics of intermediate structural states dictated by a nuanced interplay of ligand orientation, skeletal symmetry, and modes of surface adsorption. The results reveal that this isotropically flexible framework is “chaperoned” through a complex energy landscape by specific adsorbates, revealed by the reported computational approach with atomic-level insight and validated by experimentally determined structures. Thus, this work facilitates both understanding and future design of flexible materials for applications in gas storage, transport, delivery, and separation technologies.

03 NATURAL GAS↗
DOE OSTI2025

MMMnet: A Neural Network Surrogate for Real-Time Transport Prediction Based on the Updated Multi-Mode Model

The Multi-Mode Model (MMM) is a physics-based anomalous transport model integrated into TRANSP for predicting electron and ion thermal transport, electron and impurity particle transport, and toroidal and poloidal momentum transport. While MMM provides valuable predictive capabilities, its computational cost, although manageable for standard simulations, is too high for real-time control applications. MMMnet, a neural network-based surrogate model, is developed to address this challenge by significantly reducing computation time while maintaining high accuracy. Trained on TRANSP simulations of DIII-D discharges, MMMnet incorporates an updated version of MMM (9.0.10) with enhanced physics, including isotopic effects, plasma shaping via effective magnetic shear, unified correlation lengths for ion-scale modes, and a new physics-based model for the electromagnetic electron temperature gradient mode. A key advancement is MMMnet’s ability to predict all six transport coefficients, providing a comprehensive representation of plasma transport dynamics. MMMnet achieves a two-order-of-magnitude speed improvement while maintaining strong correlation with MMM diffusivities, making it well-suited for real-time tokamak control and scenario optimization.

DIII-D↗
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Tech-X contributions to the Center for Tokamak Transients (CTTS) Scientific Discovery through Advanced Computing (SciDAC) project (Final Report)

The Center for Tokamak Transients Simulations (CTTS) consists of computational plasma physicists, computer scientists, and applied mathematicians who are together developing and applying simulation software to improve our understanding of and ability to predict and control transient events in tokamak discharges, particularly those that lead to violent disruptions. We have major simulation activities concerning disruptions that are caused by violating ideal MHD stability criteria, vertical displacement events, resistive wall modes, and tearing modes. We are performing comprehensive modeling of disruption mitigation techniques, particularly by shattered pellet injection (SPI). The Center is built around two mature extended MHD codes: NIMROD and M3D-C1. These codes are being extended in a number of ways as required by the demands of disruption modeling, and are being optimized to work well for larger problem size and modern computer architectures as required to model tokamak plasmas with realistic parameters over the required timescales. The codes will be coupled to the FronTier code for state-of-the-art modeling of SPI. Validation with experimental data on DIII-D, NSTX, and JET is being facilitated by rigorous uncertainty quantification (UQ). The high-level objectives of the center are to: (1) Develop better understanding and improved predictive capability of when crossing a linear stability boundary will likely lead to a disruption, (2) Provide a tool for calculating vessel forces, both axisymmetric and non-axisymmetric, for worst case VDE and other disruptions in tokamaks, (3) Develop increased understanding and improved predictive capability for locked mode disruptions and how best to avoid them, (4) Develop, verify, and validate 3D models for disruption mitigation by shattered pellets that can be used to design and optimize a system for future tokamaks, and (5) Increase the efficiency and scalability of the M3D-C1 and NIMROD codes so that they can model disruptions in ITER scale devices using the most powerful DOE HPC hardware, and (6) Demonstrate code-coupling via FronTier coupling with NIMROD and M3D-C1 and by interfacing these codes within a WDM framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Effects of external kink and fishbone-like modes on energetic particle transport in tokamak plasmas

Transport and loss of beam injected energetic particles (EPs) due to three-dimensional perturbations, associated with the external kink (XK) instability and fishbone-like mode (FLM), are numerically investigated utilizing the guiding center following code ORBIT for static toroidal plasmas in HL-2A. The perturbation structure for the XK is computed by the MARS-F code and then mapped to the Boozer coordinates as defined in ORBIT. The simulation shows that the EP profile experiences a significant change in the middle of the plasma column, when the XK-induced radial magnetic field perturbation amplitude, normalized by the equilibrium field, exceeds a threshold value of about 10 -2 . The EP transport is found to be dominated by a diffusion process instead of convection. Furthermore, by scanning the perturbation frequency as a free parameter while maintaining the XK mode structure (thus mimicking the FLM as observed in DIII-D and JT-60U tokamaks), redistribution and loss of EPs are found to be substantially enhanced due to strong resonances between the FLM and EPs, when the mode frequency exceeds a threshold value of ~2 kHz for the case considered. For either XK or FLM, the response of passing EPs to the perturbation is dominant due to the assumed tangential neutral beam injection. Most lost EPs due to these instabilities are initially passing particles but are eventually lost through trapped orbits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Response of thermal and fast-ion transport to beam ion population, rotation and T e / T i in the DIII-D steady state hybrid scenario

The thermal and fast-ion transport properties of DIII-D steady-state hybrid discharges with normalized beta β N ≳ 3 are studied at low injected torque and an increased electron to ion temperature ratio T e / T i . Linear stability analysis performed with the TGLF turbulent code indicates that a high-k mode is usually dominant at smaller radii, whereas a low-k mode is usually dominant at larger radii in these plasmas. A reduction in the net injected torque from 8.6 to 4.3 N-m leads to reduced E × B shear and hence, an enhanced turbulence that was observed on the Doppler backscattering diagnostic and was also computed with TGLF. As T e / T i in the core was increased from 0.57 to 0.66 by adding electron cyclotron current drive (ECCD) to these plasmas, higher levels of transport are observed with increased high-k modes indicated by TGLF. The fast-ion transport level varied over an order of magnitude in these discharges depending on whether Alfvén eigenmodes, fishbones, or no instabilities were observed. Hybrid plasmas with fishbones have decreased fast-ion transport, compared to plasmas with Alfvén eigenmodes, since they are resonant with a smaller portion of phase space and their resonance is farther from the wall. This reduction in fast-ion transport with ECCD mitigates the increase in turbulent transport, resulting in higher performance than expected during strong electron heating. Similarly, the lowest fast-ion transport was observed in the low torque plasma, which also led to better than expected performance at this torque value. Finally, the thermal and fast-ion transport changes observed as the torque/rotation and T e / T i are varied indicate possible methods for transferring this scenario to a reactor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A general approach to seismic inversion with automatic differentiation

Imaging Earth structure or seismic sources from seismic data involves minimizing a target misfit function, and is commonly solved through gradient-based optimization. The adjoint-state method has been developed to compute the gradient efficiently; however, its implementation can be time-consuming and difficult. We develop a general seismic inversion framework to calculate gradients using reverse-mode automatic differentiation. The central idea is that adjoint-state methods and reverse-mode automatic differentiation are mathematically equivalent. Here, the mapping between numerical PDE simulation and deep learning allows us to build a seismic inverse modeling library, ADSeismic, based on deep learning frameworks, which supports high performance reverse-mode automatic differentiation on CPUs and GPUs. We demonstrate the performance of ADSeismic on inverse problems related to velocity model estimation, rupture imaging, earthquake location, and source time function retrieval. ADSeismic has the potential to solve a wide variety of inverse modeling applications within a unified framework.

58 GEOSCIENCES↗
DOE OSTI2023

Shedding Light on the Vibrational Signatures in Halogen‐Bonded Graphitic Carbon Nitride Building Blocks

Abstract The relative contributions of halogen and hydrogen bonding to the interaction between graphitic carbon nitride monomers and halogen bond (XB) donors containing C−X and C≡C bonds were evaluated using computational vibrational spectroscopy. Conventional probes into select vibrational stretching frequencies can often lead to disconnected results. To elucidate this behavior, local mode analyses were performed on the XB donors and complexes identified previously at the M06‐2X/aVDZ‐PP level of theory. Due to coupling between low and high energy C−X vibrations, the C≡C stretch is deemed a better candidate when analyzing XB complex properties or detecting XB formation. The local force constants support this conclusion, as the C≡C values correlate much better with the σ ‐hole magnitude than their C−X counterparts. The intermolecular local stretching force constants were also assessed, and it was found that attractive forces other than halogen bonding play a supporting role in complex formation.

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

Saddle point sampling using scaled normal coordinates

The sampling of saddle points (SPs) on a potential energy surface (PES) is critical for describing the dynamics and transport properties of materials. In this paper we propose a saddle point search (SPS) method that enables efficient sampling of practically meaningful SPs around a local minimum by converting the PES to the scaled normal coordinates (SNCs), together with the dimer method. We find that the pathway leading to a SP using SNCs is highly efficient and is independent of the system size. This results in the significantly increased SPS efficiency and the enhanced possibility of a complete catalog of the SPs. In addition, we perform SPSs using the SNCs for the diffusion of a vacancy and a dumbbell in body-centered cubic Fe and demonstrate (1) the required number of system force evaluations for a SPS decreases by at least an order of magnitude, (2) the effects of the system size on the number of force evaluations are greatly decreased, and (3) the effects of the number of atoms involved in the SPS on the probability of finding practically meaningful SPs is mostly eliminated. We also discuss the computational cost of introducing the SNCs. It should be highlighted that the SNCs could also be applied to other minimum-mode following methods, demonstrating the general versatility of the proposed method.

36 MATERIALS SCIENCE↗