Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “edge inference”

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 181 records · Page 10

Accelerating the Inference of the Exa.TrkX Pipeline

Recently, graph neural networks (GNNs) have been successfully used for a variety of particle reconstruction problems in high energy physics, including particle tracking. The Exa.TrkX pipeline based on GNNs demonstrated promising performance in reconstructing particle tracks in dense environments. It includes five discrete steps: data encoding, graph building, edge filtering, GNN, and track labeling. All steps were written in Python and run on both GPUs and CPUs. In this work, we accelerate the Python implementation of the pipeline through customized and commercial GPU-enabled software libraries, and develop a C++ implementation for inferencing the pipeline. The implementation features an improved, CUDA-enabled fixed-radius nearest neighbor search for graph building and a weakly connected component graph algorithm for track labeling. GNNs and other trained deep learning models are converted to ONNX and inferenced via the ONNX Runtime C++ API. The complete C++ implementation of the pipeline allows integration with existing tracking software. We report the memory usage and average event latency tracking performance of our implementation applied to the TrackML benchmark dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Kinetic Modeling and Experimental Studies of Capacitively Coupled Macro- and Microplasmas at Microwave Frequencies

The effect of driving frequency in the range of 13.56 MHz to 73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct temperatures for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2 eV to 10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25 % of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically-derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mass transfer from a circular cylinder: Effects of flow unsteadiness and slight nonuniformities

Experiments were performed to determine the effect of periodic variations in the angle of the flow incident to a turbine blade on its leading edge heat load. To model this situation, measurements were made on a circular cylinder oscillating rotationally in a uniform steady flow. A naphthalene mass transfer technique was developed and used in the experiments and heat transfer rates are inferred from the results. The investigation consisted of two parts. In the first, a stationary cylinder was used and the transfer rate was measured for Re = 75,000 to 110,000 and turbulence levels from .34 percent to 4.9 percent. Comparisons with both theory and the results of others demonstrate that the accuracy and repeatability of the developed mass transfer technique is about + or - 2 percent, a large improvement over similar methods. In the second part identical flow conditions were used but the cylinder was oscillated. A Strouhal number range from .0071 to .1406 was covered. Comparisons of the unsteady and steady results indicate that the magnitude of the effect of oscillation is small and dependent on the incident turbulence conditions.

Marziale, M. L.↗

Antarctic sea ice concentrations inferred from Nimbus 5 ESMR and Landsat imagery

Landsat Multispectral Scanner near-infrared band images (MSS-7) are used to calculate Antarctic ice concentration values and these results are quantitatively compared to those derived from Nimbus 5 electrically scanning microwave radiometer (ESMR) brightness temperature data. The set of images includes regions near the ice edge and near the continental boundary during late winter or spring. The percentage of ice cover is derived from the microwave data, using an algorithm that incorporates the brightness temperature emissivity and interpolated climatological physical temperatures and that is not dependent on spatial resolution. Cloud-free Landsat images are used to provide an alternative determination of large-scale ice concentrations, but this method is limited in its ability to resolve individual ice flows, especially near the ice edge. A proportional classification procedure is applied to Landsat data based on observed reflection and the results are correlated with similarly-derived ESMR ice concentration values. The corresponding ice concentration values are found to agree to + or - 15%, with the largest uncertainties occurring in regions of low ice concentration.

Comiso, J. C.↗

Persistence of Antarctic polar stratospheric clouds

The persistence of Polar Stratospheric Clouds (PSCs) observed by the Stratospheric Aerosol Measurement (SAM) 2 satellite sensor over a 9-year period is compared and contrasted. Histograms of the SAM 2 1.0 micron extinction ratio data (aerosol extinction normalized by the molecular extinction) at an altitude of 18 km in the Antarctic have been generated for three 10-day periods in the month of September. Statistics for eight different years (1979 to 1982 and 1984 to 1987) are shown in separate panels for each figure. Since the SAM 2 system is a solar occultation experiment, observations are limited to the edge of the polar night and no measurements are made deep within the vortex where temperatures could be colder. For this reason, use is made of the NMC global gridded fields and the known temperature-extinction relationship to infer additional information on the occurrence and areal coverage of PSCs. Calculations of the daily areal coverage of the 195 K isotherm will be presented for this same period of data. This contour level lies in the range of the predicted temperature for onset of the Type 1 particle enhancement mode at 50 mb (Poole and McCormick, 1988b) and should indicate approximately when formation of the binary HNO3-H2O particles begins.

Mccormick, M. Patrick↗

Do graph neural networks learn traditional jet substructure?

At the CERN LHC, the task of jet tagging, whose goal is to infer the origin of a jet given a set of final-state particles, is dominated by machine learning methods. Graph neural networks have been used to address this task by treating jets as point clouds with underlying, learnable, edge connections between the particles inside. We explore the decision-making process for one such state-of-the-art network, ParticleNet, by looking for relevant edge connections identified using the layerwise-relevance propagation technique. As the model is trained, we observe changes in the distribution of relevant edges connecting different intermediate clusters of particles, known as subjets. The resulting distribution of subjet connections is different for signal jets originating from top quarks, whose subjets typically correspond to its three decay products, and background jets originating from lighter quarks and gluons. This behavior indicates that the model is using traditional jet substructure observables, such as the number of prongs -- energetic particle clusters -- within a jet, when identifying jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry

Plasma facing components (PFCs) in the next generation of tokamak devices will operate in challenging environments, with heat loads predicted to exceed 10 MWm -2 . The magnitude of these heat loads is set by the width of the channel, the ‘scrape-off layer’ (SOL), into which heat is exhausted, and can be characterised by an e-folding length scale for the decay of heat flux across the channel. It is expected this channel will narrow as tokamaks move towards reactor relevant conditions. Understanding the processes involved in setting the SOL heat flux width is imperative to be able to predict the heat loads PFCs must handle in future devices. Measurements of the SOL width are performed on the high-field spherical tokamak, ST40, using a newly commissioned infrared thermography system. With its high on-axis toroidal magnetic field (≥1.5 T) ST40 is uniquely positioned to investigate the influence of toroidal field on the heat flux width in spherical tokamaks, whilst also extending measurements of the SOL width in spherical tokamaks to increased poloidal field (≥0.3 T). Due to the divertor on ST40 having a low degree of axisymmetry, it is necessary for a set of radial measurements of the heat flux to be taken across the divertor, made possible using an automated toolchain that fully incorporates its 3D geometry. These radial profiles are combined with the magnetic geometry of the plasma to infer the width of the SOL, with both Eich and double exponential profiles of heat flux observed. A reduction in the heat flux is observed toroidally across part of the divertor, along with increased heat loads observed locally around the edges of the tiles. Future work in characterising the impact of tile misalignment and uncertainties in the reconstructed divertor magnetic geometry is required in order to further understand the observed heat flux patterns, as are additional investigations into the role potentially being played by an inhomogeneous sheath electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Height of Shock Formation in the Solar Corona Inferred from Observations of Type II Radio Bursts and Coronal Mass Ejections

Employing coronagraphic and EUV observations close to the solar surface made by the Solar Terrestrial Relations Observatory (STEREO) mission, we determined the heliocentric distance of coronal mass ejections (CMEs) at the starting time of associated metric type II bursts. We used the wave diameter and leading edge methods and measured the CME heights for a set of 32 metric type II bursts from solar cycle 24. We minimized the projection effects by making the measurements from a view that is roughly orthogonal to the direction of the ejection. We also chose image frames close to the onset times of the type II bursts, so no extrapolation was necessary. We found that the CMEs were located in the heliocentric distance range from 1.20 to 1.93 solar radii (Rs), with mean and median values of 1.43 and 1.38 Rs, respectively. We conclusively find that the shock formation can occur at heights substantially below 1.5 Rs. In a few cases, the CME height at type II onset was close to 2 Rs. In these cases, the starting frequency of the type II bursts was very low, in the range 25-40 MHz, which confirms that the shock can also form at larger heights. The starting frequencies of metric type II bursts have a weak correlation with the measured CME/shock heights and are consistent with the rapid decline of density with height in the inner corona.

EUV↗

FPGA Acceleration of GCN in Light of the Symmetry of Graph Adjacency Matrix

Graph Convolutional Neural Networks (GCNs) are widely used to process large-scale graph data. Different from deep neural networks (DNNs), GCNs are sparse, irregular, and unstructured, posing unique challenges to hardware acceleration with regular processing elements (PEs). In particular, the adjacency matrix of a GCN is extremely sparse, leading to frequent but irregular memory access, low spatial/temporal data locality and poor data reuse. Furthermore, a realistic graph usually consists of unstructured data (e.g., unbalanced distributions), creating significantly different processing times and imbalanced workload for each node in GCN acceleration. To overcome these challenges, we propose an end-to-end hardware-software co-design to accelerate GCNs on resource-constrained FPGAs with the features including: (1) A custom dataflow that leverages symmetry along the diagonal of the adjacency matrix to accelerate feature aggregation for undirected graphs. We utilize either the upper or the lower triangular matrix of the adjacency matrix to perform aggregation in GCN to improve data reuse. (2) Unified compute cores for both aggregation and transform phases, with full support to the symmetry-based dataflow. These cores can be dynamically reconfigured to the systolic mode for transformation or as individual accumulators for aggregation in GCN processing. (3) Preprocessing of the graph in software to rearrange the edges and features to match the custom dataflow. This step improves the regularity in memory access and data reuse in the aggregation phase. Moreover, we quantize the GCN precision from FP32 to INT8 to reduce the memory footprint without losing the inference accuracy. We implement our accelerator design in Intel Stratix10 MX FPGA board with HBM2, and demonstrate 1.3x-110.5x improvement in end-to-end GCN latency as compared to the state-of the-art FPGA implementations, on the graph datasets of Cora, Pubmed, Citeseer and Reddit.

Nair, Gopikrishnan R.↗

Geomagnetic storm fields near a synchronous satellite.

An apparent early recovery of the main phase of geomagnetic storms at the distance of the synchronous satellite is examined in terms of changing electric current distributions in the magnetosphere during magnetic storms. It is suggested that a rapid recession of the edge of the plasma sheet (after the advance toward the earth during an early epoch of the main phase) is partly responsible for the early recovery. Relevant plasma sheet variations during geomagnetic storms are found to be in agreement with the inferred variations.

Kawasaki, K.↗

The radio 'lobes' of Scorpius X-1 are unrelated background sources

The VLA between 1981 and 1990 are used to produce high-resolution images of the radio emission from the region near Sco X-1. The radio proper motion of Sco X-1 was measured at 0.0148 +/-0.0011/yr, which agrees with the optical determination for the X-ray-emitting binary system. The proper motions of two nearby radio sources, juxtaposed 1 arcmin to the NE and to the SW, were measured and found to be stationary in the sky with upper limits of 0.004 arcsec/yr. A deep radio image of the 10-arcsec extended SW source shows a morphology strikingly similar to that of a typical luminous extragalactic radio source, which contains two edge-brightened lobes, a jet, and a core. The possibilities that the NE source, although nearly stationary in the sky, is associated with Sco X-1, are discussed, and it is concluded that it is an unrelated background source. It is inferred that Sco X-1 is not a miniature triple source or a 'microquasar', and its radio emission is confined to the binary system.

Fomalont, E. B.↗

Resonant X-ray emission spectroscopy from broadband stochastic pulses at an X-ray free electron laser

Abstract Hard X-ray spectroscopy is an element specific probe of electronic state, but signals are weak and require intense light to study low concentration samples. Free electron laser facilities offer the highest intensity X-rays of any available light source. The light produced at such facilities is stochastic, with spikey, broadband spectra that change drastically from shot to shot. Here, using aqueous ferrocyanide, we show that the resonant X-ray emission (RXES) spectrum can be inferred by correlating for each shot the fluorescence intensity from the sample with spectra of the fluctuating, self-amplified spontaneous emission (SASE) source. We obtain resolved narrow and chemically rich information in core-to-valence transitions of the pre-edge region at the Fe K-edge. Our approach avoids monochromatization, provides higher photon flux to the sample, and allows non-resonant signals like elastic scattering to be simultaneously recorded. The spectra obtained match well with spectra measured using a monochromator. We also show that inaccurate measurements of the stochastic light spectra reduce the measurement efficiency of our approach.

Fuller, Franklin D. (ORCID:0000000237737087)↗

Characterization of core neutrals using passive measurements of the D α spectrum near the X -point in the DIII-D tokamak

New spectroscopic measurements of deuterium Balmer-α emission are used to infer the spatial distribution of neutral particles near the X-point in a diverted high confinement mode plasma. The charge exchange neutral spectroscopy (CENS) diagnostic—recently installed on the DIII-D tokamak—uses 15 lines of sight extending from the edge of the confined region (ρ ≈ 0.7) to the X-point in lower single-null plasmas. Each CENS line of sight is spectrally resolved around the deuterium Balmer-α wavelength (6561 Å) to measure the Doppler shift, line broadening, and intensity of Dα emission. Thermal line broadening is used to identify emission from neutrals undergoing charge-exchange (CX) with high temperature ions in the confined plasma. This spectral information allows a more accurate determination of the neutral density deeper inside the confined plasma compared to traditional filter-based diagnostics. We present two methods of analyzing CENS measurements: (1) fitting the spectrum from each view-chord separately, and (2) a tomographic method for inverting the neutral density over a 2D region of space using all CENS views collectively. The neutral density profile is found to decay exponentially in the radial direction with two scale-lengths. In the pedestal region the neutral density decays at a rate approximately equal to the local mean-free-path for CX collisions, $L_{n_\mathrm{D0}}\approx \lambda_\text{CX}$. Further inside the plasma the neutral density decays at a rate equal to the mean-free-path for ionization, $L_{n_\mathrm{D0}}\approx \lambda_\text{inz.}$. The separatrix $n_\mathrm{D0}$ value is found to be approximately $2\times10^{15}$ m−3. Based on results from the 2D inversion, the density of neutrals is found to increase along the separatrix approaching the X-point by a factor of $\approx 5\times$ over the region covered by the CENS diagnostic.

X-point neutrals↗

Edge termination in vertical GaN diodes: Electric field distribution probed by second harmonic generation

We characterized the electric field distribution of GaN-on-GaN p–n diodes with partially compensated ion-implanted edge termination (ET) using an electric field induced second harmonic generation technique (EFISHG). The distributed electric field from the anode to the outer edge of the ET demonstrates the effectiveness of the ET structure. However, EFISHG also shows that its effectiveness is strongly dependent on the acceptor charge distribution in the ET's partially compensated layer (PC). A generally lower amount of acceptor charge can be inferred from the measured electric field distribution resulting from excessive ion implantation energy or dose during ET fabrication and causing lower than optimal breakdown voltage. Localized field crowding can be observed when the remaining acceptors uncompensated by the implant in the PC layer are nonuniformly distributed around the periphery of the devices. Important information can be obtained from these direct electric field measurements and used for optimizing the device design and fabrication process.

Cao, Yuke↗

3D radiation, density, and MHD structures following neon shattered pellet injection into stable DIII-D Super H-mode discharges

Six nominally repeat neon shattered pellet injection (SPI) shutdowns of stable DIII-D Super H-modes are studied to understand the 3D properties of the radiation and impurity transport. The radiation efficiency and radiation peaking determine whether first wall melting is expected following disruption mitigation in ITER. Previous studies make use of axisymmetric approximations to infer radiation efficiencies, but validating the high efficiency required by ITER necessitates improved accuracy, and this work contributes by exploring the 3D radiation and density structures that will inform forward modeling. When the neon shatter plume produced by the SPI reaches the plasma edge, m/n = 3/1 and 2/1 island O-points are observed to align with the injection trajectory in five out of six cases, suggesting that the injected material seeds the island O-points. Field aligned neon structures emitting Ne-I line radiation drift at 1 km/s in the ion diamagnetic drift direction during the pre-TQ, tracking the motion of the m/n = 2/1 island O-point. Neon fragments penetrate to the q = 2 surface by the time of the TQ. Techniques to constrain the 3D emissivity are explored, and one method constrains a 3D flux tube that is consistent with the radiation data, and when mapped to the interferometers, intersects the lasers that measure the highest density. The resulting structure derived from the radiation measurements exists near the 2/1 island X-point. In five repeatable discharges, the peak of the radiation in the toroidal direction exists in a 120° toroidal sector where the injection occurs, in contrast with the outlier discharge where the toroidal peak exists in the complementary 240° toroidal sector far from the injector, and where a 50% lower density rise is observed. In conclusion, the n = 1 phase behavior is markedly different in the outlier discharge, suggesting a possible dependence of the radiation structure and the assimilation efficiency on MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory and Development of Position-Sensitive Quantum Calorimeters

Quantum calorimeters are being developed as imaging spectrometers for future X-ray astrophysics observatories. Much of the science to be done by these instruments could benefit greatly from larger focal-plane coverage of the detector (without increasing pixel size). An order of magnitude more area will greatly increase the science throughput of these future instruments. One of the main deterrents to achieving this goal is the complexity of the readout schemes involved. We have devised a way to increase the number of pixels from the current baseline designs by an order of magnitude without increasing the number of channels required for readout. The instrument is a high energy resolution, distributed-readout imaging spectrometer called a Position-Sensitive Transition-Edge Sensor (POST). A POST is a quantum calorimeter consisting of two Transition-Edge Sensors (TESS) on the ends of a long absorber capable of one-dimensional imaging spectroscopy. Comparing rise time and energy information from the two TESS, the position of the event in the POST is determined. The energy of the event is inferred from the sum of the two pulses. We have developed a generalized theoretical formalism for distributed-readout calorimeters and apply it to our devices. We derive the noise theory and calculate the theoretical energy resolution of a POST. Our calculations show that a 7-pixel POST with 6 keV saturation energy can achieve 2.3 eV resolution, making this a competitive design for future quantum calorimeter instruments. For this thesis we fabricated 7- and 15-pixel POSTS using Mo/Au TESs and gold absorbers, and moved from concept drawings on scraps of napkins to a 32 eV energy resolution at 1.5 keV, 7-pixel POST calorimeter.

Figueroa-Feliciano, Enectali↗

Verification, Validation, and Calibration Through a Causal Lens

While typical validation and verification approaches focus on identifying the associations between data elements using statistical and machine learning methods, the novel methods in this paper focus instead on identifying causal relationships between data elements. Statistical and machine-learning-based approaches are strictly data-driven, meaning that they provide quantitative comparison measures between data sets without explicitly considering the hypotheses behind them. This can lead to the erroneous conclusion that, if two data sets are close enough, the models that generated them are similar. In addition, when experimental and simulated data differ to an extent that fails to meet the acceptance criteria, calibration techniques are used to tweak simulation model parameters to reduce the gap between the two types of data. This produces the false expectation that a simulation model will match reality. The methods presented in this paper move away from these strictly data-driven methods for validation and calibration toward more robust, model-driven methods based on causal inference. Causal inference aims to identify the possible mechanisms that might have generated data. Thus, this analysis targets the prediction of the effects when one (or more) of the identified mechanisms are altered. There are many approaches to identify, quantify, and illustrate causal relationships. For the scope of this paper, directed graphs are employed as causal models. If the directed graph lacks cycles, it is known as a directed acyclic graph. A node in such a graph represents an observed data element while a directed edge connecting two nodes represents a causal relationship between two variables. The developed causal methods are designed to extract causal models from simulation models and experimental data. Causal models capture the causal relationships between data elements (e.g., simulated and experimental data). In this context, validation and verification are performed by comparing causal models. The proposed approach does not only inform system analysts on how a simulation model matches real-world data, but also identifies elements of the simulation model that should be revised when discrepancies between simulation and experimental data are observed. Through these causal methods, analysts can identify the portion of the model equation(s) that are behind an edge connecting two variables. Hence, once the structural differences between causal models have been determined, model calibration can occur by changing only those model parameters that impact the identified causal relationships.

97 MATHEMATICS AND COMPUTING↗

VALIDATION, VERIFICATION, AND CALIBRATION THROUGH A CAUSAL LENS

This paper presents an alternative method based on causal inference to perform validation, verification, and calibration of simulation models. While classical validation and verification approaches focus on the identification of the associations between data elements using statistical and machine learning methods, the novel methods in this paper focus instead on the identification of causal relationships between data elements. Statistical and machine learning-based approaches are strictly data-driven, meaning that they provide quantitative comparison measures between datasets without explicitly considering the hypotheses behind them. This can lead to the erroneous conclusion that, if two data sets are close enough, then the models that generated them are similar. In addition, when experimental and simulated data differ to an extent that fails to meet the acceptance criteria, calibration techniques are used to tweak simulation model parameters to reduce the gap between the two types of data. This produces the false expectation that a simulation model will match reality. The methods presented in this paper move away from these strictly data-driven methods for validation and calibration toward more robust, model-driven methods based on causal inference. Causal inference aims to identify the possible mechanisms that might have generated data. Thus, this analysis targets the prediction of the effects when one (or more) of the identified mechanisms are altered. There are many approaches to identify, quantify and illustrate causal relationships. For the scope of this paper, directed graphs are employed as causal models. If the directed graph lacks cycles it is known as a directed acyclic graph (DAG). A node in such a graph represents an observed data element while a directed edge connecting two nodes represents a causal relationship between two variables. The developed causal methods are designed to extract causal models from simulation models and from experimental data. Causal models capture the causal relationships between data elements (e.g., simulated and experimental data). In this context, validation and verification are performed by comparing causal models. The proposed approach does not only inform system analysts on how a simulation model matches real-world data, but also identifies elements of the simulation model that should be revised when discrepancies between simulation and experimental data are observed. Through these causal methods, analysts have a means to identify the portion of the model equation(s) that are behind an edge connecting two variables. Hence, once the structural differences between causal models have been determined, model calibration can occur by changing only those model parameters that impact the identified causal relationships.

97 MATHEMATICS AND COMPUTING↗