Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial Temporal Model”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Policy-Driven Sustainable Saline Drainage Disposal and Forage Production in the Western San Joaquin Valley of California

Environmental policies to address water quality impairments in the San Joaquin River of California have focused on the reduction of salinity and selenium-contaminated subsurface agricultural drainage loads from westside sources. On 31 December 2019, all of the agricultural drainage from a 44,000 ha subarea on the western side of the San Joaquin River basin was curtailed. This policy requires the on-site disposal of all of the agricultural drainage water in perpetuity, except during flooding events, when emergency drainage to the River is sanctioned. The reuse of this saline agricultural drainage water to irrigate forage crops, such as ‘Jose’ tall wheatgrass and alfalfa, in a 2428 ha reuse facility provides an economic return on this pollutant disposal option. Irrigation with brackish water requires careful management to prevent salt accumulation in the crop root zone, which can impact forage yields. The objective of this study was to optimize the sustainability of this reuse facility by maximizing the evaporation potential while achieving cost recovery. This was achieved by assessing the spatial and temporal distribution of the root zone salinity in selected fields of ‘Jose’ tall wheatgrass and alfalfa in the drainage reuse facility, some of which have been irrigated with brackish subsurface drainage water for over fifteen years. Electromagnetic soil surveys using an EM-38 instrument were used to measure the spatial variability of the salinity in the soil profile. The tall wheatgrass fields were irrigated with higher salinity water (1.2–9.3 dS m -1 ) compared to the fields of alfalfa (0.5–6.5 dS m -1 ). Correspondingly, the soil salinity in the tall wheatgrass fields was higher (12.5 dS m -1 –19.3 dS m -1 ) compared to the alfalfa fields (8.97 dS mm -1 –14.4 dS mm -1 ) for the years 2016 and 2017. Better leaching of salts was observed in the fields with a subsurface drainage system installed (13–1 and 13–2). The depth-averaged root zone salinity data sets are being used for the calibration of the transient hydro-salinity computer model CSUID-ID (a one-dimensional version of the Colorado State University Irrigation Drainage Model). This user-friendly decision support tool currently provides a useful framework for the data collection needed to make credible, field-scale salinity budgets. In time, it will provide guidance for appropriate leaching requirements and potential blending decisions for sustainable forage production. This paper shows the tie between environmental drainage policy and the role of local governance in the development of sustainable irrigation practices, and how well-directed collaborative field research can guide future resource management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

2D-imaging of absolute OH and H 2 O 2 profiles in a He–H 2 O nanosecond pulsed dielectric barrier discharge by photo-fragmentation laser-induced fluorescence

We report pulsed dielectric barrier discharges (DBD) in He–H 2 O and He–H 2 O–O 2 mixtures are studied in near atmospheric conditions using temporally and spatially resolved quantitative 2D imaging of the hydroxyl radical (OH) and hydrogen peroxide (H 2 O 2 ). The primary goal was to detect and quantify the production of these strongly oxidative species in water-laden helium discharges in a DBD jet configuration, which is of interest for biomedical applications such as disinfection of surfaces and treatment of biological samples. Hydroxyl profiles are obtained by laser-induced fluorescence (LIF) measurements using 282 nm laser excitation. Hydrogen peroxide profiles are measured by photo-fragmentation LIF (PF-LIF), which involves photo-dissociating H 2 O 2 into OH with a 212.8 nm laser sheet and detecting the OH fragments by LIF. The H 2 O 2 profiles are calibrated by measuring PF-LIF profiles in a reference mixture of He seeded with a known amount of H 2 O 2 . OH profiles are calibrated by measuring OH-radical decay times and comparing these with predictions from a chemical kinetics model. Two different burst discharge modes with five and ten pulses per burst are studied, both with a burst repetition rate of 50 Hz. In both cases, dynamics of OH and H 2 O 2 distributions in the afterglow of the discharge are investigated. Gas temperatures determined from the OH-LIF spectra indicate that gas heating due to the plasma is insignificant. The addition of 5% O 2 in the He admixture decreases the OH densities and increases the H 2 O 2 densities. The increased coupled energy in the ten-pulse discharge increases OH and H 2 O 2 mole fractions, except for the H 2 O 2 in the He–H 2 O–O 2 mixture which is relatively insensitive to the additional pulses.

hydrogen peroxide↗
DOE OSTI2025

Self-Supervised T-GCN for Detection of Disturbance and Propagation in Power Grid

Urban power systems increasingly rely on dense sensing to monitor grid reliability, yet disturbance labels are scarce and events are rare. We present a self-supervised spatio-temporal method that detects, localizes, and characterizes grid frequency disturbances across urban areas using only unlabeled data. Our approach trains a tiny Temporal Graph Convolutional Network (T-GCN) to forecast per-site frequency residuals (deviation from 60 Hz). The sensor graph is constructed directly from signals using pre-event Pearson correlation with a cross-correlation lag penalty without geocoding. At inference, node-level anomalies are the model's forecast errors; region-level alarms arise from connected components of high-score nodes. We estimate disturbance propagation by computing per-node arrival times (first persistent exceedance), then fit a planar or time-of-arrival model to obtain direction, speed, and an epicenter proxy. With only three real events collected at decisecond resolution across U.S. cities, we evaluate the T-GCN and report time-to-detect, footprint size, and propagation consistency. We further show that short-window embeddings from the T-GCN's hidden states enable few-shot event-vs-background recognition via a simple prototypical classifier. Despite minimal data and no labels, our system yields fast, spatially coherent detection and interpretable propagation maps, offering a practical, lightweight pathway to city-scale grid resilience analytics.

Niu, Haoran [ORNL] (ORCID:0000000155228297)↗
DOE OSTI2023

Properties and Controlling Processes of Aerosol and Cloud Condensation Nuclei in Marine Boundary Layer Over Eastern North Atlantic (Final Report)

Remote marine low cloud systems have large spatial and temporal coverages. Because of their relatively low optical thickness and background aerosol concentrations, marine low clouds are particularly susceptible to perturbations in aerosols associated with anthropogenic emissions. Recent studies show that a large fraction of the global aerosol indirect forcing can be attributed to changes in marine low clouds, despite their relatively long distance from most anthropogenic sources. Currently, there remain large uncertainties in the magnitude of the global aerosol indirect radiative forcing. One major contribution to this large uncertainty derives from poor understanding of the aerosols under natural conditions, the perturbation by anthropogenic emissions, and the processes that drive them. In this project, we took advantage of the comprehensive, multi-dimensional characterization of meteorological parameters, trace gases, aerosol, cloud, and drizzle fields during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) and the long-term measurements at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site on Graciosa Island. By synergistically combining in situ measurements, remote sensing retrievals, and model simulations, we quantified aerosol properties in the ENA, their vertical profiles, and variations with season, synoptic state, and airmass. We investigated key aerosol processes, including entrainment of free troposphere aerosol into marine boundary layer, new particle formation, growth of Aitken mode particles, and the influences of these processes on the cloud condensation nuclei (CCN) population in marine boundary layer. The representation of key aerosol processes in the Energy Exascale Earth System Model (E3SM) was evaluated using the observations. This project has led to 26 publications, which are listed at the end of this report.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Probing nonlinear velocity profiles of shear-thinning, nematic platelet dispersions in Couette flow using x-ray photon correlation spectroscopy

In this study, we report experiments employing x-ray photon correlation spectroscopy (XPCS) to characterize the velocity profiles of complex fluids in Couette flow. The approach involves modeling the XPCS correlation functions obtained with the incident x-ray beam passing tangentially through the Couette cell gap at various distances from the inner wall. We first demonstrate the technique with measurements on a dilute colloidal dispersion in the Newtonian liquid glycerol, where the expected linear velocity profiles are recovered. We then employ the technique to map the shear-rate-dependent velocity profiles of a shear-thinning dispersion of nematically ordered Gibbsite platelets. The nonlinear velocity profiles of the Gibbsite dispersion include a narrow slip region adjacent to the outer wall and a band with a small velocity gradient in the interior of the gap that evolves into a region increasingly resembling plug flow with increasing shear rate. Variations in the velocity profile along the vorticity direction indicate an instability in the interface between this region of small velocity gradient and a region of high velocity gradient near the inner wall. The analysis of the small-angle scattering patterns provides information about the spatial and temporal variations in the nematic order of the Gibbsite dispersion and their coupling to the velocity profile. Additional potential applications of this XPCS-based technique and comparisons with established methods for characterizing velocity profiles are discussed.

42 ENGINEERING↗
DOE OSTI2025

Hestia-SWIFL: hourly anthropogenic fossil fuel CO2 and heat on the 2km WRF grid, version 1.1

The Hestia-SWIFL version 1.1 anthropogenic heat (AH) and fossil fuel CO2 (FFCO2) emissions data product represent emissions due to the combustion of fossil fuel and cement production within the state of Arizona from 2019 to 2022. This product was developed as part of the Southwest Urban Corridor Integrated Field Laboratory (SW-IFL) project, which aims to provide new knowledge and tools that address extreme heat, air quality, climate change and related urban environmental issues by integrating high-resolution observations, modeling, and civic engagement. The emissions are generated using a bottom-up/engineering approach and are tied to results generated by the Vulcan Project version 4, an effort to quantify space/time-resolved FFCO2 & AH emissions for the entire United States landscape. A large number of data sources are combined to best estimate the emissions at fine scales such as air quality emissions data, traffic flow data, building information, sociodemographic information, and fuel statistics. The AH product provides emissions for two emissions sources (transportation and point source emissions) in units of Watts per hour per square meter (W/m2) per year (annual files) or per hour (hourly files). The FFCO2 product provides emissions from nine individual emission sectors as well as the total, and in units of tons of carbon (tC) per grid cell per year or per hour. The output made available here places the native spatial resolution of the Hestia FFCO2 & AH emissions data product (points, lines, and polygons) into a regularized 2km x 2km grid at hourly and annual temporal resolutions, and stored in netCDF files. The exact spatial extent is defined by the ASU Weather Research Forecast (WRF) simulation grid. All data are processed using R/Python pm high-performance computing system. 2-27-2026 updates: Bugs in airport hourly profile (both AH and FFCO2) and building spatial patterns (FFCO2 only) were fixed. Hourly emissions are reprocessed for all years to reflect those changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

GPM IMERG V07B and V06B: Evaluation Using Ground-Based Radar Observations and Application in Global Mesoscale Convective System Tracking

This study evaluates the latest Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals for GPM (IMERG V07B) against its predecessor V06B, for studying mesoscale convective systems (MCSs). Both versions are compared using ground-based radar and rain gauge data from five meteorologically diverse regions: the contiguous United States (including eastern coastlines), Amazon rainforest, central Argentina mountains, equatorial Indian Ocean, and northern Australia across multiple temporal (0.5–6 hours) and spatial scales (0.1°–0.25°). An updated global MCS tracking dataset is developed by integrating satellite-observed infrared brightness temperature with IMERG V07B. Comparation of IMERG against radar observations reveals that IMERG demonstrates better performance in capturing the probability distribution and quantitative contributions of rainfall (from no-rain to intense-rain conditions) over tropical oceans than over land, with marked improvements in IMERG V07B for heavy-to-intense rain (> 10 mm h-1). Over land, systematic biases persist: IMERG tends to overestimate light-to-moderate rain (1–10 mm h-1) while underestimating heavy-to-intense rain. Additionally, aggregating IMERG to coarser resolutions (3-hourly or 0.25°) improves consistency with radar observations, outperforming the 1-hourly/0.1° resolution. The new IMERG V07B-based global MCS dataset exhibits consistent statistical characteristics with the V06B-based dataset, despite lower mean rain rates and reduced heavy precipitation contributions. These findings offer valuable insights for utilizing IMERG V07B in global precipitation studies, MCS characterization, and model evaluation.

Zhang, Sihan↗
DOE OSTI2024

Identification of Climatological Representative Days in the Mid-Atlantic for High-Fidelity Offshore Wind Energy Modeling

The goal of reaching 30 GW of offshore wind energy by 2030 becomes more realistic with the continued approval of offshore wind energy areas by the Biden Administration. In the Mid-Atlantic, where wind energy projects are in the most advanced stages of development, there is increased research focus on the eventual interaction of these wind farms. These interactions, in the form of wakes and cluster wakes, or wakes from multiple wind farms, could have detrimental effects on power production and forecastability for downwind wind farms (Pryor et al. 2022, Golbazi et al. 2022, Rosencrans et al. 2023). To help alleviate these issues, numerical simulations in the form of numerical weather prediction (NWP) and large eddy simulations (LES) can provide insight into when cluster wake situations may occur, but running such simulations can be expensive and difficult to run for multiple years. In this study, we leverage and build upon existing techniques in the literature (Fischereit et al. 2022) to identify climatologically representative days for wind energy areas in the Mid-Atlantic where conditions would promote cluster wake situations. We select meteorological variables (wind speed, wind direction, atmospheric stability, boundary-layer height, TKE) critical to understanding wind energy production and wake propagation. We then consider two different NWP datasets of varying spatial and temporal resolution: ERA5 provides data at hourly intervals from 1940 to present at 0.25 deg (31 km) spatial resolution (Hersbach et al. 2020), and the NOW-23 dataset provides data at 5-minute resolution for 21 years at 2-km spatial resolution (Bodini et al. 2020). Our first step is to compare these two datasets for an overlapping 21-year time period. Initial results show that the required number of days to represent the long-term climate increases with each additional variable considered. In their study of the German Bight, Fischereit et al. (2022) found that they could represent the long-term wind and wave climate in a "near-perfect" way with -180 days, by reaching a Perkins Skill Score (PSS) of 0.9; our investigation of the mid-Atlantic wind resource region with ERA5 and NOW-23 data suggests that we will need -100 days to reach a PSS of 0.9. As we expand our parameter space to include multiple variables, the number of required days will likely grow. These results will ultimately be used to select case studies to best represent cluster wake conditions that apply to this region for the lifetime of likely wind farms in this mid-Atlantic region.

clusterwakes↗
DOE OSTI2022

Lighting application efficacy: A framework for holistically measuring lighting use in buildings

Lighting consumes significant energy in buildings, but current measurements of lighting efficiency are inadequate in capturing the spatial and temporal effectiveness of light. The most widely used lighting energy efficiency measure, luminous efficacy, quantifies the amount of light generated by individual light sources. While luminous efficacy provides manufacturers and designers with a useful metric for quantifying the energy efficiency of individual lighting devices, it does not fully address the use of lighting in architectural spaces. Instead, application efficacy, the relationship between the electrical power consumed by lighting hardware and the amount of light that contributes to the visual perception of building occupants, must be measured. Here, we propose a framework for calculating lighting application efficacy (LAE) based on the primary pathway of light in architectural settings: the generation and emission of light from a luminaire, the travel of the light throughout the space and into occupants’ eyes, and the process of visual perception. The LAE framework will include the development of computational methods for estimating the proportion of emitted light that is directed to task areas and areas corresponding to occupants’ visual fields using ray-tracing lighting design software. Future modeling and validation work will quantify energy efficiency of buildings.

Durmus, Dorukalp↗
DOE OSTI2021

Novel 2D velocity estimation method for large transient events in plasmas

Dynamics of fast transient events are challenging to be analyzed with high time resolution. Such events can occur in fusion plasmas such as the filaments during edge-localized modes (ELMs). Here, we present a robust method—the spatial displacement estimation—for estimating the displacements of structures with fast dynamics from high spatial and time resolution imaging diagnostics [e.g., gas-puff imaging (GPI)] with sampling time temporal resolution. First, a background suppression method is shown, which suppresses the slowly time-evolving and spatially non-uniform background in the signal. In the second step, a two-dimensional polynomial trend subtraction method is presented to tackle the remaining polynomial order trend in the signal. After performing these pre-processing steps, the spatial displacement of the propagating structure is estimated from the two-dimensional spatial cross-correlation coefficient function calculated between consecutive frames. The method is tested for its robustness and accuracy by simulated Gaussian events and spatially displaced random noise. An example application of the method is presented on propagating ELM filaments measured by the GPI system on the National Spherical Torus Experiment spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Variational AutoEncoders Reveal Intensifying GPP Extremes in Continental United States based on CESM2 Simulations

Climate extremes significantly impact terrestrial carbon cycle dynamics, necessitating robust methods for detecting and analyzing anomalous behavior in plant productivity. This study presents a novel application of variational autoencoders (VAE) for identifying extreme events in gross primary productivity (GPP) from Community Earth System Model version 2 simulations across four AR6 regions in the Continental United States. We compare VAE-based anomaly detection with traditional singular spectral analysis (SSA) methods across three time periods: 1850-80, 1950-80, and 2050-80 under SSP5-8.5 scenario. The VAE architecture employs three dense layers and a latent space with input sequence length of 12 months, training on normalized GPP time series to reconstruct the GPP and identify anomalies based on reconstruction errors. Extreme events are defined using 5th percentile thresholds applied to both VAE and SSA anomalies. Results demonstrate strong regional agreement between VAE and SSA methods in spatial patterns of extreme event frequencies, despite VAE consistently producing higher threshold values (179-756 GgC for VAE vs. 100-784 GgC for SSA across regions and periods). Both methods reveal increasing magnitudes and frequencies of negative carbon cycle extremes toward 2050-80, particularly in Western and Central North America. The VAE approach shows comparable performance to established SSA techniques while offering computational advantages and enhanced capability for capturing non-linear temporal dependencies in carbon cycle variability. This research demonstrates the potential of deep learning approaches for extremes detection and provides a foundation for improved understanding of future carbon cycle risks under future conditions.

Sharma, Bharat [ORNL] (ORCID:0000000266982487)↗
DOE OSTI2022

In situ melt pool measurements for laser powder bed fusion using multi sensing and correlation analysis

Laser powder bed fusion is a promising technology for local deposition and microstructure control, but it suffers from defects such as delamination and porosity due to the lack of understanding of melt pool dynamics. To study the fundamental behavior of the melt pool, both geometric and thermal sensing with high spatial and temporal resolutions are necessary. This work applies and integrates three advanced sensing technologies: synchrotron X-ray imaging, high-speed IR camera, and high-spatial-resolution IR camera to characterize the evolution of the melt pool shape, keyhole, vapor plume, and thermal evolution in Ti–6Al–4V and 410 stainless steel spot melt cases. Aside from presenting the sensing capability, this paper develops an effective algorithm for high-speed X-ray imaging data to identify melt pool geometries accurately. Preprocessing methods are also implemented for the IR data to estimate the emissivity value and extrapolate the saturated pixels. Quantifications on boundary velocities, melt pool dimensions, thermal gradients, and cooling rates are performed, enabling future comprehensive melt pool dynamics and microstructure analysis. The study discovers a strong correlation between the thermal and X-ray data, demonstrating the feasibility of using relatively cheap IR cameras to predict features that currently can only be captured using costly synchrotron X-ray imaging. Such correlation can be used for future thermal-based melt pool control and model validation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

deadtrees.earth — An open-access and interactive database for centimeter-scale aerial imagery to uncover global tree mortality dynamics

Excessive tree mortality is a global concern and remains poorly understood as it is a complex phenomenon. We lack global and temporally continuous coverage on tree mortality data. Ground-based observations on tree mortality, e.g., derived from national inventories, are very sparse, and may not be standardized or spatially explicit. Earth observation data, combined with supervised machine learning, offer a promising approach to map overstory tree mortality in a consistent manner over space and time. However, global-scale machine learning requires broad training data covering a wide range of environmental settings and forest types. Low altitude observation platforms (e.g., drones or airplanes) provide a cost-effective source of training data by capturing high-resolution orthophotos of overstory tree mortality events at centimeter-scale resolution. Here, we introduce deadtrees.earth, an open-access platform hosting more than two thousand centimeter-resolution orthophotos, covering more than 1,000,000 ha, of which more than 58,000 ha are manually annotated with live/dead tree classifications. This community-sourced and rigorously curated dataset can serve as a comprehensive reference dataset to uncover tree mortality patterns from local to global scales using space-based Earth observation data and machine learning models. This will provide the basis to attribute tree mortality patterns to environmental changes or project tree mortality dynamics to the future. The open nature of deadtrees.earth, together with its curation of high-quality, spatially representative, and ecologically diverse data will continuously increase our capacity to uncover and understand tree mortality dynamics.

Citizen science↗
DOE OSTI2024

Detecting dark matter substructures on small scales with fast radio bursts

The matter power spectrum is only weakly constrained on subgalactic scales, while physics beyond the Standard Model can leave unique imprints, especially on sub-parsec scales. We propose measuring the arrival-time difference of fast radio bursts (FRBs) along two adjacent sightlines as a new probe to dark matter substructures on scales down to ∼ 1 AU . We discuss two observational scenarios in which it may be possible to place interesting constraints on such models through the monitoring of repeating FRB sources: (i) By sending radio receivers to space to form a baseline of tens of AU or more and measuring the temporal variation of the arrival-time difference between receivers. (ii) By measuring the temporal variation of the arrival-time difference between two lensed images of one strongly lensed repeater. In both scenarios, obtaining interesting constraints requires correlating the voltage time series to measure the radio signal arrival time to sub-nanosecond precision. We find that two radio dishes separated by 20 AU may be sensitive to the enhancement of small-scale structures at ∼ 10 - 8 M ⊙ masses in the QCD axion dark matter scenario, or from an early epoch of matter domination with a reheating temperature up to 60 MeV. Other dark matter models, such as those composed of ∼ 10 - 13 M ⊙ primordial black holes produced during inflation, would also be probed by this method. We further show that a strong lensing situation of multiple images provides an equivalent ∼ 2000 AU ( σ v / 10 3 km s - 1 ) ( δ t / 10 yr ) baseline, for a typical velocity of dark matter substructures σ v and an observational time span δ . This is much more sensitive, but with the uncertainty that intervening decoherence from the interstellar medium may degrade the timing precision, and that spatial variation in the FRB emission spot may result in confounding signals. We show that the lensing magnifications of Type Ia supernovae constrain a similar quantity to such FRB timing, with present limits being equivalent to ruling out the same parameter space that would be probed by a 0.14 AU baseline.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spatio-temporal electric field distributions in an atmospheric plasma jet impinging on a microchannel array surface

The electric field distribution in the ionization waves (IWs) propagating over a microchannel array dielectric surface, with the channels either empty or filled with distilled water, is measured by ps electric field induced second harmonic generation. The surface IW is initiated by the atmospheric pressure N 2 –Ar plasma jet impinging on the surface vertically and powered by ns pulse discharge bursts. The results show that the electric field inside the microchannels, specifically its horizontal component, is enhanced by up to a factor of 2. The field enhancement region is localized within the channels. The vertical electric field inside the channels lags in time compared to the field measured at the ridges, indicating the transient reversal of the IW propagation direction across the channels (toward the jet). This is consistent with the phase-locked plasma emission images and confirmed by the kinetic modeling predictions, which show that the IW ‘jumps’ over the empty channels and propagates into the channels only after the jump between the adjacent ridges. When the channels are filled with water, the wave speed increases by up to 50%, due to the higher effective dielectric constant of the surface. No evidence of a significant electric field enhancement near the dielectric surface (ceramic or water) has been detected, within the spatial resolution of the present diagnostic, ~100 μ m.

plasmas and liquids↗
DOE OSTI2022

Wave resource spatial and temporal variability dependence on WEC size

As the wave energy sector grows and looks to the Blue Economy for commercialization opportunities, there is a distinct and pressing need to clearly understand and quantify the coupled impacts of wave energy converter (WEC) size and wave resource characteristics on the annual energy production, spatial variability and temporal variability. Utilizing generic frequency domain representations of the Oscilla Power Triton WEC and spectral wave conditions at PacWave (Oregon), Los Angeles (California) and WETS (Hawaii), a series of interesting results emerge. Firstly, the ‘optimal’ WEC size, from an energy standpoint, is fundamentally dependent on the frequency distribution of the incoming wave variance density spectrum. Secondly, and from a seasonality perspective, the seasonal WEC energy generation doesn’t necessarily follow the seasonal distribution of gross wave power. Finally, from an hourly power variability perspective, a reduction in WEC size generally decreases variability. However, for each of the locations investigated, there appears to be a WEC size threshold; a threshold where further reducing WEC size results in increased power variability.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Dynamics of filaments during the edge-localized mode crash on NSTX

Edge localized modes (ELMs) are routinely observed in H-mode plasma regimes of the National Spherical Torus Experiment (NSTX). Due to the explosive nature of the instability, only diagnostics with high temporal and spatial resolution could provide a detailed insight into the dynamics associated with the ELMs. Gas-puff imaging (GPI) at NSTX provides 2D measurements of the magnetic field aligned fluctuations (e.g. ELM filaments) in the scrape-off layer and the at the plasma edge with 2.5 us temporal and 10 mm optical resolution.A novel analysis technique was developed to estimate the frame-by-frame velocities and the spatial parameters of the dominant structures associated with the ELMs. The analysis was applied to single ELM events to characterize the ELM crash dynamics, and then extended to a database of 169 ELM events.Statistical analysis was performed in order to find the characterizing dynamics of the ELM crash. The results show that on average an ELM crash consists of a filament with a circular cross-section which is propelled outwards with a characterizing peak radial velocity of ~3.3 km/s. The radial velocity was found to be linearly dependent on the distance of the filament from the separatrix, which has never been seen before. The ELM filament is characterized by propagation in the ion-diamagnetic direction poloidally with a peak velocity of 11.4 km/s. The ELM crash lasts for approximately 100us until the radial propulsion settles back to the pre-ELM level. The experimental findings were compared with analytical theory. Two possible mechanisms were identified for explaining the observations: the curvature interchange model and the current-filament interaction model.

correlation↗