Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multi channel data”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Strong-field ionization of water. II. Electronic and nuclear dynamics en route to double ionization

We investigate the role of nuclear motion and strong-field-induced electronic couplings during the double ionization of deuterated water using momentum-resolved coincidence spectroscopy. By examining the three-body dicationic dissociation channel, D$^{+}$/D$^{+}$/O, for both few- and multi-cycle laser pulses, strong evidence for intra-pulse dynamics is observed. The extracted angle- and energy-resolved double ionization yields are compared to classical trajectory simulations of the dissociation dynamics occurring from different electronic states of the dication. In contrast with measurements of single photon double ionization, pronounced departure from the expectations for vertical ionization is observed, even for pulses as short as 10~fs in duration. We outline numerous mechanisms by which the strong laser field can modify the nuclear wavefunction en-route to final states of the dication where molecular fragmentation occurs. Specifically, we consider the possibility of a coordinate-dependence to the strong-field ionization rate, intermediate nuclear motion in monocation states prior to double ionization, and near-resonant laser-induced dipole couplings in the ion. These results highlight the fact that, for small and light molecules such as D$_2$O, a vertical-transition treatment of the ionization dynamics is not sufficient to reproduce the features seen experimentally in the strong field coincidence double-ionization data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ultra-Fast Non-Volatile Resistive Switching Devices with Over 512 Distinct and Stable Levels for Memory and Neuromorphic Computing

Low-current multilevel programmability with inherent non-volatility and high stability of resistance states is required for both multi-bit memory storage and deep learning accelerators but is difficult to achieve. Here, in a resistive switching system, this work realizes >512 (>9 bits) distinct non-volatile conductance levels with stable retention for each state with current levels down to the nanoampere range, highly promising for potential integration with small processing nodes with ultra-low power consumption requirements. This is achieved by demonstrating a new thin film design concept that encompasses three key features: an ultra-thin epitaxial oxygen ionic switching layer that provides a tunable energy barrier at the bottom electrode, an overcoat amorphous layer that acts as an ion migration barrier for stable state retention, and a partial conductive filament as a localized electronic transport channel to the epitaxial switching layer. A large dynamic resistance range of up to seven orders of magnitude is achieved with reset-free transitions among intermediate states, and programmability is demonstrated with ultra-fast (20 ns) pulses. Artificial neural network (ANN) simulations, based on the experimental performance and its non-idealities, demonstrate close-to-ideal inference accuracies for various Modified National Institute of Standards and Technology (MNIST) data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Classifying thermodynamic cloud phase using machine learning models

Vertically resolved thermodynamic cloud-phase classifications are essential for studies of atmospheric cloud and precipitation processes. The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Thermodynamic Cloud Phase (THERMOCLDPHASE) value-added product (VAP) uses a multi-sensor approach to classify the thermodynamic cloud phase by combining lidar backscatter and depolarization, radar reflectivity, Doppler velocity, spectral width, microwave-radiometer-derived liquid water path, and radiosonde temperature measurements. The measured pixels are classified as ice, snow, mixed phase, liquid (cloud water), drizzle, rain, and liq_driz (liquid+drizzle). We use this product as the ground truth to train three machine learning (ML) models to predict the thermodynamic cloud phase from multi-sensor remote sensing measurements taken at the ARM North Slope of Alaska (NSA) observatory: a random forest (RF), a multi-layer perceptron (MLP), and a convolutional neural network (CNN) with a U-Net architecture. Evaluations against the outputs of the THERMOCLDPHASE VAP with 1 year of data show that the CNN outperforms the other two models, achieving the highest test accuracy, F1 score, and mean intersection over union (IOU). Analysis of ML confidence scores shows that ice, rain, and snow have higher confidence scores, followed by liquid, while mixed, drizzle, and liq_driz have lower scores. Feature importance analysis reveals that the mean Doppler velocity and vertically resolved temperature are the most influential data streams for ML thermodynamic cloud-phase predictions. Lidar measurements exhibit lower feature importance due to rapid signal attenuation caused by the frequent presence of persistent low-level clouds at the NSA site. The ML models' generalization capacity is further evaluated by applying them at another Arctic ARM site in Norway using data taken during the ARM Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) field campaign. The models demonstrated similar performance to that observed at the NSA site. Finally, we evaluate the ML models' response to simulated instrument outages and signal degradation and show that a CNN U-Net model trained with input channel dropouts performs better when input fields are missing.

ARM Aerial Facility↗
DOE OSTI2025

DOUBLY-STRANGE BARYON SPECTROSCOPY WITH THE GLUEX EXPERIMENT: GROUND STATE ¿(1320) PHOTOPRODUCTION CROSS SECTION AND PROPERTIES

In the light-quark sector, traditional quark models and, more recently, lattice QCD calculations have predicted many more baryon states than have been experimentally observed. The hyperons, baryons with strange-quark content, in particular abundantly lack these observed states, and past observations are statistically limited. For multi-strange hyperons the matter is much worse and the production mechanisms of these states in photoproduction are very poorly understood. By capitalizing on the narrow peaks, as compared to the broad and overlapping N* and ¿* states, characteristic of the lowest lying ¿ resonances, the systematic aspects of the spectrum can be analyzed. In addition, we can measure experimental observables such as the production cross section to study the production mechanisms of the ¿ resonances. By accumulating various experimental observables, this goal could be accomplished through Partial Wave Analysis (PWA). The main body of work presented in this dissertation is the physics analysis and systematic study of the production cross section for the doubly strange octet ground state baryon (¿(1320)) in the exclusive t-channel photoproduction reaction ¿p ¿ K+Y* ¿ K+(K+¿-)Y* where ¿- ¿ ¿p- and ¿ ¿ pp-. The spin quantum number of ¿(1320) is also extracted and discussed. A preliminary result to extract the mass of ¿(1320) is also explored as a result of the statistical limitations of current measurements. Furthermore, the spectrum of the excited ¿* baryons are explored and discussed in the ¿*- ¿ K-¿ decay channel. The studies on the ground state and excited ¿ baryon are performed using the GlueX Phase-I (GlueX-I) data set. We have extracted differential cross sections from seven beam photon energy bins for E¿ ¿ [6.4,11.4] GeV1 and for -t < 2.4 GeV2 in the production of K+Y*. 1Natural units (c = 1) are used for the entire body of work.

Hernandez, Jesse [Florida State Univ., Tallahassee↗
DOE OSTI2022

Leptonic cascade decays of a heavy Higgs boson through vectorlike leptons at the LHC

We demonstrate the potential of fully leptonic cascade decays of a heavy neutral Higgs boson through vectorlike leptons as a simultaneous probe for extended Higgs sectors and extra matter particles at the LHC. The processes we explore are unique in that their event topologies lead to di-boson-like leptonic final states with a lepton pair which does not reconstruct the mass of a gauge boson. By recasting existing $2\ell$ + $E^{\text{miss}}_{R}$ and $3/4\ell$ searches channels using run2 data from the LHC we obtain model independent bounds on the masses of heavy scalars and vectorlike leptons and use these results to explore future prospects at the HL-LHC. Our results can be directly applied to any kind of new physics scenarios sharing the final states and the event topology. For concreteness, we apply our results to a benchmark scenario: a two Higgs doublet model type-II augmented with vectorlike leptons. Remarkably, even with current data the sensitivity of our analysis shows a reach for masses of a heavy neutral Higgs and vectorlike leptons up to 2 TeV and 1.5 TeV, respectively. Even for low tan β ≳ 1, the analysis retains sensitivity to heavy Higgs masses slightly above 1 TeV. The future sensitivities at the HL-LHC extend the reach for heavy Higgses and new leptons to 2.7 TeV and 2 TeV, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Advanced modeling and simulation of research reactors using dynamic mode decomposition

Full text of publication follows. Due to the ever-increasing safety requirements, the current trend of nuclear reactor analysis is shifting towards high-fidelity multi-physics models, which have a very high computational cost and modelling complexity. As the cost of even a single model run makes it impossible to analyse the behaviour and performance of these models on large-scale commercial plants, it has become even more significant to provide suitable benchmarks to validate and test them extensively. In this sense, research reactors offer a promising solution for the initial validation of high-fidelity models, as they are significantly smaller than commercial reactors and their characteristics are well known. In particular, the reactors of the TRIGA family have been used to assess and validate models and methods for Generation-IV designs, as they have some similar features (such as the dominance of natural convection as cooling mechanism and the difficulties in performing sub-channel analysis using standard codes). Still, the computational requirements of high-fidelity models make them unsuitable for real-time analysis, even following their assessment on research reactors. In this sense, Model Order Reduction (MOR) techniques give an additional strategy to reduce the computational cost of high-fidelity models (whilst preserving sufficient accuracy). In particular, this work focuses on Dynamic Mode Decomposition (DMD), a non-intrusive MOR technique that aims at representing models with explicit temporal dynamics by extracting the time-varying characteristics and the governing structures based only on a set of available data, thus without needing any underlying knowledge of the governing equations. In addition, DMD also computes a low-dimensional surrogate of the dynamic matrix of the system, making it suited for stability analysis and real-time evaluations. This work focuses on the application and validation of the DMD method on the Computational Fluid-Dynamics (CFD) model TRIGA Mark II reactor, also discussing in detail the potentiality of this algorithm as an advanced modelling tool for nuclear reactor analysis. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

System Integration of Rationally Designed Dilute Alloy Catalysts for Energy-Efficient Electrochemical CO 2 -To-Fuel Conversion

Energy effiecient electrochemical conversion of CO 2 to high-demand chemicals and transportation fuels using renewable solar and wind energy is a key technology needed for a high-productivity, low-carbon future. However, the development of scalable, low-cost, active, selective, and stable electrocatalysts remains a key challenge that needs to be overcome to enable high-volume conversion of CO 2 to feedstock chemicals for the chemical industry. In previous work, Lawrence Livermore National Laboratory (LLNL) has developed a rational design platform for dilute alloy transition metal electrocatalysts that promise to make electrochemical CO 2 conversion more energy efficient and selective. In this project, we worked with our industrial partners, Twelve and TotalEnergies, to improve scale up, integration, and stability of LLNL’s dilute alloy catalyst technology into an industry-relevant zero gap electrolyzer platform. Through virtual experiments and data analysis, we designed efficient and cost-effective copper-based catalysts. The catalyst was specifically designed to streamline the slowest and most energy-intensive step of the electrochemical chemical transformation of CO 2 to multi-carbon products – that is making the carbon-carbon bond by dimerization of the reaction intermediate carbon monoxide - resulting in up to 10% improvement in energy efficiency for C 2 products while simultaneously increasing the selectivity towards C 2 products. We tested two different scalable catalyst coating technologies and down-selected magnetron sputtering as the technique that provided the best control over catalyst loading, composition, and morphology. Using this technology, we successfully demonstrated integration of our dilute alloy catalysts into a 100 cm 2 electrolyzer platform with Faradaic efficiencies for ethylene production reaching 40% at a current density of 200 mA/cm 2 . We also developed the technology to integrate a well-defined nanoscale porosity by depositing alloy compositions that were compatible with dealloying, that is, selective leaching of an alloy component to generate nanoscale porosity. We observed that integration of the dealloying-derived nanoporosity improved catalyst performance and stability by leading to a more hydrophobic catalyst/anion exchange membrane interface. Unsolved problems that still need to be addressed are corrosion of the Cu catalyst -specifically if the used catalyst is exposed to air - as well as long term stability do to salt formation/deposition and flooding of the catalyst/electrolyzer flow channels, especially at higher current densities. As we only worked on optimization of catalyst composition, coating thickness, and morphology, further performance optimization will require a system level approach that includes optimization of electrolyzer design and membrane technology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

First operation of a multi-channel Q-Pix prototype: measuring transverse electron diffusion in a gas time projection chamber

We report measurements of the transverse diffusion of electrons in P-10 gas (90% Ar, 10% CH 4 ) in a laboratory-scale time projection chamber (TPC) utilizing a novel pixelated signal capture and digitization technique known as Q-Pix. The Q-Pix method incorporates a precision switched integrating transimpedance amplifier whose output is compared to a threshold voltage. Upon reaching the threshold, a comparator sends a 'reset' signal, initiating a discharge of the integrating capacitor. The time difference between successive resets is inversely proportional to the average current at the pixel in that time interval, and the number of resets is directly proportional to the total collected charge. We developed a 16-channel Q-Pix prototype fabricated from commercial off-the-shelf components and coupled them to 16 concentric annular anode electrodes to measure the spatial extent of the electron swarm that reaches the anode after drifting through the uniform field of the TPC. The swarm is produced at a gold photocathode using pulsed UV light. The measured transverse diffusion agrees with simulations in PyBoltz across a range of operating pressures (200–1500 Torr). These results demonstrate that a Q-Pix readout can successfully reconstruct the ionization topology in a TPC.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Astronomical Spectroscopy with Skipper CCDs: First Results from a Skipper CCD Focal Plane Prototype at SIFS

We present the first on-sky results from an ultra-low-readout-noise Skipper CCD focal plane prototype for the SOAR Integral Field Spectrograph (SIFS). The Skipper CCD focal plane consists of four 6k x 1k, 15 $\mu$m pixel, fully-depleted, p-channel devices that have been thinned to ~250 $\mu$m, backside processed, and treated with an anti-reflective coating. These Skipper CCDs were configured for astronomical spectroscopy, i.e., single-sample readout noise < 4.3 e- rms/pixel, the ability to achieve multi-sample readout noise $\ll$ 1 e- rms/pixel, full-well capacities ~40,000-65,000 e-, low dark current and charge transfer inefficiency (~2 x 10$^{-4}$ e-/pixel/s and 3.44 x 10$^{-7}$, respectively), and an absolute quantum efficiency of $\gtrsim$ 80% between 450 nm and 980 nm ($\gtrsim$ 90% between 600 nm and 900 nm). We optimized the readout sequence timing to achieve sub-electron noise (~0.5 e- rms/pixel) in a region of 2k x 4k pixels and photon-counting noise (~0.22 e- rms/pixel) in a region of 220 x 4k pixels, each with a readout time of $\lesssim$ 17 min. We observed two quasars (HB89 1159+123 and QSO J1621-0042) at redshift z ~ 3.5, two high-redshift galaxy clusters (CL J1001+0220 and SPT-CL J2040-4451), an emission line galaxy at z = 0.3239, a candidate member star of the Bo\"{o}tes II ultra-faint dwarf galaxy, and five CALSPEC spectrophotometric standard stars (HD074000, HD60753, HD106252, HD101452, HD200654). We present charge-quantized, photon-counting observations of the quasar HB89 1159+123 and show the detector sensitivity increase for faint spectral features. We demonstrate signal-to-noise performance improvements for SIFS observations in the low-background, readout-noise-dominated regime. We outline scientific studies that will leverage the SIFS-Skipper CCD data and new detector architectures that utilize the Skipper floating gate amplifier with faster readout times.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Astronomical Spectroscopy with Skipper CCDs: First Results from a Skipper CCD Focal Plane Prototype at SIFS

We present the first on-sky results from an ultra-low-readout-noise Skipper CCD focal plane prototype for the SOAR Integral Field Spectrograph (SIFS). The Skipper CCD focal plane consists of four 6k $\times$ 1k, 15 $\mu$m pixel, fully-depleted, $p$-channel devices that have been thinned to $\sim$250 $\mu$m, backside processed, and treated with an anti-reflective coating. These Skipper CCDs were configured for astronomical spectroscopy, i.e., a single-sample readout noise $< 4.3$\,e$^-$\,rms/pixel, the ability to achieve multi-sample readout noise $\ll 1\,$e$^-$\,rms/pix, full-well capacities $\sim$40,000--65,000 e$^-$, low dark current and charge transfer inefficiency ($\sim 2$ $\times$ $10^{-4}$ e$^-$/pixel/s and 3.44 $\times$ $10^{-7}$, respectively), and an absolute quantum efficiency of $\gtrsim 80\%$ between 450\,nm and 980\,nm ($\gtrsim 90\%$ between 600\,nm and 900\,nm). We optimize the readout sequence timing to achieve sub-electron noise ($\sim 0.5$\,e$^-$\,rms/pix) and photon-counting ($\sim 0.22$\,e$^-$\,rms/pix) readout noise over an area of 2k $\times$ 4k and 110 $\times$ 4k pixels, respectively in a readout time of $\lesssim 17$min. We observed two Lyman-$\alpha$ emitting quasars (HB89\,1159$+$123 and QSO\,J1621$–$0042) at redshift $z \sim 3.5$, two moderate redshift galaxy clusters (CL\,J1001$+$0220 and SPT-CL\,J2040$-$4451), an emission line galaxy at $z = 0.3239$, a candidate member star of the Bo\"{o}tes~II ultra-faint dwarf galaxy, and five CALSPEC spectrophotometric standard stars (HD074000, HD60753, HD106252, HD101452, HD200654). We present charge-quantized, photon-counting observation of the quasar HB89\,1159$+$123 and show the detector sensitivity increase for faint spectral features. We demonstrate signal-to-noise performance improvement for SIFS for observations in the low-background, readout-noise-dominated regime. We outline preliminary scientific studies that will leverage the SIFS-Skipper CCD data to derive new cosmological measurements in the context of dark matter and galaxy evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Nanosecond Gated CMOS Camera (NSGCC) ICD (Rev. 2.1)

The Ultra-Fast X-ray Imager (UXI) program is an ongoing effort at Sandia National Laboratories to create high speed, multi-frame, time-gated Read Out Integrated Circuits (ROICs), and a corresponding suite of photodetectors to image a wide variety of High Energy Density (HED) physics experiments on both Sandia’s Z-Machine and LLNL’s National Ignition Facility (NIF). Several cameras have been designed over the length of the program; one of the most recent is the Icarus, which is an improvement on past imagers (Furi and Hippogriff). A second sensor that can be connected is the Daedalus sensor. The Icarus is a 1024 × 512-pixel array with either 25 μm or 8 µm spatial resolution containing four frames of storage per pixel and has improved timing generation and distribution components while achieved 2 ns time gating. The Daedalus sensor is also a 1024 x 512-pixel array with 25 µm special resolution containing three frames of storage per pixel and has an increased set of features for a wider variety of applications from interlacing of rows in each frame to configurability of all shutters. See Section 14 for details regarding the Icarus implementation of the firmware and Section 15 for details regarding the Daedalus implementation. Due to the unique test environments UXI sensors are targeted for, full custom hardware was required to physically mount an Icarus or Daedalus sensor, manage its various functions, and read out pixel data for transfer to a host computer. Beyond experimental functionality, the hardware also needed to accommodate sensor characterization requirements. Lawrence Livermore National Laboratory’s ‘Version 4.0 Board’ was the result of these efforts. It mounts all the components required to fully utilize the Icarus and Daedalus sensors including analog to digital converters to convert pixel data and various system voltages to digital form for readout and analysis, DAC channels for remote configuration of critical bias voltages, static random-access memories to buffer pixel data, RS422 and Gigabit Ethernet communications for remote access, and an FPGA to tie these components together. This document describes the FPGA electrical interfaces in detail to allow the reader a greater understanding of the device, and to facilitate implementation of custom software to control and manage it. The Version 4.0 Board is a continuation of the Nano-second Gated CMOS hardware design that retains much of the functionality of the Version 1.0 Board while adding features including a DAC instead of digital potentiometers, as well as sensors for pressure and radiation. The Version 4.0 board is intended for applications requiring tight form-factor enclosures. It is composed of two stacking boards; one holds the FPGA and regulators to power the various components of the board while the other contains the mating connector to the sensor, image-readoff ADCs, the DAC, and other components.

42 ENGINEERING↗
DOE OSTI2025

Requirements on the gain calibration for LiteBIRD polarisation data with blind component separation

The detection of primordial B modes of the cosmic microwave background (CMB) could provide information about the early stages of the Universe's evolution. The faintness of this signal requires exquisite calibration accuracy and control of instrumental systematic effects which otherwise could bias the measurements. In this work, we study the impact of an imperfect relative polarisation gain calibration on the recovered value of the tensor-to-scalar ratio r for the LiteBIRD experiment, through the application of the blind Needlet Internal Linear Combination (NILC) foreground-cleaning method. We derive requirements on the relative calibration accuracy of the overall polarisation gain (Δg ν ) for each LiteBIRD frequency channel. Our results show that minimum variance techniques, as NILC, are less sensitive to systematic gain calibration uncertainties compared to a parametric approach, if the latter is not equipped with a proper modelling of these instrumental effects. In this study, the most stringent requirements are found in the channels where the CMB signal is relatively brighter, with the tightest constraints at 166 GHz (Δg ν ≈ 0.16%). This differs from the outcome of an analogous analysis performed with a parametric method, where the tightest requirements are obtained for the foreground-dominated channels. Gain calibration uncertainties, corresponding to the derived requirements, are then simultaneously propagated into all frequency channels. By doing so, we find that the overall impact on estimated r is lower than the total gain systematic budget for LiteBIRD approximately by a factor 5, due to the correlations of the impacts of gain calibration uncertainties in different frequency channels. In order to decouple the systematic effect from the specific choice of the model, we derive the requirements assuming constant spectral parameters for the foreground emission. To assess the robustness of the obtained results against more realistic scenarios, we repeat the analysis assuming sky models of intermediate and high complexity. In these further cases, we adopt an optimised NILC pipeline, called the Multi-Clustering NILC (MC-NILC). We find that the impact of gain calibration uncertainties on r is lower than the LiteBIRD gain systematics budget for the intermediate-complexity sky model. For the high-complexity case, instead, it would be necessary to tighten the requirements by a factor 1.8.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗
DOE OSTI2023

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗
DOE OSTI2020

NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.

97 MATHEMATICS AND COMPUTING↗