Numerical study of damage in particulate composites during high-strain rate loading using novel damage model
Not provided.
SEARCH · Search NASA
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.
Not provided.
Damage in the microstructures of energetic materials (EMs), such as propellants and plastic bonded explosives (PBXs), can significantly alter their response to external loads. Both sensitization and desensitization can occur, causing concerns with safety and performance in the field; predictive models that connect damage and the sensitivity of EMs can enable design and provide confidence in their robustness and reliability. However, modeling of damage evolution is challenging for real microstructures of EMs; samples of damaged EMs are difficult to obtain, thereby hindering experiments and direct numerical simulations to determine the sensitivity of EMs at various stages of damage. Here, we develop an approach to generate synthetic, i.e., in silico produced, damaged microstructures for use in simulations to connect damage levels to sensitivity. The development of the present workflow to generate and impose varying levels of damage in microstructures, known as HEDS (Heterogeneous Energetic Material Damage Simulator), begins with a small set of images of damaged PBXs and combines a collection of deep neural network techniques to generate microstructures with varying levels of damage. By making the synthetic microstructures conform closely to those observed in available real, imaged microstructures, we develop an ensemble of damaged microstructures that can be used for in silico shock experiments. HEDS develops these microstructure ensembles as level set fields, which are directly employed in a sharp interface Eulerian hydrocode where shock simulations are performed to quantify the energy release rate from hotspot fields generated in the microstructure. These capabilities can be useful for the analysis and assessment of changes in the sensitivity of EMs and to design formulations that are less susceptible to damage-induced changes in sensitivity and performance.
A damage model is developed and validated with experimental data for the non-linear mechanical behavior of SiC/SiC composite materials in nuclear applications. Cyclic thermal and mechanical loading associated with neutron irradiation effects of these composites leads to wide-spread and progressive micro-cracking that leads to loss of thermal conductivity and further enhancement of thermo-mechanical damage. A physics-based model of wide-spread micro-cracking is developed within the thermodynamic framework of continuum damage mechanics. Evolution equations for damage parameters that describe the growth of continuum damage are developed, where the material variables are obtained from experiments. The model novelty is in coupling mechanical, thermal, and irradiation damage through a consistent thermodynamic framework, including loss of thermal conductivity due to the evolution of mechanically induced micro-cracks. A number of thermo-mechanical experiments were conducted to confirm model assumptions. The model is shown to be validated with out-of-pile experiments, and then implemented using commercial finite element code COMSOL to the fuel cladding problem with normal and high radiation dose cases.
This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.
We present a novel microcrack-damage theory for brittle solids under compression. Instead of using internal variables like zeroth, second or fourth rank damage tensors, the state of material damage is represented by an internal function that encapsulates the information of direction, density and size of microcracks. Just like other internal variables, the evolution of this state function must obey the second law of thermodynamics for arbitrary loading paths. This is done by casting the model in the framework of continuous hyperplasticity and enforcing a non-negative dissipation rate functional. The proposed framework offers predictions on the continuous evolution of microcrack density and the induced material anisotropy along with the macroscopic stress-strain curves. The use of continuous damage function grants the model significantly enhanced resolution in characterizing the direction-dependent response of cracked solids compared to classical models that are based on damage tensors. Two scenarios are considered in developing the theory, one assumes frictionless cracks and the other incorporates friction between crack surfaces. The results highlight that inelasticity, pressure dependence, and loading-unloading hysteresis exhibited by brittle solids are natural consequences of frictional microcracks. The proposed theory offers a generic and versatile framework to upscale micromechanical processes operating at individual crack scale to explain the macroscopic behavior of cracked solids.
Outdoor-aged photovoltaic modules that have experienced golf-ball sized hail are characterized by flash testing and photoluminescence, electroluminescence, thermal, and ultraviolet fluorescence imaging. These measurements were performed immediately after the hailstorm and then again after another two years and four years in the field. The modules originally experienced different ranges of hail damage with some having cracks on approximately half of the cells and some having few to no cracked cells. Flash testing shows that the hail-damaged modules with the most cracked cells have 3 to 4% less power than those with the least damage. After four more years in the outdoor climate of Golden, Colorado since the hailstorm, all of the modules still exhibit similar relative power and amounts of cell cracking, suggesting that the hail damage has only initially reduced power but has not yet increased the degradation rate.
In this study, a novel, to the best of our knowledge, method of wet chemical etching of sapphire workpieces (such as optics, wafers, windows, and cones), called the sapphire advanced mitigation process (or sapphire AMP), has been developed that exposes sub-surface mechanical damage created during the optical fabrication process and significantly enhances the surface laser damage resistance (> 2×) and mechanical strength (up to ~ 2.6×). Sapphire AMP involves first treating the workpiece with a mixture of sulfuric and phosphoric acid ([H 2 SO 4 ]:[H 3 PO 4 ]= 1 : 3) at 220°C, followed with phosphoric acid at 160°C, then with sodium hydroxide base (NaOH) and surfactant at 40°C, and finally with a high-pressure deionized water spray rinse. Sapphire AMP has been demonstrated on both A- and C-plane sapphire workpieces. The mechanism of this etch process involves the reaction of the sapphire (Al 2 O 3 ) surface with sulfuric acid (H 2 SO 4 ) forming aluminum sulfate [Al2(SO4)3], which has low solubility. The high phosphoric acid content in the first and second steps of sapphire AMP results in the efficient conversion of Al 2 (SO 4 ) 3 to aluminum phosphate (AlPO 4 ), which is very soluble, greatly reducing reaction product redeposition on the workpiece surface. Sapphire AMP is shown to expose sub-surface mechanical damage on the sapphire surface created during the grinding and polishing processes, whose etched morphology has either isotropic or anisotropic evolution depending on the nature of the initial surface damage. Sapphire AMP was also designed to remove the key known surface, laser absorbing precursors (namely, foreign chemical impurities, the fracture surface layer of preexisting sub-surface damage, and reaction product or foreign species redeposition or precipitation). Static and sliding indention induced surface microfractures on sapphire are shown after sapphire AMP to have a significant decrease in the fast photoluminescence intensity (a known metric for measuring the degree of laser damaging absorbing precursors). In addition, the onset of laser damage (at 351 nm 3 ns) on sapphire AMP treated workpieces was shown to increase in fluence from ~ 4 to >9.5J/cm 2 . Finally, biaxial ball-on-ring mechanical tests on sapphire disks showed an increase in the failure stress from 340 MPa (with pre-existing 28 µm flaws) to ~ 900MPa after sapphire AMP, which is attributed to the blunting of the surface microfractures.
Structured conical blockers manufactured on optical input surfaces are being used to shadow, and thereby, arrest growth of exit surface damage on large optics for applications in high-energy laser systems to extend their lifetime usage. However, as the size of the Shadow Cone Blocker (SCB) is increased, the deflected light must be carefully managed to prevent consequential damage from downstream intensification. Here, we design and fabricate a distinct input SCB geometry to alleviate this problem, enabling larger damage site mitigation. The demonstrated ~1 mm diameter SCB was manufactured with a convergent approach using a CO 2 pulsed laser ablative scanning sequence on a fused silica window. The SCB design was refined to maximize resistance to both input and exit surface damage initiations on 1 cm thick fused silica windows when exposed to 351 nm irradiation and validated with laser damage testing. The design showed to prevent damage onset to the exit surface for incident fluences on the SCB of 10.7 ± 1.3 J·cm -2 and is resistant to damage on the input surface exceeding 30 J·cm -2 input fluence.
Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.
Management and control of damage initiation and growth on high-value National Ignition Facility (NIF) optics are critical to its operation. Cone-shaped features are currently being used on the input surfaces to arrest growth of exit surface damage by creating “shadows.” Light refracted from the walls of the cones interferes with the incident beam at the exit surface to create an intensified annulus. Significant exit surface damage was observed at the locations of these annuli. Analysis showed that this damage is consistent with a recently reported new mechanism that degrades damage resistance by dispersing damage precursors over the surface upon exposure to UV laser radiation. This has significant implications for the design and operation of high-energy and power laser systems such as the NIF.
A high-strength dual a+ß phase titanium alloy Ti–6Al–4V is utilized as a material for beam windows in several accelerator target facilities. However, relatively little is known about how material properties of this alloy are affected by high-intensity proton beam irradiation. With plans to upgrade neutrino facilities at J-PARC and Fermilab to over 1 MW beam power, the radiation damage in the window material will reach a few displacements per atom (dpa) per year, significantly above the ~0.3 dpa level of existing data. The international RaDIATE collaboration, Radiation Damage In Accelerator Target Environments, has conducted a high intensity proton beam irradiation of various target and window material specimens at Brookhaven Linac Isotope Producer (BLIP) facility, including a variety of titanium alloys. Post-Irradiation Examination (PIE) of the specimens in the 1st capsule, irradiated at up to 0.25 dpa, is in progress. Tensile tests in a hot cell at Pacific Northwest National Laboratory (PNNL) exhibited a clear signature of radiation hardening and loss of ductility for Ti–6Al–4V, while Ti–3Al–2.5V, with less ß phase, exhibited less severe hardening. Microstructural investigations will follow to study the cause of the difference in tensile behavior between these alloys. High-cycle fatigue (HCF) performance is critical to the lifetime estimation of beam windows exposed to a periodic thermal stress from a pulsed proton beam. The first HCF data on irradiated titanium alloys are to be obtained by a conventional bend fatigue test at Fermilab and by an ultrasonic mesoscale fatigue test at Culham Laboratory. Specimens in the 2nd capsule, irradiated at up to ~1 dpa, cover typical titanium alloy grades, including possible radiation-resistant candidates. These systematic studies on the effects of radiation damage of titanium alloy materials are intended to enable us not only to predict realistic lifetimes of current beam windows made of Ti–6Al–4V, but also to extend the lifetime by choosing a more radiation and thermal shock tolerant alloy with a preferable heat treatment, or even by developing new materials.
The Facility for Rare Isotope Beams (FRIB) SRF linear heavy-ion accelerator is designed to accelerate all ions up to uranium to a maximum beam power of 400 kW. Several beam intercepting devices (BIDs) are essential to the successful operation of the accelerator, including a low power charge selector (LPCS) made of copper containing 0.15% precipitated aluminum oxide by weight called Glidcop Al-15. As FRIB ramps up the primary beam power beyond the current 20 kW level, the charge selector must withstand higher radiation damage rates, typically measured in displacements per atom (dpa). Significant beam induced radiation damage including significant deformation, like swelling, blistering, and cracking has been observed on a recently removed LPCS. We also observed physical features up to about 5.5 mm wide and appear to be deeper than the projected range of any ions. This paper presents the results of Monte Carlo simulations carried out using the Particle and Heavy Ion Transport code System (PHITS), quantifying the total damage dose and ion concentration, etc. Simulating an accurate irradiation history is essential to determining the scope of post irradiation examination (PIE) work.
Simple and complex clustered DNA damage represent the critical initial damage caused by radiation. In this paper, a multinomial probability model of clustered damage is developed with probabilities dependent on the energy imparted to DNA and surrounding water molecules. The model consists of four probabilities: (A) direct damage of sugar-phosphate moieties leading to SSB, (B) OH− radical formation with subsequent SSB and BD formation, (C) direct damage to DNA bases, and (D) energy imparted to histone proteins and other molecules in a volume not leading to SSB or BD. These probabilities are augmented by introducing probabilities for the relative location of SSB using a ≤10 bp criteria for a double-strand break (DSB) and for the possible success of a radical attack that leads to SSB or BD. Model predictions for electrons, 4He, and 12C ions are compared to the experimental data and show good agreement. Thus, the developed model allows an accurate and rapid computational method to predict simple and complex clustered DNA damage as a function of radiation quality and to explore the resulting challenges to DNA repair.
A comparative analysis is performed on the structural damage response and associated mechanisms in lanthanum aluminate and yttrium aluminate crystals under various irradiation conditions by a combination of experimental and theoretical approaches. Under low-energy Au+ irradiation, the damage accumulation curve shows a higher damage rate for LaAlO 3 crystals than YAlO 3 crystals. The relatively low irradiation tolerance of LaAlO 3 to the action of nuclear collisions is ascribed to the large amorphization cross-section and effective cross-section for defect-stimulated amorphization. Under swift Ar 12+ , Ni 19+ and Kr 17+ irradiation with different ion energies and velocities, the formed highly-disordered/amorphous latent tracks with different morphologies in pristine and predamaged crystals are discussed, and the corresponding electronic energy loss and lattice temperature thresholds are quantitatively determined. Compared to YAlO 3 , LaAlO 3 exhibits lower sensitivity and higher damage tolerance to the electronic energy loss process, attributing to its relatively high recrystallization efficiency during the rapid quenching process. Furthermore, the introduction of lattice defects into LaAlO 3 and YAlO 3 crystals considerably enhances the sensitivity and intensity of thermal spike response to the electronic energy loss, and the induced effective modification of track morphologies demonstrates the synergistic effect between the electronic energy loss and pre-existing defects created by nuclear collisions. In this case, even under the action of electronic energy loss below the threshold, the lattice temperature in the nuclear-collision damaged crystalline system could still meet the criterion for track production. The irradiation energy deposited to atoms and induced lattice temperature evolution discussed in this work provide a deeper insight into the complex processes involved in irradiation-induced latent track behaviors.
In this report, a new constitutive model is derived that combines the Wilshire equations with continuum damage mechanics (CDM) to enable the long-term prediction of creep deformation, damage, and rupture. The Wilshire equations have been demonstrated to accurately extrapolate the stress-rupture, minimum-creep-strain-rate, and time-to-strain of various alloys across decade of life. Recently, the time-to-creep-strain equation has been rearranged to predict creep deformation curves; however, the resulting equation is difficult to calibrate and does not track damage. The CDM-based Sinh constitutive model accurately predicts creep deformation, damage, and stress-rupture for a variety of alloys; however, it lacks an explicit description of stress and temperature co-dependence which limits interpolation and extrapolation ability. In this study, the Wilshire and the CDM-based Sinh equations are combined to create the “WCS” model. The WCS model combines the best features of each model while eliminating their deficiencies. Experimental data for alloy P91 is gathered and the WCS model is calibrated to the data. The WCS model accurately predicts the stress rupture, minimum-creep-strain-rate, creep deformation, and damage of alloy P91 across the experimental stress and temperature range as well as consistency in parametric simulations.
Multilayer dielectric (MLD) gratings with high diffraction efficiency and a high laser-induced damage (LID) threshold for pulse compressors are key to scaling the peak and average power of chirped pulse amplification lasers. However, surface defects introduced by manufacturing, storage, and handling processes can reduce the LID resistance of MLD gratings and impact the laser output. The underlying mechanisms of such defect-initiated LID remain unclear, especially in the femtosecond regime. In this Letter, we model dynamic processes in interactions of a 20-fs near-infrared (NIR) laser pulse and a MLD grating design in the presence of cylindrically symmetrical nodules and particle contaminants and cracks at the surface. Utilizing a dynamic model based on a 2D finite difference in time domain (FDTD) field solver coupled with photoionization, electron collision, and refractive index modification, we study the simulation results for the damage site distribution initiated by defects of various types and sizes and its impact on the LID threshold of the grating design.
The damage of approximately monolayer films of L-histidine by low-energy spin-polarized electrons (SPE) ejected from a magnetized cobalt substrate has been probed using X-ray photoelectron spectroscopy (XPS). Average damage cross sections for N-containing motifs of L-histidine are 25-30 and 2-5 Mb for zwitterions and neutral molecules, respectively. Here, the magnetization direction of the substrate, which controls the ejected SPE helicity, was reversed in situ, and statistically significant differences in the damage cross sections of 10-30% were measured between positive and negative electron helicities. This is the first measurement of spin dichroism (SD) in an amino acid. The differential cross sections suggest that inelastic scattering of SPE with chiral molecules could contribute to the persistence of one enantiomer vs the other under certain irradiation conditions, particularly for the zwitterionic species.
The Nuclear Science User Facilities (NSUF) provide researchers access to nuclear research facilities and expertise all over the United States and at the Belgian Center for Nuclear Research (SCK/CEN). One special capability of the NSUF is providing access to nuclear research and test reactors often used for conducting material testing experiments. These services are provided at no cost to the researcher through a competitive, peer-reviewed process. With so many proposals being submitted simultaneously, an NSUF information management resource called the Combined Materials Experiment Toolkit (CoMET) was implemented to speed up the review process and give users access to the information and tools necessary for designing meaningful experiments and writing the best proposals possible. One of CoMET’s tools is the Reactor Activation and Damage (RAD) Calculator. For those with little or no background in it, neutron irradiation in research and test reactors can be a challenging subject to include in a proposal. That is why the RAD Calculator was created to assist users with providing estimates of irradiation and post-irradiation conditions, enabling them to select the best neutron irradiation and post-irradiation examination facilities for reaching their experiment goals. The RAD Calculator is a general scoping tool that estimates radiation damage in terms of DPA, in addition to the resulting radioactivity of the material post-irradiation. This report evaluates the neutron damage component of the RAD Calculator. The evaluation of the activation component of the calculator is covered in report INL/EXT-20-58080.