Search NASA⌕ Search

SEARCH · Search NASA

Results for “array modeling”

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 397 records · Page 22

A [FeFe] Hydrogenase–Rubrerythrin Chimeric Enzyme Functions to Couple H 2 Oxidation to Reduction of H 2 O 2 in the Foodborne Pathogen Clostridium perfringens

[FeFe] hydrogenases are a diverse class of H 2 -activating enzymes with a wide range of utilities in nature. As H 2 is a promising renewable energy carrier, exploration of the increasingly realized functional diversity of [FeFe] hydrogenases is instrumental for understanding how these remarkable enzymes can benefit society and inspire new technologies. In this work, we uncover the properties of a highly unusual natural chimera composed of a [FeFe] hydrogenase and rubrerythrin as a single polypeptide. The unique combination of [FeFe] hydrogenase with rubrerythrin, an enzyme that functions in H 2 O 2 detoxification, raises the question of whether catalytic reactions, such as H 2 oxidation and H 2 O 2 reduction, are functionally linked. Herein, we express and purify a representative chimera from Clostridium perfringens (termed Cper HydR) and apply various electrochemical and spectroscopic approaches to determine its activity and confirm the presence of each of the proposed metallocofactors. The cumulative data demonstrate that the enzyme contains a surprising array of metallocofactors: the catalytic site of [FeFe] hydrogenase termed the H-cluster, two [4Fe-4S] clusters, two rubredoxin Fe(Cys) 4 centers, and a hemerythrin-like diiron site. The absence of an H 2 -evolution current in protein film voltammetry highlights an exceptional bias of this enzyme toward H 2 oxidation to the greatest extent that has been observed for a [FeFe] hydrogenase. Here, we demonstrate that Cper HydR uses H 2 , catalytically split by the hydrogenase domain, to reduce H 2 O 2 by the diiron site. Structural modeling suggests a homodimeric nature of the protein. Overall, this study demonstrates that Cper HydR is an H 2 -dependent H 2 O 2 reductase. Equipped with this information, we discuss the possible role of this enzyme as a part of the oxygen-stress response system, proposing that Cper HydR constitutes a new pathway for H 2 O 2 mitigation.

08 HYDROGEN↗

CTA and SWGO can discover Higgsino dark matter annihilation

Thermal Higgsino dark matter (DM), with a mass near 1.1 TeV, is one of the most well-motivated and untested DM candidates. Leveraging recent hydrodynamic cosmological simulations that give DM density profiles in Milky Way analog galaxies we show that the linelike gamma-ray signal predicted from Higgsino annihilation in the Galactic Center could be detected at high significance with the upcoming Cherenkov Telescope Array (CTA) and Southern Wide-field Gamma-ray Observatory (SWGO) for all but the most pessimistic DM profiles. We perform the most sensitive search to-date for the linelike signal using 15 years of data from the Fermi Large Area Telescope, coming within an order one factor of the necessary sensitivity to detect the Higgsino for some Milky Way analog DM density profiles. We show that H.E.S.S. has subleading sensitivity relative to Fermi for the Higgsino at present. In contrast, we analyze H.E.S.S. inner Galaxy data for the thermal wino model with a mass near 2.8 TeV; we find no evidence for a DM signal and exclude the wino by over a factor of two in cross section for all DM profiles considered. In the process, we identify and attempt to correct what appears to be an inconsistency in previous H.E.S.S. inner Galaxy analyses for DM annihilation related to the analysis effective area, which may weaken the DM cross-section sensitivity claimed in those works by around an order of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Solar PV

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and 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) research platform. This dataset represents 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 nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence data centers and other variable loads. This dataset entry describes the behavior of a 1.25-MW proton exchange membrane MC250 electrolyzer system, manufactured by Nel Hydrogen , [1] when fed historical data generated by the 430-kW, fixed-axis solar photovoltaic (PV) array located at NLR’s Flatirons Campus. (While the electrolyzer balance of plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack.) Solar PV power output data for the 2020 calendar year were categorized on a daily basis by total energy generation and standard deviation. Each day was then ranked by these metrics, and the 25th, 50th, and 100th percentiles were selected. The 75th percentile day did not exhibit sufficient variability to make for a valuable experiment. A similar process was used for the related historical wind dataset . [2] The historical days in 2020 that represented these percentiles are Dec. 19, March 29, and May 4, respectively. The entire solar day’s power profile was then fed through the MC250 electrolyzer. Due to its length, the 100th percentile day experiment was split into two parts, and the final 3 hours of the solar day were not captured. These final 3 hours contained no spikes or dips of interest and simply represented a slow decay of input solar power. Also, a single timestamp (13:13:47 on Jan. 14, 2026) was lost in the hydrogen system supervisory control and data acquisition. Finally, during the 25th percentile experiment (solar day Dec. 19, 2020) data recording was lost from 11:00:13 to 11:14:45. The roughly 15 minutes of the solar profile were rerun at the end of the experiment and spliced into this time slot during post-processing. The electrolysis system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operation 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 historical solar profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. For more details on the statistical analysis process, see the slide deck “Public Reference Data for Megawatt-Scale Hydrogen Electrolysis: NLR Historical Solar PV Analysis and Profile Generation” accessible with this data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single solar PV electrolysis experiment and is formatted as: {technology}_{percentile}_{scaling factor} For instance, “solarPV-430kW_25_2x.zip” reports the experiment using the 25th percentile solar data from the historical 2020 solar PV dataset, scaled to 200%. Scaling factors were applied to the generated solar PV power output files to more closely match the 1.25-MW capacity of the electrolyzer. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and solar power input. A PDF file detailing the historical solar data statistical analysis used to generate the solar profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all experiments combined into one dataset labeled "combined_solarPV_experiments.csv". [1] nelhydrogen.com/product/mc-series-electrolyser . [2] data.nlr.gov/submissions/316 .

08 HYDROGEN↗

Relating Hydro–Mechanical and Elastodynamic Properties of Dynamically Stressed Tensile–Fractured Rock in Relation to Applied Normal Stress, Fracture Aperture, and Contact Area

We exploit nonlinear elastodynamic properties of fractured rock to probe the micro-scale mechanics of fractures and understand the relation between fluid transport and fracture aperture under dynamic stressing. Experiments were conducted on rough, tensile-fractured Westerly granite subject to triaxial stresses. We measure fracture permeability for steady-state fluid flow with deionized water. Pore pressure oscillations are applied at amplitudes ranging from 0.2 to 1 MPa at 1 Hz frequency. During dynamic stressing we transmit ultrasonic signals through the fracture using an array of piezoelectric transducers (PZTs) to monitor evolution of interface properties. We examine the influence of fracture aperture and contact area by conducting measurements at effective normal stresses of 10–20 MPa. Additionally, the evolution of contact area with stress is characterized using pressure sensitive film. These experiments are conducted separately with the same fracture and map contact area at stresses from 9 to 21 MPa. The measurements are a proxy for “true” contact area for the fracture surface and we relate them to elastic properties using the calculated PZT sensor footprints via numerical modeling of Fresnel zones. We compare the elastodynamic response of the fracture using the stress-induced changes in ultrasonic wave velocities for transmitter-receiver pairs to image spatial variations in contact properties. We show that nonlinear elasticity and permeability enhancement decrease with increasing normal stress. Additionally, post-oscillation wave velocity and permeability exhibit quick recoveries toward pre-oscillation values. Estimates of fracture contact area (global and local) demonstrate that the elastodynamic and permeability responses are dominated by fracture topology.

58 GEOSCIENCES↗

Genetic tools for engineering Zymomonas mobilis , Cereibacter sphaeroides and Novosphingobium aromaticivorans to improve production of bioenergy compounds

Limited genetic tools for non-model bacteria are one of the limiting factors for genetic studies. This review compiles genetic tools used for three non-model alpha-proteobacteria, such as Zymomonas mobilis, Cereibacter (Rhodobacter) sphaeroides, and Novosphingobium aromaticivorans, which hold significant potential to produce industrially essential bioenergy compounds due to their distinctive metabolic pathways and resilience in extreme environments. Each of these strains has a unique genetic profile that enables them to efficiently carry out key reactions relevant to producing bioenergy compounds, such as converting sugars into bioenergy compounds and breaking down lignotoxins. Genetic tools can further optimize these strains for enhanced bioenergy compound production. This review explores the metabolic advantages of these organisms. It highlights the available array of genetic toolkits that can be shared among them to unlock their full potential for sustainable biofuel production.

Biofuel↗

Effects of the U.S. inflation reduction act on SMR economics

The U.S. Inflation Reduction Act (IRA) of 2022 provides a wide array of tax credits and other incentives for low-carbon energy. The technology-neutral clean generation production tax credit (PTC) (Section 45Y of the U.S. Internal Revenue Code) and the technology-neutral investment tax credit (ITC) (Section 48E) lower the net cost of new electricity generation projects with zero or negative greenhouse gas emission rates. We evaluate the impact of the IRA legislation—specifically the PTC and ITC—on the cost-competitiveness of small modular reactors (SMRs). We use the Argonne Low-carbon Energy Analysis Framework (A-LEAF) model to calculate the capacity factor of an SMR with a range of hypothetical variable operating and maintenance (O&M) costs in the Electric Reliability Council of Texas (ERCOT) electricity market. We selected ERCOT for market modeling because of its competitive structure, available data, and extensive use in prior literature. We use a discounted cash flow model to calculate the SMR’s net present value based on the market prices and capacity factors from A-LEAF, hypothetical ranges of capital and variable O&M costs, and other input parameters, with or without the IRA tax credits. We determine the SMR owner’s optimal choice of PTC or ITC for the hypothetical ranges of capital and variable O&M costs. We also evaluate potential shifts in the SMR owner’s optimal choice of PTC or ITC based on historical patterns of nuclear capital cost overruns in the United States. We also assess the sensitivity of our results to longer PTC period and electricity prices from the New England market, which tend to be higher than electricity prices in ERCOT. We find that even with the IRA tax credits, only SMRs with low capital and variable O&M costs would be economically feasible in the low-price ERCOT market scenario modeled. A longer PTC period and higher-price market such as New England, however, would significantly expand the economic feasibility of SMRs in the United States.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗

Combined dark matter search towards dwarf spheroidal galaxies with Fermi -LAT, HAWC, H.E.S.S., MAGIC, and VERITAS

Dwarf spheroidal galaxies (dSphs) are excellent targets for indirect dark matter (DM) searches using gamma-ray telescopes because they are thought to have high DM content and a low astrophysical background. The sensitivity of these searches is improved by combining the observations of dSphs made by different gamma-ray telescopes. We present the results of a combined search by the most sensitive currently operating gamma-ray telescopes, namely: the satellite-borne Fermi -LAT telescope; the ground-based imaging atmospheric Cherenkov telescope arrays H.E.S.S., MAGIC, and VERITAS; and the HAWC water Cherenkov detector. Individual datasets were analyzed using a common statistical approach. Results were subsequently combined via a global joint likelihood analysis. We obtain constraints on the velocity-weighted cross section 〈σv〉 for DM self-annihilation as a function of the DM particle mass. This five-instrument combination allows the derivation of up to 2-3 times more constraining upper limits on 〈σv〉 than the individual results over a wide mass range spanning from 5 GeV to 100 TeV. Depending on the DM content modeling, the 95% confidence level observed limits reach 1.5×10 -24 cm 3 s -1 and 3.2×10 -25 cm 3 s -1 , respectively, in the τ + τ - annihilation channel for a DM mass of 2 TeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Microwave microscope studies of trapped vortex dynamics in superconductors

Trapped vortices in superconductors introduce residual resistance in superconducting radio-frequency (SRF) cavities and disrupt the operation of superconducting quantum and digital electronic circuits. Understanding the detailed dynamics of trapped vortices under oscillating magnetic fields is essential for advancing these technologies. We have developed a near-field magnetic microwave microscope to study the dynamics of a limited number of trapped vortices under the probe when stimulated by a localized rf magnetic field. By measuring the local second-harmonic response (𝑃 2⁢f ) at subfemto-Watt levels, we isolate signals exclusively arising from trapped vortices, excluding contributions from surface defects and Meissner screening currents. Toy models of niobium superconductor hosting vortex pinning sites are introduced and studied with time-dependent Ginzburg-Landau (TDGL) simulations of probe/sample interaction to better understand the measured second-harmonic response. The simulation results demonstrate that the second-harmonic response of trapped vortex motion under a localized rf magnetic field shares key features with the experimental data. Here, this measurement technique provides access to vortex dynamics at the micrometer scale, such as depinning events and spatially resolved pinning properties, as demonstrated in measurements on a niobium film with an antidot flux pinning array.

43 PARTICLE ACCELERATORS↗

Demonstration of Vertical 2T-nC FeRAM Hybrid Cell and Its Scalability for High-Density 3-D Ferroelectric Capacitor Memory

In this work, we present a comprehensive experimental and modeling study on the scaling of vertical 2T-nC ferroelectric random access memory (FeRAM) hybrid cells, comprising n metal-ferroelectric–metal (MFM) capacitors, to demonstrate a high-performance and high-density 3-D capacitor memory. Our contributions include: 1) successful process integration of vertical 2T-3C FeRAM cells by stacking MFM structures on top of Si CMOS transistors; 2) experimental validation of memory cell functionality, confirming the feasibility of the vertical 2T-nC FeRAM architecture; 3) an analysis of scaling effects on parasitic capacitance in densely integrated 3-D arrays, using 3-D technology computer-aided design (TCAD) simulations; 4) exploration of aggressive stacking of write bitlines (WBLs) to enhance memory density, where ferroelectric linear capacitance ( C FE ) enables self-boosted inhibition under the V W /2 scheme, but renders the V W /3 scheme ineffective due to intolerable write disturbances; and 5) assessment of horizontal scaling, revealing significant increases in read disturbances caused by interplane capacitance between adjacent WBLs ( C Z ). This work represents an early exploration into the potential of 2T-nC FeRAM as a scalable and efficient 3-D memory solution.

42 ENGINEERING↗

Optically accessible high-finesse millimeter-wave resonator for cavity quantum electrodynamics with atom arrays

Cavity quantum electrodynamics (QED) is a powerful tool in quantum science, enabling preparation of nonclassical states of light and scalable entanglement of many atoms coupled to a single field mode. While the most coherent atom-photon interactions have been achieved using superconducting millimeter-wave cavities coupled to Rydberg atoms, these platforms so far lack the optical access required for trapping and addressing individual atomic qubits. We present a millimeter-wave Fabry-Pérot cavity with finesse 5.8⁢(1) ×10 7 at a temperature of 1 K providing generous transverse optical access (numerical aperture 0.56). Conflicting goals of strong atom-photon coupling and optical access motivate a near-confocal geometry. Close to confocality, however, postparaxial corrections to the cavity spectrum introduce unexpected degeneracies between transverse modes, leading to excess cavity loss. Modeling these corrections allows for tuning the cavity geometry to evade this loss, producing a high finesse that will enable cavity QED experiments with trapped atoms deep in the strong coupling regime.

Zhang, Tony [Stanford Univ., CA (United States); S↗

CRAGE-RB-PI-seq reveals transcriptional dynamics of plant-associated bacteria during root colonization

Plant roots release a wide array of metabolites into the rhizosphere, shaping microbial communities and their functions. While metagenomics has expanded our understanding of these communities, little is known about the physiology of their members in host environments. Transcriptome analysis via RNA sequencing is a common approach to learning more, but its use has been challenging because of low bacterial biomass and interference from plant RNA. To overcome this, we developed a randomly-barcoded promoter-library insertion sequencing (RB-PI-seq) combined with chassis-independent recombinase-assisted genome engineering (CRAGE). Using Pseudomonas simiae WCS417 as a model rhizobacterium, this method enabled targeted amplification of barcoded transcripts, bypassing plant RNA interference and allowing measurement of thousands of promoter activities during Arabidopsis root colonization. Our analysis revealed temporally resolved transcriptional regulation, including those associated with cell growth, chemotaxis, plant immune suppression, biofilm formation, and stress responses, reflecting the coordinated physiological adaptation to the root environment. Additionally, we discovered that transcriptional activation of xanthine dehydrogenase and a lysozyme inhibitor is crucial for evading plant immune systems. This framework is scalable to other bacterial species and provides new opportunities for understanding rhizobacterial gene regulation in native environments.

59 BASIC BIOLOGICAL SCIENCES↗

Hopping and crawling DNA-coated colloids

Understanding the motion of particles with multivalent ligand-receptors is important for biomedical applications and material design. Yet, even among a single design, the prototypical DNA-coated colloids, seemingly similar micrometric particles hop or roll, depending on the study. We shed light on this problem by observing DNA-coated colloids diffusing near surfaces coated with complementary strands for a wide array of coating designs. We find colloids rapidly switch between 2 modes: They hop—with long and fast steps—and crawl—with short and slow steps. Both modes occur at all temperatures around the melting point and over various designs. The particles become increasingly subdiffusive as temperature decreases, in line with subsequent velocity steps becoming increasingly anticorrelated, corresponding to switchbacks in the trajectories. Overall, crawling (or hopping) phases are more predominant at low (or high) temperatures; crawling is also more efficient at low temperatures than hopping to cover large distances. We rationalize this behavior within a simple model: At lower temperatures, the number of bound strands increases, and detachment of all bonds is unlikely, hence, hopping is prevented and crawling favored. We thus reveal the mechanism behind a common design rule relying on increased strand density for long-range self-assembly: Dense strands on surfaces are required to enable crawling, possibly facilitating particle rearrangements.

Science & Technology - Other Topics↗

Evolution of oxygen and stratification and their relationship in the North Pacific Ocean in CMIP6 Earth system models

Abstract. This study examines the linkages between the upper-ocean (0–200 m) oxygen (O2) content and stratification in the North Pacific Ocean using four Earth system models (ESMs), an ocean hindcast simulation, and an ocean reanalysis. The trends and variability in oceanic O2 content are driven by the imbalance between physical supply and biological demand. Physical supply is primarily controlled by ocean ventilation, which is responsible for the transport of O2-rich surface waters to the subsurface. Isopycnic potential vorticity (IPV), a quasi-conservative tracer proportional to density stratification that can be evaluated from temperature and salinity measurements, is used herein as a dynamical proxy for ocean ventilation. The predictability potential of the IPV field is evaluated through its information entropy. The results highlight a strong O2–IPV connection and somewhat higher (as compared to the rest of the basin) predictability potential for IPV across the tropical Pacific, where the El Niño–Southern Oscillation occurs. This pattern of higher predictability and strong anticorrelation between O2 and stratification is robust across multiple models and datasets. In contrast, IPV at mid-latitudes has low predictability potential and its center of action differs from that of O2. In addition, the locations of extreme events or hotspots may or may not differ between the two fields, with a strong model dependency, which persists in future projections. On the one hand, these results suggest that it may be possible to monitor ocean O2 in the tropical Pacific based on a few observational sites co-located with the more abundant IPV measurements; on the other, they lead us to question the robustness of the IPV–O2 relationship in the extratropics. The proposed framework helps to characterize and interpret O2 variability in relation to physical variability and may be especially useful in the analysis of new observation-based data products derived from the BGC-Argo float array in combination with the traditional but far more abundant Argo data.

Novi, Lyuba↗

The Surface Chemistry of Methanol on Cu 3 Pd(111): Effects of Metal Alloying and Reaction with Hydrogen

Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS) was used to study the adsorption and surface chemistry of methanol on a Cu 3 Pd(111) model surface. The composition and morphological properties of the pristine Cu 3 Pd(111) substrate were analyzed using a combination of low-energy electron diffraction (LEED), scanning tunneling microscopy (STM), and low-energy ion scattering (LEIS). The results of ion scattering showed segregation of Pd toward the surface, with a Pd/Cu ratio close to 0.5, a larger value than the ratio of 0.33 expected for a Cu 3 Pd bimetallic structure. The surface of the alloy exhibited a good LEED pattern with three-fold long-range periodicities. In STM, clusters of palladium with hexagonal arrays and Pd- Pd distances of 2.7-2.8 Å were detected. Bonding to copper perturbed the electronic properties of the atoms in the Pd clusters, shifting their 4d states toward higher binding energy with respect to the Fermi level. Further, the valence band spectrum of Cu 3 Pd(111) exhibited a line shape that was very different from those displayed by Cu(111) or Pd(111). At low pressures, the adsorption of methanol on Cu 3 Pd(111) at 300 K mainly produced CH 3 O, CO and CH x species. AP-XPS showed that most Pd atoms in the surface of the Cu 3 Pd(111) alloy interacted with the decomposition products of methanol. No significant changes were observed in the core levels of copper upon the adsorption and dissociation of methanol, suggesting that the molecule mainly interacted with Pd sites of the alloy. Reaction with hydrogen led to fast removal of CH x , C and PdC x species from Cu 3 Pd(111) at moderate (< 450 K) temperatures and prevented a CH x → C transformation. If the stability of adsorbed CH 3 O is used as a descriptor for the hydrogenation of CO 2 to methanol, Cu 3 Pd should be a much better catalyst than monometallic palladium. This may be a consequence of electronic and ensemble effects in the alloy that moderate the reactivity of Pd sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Data Science and Machine Learning Platform Supporting Large Particle Accelerator Control and Diagnostics Applications Final Report: SBIR Initial Phase II DE-SC0022583

The Machine Learning Data Platform (MLDP) is a product providing full-stack support for data science, Machine Learning, and Artificial Intelligence (ML/AI) applications at particle accelerator and large experimental physics facilities. It supports ML/AI applications from front-end, high-speed acquisition of heterogeneous, time-series data, through data archiving and management, to back-end analysis. The MLDP embodies a “data-science ready” platform for data analysis and ML/AI applications in diagnosis, modelling, control, and optimization of these facilities. It provides data scientists and applications a consistent, datacentric interface to archive data standardizing implementation and deployment of ML/AI algorithms to different operations configurations within the same facility, or between facilities. Being an open-source, public-domain project, the MLDP is intended for broadest possible impact by increasing accessibility and minimizing the required expertise for installation and operation. The MLDP can also be deployed at user facilities for experimental data collection, archiving, and analysis. It is capable of acquisition and archiving of heterogeneous data from experimental equipment (e.g., images, arrays, structures, etc.) along with system hardware configurations (e.g., scalars, tables), control system process variables, and any metadata required for provenance. Thus, the MLDP can manage experimental data through its entire lifecycle, from acquisition and archiving, through analysis and investigation, to release and final publication.

43 PARTICLE ACCELERATORS↗

Neural network-based classification and regression of magnetohydrodynamic modes in tokamaks

We present a machine learning-based magnetohydrodynamic (MHD) classifier and regressor that utilizes real or complex-valued 3D magnetic sensor array data to determine neoclassical tearing mode (NTM) onset times in tokamaks with millisecond accuracy. The input dataset consists of poloidal profiles of complex Fourier amplitudes with an n = 1 toroidal mode number from 144 human-labeled ITER Baseline Scenario discharges in the DIII-D tokamak, spanning both tearing-dominated and sawtooth-dominated regimes. Since m, n = 2,1 NTMs frequently emerge alongside sawteeth at the same frequency in this scenario, the focus is on isolating the m = 1 and m = 2 components of the n = 1 MHD mode near the tearing onset. To improve model regularization and prediction stability, singular value decomposition was applied to balance the sawtooth and tearing datasets. The enriched datasets facilitated training neural networks that learn the key distinguishing features of sawtooth and tearing modes in the poloidal profiles of their magnetic amplitude and phase. When the modes occur independently, the networks achieve perfect classification due to the modes’ distinct characteristics and low measurement noise. In the more experimentally relevant case where both modes coexist, the networks maintain exceptional performance across key metrics. Tests on synthetic data with known ground truth demonstrate the superior accuracy of the neural network trained on complex-valued input compared to models using real amplitude, phase, or pseudo-complex data, achieving both a mean time delay and standard deviation below 1 ms. Notably, standard linear regression methods fitting the dominant singular modes to the data closely match the neural network’s performance. Applying these methods across a broad range of H-mode scenarios will enable future studies to systematically identify dominant NTM triggers as scenario-specific variables, paving the way for more effective tearing mode avoidance strategies in future fusion reactor designs.

machine learning↗