SEARCH · Search NASA

Results for “sharp features”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Development of Novel Composite and Random Materials for Nonlinear Optics and Lasers

A qualitative model explaining sharp spectral peaks in emission of solid-state random laser materials with broad-band gain is proposed. The suggested mechanism of coherent emission relies on synchronization of phases in an ensemble of emitting centers, via time delays provided by a network of random scatterers, and amplification of spontaneous emission that supports the spontaneously organized coherent state. Laser-like emission from powders of solid-state luminophosphors, characterized by dramatic narrowing of the emission spectrum and shortening of emission pulses above the threshold, was first observed by Markushev et al. and further studied by a number of research groups. In particular, it has been shown that when the pumping energy significantly exceeds the threshold, one or several narrow emission lines can be observed in broad-band gain media with scatterers, such as films of ZnO nanoparticles, films of pi-conjugated polymers or infiltrated opals. The experimental features, commonly observed in various solid-state random laser materials characterized by different particle sizes, different values of the photon mean free path l*, different indexes of refraction, etc.. can be described as follows. (Liquid dye random lasers are not discussed here.)

Noginov, Mikhail↗

DOE OSTI2021

Investigation of hardening mechanisms and size effects in proton-irradiated HT-9 steels

Ferritic/martensitic steels, such as HT-9, are known for their complex microstructural features and mechanical properties. Here in this paper, in-situ micro-tensile tests and traditional fractography methods were utilized to study the fracture behavior of proton-irradiated HT-9 steels. First, to evaluate the viability of micro-tensile tests for nuclear material qualification process, meso-tensile tests on as-received HT-9 steels were performed. Fracture mechanisms of unirradiated HT-9 steels at both length scales were compared and underlying mechanisms discussed. The direct comparison of micro- and meso-scale data shows a distinctive size effect demonstrated by the increase in yield stress (YS). Upon completion of initial assessment, specimens were irradiated with 4 MeV + protons to three fluences, all of which were lower than 0.01 displacements per atom (dpa). As expected, the YS increases with irradiation. However, at 7 × 10 -3 dpa, the reversal of the trend was observed, and the YS exhibited sharp decline. We demonstrate that at lower length scales, grain structure has a more profound impact on the mechanical properties of irradiated materials, which provides information needed to fill in the gap in current understanding of the HT-9 fracture at different length scales.

36 MATERIALS SCIENCE↗

NASA NTRS2003

Irreversible Entropy Production in Two-Phase Mixing Layers

This report presents a study of dissipation (irreversible production of entropy) in three-dimensional, temporal mixing layers laden with evaporating liquid drops. The purpose of the study is to examine the effects of evaporating drops on the development of turbulent features in flows. Direct numerical simulations were performed to analyze transitional states of three mixing layers: one without drops, and two that included drops at different initial mass loadings. Without drops, the dissipation is essentially due to viscous effects. It was found that in the presence of drops, the largest contribution to dissipation was made by heating and evaporation of the drops, and that at large length scales, this contribution is positive (signifying that the drops reduce turbulence), while at small scales, this contribution is negative (the drops increase turbulence). The second largest contribution to dissipation was found to be associated with the chemical potential, which leads to an increase in turbulence at large scales and a decrease in turbulence at small scales. The next smaller contribution was found to be that of viscosity. The fact that viscosity effects are only third in order of magnitude in the dissipation is in sharp contrast to the situation for the mixing layer without the drops. The next smaller contribution - that of the drag and momentum of the vapor from the drops - was found to be negative at lower mass loading but to become positive at higher mass loading.

Okongo, Nora↗

NASA NTRS2016

Physical State and Distribution of Materials at the Surface of Pluto from New Horizons LEISA Imaging Spectrometer

From Earth based observations Pluto is known to be the host of N2, CH4 and CO ices and also a dark red material. Very limited spatial distribution information is available from rotational visible and near-infrared spectral curves obtained from hemispheric measurements. In July 2015 the New Horizons spacecraft reached Pluto and its satellite system and recorded a large set of data. The LEISA spectro-imager of the RALPH instruments are dedicated to the study of the composition and physical state of the materials composing the surface. In this paper we report a study of the distribution and physical state of the ices and non-ice materials on Pluto's illuminated surface and their mode and degree of mixing. Principal Component analysis as well as various specific spectral indicators and correlation plots are used on the first set of 2 high resolution spectro-images from the LEISA instrument covering the whole illuminated face of Pluto at the time of the New Horizons encounter. Qualitative distribution maps have been obtained for the 4 main condensed molecules, N2, CH4, CO, H2O as well as for the visible-dark red material. Based on specific spectral indicators, using either the strength or the position of absorption bands, these 4 molecules are found to indicate the presence of 3 different types of ices: N2-rich:CH4:CO ices, CH4-rich(:CO:N2?) ices and H2O ice. The mixing lines between these ices and with the dark red material are studied using scatter plots between the various spectral indicators. CH4 is mixed at the molecular level with N2, most probably also with CO, thus forming a ternary molecular mixture that follows its phase diagram with low solubility limits. The occurrence of a N2-rich - CH4-rich ices mixing line associated with a progressive decrease of the CO/CH4 ratio tells us that a fractionation sublimation sequence transforms one type of ice to the other forming either a N2-rich - CH4-rich binary mixture at the surface or an upper CH4-rich ice crust that may hide the N2-rich ice below. The strong CH4-rich - H2O mixing line witnesses the subsequent sublimation of the CH4-rich ice lag left behind by the N2:CO sublimation (N spring-summer), or a direct condensation of CH4 ice on the cold H2O ice (S autumn). The weak mixing line between CH4-containing ices and the dark red material and the very sharp spatial transitions between these ices and this non-volatile material are probably due to thermal incompatibility. Finally the occurrence of a H2O ice - red material mixing line advocates for a spatial mixing of the red material covering H2O ice, with possibly a small amount intimately mixed in water ice. From this analysis of the different materials distribution and their relative mixing lines, H2O ice appears to be the substratum on which other ices condense or non-volatile organic material is deposited from the atmosphere. N2-rich ices seem to evolve to CH4-dominated ices, possibly still containing traces of CO and N2, as N2 and CO sublimate away. The spatial distribution of these materials is very complex. The high spatial definition of all these composition maps, as well as those at even higher resolution that will be soon available, will allow us to compare them with Pluto's geologic features observed by LORRI panchromatic and MVIC multispectral imagers to better understand the geophysical processes in action at the surface of this astonishingly active frozen world.

Schmitt, B.↗

DOE OSTI2019

Injection of oxygen vacancies in the bulk lattice of layered cathodes

Surfaces, interfaces, and grain boundaries are, classically, known to be sinks at which defects generated within bulk lattice will preferentially migrate to and get annihilated, therefore rejuvenating a crystalline lattice to restore the intended functionalities.1, 2, 3 Here we report an inverse case by which the defects generated at the particle surface are continuously pumped into the bulk lattice and consequentially compromising the intended functionalities. We discovered that during the operation of a rechargeable battery, oxygen vacancies generated at cathode particle surfaces migrate toward inside of a perfect lattice in lithium-rich layered cathodes, featuring electrochemical driven surface vacancy creation, injection and condensation phenomenon and revealing an unrecognized surface process induced bulk lattice degradation featured with void formation and lattice reorganization. This process is associated with a high-cutoff voltage at which the anionic redox is activated. Coupling with calculations, we substantiate that triggering of anionic redox leads to sharply lowering of both formation energy of oxygen vacancy and migration barrier of oxidized oxide ions, therefore enabling migration of oxygen vacancy via oxidized oxide ion and leading to the injection of oxygen vacancy generated at the particle surface into the bulk lattice. This work unveils unpredicted behavior of oxidized oxide ion, leading to essential insight for designing of high capacity cathode for better batteries through coupled cationic and anionic redox dynamics in layered cathodes. In a broad term, this work provides insight on the intimate coupling between interfacial process and bulk lattice behavior for governing electrode functionality.

Yan, Pengfei↗

DOE OSTI2024

On the computation of moments in the Super-Transition-Arrays model for radiative opacity calculations

In the Super-Transition-Array statistical method for the computation of radiative opacity of hot dense matter, the moments of the absorption or emission features involve partition functions with reduced degeneracies, occurring through the calculation of averages of products of subshell populations. Here, in the present work, we discuss several aspects of the computation of such peculiar partition functions, insisting on the precautions that must be taken in order to avoid numerical difficulties. In a previous work, we derived a formula for supershell partition functions, which takes the form of a functional of the distribution of energies within the supershell and allows for fast and accurate computations, truncating the number of terms in the expansion. The latter involves coefficients for which we obtained a recursion relation and an explicit formula. We show that such an expansion can be combined with the recurrence relation for shifted partition functions. We also propose, neglecting the effect of fine structure as a first step, a positive-definite formula for the Super-Transition-Array moments of any order, providing an insight into the asymmetry and sharpness of the latter. The corresponding formulas are free of alternating sums. Several ways to speed up the calculations are also presented.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗

NASA NTRS2006

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗

DOE OSTI2018

Features of Short-Period Oscillations of Microwave Emission from the Solar Active Region before the Flare

The radiation of microwave sources above sunspots at a frequency of 17 GHz gives information about the parameters of solar plasma in the regions where the magneti-field strength is B ~ 2000 G in the transition region between the chromosphere and corona. Short-period oscillations (with a period of several minutes) of microwave emission from solar active regions (ARs) reflect wave processes in magnetic flux tubes of sunspots. Short-period oscillations of microwave emission from AR NOAA 12242 before two flares on December 17, 2014 are analyzed. This analysis is based on solar radio images obtained by means of the Nobeyama Radio Heliograph with a 10″−15″ two-dimensional spatial resolution. The radio maps of the whole solar disk were synthesized in a nonstandard mode with a cadence of 10 s and an averaging time of 10 s. An increase in the power of about ten-minute oscillations of microwave radiation approximately 40 to 50 min before the M1.5 flare (01: 00 UT) is found. On the same day, an increase in the power of ten-minute oscillations is observed about 60 min before the M8.7 flare (04: 42 UT). This effect is similar to the effect found earlier by two groups of authors independently for three-minute oscillations—namely, they observed a sharp increase in three-minute oscillations 15 to 20 min before the radio burst accompanying the flare. The effect in question may be interpreted as an relationship ofMHD waves propagating along the magnetic flux tube of a sunspot and the onset of the solar flare.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Overcharge investigation of large format lithium-ion pouch cells with Li(Ni 0.6 Co 0.2 Mn 0.2 )O 2 cathode for electric vehicles: Thermal runaway features and safety management method

The overcharge-induced thermal runaway features of large format commercial lithium-ion batteries with Li(Ni 0.6 Co 0.2 Mn 0.2 )O 2 (NCM622) cathode for electric vehicles under different current rates (C-rates) have been systematically studied at ambient temperature. The overcharge process is characterized as four stages. The temperature rise and the maximum temperature of the battery surface don't increase in proportion to the applied C-rates. However, with the increase of C-rates, the crest voltage of voltage curve rises linearly. When the voltage reaches approximately 5.1 V, a new voltage plateau appears in the cases below 2C. It is not sufficient that the temperature sensor is placed only near the terminal tab for most battery packs of EVs. In addition, the accumulated heat analysis demonstrates that side reactions dominate the temperature rise and contribute to most of the accumulated heat before thermal runaway. To mitigate the impact of overcharge and avoid the thermal runaway risk, a safety management method is proposed. Furthermore, the sharp drop in voltage before thermal runaway also provides a feasible approach to forewarn the users of the impending risk. These results are important for building safer batteries and providing information for the safety monitoring function of the battery management system (BMS).

25 ENERGY STORAGE↗

NASA NTRS2013

Downscaling Satellite Precipitation with Emphasis on Extremes: A Variational 1-Norm Regularization in the Derivative Domain

The increasing availability of precipitation observations from space, e.g., from the Tropical Rainfall Measuring Mission (TRMM) and the forthcoming Global Precipitation Measuring (GPM) Mission, has fueled renewed interest in developing frameworks for downscaling and multi-sensor data fusion that can handle large data sets in computationally efficient ways while optimally reproducing desired properties of the underlying rainfall fields. Of special interest is the reproduction of extreme precipitation intensities and gradients, as these are directly relevant to hazard prediction. In this paper, we present a new formalism for downscaling satellite precipitation observations, which explicitly allows for the preservation of some key geometrical and statistical properties of spatial precipitation. These include sharp intensity gradients (due to high-intensity regions embedded within lower-intensity areas), coherent spatial structures (due to regions of slowly varying rainfall),and thicker-than-Gaussian tails of precipitation gradients and intensities. Specifically, we pose the downscaling problem as a discrete inverse problem and solve it via a regularized variational approach (variational downscaling) where the regularization term is selected to impose the desired smoothness in the solution while allowing for some steep gradients(called 1-norm or total variation regularization). We demonstrate the duality between this geometrically inspired solution and its Bayesian statistical interpretation, which is equivalent to assuming a Laplace prior distribution for the precipitation intensities in the derivative (wavelet) space. When the observation operator is not known, we discuss the effect of its misspecification and explore a previously proposed dictionary-based sparse inverse downscaling methodology to indirectly learn the observation operator from a database of coincidental high- and low-resolution observations. The proposed method and ideas are illustrated in case studies featuring the downscaling of a hurricane precipitation field.

Hurricanes↗

DOE OSTI2023

Unusual magnetotransport in twisted bilayer graphene from strain-induced open Fermi surfaces

Anisotropic hopping in a toy Hofstadter model was recently invoked to explain a rich and surprising Landau spectrum measured in twisted bilayer graphene away from the magic angle. Suspecting that such anisotropy could arise from unintended uniaxial strain, we extend the Bistritzer–MacDonald model to include uniaxial heterostrain and present a detailed analysis of its impact on band structure and magnetotransport. We find that such strain strongly influences band structure, shifting the three otherwise-degenerate van Hove points to different energies. Coupled to a Boltzmann magnetotransport calculation, this reproduces previously unexplained nonsaturating B 2 magnetoresistance over broad ranges of density near filling ν = ± 2 and predicts subtler features that had not been noticed in the experimental data. In contrast to these distinctive signatures in longitudinal resistivity, the Hall coefficient is barely influenced by strain, to the extent that it still shows a single sign change on each side of the charge neutrality point—surprisingly, this sign change no longer occurs at a van Hove point. The theory also predicts a marked rotation of the electrical transport principal axes as a function of filling even for fixed strain and for rigid bands. More careful examination of interaction-induced nematic order versus strain effects in twisted bilayer graphene could thus be in order.

36 MATERIALS SCIENCE↗

NASA NTRS2007

Accounting for Uncertainties in Strengths of SiC MEMS Parts

A methodology has been devised for accounting for uncertainties in the strengths of silicon carbide structural components of microelectromechanical systems (MEMS). The methodology enables prediction of the probabilistic strengths of complexly shaped MEMS parts using data from tests of simple specimens. This methodology is intended to serve as a part of a rational basis for designing SiC MEMS, supplementing methodologies that have been borrowed from the art of designing macroscopic brittle material structures. The need for this or a similar methodology arises as a consequence of the fundamental nature of MEMS and the brittle silicon-based materials of which they are typically fabricated. When tested to fracture, MEMS and structural components thereof show wide part-to-part scatter in strength. The methodology involves the use of the Ceramics Analysis and Reliability Evaluation of Structures Life (CARES/Life) software in conjunction with the ANSYS Probabilistic Design System (PDS) software to simulate or predict the strength responses of brittle material components while simultaneously accounting for the effects of variability of geometrical features on the strength responses. As such, the methodology involves the use of an extended version of the ANSYS/CARES/PDS software system described in Probabilistic Prediction of Lifetimes of Ceramic Parts (LEW-17682-1/4-1), Software Tech Briefs supplement to NASA Tech Briefs, Vol. 30, No. 9 (September 2006), page 10. The ANSYS PDS software enables the ANSYS finite-element-analysis program to account for uncertainty in the design-and analysis process. The ANSYS PDS software accounts for uncertainty in material properties, dimensions, and loading by assigning probabilistic distributions to user-specified model parameters and performing simulations using various sampling techniques.

Nemeth, Noel↗

NASA NTRS1993

Flow of ices in the Ammonia-Water System

We have fabricated in the laboratory and subsequently deformed crystalline hydrates and partial melts of the water-rich end of the NH3-H2O system, with the aim of improving our understanding of physical processes occurring in icy moons of the outer solar system. Deformation experiments were carried out at constant strain rate. The range of experimental variables are given. Phase relationships in the NH3-H2O system indicate that water ice and ammonia dihydrate, NH3-2H2O, are the stable phases under our experiment conditions. X-ray diffraction of our samples usually revealed these as the dominant phases, but we have also observed an amorphous phase (in unpressurized samples only) and occasionally significant ammonia monohydrate, NH3-H2O. The onset of partial melting at the peritectic temperature at about 176 K appeared as a sharp transition in strength observed in samples of x(sub NH3) = 0.05 and 0.01, the effect of melt was less pronounced. For any given water ice + dihydrate alloy in the subsolidus region, we observed one rheological law over the entire temperature range from 175 K to about 140 K. Below 140 K, a shear instability similar to that occurring in pure water ice under the same conditions limited our ability to measure ductile flow. The rheological laws for the several alloys vary systematically from that of pure ice to that of dihydrate. Pure dihydrate is about 4 orders of magnitude less viscous than water ice just below the peritectic temperature, but because of a very pronounced temperature dependence in dihydrate (100 kJ/mol versus 43 kJ/mol for water ice) the viscosity of dihydrate equals or exceeds that of water ice at T less than 140 K. The large variation in viscosity of dihydrate with relatively small changes in temperature may be helpful in explaining the rich variety of tectonic and volcanic features seen on the surfaces of icy moons in the outer solar system.

Durham, W. B.↗

NASA NTRS1990

A vibrational analysis of the O2 (A 3Sigma/+/u) Herzberg I system using rocket data

An observation of the UV nightglow between 2670 and 3040 A was conducted over White Sands Missile Range on October 22, 1984. A 1/4-m spectrometer operating at 3.5-A resolution viewed the earth's limb at tangent heights between 90 and 110 km for 120 sec. A total of 41 spectral scans of the nightglow were obtained with the brightest feature being the O2 Herzberg I bands. The data were sorted into two groups, one from the top side of the layer and one containing the emission peak, and compared with synthetic spectra. The deduced vibrational distributions indicate that, at low altitudes, the higher vibrational levels (v-prime greater than 6) were relatively depleted; however, the magnitude of the vibrational shift is much less than that predicted from theories of vibrational relaxation. It is shown that increasing the electronic quenching with respect to the vibrational quenching can reduce the vibrational shift in the model and qualitatively explain the observations; however, several details of the vibrational distribution are not well reproduced.

Siskind, David E.↗

DOE OSTI2017

Mechanism and Kinetics for Reaction of the Chemical Warfare Agent Simulant, DMMP( g ), with Zirconium(IV) MOFs: An Ultrahigh-Vacuum and DFT Study

The mechanism and kinetics of interactions between dimethyl methylphosphonate (DMMP), a key chemical warfare agent (CWA) simulant, and Zr 6 -based metal organic frameworks (MOFs) have been investigated with in situ infrared spectroscopy (IR), X-ray photoelectron spectroscopy (XPS), powder X-ray diffraction (PXRD), and DFT calculations. DMMP was found to adsorb molecularly to UiO-66 through the formation of hydrogen bonds between the phosphoryl oxygen and the free hydroxyl groups associated with Zr 6 nodes on the surface of crystallites and not within the bulk MOF structure. Unlike UiO-66, the infrared spectra for UiO-67 and MOF-808, recorded during DMMP exposure, suggest that uptake occurs through both physisorption and chemisorption. The XPS spectra of MOF-808 zirconium 3d electrons reveal a charge redistribution following exposure to DMMP. In addition, analysis of the phosphorus 2p electrons following exposure and thermal annealing to 600 K indicates that two types of stable phosphorus-containing species exist within the MOF. DFT calculations, used to guide the IR band assignments and to help interpret the XPS features, suggest that uptake is driven by nucleophilic addition of an OH group to DMMP with subsequent elimination of a methoxy substituent to form strongly bound methyl methylphosphonic acid (MMPA). The rates of product formation indicate that there are likely two distinct uptake processes, requiring rate constants that differ by approximately an order of magnitude. However, the rates of molecular uptake were found to be nearly identical to the rates of reaction, which strongly suggests that the reaction rates are diffusion-limited. Here, the final products were found to inhibit further reactions within the MOFs, and these products could not be thermally driven from the MOFs prior to decomposition of the MOFs themselves.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗