Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesoscale”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the correlation between plastic strain and misorientation in polycrystalline body-centered-cubic microstructures with an emphasis on the grain size, loading history, and crystallographic orientation

In this work, the correlation between plastic strain and crystallographic misorientation, grain size, grain orientation, distance from grain boundary, and loading history were investigated experimentally and numerically for body-centered-cubic (BCC) polycrystalline microstructures based on a misorientation deviation (MD) approach. Nine monotonic tensile experiments were performed on two BCC titanium alloys inside a scanning electron microscope (SEM). The influence of reference orientation was explored both at the grain scale and at the mesoscale using electron backscattered diffraction (EBSD). The correlation between global plastic strain and the MD was quantified. The tendency for orientation change was quantified as a function of plastic strain and grain orientation for three crystallographic orientations (i.e., [100], [110], and [111]) with respect to tensile direction. The subpopulation of small grains exhibited a lower level of misorientation dispersion compared with larger grains, and this discrepancy became more pronounced at higher strains. An empirical equation was proposed to estimate the level of misorientation dispersion for individual grains as a function of grain size and global plastic strain level. Two interrupted in-situ SEM experiments were performed, and this resulted in a significantly increased misorientation compared with uninterrupted tests performed to similar plastic strain levels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Implementation and experimental validation of nonlocal damage in a large-strain elasto-viscoplastic FFT-based framework for predicting ductile fracture in 3D polycrystalline materials

Ductile materials, such as metal alloys, can undergo substantial deformation before failure. Additionally, these materials are usually of polycrystalline composition and exhibit strongly anisotropic behavior at small length scales. Previously developed fast Fourier transform (FFT)-based models can model ductile fracture of isotropic materials or the elastic–plastic behavior of anisotropic polycrystalline materials; however, there remains a need to couple both capabilities. This work extends a large-strain FFT-based crystal plasticity model to simulate ductile fracture of polycrystalline materials. In this work, a triaxiality-based continuum damage mechanics (CDM) formulation is incorporated into a large-strain elasto-viscoplastic FFT (LS-EVPFFT) framework. The CDM formulation is augmented with an integral-based nonlocal regularization approach that correctly handles gas-phase material necessary to model unconstrained surfaces. To validate the damage-enabled LS-EVPFFT framework, mesoscale copper tensile coupons were machined using microwire electrical discharge machining and experimentally characterized using electron backscatter diffraction. In-situ optical digital image correlation was performed during uniaxial testing to provide a side-by-side comparison of the experimental and computational strain fields and stress–strain responses. The damage-enabled LS-EVPFFT framework can simulate the complete macroscopic stress–strain response of ductile polycrystals to failure. The model reproduces necking behavior that qualitatively agrees with experimental observations. By leveraging the relatively low computational cost of the damage-enabled LS-EVPFFT framework, the framework presented here allows the ductile fracture response of 3D polycrystalline materials to be tractably predicted.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interaction of extended dislocations with nanovoid clusters

Voids of nanoscale dimensions in irradiated metals can act as obstacles to dislocation motion and cause strengthening. In this work, nanovoid strengthening and the influences of void size, void spacing and material properties, such as stacking fault energies, on dislocation bypass mechanisms are investigated using Phase Field Dislocation Dynamics, a three-dimensional mesoscale model that predicts the minimum energy pathway taken by discrete dislocations. A broad range of face centered cubic metals (copper, nickel, silver, rhodium, and platinum) and nanovoid sizes and spacings are treated, altogether spanning void size–to–dislocation stacking fault width ratios from less than unity to ten. Material γ-surfaces, calculated from ab initio methods, are input directly into the formulation. The analysis reveals that the critical bypass stress scales linearly with the linear void fraction, effective isotropic shear modulus, and ratio of the intrinsic to unstable stacking fault energies. With only a few exceptions, the critical stress is controlled by the stress required for the leading partial to impinge the voids (to move within range of the attractive image stress field of the void). When the void diameter is nearly an order of magnitude greater than the stacking fault width, the mechanism determining critical strength shifts to the stress for the dislocation to breakaway after partially cutting the void. Furthermore, this situation corresponds to that treated by line tension models and is realized here for Pt, with a sub-nanometer stacking fault width.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Using in-situ strain measurements to evaluate the accuracy of stress estimation procedures from fracture injection/shut-in tests

Fracture injection/shut-in tests are commonly used to measure the state of stress. Injection creates a hydraulic fracture (or in some cases, opens a preexisting fracture), and then the pressure after shut-in is monitored to identify fracture closure. Different interpretation procedures have been proposed for estimating closure, and the procedures sometimes yield significantly different results. In this study, direct, in-situ strain measurements are used to observe fracture reopening and closure. The tests were performed as part of the EGS Collab project, a mesoscale project performed at 1.25 and 1.5 km depth at the Sanford Underground Research Facility. The tests were instrumented with the SIMFIP tool, a double-packer probe with a high-resolution three-dimensional borehole displacement sensor. The measurements provide a direct observation of the fracture closure signature, enabling a high-fidelity estimate of the fracture closure stress (ie, the normal stress on the fracture). In two of the four tests, injection created an opening mode fracture, and so the closure stress can be interpreted as the minimum principal stress. In the other two tests, injection probably opened preexisting natural fractures, and so the closure stress can be interpreted as the normal stress on the fractures. The strain measurements are compared against different proposed methods for estimating closure stress from pressure transients. The shut-in transients are analyzed with two techniques that are widely used in the field of petroleum engineering – the ‘tangent’ method and the ‘compliance’ method. In three of the four tests, the tangent method significantly underestimates the closure stress. The compliance method is reasonably accurate in all four tests. Closure stress is also interpreted using two other commonly-used methods – ‘first deviation from linearity’ and the method of Hayashi and Haimson. In comparison with the SIMFIP data, these methods tend to overestimate the closure stress, evidently because they identify closure from early-time transient effects, such as near-wellbore tortuosity. In two of the tests, microseismic imaging provides an independent estimate of the size of the fracture created by injection. When combined with a simple mass balance calculation, the SIMFIP stress measurements yield predictions of fracture size that are reasonably consistent with the estimates from microseismic. In conclusion, the calculations imply an apparent fracture toughness 2-3x higher than typical laboratory-derived values.

58 GEOSCIENCES↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗
DOE OSTI2022

A mass conservative, well balanced, tangency preserving and energy decaying method for the shallow water equations on a sphere

Here a fully discrete surface finite element method is proposed for solving the viscous shallow water equations in a bounded Lipschitz domain on the sphere based on a general triangular mesh. The method consists of a modified Crank–Nicolson method in time and a Galerkin surface finite element method in space for the fluid thickness H and the fluid velocity u. A finite element space tangential to the sphere at all finite element nodes is proposed to approximate the fluid velocity u. The proposed method has second-order accuracy in time and first-order accuracy in space, and preserves mass conservation, well balancedness, tangency of velocity to the sphere, and energy decay. Numerical experiments are presented to illustrate the accuracy of the proposed method and the preservation of the physical properties, including mass conservation, well balancedness, and energy decay. A numerical simulation of ocean mesoscale activity on a circular basin with a continental shelf is provided.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Accelerating discrete dislocation dynamics simulations with graph neural networks

Discrete dislocation dynamics (DDD) is a widely employed computational method to study plasticity at the mesoscale that connects the motion of dislocation lines to the macroscopic response of crystalline materials. However, the computational cost of DDD simulations remains a bottleneck that limits its range of applicability. Here, we introduce a new DDD-GNN framework in which the expensive time-integration of dislocation motion is entirely substituted by a graph neural network (GNN) model trained on DDD trajectories. As a first application, we demonstrate the feasibility and potential of our method on a simple yet relevant model of a dislocation line gliding through an array of obstacles. We show that the DDD-GNN model is stable and reproduces very well unseen ground-truth DDD simulation responses for a range of straining rates and obstacle densities, without the need to explicitly compute nodal forces or dislocation mobilities during time-integration. Our approach opens new promising avenues to accelerate DDD simulations and to incorporate more complex dislocation motion behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Advanced magnetic X-ray spectro-microscopies to characterize mesoscopic magnetic materials

Here, this article provides a brief overview of advanced magnetic X-ray spectro-microscopies that are widely used in characterizing mesoscopic magnetic materials. Common to those techniques are various X-ray magnetic dichroism effects that are used as magnetic contrast mechanism. A particular interest is to use those techniques to image the underlying microscopic spin structures in magnetic materials with high spatial resolution and to ultimately resolve their full 3D mesoscale characteristics with elemental sensitivity, as well as their ultrafast dynamics upon excitations with field and current pulses. Recent research examples using various magnetic X-ray spectro-microscopies are presented to showcase their specific features with a focus on imaging novel topological spin textures, such as vortices, skyrmions and Hopfions, which are considered as potential building blocks towards low power, high-speed and high-density magnetic devices that could transform information technologies and potentially be used in biomedical applications. Future developments with magnetic X-ray spectro-microscopies harness the full coherence of next generation X-ray sources and could open the path towards single shot imaging with spatial and temporal resolutions down to fundamental magnetic length and time scales.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stress effects on the energy barrier and mechanisms of cross-slip in FCC nickel

The energy barrier for homogeneous cross-slip of a screw dislocation in face-centered cubic (FCC) nickel is calculated using atomistic simulations as a function of the Escaig stress on the glide plane, the Escaig stress on the cross-slip plane, and the Schmid stress on the cross-slip plane. Two cross-slip mechanisms, Friedel-Escaig (FE) and Fleischer, are examined and their energy barriers are calculated for a large number of stress combinations. For each mechanism, the energy barrier as a function of three stress components can be reduced into a one-dimensional function of an effective stress. The stress domains in which FE and Fleischer mechanisms operate respectively are determined. Additionally, the FE mechanism dominates when the Escaig stress on the glide plane (in the direction that reduces the stacking fault width) is the largest stress component. Increasing the Schmid stress and Escaig stress (in the direction that expands the stacking fault width) on the cross slip plane promotes the Fleischer mechanism. The cross slip energy barrier functions obtained here can be used as input functions for computing cross slip rates in mesoscale dislocation dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations

Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Twin nucleation and growth in hexagonal close-packed metals: The role of slip-mediated plasticity on twin embryo formation and evolution

Twinning is a key deformation mechanism in hexagonal close-packed (hcp) metals, which are typified by a lack of easily activated slip systems that can accommodate a general state of loading. Pragmatically, the nucleation and evolution of twin domains occurs concomitantly with slip, such that the eventual twin network is conditioned by both external and internal stresses resulting, among others, from the evolution of dislocations. However our understanding of the interplay between dislocations, and twin nucleation and stability remains limited. This work focuses on elucidating the influence of dislocation-mediated plasticity on the formation, growth, and stability of twin embryos. First, a new mesoscale spectral crystal plasticity-twinning framework is developed and used to quantify the change in the free energy landscape following the nucleation {$10\bar{1}2$} of twins in Mg. Representative twin morphologies are modeled under conditions of limited and profuse slip activity, which are emulative of small- and bulk-scale samples, respectively. Then, the driving traction profiles around twin embryos are investigated via a sharp interface approach to obtain insights into how concurrent slip-mediated plasticity can influence the growth/stabilization of nanometric twins. The driving traction profiles are further utilized to determine the stability of twin embryos post loading. The initial dislocation density in the samples, and within the different domains (i.e., twin vs. parent), is seen to have a significant effect on the twin nucleation stress. Namely, the activation of high levels of concomitant plasticity in the parent grain is seen to significantly drive the formation of nanometric twin nuclei at stresses as low as 250 MPa. Further, the sharp interface analysis reveals that profuse plasticity in the parent grain simultaneously alters the forward and back stresses, such that the magnitude/polarity of the driving tractions become increasingly favorable for nanometric twin growth when slip is active. Finally, prior plasticity in the parent grain is seen to result in favorable driving tractions for nanometric twin growth, even at applied stresses as low as ~ 100 MPa. In conclusion, these results are in stark contrast to a case without any dislocations, wherein applied stresses as high as ~ 500 MPa are necessary to grow the twin domains.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the role of geometrically necessary dislocations in void formation and growth in response to shock loading conditions in wrought and additively manufactured Ta

This study investigates the role of geometrically necessary dislocations (GNDs) and microstructure on void nucleation and growth in wrought and additively manufactured (AM) tantalum subjected to high-strain rate loading. Multi-modal 3D data was collected using TriBeam tomography to calculate GND densities and their spatial relationship to voids. A microstructural comparison between the wrought and AM samples identified distinct void shapes and locations, with intragranular voids and more spherical voids frequently observed in the AM dataset. Results indicate that voids preferentially form at both high-angle grain boundaries and low-angle subgrain boundaries, the latter of which are frequently observed in the AM material. Through a radial distribution analysis of all voids in the datasets, significant GND localization to near-void-surface regions was observed in both samples. 3D crystal plasticity simulations were employed to extend the experimental observations, revealing higher void growth rates in [111] oriented grains when compared to [001] grains. The simulations also suggest that GNDs can be generated as part of the void growth process, with more GND accumulation for growth in a [111] grain than a [001] grain. These findings provide valuable insights into the links between nanoscale void nucleation, mesoscale void growth, and microstructural effects in dynamically loaded tantalum.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A improved equation of state for Xe gas bubbles in gamma U-Mo fuels

A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U-Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in ?U-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.

Beeler, Benjamin↗
DOE OSTI2020

Grand potential sintering simulations of doped UO 2 accident-tolerant fuel concepts

Chromium-doped UO is a widely-studied near-term deployable accident-tolerant fuel concept because it results in a dense, large-grain structure that increases the fuel resistance to densification, swelling, and fission gas release. A new charged-interstitial mechanism was recently proposed to describe the behavior of dopants like chromium in sintered UO 2 . Based on that mechanism, manganese was suggested as an even stronger dopant than chromium. In the current work we use mesoscale sintering simulations in an effort to validate the new mechanism. We compare the relative behavior of Cr-doped and undoped UO 2 against experimental data in the literature. We also make predictions of the relative behavior of Mn-doped UO 2 . Simulations are done using the phase field-based grand potential sintering model. Dopants have two effects on sintered UO 2 . They marginally increase the densification rate and greatly increase the average grain size. Both of these effects are individually tested in small-scale simulations. Then large 3D sintering simulations are used to test the combined behavior of both effects. The results for the densification rate simulations are consistent with experiments for Cr-doped fuels. However, we note the grain growth rates are lower than what is found in the literature for Cr-doped fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Radiation driven diffusion in γU-Mo

A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). A critical phenomenon of interest with U-Mo monolithic fuel is the large amount of swelling that takes place during operation, particularly at high fission densities. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic and kinetic properties into mesoscale and continuum level fuel performance modeling codes. One such property where there exists incomplete data is the diffusion of relevant species under irradiation. Fuel performance swelling predictions rely on an accurate representation of diffusion in order to determine the rate of fission gas swelling and the local microstructural evolution. In this work, we present molecular dynamics simulations of the radiation driven diffusion of U, Mo and Xe in U-Mo nuclear fuels. Diffusion coefficients for each species are determined over a range of temperatures and compositions. In this work, updated diffusion coefficients are presented that are applicable under irradiation that incorporate both intrinsic and radiation driven diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Investigation of δ zirconium hydride morphology in a single crystal using quantitative phase field simulations supported by experiments

In light water nuclear reactors, waterside corrosion of the cladding material leads to the production of hydrogen, a fraction of which is picked up by the zirconium cladding and precipitates into brittle hydride particles. These nanoscale hydride particles aggregate into mesoscale hydride clusters. The principal stacking direction of the nanoscale hydrides precipitated in the cladding tube changes from circumferential in the absence of applied stress to radial under circumferential applied stress. A quantitative phase field model has been developed to predict the hydride morphology observed experimentally and identify the mechanisms responsible for nanoscale hydride stacking. The model focuses on nanoscale hydride precipitation in a single zirconium grain with a detailed description of the anisotropic elastic contribution. The model predictions concerning the shape, orientation, and stacking behavior of nanoscale hydride are analyzed and compared with experimental observations. The model accurately accounts for the experimentally observed elongated nanoscale hydride shape and the stacking of hydrides along the basal plane of the hexagonal zirconium matrix. When investigating the role of applied stress in hydride morphology, the model challenges some of the mechanisms previously proposed to explain hydride reorientation. Finally, although hydride reorientation has been hypothesized to be caused by a change in nanoscale hydride shape, the current study shows that these mechanisms are unlikely to occur.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

The mobility of small, over-pressurized helium bubbles in tungsten at 2000 K

Fusion reactor environments inevitably lead to the formation of high-pressure helium bubbles whose nucleation, growth, and diffusion strongly impact the performance of plasma-facing components. Here this research describes a diffusion mechanism of over-pressurized bubbles via a sequence of Frenkel pair nucleation, self-interstitial migration, and Frenkel pair annihilation. Molecular dynamics was used to simulate the diffusion of small bubbles in tungsten at 2000 K with helium-per-vacancy ratios in the range of 4.5 to 7. The diffusion coefficients are calculated and their dependence on helium content, number of vacancies, and number of attached self-interstitials is characterized. It is found that bubbles are most mobile when the nucleation/annihilation rates of Frenkel pairs are nearly equal and when the bubbles nucleate and annihilate a single self-interstitial. All bubbles experience a peak diffusivity, which can be as high as 10 -11 m 2 /s decreasing with bubble size. The calculated diffusion coefficients provide valuable insight into the mobility of small, high-pressure bubbles, and can be used as input parameters in mesoscale models to improve predictions of plasma-surface interactions. (LA-UR-21-21881)

36 MATERIALS SCIENCE↗