Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High frequency analysis”

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

One-Step Regeneration of Hairy Roots to Induce High Tanshinone Plants in Salvia miltiorrhiza

Salvia miltiorrhiza is a traditional Chinese medicinal plant of Labiatae, which has been widely utilized to treat a variety of cardiovascular and cerebrovascular diseases. However, due to the long growth cycle, low content of active ingredients, and serious quality deterioration of S. miltiorrhiza , the use of biotechnology to improve S. miltiorrhiza to meet the growing demand for clinical applications has become a research hotspot. In this study, a novel one-step hairy root regeneration method was developed, which could rapidly obtain hairy roots and regenerated plants with high tanshinone content. By optimizing the parameters of Agrobacterium rhizogenes transformation in S. miltiorrhiza , it was finally established that the explants were infected in Ar.qual (OD 600 = 0.6) for 10 min, co-cultured for 3 days, and then screened on the screening medium containing 7.5 mg/l hygromycin, the maximum transformation frequency can reach 73.85%. GFP and PCR detection yielded a total of 9 positive transgenic hairy root lines and 11 positive transgenic regenerated plants. SmGGPPS1 was successfully overexpressed in positive transgenic regenerated plants, according to the results of qRT-PCR. The content of tanshinone IIA and cryptotanshinone were dramatically enhanced in transgenic regenerated plants and hairy roots by Ultra Performance Liquid Chromatography analysis. Based on the Agrobacterium-mediated transformation of S. miltiorrhiza , this study developed a new method for regenerating plants with transgenic hairy roots. This method provides a foundation for the breeding of S. miltiorrhiza and the sustainable development of medicinal plant resources, as well as provides a useful reference for the application of other species.

Su, Yuekai↗
DOE OSTI2025

Discovering $\mu$Hz gravitational waves and ultra-light dark matter with binary resonances

In the presence of a weak gravitational wave (GW) background, astrophysical binary systems act as high-quality resonators, with efficient transfer of energy and momentum between the orbit and a harmonic GW leading to potentially detectable orbital perturbations. In this work, we develop and apply a novel modeling and analysis framework that describes the imprints of GWs on binary systems in a fully time-resolved manner to study the sensitivity of lunar laser ranging, satellite laser ranging, and pulsar timing to both resonant and nonresonant GW backgrounds. We demonstrate that optimal data collection, modeling, and analysis lead to projected sensitivities which are orders of magnitude better than previously appreciated possible, opening up a new possibility for probing the physics-rich but notoriously challenging to access $\mu\mathrm{Hz}$ frequency GWs. We also discuss improved prospects for the detection of the stochastic fluctuations of ultra-light dark matter, which may analogously perturb the binary orbits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Spectral distortions to momentum and scalar exchanges by non-turbulent motion and patchy landscape variability

Modifications to the spectra of turbulent velocity and scalars and co-spectra of vertical fluxes of momentum and scalars due to patchy landscape heterogeneity and non-stationarity are explored for a Mediterranean ecosystem. About 9 months of high frequency measurements of the three velocity components, water vapor concentration, carbon dioxide concentration, and air temperature were analyzed for different seasons (spring/summer) and prevalent wind directions (southeast/northwest). The two wind directions sampled a contrast of clumped and patchy landscape comprised of olive trees (southeast) and wall bounded flow disturbed by the presence of few upwind trees (northwest). The measured spectra and co-spectra were also compared to theoretical scaling forms from stationary, planar homogeneous flow, in the absence of subsidence as derived from the Kansas experiment. To assess the role of low frequency non-turbulent motion on the spectral and co-spectral content, a 5-min Fourier cutoff was introduced and the analysis was limited to near-neutral conditions where the boundary layer depth is shallow compared to its unstable counterpart. It was shown that the velocity statistics were not appreciably impacted by the low-frequency motion causing non-stationarity. Moreover, the turbulent scalar fluxes were also shown not to be significantly impacted by such low frequency motion. The scalar variances were impacted, especially the water vapor variance and its concomitant spectral shape. When the non-turbulent motion was filtered, the scalar spectra at low wavenumbers followed expectations from the so-called attached eddy hypothesis (i.e. exhibited a $k^{-1}_x$ scaling with $k_x$ defining the longitudinal wavenumber) applicable for near-neutral conditions. For momentum co-spectra, the canonical shapes from the Kansas experiment appear to describe well the measurements here and in both dominant directions and seasons with some adjustment to the integral time scales based on wind direction. For the scalar co-spectra, deviations from the Kansas experiment were prevalent. The most noticeable and surprising deviations were their slow decay with increased sampling frequency at inertial subrange scales. This slow decay was shown not to contribute appreciably to the overall scalar fluxes. At those fine scales, predictions from local isotropy were expected to hold. The scalar co-spectral deviations from local isotropy were then discussed using a simplified co-spectral budget model where scalar–scalar co-spectra naturally emerged and the interplay between landscape heterogeneity and a scale-dependent pressure-scalar de-correlation time was postulated. It is also envisaged that the findings here offer a preliminary template for analyzing eddy-covariance data in situations that deviate from ideal conditions, especially regarding low-frequency modulations of scalar spectra and vertical scalar flux co-spectra.

Canopy turbulence↗
DOE OSTI2025

Numerical Simulation Studies of Ultrasonic De-Icing for Heating, Ventilation, Air Conditioning, and Refrigeration Structures

Ice accumulation on heating, ventilation, air conditioning, and refrigeration (HVACR) structures presents significant operational challenges. These challenges include reduced efficiency, increased energy consumption, and potential damage to equipment. Traditional de-icing methods, such as chemical treatments, mechanical scraping, or heating-based techniques, are often labor-intensive, costly, and environmentally harmful. Here, this study uniquely investigates ultrasonic de-icing as an energy-efficient alternative for HVACR applications, focusing on the specific structural geometries found in these systems. A comprehensive numerical simulation framework was developed using finite element analysis to explore ultrasonic wave propagation across four distinct HVACR structures. Key parameters such as ultrasonic frequency, power levels, and the number and placement of actuators were examined for their impact on ice detachment efficiency. Results from simulations on a plate structure reveal that ultrasonic excitation can propagate effectively across large areas (at least 150 × 150 mm), enhancing the de-icing coverage. Lower frequency (e.g., 30 to 45 kHz) excitation results in greater displacement, improving de-icing performance, while increased actuator numbers with the same total power input also enhance effectiveness. Two actuators seem sufficient for the de-icing of a 300 × 300 mm plate. For tube-and-fin structures, specific high-power ultrasonic frequencies selectively excite the fin plates, demonstrating efficient ice removal when actuated on the tube. However, optimal performance requires careful design of actuator placement and vibration modes to accommodate the irregular shapes of these structures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Multifrequency Analysis of Favored Models for the Messier 87* Accretion Flow

Abstract The polarized images of the supermassive black hole Messier 87* (M87*) produced by the Event Horizon Telescope (EHT) provide a direct view of the near-horizon emission from a black hole accretion and jet system. The EHT theoretical analysis of the polarized M87* images compared thousands of snapshots from numerical models with a variety of spins, magnetization states, viewing inclinations, and electron energy distributions, and found a small subset consistent with the observed image. In this article, we examine two models favored by EHT analyses: a magnetically arrested disk with moderate retrograde spin and a magnetically arrested disk with high prograde spin. Both have electron distribution functions that lead to strong depolarization by cold electrons. We ray trace five snapshots from each model at 22, 43, 86, 230, 345, and 690 GHz to forecast future very long baseline interferometry (VLBI) observations and examine limitations in numerical models. We find that even at low frequencies where optical and Faraday rotation depths are large, approximately rotationally symmetric polarization persists, suggesting that shallow depths dominate the polarization signal. However, morphology and spectra suggest that the assumed thermal electron distribution is not adequate to describe emission from the jet. We find 86 GHz images show a ringlike shape determined by a combination of plasma and spacetime imprints, smaller in diameter than recent results from the Global mm-VLBI Array. We find that the photon ring becomes more apparent with increasing frequency, and is more apparent in the retrograde model, leading to large differences between models in asymmetry and polarization structure.

Astronomy & Astrophysics↗
DOE OSTI2025

Solid State Transformer Architecture and Control Compensation for Common Mode Currents

A high-altitude electromagnetic pulse (HEMP) or similar geomagnetic disturbance (GMD) has the potential to impact the operation of large-scale electric power grids. By introducing low-frequency common-mode (CM) currents, these events can degrade the performance of critical system components, such as large power transformers by introducing CM currents which can lead to magnetic saturation of the transformer core. In this work, a solid-state transformer (SST) is developed to replace susceptible equipment and improve grid resiliency by safely absorbing these CM disturbances. This device will be referred to as a common-mode solid-state transformer (CM-SST). An SST architecture based on a four-legged AC/DC converter is developed. This architecture enables active control of CM signals without disturbing the AC voltages or the real and reactive power delivery capabilities. A system-level model of this architecture is created, and time-domain simulations are performed to evaluate the SST’s performance in response to simulated CM disturbances. A control strategy for mitigating CM current is also investigated. Hamiltonian surface shaping and power flow control (HSSPFC) is used to design a nonlinear controller for the SST’s output inverter. The objectives of the controller are to suppress CM-induced AC current offsets and regulate AC currents to desired setpoints. Nonlinear system analysis is applied to design and validate the controller. Two cases are tested: (a) the proposed four-leg inverter and (b) a standard three-leg inverter. The results show that the proposed controller rapidly mitigates CM disturbances while maintaining high-quality AC current waveforms in the four-leg configuration. Finally, the hardware performance of an SST prototype is evaluated. In particular, the ability of the SST to safely redirect and absorb CM currents is demonstrated, showing how it can protect neighboring conventional transformers in the system. The study confirms that appropriate control laws allow the SST to protect both itself and adjacent transformers during a HEMP or GMD event.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Fast ion transport by sawtooth instability in the presence of ICRF-NBI synergy in JET plasmas

JET experiments have shown that the three-ion scenarios using waves in the ion cyclotron range of frequencies (ICRF) is an efficient way to build fast ion population through beam ion acceleration by radio frequency (RF) waves. Here, such a heating scheme is applied to plasmas with at least two thermal ion species. Analysis of mixed discharges with complex heating schemes requires a workflow that allows to model thermal and fast ion transport consistently. This paper is dedicated to modelling of a mixed plasma discharge with significant fraction of fast ions and contributes to development of fast ion transport models. For interpretive analysis with the TRANSP code a JET hydrogen-deuterium (HD) plasma discharge with neutral beam injection (NBI) and ICRF heating has been chosen. The task is complicated by NBI-ICRF synergy and plasma magnetohydrodynamic activity, like sawtooth crashes. D beam ions accelerated by RF waves form a high energy tail in fast ion distribution. Significant difference between the neutron rate computed by TRANSP and measured one is observed if the same diffusivity for electrons and ions is assumed. Sensitivity studies show that uncertainties in input plasma parameters and thermal ion transport models are crucial for modelling mixed plasma discharges and increased D transport is required to reach the plasma composition consistent with diagnostic measurements at the plasma edge. Fast ion redistribution by a sawtooth instability is characterised by non-resonant transport due to reconnection of magnetic field lines and resonant transport caused by resonance interaction between the instability and fast ions. With ORBIT simulations it has been shown that resonant interaction strongly affects fast ions of high energies, like beam ions accelerated by RF waves and fusion products. For the considered case, fast ion profiles simulated by ORBIT remain peaked after the sawtooth crashes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Liquid helium fluid dynamics studies. Final Technical Report

Future high energy physics accelerators depend on a number of advanced technologies to open the many doors of scientific discovery. Among these advanced technologies, superconducting magnets and superconducting radio frequency (SRF) cavities are the backbone of the accelerator and detector systems. But all these low temperature systems depend critically on successful and reliable operation of their supporting technologies, among which the liquid helium cooling system is of the utmost importance. To improve the quality of these systems both in terms of efficiency and reliability, a robust helium cryogenics research and development (R&D) effort is required. The proposed research to be conducted by the FSU cryogenics group aims to produce fundamental knowledge that meets this R&D need. The projects that we have completed over the past grant period at Florida State University consist of experimental research on liquid helium fluid dynamics and heat transfer problems relevant to the development of future superconducting particle physics accelerators. Liquid helium is the coolant used in all such facilities and in many of these facilities He II (the low temperature phase of liquid helium also known as superfluid helium) is preferred due to its outstanding heat transfer characteristics. The work consists of two main experimental studies that probe both fundamental as well as practical aspects of liquid helium cooling. The first is a broad and fundamental study of the heat and mass transfer processes that can occur during a sudden catastrophic loss of vacuum (SCLV) incident in a superconducting accelerator. SCLV refers to the remote but extremely critical accident scenario where atmospheric pressure air is allowed to flood into the insulating vacuum system and impinge on the liquid helium cooled surfaces in the accelerator. Safe performance and recovery from such accidents is essential to the reliable operation of superconducting accelerators. The dynamics of this process is quite complex and so our approach is to conduct a series of well-orchestrated experiments that probe the various physical phenomena that can occur during an SCLV event. The experiments are coupled with analytic and numerical analysis in an effort to develop a general understanding of the process and to assist with future accelerator design and development. The second activity is directed toward fundamental understanding of heat and mass transfer in He II, which is essential to the design of superconducting magnets and radio frequency cavities in accelerators. The work consists of flow visualization of the dynamics of He II using laser assisted techniques. Two complementary techniques are used to study the fundamentals of the turbulent state. The first technique uses neutrally buoyant solid hydrogen particles to probe the flow fields of the superfluid and normal fluid components. The other technique uses laser excited He2* molecules as tracers of the normal fluid motion within the He II. The activities also included an effort to use visualization techniques to locate transient hot spots in radio frequency superconducting cavities. Such work provides valuable information about the heat transfer process in He II and its impact on the performance of superconducting devices. The research effort at Florida State University is not directly in support of a specific high energy physics experiment or facility. Rather, the work is general and coordinated with HEP accelerator laboratories to provide valuable insight that can assist with future accelerator development.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Performance of Fire Danger Indices and Their Utility in Predicting Future Wildfire Danger Over the Conterminous United States

Abstract Predicting current and future wildfire frequency and size is central to wildfire control and management. Multiple fire danger indices (FDIs) that incorporate weather and fuel conditions have been developed and utilized to support wildfire predictions and risk assessment. However, the scale‐dependent performance of individual FDIs remains poorly understood, which leads to large uncertainty in the estimated fire sizes under climate change. Here, we calculate four commonly used FDIs over the conterminous United States using high‐resolution (4 km) climate and fuel data sets for the 1984–2019 period. The relationships of these four FDIs to the observed wildfire sizes show that higher values of FDIs correlate to larger total fire sizes; this correlation is more robust at larger spatial scales. Sensitivity analysis indicates that the daily minimum relative humidity and precipitation are the most important drivers of the annual mean fire danger. In the instances of extreme fire danger, wind speed becomes a critical factor and should be considered in the calculation of the FDI. To assess the impact of climate change on future fire size, we calculate the present‐day and end‐of‐century FDIs using the 12 km regional climate model simulations. The four FDIs generally predict consistent changes in future fire potential, suggesting an overall higher fire potential in conjunction with a prolonged wildfire season in future climate. Regionally, the four FDIs also reveal similar seasonal patterns as the enhancement arises mostly in spring and summer over the southwest US while in summer and fall over the northern and eastern US.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Online real-time learning of dynamical systems from noisy streaming data

Abstract Recent advancements in sensing and communication facilitate obtaining high-frequency real-time data from various physical systems like power networks, climate systems, biological networks, etc. However, since the data are recorded by physical sensors, it is natural that the obtained data is corrupted by measurement noise. In this paper, we present a novel algorithm for online real-time learning of dynamical systems from noisy time-series data, which employs the Robust Koopman operator framework to mitigate the effect of measurement noise. The proposed algorithm has three main advantages: (a) it allows for online real-time monitoring of a dynamical system; (b) it obtains a linear representation of the underlying dynamical system, thus enabling the user to use linear systems theory for analysis and control of the system; (c) it is computationally fast and less intensive than the popular extended dynamic mode decomposition (EDMD) algorithm. We illustrate the efficiency of the proposed algorithm by applying it to identify the Van der Pol oscillator, the chaotic attractor of the Henon map, the IEEE 68 bus system, and a ring network of Van der Pol oscillators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

NNSS Soils Monitoring: Plutonium Valley (CAU 366) FY2019

Desert Research Institute (DRI) conducted a field assessment of the potential for radionuclide-contaminated soil to be transported from the Plutonium Valley Dispersion Sites Contamination Area (CA) because of both wind and storm water runoff. This activity supported U.S. Department of Energy (DOE) Environmental Management Nevada Program (EM NV) efforts to establish post-closure monitoring plans for the Plutonium Valley Dispersion Sites Corrective Action Unit (CAU) 366. The DRI task was intended to identify the likely mechanism(s) of transport and determine the meteorological conditions that might cause the movement of radionuclide-contaminated soils. The emphasis of the work was on collecting sediment transported by channelized storm runoff and measuring airborne dust concentrations and associated wind conditions. These data will facilitate an appropriate post-closure monitoring program. In 2011, DRI installed meteorological monitoring stations north and south of the Plutonium Valley CA—known as Pu Valley North and Pu Valley South, respectively—as well as a fluvial sediment sampling station within the CA. Since installation, temperature, wind speed, wind direction, relative humidity, precipitation, solar radiation, barometric pressure, soil temperature, volumetric soil moisture content, and airborne particulate concentrations have been collected at both meteorological stations. The maximum, minimum, and average or total (as appropriate) for each of these parameters were recorded for each 10-minute interval. The sediment sampling station included an automatically activated sampling pump with collection bottles for suspended sediment, which was activated when sufficient flow was present in the channel, and passive traps for bedload material that was transported down the channel during runoff events. This report presents data collected from these stations during fiscal year (FY) 2019. During the FY2019 (October 1, 2018, through September 30, 2019) reporting period, the warmest month was July and the coldest month was February. Monthly total solar radiation was highest in June or July (depending on the station) and lowest in December. Monthly mean wind speeds were highest during February. At the Pu Valley North station, winds were commonly northerly, northwesterly, and southerly. At the Pu Valley South station, the wind direction was more variable and winds frequently came from most directions, excluding westerly and northwesterly winds. Winds above 15.0 miles per hour (mph) (24.1 kilometers per hour [km/hr]) were frequently southerly at both stations, and less-frequent strong winds were northerly (at the Pu Valley South station) or northwesterly (at the Pu Valley North station). Monthly average relative humidity ranged from 18 percent in August to 62 percent in January and February. Monthly total precipitation ranged from zero at both stations in September to 2.23 inches (in) (56.6 millimeters [mm]) during February at the Pu Valley South station and 2.09 in (53.1 mm) during March at the Pu Valley North station. From October 1, 2018, through September 30, 2019, the total precipitation was 8.23 in (209 mm) and 8.15 in (207 mm) at the Pu Valley South and Pu Valley North stations, respectively. The largest daily precipitation totals of 1.46 in (37.1 mm) and 1.14 in (29.0 mm) occurred at the Pu Valley South and Pu Valley North stations, respectively, on March 6, 2019. However, the ISCO autosampler was not activated during any storm event in FY2019 because the channel water depth measurements were insufficient for suspended sediment sample collection. Additionally, because no significant flow through the channel was indicated, no bedload samples were collected. Light breezes of 0.0 mph to 5.0 mph (0.0 km/hr to 8.0 km/hr) occurred most frequently (approximately 55 percent to 57 percent of the time). The frequency of occurrence diminished exponentially as the wind speed increased such that winds in excess of 25.0 mph (40.2 km/hr) occurred less than 0.01 percent of the time. Winds in excess of 15.0 mph (24.1 km/hr) were most commonly southerly, which is likely controlled by the topography of the valley because topographic highs define the east and west sides of the valley and converge toward the north.The concentrations of PM2.5 (particulate matter with an aerodynamic diameter ≤2.5 micrometers [μm]) and PM10 (particulate matter with an aerodynamic diameter ≤10 μm) in the air generally increase as wind speed increases. However, the particle profiler at the Pu Valley South station regularly malfunctioned, which resulted in poor data quality and a limited period of record within FY2019. Therefore, discussion of the annual dust observations is mostly limited to the Pu Valley North station. At the Pu Valley North station, significant increases in windblown dust concentrations were observed when wind speeds exceeded 15.0 mph (24.1 km/hr), especially during southerly winds. When all wind speeds were considered, high dust concentrations at the Pu Valley North station were frequently associated with northerly, northwesterly, or southerly winds. The majority of PM10 transport at the Pu Valley North station occurred when winds were blowing northerly and northwesterly because of the high frequency of winds up to 15.0 mph (24.1 km/hr). The ratio of PM10 to PM2.5 is a qualitative indicator of the proximity of dust sources to the observation point and the values of the PM10 to PM2.5 ratio tend to be higher nearer to the source. The PM10 to PM2.5 ratio increased from 1.2 to 2.8 with higher wind speeds at the Pu Valley North station. This increase suggests that some locally sourced PM10 was present in the air when local winds exceeded 15.0 mph (24.1 km/hr). Local dust suspension is interpreted during major dust events when peaks in the PM10 to PM2.5 ratio occur simultaneously with the maximum wind speed and the peak PM10 concentration. Analysis of the five major dust events in FY2019 revealed that northerly winds were dominant in three events and southerly winds were dominant in two events. Local suspension of dust within Plutonium Valley is suspected to have occurred during two of the five major events (April 9, 2019, and September 28 to 29, 2019), and a nonlocal source of dust is presumed for the three other major dust events (October 14, 2018; May 16, 2019; and September 2, 2019). The three suspected regional dust transport events exhibited winds that typically did not exceed the 15.0 mph (24.1 km/hr) wind speed threshold for local dust transport and featured PM10 to PM2.5 ratios that indicate large-scale regional dust transport. Similar event characteristics (e.g., wind speed and PM10 concentration) were observed at air monitoring stations approximately 70 mi (113 km) to the northwest within the Tonopah Test Range (TTR). Therefore, these three regional events may not reflect transport from the Plutonium Valley CA. Routine monitoring of meteorological and hydrologic parameters at Plutonium Valley has been discontinued at the request of EM NV, effective at the conclusion of FY2019. Environmental data acquisition has been completed, and both Pu Valley North and Pu Valley South meteorological stations within Plutonium Valley as well as the ISCO autosampler station within the CA were removed on October 14, 2019.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Method for Coupled Electromagnetic and Circuit Simulations to Evaluate Surge Arrester Performance in Protecting Equipment Against E1 HEMP

Surge arrester behavioral modeling for realistic systems embedded in an E1 high-altitude electromagnetic pulse environment inherently encompasses three interconnected complications: (1) the need to account for signal propagation across two domains, electromagnetics and electrical; (2) the need to include both linear and nonlinear circuit components in the analysis; and (3) the need to understand that the over-current and over-voltage mitigation performance is dependent not only on the properties of the surge arrester and protected load but also on the topology of the overall electrical network. This study presents a framework to address these challenges in a systematic manner to consider the effectiveness of protective measures for a common class of equipment in power generation facilities. Full-wave simulations were carried out to derive circuit-domain (i.e., lumped element–based) equivalent models for the excitation waveform and the physical components of the system. Then, these equivalent models were imported to a circuit solver and combined with a high-frequency surge arrester model to evaluate mitigation performance. The methodology outlined is general enough such that it can be applied for other electromagnetic interference problems that involve E2/E3 HEMP or microwave emissions.

42 ENGINEERING↗
DOE OSTI2022

Forced Oscillation Grid Vulnerability Analysis and Mitigation Using Inverter-Based Resources: Texas Grid Case Study

Forced oscillation events have become a challenging problem with the increasing penetration of renewable and other inverter-based resources (IBRs), especially when the forced oscillation frequency coincides with the dominant natural oscillation frequency. A severe forced oscillation event can deteriorate power system dynamic stability, damage equipment, and limit power transfer capability. This paper proposes a two-dimension scanning forced oscillation grid vulnerability analysis method to identify areas/zones in the system that are critical to forced oscillation. These critical areas/zones can be further considered as effective actuator locations for the deployment of forced oscillation damping controllers. Additionally, active power modulation control through IBRs is also proposed to reduce the forced oscillation impact on the entire grid. The proposed methods are demonstrated through a case study on a synthetic Texas power system model. The simulation results demonstrate that the critical areas/zones of forced oscillation are related to the areas that highly participate in the natural oscillations and the proposed oscillation damping controller through IBRs can effectively reduce the forced oscillation impact in the entire system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Developing Frequency Stability Constraint for Unit Commitment Problem Considering High Penetration of Renewables

As zero-carbon electricity systems become the trend of future grid, the system inertia provided by conventional synchronous generators (SGs) keeps decreasing. The resultant lower system inertia will inevitably cause frequency stability problem, especially in the first few seconds following disturbance. To tackle this challenge, this paper proposes a frequency stability constraint for power systems unit commitment problem by considering the fast frequency responses (FFRs) from inverter-based resources (IBRs). Our developed frequency stability constraint is grounded on an analytical frequency nadir estimation framework that considers both SG and IBR dynamics. The accuracy of our frequency nadir estimation framework is validated by most severe N-1 contingency simulation result in a real island system. Then, the adaptive inertia frequency stability constraint is derived by performing sensitivity analysis with our frequency nadir estimation framework. Finally, we demonstrate the effectiveness of our developed frequency stability constraint with one year day-ahead unit commitment results of the island system.

fast frequency response↗
DOE OSTI2021

The picasso map-making code: application to a simulation of the QUIJOTE northern sky survey

ABSTRACT Map-making is an important step for the data analysis of cosmic microwave background (CMB) experiments. It consists of converting the data, which are typically a long, complex, and noisy collection of measurements, into a map, which is an image of the observed sky. We present in this paper a new map-making code named picasso (Polarization and Intensity CArtographer for Scanned Sky Observations), which was implemented to construct intensity and polarization maps from the Multi Frequency Instrument (MFI) of the QUIJOTE (Q-U-I Joint TEnerife) CMB polarization experiment. picasso is based on the destriping algorithm, and is suited to address specific issues of ground-based microwave observations, with a technique that allows the fit of a template function in the time domain, during the map-making step. This paper describes the picasso code, validating it with simulations and assessing its performance. For this purpose, we produced realistic simulations of the QUIJOTE-MFI survey of the northern sky (approximately ∼20 000 deg2), and analysed the reconstructed maps with picasso, using real and harmonic space statistics. We show that, for this sky area, picasso is able to reconstruct, with high fidelity, the injected signal, recovering all the scales with ℓ > 10 in TT, EE, and BB. The signal error is better than 0.001 per cent at 20 < ℓ < 200. Finally, we validated some of the methods that will be applied to the real wide-survey data, like the detection of the CMB anisotropies via cross-correlation analyses. Despite that the implementation of picasso is specific for QUIJOTE-MFI data, it could be adapted to other experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Large-scale simulations of Floquet physics on near-term quantum computers

Abstract Periodically driven quantum systems exhibit a diverse set of phenomena but are more challenging to simulate than their equilibrium counterparts. Here, we introduce the Quantum High-Frequency Floquet Simulation (QHiFFS) algorithm as a method to simulate fast-driven quantum systems on quantum hardware. Central to QHiFFS is the concept of a kick operator which transforms the system into a basis where the dynamics is governed by a time-independent effective Hamiltonian. This allows prior methods for time-independent simulation to be lifted to simulate Floquet systems. We use the periodically driven biaxial next-nearest neighbor Ising (BNNNI) model, a natural test bed for quantum frustrated magnetism and criticality, as a case study to illustrate our algorithm. We implemented a 20-qubit simulation of the driven two-dimensional BNNNI model on Quantinuum’s trapped ion quantum computer. Our error analysis shows that QHiFFS exhibits not only a cubic advantage in driving frequency ω but also a linear advantage in simulation time t compared to Trotterization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Uncovering fast solid-acid proton conductors based on dynamics of polyanion groups and proton bonding strength

Achieving high proton conductivity in inorganic solids is key for advancing many electrochemical technologies, including low-energy nano-electronics and energy-efficient fuel cells and electrolyzers. A quantitative understanding of the physical traits of a material that regulate proton diffusion is necessary for accelerating the discovery of fast proton conductors. In this work, we have mapped the structural, chemical and dynamic properties of solid acids to the elementary steps of the Grotthuss mechanism of proton diffusion. Our approach combines ab initio molecular dynamics simulations, analysis of phonon spectra and atomic structure calculations. We have identified the donor–hydrogen bond lengths and the acidity of polyanion groups as key descriptors of local proton transfer and the vibrational frequencies of the cation framework as the key descriptor of lattice flexibility. The latter facilitates rotations of polyanion groups and long-range proton migration in solid acid proton conductors. The calculated lattice flexibility also correlates with the experimentally reported superprotonic transition temperatures. Using these descriptors, we have screened the Materials Project database and identified potential solid acid proton conductors with monovalent, divalent and trivalent cations, including Ag + , Sr 2+ , Ba 2+ and Er 3+ cations, which go beyond the traditionally considered monovalent alkali cations (Cs + , Rb + , K + , and NH 4 + ) in solid acids.

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

STEPS: A Portable Numerical Simulation Toolkit for Electrical Power System Dynamic Studies

Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.

42 ENGINEERING↗