Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Glassy Amorphous Polymer (GAP) model is a viscoelastic/plastic model developed at Los Alamos National Laboratory to accurately model a variety of polymers across a wide range of conditions and loading rates, including shock loading. In the present report we introduce and assess this model, newly implemented in the ALEGRA shock and multiphysics code, using a series of verification and application-related validation problems. We describe the mathematical and theoretical formulation of the model, as well as its implementation in ALEGRA, in detail. We provide verification results that assess the model implementation against published computational results, as well as validation results which we compare to existing experimental results when possible. These comparisons instill confidence in the implementation and indicate that the addition of the GAP model to the ALEGRA code provides users with a high-fidelity polymer modeling capability that is capable of recreating complex polymer phenomena.

A novel algorithm for explicit temporal discretization of the variable-density, low-Mach Navier-Stokes equations is presented here. Recognizing there is a redundancy between the mass conservation equation, the equation of state, and the transport equation(s) for the scalar(s) which characterize the thermochemical state, and that it destabilizes explicit methods, we demonstrate how to analytically eliminate the redundancy and propose an iterative scheme to solve the resulting transformed scalar equations. The method obtains second-order accuracy in time regardless of the number of iterations, so one can terminate this subproblem once stability is achieved. Hence, flows with larger density ratios can be simulated while still retaining the efficiency, low cost, and parallelizability of an explicit scheme. The temporal discretization algorithm is used within a pseudospectral direct numerical simulation which extends the method of Kim, Moin, and Moser for incompressible flow [17] to the variable-density, low-Mach setting, where we demonstrate stability for density ratios up to ∼25.7.

Abstract not provided.

Abstract not provided.

We investigate hydrodynamic fluctuations in the flow past a circular cylinder near the critical Reynolds number Rec for the onset of vortex shedding. Starting from the fluctuating Navier-Stokes equations, we perform a perturbation expansion around Rec to derive analytical expressions for the statistics of the fluctuating lift force. Molecular-level simulations using the direct simulation Monte Carlo method support the theoretical predictions of the lift power spectrum and amplitude distribution. Notably, we have been able to collect sufficient statistics at distances Re/Rec-1=O(10-3) from the instability that confirm the appearance of non-Gaussian fluctuations, and we observe that they are associated with intermittent vortex shedding. These results emphasize how unavoidable thermal-noise-induced fluctuations become dramatically amplified in the vicinity of oscillatory flow instabilities and that their onset is fundamentally stochastic.

Photonic computing has the potential to harness the full degrees of freedom (DOFs) of the light field, including the wavelength, spatial mode, spatial location, phase quadrature, and polarization, to achieve a higher level of computing parallelism and scalability than digital electronic processors. While multiplexing using the wavelength and other DOFs can be readily integrated on silicon photonics platforms with compact footprints, conventional mode-division multiplexed (MDM) photonic designs occupy areas exceeding tens to hundreds of microns for a few spatial modes, significantly limiting their scalability. Here, we utilize inverse design to demonstrate an ultracompact photonic computing core that calculates vector dot products based on MDM coherent mixing. Our dot-product core integrates the functionalities of two-mode multiplexers and one multimode coherent mixer within a nominal footprint of 5 μm x 3 μm . We have experimentally demonstrated computing examples on the fabricated dot-product core, including complex number multiplication and motion estimation using optical flow. The compact dot-product core design enables large-scale on-chip integration in a parallel photonic computing primitive cluster for high-throughput scientific computing and computer vision tasks.

A series of reactive-transport models of Enhanced Geothermal Systems (EGS) were constructed using the reactive transport code PFLOTRAN to examine the effect of matrix thermal contraction and mineral dissolution/precipitation on fracture flow in the context of grid cell size and model complexity. It was found that for thermal drawdown at production well, the impact of fracture zone grid cell size is negligible.

Abstract not provided.

In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.

Abstract not provided.

Abstract not provided.

I am going to present my work on multifidelity timeseries models at MS&T in Pittsburgh. We develop efficient machine learning methodologies to accelerate time-series predictions from a hierarchy of complex physics-based models.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologi-cally consistent model that honors key stratigraphic horizons and 3D fault surfaces inter-preted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection inter-val to the basement is essential to explain the destabilization of the regions where micro-seismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO2 sequestration operations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In this project we considered the initial stages of helium bubble nucleation via the proposed mechanism of self-interstitial atom nucleation. By calculating the energy barrier to self-interstitial atom nucleation in a range of Fe-Ni-Cr alloys we identified the most important energetic contributions to the phenomenon: the Frenkel-pair energy barrier in the absence of helium and the difference of insertion energy for a He cluster into a perfect lattice and vacancy. From this observation, we developed a simple model of helium-assisted self-interstitial atom nucleation.

Although there are increasing numbers of distributed energy resources (DERs) and microgrids being deployed, current IEEE and utility standards generally strictly limit their interconnection inside secondary networks. Secondary networks are low-voltage meshed (non-radial) distribution systems that create redundancy in the path from the main grid source to each load. This redundancy provides a high level of immunity to disruptions in the distribution system, and thus extremely high reliability of electric power service. There are two main types of secondary networks, called grid and spot secondary networks, both of which are used worldwide. In the future, primary networks in distribution systems that might include looped or meshed distribution systems at the primary-voltage (medium-voltage) level may also become common as a means for improving distribution reliability and resilience.

We build upon recent work on the use of machine-learning models to estimate Hamiltonian parameters using continuous weak measurement of qubits as input. We consider two settings for the training of our model: (1) supervised learning, where the weak-measurement training record can be labeled with known Hamiltonian parameters, and (2) unsupervised learning, where no labels are available. The first has the advantage of not requiring an explicit representation of the quantum state, thus potentially scaling very favorably to a larger number of qubits. The second requires the implementation of a physical model to map the Hamiltonian parameters to a measurement record, which we implement using an integrator of the physical model with a recurrent neural network to provide a model-free correction at every time step to account for small effects not captured by the physical model. We test our construction on a system of two qubits and demonstrate accurate prediction of multiple physical parameters in both the supervised context and the unsupervised context. We demonstrate that the model benefits from larger training sets, establishing that it is "learning,"and we show robustness regarding errors in the assumed physical model by achieving accurate parameter estimation in the presence of unanticipated single-particle relaxation.

Abstract not provided.

Abstract not provided.

Rapid, time-varying, three-dimensional physics underpin numerous engineering challenges. Often, these physics occur within opaque environments, internal to a component, severely limiting applicable diagnostics. Development of novel diagnostics is necessary to understand and predict transient three-dimensional (3D) phenomena within opaque environments. This report highlights progress in four key areas leading to advancements in high-speed X-ray radiography and tomography. The first area is enabling MHz-rate imaging of energetics at the Advanced Photon Source at Argonne National Laboratory. The second is modeling a high-flux, rotating-anode X-ray source to understand the heat loads on the anode. The third effort was to develop a novel reconstruction algorithm that is validated by ground experimental tomography data and synthetic tomography data. The fourth is the development of a novel approach to two-color X-ray imaging.

The SIREN (Scaling Ion Traps by Requiring iNnovative heterogenous integration) project explores the feasibility of heterogeneous integration (HI) as a transformative approach to scaling ion traps, a critical technology for advancing quantum computers and atomic clocks. Traditional ion trap architectures face significant challenges in scalability due to limitations in optical access, fabrication techniques, and material constraints. SIREN addresses these challenges by leveraging HI, which combines different materials and fabrication processes to create more complex and efficient ion trap structures. HI integrated structures can be manufactured without compromising the process to maintain compatibility to ion traps. This project focuses on integrating a separately fabricated waveguide with a fully functional ion trap. The respective alignment between the pieces needs to be accurate to less than 2 µm to ensure that the light from the waveguide can overlap with the trapping region. The fine alignment must also be maintained through an ultra-high vacuum bake, a critical step in preparing an ion trap experiment. The project's outcomes suggest that heterogeneous integration is a promising pathway for overcoming current scalability barriers, paving the way for the next generation of quantum technologies. SIREN's findings contribute significantly to the field, offering a scalable solution that could accelerate the development of practical quantum computers and highly accurate atomic clocks.

Abstract not provided.

The adoption of machine learning approaches for gamma-ray spectroscopy has received considerable attention in the literature. Many studies have investigated the deployment of various algorithm architectures to a specific task. However, little attention has been afforded to the development of the datasets leveraged to train the models. Such training datasets typically span a set of environmental or detector parameters to encompass a problem space of interest to a user. Variations in these measurement parameters will also induce fluctuations in the detector response, including expected pile-up and ground scatter effects. Fundamental to this work is the understanding that 1) the underlying spectral shape varies as the measurement parameters change and 2) the statistical uncertainties associated with two spectra impact their level of similarity. While previous studies attribute some arbitrary discretization to the measurement parameters for the generation of their synthetic training data, this work introduces a principled methodology for efficient spectral-based discretization of a problem space. A signal-to-noise ratio (SNR) respective spectral comparison measure and a Gaussian Process Regression (GPR) model are used to predict the spectral similarity across a range of measurement parameters. This innovative approach effectively showcased its capability by dividing a problem space, ranging from 5 cm to 100 cm standoff distances and 5 μCi–100 μCi of 137Cs, into three unique combinations of measurement parameters. The findings from this work will aid in creating more robust datasets, which incorporate many possible measurement scenarios, reduce the number of required experimental test set measurements, and possibly enable experimental training data collection for gamma-ray spectroscopy.

Presentation at IEEE Brain Discovery and Neurotechnology Workshop https://brain.ieee.org/2024-ieee-brain-discovery-neurotechnology-workshop/

Modern concepts for next generation pulsed power (NGPP) are slated to deliver up to ten times the energy of Z today. An increase of this magnitude is concerning insofar that Z currently exhibits sizable amounts of inner magnetically insulated transmission line (MITL) loss current on the order of 5-10%. Loss phenomenon in these systems are complex and electrode heating and subsequent thermal desorption are a leading cause. Rapid heat-driven thermal desorption of contaminants scales as the square of the current. Therefore, even a modest doubling of drive current would yield an ~ 4X in non-linear surface electrode heating, quickening thermal desorption-based current loss. Exacerbating these physics is a current inability to measure ultra fast heating rates (>20°C/ns), which are paramount to benchmarking and code validation critical to NGPP design – as an empirical approach is not viable. Therefore, Ultrafast Photoluminescent Surface Heating Optical Thermometry (UP-SHOT) was developed as a new diagnostic for measurement of GHz-scale electrode heating. The discovery of UP-SHOT leveraged expertise in Engineering Science, Material Science, Pulsed-Power, and the Center for Integrated Nanotechnologies. This report includes information on: 1) The preparation of zinc oxide (ZnO) films, characterization, post-deposition treatments 2) Time-resolved photoluminescence at elevated temperatures and thermographic sensitivity

The synthesis of Lanthanum Cobalt Oxide (LCO) via Sol-Gel method provides a potentially low-cost method of production for a P-type photocatalytic. LCO was prepped from Lanthanum and Nitrate Precursors in Aqueous solution. After a significant amount of water is evaporated, the solution is deposited onto SiO2 substrate via spin-deposition method. After annealing, the samples produce thin film LCO that display thickness of sub-500 nanometers. The samples material profile is confirmed by both Raman spectroscopy and X-Ray Diffraction spectroscopy (XRD). Additionally, two-point probing tested the conductivity of the samples. The thin film samples display characteristics of a P-type photocatalytic and may potentially be used in the formation of a P-N junction for user in water-splitting applications.

Abstract not provided.

Abstract not provided.

This paper details a computational framework to produce automated, graphical workflows, and how this framework can be deployed to support complex modeling problems like those in nuclear engineering. Key benefits of the framework include: automating previously manual workflows; intuitive construction and communication of workflows through a graphical interface; and automated file transfer and handling for workflows deployed across heterogeneous computing resources. This paper demonstrates the framework's application to probabilistic post-closure performance assessment of systems for deep geologic disposal of nuclear waste. However, the framework is a general capability that can help users running a variety of computational studies.

This report summarizes the research and capabilities developed as part of the project “Analysis and Validation of Critical Infrastructure Models using Model Order Reduction” (CI-MOR) LDRD project. CI-MOR research enables the solution of large, complex optimization models that naturally arise in national security challenges involving critical infrastructures. Specifically, CI-MOR researchers developed methods to (1) rigorously approximate complex, nonlinear optimization formulations, (2) identify alternative near-optimal solutions, (3) accelerate optimization workflows used for complex applications, and (4) rigorously integrate domain knowledge in stochastic-process models. This report provides an overview of the research done in CI-MOR, and we describe application exemplars used to illustrate CI-MOR capabilities. Furthermore, we describe the software developed by CI-MOR that researchers can leverage to analyze new applications.

Abstract not provided.

In this study, we address the challenge of enhancing image quality and spatial resolution in computed tomography (CT) imaging by introducing simulation and fabrication of high aspect ratio, point-like transmission targets. Utilizing advanced electroplating techniques, traditionally employed in the fabrication of Through Substrate Via (TSV) interconnects for CMOS circuitry, we successfully embed copper targets within silicon substrates. This method allows us to create high-aspect-ratio features specifically designed for X-ray transmission targets, resulting in micro targets that exhibit a volume increase compared to conventional evaporated surface targets. Furthermore, we present simulation results of the X-ray spectrum generated by these targets, demonstrating their potential to significantly improve both image quality and spatial resolution in CT applications. Our findings suggest that leveraging advanced fabrication techniques can open new avenues for the development of enhanced imaging technologies in medical diagnostics and beyond.

Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ∼8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2 < ML < 4) recorded at local distances (< 150 km), for which depth estimation is especially challenging. We find that differential magnitudes (Mw= ML–Mc) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡ log10 (M0)+3log10 (fc (,and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. We also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Infrasound sensing plays a critical role in the detection and analysis of bolides, offering passive, cost-effective global monitoring capabilities. Key objectives include determining the timing, location, and yield of these events. Achieving these goals requires a robust approach to detect, analyze, and interpret rapidly moving elevated sources such as bolides (also re-entry). In light of advancements in infrasonic methodologies, there is a need for a comprehensive overview of the characteristics that distinguish bolides from other infrasound sources and methodologies for bolide infrasound analysis. This paper provides a focused review of key considerations and presents a unified framework to enhance infrasound processing approaches specifically tailored for bolides. Three representative case studies are presented to demonstrate the practical application of infrasound processing methodologies and deriving source parameters while exploring challenges associated with bolide-generated infrasound. These case studies underscore the effectiveness of infrasound in determining source parameters and highlight interpretative challenges, such as variations in signal period measurements across different studies. Future research should place emphasis on improving geolocation and yield accuracy. This can be achieved through rigorous and systematic analyses of large, statistically significant samples of such events, aiming to resolve interpretative inconsistencies and explore the causes for variability in signal periods and back azimuths. The topic described here is also relevant to space exploration involving planetary bodies with atmospheres, such as Venus, Mars, and Titan.

Abstract not provided.

This report describes an approach to utilizing phasor measurement unit (PMU) data from multiple Intelligent Electronics Devices (IEDs) in a low-voltage network to produce a differential scheme for protecting the medium-voltage feeder and low-voltage network transformers. The proposed protection scheme is designed and prototyped on a real-time automation controller. Its performance is evaluated using real-time controller hardware-in-the-loop simulation. Lab testing results indicate that the proposed protection scheme allows significant distributed energy resources (DER) backfeed and enables selective and fast protection of medium voltage feeders.