SEARCH · Engineering Papers
Results for “impact modeling”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A Semi-Empirical Density Law for Ternary, Homogeneous PuCl 3 /HCl/H 2 O Solutions
Nuclear material operations pose unique hazards that are not encountered in other chemical, energy, or manufacturing industries. One of these hazards is the potential for a nuclear criticality accident when handling fissile isotopes such as 235 U and 239 Pu. These hazards are particularly high when fissile material is dissolved in solution as the neutron behaviors of the system can change rapidly with the physical and chemical changes accessible in solution. Current estimates of solution density used for criticality safety are outdated and hinder fissile material handling. Developing new estimates for these safety calculations requires experimental characterization and the derivation of empirical density models. We have derived a density law describing PuCl 3 /HCl/H 2 O solutions from experimental data characterizing solution density. Density data was treated using a Pitzer-derived eight-parameter equation, defining density as a function of analyte concentrations, temperature, and interactions between these variables. The model is predictive across the concentration and temperature ranges from which it was derived. The potential effects of varying oxidation states of plutonium, which are easily accessible in aqueous media, on the bulk solution density of the ternary system were also investigated. The resulting Pitzer-derived density law was applied to a nuclear criticality safety model, and the impact of the experimental characterization of solution density relative to previous estimates was demonstrated to be significant and suggest that the current approach to estimating density in nuclear criticality safety calculations may lead to overly conservative controls.
Off-Lattice Markov Chain Monte Carlo Simulations of Mechanically Driven Polymers
Here, we develop off-lattice simulations of semiflexible polymer chains subjected to applied mechanical forces by using Markov Chain Monte Carlo. Our approach models the polymer as a chain of fixed length bonds, with configurations updated through adaptive nonlocal Monte Carlo moves. This proposed method enables precise calculation of a polymer’s response to a wide range of mechanical forces, which traditional on-lattice models cannot achieve. Our approach has shown excellent agreement with theoretical predictions of persistence length and end-to-end distance in quiescent states as well as stretching distances under tension. Moreover, our model eliminates the orientational bias present in on-lattice models, which significantly impacts calculations such as the scattering function, a crucial technique for revealing the polymer conformation.
Assessing Implicit and Explicit Polarizable Solvation Models for Nuclear–Electronic Orbital Systems: Quantum Proton Polarization and Solvation Energetics
Accurate simulations of many chemical processes require the inclusion of both nuclear quantum effects and a solvent environment. The nuclear–electronic orbital (NEO) approach, which treats electrons and select nuclei quantum mechanically on the same level, combined with a polarizable continuum model (PCM) for the solvent environment, addresses this challenge in a computationally practical manner. In this work, the NEO-PCM approach is extended beyond the IEF-PCM (integral equation formalism PCM) and C-PCM (conductor PCM) approaches to the SS(V)PE (surface and simulation of volume polarization for electrostatics) and ddCOSMO (domain decomposed conductor-like screening model) approaches. IEF-PCM, SS(V)PE, C-PCM, and ddCOSMO all exhibit similar solvation energies as well as comparable nuclear polarization within the NEO framework. The calculations show that the nuclear density does not leak out of the molecular cavity because it is much more localized than the electronic density. Finally, the polarization of quantized protons is analyzed in both continuum solvent and explicit solvent environments described by the polarizable MB-pol model, illustrating the impact of specific hydrogen-bonding interactions captured only by explicit solvation. Furthermore, these calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.
Tracking Local pH Dynamics during Water Electrolysis via In-Line Continuous Flow Raman Spectroscopy
The performance of electrochemical devices, which play a critical role in decarbonization efforts, is often governed by proton-coupled electron transfer reactions at the electrode–electrolyte interface. These reactions are highly sensitive to the complex and dynamic microenvironment present at the electrode surface. However, characterizing this environment─particularly monitoring interfacial pH and its evolution under reaction conditions─remains challenging, necessitating the development of advanced analytical tools. Here, in this study, we introduce in-line continuous flow Raman spectroscopy (CFRS) as a spectroelectrochemical platform for quantifying interfacial pH swings generated during water-splitting. By monitoring phosphate ion speciation and controlling the hydrodynamics with a flow cell, we measure pH swings as a function of current density, flow rate, and distance from the electrode. Comparison with theoretical models reveals the impact of bulk pH, boundary layer thickness, and bubble dynamics at high current densities. Collectively, these findings establish CFRS as a platform for quantitatively investigating pH dynamics, offering critical insights for advancing electrochemical energy conversion technologies.
Microbial Community Composition in Deep‐Subsurface Reservoir Fluids Reveals Natural Interwell Connectivity
Abstract The identification of natural fractures and the wells they connect is crucial for the development of geological reservoirs because it may have an important impact on reservoir model construction and hydraulic fracture propagation. In this study we investigated the use of a novel data source, the microbial community composition in the reservoir formation fluids, for identification of interwell connectivity caused by natural fractures. We verified this concept at a newly developed mesoscale enhanced geothermal system (EGS) field testbed located 4,850 ft (1,478 m) beneath ground surface at the Sanford Underground Research Facility in Lead, SD. Fluids produced at or near the EGS testbed were sampled and subjected to high‐throughput 16S rRNA gene amplicon sequencing to analyze the microbial community profile therein. Despite the typically substantial heterogeneity across the community profiles of samples spatially distributed (10 m to 1.9 km apart) throughout the site, samples from two wells at the EGS testbed showed highly similar microbial community composition, suggesting the two wells intersected the same natural fracture. This evidence of natural connectivity between the two wells at the EGS testbed was later corroborated by core log analysis and sewer camera surveys into the boreholes. Besides the field case described in this study, microbial community composition as a reservoir diagnostic tool would be applicable in a much broader context such as unconventional hydrocarbon exploitation, groundwater reservoir characterization, and environmental remediation, adding valuable “hard” data capable of pinpointing the origins of fluids unambiguously.
Assessing the Impacts of Recent Crop Expansion on Water Quality in the Missouri River Basin Using the Soil and Water Assessment Tool
The Missouri River Basin (MORB) has experienced a resurgence of grassland conversion to crop production, which raised concerns on water quality. We applied the Soil and Water Assessment Tool (SWAT) to address how this conversion would impact water quality. We designed three crop production scenarios representing conversion of grassland to: (a) continuous corn; (b) corn/soybean rotation; and (c) corn/wheat rotation to assess the impact. The SWAT model results showed: (a) the lower MORB produced high total nitrogen (TN) and total phosphorus (TP) load before conversion (baseline) due mainly to high precipitation and high agricultural activity; (b) the greatest percentage increases of TN and TP occurred in the North and South Dakotas, coinciding with the highest amount of grassland conversion to cropland; and (c) grassland conversion to continuous corn resulted in the greatest increase in TN and TP loads, followed by conversion to corn/soybean and then conversion to corn/wheat. Although the greatest percentage increases of TN and TP occurred in the North and South Dakotas, these areas still contributed relatively low TN and TP to total basin loads after conversion. However, watersheds, predominantly in the lower MORB continued to be “hotspots” that contributed the greatest amounts of TN and TP to the total basin loads—driven by a combination of grassland conversion, high precipitation, and loading from pre-existing cropland. At the watershed outlet, the TN and TP loads were increased by 6.4% (13,800 t/yr) and 8.7% (3,400 t/yr), respectively, during the 2008–2016 period for the conversion to continuous corn scenario.
High pressure suppression of plasticity due to an overabundance of shear embryo formation
Abstract High pressure shear band formation is a critical phenomenon in energetic materials due to its influence on both mechanical strength and mechanochemical activation. While shear banding is known to occur in a variety of these materials, the governing dynamics of the mechanisms are not well defined for molecular crystals. We conduct molecular dynamics simulations of shock wave induced shear band formation in the energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) to assess shear band nucleation processes. We find, that at high pressures, the initial formation sites for shear bands, “embryos”, form in excess and rapidly lower deviatoric stresses prior to shear band formation and growth. This results in the suppression of plastic deformation. A local cluster analysis is used to quantify and contrast this mechanism with a more typical shear banding seen at lower pressures. These results demonstrate a mechanism that is reversible in nature and that supersedes shear band formation at increased pressures. We anticipate that these results will have a broad impact on the modeling and development of high-strain rate application materials such as those for high explosives and hypersonic systems.
Quantum-classical embedding via ghost Gutzwiller approximation for enhanced simulations of correlated electron systems
Simulating correlated materials on present-day quantum hardware remains challenging due to limited quantum resources. Quantum embedding methods offer a promising route by reducing computational complexity through the mapping of bulk systems onto effective impurity models, allowing more feasible simulations on pre- and early-fault-tolerant quantum devices. Here, this work develops a quantum-classical embedding framework based on the ghost Gutzwiller approximation to enable quantum-enhanced simulations of ground-state properties and spectral functions of correlated electron systems. Circuit complexity is analyzed using an adaptive variational quantum algorithm on a statevector simulator, applied to the infinite-dimensional Hubbard model with increasing ghost mode numbers from 3 to 5, resulting in circuit depths growing from 16 to 104. Noise effects are examined using a realistic error model, revealing significant impact on the spectral weight of the Hubbard bands. To mitigate these effects, the Iceberg quantum error detection code is employed, achieving up to 40% error reduction in simulations. Finally, the accuracy of the density matrix estimation and the derived spectral function is benchmarked on IBM and Quantinuum quantum hardware, featuring distinct qubit-connectivity and employing multiple levels of error mitigation techniques.
Convergence acceleration in machine learning potentials for atomistic simulations
Machine learning potentials (MLPs) for atomistic simulations have an enormous prospective impact on materials modeling, offering orders of magnitude speedup over density functional theory simulations without appreciably sacrificing accuracy of material property prediction.
Mechanisms of shape transfer and preheating in indirect-drive double shell collisions
Implosions of Hohlraum-driven double shell targets as an alternative inertial confinement fusion concept are underway at the National Ignition Facility. The double shell system relies on a series of energy transfer processes starting from thermal x-ray absorption by the outer shell, followed by collisional transfer of kinetic energy to a heavy metal inner shell, and finally, conversion to the internal energy of the deuterium-tritium fuel. During each of these energy transfer stages, low-mode asymmetries can act to reduce the ideal transfer efficiency degrading double shell performance. Mechanisms, such as hard x-ray preheat from the Hohlraum, not only decrease the efficiency of kinetic energy transfer but may also be a source of low-mode asymmetry. In this article, we evaluate the shape transfer processes through the time of shell collision using two-dimensional integrated Hohlraum and capsule computations. We find that the dominant mode of the shape transfer is well described using a “radial impulse” model from the shape of the foam pressure reservoir. To evaluate the importance of preheat on inner shell shape, we also report on first measurements of Au L-shell preheat asymmetry in a double shell with a tungsten pusher. These measurements showed a 65% higher preheat velocity at the pole of the capsule relative to the equator. We also found that the experiments provided rigorous constraints by which to test the Hohlraum model settings that impact the amount and symmetry of Au L-shell preheat via the plasma conditions inside the outer cone Au bubble.
New heat flux model for non-axisymmetric divertor infrared structures
A convective heat flux model for perturbed plasmas, based on guiding center ion drift in vacuum elds (A. Wingen, et al., Phys. Plasmas 21 (2014) 012509), has been updated. The old model only considered ion heat flux, while here also electron heat flux is included. The updated model predicts divertor heat flux distributions in non-axisymmetric (3D) plasmas with applied Resonant Magnetic Perturbation (RMP) fields, and includes electric scalar potentials. It is found that a radial electric field in the near Scrape-o Layer (SOL) can considerably shift the footprints toroidally, leading to a smearing out eff ect of the incident heat flux, while a simple model for sheath potential has little impact on footprints. Various approaches to model electron heat flux are studied. A convective electron model, based on collisionless free streaming, is found to yield the best agreement with measurements, while a conductive model requires a at temperature gradient inside lobes to yield acceptable peak heat ux values. A heuristic heat flux layer approach, based on a fixed layer width also requires a limited heat flux inside the last closed flux surface (LCFS); by selecting various locations of the LCFS, the results of the conductive or convective model can be recovered respectively. The sum of ion and electron heat fluxes, both obtained by the convective model, is compared to experimental data for multiple time slices in DIII-D. Strike point splitting is observed with peak heat fluxes and layer widths that compare well to infrared camera (IR) measurements.
Small-scale signatures of primordial non-Gaussianity in k-nearest neighbour cumulative distribution functions
ABSTRACT Searches for primordial non-Gaussianity in cosmological perturbations are a key means of revealing novel primordial physics. However, robustly extracting signatures of primordial non-Gaussianity from non-linear scales of the late-time Universe is an open problem. In this paper, we apply k-Nearest Neighbour cumulative distribution functions, kNN-CDFs, to the quijote-png simulations to explore the sensitivity of kNN-CDFs to primordial non-Gaussianity. An interesting result is that for halo samples with $M_\mathrm{ h}\langle 10^{14}$ M$_\odot$ $h^{-1}$, the kNN-CDFs respond to equilateral PNG in a manner distinct from the other parameters. This persists in the galaxy catalogues in redshift space and can be differentiated from the impact of galaxy modelling, at least within the halo occupation distribution (HOD) framework considered here. kNN-CDFs are related to counts-in-cells and, through mapping a subset of the kNN-CDF measurements into the count-in-cells picture, we show that our results can be modelled analytically. A caveat of the analysis is that we only consider the HOD framework, including assembly bias. It will be interesting to validate these results with other techniques for modelling the galaxy–halo connection, e.g. (hybrid) effective field theory or semi-analytical methods.
Dark Energy Survey Year 6 results: cell-based coadds and METADETECTION weak lensing shape catalogue
We present the metadetection weak lensing galaxy shape catalogue from the 6-yr Dark Energy Survey (DES Y6) imaging data. This data set is the final release from DES, spanning 4422 deg 2 of the southern sky. We describe how the catalogue was constructed, including the two new major processing steps, cell-based image coaddition, and shear measurements with metadetection. The DES Y6 M etadetection weak lensing shape catalogue consists of 151 922 791 galaxies detected over riz bands, with an effective number density of n eff = 8.22 galaxies per arcmin 2 and shape noise of σ e = 0.29. We carry out a suite of validation tests on the catalogue, including testing for point spread function (PSF) leakage, testing for the impact of PSF modelling errors, and testing the correlation of the shear measurements with galaxy, PSF, and survey properties. In addition to demonstrating that our catalogue is robust for weak lensing science, we use the DES Y6 image simulation suite to estimate the overall multiplicative shear bias of our shear measurement pipeline. We find no detectable multiplicative bias at the roughly half-per cent level, with m = (3.4 ± 6.1) x 10 –3 , at 3σ uncertainty. This is the first time both cell-based coaddition and Metadetection algorithms are applied to observational data, paving the way to the Stage-IV weak lensing surveys.
Dark Energy Survey Year 3 results: Measurement of the baryon acoustic oscillations with three-dimensional clustering
The three-dimensional correlation function offers an effective way to summarize the correlation of the large-scale structure even for imaging galaxy surveys. We have applied the projected three-dimensional correlation function, ξ p to measure the baryonic acoustic oscillations (BAO) scale on the first-three years Dark Energy Survey data. The sample consists of about 7 million galaxies in the redshift range 0.6< z p < 1.1 over a footprint of 4108 deg 2 . Our theory modeling includes the impact of realistic true redshift distributions beyond Gaussian photo-z approximation. ξ p is obtained by projecting the three-dimensional correlation to the transverse direction. To increase the signal-to-noise of the measurements, we have considered a Gaussian stacking window function in place of the commonly used top-hat. ξ p is sensitive to D M (z eff )/r s , the ratio between the comoving angular diameter distance and the sound horizon. Using the full sample, D M (z eff )/r s is constrained to be19.00 ± 0.67 (top-hat) and 19.15 ± 0.58 (Gaussian) at z eff = 0.835. The constraint is weaker than the angular correlation w constraint 18.84 ± 0.50), and we trace this to the fact that the BAO signals are heterogeneous across redshift. While ξ p responds to the heterogeneous signals by enlarging the error bar, w can still give a tight bound on D M s in this case. When a homogeneous BAO-signal subsample in the range 0.7 < z p <1.0 (z eff = 0.845) is considered, ξp yields 19.80 ± 0.67 (top-hat) and 19.84 ± 0.53 (Gaussian). The latter is mildly stronger than the w constraint (19.86 ± 0.55). We find that the ξ p results are more sensitive to photo-z because ξ p keeps the three-dimensional clustering information causing it to be more prone to photo-z noise. The Gaussian window gives more robust results than the top-hat as the former is designed to suppress the low signal modes. ξ p and the angular statistics such as w have their own pros and cons, and they serve an important crosscheck with each other.
Dynamic Strength of Iron at High Pressures and Strain Rates
Accurate modeling of meteorite impacts, and deformation of planetary cores require characterization of the flow strength and in-elasticity of iron in its different phases. In this letter, we investigate the flow strength of both the ambient α-phase and high-pressure ϵ-phase of iron at strain-rates of 1 × 10 5 s –1 and pressures up to 42 GPa using high pressure - pressure shear plate impact experiments. We report the strength of the ϵ-iron to be significantly higher than α-phase but consequently one order smaller than the previously reported dynamic strength at high pressures. As a result, the complete stress-strain response of the ϵ-phase is reported for the first time.
Automating the Solar Interconnection Technical Evaluation Process: PREconfiguring and Controlling Inverter SEt-Points (PRECISE)
Utilities are receiving an increasing number of interconnection requests for distributed solar photovoltaic (PV) systems from their customers, requiring a solution to quickly and intelligently perform technical assessments of these requests and to manage these new renewable energy assets. In this paper, we present a stand-alone interconnection tool, PREconfiguring and Controlling Inverter SEt-points (PRECISETM), and how it overcomes the challenges of accelerating technical evaluations and leverages smart inverter functionality to provide local voltage support on an as-needed basis. The paper details the purpose and architecture of PRECISE, including tools, data source requirements, the model-based interconnection evaluation process, and technical results. PRECISE is an integrated software tool that customizes PV inverter settings for utilities by modeling and assessing impacts on local voltages and the need for advanced inverter functions (AIFs) (e.g., volt-VAR and volt-watt). PRECISE maximizes the use of readily available utility data sets to make fast online assessments of incoming PV interconnection requests and makes accept/reject recommendations along with custom settings for AIFs.
Chemical signature characterization with hyperspectral imagery: novel deep learning model architectures and physically-motivated data augmentation techniques
The high spectral resolution afforded by Hyperspectral Imaging (HSI) sensors is poised to bring unprecedented advancements to signature characterization applications. Thus far, much of the research in the machine learning field devoted to HSI applications has focused on a few specific tasks like land-use land-cover classification. In land classification tasks, spatial information is very important, and model architectures are often designed to leverage spatial contexts. However, it is unclear how well these spatially-tuned models will translate to tasks where spectral information is critical, like the detection and characterization of chemicals. In this work, we compare spectral models (inputs are 1D spectra) and spatial-spectral models (inputs are 3D cubes) in the context of predicting chemical concentration maps. We find that spatial-spectral models perform the best, though we find a wide range in performance across the different architectures tested. Additionally, we find that model performance is impacted by the availability of training data, particularly in scenarios where the training data doesn't fully capture the true variance of real-world conditions. We find that data augmentation can help mitigate sparse coverage of observed parameter space (e.g., seasonal or geographic variability in ground cover), and present augmentation strategies that are tailored to hyperspectral data.