Results for “surface energetics” — Page 42

NASA NTRS2008

Lunar Dust Charging by Secondary Electron Emission and its Complex Role in the Lunar Environment

The lunar surface is covered with a thick layer of micron/sub-micron size dust grains formed by billions of years of meteoritic impact. With virtually no atmosphere and exposed to the solar wind plasma and solar electromagnetic radiation, the lunar surface and the dust grains are electrostatically charged. The dominant charging processes include: photoelectric emissions (UV, X-rays), impact of solar wind electrons and ions, and secondary electron emissions (SEE) induced by energetic solar wind electrons. During the Apollo missions, the astronauts found the lunar dust to be extraordinarily high in its adhesive characteristics, sticking to the suits and the mechanical equipment. Electrostatically charged lunar dust is believed to be transported over long distances by the induced electric fields, as indicated by the observed dust streamers and the horizon glow [e.g., 1-3]. The hazardous effects of dust in the lunar environment are recognized to be one of the major issues that must be addressed in planning the forthcoming missions for robotic and human exploration of the Moon. Theoretical studies are being performed along with the development of analytical models and a variety of experimental investigations, to better understand the lunar dust phenomena. [e.g., 4-6]. The lunar dust is believed to be charged negatively on the lunar night-side by interaction With solar wind electrons. However, rigorous theoretical expressions for calculation of SEE yields and the sticking efficiencies of individual micron size dust grains are not yet available, and the information has to be obtained by experiment. On theoretical considerations, however, it is well recognized that SEE yields, similar to the photoelectric yields for small-size grains, would be totally different from the corresponding bulk values [e.g., 7-9]. Some theoretical models for charging of individual small spherical particles have been developed [e.g., 10], and some limited measurements on individual metallic dust grains at keV electron energies have been made [e.g., i 1]. In this paper, we present the first measurements of the secondary electron emission yields of individual micron/sub-micron size dust grains selected from sample returns of Apollo 11 and Apollo 17 missions.

Abbas, M. M.↗

NASA NTRSDate not supplied

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa↗

NASA NTRSDate not supplied

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa↗

DOE OSTI2022

Sulfur-enhanced dynamics of coinage metal(111) surfaces: Step edges versus terraces as locations for metal-sulfur complex formation

The propensity of trace amounts of sulfur adsorbed on coinage metal(111) surfaces to dramatically enhance surface dynamics has been demonstrated by STM observations of accelerated 2D island decay for Cu and Ag. It is generally accepted that this enhancement is due to the formation of adsorbed metal-sulfur complexes, which facilitate surface mass transport of the metal. These complexes were originally proposed to form on terraces following the extraction of metal atoms from step edges and subsequent combination with sulfur on the terraces. However, even when thermodynamically feasible, this mechanism may not be kinetically viable for some complexes due to limited coupling of the complex concentration to the surface diffusion flux of metal atoms. Focusing on the case of Cu, we assess various scenarios where complexes are formed either on terraces or instead directly at step edges, the latter being a new paradigm. In this work, a new pathway is proposed for the formation on terraces. A rich variety of structures incorporating S at step edges exist, which could provide a viable source for complexes, at least from a thermodynamic perspective. However, it is necessary to also assess the activation barrier for complex formation and detachment from step edges. This is facilitated by the nudged-elastic-band analysis of the minimum energy path for this process utilizing machine-learning derived potentials based on density functional theory energetics for the metal-sulfur system.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Surface enrichment dictates block copolymer orientation

Orientation of block copolymer (BCP) morphology in thin films is critical to applications as nanostructured coatings. Despite being well-studied, the ability to control BCP orientation across all possible block constituents remains challenging. Here, in this study, we deploy coarse-grained molecular dynamics simulations to study diblock copolymer ordering in thin films, focusing on chain makeup, substrate surface energy, and surface tension disparity between the two constituent blocks. We explore the multi-dimensional parameter space of ordering using a machine-learning approach, where an autonomous loop using a Gaussian process (GP) control algorithm iteratively selects high-value simulations to compute. The GP kernel was engineered to capture known symmetries. The trained GP model serves as both a complete map of system response, and a robust means of extracting material knowledge. We demonstrate that the vertical orientation of BCP phases depends on several counter-balancing energetic contributions, including entropic and enthalpic material enrichment at interfaces, distortion of morphological objects through the film depth, and of course interfacial energies. BCP lamellae are found more resistant to these effects, and thus more robustly form vertical orientations across a broad range of conditions; while BCP cylinders are found to be highly sensitive to surface tension disparity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2021

On the structure sensitivity of and CO coverage effects on formic acid decomposition on Pd surfaces

Using density functional theory calculations, the Pd-catalyzed vapor-phase formic acid decomposition was studied, with a focus on the structure sensitivity and CO coverage effects. A comprehensive reaction network was developed on both the (111) and (100) facets of Pd, at CO coverages of 0 and 5/9 monolayer (ML). Pd(100) was determined to be more reactive than Pd(111) at both CO coverages. The introduction of 5/9 ML CO decreased the activity of both facets significantly, due to destabilization of the surface intermediates and transition states on the CO-decorated surfaces. Three reaction pathways were explored on the clean surfaces: the formate (HCOO) pathway, the carboxyl (COOH) pathway leading to the formation of CO 2 , and the COOH pathway leading to the formation of CO (COOH→CO). Based on the DFT-derived energetics alone, it appears that all three pathways contribute to the reaction on clean Pd, whereas the presence of 5/9 ML of CO inhibits the HCOO pathway on both facets and favors the COOH→CO pathway on the (111) facet, but the COOH→CO 2 one on the (100) facet. Moreover, at high CO coverages, alternative spectator CO-assisted adsorbate decomposition pathways were discovered, which could potentially play a role in formic acid decomposition on Pd catalysts under realistic reaction conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Simulations and Experiments Reveal Effect of Nanopores on Helium Diffusion in Quartz

The diffusion properties of noble gases in minerals are widely used to reconstruct the thermal histories of rocks. In this work, we combine density functional theory (DFT) calculations with laboratory experiments to investigate controls on helium diffusion in quartz. DFT calculations for perfect α-quartz predict substantially lower activation energies and frequency factors for helium diffusion than observed in laboratory experiments, especially in the [001] direction. These results imply that no helium could be retained in quartz at Earth surface temperatures, which conflicts with observations of partial cosmogenic 3 He retention over geologic time scales. Here, we implement a model of helium diffusion in α-quartz modulated by nanopore defects that disrupt energetically favorable diffusion pathways. In this model, we find that laboratory-determined diffusivities can be most closely reproduced when a helium atom encounters ~70 nanopore sites per million interstitial sites. The results of our model indicate that diffusion of helium in natural quartz, like other noble gases in other minerals, can be significantly modulated by extended defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Convergent ab initio analysis of the multi-channel HOBr + H reaction

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br–O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies −0.2 kcal mol−1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol−1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol−1 barrier. The final product energies lie −37.2, −67.9, and −7.3 kcal mol−1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin–orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol−1.

Chemistry↗

DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1993

Theoretical characterization of the reaction NH2+NO - products

The potential energy surface (PES) for NH2+NO is characterized using complete active space self-consistent-field (CASSCF)/derivative methods to locate the stationary points, followed by internally contracted configuration interaction (ICCI) calculations to determine the energetics. Production of N2+H2O is found to involve a complex mechanism, which, however, has no barrier with respect to NH2+NO. This pathway is exothermic by 124.5 kcal/mol. Production of NH2+OH can occur with no barrier other than the exothermicity from any of three isomers of HNNOH.

Walch, Stephen P.↗

NASA NTRS2010

Jovian Plasma Torus Interaction with Europa: 3D Hybrid Kinetic Simulation. First results

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa-moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements, (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy etal.,2007;Shematovichetal.,2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyro radius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions).Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization,electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma:(a) with O(++) ions; (b) with O(++) and S(++) ions. The majority of O2 atmosphere is thermal with an extended cold population (Cassidyetal.,2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.

Lipatov, A. S.↗

NASA NTRS2013

Jovian Plasmas Torus Interaction with Europa. Plasma Wake Structure and Effect of Inductive Magnetic Field: 3D Hybrid Kinetic Simulation

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa moon-magnetosphere system with respect to a variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo Orbiter mission, and for planning flyby and orbital measurements (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy et al., 2007; Shematovich et al., 2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyroradius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions). Photoionization, electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider the models with Oþ þ and Sþ þ background plasma, and various betas for background ions and electrons, and pickup electrons. The majority of O2 atmosphere is thermal with an extended non-thermal population (Cassidy et al., 2007). In this paper, we discuss two tasks: (1) the plasma wake structure dependence on the parameters of the upstream plasma and Europa's atmosphere (model I, cases (a) and (b) with a homogeneous Jovian magnetosphere field, an inductive magnetic dipole and high oceanic shell conductivity); and (2) estimation of the possible effect of an induced magnetic field arising from oceanic shell conductivity. This effect was estimated based on the difference between the observed and modeled magnetic fields (model II, case (c) with an inhomogeneous Jovian magnetosphere field, an inductive magnetic dipole and low oceanic shell conductivity).

Jovian magnetosphere↗

DOE OSTI2023

Acetylene Semi-Hydrogenation on a Perovskite Oxyhydride Surface: Insights from First Principles and Microkinetic Modeling

Although perovskite oxyhydrides (POHs) have shown promise as a support for selective hydrogenation, it is unknown if they themselves can catalyze alkyne semi-hydrogenation. Here, we use first-principles density functional theory, coupled with microkinetic modeling, to investigate acetylene semi-hydrogenation on a prototypical POH, cubic BaTiO 2.5 H 0.5 (BTOH). Two different mechanisms are examined on a representative surface of BTOH under the reaction conditions: both are based on the Horiuti–Polanyi mechanism, but the way of H 2 dissociation is different. In mechanism 1, a lattice hydride H atom and then a surface adsorbed H atom sequentially hydrogenate the adsorbed acetylene. In mechanism 2, two lattice hydride H atoms from the BTOH sequentially hydrogenate the adsorbed acetylene. In both mechanisms, the H atoms are replenished from gas phase H 2 dissociation. Using density functional theory (DFT), we have calculated pathways for these two mechanisms in hydrogenation of acetylene to ethylene, as well as further to ethane. Microkinetic modeling based on the DFT energetics indicates that at 523 K, hydrogenation to ethylene and ethane via mechanism 1 occurs at a rate up to two orders of magnitude faster than via mechanism 2. A selectivity analysis for the temperature range of 373–673 K shows that the product observed is essentially only ethylene within the more active mechanism 1. This is because at the steady state, there is a significant amount of surface anion vacancies that facilitate heterolytic H 2 dissociation and help stabilize a vinyl intermediate. In conclusion, the present findings suggest that POHs are capable of selective hydrogenation, thanks to their rich hydrogen surface-redox chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Enabling Particulate Materials Processing Science for High-Consequence, Small-Lot Precision Manufacturing

This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Effects of 2D filler on rheology of additive manufacturing polymers: Simulation and experiment on polyetherketoneketone-mica composites

Polymeric composites have attracted increasing interest in fundamental studies and industrial applications due to their often-superior characteristics over conventional pure polymers. In this work, we studied the relationship between the rheological properties, critical in processing methods such as additive manufacturing, and the structure of a polymer-2D filler system by experiment and simulation methods. Considering the high performance polymer polyetherketoneketone (PEKK) with 2D filler (mica) composite as an example, the effects of filler size, filler content, surface chemistry and shear rate were systematically investigated. Here, it is demonstrated that the interaction between the polymer and 2D filler brings two different consequences on tuning the viscosity of polymer-filler composite. The stronger interaction increases the shear stress of the polymer-2D filler composite system, leading to higher viscosity. On the other hand, the stronger interaction can promote more absorption of polymer chains on the surface, which makes the local polymer network disentangled and decreases the continuous viscosity from the polymer network. This work provides a fundamental understanding on the relationship between macroscopic rheological properties and the structural properties and energetics at molecular level of polymer-2D filler composite.

36 MATERIALS SCIENCE↗

NASA NTRS1979

On the connection of interplanetary shock wave parameters and energetic storm particle events

A series of energetic storm particle (ESP) events is investigated on the basis of observations of energetic solar protons by IMPs IV and V between 1967 and 1972. The total number and energy of particles associated with 20 ESP events are analyzed in connection with several plasma parameters at the shock surface; i.e., magnetic and thermal pressure, magnetic field and plasma mass flux. Correlation is established between the total energy of ESPs in the energy range of 1-80 MeV and the jump in the total pressure. It is pointed out that the amount of shock energy converted to the acceleration of ESP particles is not negligible and could, in principle, contribute to the deceleration of the shock wave.

Gombosi, T.↗

NASA NTRS1986

Observations of a large-scale vortex-like structure in the deep-tail plasma sheet boundary layer

ISEE-3 observations of a large-scale vortexlike structure in the deep tail of the magnetosphere at X(GSM) = -217 earth radii are reported. The structure is characterized by two clockwise rotations of the energetic-ion anisotropy vector. Variations in the magnetic-field vector approximately 180 deg out of phase with the ion variations are observed. This structure is most likely the signature within the magnetosphere of a surface wave at the magnetopause driven by a Kelvin-Helmholtz instability. Conditions inside and outside of the magnetosphere, as observed by ISEE-3 and ISEE-2, respectively, are examined; these conditions suggest that the surface wave is most likely propagating in the slow mode.

Sanderson, T. R.↗

NASA NTRS1989

Oxidation of silicon with a 5 eV O(-) beam

A silicon wafer has been oxidized at room temperature in vacuum using a pure, ground-state beam of O(-) ions. The beam was of sufficiently low energy that no displacement damage or implantation was energetically possible. The resulting SiO2 films were analyzed with X-ray photoelectron spectroscopy. A logarithmic dependence of oxide thickness on dose was observed, with an extrapolated oxidation efficiency of unity for the clean silicon surface. A distinct initial oxidation phase was observed, with an anomalously high level of silicon suboxides. In addition, the valence-band offset between the silicon and the oxide was unusually small, suggesting a large interfacial dipole.

Hecht, M. H.↗