Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability analysis”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Empirical Indicators of Transmission Value in the Southeast United States

Concurrent differences in energy price between different parts of the electric grid are a key indicator of the value of additional transmission. In areas without a wholesale electricity market, such as the Southeast, an alternative indicator to price is the Federal Energy Regulatory Commission’s (FERC) system lambda data. This economic metric represents the minimized marginal production costs of thermal generators, including fuel and other variable operation and maintenance expenses. Balancing Authorities report a single system lambda for their entire balancing area. Most Southeastern lambdas exhibit sufficient price variation to support a transmission valuation analysis, although incomplete accounting of congestion costs or scarcity rents during peak load hours may underestimate the true value of transmission capacity. With transmission value defined as the annual average hourly absolute price difference between two regions and FERC’s system lambda data used as a price proxy, we find the following results in the Southeast region during 2012-2023 (reported in $\$2024$/MWh): Intra‐regional findings: Annual averages historically span $\$2$–$\$28$/MWh and average $\$12$/MWh in SERTP and span $\$4$–$\$19$/MWh and average $\$9$/MWh in FRCC, disregarding transmission value driven by anomalous data. The ranges of transmission value reported here are large, spanning an order of magnitude in some cases. Much of this variation is driven by year-to-year changes, with 2022 having a particularly high intra-regional transmission value due to elevated natural gas prices. Inter‐regional corridors: Annual average transmission values across three broader regions range from $\$6$ to $\$28$/MWh with a long-term average of $\$11$/MWh. Much of the transmission value is concentrated in a small portion of hours. Across all regions, severe weather—particularly polar vortex events in January 2018, February 2021, and December 2022—drives the largest price spreads. Seasonal patterns also emerge, with summer afternoons and fall mornings contributing consistently to transmission value, as for example between MISO and SOCO in 2023.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Environmental DNA Detection in Marine Macrophyte Ecosystems as a Potential Blue Carbon Source in Sediments

“Blue carbon” refers to the carbon sequestered by the world’s oceanic and coastal ecosystems, particularly through coastal vegetation such as mangroves, salt marshes, seagrasses, and marine macroalgae. These ecosystems play a crucial role in the global carbon cycle by serving as significant carbon sinks, absorbing carbon dioxide from the atmosphere and storing it in biomass and sediments over long periods. This study explores the use of environmental DNA (eDNA) to detect marine macrophytes and microalgae assemblages contributing to blue carbon in sediments across various coastal ecosystems. The research addresses the challenges of traditional monitoring methods by utilizing high-throughput sequencing of the 18S-V9 region amplified using eDNA from sediment samples collected at eight locations in the United States and South Korea. The results reveal a diverse array of taxa, underscoring the variability in community composition across different conditions. Notably, sites with seagrass beds and Ulva blooms showed distinct patterns in microalgal community structure. This study underscores the potential of eDNA analysis in providing comprehensive insights into the biodiversity of marine macrophyte ecosystems, thus informing conservation efforts and enhancing the understanding of marine ecological dynamics.

Xing, Qikun (ORCID:0000000206179027)↗
DOE OSTI2023

Fermi Constraints on the Ejecta Speed and Prompt Emission Region of the Distant GRB 220101A

At redshift z = 4.618, GRB 220101A is the most distant gamma-ray burst (GRB) detected by Fermi/LAT to date. It is also a very energetic event, with an equivalent isotropic energy of 3.6 × 10 54 erg. We jointly analyzed the Fermi/GBM and LAT observations of GRB 220101A with two independent approaches and found a significant spectral break at sub-100 MeV energies during the prompt emission. The fast variability of the emission suggests that this spectral attenuation is caused by internal opacity to pair creation. Regardless of the nature of the emission processes assumed in the spectral analysis, we infer a moderate value for the jet Lorentz factor, Γ ~ 110, and find that all of the high-energy emission was produced above and near the photosphere, at a distance of ~10 14 cm from the central engine. We compare these results with the four other LAT-detected GRBs with similar properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

BAWLD-CH 4 : a comprehensive dataset of methane fluxes from boreal and arctic ecosystems

Methane (CH 4 ) emissions from the boreal and arctic region are globally significant and highly sensitive to climate change. There is currently a wide range in estimates of high-latitude annual CH 4 fluxes, where estimates based on land cover inventories and empirical CH 4 flux data or process models (bottom-up approaches) generally are greater than atmospheric inversions (top-down approaches). A limitation of bottom-up approaches has been the lack of harmonization between inventories of site-level CH 4 flux data and the land cover classes present in high-latitude spatial datasets. Here we present a comprehensive dataset of small-scale, surface CH 4 flux data from 540 terrestrial sites (wetland and non-wetland) and 1247 aquatic sites (lakes and ponds), compiled from 189 studies. The Boreal–Arctic Wetland and Lake Methane Dataset (BAWLD-CH 4 ) was constructed in parallel with a compatible land cover dataset, sharing the same land cover classes to enable refined bottom-up assessments. BAWLD-CH 4 includes information on site-level CH 4 fluxes but also on study design (measurement method, timing, and frequency) and site characteristics (vegetation, climate, hydrology, soil, and sediment types, permafrost conditions, lake size and depth, and our determination of land cover class). The different land cover classes had distinct CH 4 fluxes, resulting from definitions that were either based on or co-varied with key environmental controls. Fluxes of CH 4 from terrestrial ecosystems were primarily influenced by water table position, soil temperature, and vegetation composition, while CH 4 fluxes from aquatic ecosystems were primarily influenced by water temperature, lake size, and lake genesis. Models could explain more of the between-site variability in CH 4 fluxes for terrestrial than aquatic ecosystems, likely due to both less precise assessments of lake CH 4 fluxes and fewer consistently reported lake site characteristics. Analysis of BAWLD-CH 4 identified both land cover classes and regions within the boreal and arctic domain, where future studies should be focused, alongside methodological approaches. Overall, BAWLD-CH 4 provides a comprehensive dataset of CH 4 emissions from high-latitude ecosystems that are useful for identifying research opportunities, for comparison against new field data, and model parameterization or validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Diffusion kinetics of 3 He in pyroxene and plagioclase and applications to cosmogenic exposure dating and paleothermometry in mafic rocks

In this study, we investigate the diffusivity of cosmogenic 3 He in a variety of plagioclase and pyroxene compositions, and its application to paleothermometry and exposure dating in these minerals, through stepwise degassing experiments. While cosmogenic 3 He has been utilized for exposure dating in pyroxene for decades due to its retentivity, plagioclase, often found along with pyroxene in mafic rocks, is generally less retentive of cosmogenic noble gas. However, the diffusivity of 3 He in either plagioclase or pyroxene has not yet been measured quantitatively. A challenge in measuring diffusion kinetics by step-degassing experiments in poorly retentive minerals is the fact that significant amounts of He can be lost prior to the experiment. To address this issue, we apply a forward “multiple diffusion domain” (MDD) inversion model that includes model predictions of initial gas loss during irradiation and storage of the samples to account for this observation and add constraints to the diffusion parameters. We find that 3 He diffusivity in plagioclase appears to be highly variable. This variability can be explained by the MDD inversion models' inability to constrain the diffusion parameters when significant gas has been lost during irradiation and/or prolonged storage prior to experiment analysis, resulting in an overestimation of 3 He retentivity. Plagioclase samples that were kept frozen after irradiation, to limit this initial gas loss, yielded the most reliable estimate of diffusion kinetics. We find that 3 He in plagioclase is diffusively lost at Earth's surface temperatures on a timescale of 100 years and is therefore unsuitable for surface temperature paleothermometry. On the contrary, we find cosmogenic 3 He in pyroxene to be retentive at Earth's surface temperatures on a 1-million-year timescale.

Bergelin, Marie [Berkeley Geochronology Center (BG↗
DOE OSTI2025

Raptor

Raptor is an efficient Python-based tool for predicting the formation and morphology of stochastic lack of fusion defects in metal AM processes. A major obstacle for the qualification and certification of additively manufactured parts in critical applications continues to be performance variability caused in part by porosity-related defects. High-fidelity process models that could predict these defect features are currently too computationally expensive for component-level analysis. To address this, Raptor employs a high-performance geometric method to model the dynamic melt pool rather than relying on computationally intensive thermal fluid dynamics. This allows Raptor to rapidly identify regions of unmelted material that correspond to lack of fusion pores. The efficiency of this approach significantly reduces the time and resources needed for generating 3D defect predictions, which enables users to conduct large-scale parameter studies and evaluate how process variations affect part quality. The framework offers operational flexibility; users can execute simulations through a simple command line interface or integrate core functions as a library within larger computational workflows. Simulation outputs include 3D porosity maps for visualization and tools for quantitative morphological analysis. These results are suitable for direct comparison with experimental characterization data from methods such as X-ray computed tomography and can be used for statistical process optimization.

Subraveti, Vamsi [Vanderbilt Univ., Nashville, TN ↗
DOE OSTI2020

Fine Structures in the Main Sequence Revealed by Gaia Data Release 2

The Hertzsprung–Russell diagram (HRD) is fully examined using the Fourier analysis. This work shows more stars are above the gap in the lower main sequence than below it, and this implies that stars spend more time above the gap while they undergo variability associated with the {sup 3}He instability. The enhanced HRD also shows the width of the gap is not linear and depends on the G{sub BP}−G{sub RP} color up until G{sub BP}−G{sub RP}=2.7. Beyond this color limit, the gap is hardly seen. Besides, a new low density region is revealed for the first time centered at M {sub G} ≈ 10.7 and G {sub BP} − G {sub RP} ≈ 2.8, which is below the lower right corner of the gap. This work also shows that the main sequence appears to have fine stripes where stellar densities are relatively low or high compared to their adjacent regions on the main sequence. These stripes can be seen throughout the main sequence of stars redder than G {sub BP} − G {sub RP} = 0.8 and are not limited to any specific color or spectral type. Slopes of these features are different from the main sequence, but are pretty consistent throughout the main sequence, with a few exceptions. We are perplexed by these new features, but the complexities of stellar atmospheric features and opacities of dwarfs may have caused these patterns.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Characterizing the uppermost 100 m structure of the San Jacinto fault zone southeast of Anza, California, through joint analysis of geological, topographic, seismic and resistivity data

We present results from complementary geological, topographic, seismic and electrical resistivity surveys at the Sagebrush Flat (SGB) site along the Clark fault (CF) strand of the San Jacinto fault zone trifurcation area southeast of Anza, California. Joint interpretation of these data sets, each with unique spatiotemporal sensitivities, allow us to better characterize the shallow (<100 m) fault zone at this structurally complex site. Geological mapping at the surface shows the CF has three main subparallel strands within a <100 m zone with varying degrees of rock damage. These strands intersect units of banded gneiss and tonalite, and various sedimentary units. Near the surface, the weathered but more intact tonalite and gneiss to the southwest have relatively high VP. The low-lying flat sedimentary basins around the two southwestern-most CF strands and elevated damaged gneiss to the northeast have lowest VP <500 m s–1. The high relief of the northeast gneiss unit may in part be explained by its extensive damage and inferred increased relative rock uplift. Resistivity imaging shows the unconsolidated dry basin sediments (maximum >1300 Ohm.m) contrasted against the compacted fine-grained (potentially wet) materials within the CF core and the Bautista Formation (minimum <40 Ohm.m), which is slightly elevated above the flat basins. The inverse relationship between VP (increases) and resistivity (decreases) in the uppermost ~15 m can be characterized as log–log linear with slopes of –2.6 to –4. At depths >30 m, the velocity heterogeneity near the surface merges into larger-scale structures that are generally slower on the northeast side of the CF core compared to the southwest side (as much as ~40 per cent reduction in average VP). A previous study revealed a 20–37 per cent variability in peak ground velocities across the SGB site from local earthquakes. The upper end of that range is associated with the near-surface unconsolidated sedimentary basins and northeast damaged gneiss unit. Preliminary analysis of time-dependent topography mostly shows effects of changing vegetation and anthropogenic activity.

Geochemistry & Geophysics↗
DOE OSTI2023

Field Validation of Innovative Air-Source Heat Pumps for Cold-Climate Heating Applications

Air-source heat pumps (ASHPs) have historically found application in milder climates in the United States (US). Their application in cold climates has been hindered by reduced performance as outdoor temperatures fall below freezing and the need for backup or auxiliary electric resistance heaters to meet peak heating loads. Recent advances in cold-climate air-source heat pumps (ccASHPs) with features such as variable-speed compressors, multistage systems, and highly optimized thermal design have improved performance by increasing the coefficient of performance (COP), low-ambient-temperature capacity, and heating seasonal performance factor (HSPF) dramatically. Despite their benefits, ccASHPs are still not widely prescribed for cold climate conditions. This is due in part to a lack of verified, demonstrated performance and analysis in cold climates. Ultimately, high-efficiency ccASHPs play a key role in emerging energy/grid renovation efforts in addition to their potential energy savings and avoided carbon emissions, but their actual performance in various cold climate field settings must be better understood. The primary objective of this project was to measure the in-field performance of variable-capacity air-source heat pumps in cold climates with the goal of enabling the development of field-based performance maps. Specifically, the study looked at how the heat pumps operated in the field, the frequency of cycling, the frequency of defrost events, and the time spent in each mode of operation. The results are intended to be used by DOE to inform research and development of energy-efficient equipment and to develop guidelines for optimizing primary energy savings when using air-source heat pumps in heating-dominated regions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Multiscale and multiphysics FEA simulation and materials optimization for laser ultrasound transducers

In this study, the relationship between the nanocomposite design and the laser ultrasound transducer (LUT) characteristics was investigated through simulations in multiple scale levels for material behavior, device response, and acoustic wave propagation in media. First, the effects of the nanoparticle size and concentration on the effective properties of composites were quantitatively investigated with the finite element analysis (FEA) method. Second, the effective properties of the nanocomposite were assigned to the layer, which is modeled as a homogeneous material, in the FEA for the LUT simulating the energy conversion from the incident laser to the acoustic wave. Finally, the ultrasound propagation in the water was calculated by a theoretical wave propagation model. The FEA-based prediction was compared with the experimental data in the literature and a theoretical analysis for LUT based on Thermal-Acoustic coupling. As a result, the ultrasound waves on the transducer surface and at a distance in the water could be predicted. Based on the hierarchically integrated prediction procedure, the optimal conditions of the photoacoustic nanocomposites were investigated through the parametric study with the particle size and concentration as variables. The results guide the material designs optimized for different device characteristics, such as high pressure and broad bandwidth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Generic and ML Workloads in an HPC Datacenter: Node Energy, Job Failures, and Node-Job Analysis

HPC datacenters offer a backbone to the modern digital society. Increasingly, they run Machine Learning (ML) jobs next to generic, compute-intensive workloads, supporting science, business, and other decision-making processes. However, understanding how ML jobs impact the operation of HPC datacenters, relative to generic jobs, remains desirable but understudied. In this work, we leverage long-term operational data, collected from a national-scale production HPC datacenter, and statistically compare how ML and generic jobs can impact the performance, failures, resource utilization, and energy consumption of HPC datacenters. Our study provides key insights, e.g., ML-related power usage causes GPU nodes to run into temperature limitations, median/mean runtime and failure rates are higher for ML jobs than for generic jobs, both ML and generic jobs exhibit highly variable arrival processes and resource demands, significant amounts of energy are spent on unsuccessfully terminating jobs, and concurrent jobs tend to terminate in the same state. We open-source our cleaned-up data traces on Zenodo (https://doi. org/10.5281/zenodo.13685426), and provide our analysis toolkit as software hosted on GitHub (https://github.com/atlarge-research/2024-icpads-hpc-workload-characterization). This study offers multiple benefits for data center administrators, who can improve operational efficiency, and for researchers, who can further improve system designs, scheduling techniques, etc.

crossanalysis↗
DOE OSTI2023

Variability in Diurnal and Seasonal Ambient Conditions on Geothermal Plant Performance and Cost: Preprint

Geothermal plant performance is bounded by the second law efficiency, which accounts for the quantity of exergy that can be converted into useful work. This in turn is dependent on the geothermal resource temperature and the temperature of the heat sink (i.e., the ambient temperature). In this study we show that ambient temperature variability on a diurnal and seasonal basis can affect performance and cost estimations for geothermal plants. We have utilized the updated System Advisor Model (SAM) to assess nine geothermal sites with existing resource capacities across three climate zones. Our analysis shows that both evaporatively-cooled flash and air-cooled binary cycle plants are affected by temperature, with a slightly higher effect in EGS binary sites. By assuming an ambient (wet bulb) temperature baseline of 15.6 degrees C (60 degrees F) and comparing baseline results to those from site-specific data we observe up to 15% underestimation of plant performance and up to 20% overestimation of cost.

ambient temperature↗
DOE OSTI2020

Sensitivity of meteorological-forcing resolution on hydrologic variables

Projecting the spatiotemporal changes in water resources under a no-analog future climate requires physically based integrated hydrologic models which simulate the transfer of water and energy across the earth's surface. These models show promise in the context of unprecedented climate extremes given their reliance on the underlying physics of the system as opposed to empirical relationships. However, these techniques are plagued by several sources of uncertainty, including the inaccuracy of input datasets such as meteorological forcing. These datasets, usually derived from climate models or satellite-based products, are typically only resolved on the order of tens to hundreds of kilometers, while hydrologic variables of interest (e.g., discharge and groundwater levels) require a resolution at much smaller scales. In this work, a high-resolution hydrologic model is forced with various resolutions of meteorological forcing (0.5 to 40.5 km) generated by a dynamical downscaling analysis from the regional climate model Weather Research and Forecasting (WRF). The Cosumnes watershed, which spans the Sierra Nevada and Central Valley interface of California (USA), exhibits semi-natural flow conditions due to its rare undammed river basin and is used here as a test bed to illustrate potential impacts of various resolutions of meteorological forcing on snow accumulation and snowmelt, surface runoff, infiltration, evapotranspiration, and groundwater levels. Results show that the errors in spatial distribution patterns impact land surface processes and can be delayed in time. Localized biases in groundwater levels can be as large as 5–10 m and 3 m in surface water. Most hydrologic variables reveal that biases are seasonally and spatially dependent, which can have serious implications for model calibration and ultimately water management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simultaneous isotopic analysis of fission product Sr, Mo, and Ru in spent nuclear fuel particles by resonance ionization mass spectrometry

Abstract Fission product Sr, Mo, and Ru isotopes in six 10-μm particles of spent fuel from a pressurized water reactor were analyzed by resonance ionization mass spectrometry (RIMS) and evaluated for utility in nuclear material characterization. Previous measurements on these same samples showed widely varying U, Pu, and Am isotopic compositions owing to the samples’ differing irradiation environments within the reactor. This is also seen in Mo and Ru isotopes, which have the added complication of exsolution from the UO 2 fuel matrix. This variability is a hindrance to interpreting data from a collection of particles with incomplete provenance since it is not always possible to assign particles to the same batch of fuel based on isotopic analyses alone. In contrast, the measured 90 Sr/ 88 Sr ratios were indistinguishable across all samples. Strontium isotopic analysis can therefore be used to connect samples with otherwise disparate isotopic compositions, allowing them to be grouped appropriately for interpretation. Strontium isotopic analysis also provides a robust chronometer for determining the time since fuel irradiation. Because of the very high sensitivity of RIMS, only a small fraction of material in each of the 10 μm samples was consumed, leaving the vast majority still available for other analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Scaling Turbulent Combustion Fields in Explosions

We considered the topic of explosions from spherical high-explosive (HE) charges. We studied how the turbulent combustion fields scale. On the basis of theories of dimensional analysis by Bridgman and similarity theories of Sedov and Barenblatt, we found that all fields scaled with the explosion length scale r0. This included the blast wave, the mean and root mean squared (RMS) profiles of thermodynamic variables, combustion variables, velocities, vorticity, and turbulent Reynolds stresses. This was a consequence of the formulation of the problem and our numerical method, which both satisfied the similarity conditions of Sedov. We performed numerical simulations of 1 g charges and 1 kg charges; the solutions were identical (within roundoff error) when plotted in scaled variables. We also explored scaling laws related to three-phase pyrotechnic explosions. We show that although the scaling formally broke down, the fireball still essentially scaled with the explosion length scale r0. However, the discrete Lagrange particles (DLP) (phase 2) and the heterogeneous continuum model (HCM) of the DLP wakes (phase 3) did not scale with r0, and mean and RMS profiles could differ by a factor of 10 in some regions. This was because the DLP particles and wakes introduced an additional scale that broke the similarity conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Behavioral and neuroinflammatory changes caused by glyphosate: Base herbicide in mice offspring

Abstract Background Glyphosate is a pesticide considered of low toxicity, but scientific evidences show it can be harmful to health. This study aimed to evaluate the toxicity in mice offspring exposed to glyphosate‐based herbicide (GBH) during the intrauterine period. Methods Female matrices received glyphosate 0.3 mg/kg daily per oral throughout the gestational period, which was variable between 18 and 22 days. From the 25th until the 28th days post‐birth, mice offspring were subjected to behavioral tests, and the prefrontal cortex was processed for immunohistochemical analysis. Results Two significant behavioral changes were observed: anxiety in the GLIF0.3 group, increase in the behavior burying marbles in the marble‐burying test and hyperactivity, expressed by the significant increase of the crossing number in the open field test. The increased microglia, TNF‐alpha, and astrocyte expression were also observed in the prefrontal cortex of offspring treated with GLIF0.3. Conclusion Exposure to GBH during mice intrauterine development induces hyperactive and anxious behavior, evidencing neuroinflammation.

de Castro Vieira Carneiro, Carla Larissa↗
DOE OSTI2022

The value of in-reservoir energy storage for flexible dispatch of geothermal power

Geothermal systems making use of advanced drilling and well stimulation techniques have the potential to provide tens to hundreds of gigawatts of clean electricity generation in the United States by 2050. With near-zero variable costs, geothermal plants have traditionally been envisioned as providing “baseload” power, generating at their maximum rated output at all times. However, as variable renewable energy sources (VREs) see greater deployment in energy markets, baseload power is becoming increasingly less competitive relative to flexible, dispatchable generation and energy storage. Herein we conduct an analysis of the potential for future geothermal plants to provide both of these services, taking advantage of the natural properties of confined, engineered geothermal reservoirs to store energy in the form of accumulated, pressurized geofluid and provide flexible load-following generation. We develop a linear optimization model based on multi-physics reservoir simulations that captures the transient pressure and flow behaviors within a confined, engineered geothermal reservoir. We then optimize the investment decisions and hourly operations of a power plant exploiting such a reservoir against a set of historical and modeled future electricity price series. Further, we find that operational flexibility and in-reservoir energy storage can significantly enhance the value of geothermal plants in markets with high VRE penetration, with energy value improvements of up to 60% relative to conventional baseload plants operating under identical conditions. Across a range of realistic subsurface and operational conditions, our modeling demonstrates that confined, engineered geothermal reservoirs can provide large and effectively free energy storage capacity, with round-trip storage efficiencies comparable to those of leading grid-scale energy storage technologies. Optimized operational strategies indicate that flexible geothermal plants can provide both short- and long-duration energy storage, prioritizing output during periods of high electricity prices. Sensitivity analysis assesses the variation in outcomes across a range of subsurface conditions and cost scenarios.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

A molecular dynamics survey study of impact-induced amorphization in yttria nanoparticles

The impact of a single yttria (Y 2 O 3 ) nanoparticle onto a yttria substrate was studied as a function of particle velocity (300–1200 m/s) and diameter (5–12 nm) for three particle orientations and two substrate orientations using molecular dynamics simulations. Impacts resulted in minimal dislocation motion but extensive disordering of the crystal lattice without approaching the melting temperature. The resulting amorphous microstructure was both precisely coincident with plastic strain and retained until the end of the simulation, suggesting a deformation mechanism of amorphization-facilitated viscous flow. Further analysis revealed that the degree of amorphization depended primarily on the impact velocity and secondarily on the orientations of the particle and substrate, with significant second-order effects from variable interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗