Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inelasticity”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Phonons reveal coupled cholesterol-lipid dynamics in ternary membranes

Experimental studies of collective dynamics in lipid bilayers have been challenging due to the energy resolution required to observe these low-energy phonon-like modes. However, inelastic X-ray scattering (IXS) measurements — a technique for probing vibrations in soft and biological materials — are now possible with sub-meV resolution, permitting direct observation of low energy, phonon-like modes in lipid membranes. Here, in this study, IXS measurements with sub-meV energy resolution reveal a low-energy optic-like phonon mode at roughly 3 meV in the liquid-ordered (L o ) and liquid-disordered (L d ) phases of a ternary lipid mixture. This mode is only observed experimentally at momentum transfers greater than 5 nm -1 in the L o system. A similar gapped mode is also observed in all-atom molecular dynamics (MD) simulations of the same mixture, indicating that the simulations accurately represnt the fast, collective dynamics in the L o phase. It’s optical nature and the Q range of the gap together suggest that the observed mode is due to the coupled motion of cholesterol-lipid pairs, separated by several hydrocarbon chains within the membrane plane. Analysis of the simulations provide molecular insight into the origin of the mode in transient, nanoscale substructures of hexagonally packed hydrocarbon chains. This nanoscale hexagonal packing was previously reported based on molecular dynamics simulations and later by NMR measurements. Here, however, the integration of IXS and MD simulations identifies a new signature of the L o substructure in the collective lipid dynamics, thanks to the recent confluence of IXS sensitivity and MD simulation capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vacancy-driven variations in the phonon density of states of fast neutron irradiated nuclear graphite

Research examining radiation damage and its effects in graphite began in the 1940’s during the development of moderated nuclear reactors. Interest in this topic is expanding because of emerging applications associated with fullerenes and carbon nanostructures as well as its long-standing use as a fission moderator. In this work, we report the measurements of the full phonon density of states of irradiated nuclear graphite using inelastic neutron scattering experiments at room temperature. The samples were previously exposed to different levels of neutron doses and irradiation temperatures. The phonon density of states of perfect and defected graphite supercells including different configurations of vacancies and interstitials were calculated using first-principles direct method. Despite high neutron irradiation doses, the scattering measurements show that the induced damage is localized, and the layered structure is preserved. In conclusion, a comparison of the measured phonon densities of states of irradiated samples with those calculated of defected supercells indicates that the main changes observed are mainly attributed to formation of vacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Work function measurements of clean and modified carbon nanospikes

Carbon surfaces incorporating densely packed spikes with nanometer sharpness have shown exceptional electrocatalytic properties attributed to the enhanced electric field from the topography of the spikes. In this article, the composition of nanospike surfaces is identified by x-ray photoelectron spectroscopy (XPS) and their absolute work function determined using ultraviolet photoemission spectroscopy (UPS). Low temperature annealing of as-grown samples above 275 °C is required to produce a clean surface which has a 4.13 eV work function, a half volt lower than that of flat graphite. Treatments that contaminate the spiked surface, including exposure to air, oxidation at elevated temperature, or immersion in water or hydrocarbons, all increase the work function. Processing that blunts the spikes, including exposure to an oxygen plasma, argon sputtering, or high temperature annealing (800 °C) result in a work function close to that of flat graphite. An unusual double onset in the UPS secondary electron intensity is observed on as-grown nanospike samples and is reproduced by absorbing hydrocarbons on the surface. In conclusion, this double onset appears to be unique to carbon substrates and may originate in inelastic scattering of photoelectrons.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A multiscale cohesive law for carbon fiber networks

Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. In this study, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Cam-Clay plasticity. Part IX: On the anisotropy, heterogeneity, and viscoplasticity of shale

Here, we investigate three aspects of material behavior in this ninth installment of the Cam-Clay series, namely, anisotropy, heterogeneity, and viscoplasticity. The main focus of the paper is creep in shale and how this time-dependent deformation behavior may be quantified across the scales, from nanometers to millimeters. Recognizing the highly heterogeneous nature of shale, we adopt a simplified representation of this material as a mixture of softer matter representing organics and clay, and harder matter representing the inorganic rock matrix. Due to the presence of bedding planes in the rock matrix, anisotropy in both the elastic and inelastic responses is assumed; yet, the superimposed softer matter may be taken to be isotropic unless experimental evidence indicates otherwise. Viscoplasticity is considered for the first time in this series of work, in which both the Duvaut-Lions and Perzyna formulations are utilized. The two viscoplastic formulations are shown to predict very similar time-dependent deformation responses. The framework is used to interpret the results of multiscale triaxial laboratory creep tests in organic-rich Barnett shale. Time-dependent strain localization in the form of dilative shear bands is also shown to result from loading with varying strain rates, as well as from creep processes that accommodate the multiscale heterogeneity of shale.

42 ENGINEERING↗
DOE OSTI2022

Role of electronic energy loss on defect production and interface stability: Comparison between ceramic materials and high-entropy alloys

High-entropy alloys (HEAs) and some complex alloys exhibit desirable properties and significant structural stability in harsh environments, including possible applications in advanced reactors. Energetic ion irradiation is often used as a surrogate for neutron irradiation; however, the impact of ion electronic energy deposition and dissipation is often neglected. Moreover, differences in recoil energy spectrum and density of cascade events on damage evolution must also be considered. In many chemically complex alloys, the mean free path of electrons is reduced significantly, thus their decreased thermal conductivity and slow dissipation of localized radiation energy can have noticeable effects on displacement cascade evolution that is greatly different from metals with high thermal conductivity. In this work, nanocrystalline HEAs of Ni 20 Fe 20 Co 20 Cr 20 Cu 20 and nonequiatomic (NiFeCoCr) 97 Cu 3 , both having much lower room-temperature thermal conductivity than pure Ni or Fe, are chosen as model HEAs to reveal the role that electronic energy loss during ion irradiation has in complex alloys. The response of nanocrystalline HEAs is investigated under irradiation at room temperature using MeV Ni and Au ions that have different ratios of electronic energy to damage energy, which is the energy dissipated in displacing atoms. Different from previously reported amorphization of nanocrystalline SiC, experimental results on these HEAs show that, similar to the process in nanocrystalline oxide materials, both inelastic thermal spikes via electron–phonon coupling and elastic thermal spikes via collisions among atomic nuclei contribute to the overall grain growth. The growth follows a power law dependence with the total deposited ion energy, and the derived value of the power-exponent suggests that the irradiation-induced instability at and near grain boundaries leads to local rapid atomic rearrangements and consequently grain growth. The high power-exponent value can be attributed to the sluggish diffusion and delayed defect evolution arising from the chemical complexity intrinsic to HEAs. Here, this work calls attention to quantified fundamental understanding of radiation damage processes beyond that of simplified displacement events, especially in simulating neutron environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Efficient Fourier transforms for transverse momentum dependent distributions

Hadron production at low transverse momenta in semi-inclusive deep inelastic scattering can be described by transverse momentum dependent (TMD) factorization. This formalism has also been widely used to study the Drell-Yan process and back-to-back hadron pair production in $e^+e^-$ collisions. These processes are the main ones for extractions of TMD parton distribution functions and TMD fragmentation functions, which encode important information about nucleon structure and hadronization. One of the most widely used TMD factorization formalism in phenomenology formulates TMD observables in coordinate $b_\perp$-space, the conjugate space of the transverse momentum. The Fourier transform from $b_\perp$-space back into transverse momentum space is sufficiently complicated due to oscillatory integrands that it requires a careful and computationally intensive numerical treatment in order to avoid potentially large numerical errors. Within the TMD formalism, the azimuthal angular dependence is analytically integrated and the two-dimensional $b_\perp$ integration reduces to a one-dimensional integration over the magnitude $b_\perp$. In this paper we develop a fast numerical Hankel transform algorithm for such a $b_\perp$-integration that improves the numerical accuracy of TMD calculations in all standard processes. Libraries for this algorithm are implemented in Python 2.7 and 3, C++, as well as FORTRAN77. All packages are made available open source.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

LeptonInjector and LeptonWeighter: A neutrino event generator and weighter for neutrino observatories

In this work, we present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples. Program summary: Program Titles: LeptonInjector and LeptonWeighter CPC Library link to program files: https://doi.org/10.17632/662gkpjfd9.1 Developer's repository links: https://github.com/icecube/LeptonInjector and https://github.com/icecube/LeptonWeighter Licensing provisions: GNU Lesser General Public License, version 3. Programming Language: C++11 External Routines: • Boost • HDF5 • nuflux (https://github.com/icecube/nuflux) • nuSQuIDS (https://github.com/arguelles/nuSQuIDS) • Photospline (https://github.com/icecube/photospline) • SuiteSparse (https://github.com/DrTimothyAldenDavis/SuiteSparse) Nature of problem: LeptonInjector: Generate neutrino interaction events of all possible topologies and energies throughout and around a detector volume. LeptonWeighter: Reweight Monte Carlo events, generated by a set of LeptonInjector Generators, to any desired physical neutrino flux or cross section. Solution method: LeptonInjector: Projected ranges of generated leptons and the extent of the detector, in terms of column depth, are used to inject events in and around the detector volume. Event kinematics follow distributions provided in cross section files. LeptonWeighter: Event generation probabilities are calculated for each Generator, which are then combined into a generation weight and used to calculate an overall event weight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Electronic structure of reduced CeO 2 (111) surfaces interacting with hydrogen as revealed through electron energy loss spectroscopy in comparison with theoretical investigations

Based on both, ab-initio cluster calculations as well as periodic density functional based loss function calculations, we have assigned the origins of the valence electron excitation regime in electron energy loss spectra of well-ordered ceria films in (111) orientation at various states of reduction as well as after exposing the reduced films to hydrogen from the gas phase (Li et al, Angew. Chem. Int. Ed., 58 (2019) 14686–14693). The explicit calculation of intensity distributions using the dipole approximation allow us to draw conclusions about the nature of oxygen vacancies, which occur upon reduction and how those interact with hydrogen to form hydride species. Herein, we also reports a brief discussion of vibrational excitations of hydrogen loaded ceria and corroborates previous interpretations, based on inelastic neutron scattering.

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

Real-space Green’s function approach to photoelectron diffraction

We discuss the real-space Green’s function (RSGF) approach for core-level photoelectron diffraction (PD) and its comparison to x-ray absorption spectra (XAS). We focus especially on the separable Green’s function formalism for efficient calculations of high-order multiple scattering. Computational details such as scattering potentials, self-energy effects, Debye–Waller factors, and inelastic losses are discussed. Finally, we consider prospects for improving ab initio PD calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Measurement of mechanical and fracture properties of solid electrolyte interphase on lithium metal anodes in lithium ion batteries

Mechanical integrity of the solid electrolyte interphase (SEI) plays an essential role in determining the life and performance of lithium-ion batteries. Fracture and continued formation of the SEI contribute to consumption of lithium, drying of electrolyte, increase in impedance, and growth of dendrites resulting in capacity fade and premature failure. Electrolyte additives such as fluoroethylene carbonate (FEC) have been known to improve performance, but the underlying reasons have been elusive. Despite its importance, reliable methods for mechanical characterization of SEI have been lacking. In this paper, we present a new experimental technique that combines atomic force microscopy and membrane-bulge configuration to accurately measure the stress-strain behavior of SEI, including the onset of inelastic response and evolution of fracture. We characterize the SEI formed with two ethylene carbonate-based electrolytes, without and with fluoroethylene carbonate (FEC) additive. The measurements show a striking contrast; SEI with FEC additive has 80% higher elastic modulus and a vastly higher resistance to fracture. These findings offer a mechanical-behavior based rationale to understand how SEI controls battery performance. Moreover, the experimental technique offers a robust diagnostic tool to design electrolytes that can form SEI with the desired mechanical properties for optimal battery performance.

25 ENERGY STORAGE↗
DOE OSTI2020

Fluorination effect for stabilizing cationic and anionic redox activities in cation-disordered cathode materials

We see that cation-disordered Li-excess cathodes with oxygen redox reactions are promising candidates for high-energy-density Li ion batteries. Nevertheless, the oxygen redox process that is required for the high capacity often comes with the oxygen loss, which leads to severe capacity degradation and voltage decay. In this work, we have successfully synthesized a series of Li-excess cation-disordered cathodes (Li 1.2 Mn 0.4+x Ti 0.4-x O 2-x Fx) (0 ≤ x ≤ 0.2) with different fluorine (F) contents. The electrochemical performance results show that the Li 1.2 Mn 0.55 Ti 0.25 O 1.85 F 0.15 (LMTOF0.15) exhibits the highest reversible capacity (275 mAh g -1 , under 30 mA g -1 ), cyclability, and voltage retentions. The mapping of resonant inelastic X-ray scattering (mRIXS) and differential electrochemical mass spectroscopy (DEMS) results reveal that the fluorination enhances the reversible lattice oxygen redox reaction while suppressing irreversible gas release and surface reactions. The X-ray Absorption Spectroscopy (XAS) during the initial two cycles shows that F-substitution alleviates the reduction of the Mn valence state during the whole (dis)charge processes in the bulk and at the surface of the material, results in higher average discharge voltage. In addition, the introduction of F improves the structural stability and suppresses local lattice distortion of the material. Therefore, LMTOF 0.15 is able to cycle with smaller polarization, less interfacial side reaction and Mn dissolution, and therefore results in enhanced cyclability. This work provides a comprehensive understanding of the fluorination effect on the cationic and anionic redox activities in cation-disordered Li-excess cathodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Assessment of thermal and radiation induced creep in the dual cooled lead lithium blanket

The creep in recently designed blanket for the fusion nuclear science facility is described by decomposing it into thermal, irradiation and cavity swelling creep. Initial estimates of thermal creep using elastic high temperature rules showed that the recent heat transfer improvements at the first wall (FW) is feasible, and the FW can withstand heat flux of 0.35 MW/m 2 at 3000 hrs without rupture. A viscoplastic model that is based on Norton's power law and relies on Multiphysics coupling of solid mechanics and heat transfer modules is used to capture the inelastic deformation in the blanket. Results showed that the design peak heat flux of 0.25 MW/m 2 produced maximum thermal creep strain of 0.45 % and the relaxation of the thermal stress at the FW. Irradiation creep is prescribed to be proportional to the displacement damage dose and applied stress. The displacements from irradiation creep radially decrease from the FW to the back wall, but the maximum deformation is found at the back wall that is connected to the Helium manifold due to the high stresses at the region. In conclusion, cavity swelling creep is dominant at steady state and combined with radiation creep to produce displacement of about 8 mm at the FW.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Redox and structural controls on tin isotopic fractionations among magmas

Recent analytical developments have made possible the determination of the isotopic composition of tin in igneous rocks. In order to establish a framework to interpret the mass-dependent tin isotopic signatures of planetary materials, in this work, seven geologically-relevant silicate glasses (basalt, rhyolite, enstatite and anorthite glasses) were synthesized with moderate amounts of 119 Sn (on the order of a weight percent). Redox conditions were controlled during sample synthesis to set the redox ratio (Sn 2+ /Sn tot ) from stannous (Sn 2+ ) to stannic (Sn 4+ ) glasses. The mean force constants of tin bonds in these glasses were determined by synchrotron nuclear resonant inelastic X-ray scattering (NRIXS) in order to determine the reduced isotopic partition function ratios (β-factors) of these glasses. Clues on the coordination chemistry and the valence state of tin in these glasses were also derived from synchrotron Mossbauer spectroscopy (SMS). The force constants of tin drastically increases from Sn 2+ -bearing to Sn 4+ -bearing glasses and varies significantly with the glass composition at a given redox state. The average coordination number of tin likely controls these variations with glass composition as suggested by SMS results. It is concluded that large isotope fractionation is expected between materials containing Sn 2+ and Sn 4+ respectively even at magmatic temperatures and that the coordination chemistry of tin in silicates strongly affect its isotope partitioning behavior. Our experimental data are finally used to interpret available Sn isotope data collected in terrestrial rocks. The incompatible behavior of Sn 4+ in mantle minerals leads to a enrichment in heavy isotopes in mantle melts and to the depletion in heavy isotopes in solid residues of melting with a magnitude consistent with the isotope fractionation between Sn 2+ and Sn 4+ predicted by NRIXS data. Finally, we show that during fractional crystallization of basalt and considering the effect of tin coordination number in minerals and melts, the partitioning of Sn 4+ into ilmenite leads to an enrichment in light isotopes in the residual melt.

58 GEOSCIENCES↗
DOE OSTI2021

Iron, magnesium, and titanium isotopic fractionations between garnet, ilmenite, f ayalite, biotite, and tourmaline: Results from NRIXS, ab initio , and study of mineral separates from the Moosilauke metapelite

Interpreting isotopic signatures documented in natural rocks requires knowledge of equilibrium isotopic fractionation factors. Here, we determine equilibrium Fe isotope fractionation factors between several common rock-forming minerals using a comparative approach involving three independent methods: (i) isotopic analyses of natural minerals from a metapelite from Mt. Moosilauke, New Hampshire, for which equilibration temperature and pressure are well constrained to be near the aluminosilicate triple point (T ≃ 500 °C, P ≃ 4 kbar), (ii) Nuclear Resonant Inelastic X-ray Scattering (NRIXS) measurements of Fe force constants of minerals, and (iii) Density Functional Theory (DFT) ab initio calculations of Fe force constants of minerals. The minerals studied for Fe isotopes include, in increasing order of their β-factors: garnet < ilmenite ≈ fayalite < biotite < tourmaline < muscovite ≈ plagioclase. Some of this ordering is affected by the presence of Fe 3+ in the minerals, which tends to form stiffer bonds and be associated with heavy Fe isotope enrichments relative to Fe 2+ . We are, however, able to assess the magnitude of the effect of the ratio Fe 3+ /ΣFe on equilibrium fractionation factors, notably on the ilmenite-hematite solid solution. Equilibrium Fe isotopic fractionation factors between garnet, ilmenite, biotite, tourmaline and fayalite are determined. We also report Mg and Ti isotopic compositions of selected Moosilauke minerals that allow us to better constrain the equilibrium fractionation factors for garnet-biotite-tourmaline (Mg isotopes) and biotite-ilmenite (Ti isotopes). We show how the newly determined equilibrium fractionation factors can be used to address diverse problems in Earth and planetary sciences, notably (i) Fe and Mg isotopic fractionation during anatexis, (ii) Fe isotopic fractionation in lunar ilmenite, and (iii) Ti isotopic fractionation during fluvial transport of minerals.

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

Equilibrium Fe isotope fractionation between olivine, pyroxene, spinel and MORB glass: Implications for mantle partial melting to generate MORBs

Primitive mid-ocean ridge basalts (MORBs) exhibit Fe isotopic compositions heavier than the upper mantle by +0.074 ± 0.028 ‰ for δ 56 Fe. The processes responsible for this isotopic difference remain unclear. Modeling of Fe isotope fractionation during mantle partial melting requires reliable equilibrium Fe isotope fractionation factors between minerals and melts, for which consistent data are still lacking. Here, in this study, we used Nuclear Resonant Inelastic X-ray Scattering (NRIXS) technique to measure Fe force constants for a MORB glass (ALV 519-4-1) and natural mantle minerals (olivine, orthopyroxene, clinopyroxene, and spinel) to determine the equilibrium Fe isotope fractionation factors between them. The force constants determined in this study, in increasing order, are 167 ± 26 N/m for spinel, 175 ± 17 N/m for olivine, 176 ± 20 N/m for MORB glass, 205 ± 26 N/m for clinopyroxene, and 219 ± 36 N/m for orthopyroxene. We evaluated the previously proposed mechanisms for the heavy Fe isotopic composition of MORBs, including (i) mantle partial melting, (ii) mantle lithological heterogeneity, with pyroxenite in the source, (iii) mantle metasomatism by low-degree melts, and (iv) fractional crystallization of olivine from melts. For (i), we used the pMELTS program to simulate adiabatic decompression melting of mantle peridotites, and calculated Fe isotope fractionation based on Fe 3+ –Fe 2+ equilibrium-controlled fractionation, where Fe 3+ forms stronger bonds and is more incompatible than Fe 2+ . At 10 wt% peridotite melting, corresponding to MORB generation, only +0.03 ‰ Fe isotope fractionation between the melt and the original bulk composition (Δ 56 Fe = δ 56 Fe melt - δ 56 Fe 0 ) was produced, insufficient to account for the observed MORB-upper mantle difference. For (ii), melting of pyroxenites yields smaller Fe isotope fractionation than melting of peridotites, making it unlikely the cause for the MORB-upper mantle isotopic difference. For (iii), both the Fe 3+ /ΣFe ratio and the δ 56 Fe of melts increase with the degree of partial melting, indicating that low-degree melts are not isotopically heavy enough to significantly alter the isotopic composition of lithospheric mantle through metasomatism. For (iv), equilibrium isotope fractionation between olivine and melt is near zero. These results suggest that equilibrium Fe isotope fractionation alone cannot explain the MORB isotopic signature, highlighting the potential role of kinetic isotope fractionation. Using a diffusion model, we calculated kinetic Fe and Mg isotope fractionations associated with (iv) olivine crystallization from a melt, and found that the predicted Fe and Mg isotope fractionations were inconsistent with observations in MORBs. Qualitatively, two processes could have induced kinetic Fe isotope fractionation during MORB generation: (a) Fe-Mg interdiffusion between melt and solid during melt migration and (b) reactive melt-rock interactions during melt focusing. However, a quantitative understanding of their role in modifying the melt isotopic composition remains limited and requires further investigation.

Fe isotopes↗
DOE OSTI2020

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets—Part I: Mesoscale simulations

Recently a mesoscale model was developed that includes thermomechanically consistent constitutive equations with dissipation due to porous compaction, inelastic distortional deformation rate, and slipping on persistent dry joint surfaces. In this Part I of a two-part paper, we describe enhancements of the mesoscale model which include: normal compliance with hysteresis, non-persistent joints with finite area, joint sets with angular deviation, and the weakening effect of fluid pressure in the joints. Examples study the influences of these effects on average strength. These computationally expensive simulations are used to inspire functional forms in a continuum model, described in Part II, which is needed for large scale simulations of real applications.

42 ENGINEERING↗
DOE OSTI2020

A Physically-Based Model for Thermo-Oxidative and Hydrolytic Aging of Elastomers

A computationally efficient model is proposed to capture the loss of mechanical performance due to chemical aging that are formed as the competition of chain scission and cross-link formation/dissolution, such as thermo-oxidative aging or hydrolytic aging. The model should be considered an extension of our recent models [1, 2, 3] which further simplifies the matrix behavior based on the assumption of independence of environmental and mechanical damage. The model uses this assumption to reduce the necessary material parameters needed to model constitutive and inelastic behavior of elastomers during aging. To this end, the model can provide accurate predictions of the material performance with the significantly fewer number of fitting parameters. The model is relevant for all decay mechanisms formed by the occurrence of two simultaneous micro-mechanisms; (i) formation/reduction of the cross-links, and (ii) chain scission, both of which are present in thermo-oxidation and hydrolytic aging. Assuming the alteration of the chain density along the aging trajectory is identical to the peroxide cross-link density for thermo-oxidation, and the change of the average molecular weight for hydrolysis, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time and aging temperature. Next, the modified network alteration model is formulated for implementation into Finite Element (FE) simulations. The model is built on the presumption of homogeneous and consistent oxygen/water absorption and thus is mainly relevant for relatively thin samples exposed to environmental loads for a long time. The proposed model includes only six physically inspired material parameters. Therefore, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is the right choice for advanced implementations in FE programs.

42 ENGINEERING↗