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 163 records · Page 9

Waveform Simulation Framework: User Manual with Tutorials

This manuscript describes the Waveform Simulation Framework (WSF), a Python-based framework that provides a unified, programmable interface for generating synthetic seismograms for applications such as seismic array design, method development, and special event analysis. WSF standardizes how users define sources, receivers, and velocity models while abstracting simulator-specific configuration details, enabling workflows that are largely independent of the underlying numerical engine. The document provides installation guidance and tutorial-driven examples for three WSF simulator wrappers—WSF PyFK, WSF SW4, and WSF SPECFEM2D—illustrating end-to-end workflows from forward waveform simulation to common post-processing tasks (e.g., visualization and backprojection) using consistent data products (e.g., ObsPy Stream objects and SAC files).

97 MATHEMATICS AND COMPUTING↗

SymbolNet: neural symbolic regression with adaptive dynamic pruning for compression

Abstract Compact symbolic expressions have been shown to be more efficient than neural network (NN) models in terms of resource consumption and inference speed when implemented on custom hardware such as field-programmable gate arrays (FPGAs), while maintaining comparable accuracy (Tsoi et al 2024 EPJ Web Conf. 295 09036). These capabilities are highly valuable in environments with stringent computational resource constraints, such as high-energy physics experiments at the CERN Large Hadron Collider. However, finding compact expressions for high-dimensional datasets remains challenging due to the inherent limitations of genetic programming (GP), the search algorithm of most symbolic regression (SR) methods. Contrary to GP, the NN approach to SR offers scalability to high-dimensional inputs and leverages gradient methods for faster equation searching. Common ways of constraining expression complexity often involve multistage pruning with fine-tuning, which can result in significant performance loss. In this work, we propose S y m b o l N e t , a NN approach to SR specifically designed as a model compression technique, aimed at enabling low-latency inference for high-dimensional inputs on custom hardware such as FPGAs. This framework allows dynamic pruning of model weights, input features, and mathematical operators in a single training process, where both training loss and expression complexity are optimized simultaneously. We introduce a sparsity regularization term for each pruning type, which can adaptively adjust its strength, leading to convergence at a target sparsity ratio. Unlike most existing SR methods that struggle with datasets containing more than O ( 10 ) inputs, we demonstrate the effectiveness of our model on the LHC jet tagging task (16 inputs), MNIST (784 inputs), and SVHN (3072 inputs).

Tsoi, Ho Fung (ORCID:0000000225502184)↗

Stairway to discovery: A report on the CMS programme of cross section measurements from millibarns to femtobarns

The Large Hadron Collider at CERN, delivering proton-proton collisions at much higher energies and far higher luminosities than previous machines, has enabled a comprehensive programme of measurements of the standard model (SM) processes by the CMS experiment. These unprecedented capabilities facilitate precise measurements of the properties of a wide array of processes, the most fundamental being cross sections. The discovery of the Higgs boson and the measurement of its mass became the keystone of the SM. Knowledge of the mass of the Higgs boson allows precision comparisons of the predictions of the SM with the corresponding measurements. These measurements span the range from one of the most copious SM processes, the total inelastic cross section for proton-proton interactions, to the rarest ones, such as Higgs boson pair production. They cover the production of Higgs bosons, top quarks, single and multibosons, and hadronic jets. Associated parameters, such as coupling constants, are also measured. These cross section measurements can be pictured as a descending stairway, on which the lowest steps represent the rarest processes allowed by the SM, some never seen before.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

HDBind: encoding of molecular structure with hyperdimensional binary representations

Traditional methods for identifying “hit” molecules from a large collection of potential drug-like candidates rely on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have a significant limitation in that they require exceptional computing capabilities for even relatively small collections of molecules. Increasingly large and complex state-of-the-art deep learning approaches have gained popularity with the promise to improve the productivity of drug design, notorious for its numerous failures. However, as deep learning models increase in their size and complexity, their acceleration at the hardware level becomes more challenging. Hyperdimensional Computing (HDC) has recently gained attention in the computer hardware community due to its algorithmic simplicity relative to deep learning approaches. The HDC learning paradigm, which represents data with high-dimension binary vectors, allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas (computer vision, bioinformatics, mass spectrometery, remote sensing, edge devices, etc.). To the best of our knowledge, our work is the first to consider HDC for the task of fast and efficient screening of modern drug-like compound libraries. We also propose the first HDC graph-based encoding methods for molecular data, demonstrating consistent and substantial improvement over previous work. We compare our approaches to alternative approaches on the well-studied MoleculeNet dataset and the recently proposed LIT-PCBA dataset derived from high quality PubChem assays. We demonstrate our methods on multiple target hardware platforms, including Graphics Processing Units (GPUs) and Field Programmable Gate Arrays (FPGAs), showing at least an order of magnitude improvement in energy efficiency versus even our smallest neural network baseline model with a single hidden layer. Our work thus motivates further investigation into molecular representation learning to develop ultra-efficient pre-screening tools. We make our code publicly available at https://github.com/LLNL/hdbind.

59 BASIC BIOLOGICAL SCIENCES↗

A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots

While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.

97 MATHEMATICS AND COMPUTING↗

Integrated modeling of alpha particle losses in JET DT plasmas

Alpha particle confinement is crucial for sustaining burning plasmas and designing future reactor concepts. Along with classical/prompt losses, various magnetohydrodynamic instabilities can lead to wave-particle interactions which can transport alpha particles outward from the plasma. This can result in a reduction in plasma heating/performance, and, at worst, damage in-vessel components. Joint European Torus’s recent deuterium–tritium campaigns in 2021–2023 have produced numerous alpha particle loss measurements with its scintillator probe and Faraday cup array fast ion loss detectors as discussed in Bonofiglo et al (2024 Nucl. Fusion 64 096038). This paper will report on integrated energetic particle transport modeling in support of those measurements. The modeling is accomplished with the TRANSP and ORBIT-kick codes with the use of recently developed reduced models which calculate mode structure, amplitude, and the evolving dynamics. When possible, constraints and comparisons to experiment are conducted. Case studies are performed on a variety of magnetohydrodynamic activity, including: fishbones, tearing modes (TMs), and sawtooth crashes. Additionally, a special case of an alpha-driven toroidal Alfvén eigenmode is briefly discussed, where modeling showed marginally weak alpha losses and was unable to support experimental observations. Coupled effects between a TM and toroidal field ripple are presented and were unable to replicate the observations in lost particle pitch but did duplicate the localized flattening of the measured neutron profile. Additional modeling results compare the magnitude of losses and energy/velocity-space sensitivities against experimental observations/measurements for each scenario. This work corroborates numeric alpha transport modeling while also identifying model deficiencies. While this report details alpha transport, it also presents open issues for discussion in assessing the validity of our numerical models towards burning plasmas.

Bonofiglo, Phillip J. [Princeton Plasma Physics La↗

Towards a new generation of solid total-energy detectors for neutron-capture time-of-flight experiments with intense neutron beams

Challenging neutron-capture cross-section measurements of small cross sections and samples with a very limited number of atoms require high-flux time-of-flight facilities. In turn, such facilities need innovative detection setups that are fast, have low sensitivity to neutrons, can quickly recover from the so-called γ-flash, and offer the highest possible detection sensitivity. In this paper, we present several steps towards such advanced systems. Specifically, we describe the performance of a high-sensitivity experimental setup at CERN n_TOF EAR2. It consists of nine sTED detector modules in a compact cylindrical configuration, two conventional used large-volume C 6 D 6 detectors, and one LaCl 3 (Ce) detector. The performance of these detection systems is compared using 93 Nb(n, γ) data. We also developed a detailed G EANT small4 Monte Carlo model of the experimental EAR2 setup, which allows for a better understanding of the detector features, including their efficiency determination. This Monte Carlo model has been used for further optimization, thus leading to a new conceptual design of a γ detector array, STAR, based on a deuterated-stilbene crystal array. Finally, the suitability of deuterated-stilbene crystals for the future STAR array is investigated experimentally utilizing a small stilbene-d12 prototype. The results suggest a similar or superior performance of STAR with respect to other setups based on liquid-scintillators, and allow for additional features such as neutron-gamma discrimination and a higher level of customization capability.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Supercooled phase transitions in conformal dark sectors explain NANOGrav data

According to recent lore, it is difficult to explain the evidence for a stochastic gravitational wave background obtained by pulsar timing arrays with supercooled first-order phase transitions (FOPTs). We demonstrate that supercooled FOPTs in dark U(1) models with a conformal dark sector easily explain the nHz signal at NANOGrav.

Conformal models↗

Reaction Diffusion Modelling of 3D Pillar Electrodes in Single-Catalyst CO 2 Reduction Cascades

Effective electrochemical CO 2 reduction to liquid fuels requires that the local catalytic environment facilitates the desired reactivity, yet a microscopic understanding of this environment is difficult to achieve from experiment alone. In this work, a 3D reaction-diffusion model was developed to explore the effects of electrode surface area and local geometry on the performance of a heterogeneous catalyst that performs a two-step CO 2 reduction cascade reaction to CO and then CH 3 OH under aqueous conditions. Kinetic parameters for the model were inspired by experimental results using a cobalt phthalocyanine (CoPc) catalyst. Three-dimensional architectures composed of arrays of square pillars with varying dimensions and either smooth or periodically modulated surfaces were tested, revealing the extent to which geometry modulates the performance of the cascade reactions. Although structural variations modulate local concentration gradients, we find that electrochemically active surface area predominantly governs the overall cascade reaction. Moreover, the results suggest that supersaturation of CO, with concentrations up to ten-fold higher than the equilibrium solubility limit, might be critical for more efficient conversion to CH 3 OH. For any given geometry, the spatially averaged ratio of [CO] to [CO 2 ] is dictated by the electrochemically active surface area and determines the yield of CH 3 OH. For a fixed surface area, geometries that spatially confine the electrolyte yield moderate local [CO] to [CO 2 ] ratios within small volumes. In contrast, less confining geometries result in a broader distribution of local ratios spread over larger volumes, with both configurations yielding the same spatially averaged [CO] to [CO 2 ] ratio. These insights provide valuable design principles—highlighting the critical importance of surface area and possibly CO supersaturation—for engineering advanced electrode architectures that leverage intermediate trapping and CO supersaturation to enhance overall performance in tandem CO 2 reduction systems.

COMSOL↗

The Functor system: a new on-the-fly take on Material Properties based on C++ functions

In the context of solving multiphysics problems, the discretization of the partial differential equations (PDE) at hand often takes the spotlight. However, for most engineering users and even application developers, the discretization of the equations has already been performed. Instead, they are tasked with implementing specific closure relations and material properties. MOOSE has long enabled this using the Materials system. This system relied on the pre-computation of all properties before they are used in the PDE or in postprocessing. In this talk we will introduce the Functor system, which was deployed in MOOSE in 2021, then present a few applications of functors in flow modeling simulations by the NEAMS program. Functors first offer great flexibility in their evaluation. Rather than storing various arrays for material properties, they are evaluated on the fly at the location and state, e.g. current or old value, requested. Unlike regular material properties, several operations such as the time derivative, the divergence and the curl can be requested from a functor. Similar to material properties, functors can be made to depend on arbitrary combinations of variables, functions, postprocessors and other properties. However, unlike material properties, any of these can be substituted for a functor material property. Thanks to this, objects no longer need to be duplicated based on the types of their parameters.

97 - MATHEMATICS AND COMPUTING↗

The VIBES Are Shifting: Assessing Emergent Capabilities in Multi-Modal Models

Researchers assessing open-source domains such as the internet, and particularly those studying information conflict, often have no single prescribed workflow. In the course of their research, they may need to perform a diverse array of tasks far beyond simply identifying an ever-changing set of objects. These analytical tasks can include ascertaining the provenance of images, understand an image in the context of accompanying text-based data, or identifying indicators of digital image manipulation. They must further be able to do this at the scale of tens of thousands of images or more. Traditional machine vision models have typically lacked the flexibility and breadth of performance sufficient for these needs. The research team from Pacific Northwest National Laboratory assessed the performance of a single baseline CLIP ViT-L model against a series of analytical tasks relevant for the study of online information conflict.

97 MATHEMATICS AND COMPUTING↗

The U.S. Agrivoltaic Shading Tool: A National-Scale Interface for Modeling Light and Shade Patterns in Ten Common Agrivoltaic Configurations

Agrivoltaic systems are dual-use configurations that co-locate agriculture and photovoltaic (PV) infrastructure and require careful design to balance crop performance and energy generation. A critical element of agrivoltaic design is the spatial and temporal distribution of irradiance and shade within and around PV arrays. To support research, planning, and stakeholder decision-making, we introduce the U.S. Agrivoltaic Shading Tool, a novel web-based application that delivers high-resolution irradiance and photosynthetically active radiation (PAR) modeling for ten standardized PV configurations across the conterminous United States. The tool leverages the National Laboratory of the Rockies (NLR) System Advisor Model (SAM) to perform detailed irradiance simulations, using meteorological data from the National Solar Radiation Database (NSRDB). Outputs include seasonal, monthly, weekly, and diurnal patterns of available sunlight, amount of shade, irradiance, and PAR at ground level within agrivoltaic system footprints. For a user's selected location, these results are visualized through interactive visualizations, heatmaps, and time-series plots, designed to be accessible to both technical and non-technical users. In addition to facilitating rapid spatial exploration of agrivoltaic light environments, the tool will offer seamless integration with the InSPIRE Agrivoltaics Design and Analysis Model (ADAM). This optional workflow will allow users to port selected site and configuration parameters into a more advanced modeling environment for further customization of structural layouts, crop-system compatibility, power generation, and technoeconomic performance. Finally, to promote open science, the entire dataset will be hosted and available for open access through the OpenEI platform. By standardizing and disseminating high-quality irradiance data and design tools, the U.S. Agrivoltaic Shading Tool supports a wide range of users, including researchers, landowners, energy developers, and policymakers, in evaluating the agronomic and energetic feasibility of agrivoltaic systems across the United States.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Fast and flexible neutrino decoupling. Part I. The Standard Model

Cosmological determinations of the number of relativistic neutrino species, N eff , are becoming increasingly accurate, and further improvements are expected both from CMB and BBN data. Given this context, we update the evaluation of N eff and the current entropy density via the momentum-averaged approach. This allows for a numerically fast description of neutrino decoupling, easily portable to an array of new physics scenarios. We revisit all aspects of this approach, including collision terms with full electron mass dependence, finite temperature QED corrections to the equation of state, neutrino oscillations, and the modelling of neutrino ensembles with effective chemical potentials. For integrated observables, our results differ by less than 0.04% from the solution of the momentum-dependent evolution equation. We outline how to extend the approach to BSM settings, and will highlight its power in Part II. To facilitate the practical implementation, we release a Mathematica and python code within nudec_BSM_v2, easily linkable to BBN codes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Runtime performance of a GAMESS quantum chemistry application offloaded to GPUs

Summary Computational chemistry is at the forefront of solving urgent societal problems, such as polymer upcycling and carbon capture. The complexity of modeling these processes at appropriate length and time scales is mainly manifested in the number and types of chemical species involved in the reactions and may require models of several thousand atoms and large basis sets to accurately capture the chemical complexity and heterogeneity in the physical and chemical processes. The quantum chemistry package General Atomic and Molecular Electronic Structure System (GAMESS) has a wide array of methods that can efficiently and accurately treat complex chemical systems. In this work, we have used the GAMESS Effective Fragment Molecule Orbital (EFMO) method for electronic structure calculation of a challenging mesoporous silica nanoparticle (MSN) model surrounded by about 4700 water molecules to investigate the strong scaling and GPU offloading on hybrid CPU‐GPU nodes. Experiments were performed on the Perlmutter platform at the National Energy Research Scientific Computing Center. Good strong scaling and load balancing have been observed on up to 88 hybrid nodes for different settings of the execution parameters for the calculation considered here. When GPUs are oversubscribed by offloading work from multiple CPU processes, using the NVIDIA multi‐process service (MPS) has consistently reduced time to solution and energy consumed. Additionally, for some configuration parameter settings, oversubscription with MPS improved performance by up to 5.8% over the case without oversubscription.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Multichannel Silicon Package for Large-Scale Skipper-CCD Experiments

The next generation of experiments for light-dark-matter and neutrino searches based on skipper charge-coupled devices (skipper-CCDs) introduces new challenges for the sensor packaging and readout architecture. Scaling the active mass while simultaneously reducing the experimental backgrounds in orders of magnitude requires a novel high-density silicon-based package that must be massively produced and tested. In this work, we present a silicon multichip module design capable of hosting up to 16 skipper-CCDs, along with the fabrication process for the first prototypes. We thoroughly tested and characterized the assembled prototypes to build an empirical model for the video output signal of skipper-CCDs integrated into the silicon package. We then used this model in simulations to optimize the fabrication process and achieve the robust performance required for the full-scale array, which we validated through a new round of prototype fabrication. We outline the specifications selected for the ongoing production of 1500 silicon wafers that will ultimately add up to a 10-kg skipper-CCD array with 24000 readout channels.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Mechanical sensors for ultraheavy dark matter searches via long-range forces

Dark matter candidates with masses around the Planck scale are theoretically well-motivated, and it has been suggested that it might be possible to search for dark matter solely via gravitational interactions in this mass range. In this work, we explore the pathway towards searching for dark matter candidates with masses around the Planck scale using mechanical sensors while considering realistic experimental constraints, and develop analysis techniques needed to conduct such searches. These dark matter particles are expected to leave tracks as their signature in mechanical sensor arrays, and we show that we can effectively search for such tracks using statistical approaches to track-finding. We analyze a range of possible experimental setups and compute sensitivity projections for searches for ultraheavy dark matter coupling to the Standard Model via long-range forces. We find that while a search for Planck-scale dark matter purely via gravitational couplings would be exceedingly difficult, requiring ∼80 dB of quantum noise reduction with a 100 3 array of devices, there is a wide range of currently unexplored dark matter candidates which can be searched for with already existing or near-term experimental platforms.

Qin, Juehang [Rice Univ., Houston, TX (United Stat↗

Scalable Thin Light-Emitting Diode (LED) Light Sheet Platform

The goal of the Scalable Thin Light-Emitting Diode (LED) Light Sheet Platform project is to derisk the manufacturing scalability for a high chip count heterogeneous smart lighting platform, with the potential to enable significant energy savings through dynamic directional light control and custom. Specifically, the project derisks the scalability of a disruptive new computer controlled microassembly fabrication process SRI International is developing, to address fundamental cost barriers to mass production of high chip count systems. This is key for enabling mass adoption and thus maximum societal energy savings impact. The LED light sheet technology is a smart illumination platform that can be applied to lighting, signage, and display. It has features and a form factor similar to bendable or conformal OLED light sheets but uses more efficient LEDs with a remote phosphor layer that is very close to the LED to deliver a luminous efficacy that exceeds 125 lumen/Watt and enable > 50% Lighting Application Efficiency (LAE) energy savings. During Budget Period 1 (BP1), all of the BP1 milestones were successfully completed: (M2.0.1) Design demonstrator details that will meet the final project goals, (M6.1.1) Show lighting output model can predict illuminance, spatial, spectral distribution for specific light-sheet designs and use case, (M3.2.1) Automated loading supports 200 chip arrays, and (M3.3.1) Finish first interconnect process run for 200 chip array. The BP1 Go/No-Go Decision Point (G/NG 1) 200 chip demonstrator sample was in the process of final assembly and test at the end of BP1 on December 31, 2023. The objective was successfully achieved on January 16, 2024, when a wired 200-chip sample underwent confirmation testing that demonstrated a 93.75% first pass electrical yield that successfully exceeded the 75% first-pass electrical yield BP1 G/NG requirement. In addition, all 93.75% of the sample devices lit up, further demonstrating the ability to assemble and transfer small LED chips without damaging their functionality.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗