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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Advanced modeling and simulation of research reactors using dynamic mode decomposition

Full text of publication follows. Due to the ever-increasing safety requirements, the current trend of nuclear reactor analysis is shifting towards high-fidelity multi-physics models, which have a very high computational cost and modelling complexity. As the cost of even a single model run makes it impossible to analyse the behaviour and performance of these models on large-scale commercial plants, it has become even more significant to provide suitable benchmarks to validate and test them extensively. In this sense, research reactors offer a promising solution for the initial validation of high-fidelity models, as they are significantly smaller than commercial reactors and their characteristics are well known. In particular, the reactors of the TRIGA family have been used to assess and validate models and methods for Generation-IV designs, as they have some similar features (such as the dominance of natural convection as cooling mechanism and the difficulties in performing sub-channel analysis using standard codes). Still, the computational requirements of high-fidelity models make them unsuitable for real-time analysis, even following their assessment on research reactors. In this sense, Model Order Reduction (MOR) techniques give an additional strategy to reduce the computational cost of high-fidelity models (whilst preserving sufficient accuracy). In particular, this work focuses on Dynamic Mode Decomposition (DMD), a non-intrusive MOR technique that aims at representing models with explicit temporal dynamics by extracting the time-varying characteristics and the governing structures based only on a set of available data, thus without needing any underlying knowledge of the governing equations. In addition, DMD also computes a low-dimensional surrogate of the dynamic matrix of the system, making it suited for stability analysis and real-time evaluations. This work focuses on the application and validation of the DMD method on the Computational Fluid-Dynamics (CFD) model TRIGA Mark II reactor, also discussing in detail the potentiality of this algorithm as an advanced modelling tool for nuclear reactor analysis. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Harnessing High‐Throughput Computational Methods to Accelerate the Discovery of Optimal Proton Conductors for High‐Performance and Durable Protonic Ceramic Electrochemical Cells

Abstract The pursuit of high‐performance and long‐lasting protonic ceramic electrochemical cells (PCECs) is impeded by the lack of efficient and enduring proton conductors. Conventional research approaches, predominantly based on a trial‐and‐error methodology, have proven to be demanding of resources and time‐consuming. Here, this work reports the findings in harnessing high‐throughput computational methods to expedite the discovery of optimal electrolytes for PCECs. This work methodically computes the oxygen vacancy formation energy (E V ), hydration energy (E H ), and the adsorption energies of H 2 O and CO 2 for a set of 932 oxide candidates. Notably, these findings highlight BaSn x Ce 0.8‐x Yb 0.2 O 3‐δ (BSCYb) as a prospective game‐changing contender, displaying superior proton conductivity and chemical resilience when compared to the well‐regarded BaZr x Ce 0.8‐x Y 0.1 Yb 0.1 O 3‐δ (BZCYYb) series. Experimental validations substantiate the computational predictions; PCECs incorporating BSCYb as the electrolyte achieved extraordinary peak power densities in the fuel cell mode (0.52 and 1.57 W cm −2 at 450 and 600 °C, respectively), a current density of 2.62 A cm −2 at 1.3 V and 600 °C in the electrolysis mode while demonstrating exceptional durability for over 1000‐h when exposed to 50% H 2 O. This research underscores the transformative potential of high‐throughput computational techniques in advancing the field of proton‐conducting oxides for sustainable power generation and hydrogen production.

08 HYDROGEN↗
DOE OSTI2024

Communication Lower Bounds and Optimal Algorithms for Multiple Tensor-Times-Matrix Computation

Multiple tensor-times-matrix (Multi-TTM) is a key computation in algorithms for computing and operating with the Tucker tensor decomposition, which is frequently used in multidimensional data analysis. Here, we establish communication lower bounds that determine how much data movement is required (under mild conditions) to perform the Multi-TTM computation in parallel. The crux of the proof relies on analytically solving a constrained, nonlinear optimization problem. We also present a parallel algorithm to perform this computation that organizes the processors into a logical grid with twice as many modes as the input tensor. We show that, with correct choices of grid dimensions, the communication cost of the algorithm attains the lower bounds and is therefore communication optimal. Finally, we show that our algorithm can significantly reduce communication compared to the straightforward approach of expressing the computation as a sequence of tensor-times-matrix operations when the input and output tensors vary greatly in size.

HBL-inequalities↗
DOE OSTI2024

Resistive hose modes in tokamak runaway electron beams

Beams of energetic runaway electrons are generated during disruptions in tokamaks, and fluid models are used to study their effects on macroscale dynamics. Linear computations of a massless, runaway electron beam coupled to MHD plasma show that resistive hose instabilities grow faster than tearing modes at large resistivity. Eigenvalue results with reduced models of the resistive hose instability are compared with results from the full MHD and beam system, showing that the resistive hose decouples from any plasma response. An estimate of plasma temperature at which growth of the resistive hose dominates tearing for post-disruption DIII-D plasma parameters is in a physically relevant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Experimental Insights into the Coupling of Methane Combustion and Steam Reforming in a Catalytic Plate Reactor in Transient Mode

The microstructured reactor concept is very promising technology to develop a compact reformer for distributed hydrogen generation. In this work, a catalytic plate reactor (CPR) is developed and investigated for the coupling of methane combustion (MC) and methane steam reforming (MSR) over Pt/Al 2 O 3 -coated microchannels in cocurrent and counter-current modes in transient experiments during start-up. A three-dimensional (3D) computational fluid dynamics (CFD) simulation shows uniform velocity and pressure distribution profiles in microchannels. For a channel velocity from 5.1 to 57.3 m/s in the combustor, the oxidation of methane is complete and self-sustainable without explosion, blow-off, or extinction; nevertheless, flashbacks are observed in counter-current mode. In the reformer, the maximum methane conversion is 84.9% in cocurrent mode, slightly higher than that of 80.2% in counter-current mode at a residence time of 33 ms, but at the cost of three times higher energy input in the combustor operating at ~1000 °C. Nitric oxide (NO) is not identified in combustion products, but nitrous oxide (N 2 O) is a function of coupling mode and forms significantly in cocurrent mode. Furthermore, this research would be helpful to establish the start-up strategy and environmental impact of compact reformers on a small scale.

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

Quantification of the carbon bonding state in amorphous carbon materials: A comparison between EELS and NEXAFS measurements

In this report the quantitative determination of the carbon hybridization is critical for establishing processing-structure-properties relationships for carbon-based materials, including amorphous carbon coatings. While several techniques have been employed to characterize the amount of sp 2 and sp 3 carbon in these materials, direct comparisons between analytical results are limited. Here, we compare near edge X-ray absorption fine structure (NEXAFS) spectra of a silicon- and oxygen-containing hydrogenated amorphous carbon (a-C:H:Si:O) coating acquired in synchrotron-based scanning transmission X-ray microscopy (STXM) mode with electron energy loss spectra (EELS) obtained from the same a-C:H:Si:O lamella. While the fractions of sp 2 carbon computed from STXM and EELS spectra are in close agreement, the comparison of NEXAFS spectra acquired in STXM mode with NEXAFS spectra collected in partial electron yield mode on a flat a-C:H:Si:O surface indicated that the destructive preparation of thin lamellae for STXM analyses induces variations in the structure of a-C:H:Si:O, namely the breakage of carbon-silicon and carbon-hydrogen bonds, a change in ordering of sp 2 -bonded carbon, and an increase in the sp 2 carbon fraction. These findings can help scientists in the careful interpretation of spectroscopic results obtained from the analysis of samples made of metastable materials after the destructive preparation of specimens for analytical purposes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shape restricted additive hazards models: Monotone, unimodal, and U‐shape hazard functions

We consider estimation of the semiparametric additive hazards model with an unspecified baseline hazard function where the effect of a continuous covariate has a specific shape but otherwise unspecified. Such estimation is particularly useful for a unimodal hazard function, where the hazard is monotone increasing and monotone decreasing with an unknown mode. A popular approach of the proportional hazards model is limited in such setting due to the complicated structure of the partial likelihood. Our model defines a quadratic loss function, and its simple structure allows a global Hessian matrix that does not involve parameters. Thus, once the global Hessian matrix is computed, a standard quadratic programming method can be applicable by profiling all possible locations of the mode. However, the quadratic programming method may be inefficient to handle a large global Hessian matrix in the profiling algorithm due to a large dimensionality, where the dimension of the global Hessian matrix and number of hypothetical modes are the same order as the sample size. We propose the quadratic pool adjacent violators algorithm to reduce computational costs. The proposed algorithm is extended to the model with a time‐dependent covariate with monotone or U‐shape hazard function. In simulation studies, our proposed method improves computational speed compared to the quadratic programming method, with bias and mean square error reductions. We analyze data from a recent cardiovascular study.

Mathematical & Computational Biology↗
DOE OSTI2025

Solution of the linear wave-particle kinetic equation for global modes of arbitrary frequency in a tokamak

The linear response of a plasma to perturbations of arbitrary frequency and wavelength is derived for any axisymmetric magnetized toroidal plasma. An explicit transformation to action-angle coordinates is achieved using orthogonal magnetic coordinates and the Littlejohn Lagrangian, establishing the validity of this result to arbitrary order in normalized Larmor radius. The global resonance condition for compressional modes is clarified in more detail than in previous works, confirming that the poloidal orbit-average of the cyclotron frequency gives the desired result at lowest order in Larmor radius. The global plasma response to the perturbation at each resonance is captured by a poloidal and gyroaverage of the perturbing potential. A “global gyroaveraging” of the potential is a natural by-product of this analysis which takes into account the changing of the magnetic field over an orbit. The resonance condition depends on two arbitrary integers which completely separately capture the effects poloidal non-uniformity and finite Larmor radius in generating sidebands. We learn that poloidal sidebands generated for compressional modes are dominated by the change in gyrofrequency over the orbit, which is very different to shear modes where the gyrofrequency only contributes via a finite Larmor radius effect. This increases the number of bounce harmonics required to compute the linear drive, giving a more complicated resonance map. An example calculation is given comparing resonance of shear and compressional modes in a published DIII-D case.

Compressional↗
DOE OSTI2025

Hybrid quantum simulations with qubits and qumodes on trapped-ion platforms

We explore the feasibility of gate-based hybrid quantum computing using both discrete (qubit) and continuous (qumode) variables on trapped-ion platforms. Trapped-ion systems have demonstrated record one- and two-qubit gate fidelities and long qubit coherence times, while qumodes, which can be represented by the collective vibrational modes of the ion chain, have remained relatively unex- plored for their use in computing. Using numerical simulations, we show that high-fidelity hybrid gates and measurement operations can be achieved for existing trapped-ion quantum platforms. As an exemplary application, we consider quantum simulations of the Jaynes-Cummings-Hubbard model, which is given by a one-dimensional chain of interacting spin and boson degrees of free- dom. Using classical simulations, we study its real-time evolution and develop a suitable variational quantum algorithm for ground state preparation. Furthermore, our results motivate further studies of hybrid quantum computing in this context, which may lead to direct applications in condensed matter and fundamental particle and nuclear physics.

Lower-dimensional field theories↗
DOE OSTI2020

Unitary Subharmonic Response and Floquet Majorana Modes

Detection and manipulation of excitations with non-Abelian statistics, such as Majorana fermions, are essential for creating topological quantum computers. To this end, we show the connection between the existence of such localized particles and the phenomenon of unitary subharmonic response (SR) in periodically driven systems. In particular, starting from highly nonequilibrium initial states, the unpaired Majorana modes exhibit spin oscillations with twice the driving period, are localized, and can have exponentially long lifetimes in clean systems. While the lifetime of SR is limited in translationally invariant systems, we show in this work that disorder can be engineered to stabilize the subharmonic response of Majorana modes. A viable observation of this phenomenon can be achieved using modern multiqubit hardware, such as superconducting circuits and cold atomic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mesoscopic Modeling and Rapid Simulation of Incremental Changes in Epidemic Scenarios on GPUs

In simulation-based studies and analyses of epidemics, a major challenge lies in resolving the conflict between fidelity of models and the speed of their simulation. Another related challenge arises in dealing with the large number of what–if scenarios that need to be explored. Here, we describe new computational methods that together provide an approach to dealing with both challenges. A mesoscopic modeling approach is described that strikes a middle ground between macroscopic models based on coupled differential equations and microscopic models built on fine-grained behaviors at the individual entity level. The mesoscopic approach offers the ability to incorporate complex compositions of multiple layers of dynamics even while retaining the potential for aggregate behaviors at varying levels. It also is an excellent match to the accelerator-based architectures of modern computing platforms in which graphical processing units (GPUs) can be exploited for fast simulation via the parallel execution mode of single instruction multiple thread (SIMT). The challenge of simulating a large number of scenarios is addressed via a method of sharing model state and computation across a tree of what–if scenarios that are localized, incremental changes to a large base simulation. A combination of the mesoscopic modeling approach and the incremental what–if scenario tree evaluation has been implemented in the software on modern GPUs. Synthetic simulation scenarios are presented to demonstrate the computational characteristics of our approach. Results from the experiments with large population data, including USA, UK, and India, illustrate the modeling methodology and computational performance on thousands of synthetically generated what–if scenarios. Execution of our implementation scaled to 8192 GPUs of supercomputing platforms demonstrates the ability to rapidly evaluate what–if scenarios several orders of magnitude faster than the conventional methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Improving Estimation of the Koopman Operator with Kolmogorov–Smirnov Indicator Functions

It has become common to perform kinetic analysis using approximate Koopman operators that transform high-dimensional timeseries of observables into ranked dynamical modes. The key to the practical success of the approach is the identification of a set of observables that form a good basis on which to expand the slow relaxation modes. Good observables are, however, difficult to identify a priori and suboptimal choices can lead to significant underestimations of characteristic time scales. Leveraging the representation of slow dynamics in terms of Hidden Markov Models (HMM), we propose a simple and computationally efficient clustering procedure to infer surrogate observables that form a good basis for slow modes. Here, we apply the approach to an analytically solvable model system as well as on three protein systems of different complexities. We consistently demonstrate that the inferred indicator functions can significantly improve the estimation of the leading eigenvalues of Koopman operators and correctly identify key states and transition time scales of stochastic systems, even when good observables are not known a priori.

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

Improving ideal MHD equilibrium accuracy with physics-informed neural networks

We present a novel approach to compute three-dimensional magnetohydrodynamic equilibria with isotropic pressure profiles and nested surfaces by parametrizing Fourier modes with artificial neural networks (NNs). The full nonlinear global force residual of single equilibria across the volume in real space is then minimized with first order optimizers and compared to equilibria computed by conventional solvers. Already, we observe competitive computational cost to arrive at the same minimum residuals computable with existing codes. With increased computational cost, lower minima of the residual are computable with the NNs than with any other tested solver, establishing a new lower bound for the force residual. We use minimally complex NNs, and we expect significant improvements for solving not only single equilibria with NNs, but also for creating NN models valid over continuous distributions of equilibria.

ideal magnetohydrodynamics↗
DOE OSTI2021

Redshift space power spectrum beyond Einstein-de Sitter kernels

Abstract We develop a framework to compute the redshift space power spectrum (PS), with kernels beyond Einstein-de Sitter (EdS), that can be applied to a wide variety of generalized cosmologies. We build upon a formalism that was recently employed for standard cosmology in Chen, Vlah & White (2020), and utilize an expansion of the density-weighted velocity moment generating function that explicitly separates the magnitude of the k-modes and their angle to the line-of-sight direction dependencies. We compute the PS for matter and biased tracers to 1-loop Perturbation Theory (PT) and show that the expansion has a correct infrared and ultraviolet behavior, free of unwanted divergences. We also add Effective Field Theory (EFT) counterterms, necessary to account for small-scale contributions to PT, and employ an IR-resummation prescription to properly model the smearing of the BAO due to large scale bulk flows within Standard-PT. To demonstrate the applicability of our formalism, we apply it on the ΛCDM and the Hu-Sawicki f(R) models, and compare our numerical results against the elephant suite of N-body simulations, finding very good agreement up to k = 0.27 Mpc -1 h at z = 0.5 for the first three non-vanishing Legendre multipoles of the PS. To our knowledge, the model presented in this work is the most accurate theoretical EFT-PT for modified gravity to date, being the only one that accounts for beyond linear local biasing in redshift-space. Hence, we argue our RSD modeling is a promising tool to construct theoretical templates in order to test deviations from ΛCDM using real data obtained from the next stage of cosmological surveys such as DESI and LSST.

Astronomy & Astrophysics↗
DOE OSTI2022

Topological surface superconductivity in FeSe 0.45 Te 0.55

The engineering of Majorana zero modes in topological superconductors, a paradigm for the realization of topological quantum computing and topology-based devices, has been hampered by the absence of materials with sufficiently large superconducting gaps. Recent experiments, however, have provided enthralling evidence for the existence of topological surface superconductivity in the iron-based superconductor FeSe 0.45 Te 0.55 possessing a full s> ± -wave gap of a few meV. Here, we propose a mechanism for the emergence of topological superconductivity on the surface of FeSe 0.45 Te 0.55 by demonstrating that the interplay between the s ± -wave symmetry of the superconducting gap, surface magnetism, and a Rashba spin–orbit interaction gives rise to robust topological superconducting phases. Moreover, the proposed mechanism explains a series of experimentally observed hallmarks of topological superconductivity, such as the emergence of Majorana zero modes in the center of vortex cores and at the end of line defects, as well as of chiral Majorana edge modes along domain walls. We also propose that the spatial distribution of supercurrents near a domain wall is a characteristic signature measurable via a scanning superconducting quantum interference device that can distinguish between chiral Majorana edge modes and trivial in-gap states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Are Marine Low Cloud Droplet Concentrations Buffered by Entrained Aitken‐Mode Aerosol (Final Technical Report)

During the summertime, the high-latitude oceans come to life with green phytoplankton, which gain their energy from sunlight and are food for sea creatures small and large. When the phytoplankton are eaten or die, sulfur-rich gases are released under the ocean surface and mix into the air. Observations suggest that, over the Southern Ocean, frequent storms lift this air high into the atmosphere while raining out particulates like salt. As a result, the sulfur-rich air then spawns high concentrations of small ‘Aitken-mode’ aerosol particles. We hypothesize that these particles work their way down into the marine boundary layer, where they can replenish the supply of cloud-condensation nuclei scavenged by frequent precipitation. Further, this process maintains high concentrations of liquid cloud droplets in austral summer, promoting more sunlight to be reflected to space. We call this ‘Aitken buffering’. The primary objective of this project has been to document and test what role Aitken-mode aerosols play in clouds over the Southern Ocean and elsewhere. This effort has included three main components: 1) developing a computer model that realistically simulates the aerosol processes and the small-scale turbulent air motions that move aerosols around and create the clouds, 2) using that model to interpret and extend these observations for process understanding, by allowing different factors that contribute to the aerosol budget, such as surface wind speed, precipitation, surface gas exchange, etc. to be separated, and 3) studying Aitken-mode aerosol and its variability with a focus over the Southern Ocean and Antarctica, using data from Atmospheric Radiation Measurement (ARM) sites and other available observations. Initial computer studies in more idealized conditions found that elevated concentrations of Aitken-mode aerosols above the clouds could help prevent the breakup of those clouds by acting as cloud-condensation nuclei after they were entrained into the cloudy boundary layer. Simulations of a day during the ACE-ENA field campaign showed that Aitken-mode aerosols could also prevent cloud breakup under more realistic conditions. A new method that extracts information about Aitken-mode aerosols from measurements of aerosols onto which cloud droplets can form finds that Aitken-mode aerosols do vary seasonally over the Southern Ocean, with a peak in summertime, as described above. Other work during this project has focused on understanding how patterns of water vapor, clouds and precipitation are coupled within low-lying clouds over the oceans, and also on how cloud droplets cluster within clouds and how the distribution of cloud droplet sizes change as dry air is mixed into clouds, with the latter studies also using observations from ACE-ENA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Photodetachment spectroscopy and resonant photoelectron imaging of the 2-naphthoxide anion via dipole-bound excited states

We report a photodetachment spectroscopy and high-resolution resonant photoelectron imaging study of cryogenically-cooled 2-naphthoxide anions (C 10 H 7 O - ). The photodetachment spectrum revealed a dipole-bound state (DBS) 202(4) cm -1 below the detachment threshold and thirty-eight resonances corresponding to the vibrational levels of the DBS. By tuning the detachment laser to these above-threshold resonances, we obtained thirty-eight resonantly-enhanced photoelectron spectra, which were highly non-Franck-Condon as a result of mode-selective vibrational auto-detachment from the DBS. The resonances were assigned by comparing the resonant and non-resonant photoelectron spectra, assisted by the computed vibrational frequencies. Specifically, vibrational features with low Franck-Condon factors or from Franck-Condon-forbidden vibrational modes were significantly enhanced in the resonant photoelectron spectra, resulting in much richer spectroscopic information. The electron affinity of the 2-naphthoxy radical was measured accurately to be 19,387(4) cm -1 or 2.4037(5) eV. In addition, a total of seventeen vibrational frequencies were obtained for the 2-naphthoxy radical. In particular, seven Franck-Condon-forbidden out-of-plane bending modes, including the two lowest frequency modes (ν 48 at 102 cm -1 and ν47 at 171 cm -1 ), were observed, demonstrating the advantages of combining photodetachment spectroscopy and resonant photoelectron spectroscopy in obtaining vibrational information for polar radical species via DBS.

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

Open quantum systems and Schwinger-Keldysh holograms

We initiate the study of open quantum field theories using holographic methods. Specifically, we consider a quantum field theory (the system) coupled to a holographic field theory at finite temperature (the environment). We investigate the effects of integrating out the holographic environment with an aim of obtaining an effective dynamics for the resulting open quantum field theory. The influence functionals which enter this open effective action are determined by the real-time (Schwinger-Keldysh) correlation functions of the holographic thermal environment. To evaluate the latter, we exploit recent developments, wherein the semiclassical gravitational Schwinger-Keldysh saddle geometries were identified as complexified black hole spacetimes. We compute real-time correlation functions using holographic methods in these geometries, and argue that they lead to a sensible open effective quantum dynamics for the system in question, a question that hitherto had been left unanswered. In addition to shedding light on open quantum systems coupled to strongly correlated thermal environments, our results also provide a principled computation of Schwinger-Keldysh observables in gravity and holography. In particular, these influence functionals we compute capture both the dissipative physics of black hole quasi- normal modes, as well as that of the fluctuations encoded in outgoing Hawking quanta, and interactions between them. We obtain results for these observables at leading order in a low frequency and momentum expansion in general dimensions, in addition to determining explicit results for two dimensional holographic CFT environments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗