Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-varying 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Physics-constrained superresolution diffusion for six-dimensional phase space diagnostics

Adaptive physics-constrained superresolution diffusion is developed for noninvasive virtual diagnostics of the six-dimensional (6D) phase space density of charged particle beams. An adaptive variational autoencoder embeds initial beam condition images and scalar measurements to a low-dimensional latent space from which a 32 6 pixel 6D tensor representation of the beam's 6D phase space density is generated. Projecting from a 6D tensor generates physically consistent two-dimensional projections. Physics-guided superresolution diffusion transforms low-resolution images of the 6D density to high resolution 256 × 256 pixel images. Unsupervised adaptive latent space tuning enables tracking of time-varying beams without knowledge of time-varying initial conditions. The method is demonstrated with experimental data and multiparticle simulations at the HiRES UED. The general approach is applicable to a wide range of complex dynamic systems evolving in high-dimensional phase space. The method is shown to be robust to distribution shift without retraining. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Optical emission spectroscopy and imaging of low-pressure N2 plasmas generated by intense fast-pulsed electron beams

An optical emission spectroscopic (OES) and imaging characterization is conducted on N2 plasmas generated by a 100 keV fast-pulsed electron beam. The electron beams are injected into an N2 gas filled volume with a current of 4.5 kA (300 A/cm2) and a 100 ns pulse width. The characterization is conducted at the pressures, 1 Torr and 0.1 Torr, corresponding to two distinct regimes that exhibit significantly different plasma dynamics. Beam impact ionization is shown to be a primary mechanism for producing low temperature plasmas at 1 Torr during beam output. After beam termination, ionization by an inductive electric field becomes the primary mechanism for plasma formation later in time for both pressures. OES and plasma imaging are used in this work as a diagnostic tool to track the distribution of electronic, vibrational, and rotational state transitions and ionized species. This is achieved with the use of a multi-resolution suite of spectrometers capable of acquiring time-resolved spectra. Vibrational and rotational bands of the N2 second positive system (C3Πu→B3Πg) and the N2+ first negative system (B2Σu+→X2Σg+) are identified in both regimes as well as N+ states exclusively in the lower pressure regime. Vibrational and rotational spectra are shown to track the evolution of the time-varying plasma current. Plasma imaging also reveals spatially nonuniform plasma emission at 0.1 Torr. A few hundred shots are recorded to fully characterize the emissions, and results are shown to be highly reproducible (≤±1% with 2σ confidence).

Kaiser, E. R. (ORCID:0000000336493938)↗
DOE OSTI2026

Heat exposure and maternal stress: evidence from the GRAPHS pregnancy cohort in Ghana

Heat exposure has been linked to psychosocial stress, an established antecedent of perinatal depression; however, evidence on heat-related stress during pregnancy in sub-Saharan Africa remains limited. We analyzed psychosocial stress scores and covariate data from the Ghana Randomized Air Pollution and Health Study, linking daily maximum and minimum shaded wet bulb globe temperature (WBGT) metrics to participants’ stress scores derived from the Crisis in Family Systems-Revised Life Events Questionnaire. We evaluated associations using ordinal logistic regression of pregnancy-average and trimester-average exposures and distributed lag non-linear models (DLNMs) to assess time-varying associations across gestation. Higher average maximum WBGT exposure across pregnancy was associated with increased odds of higher psychosocial stress; each 1 °C increase in maximum WBGT was associated with 64% higher odds of belonging to a higher stress category (OR = 1.64; 95% CI = 1.17–2.31; p = 0.0040). In trimester-average models, higher first-trimester maximum WBGT was also associated with higher stress (OR = 1.44; 95% CI = 1.15–1.81; p = 0.0014). DLNMs suggested that relatively cooler daily maximum WBGT values (25th percentile) were associated with decreased odds of stress in early pregnancy, whereas extreme daily maximum WBGT values (99th percentile) showed a pattern consistent with increased odds of stress in mid-to-late gestation. These findings highlight gestational windows in which heat exposure may influence stress, emphasizing the need for further research into underlying mechanisms and effective interventions to protect maternal mental health in heat-vulnerable settings.

White, Lewis [Columbia University] (ORCID:00090005↗
DOE OSTI2026

Nanoscale strain wave generation by a piezoelectric grating from polar vortices

Nanostructures formed by spontaneously broken symmetry have provided new ways to manipulate quantum states. Specifically, topological structures with periodic spatial ordering, such as polar vortices and skyrmions, can be ideal hosts for creating engineered responses in both spatial and frequency domains. So far, however, only a few examples of such hierarchical engineering have been reported in the literature. Here we demonstrate that the spatially modulated piezoelectric response of a polar vortex structure can create strain waves with a characteristic nanoscale wavefront. Using time-resolved pump–probe resonant X-ray scattering and diffraction measurements, coupled with dynamical phase-field simulations, we show that the piezoelectric modulation of the spontaneously formed polar vortex crystal functions as an acoustic diffraction grating. This system converts incoming laterally uniform strain waves into outgoing waves with a characteristic sub-terahertz frequency, driven by an intrinsic excitation of the polar vortex crystal. Moreover, our phase-field simulations suggest that the dynamic mechanical displacements exhibiting vortex textures are generated from both space- and time-varying piezoelectric responses. Our findings illustrate a new method for generating nanoscale strain waves with unique spatial textures by tuning the hierarchical order of polar topologies to engineer new collective modes, allowing for a wide range of control through the topological lattice.

ferroelectrics↗
DOE OSTI2024

Stochastic finite volume method for uncertainty quantification of transient flow in gas pipeline networks

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. In conclusion, we demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Planar shock-induced bubble collapse and jetting in water captured via X-ray phase contrast imaging

Shock wave–bubble interactions in water manifest rich dynamics driven by a combination of strong pressure and density mismatches. They have a wide variety of applications, including the injection of pharmaceuticals, and through scaling, enable the exploration of various aspects of high-energy-density systems such as inertial confinement fusion. In this work, the interaction between a micrometric nitrogen bubble and a planar shock wave, characterized by a Mach number of M = 1.24 and a peak pressure of p max = 0.57 , is experimentally recorded using ultra-high-speed x-ray phase contrast imaging. Highly resolved radiographs provide access to all phase discontinuities along the beam path, offering quantities such as the time-varying bubble size, the speed of a jet produced during the bubble collapse, and the time evolution of the shock wave front, which are critical benchmark data for numerical scheme validation. This study addresses the lack of well-characterized, repeatable, and high spatiotemporal resolution experiments at negative Atwood numbers by providing shock–bubble visualization and corresponding numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Influence of background sources and topographic resolution in the Weather Research and Forecasting Model on xenon plume characteristics at monitoring stations

For many atmospheric monitoring applications, networks of measurement sites—such as the radionuclide stations of the International Monitoring System—can be sparse. With measurement locations potentially hundreds to thousands of kilometers from a release it is important to quantify the effects of physical processes on transport and dispersion of plumes between source and measurement locations. This study addresses the effects of background sources and topography resolution near the release location of radionuclides. We use the Weather Research and Forecasting (WRF) model with inline chemistry to investigate (1) how an additional, time-varying source of 133 Xe, such as an operational medical isotope production facility, contributes to activity concentration measurements at monitoring sites, and (2) how complex topography influences on atmospheric conditions near emission sources impact plume concentrations at varying distances from the source. Two 133 Xe emission sources, including (1) a high flux rate of short duration representing an explosive event, and (2) a variable and continuous background source, are simulated. The continuous background source contributes significantly to total 133 Xe concentrations at several monitoring stations. Further, a WRF simulation at 9 km horizontal resolution is compared with a nested grid simulation, where the innermost domain has a resolution of 1 km. Increased topographic resolution leads to an improved representation of plume responses to local winds, with topographic influences greatest at locations closest to the sources. Differences between the two domain resolutions decrease at greater distances from the sources, as plumes have time to spread and mix and are influenced by synoptic scale circulation patterns that are represented similarly in both simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗
DOE OSTI2024

Implementation and validation of optimal start control strategy for air conditioners and heat pumps

Commercial buildings are responsible for approximately 20 % of the total energy consumption and greenhouse gas emissions in the United States. Over 85 % of these buildings lack building automation systems, and many are small (<50,000 square feet), underserved, and use rooftop units (RTUs) for heating, ventilation, and air-conditioning needs. Because these buildings lack proper energy management systems, several operational deficiencies lead to excess energy consumption. Studies have shown that managing the RTUs’ heating and cooling set points, schedules, setbacks, and optimal start can result in a 20 % to 25 % reduction in electricity consumption in small commercial buildings. These buildings typically use fixed schedules to start the RTUs 60 to 120 min before occupancy begins, which results in excess energy consumption. This paper presents research that demonstrates and evaluates the performance of four optimal start methods, which utilize data-based modeling as a key element in facilitating adaptive control in response to time-varying inputs while requiring minimal sensor inputs. The evaluation found energy savings in two commercial buildings equipped with RTUs by periodically alternating four different optimal start models during the cooling and heating season. The resulting energy savings are positive for all models and range from 2 to 5 kWh/day/unit. The units on the east side of the building showed higher savings, while interior units showed greater variability in savings due to the differences in capacities and room sizes. Savings were considerably greater during the heating season compared to the cooling season. The performance of all four models on Mondays was poor; models suggested a shorter optimal start time, which resulted in relatively larger errors. Finally, the future work will look at using a different model for the days after weekends and holidays.

42 ENGINEERING↗
DOE OSTI2025

Incorporating civilian radioxenon background estimates in anomaly detection

A nuclear explosion screening exercise in 2023 (Maurer et al., 2023) found challenges with discerning anomalous radioxenon activity concentrations relative to elevated background concentrations. Research has continued into methods to detect anomalous radioxenon concentrations by comparing samples to estimates of atmospheric radioxenon background concentrations caused by releases at nuclear reactors or medical isotope production facilities. A new approach estimates the sample concentrations using time-varying radioxenon release rates obtained using optimization techniques that constrain the facility release rates to plausible amounts based on historical data or facility knowledge. The purpose of the optimization is to determine whether any combination of plausible release rates from emitting facilities can explain a series of radioxenon measurements at one or more sampling stations. A case study uses radioxenon data collected at three locations in western Europe for a month in 2021 and considers releases from 77 locations. Fewer samples are identified as being anomalous using a simplistic flagging rule than from an application of the current International Monitoring System (IMS) activity concentration-level rule.

Environmental sciences↗
DOE OSTI2024

Accretion onto a Supermassive Black Hole Binary before Merger

While supermassive binary black holes (SMBBHs) inspiral toward merger they may also accrete matter from a surrounding disk. To study the dynamics of this system requires simultaneously describing the evolving spacetime and the magnetized plasma. We present the first relativistic calculation simulating two equal-mass, nonspinning black holes as they inspiral from a 20 M (G = c = 1) initial separation almost to merger. Our results imply important observational consequences: for instance, the accretion rate $\dot{M}$ onto the black holes first decreases and then plateaus, dropping by only a factor of ~3 despite the rapid inspiral. An estimated bolometric light curve follows the same profile, suggesting some merging SMBBHs may be significantly luminous past the predicted circumbinary disk decoupling. The minidisks are nonstandard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, "sloshing" streams transfer mass between minidisks, carrying kinetic energy at a rate sometimes as high as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives time-varying minidisk tilts. These complex dynamics all contribute to unique cyclical behavior in the light curves of late-time inspiraling SMBBHs. The poloidal magnetic flux on the black holes is roughly constant at a dimensionless level $\phi$ ~ 2–3, but doubles just before merger; for significant black hole spin, this flux predicts powerful jets with variability driven by binary dynamics, another potentially unique electromagnetic signature. This simulation is the first to employ our multipatch infrastructure PATCHWORK MHD, decreasing the computational expense to ~3% of conventional single-grid methods' cost.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nuclear Materials Packaging, Transportation, and Systems Analysis Group Software Quality Assurance Plan: ANSYS Mechanical Finite Element Analysis Software Version 2023R1

ANSYS Inc. develops and markets engineering simulation software and services used in the aerospace, automotive, manufacturing, electronics, biomedical, energy, defense, and many other industries. ANSYS is dedicated to engineering simulation and is the world’s leading software provider. ANSYS was founded in 1970 and is headquartered in Canonsburg, Pennsylvania. ANSYS provides an engineering analysis tool combining structural, thermal, computational fluid dynamics, acoustic, and electromagnetic simulation capabilities. ANSYS has two main programs, which use the same solvers: (1) Mechanical APDL (ANSYS Design Parametric Language), a Fortran-based coding platform, and (2) ANSYS Workbench, which uses a graphical user interface to aid in finite element analysis implementation. This plan covers both APDL and Workbench. The ANSYS computer program is a large-scale, multipurpose finite element program that can be used to solve several classes of engineering analyses. The analysis capabilities of ANSYS include the ability to solve static and dynamic structural analyses, steady-state and transient heat transfer problems, mode-frequency and buckling eigenvalue problems, static or time-varying magnetic analyses, and various types of field and coupled-field applications. The program contains many special features that allow nonlinearities or secondary effects such as plasticity, large strain, hyperelasticity, creep, swelling, large deflections, contact, stress stiffening, temperature dependency, material anisotropy, and radiation to be included in the solution. As ANSYS has been developed, other special capabilities such as substructuring, submodeling, random vibration, kinetostatics, kinetodynamics, free convection fluid analysis, acoustics, magnetics, piezoelectrics, coupled-field analysis, and design optimization have been added to the program. These capabilities contribute further to making ANSYS a multipurpose analysis tool for varied engineering disciplines. The ANSYS program has been in commercial use for over 50 years and has been used extensively in the aerospace, automotive, construction, electronic, energy services, manufacturing, nuclear, plastics, oil, and steel industries. Additionally, many consulting firms and hundreds of universities have used ANSYS for analysis, research, and educational purposes. ANSYS is recognized worldwide as one of the most widely used and capable programs of its type. Ansys design analysis software is the first created within a quality system with ISO 9001 certification, the internationally accepted quality standard. Product development, testing, maintenance and support processes also meet the United States Nuclear Regulatory Commission's quality requirements, as they have for nearly four decades. The Quality Assurance Service Agreement is suitable for the customers working in the nuclear industry who need to meet specific federal regulations including 10CRF50 Appendix B and provisions of 10CFR21. ANSYS has retained its original International Organization for Standardization (ISO) 9001 accreditation certificate since1995-05-04, It’s current certificate is valid until 2027-05-29.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Multifrequency eddy-current detection of fast transient thermal signatures for in situ monitoring applications

Eddy-current (EC) nondestructive evaluation has a long history of use in a variety of ex situ defect monitoring applications because of its exquisite sensitivity to local material variations. Due to the relationship between a material's conductivity and its temperature, EC methods have also been used to investigate quasistatic, long-range temperature variations in casting applications. However, these techniques remain underutilized for the measurement of rapidly varying, spatially nonuniform temperature distributions. In this work, we construct a model system capable of generating repeatable temperature transients in steel plates and measure real-time eddy-current signals with millisecond time resolution and spatial resolution of the order of 1 mm. Using a combination of Multiphysics simulations, fast thermal imaging, and time-resolved holographic interferometry, we tease apart contributions to the eddy-current signals arising from temperature variations and transient plate deformation. Finally, we perform a systematic study in which we vary the plate thickness and the eddy-current excitation frequency to demonstrate that eddy-current techniques can provide information about a three-dimensional, time-varying, subsurface thermal distribution, which is inaccessible to the traditional thermal imaging techniques.

Rosenberg, Ethan R. [Lawrence Livermore National L↗
DOE OSTI2025

A Two-Stage Approach for PV Inverter Engagement in Power Factor Correction and Voltage Regulation

The rapid integration of distributed energy resources, like solar photovoltaics (PVs), can lead to overvolt-age challenges due to reverse power flow and a noticeable decrease in power factor at the substation interface. While existing literature extensively explores utilizing smart inverter capabilities for reactive power flexibility using a volt-var curve (VVC), obtaining time-varying operating points of such curves in real-time is challenging due to computational demands and communication requirements. Similarly, employing optimization-based approaches for reactive power control and active voltage regulation in large-scale distribution feeders is difficult due to the complexity of the problem and the challenges in effectively engaging customer-owned resources. This paper proposes a two-stage strategy to harness smart inverters for reactive power support. The first stage formulates short-term planning by optimally designing VVCs (on a daily or hourly basis) for large-scale solar PVs based on projected system needs and communicating optimal curves to smart inverters in advance. Subsequently, the second stage employs a transactive-based method to involve customer-owned PVs for reactive power support, effectively enhancing overall system performance and addressing real-time demands. In conclusion, the efficacy of this approach will be demonstrated using real-world distribution circuits provided by Vermont Electric Power Company (VELCO) and Vermont Electric Cooperative (VEC).

Poudel, Shiva [Pacific Northwest National Laborato↗
DOE OSTI2025

Cooperative On-Ramp Merging with Time-Varying Vehicle-to-Vehicle Communication Delay Compensation via a Model-Free Approach

Cooperative merging strategies enabled by vehicle-to-vehicle (V2V) communication have shown promise in addressing congestion, fuel inefficiency, and collision risks. However, their performance can be severely degraded by time-varying and uncertain communication delays-an issue often overlooked in existing research, which primarily focuses on merging sequence determination and trajectory planning. Furthermore, practical considerations such as heterogeneous vehicle dynamics, varying road conditions, and real-time implementation complexities are frequently neglected. This paper presents a model-free, online planning framework for cooperative on-ramp merging of connected and automated vehicles (CAVs), explicitly accounting for time-varying V2V communication delays. Without relying on detailed vehicle dynamics, the proposed method introduces a data-driven delay compensation scheme. A co-simulation platform integrating high-fidelity vehicle dynamics, traffic simulation (SUMO), and V2V communication within MATLAB/Simulink is developed to evaluate the proposed method. Simulation results demonstrate that unaddressed V2V communication delays significantly impair merging performance. In contrast, the proposed framework enhances intervehicle distance tracking and maintains low CO2 emissions and fuel consumption, under communication delay across different communication frequencies. In conclusion, its lightweight design also facilitates real-time implementation, making it well-suited for deployment in practical CAV systems.

Accounting↗
DOE OSTI2024

Implementing Superresolution of Nonstationary Tides with Wavelets: An Introduction to CWT_Multi

Abstract Tides are often nonstationary due to nonastronomical influences. Investigating variable tidal properties implies a trade-off between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short analysis windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a MATLAB code that 1) uses CWT linearity (via the “response coefficient method”) to implement superresolution, i.e., resolving tidal constituents beyond the Rayleigh criterion; 2) provides a Munk–Hasselmann constituent selection criterion appropriate for superresolution; and 3) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time scales of days, and multiple constituents per species with fortnightly filters. It outputs astronomical phase lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s), and reconstruction(s); confidence limits; and signal-to-noise ratios. Artificial data and water levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K 1 , S 1 , and P 1 ) on fortnightly time scales, quantifying impacts of diel power peaking (with a 24-h period, like S 1 ) on water levels in the LCRE. CWT_Multi also helps quantify the impacts of high flows and a salt barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to nonstationary data. Significance Statement Ocean tides, especially in coastal and estuarine systems, are often nonstationary, in the sense that the mean and standard deviation of tidal properties vary over time, usually in response to some nontidal process. We introduce here a MATLAB code, CWT_Multi, that uses wavelet transforms to resolve both tidal species and constituents on time scales from a few days to months. Our code accommodates multiple scalar time series and has typical tidal analysis features like constituent selection and inference, plus two forms of uncertainty analyses. It is flexible, allowing the user to adapt analysis properties to diverse datasets. CWT_Multi is applicable to many problems involving time-variable tides, including sea level rise, compound flooding, sediment transport, and wetland habitat analyses. Application to vector data is a straightforward extension, but further development of our uncertainty analysis is merited. Because nonstationary tidal analysis is rapidly advancing, we also define the features of a “well-formed” analysis code.

Lobo, Matthew↗
DOE OSTI2024

A Roadmap for a Lightning Modeling Grand Challenge

This document is a roadmap for building an interconnected model of the physical processes that produce a lightning discharge, and its observable optical and radio signals. We call this a Lightning Modeling Grand Challenge, recognizing that significant effort and coordination of human and financial resources is required to realize the capability. The roadmap serves to outline the coordination of resources necessary to enable stitching together existing knowledge and model components to make a lightning prediction, and to test these predictions with observations. Such a capability does not currently exist. The roadmap is motivated not only by a spirit of scientific inquiry, but by practical challenges faced by US Federal and societal stakeholders. Advancements in lightning observations have outpaced our tests of integrated understanding, leaving many stakeholders unsure how to design their missions to properly detect and discriminate lightning, and unsure how to apply the sometimes-disagreeing lightning signals from diverse instruments. The time is right to connect existing theories and models to support stakeholders in understanding the signals they observe, for needs as diverse as climate monitoring, national security, weather forecasting, public safety, and protection of natural and built environments. The roadmap’s two main technical sections describe the components of a linked physical model, followed by a description of models of lightning signals and sensors that are driven by outputs from the physical model. The goal is to predict the time-varying physical properties of lightning that are self-consistent with the thunderstorm’s structure and dynamics. These lightning signals then propagate through the storm, with realistic dispersion and attenuation, to receivers on the ground or in space. At a high level, the model begins with weather (cloud) model output, including explicit prediction of the electrification of cloud particles. The cloud’s electrical structure drives a model of lightning physics, from initiation, through channel development, and discharges along those channels. Key lightning parameters, such as the temperature and currents in the channel, and their space and time distribution, are then used to produce optical and electromagnetic signal sources that propagate to modeled receivers. This architecture therefore generates a dataset suitable for comparison to existing and envisioned observing systems. The need for additional measurements and field campaigns to support model development is described. In each model sub-component, inputs, outputs, uncertainties, evaluation methods, and next steps are summarized, interleaved with references to the scientific literature. Identifying boundaries between the model sub-components aids in segmenting an integrated, complex model into practical work packages and system sub-components, allowing a diverse team to contribute and maintain the system. We estimate that at least five years of effort and a $\$$10M initial investment is necessary to make a significant step forward. Mechanisms to facilitate community coordination, including annual workshops and open-source code repositories, are described.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Delivery Ring Spill Characterization and Impulse Study

High-intensity particle physics experiments require uniform beam extraction to prevent instantaneous rate spikes from overwhelming detector systems. By analyzing accelerator parameters and extracted beam dynamics, we directly inform spill regulation systems that make real-time adjustments to minimize non-uniformity. This Department of Energy Visiting Faculty Program project transitioned from characterizing Main Injector half-integer slow extraction for SpinQuest to Delivery Ring third-integer slow extraction for Mu2e. Working alongside the Fast Adaptive Neural Control (FANC) group, we developed an automated pipeline that aligns asynchronous instrument channels, embeds quality metrics, and isolates clean spill populations. Analyzing baseline spills alongside a dedicated quadrupole impulse study allowed us to quantify noise structures while mapping time-varying beam response and transit-delay dynamics. These empirical measurements directly ground digital twin models, supporting FANC’s deployment of real-time, FPGA-based neural network controllers in the Mu2e Spill Regulation System.

Dolen, James William [Purdue U., West Lafayette] (↗