Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress field”

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 145 records · Page 8

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↗

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

A parallel discrete dislocation dynamics/kinetic Monte Carlo method to study non-conservative plastic processes

Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.

36 MATERIALS SCIENCE↗

Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions

Understanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo–Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High–resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first–order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. Here, the results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.

58 GEOSCIENCES↗

Large deformation and brittle failure calculated using the dual-domain material point method

The dual domain material point (DDMP) method is explored as a candidate to be implemented in a general purpose code to perform simulations of materials with complex geometry that undergo large history-dependent deformation and failure. To test its candidacy, we study its mesh convergence, its sensitivity to mesh orientation, and its ability to handle softening and failure of a material. Simulations of large deformation and simulations of mechanical failure are performed using both DDMP and the material point method (MPM). When cell-crossing of material points is not an issue and when there are a sufficient number of material points in each computation cell, the numerical error decreases with the square of the cell size as expected for both MPM and DDMP. DDMP has reduced error compared with MPM when there are many instances of material points crossing cell boundaries due to the continuous nature of the modified gradient of the shape functions. Simulations of a specimen under tension are also performed where the background mesh is aligned and misaligned with the tension direction. MPM displays a significant mesh-dependent stress field, DDMP shows negligible mesh dependency. Despite a mesh orientation-dependent stress field from MPM, the critical tension and failure mode from both MPM and DDMP calculations have negligible mesh dependency when using a non-local failure model. If only the failure mode is important (i.e., local stresses are unimportant), MPM with a non-local failure model is a suitable method for modeling failure with small deformations. However, if local stresses are also important or if there are large deformations with many cell-crossings before failure, DDMP should be the method that is used. A needed improvement for DDMP is identified from our numerical simulations.

36 MATERIALS SCIENCE↗

Unraveling the dislocation–precipitate interactions in high-entropy alloys

The precipitates play a significant role in not only enhancing the strength, but also maintaining the high toughness in alloys. However, the interactions of the nanoscale precipitates with dislocations in the high entropy alloys (HEAs) are difficult to observe directly by in-situ TEM experiments due to the limits of the resolution and time. Here, using atomic simulations we report the synergistic strengthening of the coherent precipitate and atomic-scale lattice distortion in the HEAs at cryogenic/elevated temperatures. The effects of temperature, chemical disorder, precipitate spacing, precipitate size, elemental segregation, and dislocation-cutting number on the critical stress for the dislocation to overcome a row of precipitates are studied. A random stacking fault energy landscape along the slip plane, the lattice distortion at different temperatures, and the interface/surface energy at various crystallographic orientations are obtained. Compared with the traditional metals and alloys, HEAs have the severe atomic-scale lattice distortions to generate the local high tensile/compressive stress fields. This complex stress causes the dislocation line to bend, and thus improves the dislocation slip resistance, resulting in the strong solid-solution strengthening. The stacking fault strengthening induced by the obvious difference of the stacking fault energies between the HEA matrix and precipitate (within the inner of the HEA matrix), and the formation of the antiphase domain boundary contribute to the high strength. The precipitate embedded by the solute atoms produces the strong lattice distortion to enhance the dislocation slip resistance at high temperatures. Hence, the current results provide the mechanistic insight into the phenomenon that the coherent precipitate combined with the severe atomic-scale lattice distortion can enhance the strength at cryogenic/elevated temperatures to further broaden the scope of applications of advanced HEAs.

36 MATERIALS SCIENCE↗

Sound Assisted Low Temperature Spalling: Upscaling and Throughput

In evaluating the state-of-the-art spalling techniques, the key factor that is missing is the control of crack propagation, which is responsible for low quality wafers due to high surface roughness and Total Thickness Variation, TTV. This project has shown that acoustic waves can alter the conditions at the crack tip during thermal spalling, modifying the propagation speed and the stress field around it (stress intensity factor, K). Both factors are directly related to the roughness of the spalled surface, so as we control speed and stress, we can fine-tune the surface roughness through the entire wafer as suggested by Arakawa.

42 ENGINEERING↗

Nitrogen Vacancy Center in Diamond for the Stress and Field Sensing Applications

The nitrogen-vacancy (NV) center in a nanodiamond (ND) crystal is a promising material for quantum information processing, sensing, and computing applications. It is one of the best candidate materials for quantum sensing and metrology expected to work at elevated temperatures and pressures conditions. We computationally show the effect of strain on the defect band edges and band gaps in the NV center diamond. A low energy Hamiltonian is developed for the ±1 spin manifold at the ground state. We show the quantum sensing device is a few orders of magnitude superior in sensing than the traditional optical sensing devices. We also discuss experimental results from the optically detected magnetic resonance (ODMR) and the spin relaxometry for the field sensing applications. The presentation concludes by providing a model for free spins detection of the rare earth ions.

field sensing applications↗

Ethiopian Tertiary dike swarms

Mapping of the Ethiopian rift and Afar margins revealed the existence of Tertiary dike swarms. The structural relations of these swarms and the fed lava pile to monoclinal warping of the margins partly reflect a style of continental margin tectonics found in other parts of the world. In Ethiopia, however, conjugate dike trends appear to be unusually strongly developed. Relation of dikes to subsequent margin faulting is ambiguous, and there are instances where the two phenomena are spatially separate and of differing trends. There is no evidence for lateral migration with time of dike injection toward the rift zone. No separate impingement of Red Sea, Gulf of Aden, and African rift system stress fields on the Ethiopian region can be demonstrated from the Tertiary dike swarms. Rather, a single, regional paleostress field existed, suggestive of a focus beneath the central Ethiopian plateau. This stress field was dominated by tension: there is no cogent evidence for shearing along the rift margins. A gentle compression along the rift floor is indicated. A peculiar sympathy of dike hade directions at given localities is evident.

Mohr, P. A.↗

Separating the Influence of Environment from Stress Relaxation Effects on Dwell Fatigue Crack Growth in a Nickel-Base Disk Alloy

Both environmental embrittlement and crack tip visco-plastic stress relaxation play a significant role in determining the dwell fatigue crack growth (DFCG) resistance of nickel-based disk superalloys. In the current study performed on the Low Solvus High Refractory (LSHR) disk alloy, the influence of these two mechanisms were separated so that the effects of each could be quantified and modeled. Seven different microstructural variations of LSHR were produced by controlling the cooling rate and the subsequent aging and thermal exposure heat treatments. Through cyclic fatigue crack growth testing performed both in air and vacuum, it was established that four out of the seven LSHR heat treatments evaluated, possessed similar intrinsic environmental resistance to cyclic crack growth. For these four heat treatments, it was further shown that the large differences in dwell crack growth behavior which still persisted, were related to their measured stress relaxation behavior. The apparent differences in their dwell crack growth resistance were attributed to the inability of the standard linear elastic fracture mechanics (LEFM) stress intensity parameter to account for visco-plastic behavior. Crack tip stress relaxation controls the magnitude of the remaining local tensile stresses which are directly related to the measured dwell crack growth rates. It was hypothesized that the environmentally weakened grain boundary crack tip regions fail during the dwells when their strength is exceeded by the remaining local crack tip tensile stresses. It was shown that the classical creep crack growth mechanisms such as grain boundary sliding did not contribute to crack growth, but the local visco-plastic behavior still plays a very significant role by determining the crack tip tensile stress field which controls the dwell crack growth behavior. To account for the influence of the visco-plastic behavior on the crack tip stress field, an empirical modification to the LEFM stress intensity parameter, Kmax, was developed by incorporating into the formulation the remaining stress level concept as measured by simple stress relaxation tests. The newly proposed parameter, Ksrf, did an excellent job in correlating the dwell crack growth rates for the four heat treatments which were shown to have similar intrinsic environmental cyclic fatigue crack growth resistance.

Dwell crack growth↗

Separating the Influence of Environment from Stress Relaxation Effects on Dwell Fatigue Crack Growth

Seven different microstructural variations of LSHR were produced by controlling the cooling rate and the subsequent aging and thermal exposure heat treatments. Through cyclic fatigue crack growth testing performed both in air and vacuum, it was established that four out of the seven LSHR heat treatments evaluated, possessed similar intrinsic environmental resistance to cyclic crack growth. For these four heat treatments, it was further shown that the large differences in dwell crack growth behavior which still persisted, were related to their measured stress relaxation behavior. The apparent differences in their dwell crack growth resistance were attributed to the inability of the standard linear elastic fracture mechanics (LEFM) stress intensity parameter to account for visco-plastic behavior. Crack tip stress relaxation controls the magnitude of the remaining local tensile stresses which are directly related to the measured dwell crack growth rates. It was hypothesized that the environmentally weakened grain boundary crack tip regions fail during the dwells when their strength is exceeded by the remaining local crack tip tensile stresses. It was shown that the classical creep crack growth mechanisms such as grain boundary sliding did not contribute to crack growth, but the local visco-plastic behavior still plays a very significant role by determining the crack tip tensile stress field which controls the dwell crack growth behavior. To account for the influence of the visco-plastic behavior on the crack tip stress field, an empirical modification to the LEFM stress intensity parameter, Kmax, was developed by incorporating into the formulation the remaining stress level concept as measured by simple stress relaxation tests. The newly proposed parameter, Ksrf, did an excellent job in correlating the dwell crack growth rates for the four heat treatments which were shown to have similar intrinsic environmental cyclic fatigue crack growth resistance.

environmental damage↗

The mode 3 crack problem in bonded materials with a nonhomogeneous interfacial zone

The mode 3 crack problem for two bonded homogeneous half planes was considered. The interfacial zone was modelled by a nonhomogeneous strip in such a way that the shear modulus is a continuous function throughout the composite medium and has discontinuous derivatives along the boundaries of the interfacial zone. The problem was formulated for cracks perpendicular to the nominal interface and was solved for various crack locations in and around the interfacial region. The asymptotic stress field near the tip of a crack terminating at an interface was examined and it was shown that, unlike the corresponding stress field in piecewise homogeneous materials, in this case the stresses have the standard square root singularity and their angular variation was identical to that of a crack in a homogeneous medium. With application to the subcritical crack growth process in mind, the results given include mostly the stress intensity factors for some typical crack geometries and various material combinations.

Erdogan, Fazil↗

The mode III crack problem in bonded materials with a nonhomogeneous interfacial zone

The mode 3 crack problem for two bonded homogeneous half planes was considered. The interfacial zone was modelled by a nonhomogeneous strip in such a way that the shear modulus is a continuous function throughout the composite medium and has discontinuous derivatives along the boundaries of the interfacial zone. The problem was formulated for cracks perpendicular to the nominal interface and was solved for various crack locations in and around the interfacial region. The asymptotic stress field near the tip of a crack terminating at an interface was examined and it was shown that, unlike the corresponding stress field in piecewise homogeneous materials, in this case the stresses have the standard square root singularity and their angular variation was identical to that of a crack in a homogeneous medium. With application to the subcritical crack growth process in mind, the results given include mostly the stress intensity factors for some typical crack geometries and various material combinations.

Erdogan, F.↗

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗

Thermoelastic stress - How important as a cause of earthquakes in young oceanic lithosphere?

Thermoelastic or thermal stress is a potentially important contributor to the state of stress in the oceanic lithosphere. The present paper provides several simple models for the state of thermoelastic stress in a young oceanic lithosphere, taking into account a comparison of the predictions of these models with the characteristics of near-ridge earthquakes. Attention is given to the characteristics of near-ridge earthquakes, sources of stress in an oceanic lithosphere, previous models of thermal stress, the calculation of thermal stress, and thermal stress models. A test is conducted of the hypothesis that thermoelastic stress is a significant component of the stress field in a young oceanic lithosphere. The considered models support the hypothesis that thermoelastic stress is a significant component of the stress field in a young oceanic lithosphere.

Bratt, S. R.↗

Ultrasonic measurement of stress in 2219-T87 aluminum plate

The basic relationship of ultrasonic signal velocity to directional subsurface stress is reviewed. Inappropriateness of dependency on a single correlative value of constant for a three dimensional stress field in metallic materials is discussed. Implementation of conventional ultrasonic nondestructive testing capabilities integrated to provide a composite technique for the measurement of orthogonal stress components is described, and the procedures for performing the preparatory calibration and subsequent stress field measurements are presented. In conclusion, the prime effect of stress on ultrasonic signal velocity occurs only in the direction of material excitation or particle motion.

Clotfelter, W. N.↗

Wavy and Cycloidal Lineament Formation on Europa from Combined Diurnal and Nonsynchronous Stresses

In a companion abstract, we show that fractures propagated into combined diurnal and nonsynchronous rotation (NSR) stress fields can be cycloidal, "wavy," or arcuate in planform as the relative proportion of NSR stress in increased. These transitions occur as NSR stress accumulates over approx. 0 to 10 deg of ice shell rotation, for average fracture propagation speeds of approx. 1 to 3 m/s. Here we consider the NSR speed parameter space for these morphological transitions, and explore the effects on cycloids of adding NSR to diurnal stress. Fitting individual Europan lineaments can constrain the combined NSR plus diurnal stress field at the time of formation.

Gleeson, Damhnait↗

Green's function fast multipole method for continuum mechanics (SM-FMM)

Solid Mechanics Fast Multipole Method based on elastic Green's function and accelerated using FFTs The code calculates the mechanical fields (stress and strain) for a heterogeneous elasto-plastic material unit cell under quasi-static conditions (no dynamic effects). The fields are calculated using a Green's function method, where the strain field is given by a discrete convolution of Green's operator with an auxiliary stress field. The convolution is calculated using the fast multipole method.

Zecevic, Miroslav↗