Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response simulation”

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Comparison of irradiation tolerance of two MAX phases-Ti 4 AlN 3 and Ti 2 AlN

The microstructures of two MAX phases, Ti 2 AlN and Ti 4 AlN 3 , irradiated at room temperature with 70 keV He ions over a wide fluence range have been studied to reveal the dependence of MAX phase irradiation tolerance on structure. Grazing incident x-ray diffraction (GIXRD) reveals that a new set of diffraction peaks appears after ion irradiation for both MAX phases and the positions of diffraction peaks gradually shift towards lower 2θ values with increasing fluence. Transmission electron microscope characterizations show both MAX phases experience a similar structure transformation at low fluences, forming an intermediate structure. The d-spacings derived from electron diffraction analysis of intermediate structures are basically consistent with that of new peaks in GIXRD. At high fluences a nano-twinned microstructure is found in Ti 2 AlN, but no such microstructure is observed in Ti 4 AlN 3 . Based on the simulation results of electron diffractions and the irradiation responses of other MAX phases, the structure characteristics and formation mechanisms of the intermediate structure are derived and the discrepancy in the phase stability of both MAX phases at high fluences is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Strength of metals in liquid and solid states at extremely high tension produced by femtosecond laser heating

We will discuss results of combined experimental and theoretical investigations of ablation and laser-driven shock-wave phenomena in metal films irradiated by femtosecond laser pulses. The femtosecond interferometric microscopy technique was used to make time-resolved measurements of optical properties as well as record the deformation dynamics at both the rear and frontal surfaces during initial two-temperature electron-ion relaxation and subsequent hydrodynamic expansion. In conjunction with experiment, the formation and propagation of strong tensile and compression waves were investigated by a combination of two-temperature hydrodynamic modeling and molecular dynamics simulations. The experimental tensile strengths of aluminum and nickel in solid and liquid states at extremely high strain rates in range 10{sup 8} Division-Sign 10{sup 9}s{sup -1} were obtained from the time evolution of rear and frontal surface velocities. Theoretical tensile strengths calculated by atomistic simulations of ablation and spallation using micron-sized films agree well with experiment. Elastic-plastic response of metallic films to shock compression investigated by both experiment and theory/modeling will also be discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Smooth periodic gauge satisfying crystal symmetry and periodicity to study high-harmonic generation in solids

Intense lasers can easily drive nonadiabatic transitions of excited electron wave packets across the Brillouin zones, thus transition dipole moments (TDM) between energy bands of solids should be continuous, satisfying crystal symmetry, and periodic at zone boundaries. While current ab initio algorithms are powerful in calculating band structures of solids, they all introduced random phases into the eigenfunctions at each crystal momentum k. In this work, we show how to choose a “smooth-periodic” gauge where TDMs can be smooth versus k, preserving crystal symmetry, as well as maintaining periodic at boundaries. The symmetry properties of TDMs with respect to k ensure the absence of even-order harmonics from MgO with inversion symmetry, while the TDM in the “smooth-periodic” gauge for broken-symmetry ZnO is responsible for even harmonics that were underestimated in previous simulations. These results reveal the importance of correctly treating the complex TDMs that satisfy crystal symmetry and continuous across zone boundaries in nonlinear laser-solid interactions, which has been elusive in most theories so far.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Embedded Fiber Bragg Grating (FBG) Sensors Fabricated by Ultrasonic Additive Manufacturing for High-Frequency Dynamic Strain Measurements

This article demonstrates high-frequency dynamic strain measurements using fiber Bragg grating (FBG) sensors embedded in metal parts. Using an ultrasonic additive manufacturing (UAM) process, FBGs inscribed in polyimide coated optical fibers were embedded in aluminum parts. An electromagnetic shaker was used to exert dynamic events on the embedded FBG sensors with frequencies from 1 to 10 kHz. The high-speed interrogation of FBG sensors was accomplished using a tunable vertical-cavity surface-emitting laser (VCSEL) and a high-speed interrogation system sampling at 120 kHz. The strain response measured by the FBG sensors was compared with real-time measurements using a laser velocimeter. A finite-element analysis (FEA) was performed to simulate responses to both static strain and high-frequency dynamic strain. In conclusion, results show that strains as small as 2.5$με$ can be resolved at frequencies up to 10 kHz.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Mechanochemical Genesis of Prebiotics (18-LW-36) (Final Report)

In this project we explored the athermal condensation of simple amino acids towards the formation of peptides under high-pressure conditions. Our main aim was to address the provocative question that science continues explore: “Does life exists beyond Earth?”. We followed a combined experimental and theoretical investigation of simple amino acids under static high pressure and also under simultaneous application of high pressure and plastic shear. Surprisingly, we find that the ambient crystal structure of a-glycine has an extensive pressure resilience and stability up to at least 50 GPa. In the case glycylglycine (the simplest polypeptide) an isosymmetric phase transition was observed at ~7GPa. In the framework of this project a novel QMD approach was developed aiming to simulate the response of organic, molecules (with focus on glycine) under the simultaneous application of high pressure and plastic shear (V-RDAC). The V-RDAC simulations predict that glycine readily reacts when subject to high strain rate shearing under applied compressive loads greater than 7 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Strain-path change tests and physically based polycrystalline modeling of the behavior of recrystallized zirconium alloys - 2016-0092

Zirconium alloys used as cladding tubes for the fuel of pressurized water reactors can undergo high applied stress during power transients. In these conditions, biaxial loading may lead to plastic deformation of the cladding. A comprehensive understanding of the material mechanical response during loading-path change tests is an important step toward predicting the behavior in these specific conditions. Using a nonstandard mechanical testing machine, we conducted loading-path change tests at 623 K on as-received recrystallized Zircaloy-4 tubes. These tests consisted of an axial tensile loading and unloading followed by an internal pressure (or pure hoop tension) loading and unloading. These tests examined the kinematic and isotropic hardening components of the strain hardening behavior of the thin cladding tubes. The isotropic hardening was attributed to dislocation multiplication and dislocation-dislocation interactions. The kinematic hardening was attributed to the interaction of the grains with each other. A polycrystalline model was also used and improved in order to simulate the tests. A good prediction of the isotropic and kinematic hardening was provided by the modeling. (authors)

36 MATERIALS SCIENCE↗
NASA NTRS1991

Long discontinuous fiber composite structure: Forming and structural mechanics

Cost effective composite structure has motivated the investigation of several new approaches to develop composite structure from innovative material forms. Among the promising new approaches is the conversion of planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. In the present study, the authors have established a framework which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. The initial study focuses upon the establishment of micromechanics models for prediction of the effective anisotropic viscosities of the oriented fiber assembly in a viscous matrix. Next, the developed constitutive relation is employed through an analogy with incompressible elasticity to exercise the finite element technique for determination of local fiber orientation and laminate thickness after forming. Results are presented for the stretch bending of a curved beam from an arbitrary composite laminate and the bulging of a clamped sheet. Structural analyses are conducted to determine the effect of microstructure on the performance of curved beams manufactured from long discontinuous fiber composites. For the purposes of this study, several curved beams with ideal and non-ideal microstructures are compared for response under pure bending. Material parameters are determined from a separate microstructural analysis.

Pipes, R. B.↗
DOE OSTI2024

Tailoring Carbide Dispersed Steels: A Path to Increased Strength and Hydrogen Tolerance

The use of transition metal carbides is reported for use as a hydrogen trapping mechanism for ferritic and austenitic steel materials. The program combined computational modeling and simulations to guide experiments towards candidate metal carbide traps, both for interfacial and interior trapping. It was found that interfacial trapping is less effective than interior trapping, with the group IVB transition metal carbides being the most effect internal traps with a loss of carbon. The sub-stoichiometric rocksalt structure accommodate the hydrogen atoms in its octahedral interstices. Using percolation theory, carbon loss of approximately 25% or more was sufficient to ensure an interconnected network of vacancies for such trapping from the surface to the internal sites within the carbide. Using this as a guide, the program developed a means to provide a uniform dispersion of ZrC nanoparticles with either Fe or 304L micron-scale powders which was then consolidated by direct current sintering. Electrolytic hydrogen diffusivity studies confirmed the reduction of hydrogen diffusion in the matrix with increasing ZrC content, which was a linear response over the sample range studied (0.01 to 1.0 wt.%). The consolidated material was micro-tensile tested in either a non-hydrogen or hydrogen charge condition and compared to a control with no carbides. Additions up to 0.05 wt.% ZrC increased the yield strength with no loss in ductility in either the non-hydrogen or hydrogen tested condition. ZrC concentrations above this amount further increased the yield strength at the expense of ductility. While these samples had a lower absolute ductility value prior to failure, the relative change in ductility between the non-hydrogen and hydrogen charge states was less for the carbides than that of the control. Metal-rich ZrC nanoparticles were fabricated through a conformal coating process yielding ZrC0.66 particles that were then incorporated into a metal matrix. Notch fatigue testing in a hydrogen environment was conducted where the number of cycles to failure was found to be less in the control than that of the carbide addition. However, the spread in experimental data and the number of samples tested limits a conclusive outcome based on defects noticed in the gauge section of all the powder processed samples. The collective outcomes of this report provide further insight into the mechanisms by which carbides act as hydrogen traps; a means to process such carbides through powder metallurgy; and their associated mechanical performance in either a non-hydrogen or hydrogen-charged condition.

08 HYDROGEN↗
DOE OSTI2023

Grain boundary softening from stress assisted helium cavity coalescence in ultrafine-grained tungsten

The formation of helium cavities in coarse-grained materials produces hardening proportional to the number density and size of the cavities and due to the interaction of dislocations with intragranular helium defects. In nanostructured metals containing a high density of interfacial sinks, preferential cavity formation in the grain boundaries instead produces softening that is often attributed to enhanced interfacial plasticity. Here, employing two grades of ultrafine-grained tungsten, we explore this effect using targeted implantation studies to map cavity evolution as a function of the irradiation conditions and quantify its impact on the mechanical response through nanoindentation. Softening is reported at implantation temperatures above the threshold for preferential grain boundary cavity formation but at a sufficiently low fluence prior to the growth of intragranular cavities. Collective changes in the mean cavity size, density, and morphology beneath a residual impression on an implanted surface indicate that cavity coalescence accompanied the reduction in hardness. Complementary atomistic simulations demonstrate that, in tungsten grain structures exhibiting softening, grain boundary bubble coalescence is driven by stress concentrations that further act to localize strain in the grain boundaries through cooperative deformation processes involving local atomic shuffling and sliding, dislocation emission, and even the nucleation of unstable twinning events.

36 MATERIALS SCIENCE↗
DOE OSTI2020

LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report

In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Summary of Geometric Parameters and Their Effects on Performance of U-10Mo Fuel Plates

A monolithic plate-type fuel system has been under development to convert high-performance test reactors from highly enriched uranium to low-enrichment uranium fuels and is now moving into the qualification phase, a predecessor to the timely conversion of the target reactors. To qualify this fuel system, the plates must meet safety standards and perform well in a reactor. The plates must maintain mechanical integrity, exhibit geometric stability, and have stable and predictable in-reactor behavior. The requirement to maintain mechanical integrity under normal operating conditions is primarily demonstrated by successful testing. However, each high-performance reactor employs a distinct design, resulting in distinct plate geometries, with unique features, attributes, irregularities, and tolerances. Due to the abundance of such distinct geometric varieties, a single “generic” plate geometry capturing all extremes is not achievable. It is also impractical to test each of these proposed designs in a reactor. This limitation necessitates cautious evaluations since the thermomechanical response of a plate with a certain geometry may not be representative for a plate with a significantly different geometry. To address concerns related to in-reactor plate performance, large set of sensitivity studies were performed. These parametric studies aimed to better understand irradiation performance, while evaluating the sensitivity of results to various modeling inputs, including geometric, operational, and material parameters. This work studied selected geometric parameters based on provided fuel specifications and performed a series of parametric simulations. The resulting temperature, displacement, and stress strains were comparatively evaluated to determine the effects of various geometric parameters. This draft provides a “high-level summary” of our parametric sensitivity studies and summarizes the key findings.

42 ENGINEERING↗
NASA NTRSDate not supplied

Ares: A Coupling Methodology for Ablation Modeling

To enable modeling of complex and coupled ablation problems, a multi-physics framework is developed. A methodology for modeling shape change in coupled systems is presented. The approach taken to model gas-surface interactions and translate coupled surface phenomena to physically meaningful boundary conditions in the distinct solvers is discussed. Particular emphasis is placed on the nature of coupled boundary conditions pertaining to surface energy and mass balances as well as surface chemistry modeling. The developed methodology is used to simulate a shear test in arc-jet conditions in order to assess the validity of the coupled approach as well as the implementation of the relevant physical processes.

EDL↗
DOE OSTI2025

Modeling of Stress and Temperature Effects on Creep of Reduced Activation Ferritic-Martensitic Steel Alloy F82H (Tertiary Creep Modeling of RAFM Steel)

A Bayesian optimization procedure is presented for calibrating a multi-mechanism micromechanical model for creep to experimental data of F82H steel. Reduced activation ferritic martensitic (RAFM) steels based on are the most promising candidates for some fusion reactor structures. Although there are indications that RAFM steel could be viable for fusion applications at temperatures up to 600 °C, the maximum operating temperature will be determined by the creep properties of the structural material and the breeder material compatibility with the structural material. Due to the relative paucity of available creep data on F82H steel compared to other alloys such as Grade 91 steel, micromechanical models are sought for simulating creep based on relevant deformation mechanisms. As a point of departure, this work recalibrates a model form that was previously proposed for Grade 91 steel to match creep curves for F82H steel. Due to the large number of parameters (9) and cost of the nonlinear simulations, an automated approach for tuning the parameters is pursued using a recently developed Bayesian optimization for functional output (BOFO) framework [1]. Incorporating extensions such as batch sequencing and weighted experimental load cases into BOFO, a reasonably small error between experimental and simulated creep curves at two load levels is achieved in a reasonable number of iterations. Validation with an additional creep curve provides confidence in the fitted parameters obtained from the automated calibration procedure to describe the creep behavior of F82H steel at 600 °C. The model is further extended using a temperature dependent scaling law approach to simulate creep response between 550 °C and 650 °C. The efficacy of this extension is compared with the previously used scaling law approach for Grade 91 steel.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A New Look at Historic CRAC 'Divergence' Experiments - Implications for a CAAS Minimum Accident of Concern

This paper examines power history parameters and results based on historic fissile solution criticality 'divergence' experiments (referenced here as Divergence Experiments) from the 1960's-70's' Consequences Radiologiques d'un Accident de Criticite' (CRAC) campaign. The results from specific Divergence Experiments have implications for an ANSI/ANS-8.3 (1997) Criticality Accident Alarm System (CAAS) Minimum Accident of Concern (MAC). Specifically, this paper examines select CRAC Divergence Experiments in conjunction with 'qualifying reactivity additions' for a MAC defined by an absolute time constraint. A previous paper presented at the 2011 International Conference on Nuclear Criticality (ICNC 2011) illustrated the differences in 'qualifying reactivity additions' above delayed critical (and below prompt) for a variety of MAC constructs. That paper was the first to examine delayed critical kinetics in greater detail over a range of process and MAC considerations. The ICNC2011 work did not consider an absolute MAC time constraint from the excursion onset and assumed an initial excursion power. That study evaluated MAC constructs based on 'relative time differences' between reaching a specific fission rate and total fissions. The work presented here thus picks up where the ICNC2011 study left off. The primary focus of this work provides simulated power histories for specific Divergence Experiments that served as a basis for justifying a MAC noted in ANS-8.3 Appendix A. This paper documents computational modeling and dynamic multiphysics simulation with known Divergence Experiment parameters for reactivity additions at, below, and above the reported terminal reactivity. The initial step considers replicating what could be a reasonable power history for the terminal reactivity addition to the extent the same peak power and total fissions are obtained over the total time considered. Power histories for all other reactivity additions are simulated by varying the input reactivity and solution heights for the same model. A relationship is obtained between reactivity addition and the time to reach specific MAC power and fission levels within a given time constraint (e.g. 60 s). The 'qualifying reactivity addition' is determined for the total MAC fissions within that time. This value is then compared to the historic terminal value (e.g. $0.33) or any other values (e.g. 'excess reactivity of only a few cents') noted in ANS-8.3 Appendix A. In general, calculated MAC power levels have ranged from 1.8 kW to 173 kW depending on material type and form, H/U, enrichment, critical volume, and dose response function. A minimum accident consideration noted a SHEBA experiment with an excess reactivity 7 cents above delayed critical (TD ∼110 s). This excursion provided an additional 10{sup 15} fissions after an undetermined buildup time resulting in a 'few times 10{sup 15} fissions'. ICNC 2011 LEU solutions calculations show 0.2 Gy in air corresponds to 1.4 x10{sup 16} fissions. Thus, while historically, 0.2 Gy has been loosely associated with 'several times 10{sup 15} fissions', in reality, the actual total fissions could be an order of magnitude larger, and the power associated with 0.2 Gy/min (7.5 kW) larger as well. The reactivity addition required to obtain the actual 0.2 Gy total fissions over 60 s could be much larger than considered. Another reason the qualifying reactivity level could be larger than previously noted, involves temperature reactivity feedback. Compensated reactivity results in a departure from a pure exponential power rise, increasing the doubling time (TD) as the excursion approaches peak power. Uncompensated reactivity additions and constant TD have been traditionally assumed for a pure exponential power rise through the kW range to determine either total MAC fissions, power, or arrival time to these values.

42 ENGINEERING↗
NASA NTRS1996

The Response of the BATSE LADs to Radiation From the Crab Nebula and Plans for Radioactivity Studies on Space Station

The Burst And Transient Source Experiment (BATSE) onboard the Compton Gamma-Ray Observatory (CGRO) was designed to measure X-rays and gamma rays with energies from about 50 keV to above 2 MeV. As with many scientific investigations, the success of the original experiment lead to additional areas of research interest. In the case of BATSE the ability to observe the radiation from sources down to about 20 keV became readily apparent. This lead to a continuing program of measuring the spectrum of radiation from stellar objects at these lower energies. One of these, the Crab Nebula, has a very steady radiation flux and, thus, has become a "standard candle" for such measurements. The Large Area Detectors (LADS) on BATSE contain a 1.27-cm thick, 25.4-cm radius NaI(Tl) detector behind a 6.35-mm thick polystyrene Charged Particle Detector (CPD) used to "veto" charged particles signals. The detectors have been calibrated with a series of gamma and X-ray sources and the results carefully simulated with a Monte Carlo code. In the calibration process the computer simulation accounts for scattering from material in the counting room as well as the BATSE structure. For an orbiting detector, scattering from the entire spacecraft must be modeled as well as for all covering material over the detectors. Five years after CGRO was launched on April 5, 1991, a large body of observational data has been taken of the Crab Nebula. The technique used for these observations, and for many other X-ray sources, is Earth occultation. From the perspective of the spacecraft, the Earth occults most stellar objects once in orbit, i.e., the signal is lost as the source sets and is regained as the source rises. A careful analysis of the continuing signals from all sources measured allows for an accurate measurement of the spectrum of a given source. An analysis of this data from the Crab has indicated that the LADs are very responsive at energies as low as 20 keV-at energies below the range of calibration. While the model accounts for many of the interactions of the photons with the detectors, the observation of nonstatistical deviations at low energy and at small angles has suggested a need to recalibrate at energies where the attenuation effects are increasing exponentially.

Laird, C. E.↗
DOE OSTI2024

Mechanism of Antiferroelectricity in Polycrystalline ZrO 2

The size and electric field dependent induction of polarization in antiferroelectric ZrO 2 is the key to several technological applications that are unimaginable a decade ago. However, the lack of a deeper understanding of the mechanism hinders progress. Molecular dynamics simulations of polycrystalline ZrO 2 , based on machine-learned interatomic forces with near ab initio quality, shed light on the fundamental mechanism of the size effect on the transition fields. Stress in the oxygen sublattice is the most important factor. The so constructed interatomic forces allow the calculation of the transition fields as a function of the ZrO 2 film thickness and predict the ferroelectricity at large thickness. The simulation results are validated with electrical and piezo response force microscopy measurements. The results allow a clear interpretation of the properties of the double-hysteresis loops as well as the construction of the free energy landscape of ZrO 2 grains.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A data driven approach for cross-slip modelling in continuum dislocation dynamics

Cross-slip is a thermally activated process by which screw dislocation changes its glide plane to another slip plane sharing the same Burgers vector. The rate at which this process happens is determined by a Boltzmann type expression that is a function of the screw segment length and the stress acting on the dislocation. In continuum dislocation dynamics (CDD), the information regarding the length of the screw dislocation segment and local stress state on dislocations are lost due to the coarse-grained representation of the density. Here, in this work, a data driven approach to characterize the lost information by analyzing the discrete dislocation configurations is proposed to enable cross-slip modeling in the CDD framework in terms of the coarse-grained dislocation density and stress fields. The analysis showed that the screw segment length follows an exponential distribution, and the stress fluctuations, defined as the difference between the stress on the dislocations and the mean field stress in CDD, follows a Lorentzian distribution. A novel approach for cross slip implementation in CDD employing the screw segment length and stress fluctuation statistics was proposed and rigorously tested by comparing the CDD cross-slip rates with discrete dislocation dynamics (DDD) rates. This approach has been applied in conjunction with three cross-slip models used in DDD simulations differing mainly in the functional form of cross slip activation energy. It was found that different cross-slip activation energy formulations yielded different cross-slip rates, yet the effect on mechanical stress-strain response and dislocation density evolution was minimal for the [001] type loading.

42 ENGINEERING↗
DOE OSTI2022

Investigating laser ablated plume dynamics of carbon and aluminum targets

Recently acquired high-resolution images of nanosecond laser ablation plumes suggest a strong correlation between the internal plume structure and the type of material being ablated. However, the details of this relation are currently not well understood. In this work, we attempt to explore this correlation using a 2D radiation hydrodynamics model to study the dependence of internal plume structure formation on the ablation material. Here spatio-temporal emission maps and plume expansion velocities from experimental measurements are compared with the model predictions, including synthetic emission maps. The shape and expansion rate of an outer air plume region are found to be in good agreement for both carbon and aluminum, as are the inner material plume dynamics for carbon ablation. The largest disagreement is observed in the case of a polished aluminum target, where the chaotic inner plume features seen in the experimental images are not observed in the model. The possible physical mechanisms responsible for this discrepancy are discussed. This effort constitutes a continued development toward a predictive model of ablation plume dynamics and chemistry for various materials in extreme environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗