Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Periodic systems”

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

Multi-Period Optimization of Multi-Timescale Energy Systems: Application to Solid-Oxide Electrolysis Cells

A presentation detailing the recent progress in the optimization of Solid-Oxide Cells under Chemical Degradation. The existing methods typically involve a quasi-steady state assumption to make this large-scale problem tractable. Here, we introduce an extension to this method that allows for adaptive coupling of degradation with the fast-timescale process based on a variety of error thresholds. This reduces the error accumulated during long-term optimization of the SOC under degradation and allows for decision making at multiple timescales. The results presented include long-term operating profiles of the SOC under degradation with both steady-state and fluctuating day to day operation.

Giridhar, Nishant↗
DOE OSTI2020

Classical discrete time crystals

The spontaneous breaking of time-translation symmetry in periodically driven quantum systems leads to a new phase of matter: the discrete time crystal (DTC). This phase exhibits collective subharmonic oscillations that depend upon an interplay of non-equilibrium driving, many-body interactions and the breakdown of ergodicity. However, subharmonic responses are also a well-known feature of classical dynamical systems ranging from predator–prey models to Faraday waves and a.c.-driven charge density waves. This raises the question of whether these classical phenomena display the same rigidity characteristic of a quantum DTC. In this work, we explore this question in the context of periodically driven Hamiltonian dynamics coupled to a finite-temperature bath, which provides both friction and, crucially, noise. Focusing on one-dimensional chains, where in equilibrium any transition would be forbidden at finite temperature, we provide evidence that the combination of noise and interactions drives a sharp, first-order dynamical phase transition between a discrete time-translation invariant phase and an activated classical discrete time crystal (CDTC) in which time-translation symmetry is broken out to exponentially long timescales. Power-law correlations are present along a first-order line, which terminates at a critical point. We analyse the transition by mapping it to the locked-to-sliding transition of a d.c.-driven charge density wave. To conclude, building upon results from the field of probabilistic cellular automata, we conjecture the existence of classical time crystals with true long-range order, where time-translation symmetry is broken out to infinite times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dissipative Floquet Dynamics: from Steady State to Measurement Induced Criticality in Trapped-ion Chains

Quantum systems evolving unitarily and subject to quantum measurements exhibit various types of non-equilibrium phase transitions, arising from the competition between unitary evolution and measurements. Dissipative phase transitions in steady states of time-independent Liouvillians and measurement induced phase transitions at the level of quantum trajectories are two primary examples of such transitions. Investigating a many-body spin system subject to periodic resetting measurements, we argue that many-body dissipative Floquet dynamics provides a natural framework to analyze both types of transitions. We show that a dissipative phase transition between a ferromagnetic ordered phase and a paramagnetic disordered phase emerges for long-range systems as a function of measurement probabilities. A measurement induced transition of the entanglement entropy between volume law scaling and sub-volume law scaling is also present, and is distinct from the ordering transition. The two phases correspond to an error-correcting and a quantum-Zeno regimes, respectively. The ferromagnetic phase is lost for short range interactions, while the volume law phase of the entanglement is enhanced. An analysis of multifractal properties of wave function in Hilbert space provides a common perspective on both types of transitions in the system. Our findings are immediately relevant to trapped ion experiments, for which we detail a blueprint proposal based on currently available platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Automated, reliable, and efficient continental-scale replication of 7.3 petabytes of computational simulation data: A case study

We report on our experiences replicating 7.3 petabytes (PB) of Earth System Grid Federation (ESGF) computational simulation data from Lawrence Livermore National Laboratory (LLNL) in California to Argonne National Laboratory (ANL) in Illinois and Oak Ridge National Laboratory (ORNL) in Tennessee—a task motivated by a need for increased reliability, capacity, and performance. This task presented significant challenges: the need to move 29 million files twice under time pressure from aging storage hardware; a source file system bottleneck limiting throughput to 1.5 GB/s; frequent site maintenance windows; and the need for complete reliability at scale. We addressed these challenges using a simple replication tool that invoked Globus to transfer large bundles of files while tracking progress in a database, dynamically rerouting transfers to work around maintenance periods and file system limitations. Under the covers, Globus organized transfers to make efficient use of the high-speed Energy Sciences network (ESnet) and the data transfer nodes deployed at participating sites, and also addressed security, integrity checking, and recovery from a variety of transient failures. This success demonstrates the considerable benefits that can accrue from the adoption of performant data replication infrastructure. The replication tool is available at https://github.com/esgf2-us/data-replication-tools.

Globus↗
DOE OSTI2020

Thermal Insulation for a High Temperature Molten Salt Storage Tank in a CSP Plant

A ternary molten chloride has been suggested as a high temperature (550-720 °C) sensible heat storage medium for a two-tank system. An effective thermal insulation is proposed to provide an opportunity for the fabrication of both the hot and cold storage tanks from lower cost metals, e.g. A240-347H, by maintaining the tank wall below 500 °C. To achieve this, the hot storage tank mainly comprises of inner shell (11.4 cm of a refractory brick as liner (KX-99®), 30 cm of insulation brick) and outer metal shell, where 5 cm of frozen layer of PCM800 (NaCl), and 5 cm of ceramic wool fill the annular space between the shells. Another layer of ceramic wool was introduced at the outer side of the metal shell to keep the temperature lower than 60 °C. ANSYS Fluent was used for the numerical study of the cool down behavior of the system during idle period when tank is quarter-filled. The heat transfer modes were considered as radiation and natural convection heat transfer from the hot molten salt at 720 °C, conduction through the solid layers and convection heat transfer from the outer shell to the environment at 15 °C. The results showed that the proposed insulation system including a frozen layer of PCM800 maintains the temperature at the metal wall lower than 500 °C, and lower than 60°C at the outer side of the tank. Apart from the insulation property of the frozen layer of PCM800, it absorbs any leak of hot molten salt, acting as a barrier between salt and the metal shell. The thermal performance analysis provide an insight into the benefits of the multilayer insulation.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2020

Thermal Insulation for a High Temperature Molten Salt Storage Tank in a CSP Plant: Preprint

A ternary molten chloride has been suggested as a high temperature (550-720 °C) sensible heat storage medium for a two-tank system. An effective thermal insulation is proposed to provide an opportunity for the fabrication of both the hot and cold storage tanks from lower cost metals, e.g. A240-347H, by maintaining the tank wall below 500 °C. To achieve this, the hot storage tank mainly comprises of inner shell (11.4 cm of a refractory brick as liner (KX-99 ®), 30 cm of insulation brick) and outer metal shell, where 5 cm of frozen layer of PCM800 (NaCl), and 5 cm of ceramic wool fill the annular space between the shells. Another layer of ceramic wool was introduced at the outer side of the metal shell to keep the temperature lower than 60 °C. ANSYS Fluent was used for the numerical study of the cool down behavior of the system during idle period when tank is quarter-filled. The heat transfer modes were considered as radiation and natural convection heat transfer from the hot molten salt at 720 °C, conduction through the solid layers and convection heat transfer from the outer shell to the environment at 15 °C. The results showed that the proposed insulation system including a frozen layer of PCM800 maintains the temperature at the metal wall lower than 500 °C, and lower than 60 °C at the outer side of the tank. Apart from the insulation property of the frozen layer of PCM800, it absorbs any leak of hot molten salt, acting as a barrier between salt and the metal shell. The thermal performance analysis provide an insight into the benefits of the multilayer insulation.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

Renewable hydrogen and ammonia for combined heat and power systems in remote locations: Optimal design and scheduling

Abstract Using hydrogen (H ) and ammonia (NH ) for renewable energy storage has the potential to enable economical power and heat supply with high renewable penetrations, especially in remote locations which are characterized by high energy costs. In this work we assess the economic competitiveness of renewable combined heat and power (CHP) systems in Mahaka HI, Nantucket MA, and Northwest Arctic Borough (NWAB) AK by optimally designing these systems for scenarios in which power and heat can be purchased over a range of historical energy prices as well as when 100% renewable supply is required. We use a combined optimal design and scheduling model which minimizes annualized net present cost by determining optimal technology selection and size simultaneously with optimal schedules for each period of a system operating horizon aggregated from full year hourly resolution data via a consecutive temporal clustering algorithm. We find that renewable generation meets at least 85% of power demands and 75% of heat demands under the lowest energy prices investigated. Higher conventional energy prices lead to increased renewable penetration which is facilitated by renewable NH as a seasonal energy storage medium, as are 100% renewable CHP systems. NH is used for power generation with heat cogeneration in all three locations, as well as directly for heating in NWAB. On an annual cost basis, NH ‐enabled 100% renewable CHP is only 3% more expensive in Mahaka and NWAB than systems which can purchase energy at the lowest prices, while it is 15% more expensive in Nantucket.

Palys, Matthew J.↗
DOE OSTI2026

Quantum Computing Algorithms and Applications for Coherent and Strongly Correlated Chemical Systems

This project advanced quantum algorithms, quantum information theory, strongly correlated electronic structure methods, molecular quantum materials, and exciton transport imaging in coherent condensed phase systems. Across the award period, the team developed new methods for open-quantum-system simulation, Hamiltonian learning, state tomography, reduced-density-matrix and contracted-quantum eigensolver approaches, and quantum diagnostics for device capability and openness.

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

Modeling Flexible Generator Operating Regions via Chance- constrained Stochastic Unit Commitment

Here, we introduce a novel chance-constrained stochastic unit commitment model to address uncertainty in renewables' production uncertainty in power systems operation. For most thermal generators,underlying technical constraints that are universally treated as "hard" by deterministic unit commitment models are in fact based on engineering judgments, such that system operators can periodically request operation outside these limits in non-nominal situations, e.g., to ensure reliability. We incorporate this practical consideration into a chance-constrained stochastic unit commitment model, specifically by in-frequently allowing minor deviations from the minimum and maximum thermal generator power output levels. We demonstrate that an extensive form of our model is computationally tractable for medium-sized power systems given modest numbers of scenarios for renewables' production. We show that the model is able to potentially save significant annual production costs by allowing infrequent and controlled violation of the traditionally hard bounds imposed on thermal generator production limits. Finally, we conduct a sensitivity analysis of optimal solutions to our model under two restricted regimes and observe similar qualitative results.

42 ENGINEERING↗
DOE OSTI2022

Volume extrapolation via eigenvector continuation

Here, we develop an extension of eigenvector continuation (EC) that makes it possible to extrapolate simulations of quantum systems in finite periodic boxes across large ranges of box sizes. The formal justification for this approach, which we call finite-volume eigenvector continuation (FVEC), is provided by matching periodic functions at different box sizes. As concrete FVEC implementation we use a discrete variable representation based on plane-wave states and present several applications calculated within this framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Complex scaling in finite volume

Quantum resonances, i.e., metastable states with a finite lifetime, play an important role in nuclear physics and other domains. Describing this phenomenon theoretically is generally a challenging task. In this work, we combine two established techniques to address this challenge. Complex scaling makes it possible to calculate resonances with bound-state-like methods. Finite-volume simulations exploit the fact that the infinite-volume properties of quantum systems are encoded in how discrete energy levels change as one varies the size of the volume. Herein we apply complex scaling to systems in finite periodic boxes and derive the volume dependence of states in this scenario, demonstrating with explicit examples how one can use these relations to infer infinite-volume resonance energies and lifetimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Plutonium in Manhattan Project workers: Using autopsy data to evaluate organ content and dose estimates based on urine bioassay with implications for radiation epidemiology

Purpose Radiation dose estimates in epidemiology typically rely on intake predictions based on urine bioassay measurements. The purpose of this article is to compare the conventional dosimetric estimates for radiation epidemiology with the estimates based on additional post-mortem tissue radiochemical analysis results. Methods The comparison was performed on a unique group of 11 former Manhattan Project nuclear workers, who worked with plutonium in the 1940s, and voluntarily donated their bodies to the United States Transuranium and Uranium Registries. Results Post-mortem organ activities were predicted using different sets of urine data and compared to measured activities. Use of urinalysis data collected during the exposure periods overestimated the systemic (liver+skeleton) deposition of 239 Pu by 155±134%, while the average bias from using post-exposure urinalyses was –4±50%. Committed effective doses estimated using early urine data differed from the best estimate by, on average, 196±193%; inclusion of follow-up urine measurements in analyses decreased the mean bias to 0.6±36.3%. Cumulative absorbed doses for the liver, red marrow, bone surface, and brain were calculated for the actual commitment period. Conclusion On average, post-exposure urine bioassay results were in good agreement with post-mortem tissue analyses and were more reliable than results of urine bioassays collected during the exposure.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Shapiro steps and stability of skyrmions interacting with alternating anisotropy under the influence of ac and dc drives

Here we use atomistic simulations to examine the sliding dynamics of a skyrmion in a two-dimensional system containing a periodic one-dimensional stripe pattern of variations between low and high values of the perpendicular magnetic anisotropy. The skyrmion changes in size as it crosses the interface between two anisotropy regions. On applying combined dc and ac driving in either parallel or perpendicular directions, we observe a wide variety of Shapiro steps, Shapiro spikes, and phase-locking phenomena. The phase-locked orbits have two-dimensional dynamics due to the gyrotropic or Magnus dynamics of the skyrmions and are distinct from the phase-locked orbits found for strictly overdamped systems. Along a given Shapiro step when the ac drive is perpendicular to the dc drive, the velocity parallel to the ac drive is locked while the velocity in the perpendicular direction increases with increasing drive to form Shapiro spikes. At the transition between adjacent Shapiro steps, the parallel velocity jumps up to the next step value, and the perpendicular velocity drops. The skyrmion Hall angle shows a series of spikes as a function of increasing dc drive, where the jumps correspond to the transition between different phase-locked steps. At high drives, the Shapiro steps and Shapiro spikes are lost. When both the ac and dc drives are parallel to the stripe periodicity direction, Shapiro steps appear, while if the dc drive is parallel to the stripe periodicity direction and the ac drive is perpendicular to the stripe periodicity, then there are only two locked phases, and the skyrmion motion consists of a combination of sliding along the interfaces between the two anisotropy values and jumping across the interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory

Here, we investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet–triplet excitation energy is 2.31 eV and that for the optically inactive singlet–singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV – defect using NEVPT2 for periodic supercell models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5

Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.

CATHARE↗
DOE OSTI2023

Drude weights in one-dimensional systems with a single defect

Ballistic transport of a quantum system can be characterized by Drude weight, which quantifies the response of the system to a uniform electric field in the infinitely long timescale. The Drude weight is often discussed in terms of the Kohn formula, which gives the Drude weight by the derivative of the energy eigenvalue of a finite-size system with the periodic boundary condition in terms of the Aharonov-Bohm flux. Recently, the Kohn formula is generalized to nonlinear responses. However, the nonlinear Drude weight determined by the Kohn formula often diverges in the thermodynamic limit. In order to elucidate the issue, in this work we examine a simple example of a one-dimensional tight-binding model in the presence of a single defect at zero temperature. We find that its linear and nonlinear Drude weights given by the Kohn formula (i) depend on the Aharonov-Bohm flux and (ii) diverge proportionally to a power of the system size. Here, we argue that the problem can be attributed to different order of limits. The Drude weight according to the Kohn formula (“Kohn-Drude weight”) indicates the response of a finite-size system to an adiabatic insertion of the Aharonov-Bohm flux. While it is a well-defined physical quantity for a finite-size system, its thermodynamic limit does not always describe the ballistic transport of the bulk. The latter should be rather characterized by a “bulk Drude weight” defined by taking the thermodynamic limit first before the zero-frequency limit. While the potential issue of the order of limits has been sometimes discussed within the linear response, the discrepancy between the two limits is amplified in nonlinear Drude weights. We demonstrate the importance of the low-energy excitations of O(1/L), which are excluded from the Kohn-Drude weight, in regularizing the bulk Drude weight.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Temporal entanglement transitions in the periodically driven Ising chain

Periodically driven quantum systems can host non-equilibrium phenomena without static analogs, including in their entanglement dynamics. Here, we discover (TET) in a Floquet spin chain, which correspond to a quantum phase transition in the spectrum of the entanglement Hamiltonian and are signaled by dynamical spontaneous symmetry breaking. We identify the symmetry principles underlying these transitions: they appear when the driven Hamiltonian preserves global symmetry (here, ℤ2), the initial state respects this symmetry, and the reduced density matrix carries weight in both subsystem-parity sectors, with TET occurring precisely when the sector weights become equal (given the previous two conditions are also satisfied). Intriguingly, we find these transitions across a broad range of driving frequencies (from adiabatic to high-frequency regime) and independently of drive details, where they manifest as periodic, sharp entanglement spectrum reorganizations marked by the Schmidt-gap closure, a vanishing entanglement echo, and symmetry-quantum-number flips, while remaining invisible to conventional local observables. At high frequencies, the entanglement Hamiltonian acquires an intrinsic timescale decoupled from the drive period, rendering the transitions genuine steady-state features. Finite-size scaling reveals universal critical behavior with correlation-length exponent 𝜈 = 1, matching equilibrium Ising universality despite its emergence from purely dynamical mechanisms decoupled from static criticality. Furthermore, our work establishes TET as novel features in Floquet quantum matter.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evolving Metrics for Resource Adequacy Assessment

Resource adequacy analysis quantifies the likelihood of capacity shortfall on a power system in a probabilistic manner. Using a combination of statistical techniques and power system fundamentals, the analysis typically evaluates hundreds or thousands of stochastic random samples (replications) of varying load, generator outages, variable renewable energy availability, and other aspects of power system uncertainty. In this range of uncertainty, there are - at times - periods where the power system's available resources are insufficient to meet system demand, referred to as a shortfall event. Today's power systems' rapidly evolving generation mix is changing the types of data needed by system planners and regulators, which can often render traditional resource adequacy metrics insufficient for ensuring resource adequacy for tomorrow's grid. In this paper we provide a critical assessment of traditional measures of shortfall risk in power systems, discussing their shortcomings and how they compare to metrics used in other domains. From this analysis we propose four steps forward for improving power system resource adequacy risk metrics in the future.

ENERGY PLANNING, POLICY, AND ECONOMY,POWER TRANSMI↗