Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “exponential time integration”

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 OSTI2023

Generation of Near-Equipartition Magnetic Fields in Turbulent Collisionless Plasmas

The mechanisms that generate “seed” magnetic fields in our Universe and that amplify them throughout cosmic time remain poorly understood. By means of fully kinetic particle-in-cell simulations of turbulent, initially unmagnetized plasmas, here we study the genesis of magnetic fields via the Weibel instability and follow their dynamo growth up to near-equipartition levels. In the kinematic stage of the dynamo, we find that the rms magnetic field strength grows exponentially with rate γ B ≃0.4⁢u rms /L, where L/2⁢π is the driving scale and u rms is the rms turbulent velocity. In the saturated stage, the magnetic field energy reaches about half of the turbulent kinetic energy. Here, magnetic field growth is balanced by dissipation via reconnection, as revealed by the appearance of plasmoid chains. At saturation, the integral-scale wave number of the magnetic spectrum approaches k int ≃12⁢π/L. Our results show that turbulence—induced by, e.g., the gravitational buildup of galaxies and galaxy clusters—can magnetize collisionless plasmas with large-scale near-equipartition fields.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Modeling of Turbulence Generated Noise in Jets

A numerically calculated Green's function is used to predict jet noise spectrum and its far-field directivity. A linearized form of Lilley's equation governs the non-causal Green s function of interest, with the non-linear terms on the right hand side identified as the source. In this paper, contributions from the so-called self- and shear-noise source terms will be discussed. A Reynolds-averaged Navier-Stokes solution yields the required mean flow as well as time- and length scales of a noise-generating turbulent eddy. A non-compact source, with exponential temporal and spatial functions, is used to describe the turbulence velocity correlation tensors. It is shown that while an exact non-causal Green's function accurately predicts the observed shift in the location of the spectrum peak with angle as well as the angularity of sound at moderate Mach numbers, at high subsonic and supersonic acoustic Mach numbers the polar directivity of radiated sound is not entirely captured by this Green's function. Results presented for Mach 0.5 and 0.9 isothermal jets, as well as a Mach 0.8 hot jet conclude that near the peak radiation angle a different source/Green's function convolution integral may be required in order to capture the peak observed directivity of jet noise.

Khavaran, Abbas↗
NASA NTRS2012

SGR J1550-5418 Bursts Detected with the Fermi Gamma-Ray Burst Monitor during Its Most Prolific Activity

We have performed detailed temporal and time-integrated spectral analysis of 286 bursts from SGR J1550-5418 detected with the Fermi Gamma-ray Burst Monitor (GBM) in 2009 January, resulting in the largest uniform sample of temporal and spectral properties of SGR J1550-5418 bursts. We have used the combination of broadband and high time-resolution data provided with GBM to perform statistical studies for the source properties.We determine the durations, emission times, duty cycles, and rise times for all bursts, and find that they are typical of SGR bursts. We explore various models in our spectral analysis, and conclude that the spectra of SGR J15505418 bursts in the 8-200 keV band are equally well described by optically thin thermal bremsstrahlung (OTTB), a power law (PL) with an exponential cutoff (Comptonized model), and two blackbody (BB) functions (BB+BB). In the spectral fits with the Comptonized model, we find a mean PL index of -0.92, close to the OTTB index of -1. We show that there is an anti-correlation between the Comptonized E(sub peak) and the burst fluence and average flux. For the BB+BBfits, we find that the fluences and emission areas of the two BB functions are correlated. The low-temperature BB has an emission area comparable to the neutron star surface area, independent of the temperature, while the high temperature BB has a much smaller area and shows an anti-correlation between emission area and temperature.We compare the properties of these bursts with bursts observed from other SGR sources during extreme activations, and discuss the implications of our results in the context of magnetar burst models.

bursts↗
DOE OSTI2022

EVSE Cybersecurity and Resilience

Consequence-driven Cybersecurity Analysis for Extreme Fast Charging Electric Vehicle Infrastructure Electric vehicle (EV) development and associated charging infrastructure are expected to advance rapidly. Thirty percent of all global vehicle sales may be EVs and hybrid EVs by 2025, and they will rely on increasingly sophisticated strategies for grid integration. Next-generation EV charging infrastructure is expected to include interconnected renewable resources, such as photovoltaic (PV) arrays and battery storage systems, along with grid-edge devices. Although distributed energy resources (DERs) are useful in several ways, such as peak shaving at high demand times and backup supply for added resilience, the integration of vehicle charging and DERs could create more avenues for cyberattack. Physical and/or remote access to EV charging station components, including charge ports, power electronics, controllers, and local generation (e.g., PV and energy storage) could be paths to cause power fluctuations, leading to altered operations at the charging station, escalated privileges to administrative systems, exfiltration of financial information (including personally identifiable information), and reduced grid stability. One compromised EV supply equipment component can open the door to a variety of exploitable vulnerabilities. Cloud computing and mobile application control have the potential to expand the threat surface to non-repudiation and firmware integrity challenges. Vendor clouds have access to hundreds of chargers, and if compromised, can scale the attack surface exponentially. The high power and voltage levels of xFC infrastructure (e.g., 400 kW at 1000- V DC) increase the hazards and ability to impact the grid and vehicles more than lower-power charging systems. Legacy communications systems and protocols could also put EV infrastructure at risk of cyberattacks requiring a robust patch management process. Communications networks link EVs and chargers to several stakeholders - including charging station operators, grid operators, vendors/manufacturers, and aggregators - who have both physical and network access to share information for control, monitoring, and analytics. Information in these networks that is vulnerable to compromise includes the state of charge, charging duration, payment information, electricity price, and load control. Analyzing and prioritizing these interconnections risks could help address cybersecurity related to data leakage and manipulation.

charging↗
NASA NTRS1991

Asymptotic integration algorithms for nonhomogeneous, nonlinear, first order, ordinary differential equations

New methods for integrating systems of stiff, nonlinear, first order, ordinary differential equations are developed by casting the differential equations into integral form. Nonlinear recursive relations are obtained that allow the solution to a system of equations at time t plus delta t to be obtained in terms of the solution at time t in explicit and implicit forms. Examples of accuracy obtained with the new technique are given by considering systems of nonlinear, first order equations which arise in the study of unified models of viscoplastic behaviors, the spread of the AIDS virus, and predator-prey populations. In general, the new implicit algorithm is unconditionally stable, and has a Jacobian of smaller dimension than that which is acquired by current implicit methods, such as the Euler backward difference algorithm; yet, it gives superior accuracy. The asymptotic explicit and implicit algorithms are suitable for solutions that are of the growing and decaying exponential kinds, respectively, whilst the implicit Euler-Maclaurin algorithm is superior when the solution oscillates, i.e., when there are regions in which both growing and decaying exponential solutions exist.

Walker, K. P.↗
DOE OSTI2025

Anatomy of information scrambling and decoherence in the integrable Sachdev-Ye-Kitaev model

The growth of information scrambling, captured by out-of-time-order correlation functions (OTOCs), is a central indicator of the nature of many-body quantum dynamics. Here, we compute analytically the complete time dependence of the OTOC for an integrable Sachdev-Ye-Kitaev (SYK) model, 𝑁 Majoranas with random two-body interactions of infinite range, coupled to a Markovian bath at finite temperature. In the limit of no coupling to the bath, the time evolution of scrambling experiences different stages. For 𝑡≲√𝑁, after an initial polynomial growth, the OTOC approaches saturation in a power-law fashion with oscillations superimposed. At 𝑡∼√𝑁, the OTOC reverses trend and starts to decrease linearly in time. The reason for this linear decrease is that, despite being a subleading 1/𝑁 effect, the OTOC in this region is governed by the spectral form factor of the antisymmetric couplings of the SYK model. The linear decrease stops at 𝑡∼2⁢𝑁, the Heisenberg time, where saturation occurs. The effect of the environment is an overall exponential decay of the OTOC for times longer than the inverse of the coupling strength to the bath. Here, the oscillations at 𝑡≲√𝑁 indicate lack of thermalization—a desired feature for better performance of quantum information devices.

Open quantum systems & decoherence↗
DOE OSTI2021

Dimensionality reduction of the many-body problem using coupled-cluster subsystem flow equations: classical and quantum computing perspective

We discuss reduced-scaling strategies employing recently introduced sub-system embedding sub-algebras coupled-cluster formalism (SES-CC) to describe many-body systems. These strategies utilize properties of the SES-CC formulations where the equations describing certain classes of sub- systems can be integrated into a computational flows composed coupled eigenvalue problems of reduced dimensionality. Additionally, these flows can be defined at the level of the CC Ansatz defined by selected classes of cluster amplitudes, which define the wave function ”memory” of possible partitionings of the many-body system into constituent sub-systems. One of the possible ways of solving these coupled problems is through implementing procedures, where the information is passed between the sub-systems in a self-consistent manner. As a special case, we consider local flow formulations where the so-called local character of correlation effects can be closely related to properties of sub-system embedding sub-algebras employing localized molecular basis. We also generalize flow equations to the time domain and to downfolding methods utilizing double exponential unitary CC Ansatz (DUCC), where reduced dimensionality of constituent sub-problems offer a possibility of efficient utilization of limited quantum resources in modeling realistic systems.

Electron correlation, quantum chemistry, quantum c↗
NASA NTRS1994

O(+)-O collision frequency in high-speed flows

Throughout much of the terrestrial thermosphere and ionosphere, the motions of the neutral and ionized constituents are closely coupled and relative velocities are small, of the order of 100 m/s or less. This is particularly true at midlatitudes to low latitudes where typical velocities in the neutral gas due to tidal forcing are only 20-50 m/s. However, the solar wind-magnetosphere interaction drives a large-scale convection pattern in the polar ionosphere. When the rapid adjustment of the plasma to changes in the solar wind is combined with the slower response of the more massive neutral gas, large relative velocities on the order of 1 km/s can exist for substantial lengths of time. This will be more common during periods of high geomagnetic activity, as a result of the greater number of magnetic substorms and other particle precipitation events. When a significant relative velocity is present, the calculation of interaction parameters of the two gases passing through each other, such as collision frequency, must include that velocity. These effects are usually neglected when interpreting wind and ion drift observations. We show how the collision frequency is affected by a directed velocity between any two gases interacting with a power law or exponential potential energy curved. The directed velocity increases the collision frequency at all temperatures for most ion-neutral interactions. For certain power law potentials, such as the charge quadrupole, the collision frequency is decreased. We present an analytic solution for the high-speed collision integral using the resonance charge exchange cross section.

Pesnell, William Dean↗
NASA NTRS2013

Tuning and Robustness Analysis for the Orion Absolute Navigation System

The Orion Multi-Purpose Crew Vehicle (MPCV) is currently under development as NASA's next-generation spacecraft for exploration missions beyond Low Earth Orbit. The MPCV is set to perform an orbital test flight, termed Exploration Flight Test 1 (EFT-1), some time in late 2014. The navigation system for the Orion spacecraft is being designed in a Multi-Organizational Design Environment (MODE) team including contractor and NASA personnel. The system uses an Extended Kalman Filter to process measurements and determine the state. The design of the navigation system has undergone several iterations and modifications since its inception, and continues as a work-in-progress. This paper seeks to show the efforts made to-date in tuning the filter for the EFT-1 mission and instilling appropriate robustness into the system to meet the requirements of manned space ight. Filter performance is affected by many factors: data rates, sensor measurement errors, tuning, and others. This paper focuses mainly on the error characterization and tuning portion. Traditional efforts at tuning a navigation filter have centered around the observation/measurement noise and Gaussian process noise of the Extended Kalman Filter. While the Orion MODE team must certainly address those factors, the team is also looking at residual edit thresholds and measurement underweighting as tuning tools. Tuning analysis is presented with open loop Monte-Carlo simulation results showing statistical errors bounded by the 3-sigma filter uncertainty covariance. The Orion filter design uses 24 Exponentially Correlated Random Variable (ECRV) parameters to estimate the accel/gyro misalignment and nonorthogonality. By design, the time constant and noise terms of these ECRV parameters were set to manufacturer specifications and not used as tuning parameters. They are included in the filter as a more analytically correct method of modeling uncertainties than ad-hoc tuning of the process noise. Tuning is explored for the powered-flight ascent phase, where measurements are scarce and unmodelled vehicle accelerations dominate. On orbit, there are important trade-off cases between process and measurement noise. On entry, there are considerations about trading performance accuracy for robustness. Process Noise is divided into powered flight and coasting ight and can be adjusted for each phase and mode of the Orion EFT-1 mission. Measurement noise is used for the integrated velocity measurements during pad alignment. It is also used for Global Positioning System (GPS) pseudorange and delta- range measurements during the rest of the flight. The robustness effort has been focused on maintaining filter convergence and performance in the presence of unmodeled error sources. These include unmodeled forces on the vehicle and uncorrected errors on the sensor measurements. Orion uses a single-frequency, non-keyed GPS receiver, so the effects due to signal distortion in Earth's ionosphere and troposphere are present in the raw measurements. Results are presented showing the efforts to compensate for these errors as well as characterize the residual effect for measurement noise tuning. Another robustness tool in use is tuning the residual edit thresholds. The trade-off between noise tuning and edit thresholds is explored in the context of robustness to errors in dynamics models and sensor measurements. Measurement underweighting is also presented as a method of additional robustness when processing highly accurate measurements in the presence of large filter uncertainties.

Holt, Greg N.↗
DOE OSTI2023

Robust highly durable solid oxide fuel cell cathodes – Improved materials compatibility & self-regulating surface chemistry

Solid oxide fuel cells (SOFCs) are electrochemical conversion devices that directly transform hydrogen or hydrocarbon fuels to electricity, with energy efficiencies as high as 90%, coupled with reduced emissions. Several factors, however, remain to be addressed when considering scale-up of SOFC technology, including the need to overcome decreased performance due to sluggish rates of the oxygen reduction reaction (ORR) at the cathode under reduced temperatures and susceptibility to degradation in performance from surface poisoning e.g. from chromia, while limiting the use of critical raw materials (lanthanides and transition metals) present in high performing mixed ionic electronic conducting electrodes like (La,Sr)CoO 3 (LSC). In this project we explored the key descriptors for determining ORR activity in SOFC electrodes and tried to recover performance degradation by applying them to SOFC electrodes. In order to do this, we first selected a model mixed ionic electronic conducting (MIEC) oxide, Pr-doped CeO 2 (Pr 0.1 Ce 0.9 O 2-δ , PCO), which is a chemically stable fluorite and free of inherent poison sources (e.g. Sr segregation in LSC) that potentially react with external impurities such as Cr-species vaporized from the interconnect. The three approaches originally planned in this project are as follows: 1) evaluation of scavenger exsolution characteristics, 2) study of scavengers gettering efficacy following Cr and Si poisoning and 3) integration of new compositions into porous electrodes. Among them, exceptional progress has been made in 2) and 3), especially understanding the role of surface infiltrants in impacting electrode performance and degradation of PCO materials. We found that the Smith acidity scale for binary oxides serves as a powerful descriptor for tuning and predicting the oxygen exchange kinetics on MIEC PCO surfaces. As a result, with infiltration with binary oxides, ranging from strongly basic (Li 2 O) to strongly acidic (SiO 2 ) onto the surface of porous PCO, it was possible to systematically vary the chemical surface exchange coefficient (k chem ) by 6 orders of magnitude! L i2 O increased k chem by nearly 1,000 times over that of pristine PCO, while SiO 2 decreased k chem by nearly the same factor. Strikingly, although the pre-exponential of k chem scales linearly with the acidity of the infiltrated binary oxide, there is nearly no change in the activation energy. With this insight, we attributed the origin of these dramatic changes in k chem values to the systematic increase and decrease in the surface electron density induced by infiltrated binary oxides. More interestingly, although both Cr 2 O 3 and SiO 2 were determined to be acidic by Smith, suggesting that this feature could likely be the primary reason that these compounds serve to poison the ORR on SOFC cathodes, the effect of poisoning could be subsequently tuned by adding multiple infiltrants and controlling their relative surface acidities. We also systematically examined the effect of serial infiltration of basic and acidic oxides. It turned out that serial infiltration of Li not only recovers approximately 20-fold degraded k chem of PCO by acidic Cr 2 O 3 but its k chem is enhanced even beyond that of the non-infiltrated PCO by more than three orders of magnitude. This was further verified with a screen-printing PCO symmetric cell in terms of the electrode performance (area-specific resistance, ASR) related to approach 3). These observations point to acidity as a key descriptor not only in tuning and predicting the ORR activity of SOFC cathodes that largely determines the overall performance of SOFC, but in mitigating and reactivating poisoned electrode performance. This work provides novel guidelines for making the electrode performance much more active and robust in SOFCs, which can further be applied to all applications requiring oxygen exchange reaction, such as electrolyzers, permeation membranes and gas sensors.

30 DIRECT ENERGY CONVERSION↗
NASA NTRSDate not supplied

Developing Mars-Based Clinical Scenarios for an Earth Independent Medical Operations (EIMO) – Based Decision Support Service

As crewed missions move beyond Low-Earth Orbit, pre-mission planning cannot fully buy down the medical risks of exploration-class missions. Martian missions, where increased hazards exist, (such as long-duration spaceflight, surface-level EVA operations, and communications delays) will require a paradigm shift in the structure of a medical system. An Earth-Independent Medical Operations-based Medical System (EIMO-MS) will need to optimize four critical domains to help provide medical care: utilization of Pre-Mission Planning, augmentation of Acute and Prolonged Medical Decision Making, automated tracking of Resource Management, and assistance in Task Load Balance. The ideal EIMO-MS will be able to accomplish this goal by having an interactive, adaptable interface that will be able to provide real-time medical services. It must respond based on the level of crewmember training, medical situation, and available medical and non-medical resources. To showcase the capabilities and requirements of such a sophisticated automated MS, a series of clinical scenarios of escalating complexity were developed with clinical and systems engineering input. These scenarios describe in clinical detail what a theoretical future medical system, enhanced with multiple information streams (such as a medical database, an AI-based Decision Support System, real-time monitoring, enhanced in-situ laboratory imaging, etc.) can achieve in conjunction with a trained and experienced crew. Scenarios are comprised of: a context section including objectives and applicable spaceflight environment, a highlighted assumptions section, a clinical narrative section, and a systems engineering activity diagram demonstrating the integrated Medical System (MS). The “swim lanes” of the activity diagram act as the logistical core of each scenario and show how the MS will interact with the crew, ground support, and other in-flight systems. The Design Reference Mission that is used for the scenarios is based on existing reference mission profiles [1] with a projected 30-sol stay on the Martian surface. Scenarios span the spectrum from planned evaluations, minor medical care, urgent care, surgical guidance, critical and expectant management, and behavioral health care. Mission complexity will exponentially increase during deep space and Mars exploration-class missions, and medical support for these missions will likewise need to increase in autonomy and adaptability. The integrated system that will support these missions will need to provide assistance in a variety of anticipated and unforeseen scenarios. These medical scenarios, guided by clinician input, are initial steps in crafting the requirements for an EIMO-based medical system. By working in a systems engineering framework, requirements and capabilities can be extracted and mapped while maintaining a clinical core.

Prashant Parmar↗
DOE OSTI2024

Future of plasma etching for microelectronics: Challenges and opportunities

Plasma etching is an essential semiconductor manufacturing technology required to enable the current microelectronics industry. Along with lithographic patterning, thin-film formation methods, and others, plasma etching has dynamically evolved to meet the exponentially growing demands of the microelectronics industry that enables modern society. At this time, plasma etching faces a period of unprecedented changes owing to numerous factors, including aggressive transition to three-dimensional (3D) device architectures, process precision approaching atomic-scale critical dimensions, introduction of new materials, fundamental silicon device limits, and parallel evolution of post-CMOS approaches. The vast growth of the microelectronics industry has emphasized its role in addressing major societal challenges, including questions on the sustainability of the associated energy use, semiconductor manufacturing related emissions of greenhouse gases, and others. The goal of this article is to help both define the challenges for plasma etching and point out effective plasma etching technology options that may play essential roles in defining microelectronics manufacturing in the future. The challenges are accompanied by significant new opportunities, including integrating experiments with various computational approaches such as machine learning/artificial intelligence and progress in computational approaches, including the realization of digital twins of physical etch chambers through hybrid/coupled models. These prospects can enable innovative solutions to problems that were not available during the past 50 years of plasma etch development in the microelectronics industry. To elaborate on these perspectives, the present article brings together the views of various experts on the different topics that will shape plasma etching for microelectronics manufacturing of the future.

Engineering↗
NASA NTRS1977

Control logic to track the outputs of a command generator or randomly forced target

A procedure is presented for synthesizing time-invariant control logic to cause the outputs of a linear plant to track the outputs of an unforced (or randomly forced) linear dynamic system. The control logic uses feed-forward of the reference system state variables and feedback of the plant state variables. The feed-forward gains are obtained from the solution of a linear algebraic matrix equation of the Liapunov type. The feedback gains are the usual regulator gains, determined to stabilize (or augment the stability of) the plant, possibly including integral control. The method is applied here to the design of control logic for a second-order servomechanism to follow a linearly increasing (ramp) signal, an unstable third-order system with two controls to track two separate ramp signals, and a sixth-order system with two controls to track a constant signal and an exponentially decreasing signal (aircraft landing-flare or glide-slope-capture with constant velocity).

Trankle, T. L.↗
NASA NTRS1999

Control Design Strategies to Enhance Long-Term Aircraft Structural Integrity

Over the operational lifetime of both military and civil aircraft, structural components are exposed to hundreds of thousands of low-stress repetitive load cycles and less frequent but higher-stress transient loads originating from maneuvering flight and atmospheric gusts. Micro-material imperfections in the structure, such as cracks and debonded laminates, expand and grow in this environment, reducing the structural integrity and shortening the life of the airframe. Extreme costs associated with refurbishment of critical load-bearing structural components in a large fleet, or altogether reinventoring the fleet with newer models, indicate alternative solutions for life extension of the airframe structure are highly desirable. Increased levels of operational safety and reliability are also important factors influencing the desirability of such solutions. One area having significant potential for impacting crack growth/fatigue damage reduction and structural life extension is flight control. To modify the airframe response dynamics arising from command inputs and gust disturbances, feedback loops are routinely applied to vehicles. A dexterous flight control system architecture senses key vehicle motions and generates critical forces/moments at multiple points distributed throughout the airframe to elicit the desired motion characteristics. In principle, these same control loops can be utilized to influence the level of exposure to harmful loads during flight on structural components. Project objectives are to investigate and/or assess the leverage control has on reducing fatigue damage and enhancing long-term structural integrity, without degrading attitude control and trajectory guidance performance levels. In particular, efforts have focused on the effects inner loop control parameters and architectures have on fatigue damage rate. To complete this research, an actively controlled flexible aircraft model and a new state space modeling procedure for crack growth have been utilized. Analysis of the analytical state space model for crack growth revealed the critical mathematical factors, and hence the physical mechanism they represent, that influenced high rates of airframe crack growth. The crack model was then exercised with simple load inputs to uncover and expose key crack growth behavior. To characterize crack growth behavior, both "short-term" laboratory specimen test type inputs and "long-term" operational flight type inputs were considered. Harmonic loading with a single overload revealed typical exponential crack growth behavior until the overload application, after which time the crack growth was retarded for a period of time depending on the overload strength. An optimum overload strength was identified which leads to maximum retardation of crack growth. Harmonic loading with a repeated overload of varying strength and frequency again revealed an optimum overload trait for maximizing growth retardation. The optimum overload strength ratio lies near the range of 2 to 3 with dependency on frequency. Experimental data was found to correlate well with the analytical predictions.

Newman, Brett A.↗
DOE OSTI2020

Integration of Total-Sky Imager Data with a Physics-Based Smart Persistence Model for Intra-Hour Forecasting of Solar Radiation

Short-term solar forecasting models based solely on global horizontal irradiance (GHI) measurements are often unable to discriminate the forecasting of the factors affecting GHI from those that can be precisely computed by atmospheric models. Our previous study introduced a Physics-based Smart Persistence model for Intra-hour forecasting of solar radiation (PSPI) that decomposed the forecasting of GHI into the computation of extraterrestrial solar radiation and solar zenith angle and the forecasting of cloud albedo and cloud fraction. The extraterrestrial solar radiation and solar zenith angle were accurately computed by the Solar Position Algorithm (SPA) developed at the National Renewable Energy Laboratory (NREL). A cloud retrieval technique was used to estimate cloud albedo and cloud fraction from surface-based observations of GHI. With the assumption of persistent cloud structures, the cloud albedo and cloud fraction were predicted for future time steps using a two-stream approximation and a 5-minute exponential weighted moving average, respectively. The model evaluation indicated the estimation and forecast of cloud fraction mostly contributed to the uncertainty of the PSPI though it overcame the persistence and smart persistence models in all forecast time horizons between 5 and 60 minutes. This study aims to enhance the PSPI by ingesting surface-based observations of cloud fraction from a total sky imager (TSI). The estimation and forecast of cloud albedo is correspondingly improved by utilizing the cloud fraction observations and thus leads to more accurate GHI forecast. Various time-series analysis methods are also investigated on the forecasting of cloud fraction and cloud albedo for further improving the GHI forecast. These improvements are valuable for many applications, such as forecasting energy use for buildings, grid operations, and ultimately bringing down the cost of solar energy.

14 SOLAR ENERGY↗
DOE OSTI2021

Linear stability analysis of two fluid columns of different densities and viscosities in a gravity field

The linear stability of a vertical interface separating two miscible fluid columns of different densities and viscosities under the influence of gravity is investigated. This flow possesses a time-dependent reference state (each column accelerates at different rates owing to their different densities) and the interface thickness grows as the square root of time (by diffusion). Numerical integration of the linear initial-value problem is carried out and discussed in detail as a function of vertical and spanwise wavenumbers and the flow parameters. Adjoint-based optimization is performed in order to determine initial conditions that lead to maximum growth of disturbances in finite time. Results indicate that the rate of growth of the perturbation energy at small wavenumbers (less affected by viscosity initially) is dominated by two-dimensional modes (no spanwise variation). Substantial transient growth is observed at higher wave modes initially, followed by asymptotic decay of the perturbations at large time. Sensitivity of perturbation growth with respect to initial time, density and viscosity ratios is investigated. This work is complementary to previous inviscid analysis of this configuration, which showed that the interface was unconditionally unstable at all wave modes, even in the presence of surface tension, and that instability grew as the exponential of time squared.

Mechanics↗
DOE OSTI2022

Combined Cycle Integrated Thermal Energy Storage “CiTES” (Final Scientific/Technical Report)

The Phase I of this project confirmed the technical feasibility of a Combined Cycle integrated Thermal Energy Storage “CiTES” system, calculated the key performance parameters like power efficiency and costs, and proved its commercial value with full-year simulations for several US electricity markets with high degree of variable renewable generation and volatile hourly electricity prices. The core element of this project is the Electro Thermal Energy Storage (ETES) technology from Siemens Gamesa Renewable Energy GmbH, using thermally stable and inexpensive volcanic rocks as storage material and air as heat transfer medium. This technology is backed by more than 10 years of experience and a 440MMBTU (130MWh-th) pilot plant in Hamburg, Germany, which is in operation since 2019. The integration of this thermal storage in an existing combined cycle power plant (CCPP) is typical power plant technology without any major technology risks. It allows the storage of inexpensive renewable energy during times of surplus renewable generation and the discharge of this energy in times of high energy demand when the fossil plant is in operation. This supplements the fossil power generation with CO2-emission-free energy. The secondary effect of the CiTES system is that a small part of the stored thermal energy is used to keep the heat recovery steam generator (HRSG) and steam turbine (ST) of the combined cycle power plant in hot and ready-to-start condition. This enables the plant to start rapidly when fossil generation is required to satisfy demand as soon variable generation drops off in the evenings or during cloud cover and calm wind periods. Without pre-warming of the HRSG and ST, the CCPP would need several hours for a cold or warm start, burn a lot of gas and release high NOx emissions during start and wouldn’t be able to use the short times of high energy prices in an efficient or economical manner. The economic parameters of CiTES were determined by a full year “8760” simulation using a data set calculation for each of the hours of the year, and historical electricity and gas prices. For consistency, the simulations were focused on the pre-COVID year 2019. The financially most attractive markets were in the Energy Reliability Council of Texas (ERCOT) region, which allowed substantial value generation with arbitrage (charge with cheap energy during renewable surplus times and discharge when energy is needed and expensive). The improvement of flexibility with the CiTES system by pre-heating and warm-keeping of the CCPP allowed for additional power generation during short time periods when demand is high but renewable generation is down; when the hourly energy prices are highest in these markets. The simulations are based on 2019 data, when ERCOT had 27GW of installed photovoltaic (PV) and wind generation. They showed that the created revenue with the prototypically sized CiTES system of 1,000MMBTU (300MWh-th) falls a little bit short of what is expected from a commercially viable investment. The system has optimization opportunities for cost reduction and increased effectiveness which will be realized during a potential Phase II Pre-FEED study following this project. Furthermore, it is safe to assume that a lot of renewable generation capacity will be added all over the US in the coming years. As an example, ERCOT is predicting to more than double its renewable generation from 27GW in 2019 to a forecasted 63GW in 2023. This will increase the amount of renewable overproduction exponentially. This rapid increase of local overproduction and the need to curtail renewable generation is well documented by the California ISO (www.CAISO.com / managing oversupply). However, the simulations also revealed a weakness in the structure of the electricity markets in the US. More specifically, when electricity prices are very low and approaching negative levels, the owners of Variable Renewable Energy (VRE) will curtail a part of their facility to stabilize the price by reducing supply. This results in a situation in which storage facilities, which are integrated in existing fossil assets and don’t have the behind-the-meter benefit of a VRE, won’t be able to purchase low cost – otherwise curtailed – renewable energy off the grid. A special tariff, which motivates VRE owners to sell otherwise curtailed renewable energy to storage facilities (Hydrogen, thermal, pumped hydro, etc.) can solve this issue. The implementation of such a regulating tariff by Independent System Operators, thus avoiding renewable curtailment, is a pre condition for successful commercialization for renewable energy storage technologies. With this advancement of design and technology and improvements in the market environment, it can be expected that the Combined Cycle integrated Thermal Energy Storage proves itself as an important innovation to keep highly efficient, natural gas-based power generation economically successful and relevant for the power industry in the United States of America.

Wolf, Thorsten↗
DOE OSTI2025

Nonperturbative quantum gravity in a closed Lorentzian universe

We study how meaningful physical predictions can arise in nonperturbative quantum gravity in a closed Lorentzian universe. In such settings, recent developments suggest that the quantum gravitational Hilbert space is one-dimensional and real for each α-sector, as induced by spacetime wormholes. This appears to obstruct the conventional quantum-mechanical prescription of assigning probabilities via projection onto a basis of states. While previous approaches have introduced external observers or augmented the theory to resolve this issue, we argue that quantum gravity itself contains all the necessary ingredients to make physical predictions. We demonstrate that the emergence of classical observables and probabilistic outcomes can be understood as a consequence of partial observability: physical observers access only a subsystem of the universe. Tracing out the inaccessible degrees of freedom yields reduced density matrices that encode classical information, with uncertainties exponentially suppressed by the environment’s entropy. We develop this perspective using both the Lorentzian path integral and operator formalisms and support it with a simple microscopic model. Our results show that quantum gravity in a closed universe naturally gives rise to meaningful, robust predictions without recourse to external constructs.

AdS-CFT Correspondence↗