Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “slip”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2024

Quantification of grain boundary effects on the geometrically necessary dislocation density evolution and strain hardening of polycrystalline Mg 4Al using in situ tensile testing in scanning electron microscope and HR-EBSD

In situ tensile testing in a scanning electron microscope (SEM) in conjunction with high-resolution electron backscatter diffraction (HR-EBSD) under load was used to characterize the evolution of geometrically necessary dislocation (GND) densities at individual grain boundaries as a function of applied strain in a polycrystalline Mg—4Al alloy. The increase in GND density was investigated at plastic strains of 0 %, 0.6 %, 2.2 %, 3.3 % from the area including 76 grains and correlated with (i) geometric compatibility between slip systems across grain boundaries, and (ii) plastic incompatibility. We develop expressions for the grain boundary GND density evolution as a function of plastic strain and plastic incompatibility, from which uniaxial tensile stress-strain response of polycrystalline Mg—4Al are computed and compared with experimental measurement. The findings in this study contribute to understanding the mechanisms governing the strain hardening response of single-phase polycrystalline alloys and more reliable prediction of mechanical behaviors in diverse microstructures.

Geometrically necessary dislocations↗
DOE OSTI2022

Crystal plasticity including a phase-field deformation twinning model for the high-rate deformation of cyclotetramethylene tetranitramine

A new finite strain thermomechanical model for the high-rate deformation of the β-polymorph of cyclotetramethylene tetranitramine (β-HMX) has been developed and applied to simulations of plate impact experiments. The crystal plasticity model is based on a model developed previously for RDX (Luscher et al. 2017), which is extended to incorporate deformation twinning. Twinning during normal plate impacts is simulated with a phase-field twin model. First, material parameters governing the kinetics of dislocation slip are calibrated on the subset of simulations which had negative Schmid factors for the twin system. Second, a parametric study of the twin material parameters was performed to find suitable values. Overall, the results of the simulations with the phase-field twinning model are reported for impacts on several crystal orientations. We find that the twin growth decreases with increasing distance from the impact surface because of dissipation of the shock front via dislocation-mediated plasticity, and that the simulated interface velocity with and without phase-field twinning do not show appreciable differences. These modeling results suggest that the significance of twinning in β-HMX cannot be determined with traditional loading configurations and diagnostics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Non-local large-strain FFT-based formulation and its application to interface-dominated plasticity of nano-metallic laminates

Here, this paper presents a novel formulation and its robust numerical implementation of strain-gradient (SG) crystal plasticity within a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. The resulting non-local SG-LS-EVPFFT formulation is used to model and understand the process of kink band formation during layer-parallel compression of nano-metallic laminates (NMLs). NMLs are layered composites with nanoscale thicknesses, thus requiring consideration of the interaction between dislocations and interfaces within the micromechanical model. The length-scale parameter of the SG model is calibrated by simulating a double pile-up and comparing predictions to analytical solution. This required new expressions for the defect energy, resulting in more accurate double pile-up predictions. The calibrated SG-LS-EVPFFT model is then used to simulate layer-parallel compression of copper-niobium NML. Formation of kink bands is predicted, and the model is used to rationalize the microscopic mechanisms enabling the formation process. It is found that accumulation of dislocations at interfaces leads to activation of layer-parallel slip, which in turn leads to kink band formation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of Zn and Ca on the deformation texture evolution of Mg–Zn–Ca alloys at elevated temperatures

This work utilized crystal plasticity simulations and experimental thermo-mechanical processing to examine how Zn and Ca content affects texture evolution in Mg–Zn–Ca alloys during simulated hot rolling. Four different compositions were studied: unalloyed Mg and three Mg alloys of ZX0p50 (Mg- 0.5 wt% Zn- 0.1 wt% Ca), ZX30 (Mg-3 wt% Zn- 0.1 wt% Ca), and ZX31 (Mg-3 wt% Zn- 0.3 wt% Ca). Multi-pass Gleeble experiments simulated hot rolling using plane strain conditions. Unalloyed Mg developed a strong basal texture, while the addition of Zn and Ca significantly weakened the texture in ZX30 and ZX31; this reduction was not seen in ZX0p50. Notably, ZX31 exhibited a split basal texture aligned with the rolling direction. Crystal plasticity simulations revealed that the weaker basal texture in ZX30 and ZX31 resulted from reduced activation of basal and twinning modes compared to unalloyed Mg and Zx0p50. Increased Zn and Ca content raised the ratios of basal to pyramidal critical resolved sheer stress (CRSS) and twin to pyramidal CRSS, enhancing pyramidal slip activity and weakening the basal texture. Extension twinning played a critical role in texture development during the first deformation pass. In unalloyed Mg and ZX0p50, twinning led to a strong basal texture, while ZX30 and ZX31 had weak basal textures that increased only slightly with further deformation. The split basal poles simulated in ZX31 were consistent with the experimental findings, highlighting the interplay between alloy composition and texture evolution.

Crystal plasticity finite element↗
DOE OSTI2022

Stress-controlled fatigue of HfNbTaTiZr high-entropy alloy and associated deformation and fracture mechanisms

The stress-controlled fatigue tests are carried out at a stress ratio of 0.1 and a frequency of 10 Hz, and span both low-cycle and high-cycle regimes by varying the applied stress amplitudes. The high-cycle fatigue regime gives a fatigue strength of 497 MPa and a fatigue ratio of 0.44. At equivalent conditions, the alloy's fatigue strength is greater than all other high-entropy alloys (HEAs) with reported high-cycle fatigue data, dilute body-centered cubic alloys, and many structural alloys such as steels, titanium alloys, and aluminum alloys. Through in-depth analyses of crack-propagation trajectories, fracture-surface morphologies and deformation plasticity by means of various microstructural analysis techniques and theoretical frameworks, the alloy's remarkable fatigue resistance is attributed to delayed crack initiation in the high-cycle regime, which is achieved by retarding the formation of localized persistent slip bands, and its good resistance to crack propagation in the low-cycle regime, which is accomplished by intrinsic toughening backed up by extrinsic toughening. Additionally, the stochastic nature of the fatigue data is neatly captured with a 2-parameter Weibull model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Shear band velocity and activation volume during shear deformation by acoustic emission in a Zr-based bulk metallic glass

Here, an analysis is presented for the shear band velocity of plastic deformation and activation volume of the shear transformation zone (STZ) in a Zr 60 Cu 30 Al 10 bulk metallic glass (BMG). Compared to the values calculated by the stick-slip model, the values of shear band velocity calculated after optimization by the acoustic emission technique are much larger than those reported previously. Plastic flow in BMGs is found to be described by a power-law correlation between the maximum velocity of the shear band and the strain rates. Combining the free volume theory and the STZ theory, the relationship between the volume change in the STZ and the shear band velocity is derived as $Ω=A$(27–ln$v/l$). The value of A, which is a material dependent constant, is 5.43 × 10 –3 for the current BMG at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Global and local texture development during initial plastic deformation of cold-pilgered Zircaloy-4 tubing

In this work we present a thorough characterization of the texture development during the initial step of a cold pilgering process of a Zircaloy-4 TREX tube, on which the axial cross section is reduced by 80% over a 40cm length with an overall Q-ratio of 2.34, followed by a recrystallization thermal treatment. Texture across a tapered specimen representing this working zone was characterized by neutron diffraction and high energy synchrotron x-ray diffraction at 28 positions (4 locations along its length and 7 locations across its thickness with 0.3 mm spatial resolution), and on the deformed end after the thermal treatment. Results show that the starting material (TREX tube) presents a texture gradient through its wall thickness, with the c poles aligned to the hoop direction at the outer surface (with radial and hoop Kearns factors: f R =0.36, f H = 0.50), and gradually rotating around the axial direction to become a quasi-radial orientation (~35 o out of the radial direction) at the inner surface ( f R =0.43, f H =0.35). The overall texture becomes more radial due to the cold pilgering, but the texture gradient remains during the whole deformation process, being always more radial at the inner surface, ( f R =0.53, f H =0.37 at 80%) than at the outer surface ( f R =0.41, f H =0.50at 80%). The most remarkable changes in texture produced by the rolling process is a 30 o rotation of the grains around their c-axis, producing a marked increase of the <11 2 ¯ 0> // axial fiber texture at expenses of the original <10 1 ¯ 0> // axial fiber texture observed in the TREX tube. The observed changes in texture during deformation were correctly explained by visco-plastic self-consistent (VPSC) modelling. This was achieved after the inclusion of the basal <a> slip system, which plays a central role to explain the 30 o rotation around the c-axis during deformation. The texture evolution was correctly described after the inclusion of a drag effect to grain rotation due to the interaction with neighbouring grains. The through thickness texture differences are greatly reduced by the recrystallization thermal treatment, with similar Kearns factors at the inner surface ( f R =0.50, f H =0.28) and the outer surface ( f R =0.48, f H =0.39). The predominance of the <10 1 ¯ 0> // axial fiber texture is recovered after recrystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deformation localisation in ion-irradiated Fe and Fe10Cr

Determining the mechanisms for irradiation-induced ductility loss is crucial for the design of reactor structural components. Here, the deformation characteristics around nanoindents in Fe and Fe10Cr irradiated with Fe ions to ~1 displacement-per-atom at 313 K are non-destructively studied. Slip steps surrounding the nanoindents indicate that deformation is localised in the irradiated materials. Lattice rotation and strain fields near the indent site show over 87% confinement of plasticity in the irradiated material. Cr has little effect on the irradiation-induced changes in pile-up topography and deformation fields, suggesting it has limited impact on retaining strain hardening capacity and reducing irradiation-induced embrittlement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Assessing the viability of a new subsize tensile specimen geometry for evaluation of structural nuclear and additively manufactured materials

The use of subsize specimens in nuclear materials testing has been a subject of ongoing interest due to radiation safety concerns and resource conservation needs. It is also of interest in additive manufacturing (AM), also known as 3D-printing, as subsize specimens can more accurately represent the behavior of the small, intricate geometries often produced using AM. A novel, extremely small geometry, called the Subsize Teeny (SST) was recently developed and is of interest for implementation. Given its extraordinarily small size and the complexities associated with subsize specimen testing, adequate vetting of this geometry is necessary to ensure data quality. Here, a variety of unirradiated nuclear structural materials were tested in the SST geometry and compared against the well-established SSJ3 geometry. In addition, two case studies implementing the SST as a screening geometry for AM materials were also conducted. The question of SST viability was found to be highly nuanced and will often be dependent on the context or application in question. It was determined, however, that the SST is a largely invalid geometry for exceptionally coarse-grained materials or in cases where the physical defect volume equals or exceeds 0.1 % or where the specimen machining parameters result in significant surface alterations. On the other hand, it was determined that the SST may be employed with confidence if the test material is nearly or totally free of physical defects, isotropic, demonstrates homogeneous plastic deformation, and possesses a fine-grained, nearly or totally homogeneous microstructure with at least twelve slip systems.

Additive manufacturing↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
DOE OSTI2021

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗
DOE OSTI2021

Grain-resolved temperature-dependent anisotropy in hexagonal Ti-7Al revealed by synchrotron X-ray diffraction

Hexagonal metals have anisotropic coefficients of thermal expansion causing grain-level internal stresses during heating. High energy x-ray diffraction microscopy, a non-destructive, in situ, micromechanical and microstructural characterization technique, has been used to determine the anisotropic coefficients of thermal expansion (CTEs) for Ti-7Al. Two samples of polycrystalline α-phase Ti-7Al were continuously heated from room temperature to 850 °C while far-field HEDM scans were collected. The results showed a change in the ratio of the CTEs in the ‘a’ and ‘c’ directions which explains discrepancies found in the literature. The CTE additionally appears to be affected by the dissolution of α 2 precipitates. Analysis of the grain-resolved micromechanical data also shows reconfiguration of the grain scale stresses likely due to anisotropic expansion driving crystallographic slip.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Strengthening mechanisms, hardening/softening behavior, and microstructure evolution in an LPSO magnesium alloy at elevated temperatures

Here, rapid degradation of the mechanical properties of conventional Mg alloys with temperature precludes their wide application in industry. The LPSO (long period stacking ordered) Mg alloys, on the other hand, show considerable potential to achieve excellent high-temperature mechanical performance. In this study, a Mg-5.1Y-2.4Zn (wt.%) extrusion alloy with the microstructure comprising of α-Mg matrix and LPSO phase was prepared. It is shown that the yield strength of the alloy at 200 °C retained ~83% of the yield strength of 198 MPa at room temperature (RT). Importantly, the ductility also improved, owing to the increased non-basal slip activities at elevated temperatures. In contrast to RT and 200 °C, the alloy exhibited softening tensile behavior with an apparent loss of yield strength at 300 °C. In-situ lattice strain analysis revealed that load transfer from α-Mg to LPSO remains effective at temperatures up to 200 °C, governing the high yield strength of the alloy. The post-mortem analysis by means of electron microscopy confirmed that the LPSO phase is structurally stable at 200 °C. In particular, the high strain hardenability is rationalized by the sandwiched LPSO structures together with stacking faults serving as effective obstacles to the motion of non-basal dislocations. The load transfer effect disappears at the temperature of 300 °C, resulting in the fast deterioration of strength of the alloy. The cracking of the LPSO phase along with the recrystallizing of α-Mg grains were found to be the main causes for the softening behavior of the alloy at 300 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗
DOE OSTI2021

Ultrasonic additive manufacturing of zirconium: Pilot results

Ultrasonic additive manufacturing (UAM) was successfully applied to the Zirconium material system to create a three-dimensional component prototype. The UAM process resulted in grain size refinement and significant localized deformation. Electron backscatter diffraction analysis revealed activated slip systems only applicable at high deformation temperatures, while texture analysis showed a decrease in measurable texture of the UAM build in comparison to the initial Zr-foil. Average tensile strengths in the X (sonotrode travel), Y (vibration), and Z (build) directions averaged 435 MPa, 458 MPa, and 359 MPa, respectively, demonstrating the viability of UAM for Zr-based materials. Furthermore, the delamination of some Z specimens along foil boundaries during loading suggests an interplay between interfacial Ti impurities introduced during welding and the spatial dependence of weld quality for this material system.

36 MATERIALS SCIENCE↗