Search NASA⌕ Search

SEARCH · Search NASA

Results for “time evolution operators”

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 487 records · Page 27

NBI optimization on SMART and implications for scenario development

Abstract The SMall Aspect Ratio Tokamak (SMART) under commissioning at the University of Seville, Spain, aims to explore confinement properties and possible advantages in confinement for compact/spherical tokamaks operating at negative vs. positive triangularity. This work explores the benefits of auxiliary heating through Neutral Beam Injection (NBI) for SMART scenarios beyond the initial Ohmic phase of operations, in support of the device’s mission. Expected values of electron and ion temperature achievable with NBI heating are first predicted for the current flat-top phase, including modeling to optimize the NBI injection geometry to maximize NBI absorption and minimize losses for a given equilibrium. Simulations are then extended for a selected case to cover the current ramp-up phase. Differences with results obtained for the flat-top phase indicate the importance of determining the plasma evolution over time, as well as self-consistently determining the edge plasma parameters for reliable time-dependent simulations. Initial simulation results indicate the advantage of auxiliary NBI heating to achieve nearly double values of pressure and stored energy compared to Ohmic discharges, thus significantly increasing the device’s performance. The scenarios developed in this work will also contribute to diagnostic development and optimization for SMART, as well as providing test cases for initial predictions of macro- and micro-instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutronics analyses for the conceptual design of the SNS Second Target Station

The Spallation Neutron Source, in operation at the Oak Ridge National Laboratory since 2006, was designed to allow the addition of a second target station and an upgrade of the accelerator proton power. Both upgrades are now underway. Here we describe the evolution of the design of the target of the second target station with the emphasis on the effects of the proton beam footprint on the energy deposition in the target, stresses induced by the pulsed operation, and the importance of the residual heat. The moderator configurations and their optimization are discussed. With the utilization of pure parahydrogen moderators, small neutron beam cross-sections, and the specific optimization, neutron beams of the second target station will achieve exceptionally high peak brightness and time-averaged brightness.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling the evolution of the US opioid crisis for national policy development

The opioid crisis remains one of the greatest public health challenges in the United States. The crisis is complex, with long delays and feedbacks between policy actions and their effects, which creates a risk of unintended consequences and complicates policy decision-making. We present SOURCE (Simulation of Opioid Use, Response, Consequences, and Effects), an operationally detailed national-level model of the opioid crisis, intended to enhance understanding of the crisis and guide policy decisions. Drawing on multiple data sources, SOURCE replicates how risks of opioid misuse initiation and overdose have evolved over time in response to behavioral and other changes and suggests how those risks may evolve in the future, providing a basis for projecting and analyzing potential policy impacts and solutions.

59 BASIC BIOLOGICAL SCIENCES↗

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Multi-modal Approach to Understanding Degradation of Organic Photovoltaic Materials

State-of-the-art organic photovoltaic (OPV) materials are composed of complex, chemically diverse polymeric and molecular structures that form highly intricate solid-state interactions, collectively yielding exceptional tunability in performance and aesthetics. These properties are especially attractive for semitransparent power-generating windows or shades in living environments, greenhouses, or other architectural integrations. However, before such a future is realized, a broader and deeper understanding of property stability must be acquired. Stability during operating and environmental conditions is critical, namely, material color steadfastness, optoelectronic performance retention, morphological rigidity, and chemical robustness. To date, no single investigation encompasses all four distinct, yet interconnected, metrics. Here, we present a multimodal strategy that captures a dynamic and interconnected evolution of each property during the course of an accelerated photobleaching experiment. We demonstrate this approach across relevant length scales (from molecular to visual macroscale) using X-ray photoelectron spectroscopy, grazing-incidence X-ray scattering, microwave conductivity, and time-dependent photobleaching spectroscopies for two high-performance semitransparent OPV blends—PDPP4T:PC 60 BM and PDPP4T:IEICO-4F, with comparisons to the stabilities of the individual components. We present direct evidence that specific molecular acceptor (fullerene vs nonfullerene) designs and the resulting donor–acceptor interactions lead to distinctly different mechanistic routes that ultimately arrive at what is termed “OPV degradation.” We directly observe a chemical oxidation of the cyano endcaps of the IEICO-4F that coincides with a morphological change and large loss in photoconductivity while the fullerene acceptor-containing blend demonstrates a significantly greater fraction of oxygen uptake but retains 55% of the photoconductivity. Finally, this experimental roadmap provides meaningful guidance for future high-throughput, multimodal studies, benchmarking the sensitivity of the different analytical techniques for assessing stability in printable active layers, independent of complete device architectures.

14 SOLAR ENERGY↗

Semiclassical initial value representation for the quantum propagator in the Heisenberg interaction representation

One of the challenges facing on-the-fly ab initio semiclassical time evolution is the large expense needed to converge the computation. In this paper, we suggest that a significant saving in computational effort may be achieved by employing a semiclassical initial value representation (SCIVR) of the quantum propagator based on the Heisenberg interaction representation. We formulate and test numerically a modification and simplification of the previous semiclassical interaction representation of Shao and Makri [J. Chem. Phys. 113, 3681 (2000)]. The formulation is based on the wavefunction form of the semiclassical propagation instead of the operator form, and so is simpler and cheaper to implement. The semiclassical interaction representation has the advantage that the phase and prefactor vary relatively slowly as compared to the “standard” SCIVR methods. This improves its convergence properties significantly. Using a one-dimensional model system, the approximation is compared with Herman-Kluk’s frozen Gaussian and Heller’s thawed Gaussian approximations. The convergence properties of the interaction representation approach are shown to be favorable and indicate that the interaction representation is a viable way of incorporating on-the-fly force field information within a semiclassical framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the stability and stationarity of the Super H-mode combined with an ion transport barrier in the core

The Super H-mode, a regime with high pedestal pressure and stored energy, is explored on DIII-D and combined with an ion transport barrier in the plasma core to increase performance. A significant improvement of ion temperatures and confinement is facilitated by favorable conditions such as high rotational shear and high ion pedestal temperatures. As a result of a rise in density and simultaneous decrease in rotation, the ion transport barrier disappears during the discharge evolution, leading to a transition from a very high confinement state at early times, to a reduced but still high confinement phase. Additionally, in many discharges, a global MHD event consistent with the coupling of a destabilized internal mode to an edge localized mode (ELM) causes a large energy loss and leads to a reorganization of the plasma into a lower temperature, higher density state. Depending on the magnitude of the global MHD event, the plasma edge collisionality can increase significantly and shift the operational boundary from the peeling to the ballooning side, which can be understood as a drop out of the Super H-mode channel into standard H-mode. Hence, in Super H-mode discharges with ion transport barriers, both the improved pedestal height and rotational shear contribute to the high stored energy. Here, at very low levels of rotation, the confinement factor for SH modes is still expected to exceed standard H-mode by 20-30 %. With their overall stationarity and high-performance levels, Super H-mode discharges provide an attractive regime for ITER and may enable a more compact design of future fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microstructure Investigations of U{sub 3}Si{sub 2} Irradiated by Heavy Ions at LWR Temperatures

Searching for a nuclear fuel material with enhanced accident tolerance compared to the conventional UO{sub 2} fuel is one of the key objectives of the accident tolerant fuel (ATF) campaign. Having both high thermal conductivity and high uranium load, U{sub 3}Si{sub 2} is expected to provide lower operating temperature/stored energy and advantageous neutron economy as a light water reactor (LWR) fuel material. Therefore, U{sub 3}Si{sub 2} has been regarded as a promising candidate fuel material by the ATF community. To evaluate the qualification of U{sub 3}Si{sub 2} as an ATF fuel in LWRs, its fuel performance during both steady-state operation and power transients must be well understood. As a typical nuclear fuel material applied in research reactors, the fuel performance of U{sub 3}Si{sub 2} at low temperatures (<300 deg. C) has been extensively studied through both experimental and computational approaches. However, under LWR conditions, the temperature of U{sub 3}Si{sub 2} can go beyond 600 deg. C. In that case, the references collected in U{sub 3}Si{sub 2}'s applications in research reactors are inadequate for the validation of U{sub 3}Si{sub 2} as an LWR ATF. Experimental efforts that are dedicated to the LWR conditions must be made to help examine U{sub 3}Si{sub 2}'s qualifications as well as provide valuable references for the development related models in advanced fuel performance code. Gaseous fission products accumulate in nuclear fuel during operation and form fission gas bubbles. The evolution of those fission gas bubbles results in increase of volume (swelling), degradation of thermal conductivity, and eventually fission gas release, compromising the integrity and performance of fuel elements in multiple aspects. Hence, fission gas behavior is one of the most important components of fuel performance. In-pile irradiation and corresponding post-irradiation examinations (PIEs) undoubtedly provide most reliable information on the actual performance of a nuclear fuel. However, the high cost and long lead time involved in in-pile irradiation limit its wide applications. On the other hand, heavy ion irradiation is capable of simulating the neutron effects on materials, especially the behavior of energetic fission fragments in nuclear fuel materials, to a considerable extent. Therefore, heavy ion irradiation was used in this study to investigate the microstructural evolution of U{sub 3}Si{sub 2} under LWR-temperature irradiation, focusing on the difference from low temperature irradiation and fission gas behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Carbon Abundance Inhomogeneities and Deep Mixing Rates in Galactic Globular Clusters

Among stars in Galactic globular clusters the carbon abundance tends to decrease with increasing luminosity on the upper red giant branch, particularly within the lowest metallicity clusters. While such a phenomena is not predicted by canonical models of stellar interiors and evolution, it is widely held to be the result of some extra mixing operating during red giant branch ascent which transports material exposed to the CN(O)-cycle across the radiative zone in the stellar interior and into the base of the convective envelope, whereupon it is brought rapidly to the stellar surface. Here we present measurements of [C/Fe] abundances among 67 red giants in 19 globular clusters within the Milky Way. Building on the work of Martell et al., we have concentrated on giants with absolute magnitudes of M {sub V} ∼ −1.5 within clusters encompassing a range of metallicity (−2.4 < [Fe/H] < −0.3). The Kitt Peak National Observatory (KPNO) 4 m and Southern Astrophysical Research (SOAR) 4.1 m telescopes were used to obtain spectra covering the λ4300 CH and λ3883 CN bands. The CH absorption features in these spectra have been analyzed via synthetic spectra in order to obtain [C/Fe] abundances. These abundances and the luminosities of the observed stars were used to infer the rate at which C abundances change with time during upper red giant branch evolution (i.e., the mixing efficiency). By establishing rates over a range of metallicity, the dependence of deep mixing on metallicity is explored. We find that the inferred carbon depletion rate decreases as a function of metallicity, although our results are dependent on the initial [C/Fe] composition assumed for each star.

79 ASTRONOMY AND ASTROPHYSICS↗

Stress evolution and creep deformation in solid-oxide electrolysis cell systems – Dynamic modeling and multi-objective optimization to maximize stack life and efficiency

Here, this study develops a thermal stress model of solid-oxide electrolysis cells (SOECs) including a model for creep strain and failure probability that is integrated with a dynamic plant-wide model of a hydrogen production process. Uncertainties in key material properties of the cell are quantified to assess their impact on stress profile variability. The oxygen electrode is found to have about 10 times higher failure probability compared to the fuel electrode. The study shows that if the stack operation is not optimized, cycling operation would lead to stress build-up eventually leading to catastrophic failure. A dynamic optimization problem is set up for obtaining the optimal operational profile considering a variable hydrogen production rate. Due to the tradeoff between the efficiency and stress build-up, the dynamic optimization problem is multi-objective. It is observed that the optimizer can considerably reduce the stress build-up (i.e., can increase the stack life) albeit at the cost of a lower efficiency thus exhibiting strong tradeoffs between capital and operating costs. For example, if the stack would be replaced in 0.5 yr, specific energy requirement would be 48.5 kWh/kg H 2 while for a stack replacement time of about 6 yr, the specific energy requirement rises by about 4.2 %.

SOEC↗

Study of MHD Stability, Active Mode Control, and Disruption Avoidance in NSTX-U

Learning how to achieve sustained operation at high plasma pressure is an important goal of the National Spherical Torus Upgrade (NSTX-U) program and is critically important to any future fusion power plant. This research addresses high priority stability and control issues outlined by the NSTX-U Program. Reduced models of the stabilization physics learned continue to be created and tested against actual plasma operation in NSTX, NSTX-U, and other tokamak devices. This research is forming physics-based models for disruption forecasting that aims to unify results across different fusion machines for more confident extrapolation of the techniques studied to future devices including ITER. The Disruption Event Characterization and Forecasting (DECAF) code was expanded to analyze data from the TCV tokamak, including real-time MHD analysis signals. The presence of a mode and likelihood of its m/n numbers can be tracked, as well as the mode frequency, although a limitation at low frequency does not yet allow tracking modes to the locking point in TCV. Initial analyses of NSTX-U operation were conducted using the DECAF code. The timing of sawtooth oscillations was investigated taking advantage of the more extensive L-mode operation of the device in its first run campaign. An initial investigation of thermal collapses in NSTX-U was also conducted with specific comparison to the Greenwald density limit. A local island power balance model to explain the density limit in tokamaks has been tested with NSTX data and is being evaluated as a potential disruption forecaster in the DECAF code, presently showing similarities to the global Greenwald density evolution and limit. A module was written for the DECAF code that generalizes magnetic spectrogram analysis of MHD modes. This new portable code, written in Python, discriminates the toroidal mode number of instabilities and can track each mode as a DECAF event, monitoring the evolution of the mode through bifurcation and locking points. This is a critical new capability to allow disruption event characterization for general tokamak databases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reactor simulations with nuclear data uncertainties

This work demonstrates the influence of uncertainties in microscopic nuclear data on the results of reactor simulations, both for stationary and transient states. It gives an overview on the methods in use for uncertainty and sensitivity analyses with respect to nuclear data, and discusses the pros and cons. For full-scale reactor simulations, in particular with coupled neutron transport and thermo-hydraulics, random sampling provides a powerful means to propagate nuclear data uncertainties through the complete calculation sequence. Results of uncertainty analyses performed with the GRS XSUSA – “cross section (XS) Uncertainty and Sensitivity Analysis” methodology are shown for radial power distributions from steady-state PWR calculations and for the time evolution of the reactor power in the course of a control rod withdrawal from a PWR mini-core. In all cases, the output uncertainties are considerable. For the radial power distributions, relative 1σ uncertainties of up to 10% are observed, and for the power peak during the transient, the relative 1σ uncertainty reaches 20%. These large uncertainties strongly suggest to routinely accompany best-estimate simulations by uncertainty analyses with respect to nuclear data, in particular for systems beyond LWR for which much less operation experience is available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multimodal Analysis of Reaction Pathways of Cathode Materials for Lithium Ion Batteries

Conversion mechanism in lithium ion batteries provides higher capacity than intercalation mechanism since multiple numbers of electrons and lithium ions are associated. However, poor cycling stability, large voltage hysteresis, and low energy efficiency have been great challenges of conversion reaction. To address those issues, conversion-type electrode materials have been reformed via doping or substituting other elements. For example, iron fluorides (FeF 2 , FeF 3 ) have modified as iron oxyfluorides (FeF 1-x O x ), showing enhanced long-term stability. Furthermore, co-substituted (both anion and cation substituted) Fe 0.9 Co 0.1 OF (FeCoOF) was demonstrated excellent cycling stability. (capacity of 350 mAh g -1 at a current of 500 mA g -1 for 1000 cycles). Substituting anion and cation in iron fluoride has been suggested as an effective method to achieve better reversibility but understanding of lithiation reactions in co-substituted FeCoOF is not clear. This work takes advantage of ex-situ/ in-situ synchrotron X-ray based techniques and transmission electron microscopy to elucidate structural changes with lithium ion, which may provide fundamental insights into modifying conversion-type materials. Figure 1 presents discharge-charge curves and pair distribution function patterns acquired at each potential. As lithium ions were inserted, structural changes were noticed both at short-range and long-range. However, long-range ordering was nearly maintained even at 1 V, indicating absence of conversion reaction. Figure 2 shows lithiation induced structural evolution of a single FeCoOF nanorod observed in real time. As lithiation proceeds, the width of nanorods shows a stepwise increase, particularly A in figure 2d, which may indicate multiple steps of lithiation occur. Considering that conversion reaction takes place around 2 V in FeF 3 , co-substitution Co and O into iron fluoride may change thermodynamic features of lithiation reactions by lowering the initiation potential for conversion reaction. Instead, phase transformations occur at long-range order, which may help maintaining structural integrity during operation, eventually, achieving cycling stability.

25 ENERGY STORAGE↗

Calculating Time Eigenvalues of the Neutron Transport Equation with Dynamic Mode Decomposition

A novel method to compute time eigenvalues of neutron transport problems is presented based on solutions to the time-dependent transport equation. Using these solutions, we use the dynamic mode decomposition to form an approximate transport operator. This approximate operator has eigenvalues that are mathematically related to the time eigenvalues of the neutron transport equation. This approach works for systems of any level of criticality and does not require the user to have estimates for the eigenvalues. Numerical results are presented for homogeneous and heterogeneous media. The numerical results indicate that the method finds the eigenvalues that contribute the most to the change in the solution over a given time range, and the eigenvalue with the largest real part is not necessarily important to the system evolution at short and intermediate times.

42 ENGINEERING↗

Triviality of quantum trajectories close to a directed percolation transition

Here we study quantum circuits consisting of unitary gates, projective measurements, and control operations that steer the system toward a pure absorbing state. Two types of phase transition occur as the rate of these control operations is increased: a measurement-induced entanglement transition, and a directed percolation transition into the absorbing state (taken here to be a product state). In this work, we show analytically that these transitions are generically distinct, with the quantum trajectories becoming disentangled before the absorbing state transition is reached, and we analyze their critical properties. We introduce a simple class of models where the measurements in each quantum trajectory define an effective tensor network (ETN)—a subgraph of the initial spacetime graph where nontrivial time evolution takes place. By analyzing the entanglement properties of the ETN, we show that the entanglement and absorbing-state transitions coincide only in the limit of the infinite local Hilbert-space dimension. Focusing on a Clifford model which allows numerical simulations for large system sizes, we verify our predictions and study the finite-size crossover between the two transitions at large local Hilbert space dimension. We give evidence that the entanglement transition is governed by the same fixed point as in hybrid circuits without feedback.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Clean Grid Vision: A U.S. Perspective - Chapter 5. Global Power Market Trends

The evolution of the electricity sector toward competitive modern markets is occurring over a long period and responding to rapid changes in technologies, business models, and other market factors. No reformed power market anywhere in the world looks the same today as it did a decade ago. Technology has enabled some of the evolution, but perhaps it is more important that people and institutions learn from experience how to make markets better over time. Governments still have an important role—not as the authority controlling prices, but as the authority safeguarding the integrity of the economic processes that create prices in response to supply and demand and ensure reliability and other objectives are met. More VRE resources will change market conditions, but they will not change market fundamentals. Prices provide important signals for new investments and how the grid operates—including new technologies and processes discussed elsewhere in this report. If poor market design or the exercise of market power by a large entity distort these signals, outcomes will cease to be efficient and the public will bear the additional cost. Fundamentally, prices reflect supply, demand, and the cost of technological options.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Development of High Capacity Energy Storage Materials

Hydrogen fuel cells have the potential to offer energy and power density advantages over lithium ion batteries in automotive and portable power applications when paired with an appropriate hydrogen storage system. Development of the ideal hydrogen storage material has been immensely sought after but plagued by limitations present in each type of material. These limitations typically include nonidealistic operational temperatures, low capacities, excessive costs, lack of reversibility, or evolution of impurities which irreversibly damage fuel cell performance. Many complex metal hydrides possess suitable hydrogen capacities but unfortunately suffer from impurity release. In order to address this problem, our team has explored in-situ catalyst doping of alkali metal amides to demonstrate methodologies to reduce ammonia release. The incorporation of iridium metal in lithium amide significantly reduces the ammonia release and decomposes the ammonia to provide additional hydrogen. Furthermore, we demonstrate, for the first time, high resistance of low temperature fuel cells to ammonia in comparison with typical proton exchange membrane fuel cells. The utilization of these two approaches in tandem provides a novel pathway for the development and implementation of high capacity energy storage materials for fuel cell applications.

25 ENERGY STORAGE↗

Harnessing Virtual Power Plants Reliably: Enabling tools for increased observability, controllability, operation, and aggregation of distributed energy resources

Harnessing virtual power plants enhances the integration of distributed energy resources into utility grids for a sustainable energy future. Virtual power plants (VPPs) aggregate DERs to enhance resource adequacy and reduce emissions. U.S. utilities are exploring various technologies to manage DERs effectively. FERC Order 2222 allows DERs to participate in both wholesale and retail markets. Enhancing observability and controllability of behind-the-meter (BTM) DERs is essential for reliable grid operations. A hierarchical control architecture can improve coordination among residential energy resources. Field tests showed nearly 20% energy savings and 30% peak power reduction during grid events. Effective management of DERs requires enhanced situational awareness to prevent grid congestion. Integrating DER management systems (DERMS) with existing planning tools can improve operational security. Near-real-time grid models can validate optimal resource set points against resource uncertainty. Traditional uninterruptible power supplies (UPS) can be upgraded to support grid services and become part of VPPs. Upgrading UPS systems can reduce costs by 75% and unlock significant battery capacity. New battery management systems and grid-aware controllers are essential for optimizing UPS performance. Continued research and development are necessary to address challenges in integrating DERs into utility grids. Encouraging customer participation in pilot programs is vital for the evolution of VPPs. Here, the shift towards price-responsive DERs and VPPs is expected to enhance energy distribution efficiency.

24 POWER TRANSMISSION AND DISTRIBUTION↗