SEARCH · Engineering Papers

Results for “RANS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Mechanism of RanGTP priming H2A-H2B release from Kap114 in an atypical RanGTP•Kap114•H2A-H2B complex

Previously, we showed that the nuclear import receptor Importin-9 wraps around the H2A-H2B core to chaperone and transport it from the cytoplasm to the nucleus. However, unlike most nuclear import systems where RanGTP dissociates cargoes from their importins, RanGTP binds stably to the Importin-9•H2A-H2B complex, and formation of the ternary RanGTP•Importin-9•H2A-H2B complex facilitates H2A-H2B release to the assembling nucleosome. It was unclear how RanGTP and the cargo H2A-H2B can bind simultaneously to an importin, and how interactions of the three components position H2A-H2B for release. Here, we show cryo-EM structures of Importin-9•RanGTP and of its yeast homolog Kap114, including Kap114•RanGTP, Kap114•H2A-H2B, and RanGTP•Kap114•H2A-H2B, to explain how the conserved Kap114 binds H2A-H2B and RanGTP simultaneously and how the GTPase primes histone transfer to the nucleosome. In the ternary complex, RanGTP binds to the N-terminal repeats of Kap114 in the same manner as in the Kap114/Importin-9•RanGTP complex, and H2A-H2B binds via its acidic patch to the Kap114 C-terminal repeats much like in the Kap114/Importin-9•H2A-H2B complex. Ran binds to a different conformation of Kap114 in the ternary RanGTP•Kap114•H2A-H2B complex. Here, Kap114 no longer contacts the H2A-H2B surface proximal to the H2A docking domain that drives nucleosome assembly, positioning it for transfer to the assembling nucleosome or to dedicated H2A-H2B chaperones in the nucleus.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Droplet and Particle Methods to Investigate Turbulent Particle Laden Jets

The SARS-CoV-2 pandemic has heightened the interest in particle-laden turbulent jets generated by breathing, talking, coughing and sneezing, and how these can contribute to disease transmission. We present quantitative measurement methods for such flows, while exploring and offering improvements for common shortcomings. We generate jets consisting of either liquid droplets or solid particles in an isothermal, quiescent and electrically isopotential experimental chamber that was constructed to control the effects of ambient forcing on jet behavior. For liquid droplets, we find promise in surface deposition analysis based on fluorescent tracer use. For particles, we explore the performance of commercially available adhesive sampling strips and develop conductive grounded carbon tape based sampling strips. We explore ways in which the smallest of thermal gradients or electrostatic charge issues can affect particle dispersion, and suggest practical methods to address these issues. The developed methods are applied to study the simultaneous deposition of (Formula presented.) 25, 50 and 200 μm solid particles from a particle laden turbulent jet with a mean velocity of 33.2 m/s. The deposition location as a function of particle size was compared to results from a simple numerical RANS model, and illustrates ways in which imprecise initial or boundary conditions can lead to a notable deviation from experimental results. The differences in deposition pattern seen in experimental and numerical results despite a carefully controlled environment and characterized particle ejection indicate the need for a more stringent numerical model validation, especially when studying fate and transport of mid-range (neither purely aerosol or ballistic) sized particles.

42 ENGINEERING↗

DOE OSTI2022

Mechanical loading prediction through accelerometry data during walking and running

ABSTRACT Currently, there is no way to assess mechanical loading variables such as peak ground reaction forces (pGRF) and peak loading rate (pLR) in clinical settings. The purpose of this study was to develop accelerometry‐based equations to predict both pGRF and pLR during walking and running. One hundred and thirty one subjects (79 females; 76.9 ± 19.6 kg) walked and ran at different speeds (2–14 km·h −1 ) on a force plate–instrumented treadmill while wearing accelerometers at their ankle, lower back and hip. Regression equations were developed to predict pGRF and pLR from accelerometry data. Leave‐one‐out cross‐validation was used to calculate prediction accuracy and Bland–Altman plots. Our pGRF prediction equation was compared with a reference equation previously published. Body mass and peak acceleration were included for pGRF prediction and body mass and peak acceleration rate for pLR prediction. All pGRF equation coefficients of determination were above 0.96, and a good agreement between actual and predicted pGRF was observed, with a mean absolute percent error (MAPE) below 7.3%. Accuracy indices from our equations were better than previously developed equations. All pLR prediction equations presented a lower accuracy compared to those developed to predict pGRF. Walking and running pGRF can be predicted with high accuracy by accelerometry‐based equations, representing an easy way to determine mechanical loading in free‐living conditions. The pLR prediction equations yielded a somewhat lower prediction accuracy compared with the pGRF equations.

Veras, Lucas↗

DOE OSTI2019

A real-time machine learning-based disruption predictor in DIII-D

A disruption prediction algorithm, named DPRF (Disruption Prediction using Random Forests), has run in real-time in the DIII-D plasma control system (PCS) over more than 900 discharges. DPRF naturally provides a probability mapping associated to its predictions, i.e. the disruptivity signal, now incorporated in the DIII-D PCS. Here, we report on disruption prediction accomplishments in terms of shot-by-shot performances, by simulating alarms on each discharge as in the PCS framework. Due to the optimised performance metric chosen to evaluate DPRF, we find that almost all disruptive discharges are detected in average with few hundred of milliseconds of warning time, but this comes at a high cost of false alarms produced. Performances are not satisfying ITER requirements, where the success rate has to be higher than 95%, but this is not completely unexpected. DPRF is trained on many years of major disruptions occurring during the flattop phase of the plasma current on DIII-D, but without any differentiation by cause. Furthermore, we find that DPRF is affected by a relatively high fraction of false alarms occurring during the first 500 milliseconds from the flattop onset. This subtle effect, more evident on the discharges where DPRF ran in real-time, can be marginalized by taking specific precautions on the validity range of the predictions, and performances do improve. Even if presently burdened by some limitations, DPRF provides an incredible and novel advantage. Thanks to the feature contribution analysis (e.g., the identification of which signals contributed to triggering an alarm), it is possible to interpret and explain DPRF predictions. It is the first time that such interpretability features is exploited by a disruption predictor: by uncovering the disruption events' causes, a better understanding of disruption dynamics is achieved, and a clear path toward the design of disruption avoidance strategies can be provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Verification of induction zone models for wind farm annual energy production estimation

Fast engineering wake models are the backbone of wind farm annual energy production (AEP) estimators, whereas the addition of induction zone models in existing tools is a more recent response to rising concerns over wind farm blockage associated losses. Here, "blockage" describes the combined induction fields of all wind turbines inside a farm. Unlike the term might suggest, blockage not only reduces flow speeds, but also increases them; for instance along the outer edges of wind farms. Evaluating the overall impact on AEP of these gains and losses necessitates accurate wind farm induction models. Whilst engineering wake models are tuned to predict wind farm performance, existing induction zone models are all derived for accurately predicting the near field induction not the far field value—important for wind farm simulations. This paper presents the induction models implemented in the wind farm simulation tool PyWake, as well as results from novel analytical models and compares their far field predictions to RANS-AD simulations of different turbines. We demonstrate that when including blockage in AEP simulations, the downstream speed-ups need to be included to avoid an unrealistic bias toward AEP loss and that wake expansion significantly impacts induction at rated thrust levels.

17 WIND ENERGY↗

DOE OSTI2026

EPICS for small-scale laboratories with Python soft IOCs

While the Experimental Physics and Industrial Control System (EPICS) is widely used at large laboratories for slow controls and instrumentation, the deployment of a full EPICS installation can be difficult, with a steep learning curve to new users. Taking advantage of the pythonSoftIOC module, we developed an EPICS slow controls implementation for Jefferson Lab's Hall B cryotarget written entirely in Python and based on software IOCs that communicate with instruments over Ethernet. Here, this system ran successfully, interfacing with Jefferson Lab's full EPICS network, and we offer it as an example of the capabilities of pythonSoftIOC to build lightweight, yet robust and flexible instrumentation platforms that would be easily adapted for use at a small-scale laboratory. University groups can use these examples to build complete slow controls systems, from device communication to data archiving and display, using open-source, mature EPICS tools and student-friendly Python as an alternative to expensive and proprietary systems such as LabVIEW.

Computing↗

DOE OSTI2025

Plasma phosphorylated tau217 strongly associates with memory deficits in the Alzheimer’s disease spectrum

Abstract Plasma phosphorylated tau (p-tau) biomarkers open unprecedented opportunities for identifying carriers of Alzheimer’s disease pathophysiology in early disease stages using minimally invasive techniques. Plasma p-tau biomarkers are believed to reflect tau phosphorylation and secretion. However, it remains unclear to what extent the magnitude of plasma p-tau abnormalities reflects neuronal network disturbance in the form of cognitive impairment. To address this question, we included 103 cognitively unimpaired elderly and 40 cognitively impaired, amyloid-β-positive individuals from the TRIAD cohort, in addition to 336 cognitively unimpaired and 216 cognitively impaired, amyloid-β-positive older adults from the BioFINDER-2 cohort. Participants had tau PET scans, amyloid PET scans or amyloid CSF, p-tau217, p-tau181 and p-tau231 blood measures, structural T1-MRI and cognitive assessments. In this cross-sectional study, we used regression models and correlation analyses to assess the relationship between plasma biomarkers and cognitive scores. Furthermore, we applied receiver operating characteristic curves to assess cognitive impairment across plasma biomarkers. Finally, we categorized participants into amyloid (A), p-tau (T1) and tau PET (T2) positive (+) or negative (−) profiles and ran non-parametric comparisons to assess differences across cognitive domains. We found that plasma p-tau217 was more associated with cognitive performance than p-tau181 and p-tau231 and that this relationship was particularly strong for memory scores (TRIAD: βp-tau217 = −0.53, βp-tau181 = −0.35 and βp-tau231 = −0.24; BioFINDER-2: βp-tau217 = −0.52, βp-tau181 = −0.24 and βp-tau231 = −0.29). Associations in amyloid-β-positive participants resembled these results, but other cognitive scores also showed strong associations in cognitively impaired individuals. Moreover, plasma p-tau217 outperformed plasma p-tau181 and plasma p-tau231 in identifying memory impairment (area under the curve values for TRIAD: p-tau217 = 0.86, p-tau181 = 0.77 and p-tau231 = 0.75; and for BioFINDER-2: p-tau217 = 0.86, p-tau181 = 0.76 and p-tau231 = 0.81) and in identifying executive function impairment only in the BioFINDER-2 cohort (p-tau217 = 0.82, p-tau181 = 0.76 and p-tau231 = 0.76). Lastly, we showed that subtle memory deficits were present in A+T1+T2− participants for plasma p-tau217 (P = 0.007) and plasma p-tau181 (P = 0.01) in the TRIAD cohort and for all biomarkers across cognitive domains in A+T1+T2− and A+T1+T2− individuals (P < 0.001 in all) in the BioFINDER-2 cohort. The A+T1+T2− individuals showed cognitive deficits in both cohorts (P < 0.001 in all). Together, our results suggest that plasma p-tau217 stands out as a biomarker capable of identifying memory deficits attributable to Alzheimer’s disease and that memory impairment certainly occurs in amyloid-β- and plasma p-tau-positive individuals who have no significant amounts of tau in the neocortex.

Neurosciences & Neurology↗

DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗

DOE OSTI2023

Measurement of the branching fraction and search for CP violation in D 0 → K $^{0}_{S}$ K $^{0}_{S}$ ⁢π + ⁢π - decays at Belle

We measure the branching fraction for the Cabibbo-suppressed decay D 0 → K $^{0}_{S}$ K $^{0}_{S}$ ⁢π + ⁢π - and search for CP violation via a measurement of the CP asymmetry A CP as well as the T-odd triple-product asymmetry a$^{T}_{CP}$. We use 922 fb -1 of data recorded by the Belle experiment, which ran at the KEKB asymmetric-energy e + ⁢e - collider. The branching fraction is measured relative to the Cabibbo-favored normalization channel D 0 → K $^{0}_{S}$ ⁢π + ⁢π - ; the result is $\mathscr{B}$⁡( D 0 → K $^{0}_{S}$ K $^{0}_{S}$ ⁢π + ⁢π - )=[4.79±0.08⁢(stat)±0.10⁢(syst)±0.31⁢(norm)]×10 -4 , where the first uncertainty is statistical, the second is systematic, and the third is from uncertainty in the normalization channel. We also measure A CP =[-2.51±1.44⁢(stat)$^{+0.11}_{-0.10}$⁢(syst)]%, and a$^{T}_{CP}$=[-1.95±1.42⁢(stat)$^{+0.14}_{-0.12}⁢$(syst)]%. These results show no evidence of CP violation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of three Rayleigh-Taylor mixing configurations with gravity reversal

In this work, high-fidelity large-eddy simulation (LES) is performed of Rayleigh-Taylor (RT) mixing in three different configurations involving gravity reversal. In each configuration, LES results are compared with one-dimensional Reynolds-averaged Navier-Stokes (RANS) results, and a deficiency in a commonly used transport equation for the mass-flux velocity, $a_j$, is identified. In the first configuration, a classical two-component RT mixing layer is allowed to develop before it is subjected to rapid acceleration reversal. In the second configuration, a three-component RT mixing layer with an intermediate density layer is allowed to develop before being subjected to rapid acceleration reversal. Finally, in the third configuration, a light layer is interposed between two heavy layers; in this configuration, only one interface is RT-unstable at a time as it undergoes rapid acceleration reversal. In all cases, a commonly used buoyancy production closure in the $a_j$ transport equation is shown to lead to significant over-prediction of mixing layer growth after gravity reversal. An alternative formulation for this closure is then presented which is shown to more accurately capture the stabilization effect of gravity reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Analysis of real-fluid thermodynamic effects on turbulent statistics in transcritical channel flows

Wall-bounded turbulence at high-pressure transcritical conditions with intense density fluctuations are encountered in many technical applications. In this study, we analyze the turbulent energy transport in transcritical channel flows specifically focusing on dissipation rate, turbulent kinetic energy budgets, heat fluxes, and momentum-fluctuation statistics; results from this analysis are used to guide the development of turbulent scaling laws. Furthermore, we find that the dissipation rate of turbulent kinetic energy is dominated by the enstrophy in the logarithmic layer, and the fluctuating viscosity results in the reduced tilting of the vortical structures and the attenuation of streamwise vorticity in the near-wall layer; the fluctuating viscosity attenuates the dissipation rate by reducing the shear strain and the enstrophy production. Local equilibrium of the turbulent kinetic energy exists in the logarithmic layer. We show that the real-fluid thermodynamic effects significantly change the turbulent heat flux correlated with the sweep and the ejection events; the density changes alter the turbulent transport and result in noticeable magnitudes of density-fluctuation-related momentum-fluctuation statistics. From these results, scaling laws for the turbulent length scales and turbulent kinetic energy budgets are proposed, thereby contributing to improvement of the wall models in large-eddy simulations and Reynolds-averaged Navier-Stokes (RANS) simulations.

42 ENGINEERING↗

DOE OSTI2021

Systematic Operating Temperature Differences Between Al-BSF, PERC, and PERT-With-Optimized-Rear-Reflector Solar Mini-Modules Due to Rear Reflectance

Reflecting sub-bandgap light from photovoltaic modules has the potential to improve lifetime energy generation of fielded systems by reducing operating temperature. In this article, the temperature of fielded aluminum back-surface field (Al-BSF) and passivated emitter and rear contact (PERC) mini-modules was monitored every 5 minutes for 75 days along with corresponding meteorological data. Additionally, passivated emitter rear totally diffused (PERT) mini-modules with high-performance sub-bandgap rear reflectors were tested and compared to the state-of-the-art industrial modules. These reflectors consisted of a >300-nm-thick silicon dioxide nanoparticle film with a low refractive index. The impact of reflectance on measured operating temperature was isolated with a previously developed thermal model and quantified as the reflectance-induced median temperature difference between each tested module at representative outdoor conditions (1000 W m -2 , 25 degrees C ambient temperature, and 1.43 m s -1 wind speed). We found that, because of their reflectance differences, PERC modules ran systematically cooler than Al-BSF modules by 1.0 degrees C, whereas the PERT-with-optimized-rear-reflector systematically operated 1.4 degrees C cooler than the Al-BSF module and 0.4 degrees C than the PERC module. We also found that the rear reflector provided the greatest temperature benefit during periods of highest irradiance.

14 SOLAR ENERGY↗

DOE OSTI2025

Traffic Control via Connected and Automated Vehicles (CAVs): An Open-Road Field Experiment with 100 CAVs

The CIRCLES project aims to reduce instabilities in traffic flow, which are naturally occurring phenomena due to human driving behavior. Also called “phantom jams” or “stop-and-go waves,” these instabilities are a significant source of wasted energy. Toward this goal, the CIRCLES project designed a control system, referred to as the MegaController by the CIRCLES team, that could be deployed in real traffic. Our field experiment, the MegaVanderTest (MVT), leveraged a heterogeneous fleet of 100 longitudinally controlled vehicles as Lagrangian traffic actuators, each of which ran a controller with the architecture described in this article. The MegaController is a hierarchical control architecture that consists of two main layers. The upper layer is called the Speed Planner and is a centralized optimal control algorithm. It assigns speed targets to the vehicles, conveyed through the LTE cellular network. The lower layer is a control layer, running on each vehicle. It performs local actuation by overriding the stock adaptive cruise controller, using the stock onboard sensors. The Speed Planner ingests live data feeds provided by third parties as well as data from our own control vehicles and uses both to perform the speed assignment. The architecture of the Speed Planner allows for the modular use of standard control techniques, such as optimal control, model predictive control (MPC), kernel methods, and others. The architecture of the local controller allows for the flexible implementation of local controllers. Corresponding techniques include deep reinforcement learning (RL), MPC, and explicit controllers. Depending on the vehicle architecture, all onboard sensing data can be accessed by the local controllers or only some. Likewise, control inputs vary across different automakers, with inputs ranging from torque or acceleration requests for some cars to electronic selection of adaptive cruise control (ACC) setpoints in others. The proposed architecture technically allows for the combination of all possible settings proposed previously, that is {Speed Planner algorithms} × {local Vehicle Controller algorithms} × {full or partial sensing} × {torque or speed control}. As a result, most configurations were tested throughout the ramp up to the MegaVandertest (MVT).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2019

Quantum computation of nanosheets in a background magnetic field for external control of nanosystems

We ran quantum simulations on lattices of interacting fermions in a background magnetic field [1]. We created discretized Hamiltonian matrices for each system and pre-processed them so that they could be read by the quantum simulator implemented in IBM's Qiskit package [2]. The simulator uses the variational quantum eigensolver (VQE) algorithm [3] to approximate the ground state energy of the input Hamiltonian, which is then compared to the classically-calculated exact value. This classical evaluation is possible because the systems studied are still within the grasp of classical computers. We also investigated ways to simplify the system before the pre-processing step. We found that in almost all cases these simplifications made our quantum calculations converge to better answers using simpler starting conditions. We are exploring other ways to interface with IBM's quantum computing software so that we no longer need to explicitly define the Hamiltonian matrices for our systems, a problem which will itself soon become intractable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Scalable Knowledge Graph Analytics at 136 Petaflop/s

We are motivated by newly proposed methods for data mining large-scale corpora of scholarly publications, such as the full biomedical literature, which may consist of tens of millions of papers spanning decades of research. In this setting, analysts seek to discover how concepts relate to one another. They construct graph representations from annotated text databases and then formulate the relationship-mining problem as one of computing all-pairs shortest paths (APSP), which becomes a significant bottleneck. In this context, we present a new high-performance algorithm and implementation of the Floyd-Warshall algorithm for distributed-memory parallel computers accelerated by GPUs, which we call DSNAPSHOT (Distributed Accelerated Semiring All-Pairs Shortest Path). For our largest experiments, we ran DSNAPSHOT on a connected input graph with millions of vertices using 4, 096nodes (24,576GPUs) of the Oak Ridge National Laboratory's Summit supercomputer system. We find DSNAPSHOT achieves a sustained performance of 136×1015 floating-point operations per second (136petaflop/s) at a parallel efficiency of 90% under weak scaling and, in absolute speed, 70% of the best possible performance given our computation (in the single-precision tropical semiring or “min-plus” algebra). Looking forward, we believe this novel capability will enable the mining of scholarly knowledge corpora when embedded and integrated into artificial intelligence-driven natural language processing workflows at scale.

Kannan, Ramakrishnan {ramki}↗

DOE OSTI2024

Toward Automated Detection of Portability Bugs in Kokkos Parallel Programs

Performance-portable programming frameworks provide abstractions for parallel execution to allow easily porting an application to multiple backend programming models, such as CUDA, HIP, and OpenMP. However, programs may still have portability bugs that manifest only on specific backends. Traditional testing is ineffective in discovering these bugs, as it would require concrete execution on all supported hardware configurations for a potentially infinite set of inputs. To mitigate this issue, we focused on a specific programming framework, Kokkos, and identified several categories of common portability bugs. We then developed Klokkos, a static analysis approach based on symbolic execution that can run on commodity hardware, before execution on supercomputers. As a proof-of-concept, we ran Klokkos on examples encoding the identified bugs. Our results show that Klokkos is effective, efficient, and precise: it detected all the considered bugs, quickly, and without any false positives. Although preliminary, our results motivate further research and development in this direction.

Kale, Vivek↗

DOE OSTI2023

Material Identification From Radiographs Without Energy Resolution

We propose a method for performing material identification from radiographs without energy-resolved measurements. Material identification has a wide variety of applications, including in biomedical imaging, nondestructive testing, and security. While existing techniques for radiographic material identification make use of dual energy sources, energy-resolving detectors, or additional (e.g., neutron) measurements, such setups are not always practical— requiring additional hardware and complicating imaging. We tackle material identification without energy resolution, allowing standard X-ray systems to provide material identification information without requiring additional hardware. Assuming a setting where the geometry of each object in the scene is known and the materials come from a known set of possible materials, we pose the problem as a combinatorial optimization with a loss function that accounts for the presence of scatter and an unknown gain and propose a branch and bound algorithm to efficiently solve it. We present experiments on both synthetic data and real, experimental data with relevance to security applications— thick, dense objects imaged with MeV X-rays. We show that material identification can be efficient and accurate, for example, in a scene with three shells (two copper, one aluminum), our algorithm ran in six minutes on a consumer-level laptop and identified the correct materials as being among the top 10 best matches out of 8,000 possibilities.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Lightning Current Propagation in Electrical Conduit

Investigation of lightning strikes to conductors ran through long spans of rigid steel conduit was performed. An overdamped-exponential current waveform with controlled peaks and rise rates was used to inject simulated lightning strikes. The impact of the length of wire, length of conduit, grounding location/s, and load type was investigated. Breakdown of 600 V, 12 AWG, THHN insulated wire (3.23 mm OD, 2 mm conductor diameter) was observed for voltages above 45 kV. The presence of resistive loads (between wire and ground) in excess of 20 ohms, or current rise times in excess of 5 kA/μs were found to consistently produce breakdown between wire and conduit. Practical power circuit elements such as outlets and splices were found to breakdown at voltage levels much below the wire insulation failure threshold.

47 OTHER INSTRUMENTATION↗