Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Aggregates”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Impacts of advanced diesel combustion operation and fuel formulation on soot nanostructure and reactivity

Advanced diesel combustion, accomplished via a single pulse fuel injection and high levels of exhaust gas recirculation (referred to as “PCCI”, partially-premixed charge compression ignition), is shown to be a path to reduce oxides of nitrogen and particulate matter simultaneously. This is well established in the literature. Less established is the extent to which such dilute combustion processes influence soot formation and affect soot that is emitted from diesel engines under such combustion modes. This work focuses on characterization of the nanostructure and oxidative reactivity of soot generated by a light-duty turbodiesel engine operating under a PCCI combustion mode, a dilute, low-temperature combustion process. Previous work on a type of PCCI combustion, referred to as high-efficiency clean combustion (HECC), showed soot samples having a fullerenic nanostructure, characterized by high levels of tortuosity of the fringe layers as seen in transmission electron micrograph images and as quantified using an image processing algorithm. Thermogravimetric analysis of the HECC mode soot samples showed that they displayed higher rates of oxidation than soot samples from a conventional diesel combustion mode. The present work returns to PCCI combustion, considering the timing of the main fuel injection, and the effects of operating on fuels rich in n-alkanes, particularly a fuel produced from a low temperature Fischer-Tropsch process (LTFT) and a renewable diesel fuel (RD) produced via hydrodeoxygenation of a plant oil. PCCI combustion conditions yield soot that shows higher reactivity compared to soot from the conventional combustion mode regardless of fuel type. LTFT and RD fuels produce soots with lower reactivity compared to ULSD. Soots produced from PCCI combustion have higher surface oxygen concentration and higher proportion of amorphous carbon. In addition, TEM images show that PCCI soots from all three fuels have smaller primary particle and particle aggregate sizes, and smaller graphene layers. These properties explain the higher reactivity of soot from PCCI combustion. The less reactive soots, which are produced from LTFT and RD fuel under conventional combustion, show internal burning during oxidation. However, soots with higher reactivity which are produced from late injection PCCI combustion and ULSD show shrinking core oxidation, most likely because of their overall amorphous structure.

42 ENGINEERING↗
DOE OSTI2021

Analysis of digitized herbarium records and community science observations provides a glimpse of downy mildew species diversity of North America, reveals potentially undescribed species, and documents the need for continued digitization and collecting

Downy mildew diseases caused by Peronosporaceae cause significant crop losses globally, with several emerging and resurgent threats in recent decades. Biodiversity data from digitized herbarium specimens provide an opportunity to develop a baseline census of species diversity, however, these resources may represent aggregations of nonrandom and opportunistic collecting efforts, which could lead to spurious results. Here, the MyCoPortal census of digitized herbarium records for downy mildew species collected from North America 1800 to present were analyzed. From 9838 unique records, 196 species were identified, reflecting ~28% of known species diversity. Temporal and geographic collecting biases were observed, with 90% of the collections made prior to 1960 and the efforts of six “super-collectors” accounting for 25% of the collections. The presence of 50–100 undescribed species in North America was inferred from the records. Together, these results highlight the need for continued downy mildew collections, taxonomic research and digitization efforts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Co-precipitation induces changes to iron and carbon chemistry and spatial distribution at the nanometer scale

Association of organic matter (OM) with mineral phases via co-precipitation is expected to be a widespread process in environments with high OM input and frequent mineral dissolution and re-precipitation. In contrast to surface area-limited adsorption processes, co-precipitation may allow for greater carbon (C) accumulation. However, the potential sub-micrometer scale structural and compositional differences that affect the bioavailability of co-precipitated C are largely unknown. In this work, we used a combination of high-resolution analytical electron microscopy and bulk spectroscopy to probe interactions between a mineral phase (ferrihydrite, nominally Fe 2 O 3 •0.5H 2 O) and organic soil-derived water-extractable OM (WEOM). In co-precipitated WEOM-Fe, nanometer-scale scanning transmission electron microscopy with electron energy loss spectroscopy (STEM-EELS) revealed increased Fe(II) and less Fe aggregation relative to adsorbed WEOM-Fe. Spatially distinct lower- and higher-energy C regions were detected in both adsorbed and co-precipitated WEOM-Fe. In co-precipitates, lower-energy aromatic and/or substituted aromatic C was spatially associated with reduced Fe(II), but higher-energy oxidized C was enriched at the oxidized Fe(III) interface. Therefore, we show that co-precipitation does not constitute a non-specific physical encapsulation of C that only affects Fe chemistry and spatial distribution, but may cause a bi-directional set of reactions that lead to spatial separation and transformation of both Fe and C forms. In particular, we propose that abiotic redox reactions between Fe and C via substituted aromatic groups (e.g., hydroquinones) play a role in creating distinct co-precipitate composition, with potential implications for its mineralization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation

Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.

58 GEOSCIENCES↗
DOE OSTI2022

On the provenance of GEMS, a quarter century post discovery

The provenance of GEMS (glass with embedded metal and sulfides) in cometary type interplanetary dust particles is investigated using analytical scanning transmission electron microscopy and secondary ion mass spectrometry. We review the current state of knowledge and closely examine the densities, elemental compositions and distributions, iron oxidation states, and isotopic compositions of a subset of GEMS in chondritic porous interplanetary dust. We find that GEMS are underdense with estimated densities that are 35–65% of compositionally equivalent crystalline aggregates. GEMS low densities result in a lower contribution to the bulk compositions of IDPs than has been assumed based on their volume fraction. We also find that element/Si ratios, assumed to be primary (indigenous), are instead perturbed by contamination and secondary alteration, including pulse heating during atmospheric entry. Fe in pyrrhotite inclusions was oxidized and Mg, S, Ca, and Fe were depleted relative to lithophile Al and Si, resulting in reduction in element/Si ratios. Because they trap outgassing elements, Fe-rich oxide rims that formed on the surface of GEMS are serendipitous “witness plates” to the changes in composition that accompany atmospheric entry. As a result of alteration, GEMS elemental compositions cannot reliably inform about their provenance. Except for highly anomalous oxygen isotope ratios measured in some large GEMS grains that indicate a contribution from circumstellar dust, oxygen isotope compositions are generally poor indicators of provenance. Prior work indicates that most GEMS fall close to the terrestrial oxygen isotope composition, which, however, does not exclude a presolar interstellar origin. Nitrogen isotopic compositions are more diagnostic. Elevated 15 N/ 14 N ratios indicate that GEMS accreted in conjunction with formation of organic matter by ion–molecule reactions in a cold (<50 K) presolar environment like the extreme outer nebula or interstellar medium. Considering all observations, we conclude that GEMS are most likely processed interstellar silicates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Composition and persistence of soil organic matter along eroding and depositional transects in buried vs. modern soil layers: A case of the Brady paleosol at Wauneta, Nebraska

Paleosols form when soils are buried through deposition by aeolian, colluvial, alluvial or other processes. Burial of former topsoil isolates soil organic matter (SOM) from surface conditions, allowing carbon to accumulate and potentially remain stable for millennia. In this study, SOM composition, distribution, and persistence were analyzed in the Brady Soil of Nebraska, USA to compare SOM spatial variability in modern and buried soils, as well as the impact of erosional exposure on SOM stability. The Brady Soil, formed as a surface soil during the Pleistocene-Holocene transition and now a paleosol buried up to 6 m deep (or more) by loess deposition during the Holocene, was sampled along burial (up to 5.8 m depth) and erosional (up to 1.8 m depth) transects to compare SOM dynamics in different geomorphic settings. Fourier Transform Infrared Spectroscopy (FTIR) and Fourier Transform ion cyclotron resonance mass spectrometry (FTICR-MS) were used to analyze SOM composition, while δ 13 C isotope analyses identified SOM sources and radiocarbon values were used to estimate turnover rates. Results confirmed a vegetation shift from C3 to C4 plants after Brady Soil formation, reflecting warming climatic conditions. Increasing SOM age and decreasing δ 13 C and δ 15 N values with depth indicated slowing of decomposition rate in buried soils. Higher pH in the Brady Soil suggested greater base cation content, supporting SOM stabilization through organo-mineral associations and aggregate formation. However, exposure of the Brady Soil due to surface erosion caused faster SOM turnover. This result suggested susceptibility of buried SOM to losses via decomposition upon erosional exposure, possibly accelerated by priming in response to modern SOM inputs. These findings highlight the potential loss of carbon stocks in buried soils under future climate change, as shifts in soil physicochemical properties may destabilize long-preserved SOM.

Environmental sciences↗
DOE OSTI2022

Induced seismicity red-light thresholds for enhanced geothermal prospects in the Netherlands

Prospects for geothermal energy in the Netherlands have renewed concerns around induced earthquakes. Risks from induced earthquakes are managed by traffic light protocols (TLPs), where the red-light is chosen as the stop-point before exceeding a tolerance to risk. Here, we simulate post-shut-in earthquake scenarios based on realistic information for the Netherlands. We focus on three risk metrics: aggregates like nuisance and damage impacts and also local personal risk (LPR) – a likelihood of building collapse fatality for an individual. Our results show that the severity of these risks varies spatially by orders of magnitude. Prior induced seismicity (e.g., the 2012 Huizinge event) provides a reference baseline to calibrate the Dutch earthquake risk tolerances. We find that these calibrated risk tolerances are similar to those observed in North America, suggesting an underlying sociological ‘license to operate.’ Furthermore, the use of calibrated risk tolerances results in nuisance concerns completely eclipsing the other two metrics. We compare our results to a hypothetical Groningen geothermal operation and find that our approach sets red-light thresholds approximately one magnitude unit below the ML 3.6 Huizinge event. Overall, our results provide a first-order recommendation for red-light thresholds and proactive management of Dutch enhanced geothermal induced seismicity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Real-time streamflow forecasting: AI vs. Hydrologic insights

In this paper, we propose a set of simple benchmarks for the evaluation of data-based models for real-time streamflow forecasting, such as those developed with sophisticated Artificial Intelligence (AI) algorithms. The benchmarks are also data-based and provide context to judge incremental improvements in the performance metrics from the more complicated approaches. The benchmarks include temporal and spatial persistence, persistence corrected for baseflow and streamflow, as well as river distance weighted runoff obtained from space-time distributed rainfall. In the development of the benchmarks, we use basic hydrologic insights such as flow aggregation by the river network, scale-dependence in basin response, streamflow partitioning into quick flow and baseflow, water travel time, and rainfall averaging by the basin width function. The study uses 140 streamflow gauges in Iowa that cover a range of basin scales between 7 and 37,000 km 2 . The data cover 17 years. This work demonstrates that the proposed benchmarks can provide good performance according to several commonly used metrics. For example, streamflow forecasting at half of the test locations across years achieves a Kling-Gupta Efficiency (KGE) score of 0.6 or higher at one-day ahead lead time, and 20% of cases reach the KGE of 0.8 or higher. The proposed benchmarks are easy to implement and should prove useful for developers of data-based as well as physics-based hydrologic models and real-time data assimilation techniques.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Robust Adaptive Control for Large-scale Inverter-based Resources with Partial and Complete Loss of Inverters

This article proposes an approach to address the current and the aggregated active power control challenge for large-scale inverter-based resources subjected to partially or completely loss of inverters and grid voltage variations. To address this problem, a distributed active power mechanism is proposed which generates desired inverters current. Then an adaptive mechanism distributes the voltage control input automatically between inverters in response to the partial or complete loss of inverters. An L2-gain-based controller is designed for each inverter to track the desired current and rejects the grid voltage disturbance. Here, simulation results show a significant robust tracking for collective active power and current.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Investigating the Determinants of Household Capabilities Burden During Power Outages: The Case of Winter Storm Uri

Existing research primarily uses census data to identify the vulnerability of communities to hazards. These vulnerability indices provide aggregated data and are not hazard-specific nor well-validated with post-event data. In contrast, our study uses household survey data (n=1065) to understand which Texan households suffered the greatest loss of their capabilities due to power outages and other utility service disruptions during Winter Storm Uri. Inspired by the Capabilities Approach, our measures of burden include the number of household capability types disrupted during the outages (e.g., cooking, heating, refrigeration), the severity of impact for each disrupted capability, and the additional time and financial costs of coping with these disruptions. We perform a clustering analysis, and find two distinct groups in our data, consisting of ‘lesser burden' and ‘heavier burden' households. Results indicate that the households experiencing the heaviest capabilities burden were most likely to experience longer power outages and the loss of water services. They were also more likely to have a Hispanic-Latino household member, lack access to a generator, live in a rented home, have larger households with more young children, fewer adults over 65, lower household incomes, been impacted by the COVID-19 pandemic, and more family characteristics that made life harder. We also fit a logistic regression model to assess the role of outage, household, and community characteristics in predicting differences in capabilities burden. Our results offer insights into enumerating the consequences of utility service disruptions on households, which can inform more targeted and equitable resilience strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Power system inertia estimation: Review of methods and the impacts of converter-interfaced generations

Understanding and quantifying the inertia of power systems with the integration of converter-interfaced generation (CIG) plays an essential role in the safe transition to a future low-inertia scenario. Here we provide a comprehensive summary of inertia definitions for both synchronous generators and CIGs as well as their corresponding estimation methods. In particular, the estimation methods are categorized as model-based and measurement-based approaches considering both small and large disturbances. The advantages and disadvantages of different methods are carefully discussed. This paper also offers for the first time a framework to quantify the virtual inertia of CIGs at the component and aggregation levels, an open problem in the literature. Finally, future directions for inertia estimation are identified and discussed. This significantly benefits the design of appropriate control and protection schemes in achieving a more reliable, secure, and resilient power system.

42 ENGINEERING↗
DOE OSTI2022

Synergistic effect of microstructure and defects on the initiation of fatigue cracks in additively manufactured Inconel 718

Fatigue cracks in additively manufactured (AMed) Inconel 718 (IN-718) in machined surface condition often initiate from persistent slip bands (PSBs) unlike other popular AM alloys such as 17-4 PH stainless steel or Ti-6Al-4V, where fatigue crack initiation is exclusively from volumetric defects; therefore, a competition between PSB- vs. defect- mediated crack initiation clearly exists. To shed light on the factors governing the competition, this study investigates the characteristics of cyclic strain localization, PSB formation, and crack initiation via crystal plasticity (CP) modeling of cyclic loading on polycrystalline aggregates which are then validated by experiments. In this work, a physics-based, free slip distance (FSD) dependent slip strength evolution law is proposed, which is shown to enable the CP model to simulate the heterogeneous strain distribution in IN-718. Implementing a crack initiation criterion based on strain contrast, the locations and lives for crack initiation can also be calculated. It is shown that both FSD and resolve shear stress influence the strain localization and crack initiation behaviors. The distribution of a localization parameter calculated based on the multiplication of FSD and Schmid factor within a grain is found to correlate well with the locations of PSBs. The maximum values of the localization parameters within a microstructure are shown to correlate well with the experimentally obtained crack initiation lives. The presence of volumetric defects in IN-718 generally do not impact strain localization behavior unless their size is large compared to the grain size.

42 ENGINEERING↗
DOE OSTI2022

A multi-criteria CCUS screening evaluation of the Gulf of Mexico, USA

Continued research into reservoir characterization along with offshore carbon dioxide (CO 2 ) transportation and infrastructure assets is needed to facilitate development of safe and successful carbon capture, utilization, and storage (CCUS) projects. This paper outlines a multi-criteria evaluation methodology that incorporates disparate sets of quantitative, spatially variable data into a decision-making framework for screening the Gulf of Mexico (GOM) outer continental shelf (OCS) for potentially viable CO 2 storage and enhanced oil recovery (EOR) sites. Criteria categories include favorable geologic characteristics, logistics, and potential risks. Data compiled for 14 criteria from several publicly available geographic information system (GIS) layers was aggregated over 2559 spatially balanced points across the study area using the National Energy Technology Laboratory (NETL)-developed Cumulative Spatial Impact Layers™ (CSIL) GIS tool. Criteria are weighted by qualitative expert opinion relative to their perceived importance to given scenarios— the output of combined criteria values and weights enables regional CO 2 storage suitability differentiation. The methodology considers both technical and non-technical factors impacting CCUS decision-making. The flexible methodology enables a systematic approach to regional ranking at high spatial resolution over a large study domain. Additionally, the framework enables high-grading of priority sites that warrant further characterization and follow-on analysis. Areas along the Louisiana coast and Mississippi River Delta consistently rank high for all scenarios largely a result of the favorable geology with the potential for stacked storage, as well as the density of existing pipelines and platforms, and proximity to several onshore CO 2 sources. High-graded regions for the CO 2 EOR-related scenarios are typically located further offshore towards the middle and edge of the OCS compared to higher priority regions for the geologic storage scenarios which fall closer to the Louisiana coastline.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An experimental investigation on the CO 2 storage capacity of the composite confining system

Assuring secure containment of stored CO 2 is of paramount importance—for climate change mitigation, for permitting, and for reassuring the public. Regional seals, such as those sealing petroleum accumulations, have proven to be effective for securing CO 2 . However, the goal of permanent sequestration can also be satisfied with “composite confining systems” consisting of multiple, possibly discontinuous flow barriers that, in aggregate, create a system with very high permeability anisotropy and effectively retard the vertical migration of CO 2 . This study focuses on investigating the barrier characteristics necessary for effective containment of CO 2 , using both physical flow experiments and modified invasion percolation simulations. The simulations are calibrated to the physical experiments and are used to further extend the analysis. Results show that for a composite confining system, a) even barriers with low capillary entry pressure contrast to the underlying flow unit can divert rising CO 2 ; b) curved or anticlinal barrier topography can enhance CO 2 trapping; c) fining-upward gradations make little difference to barrier effectiveness or CO 2 retention; d) longer barriers retain more CO 2 regardless of barrier topography. Finally, field-scale simulations have demonstrated the importance of barrier length for increasing CO 2 storage capacity. Furthermore, the results presented can be used to inform the development of new screening criteria for characterization and effectiveness of composite confining systems.

58 GEOSCIENCES↗
DOE OSTI2026

Carbon dioxide pipeline network transportation cost model: evaluating economic and geographic factors for efficient carbon capture, storage, and utilization

This study presents a comprehensive pipeline network modeling framework to estimate the CO 2 delivery cost for CO 2 utilization and geologic CO 2 storage across the United States. We developed a Python-based CO 2 pipeline transportation cost model leveraging Argonne National Laboratory’s pipeline engineering expertise and detailed natural gas transmission pipeline cost data across U.S. regions. Using existing road corridors as practical routing guides, the model designs pipeline networks that aggregate CO 2 from one or multiple sources and deliver it to selected destinations. It then minimizes the total transportation cost by optimizing pipeline diameters and incorporating booster pumps. A key contribution is the incorporation of up-to-date, region-specific cost factors with itemized components for materials, labor, miscellaneous construction expenses, and right-of-way acquisition. Results emphasize that regional variation and economies of scale associated with CO 2 pipeline costs are significant and should be explicitly accounted for in screening and planning studies. By combining realistic routing constraints with regionalized cost inputs, the model provides transparent design methodology and location-specific insights into source–destination delivery costs, including the effects of routing complexity along existing road networks. We demonstrate the model with two illustrative case studies – one for CO 2 storage and one for CO 2 utilization – in which the model designs pipeline networks spanning hundreds of miles across the states, collecting CO 2 from multiple sources and delivering it to designated endpoints while minimizing levelized cost of delivery via diameter and compression optimization. The model offers a practical, scalable approach for alternative design option screening and early-stage CO 2 transportation planning.

CCS↗
DOE OSTI2026

System-of-systems optimization of hydrogen infrastructure for heavy-duty freight corridors: The interstate 10 case study

Medium and heavy-duty freight transportation requires hydrogen energy infrastructure that is cost-effective, operationally reliable, spatially coherent, and resilient to demand variability along major corridors. This paper presents an integrated hydrogen corridor planning framework using Oak Ridge National Laboratory's OR-AGENT that couples freight-driven, route-resolved hydrogen demand modeling with optimized station siting, sizing, and station-level techno-economic analysis. The framework is demonstrated for the Interstate 10 freight corridor and the Houston-to-Los-Angeles region. Hydrogen demand is derived from high-resolution origin–destination freight data, duty-cycle characterization, and physics-based energy consumption modeling. Candidate refueling sites are selected from existing heavy-duty diesel fueling locations and optimized subject to onboard storage and station capacity constraints. Resulting station throughputs are evaluated using established techno-economic models for electrolytic hydrogen production and dispensing infrastructure. Results show that a regional, portfolio-level aggregation, average dispensed electrolytic hydrogen cost of $6.87–$7.26/kg is currently feasible, and is strongly influenced by demand density and utilization.

Sujan, Vivek [ORNL] (ORCID:0000000269882342)↗
DOE OSTI2020

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗