Dual-Domain Material Point method with Local Stress Difference for Fragmentation
Explore the source record for details and available documents.
36 MATERIALS SCIENCE↗
SEARCH · Search NASA
Results for “local stress difference”
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.
Micro-strains, local stresses, and coherently diffracting domain size in shock compressed Al(100) single crystals
In situ x-ray diffraction (XRD) measurements and their analysis in Al single crystals shock compressed along the [100]-direction were utilized to examine shock wave induced microstructural heterogeneities. High-resolution XRD line profiles for the 200, 400, and 600 Al peaks were measured in uniaxial strain compression states to either 5.6 or 11.7 GPa and partial stress release to 3.5 or 6.6 GPa, respectively. Broadening of the XRD line profiles was analyzed to determine the magnitude of the longitudinal micro-strain distribution (0.195% and 0.28% full width at half maximum for 3.5 and 6.6 GPa stresses, respectively) and the size of coherently diffracting domains (CDDs) (0.125 and 0.07 μm for 3.5 and 6.6 GPa stresses, respectively). From the longitudinal micro-strain distributions, the distribution of local stress differences (or stress deviators) was obtained in the shocked state. The full width at half maximum of this distribution, a measure of the local stress inhomogeneities, is greater than half of the macroscopic stress difference for both 3.5 and 6.6 GPa peak stresses, suggesting considerable variation in local stress deviators. The CDD sizes determined here are comparable to characteristic length scales for defect-free regions determined from defect density measurements in post-shock recovery experiments. As a result, the present work represents an important step in understanding material microstructure and inhomogeneities in the shock-compressed state.
36 MATERIALS SCIENCE↗
Rapid particle generation from an STL file and related issues in the application of material point methods to complex objects
Abstract In this paper, we focus on three issues related to applications of material point methods (MPMs) to objects with complex geometries. They are material point generation, compatibility of material points with a mesh, and sensitivity to mesh orientation. An efficient method of generating material points from a stereolithography (STL) file is introduced. This material point generation method is independent of the mesh used in MPM calculations. The compatibility between the material points and the mesh is then studied. We also show that the original MPM and the dual domain material point (DDMP) method are sensitive to mesh orientation. These issues are related to the calculation of the internal force and are concerns of the MPMs. They become more prominent when MPMs are applied to complex geometries. Our numerical results show that the recently developed local stress difference (LSD) algorithm (Perez et al. in J Comp Phys 498:112681, 2024) can be used to effectively address them.
36 MATERIALS SCIENCE↗
Statistics of internal stress fluctuations in dislocated crystals and relevance to density-based dislocation dynamics models
A statistical analysis of internal stress fluctuations, defined as the difference between the local mean stress and stress on dislocations, is presented for deforming crystals with 3D discrete dislocation systems. Dislocation realizations are generated using dislocation dynamics simulations and the associated stress field is computed as a superposition of a regularized stress field of dislocation lines within the domain of the solution and a complementary stress field computed via a finite-element boundary value problem. The internal stress fluctuations of interest are defined by an ensemble of the difference between the stress on dislocation lines and the local mean field stress in the crystal. The latter is established in a piecewise fashion over small voxels in the crystal thus allowing the difference between the local average stress and stress on segments to be easily estimated. The results show that the Schmid stress (resolved shear stress) and Escaig stress fluctuations on various slip systems sampled over a random set of points follow a Cauchy (Lorentz) distribution at all strain levels, with the amplitude and width of the distribution being dependent on the strain. Finally, the implications of the Schmid and Escaig internal stress fluctuations are discussed from the points of view of dislocation cross-slip and the dislocation motion in continuum dislocation dynamics.
36 MATERIALS SCIENCE↗
Intermittent Criticality Multi‐Scale Processes Leading to Large Slip Events on Rough Laboratory Faults
Abstract We discuss data of three laboratory stick‐slip experiments on Westerly Granite samples performed at elevated confining pressure and constant displacement rate on rough fracture surfaces. The experiments produced complex slip patterns including fast and slow ruptures with large and small fault slips, as well as failure events on the fault surface producing acoustic emission bursts without externally‐detectable stress drop. Preparatory processes leading to large slips were tracked with an ensemble of ten seismo‐mechanical and statistical parameters characterizing local and global damage and stress evolution, localization and clustering processes, as well as event interactions. We decompose complex spatio‐temporal trends in the lab‐quake characteristics and identify persistent effects of evolving fault roughness and damage at different length scales, and local stress evolution approaching large events. The observed trends highlight labquake localization processes on different spatial and temporal scales. The preparatory process of large slip events includes smaller events marked by confined bursts of acoustic emission activity that collectively prepare the fault surface for a system‐wide failure by conditioning the large‐scale stress field. Our results are consistent overall with an evolving process of intermittent criticality leading to large failure events, and may contribute to improved forecasting of large natural earthquakes.
Geochemistry & Geophysics↗
Regional cutaneous microvascular flow responses during gravitational and LBNP stresses
Due to the regional variability of local hydrostatic pressures, microvascular flow responses to gravitational stress probably vary along the length of the body. Although these differences in local autoregulation have been observed previously during whole-body tilting, they have not been investigated during application of artificial gravitational stresses, such as lower body negative pressure or high gravity centrifugation. Although these stresses can create equivalent G-levels at the feet, they result in distinct distributions of vascular transmural pressure along the length of the body, and should consequently elicit different magnitudes and distributions of microvascular response. In the present study, the effects of whole-body tilting and lower body negative pressure on the level and distribution of microvascular flows within skin along the length of the body were compared.
Breit, Gregory A.↗
Isostatic geoid anomalies on a sphere
Long-wavelength geoid anomalies due to lateral variations in the density or thickness of a thin layer in local isostatic equilibrium on the surface of a sphere depend sensitively on the assumed state of stress within the layer. A number of common intuitive definitions of local isostasy generally believed to be essentially identical are associated with quite different states of isostatic stress, and the corresponding theoretical geoid anomalies can vary by more than a factor of 2. This sensitivity of the theoretical anomaly to the exact definition of local isostasy constitutes an obstacle to any proposed program of interpreting or exploiting observed global isostatic geoid anomalies such as those associated with the aging of the oceanic lithosphere.
Dahlen, F. A.↗
Microstructural changes of proton irradiated Hastelloy-N and in situ micropillar compression testing of one single grain at different local damage levels
Here, in situ micropillar compression was used to study the deformation of proton-irradiated Hastelloy-N at different damage levels. Multiple pillars were prepared from a single domain along the cross-section of 2.5 MeV proton-irradiated Hastelloy-N. Depending on the location of micropillars, the critical resolved shear stress was obtained as a function of local damage levels. Such an approach eliminates the variation of yield stress due to the difference in the Schmid factor. Microstructural characterization showed that irradiation hardening is contributed by complicated defect structures, including (a) dislocation loops with many in corduroy-like alignments, (2) dislocations pile up, and (3) Cottrell atmosphere. Silicon atoms are found to segregate at dislocation lines and loops and form complicated Cottrell cloud patterns. The critical resolved shear stress, $t_{crss}$, and hardening exponents were both extracted as a function of displacements per atom values up to 2.3. There was a 60% increase in $t_{crss}$ at the highest damage level.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
[Presentation Slides] Code-to-Code Benchmark Study for Thermal Stress Modeling and Preliminary Analysis of the High-temperature Single Heat-Pipe Experiment
In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios. This presentation slides discuss finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL), using two commercial software packages, Abaqus and Star-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain preliminary insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the inner surface of the heat pipe hole surrounded by electric heaters where the largest temperature gradients appear. The temperature fields have good agreements between Abaqus and Star-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes. It is found that the local maximum thermal stresses may reach close to the ultimate tensile strength and yield strength of structural material depending on the heater power. Ultimately, the coupled thermal-structural analysis will help guide the current experimental plan and ensure the facility safety for future experimental study.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
A dynamic localization model with stochastic backscatter
The modeling of subgrid scales in large-eddy simulation (LES) has been rationalized by the introduction of the dynamic localization procedure. This method allows one to compute rather than prescribe the unknown coefficients in the subgrid-scale model. Formally, the LES equations are supposed to be obtained by applying to the Navier-Stokes equations a 'grid filter' operation. Though the subgrid stress itself is unknown, an identity between subgrid stresses generated by different filters has been derived. Although preliminary tests of the Dynamic Localization Model (DLM) with k-equation have been satisfactory, the use of a negative eddy viscosity to describe backscatter is probably a crude representation of the physics of reverse transfer of energy. Indeed, the model is fully deterministic. Knowing the filtered velocity field and the subgrid-scale energy, the subgrid stress is automatically determined. We know that the LES equations cannot be fully deterministic since the small scales are not resolved. This stems from an important distinction between equilibrium hydrodynamics and turbulence. In equilibrium hydrodynamics, the molecular motions are also not resolved. However, there is a clear separation of scale between these unresolved motions and the relevant hydrodynamic scales. The result of molecular motions can then be separated into an average effect (the molecular viscosity) and some fluctuations. Due to the large number of molecules present in a box with size of the order of the hydrodynamic scale, the ratio between fluctuations and the average effect should be very small (as a result of the 'law of large numbers'). For that reason, the hydrodynamic balance equations are usually purely deterministic. In turbulence, however, there is no clear separation of scale between small and large eddies. In that case, the fluctuations around a deterministic eddy viscosity term could be significant. An eddy noise would then appear through a stochastic term in the subgrid-scale model and could be the source of backscatter.
Carati, Daniele↗
Utah FORGE: 2024 Discrete Fracture Network Model Data
The Utah FORGE 2024 Discrete Fracture Network (DFN) Model dataset provides a set of files representing discrete fracture network modeling for the FORGE site near Milford, Utah. The dataset includes four distinct DFN model file sets, each corresponding to different time frames and modeling approaches in 2024. These models characterize both natural and induced fractures in the geothermal reservoir, which consists of crystalline granitic and metamorphic rock approximately 8,000 feet below the ground surface. The dataset includes a reference DFN model from February 2024 that incorporates planar fractures and well trajectories, as well as upscaled permeability, porosity, compressibility, and storage values on specified grids. Additionally, there are models based on new microseismic (MEQ) data from May and July 2024, including fracture planes fitted to the latest MEQ catalog datasets, tensile fractures from hydraulic stimulation, and an alternative connected DFN for modeling purposes. Coordinate data is provided in both global and local frames, with detailed instructions on the transformations used to align with principal stress orientations. The dataset also includes notes and calculation files for estimating fracture sizes and differences between various fracture sets. There are subfolders for Global Coordinates and Local Coordinates. To move from the global to the local coordinate frame, fractures and wells were a) rotated 20 degrees counterclockwise looking down about the global point (335376.400482041, 4263189.99998761, 250.093546450195) to better align with the principal stresses; and b) translated by (-335408.68, -4263010.9, 1150). Upscaled permeability values using the _XYZ suffix show directions with respect to the global XYZ coordinate frame, while those using the _IJK suffix are aligned with local coordinate frame.
15 GEOTHERMAL ENERGY↗
Predicting Fracture Network Development in Crystalline Rocks
The geometric properties of fractures influence whether they propagate, arrest, or coalesce with other fractures. Thus, quantifying the relationship between fracture network characteristics may help predict fracture network development, and perhaps precursors to catastrophic failure. To constrain the relationship and predictability of fracture characteristics, we deform eight one centimeter tall rock cores under triaxial compression while acquiring in situ X-ray tomograms. The tomograms reveal precise measurements of the fracture network characteristics above the spatial resolution of 6.5 µm. We develop machine learning models to predict the value of each characteristic using the other characteristics, and excluding the macroscopic stress or strain imposed on the rock. The models predict fracture development more accurately in the experiments performed on granite and monzonite, than the experiments on marble. Fracture network development may be more predictable in these igneous rocks because their microstructure is more mechanically homogeneous than the marble, producing more systematic fracture development that is not strongly impeded by grain contacts and cleavage planes. The varying performance of the models suggest that fracture volume, length, and aperture are the most predictable of the characteristics, while fracture orientation is the least predictable. Orientation does not correlate with length, as suggested by the idea that the orientation evolves with increasing differential stress and thus fracture length. This difference between the observed and expected relationship between orientation and length highlights the influence of mechanical heterogeneities and local stress perturbations on fracture growth as fractures propagate, link, and coalesce.
58 GEOSCIENCES↗
The Role of Stratigraphy and Loading History in Generating Complex Compaction Bands in Idealized Field-Scale Settings
The Buckskin Gulch locality in Utah is a landmark example of compaction localization. The outcrop of this locality involves distinct stratigraphic heterogeneity and was exposed to complex loading history. It features multiple sets of deformation bands with different kinematics and orientation. Similar formations were seen in the Valley of Fire, Nevada, and the Orange quarry, France, among other localities. The formation of such complex structures, their propagation mechanisms, and frequency is affected by numerous local and ambient factors whose impacts are not yet fully understood. The simulation of the above-mentioned localities is not feasible because of the limited amount of available information. This work, instead, investigates from a geomechanics standpoint how the interplay among material nonlinearity, outcrop stratigraphy, and loading history interconnects with specific spatiotemporal patterns of compaction band propagation. Our study shows that the system stratigraphy can be responsible for the emergence of coexisting compaction bands with different inclination and kinematics. Specifically, we show that stiffness contrasts induce nonlocal stress changes which may favor the initiation of secondary structures with different compaction localization characteristics. Furthermore, systems of inclined compaction bands induced by burial increase display secondary, noncontemporaneous sets of vertical compaction bands under the effects of postburial shortening. Our results indicate that stages of intermediate burial decrease prior to tectonic shortening can promote the formation of such complex systems. Despite the simplifications involved in our analyses, these findings show how geomechanics computations complement field observations and could provide a mechanics-based validation of site-specific reconstruction hypothesis.
58 GEOSCIENCES↗
The 3-D inelastic analysis methods for hot section components: Brief description
Advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate yet more cost-effective analysis of components subjected to severe thermal gradients and loads in the presence of mechanical loads, with steep stress and strain gradients are being developed. Anisotropy, time and temperature dependent plasticity and creep effects are also addressed. The approach is to develop four different theories, one linear and three higher order theories (polynomial function, special function, general function). The theories are progressively more complex from linear to general function in order to provide streamlined analysis capability with increasing accuracy for each hot section component and for different parts of the same component according to the severity of the local stress, strain and temperature gradients associated with hot spots, cooling holes and surface coating cracks. To further enhance the computational effectiveness, the higher order theories will have embedded singularities (cooling passages, for example) in the generic modeling region. Each of the four theories consists of three formulation models derivable from independent theoretical formulations. These formulation models are based on: (1) mechanics of materials; (2) special finite elements; and (3) an advanced formulation to be recommended by the contractor.
Chamis, C. C.↗
The Influence of Residual Stress on Fatigue Crack Growth Rates in Stainless Steel Processed by Different Additive Manufacturing Methods
The properties and microstructure of Type 304L stainless steel produced by two additive manufacturing (AM) methods—directed energy deposition (DED) and powder bed fusion (PBF)—are evaluated and compared. Localized heating and steep temperature gradients of AM processes lead to significant residual stress and distinctive microstructures, which may be process-specific and influence mechanical behavior. Test data show that materials produced by DED and PDF have small differences in tensile strengths but clear differences in residual stress and microstructural features. Measured fatigue crack growth rates (FCGRs) for cracks propagating parallel to and perpendicular to the build directions differ between the two AM materials. To separate the influences of residual stress and microstructure, K-control test procedures with decreasing and constant stress intensity factor ranges are used to measure FCGRs in the near-threshold regime (crack growth rates ≤ 1 × 10 −8 m/cycle). Residual stress is quantified by the residual stress intensity factor, K res , measured by the online crack compliance method. Correcting the FCGR data for differences in K res brings results for specimens of the two AM materials into agreement with each other and with results for wrought specimens, when the latter are corrected for crack closure. Differences in microstructure and tensile strength have an insignificant influence on FCGRs in these tests.
additive manufacturing↗
Chemo‐Mechanics Interplay Dictates Interface Instability and Asymmetry in Plating and Stripping of Sodium Metal Electrodes
Abstract The development of practical sodium (Na) metal batteries is hindered by key challenges including dendrite growth, dead metal formation, and unstable solid electrolyte interphase (SEI) growth. A fundamental understanding of the chemo‐mechanical interactions at the Na/SEI interface is critical for designing stable Na metal electrodes. In this work, the coupled electrochemical‐mechanical processes governing the morphological evolution and stability of Na metal during plating and stripping are investigated. The heterogeneous nature of transport and morphological interactions at the Na/SEI interface is revealed to result in nonuniform mechanical overpotentials and reaction fronts, eventually leading to Na filaments or pits. The spatio‐temporal evolution of stress heterogeneities during Na plating and stripping is shown to be asymmetric, manifesting in varying morphological nonuniformities and instability modes. The crucial role of external pressure in modulating the electrochemical‐transport interactions, the mechanical response of Na, and the localized reaction currents and stresses at the Na/SEI interface is demonstrated. At different external pressure conditions, the correlation between transport heterogeneities in the SEI and the onset and propagation of interface instability has been delineated. This work highlights the need for synergistic tailoring of external pressure and SEI heterogeneity toward achieving stable electrodeposition and dissolution in Na metal electrodes.
Singla, Aditya [School of Mechanical Engineering P↗
Directed aging, memory, and nature’s greed
Disordered materials are often out of equilibrium and evolve very slowly in a rugged and tortuous energy landscape. This slow evolution, referred to as aging, is deemed undesirable as it often leads to material degradation. However, we show that aging also encodes a memory of the stresses imposed during preparation. Because of inhomogeneous local stresses, the material itself decides how to evolve by modifying stressed regions differently from those under less stress. Because material evolution occurs in response to stresses, aging can be “directed” to produce sought-after responses and unusual functionalities that do not inherently exist. Aging obeys a natural “greedy algorithm” as, at each instant, the material simply follows the path of most rapid and accessible relaxation. Our experiments and simulations illustrate directed aging in examples in which the material’s elasticity transforms as desired because of an imposed deformation.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Impacts of Substrate Thinning on FPGA Performance and Reliability [Slides]
Substrate thinning is necessary in devices with flip-chip BGA packages to enable both radiation testing and component qualification and high-spatial resolution beam-based failure analysis methods. We investigated three factors affecting device performance: subsurface damage from the thinning process, reduced heat spreading in thin substrates, and changes in device switching speed. We conclude subsurface damage to crystalline Si caused by the thinning process is removable with sufficient SiO 2 slurry polishing. Local temperature differences increase minimally in devices thinned to 3 μm. Compressive stress in the Si increases globally after device thinning and leads to slowing of ring oscillator frequency by about 0.5% compared to full-thickness devices. Future work will include extending the results to submicron Si thickness values, which also has important benefits for failure analysis, debug, and security assessments. We also plan to extend this type of work to other FPGAs and other devices like memory and processors.
36 MATERIALS SCIENCE↗