Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS CALCULATION”

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

High-burnup fuel stress analysis prior to and during a LOCA transient

High-burnup fuel fragmentation and pulverization have been experimentally observed during simulated loss-of-coolant accident (LOCA) transient testing of high-burnup (>62 GWd/tU) nuclear fuel. The cause of this phenomenon is generally accepted to be the formation of stresses arising from the off-normal conditions. However, the primary mechanism(s) contributing to the stress formation is/are still being debated. Leading hypotheses suggest that the concurrence of fission gas bubbles and high burnup structure on the pellet periphery (i.e., the rim) leave the fuel susceptible to pulverization. Other hypotheses suggest that pulverization is artificially inflated during externally heated tests that do not reproduce prototypic fuel stresses and radial temperature profiles as the fuel transitions from full-power nuclear heating to decay heating. This work performs finite element analyses of fuel pellet temperatures and macroscopic stresses prior to and during a LOCA transient to identify (1) operating conditions that may exacerbate the fuel pellet stress state, (2) differences in predicted stresses with nuclear vs. furnace heating, and (3) how fuel cracking impacts the subsequent stress state in the pellet fragments. Analyses predict that macroscopic transient fuel stresses resulting from constrained thermal expansion will remain below the fracture stress of fresh UO 2 once the fragment size reaches ~4 mm or less. Calculated fuel stresses from simulated LOCA tests performed with representative nuclear vs. electrical heating rates are shown to be equal within 5% and essentially unchanged for all times after the first 27 seconds of the transients. Because the fuel stresses do not change once the temperature profile becomes uniform, macroscopic transient stresses would be no different than those following a normal shutdown. Because pulverization of high burnup fuel does not occur prior to a LOCA test, pulverization of high burnup UO 2 cannot be initiated by differential thermal stresses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Determination of In-Situ Stress and Geomechanical Properties from Borehole Deformation

This paper proposes a cost-effective technique to determine geomechanical properties and in-situ stress from borehole deformation data. In this approach, an artificial neural network (ANN) is applied to map the relationship among in-situ stress, borehole size, geomechanical properties, and borehole displacements. The genetic algorithm (GA) searches for the set of unknown stresses and geomechanical properties that matches the objective borehole deformation function. Probabilistic recapitulation (PR) analysis is conducted after each ANN-GA modeling cycle and will be repeated with a reduced number of unknowns for the next ANN-GA modeling cycle until unequivocal results are achieved. The PR-GA-ANN method has been demonstrated by a field case study to estimate borehole size, Young’s modulus, Poisson’s ratio, and the two horizontal stresses using borehole deformation information reported from four-arm caliper log of a vertical borehole. The methodology effectively solves the issue of the multiple solutions (various rock mechanical properties and in-situ stresses combinations) for a certain borehole deformation. The case study also indicated that the calculated horizontal stresses are in reasonable agreement with the filed hydraulic fracture treatment observations and the reported regional stress study of the area.

58 GEOSCIENCES↗
DOE OSTI2021

Disentangling the Effects of Vapor Pressure Deficit and Soil Water Availability on Canopy Conductance in a Seasonal Tropical Forest During the 2015 El Niño Drought

In land models used in Earth system models, vegetation water stress is heuristically related to soil water matric potential, and root water uptake in soil layers is simply proportional to root density to meet the transpiration demand. The efficiency of roots utilizing deeper soil water under water stress cannot be well represented, but it is of critical importance for transpiration regulation and carbon flux. To improve modeling of tropical forest response to drought, a plant hydrodynamic scheme (HYDRO), adapted from a trait-based model in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), has been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to represent the hydraulic dynamics and root water uptake strategy. The adapted model conforms to the big-leaf photosynthesis model in ELM. Vegetation water stress in the model is a function of prognostic leaf water potential and is used to regulate stomatal conductance. Point simulations were performed at two tropical forest sites in the Amazon and Panama. Model sensitivity to the numerical solver, plant hydraulic traits, soil properties, and atmospheric forcing has been evaluated. Our results showed that, even though water stress factors calculated at our study sites using HYDRO are comparable to that simulated by the heuristic water stress model, soil moisture prediction is improved with soil water redistributions through roots. Simulated results using HYDRO is not as sensitive to plant hydraulic traits measured for different tree species as to the soil properties and atmospheric forcing. Simplified numerical solver may cause inaccuracy as to where roots uptake water. This study suggests that plant hydrodynamic representation is critical for improved simulations of soil moisture under water stresses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Atomistic Investigation of Plastic Deformation and Dislocation Motion in Uranium Mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. However, many aspects of its mechanical behavior, particularly at reactor-relevant conditions, remain unclear. In this study, molecular dynamics (MD) simulations were employed to investigate the deformation behavior and dislocation motion in UN. We found that the Kocevski potential predicts the principal slip system as $\frac{1}{2}$ $\langle110\rangle${110}, aligning with experimental data. On the other hand, the Tseplyaev potential predicts slip to primarily occur on $\frac{1}{2}$ $\langle110\rangle${111}. MD simulations of stress–strain behavior were used to estimate the nanoindentation hardness, revealing that the Kocevski potential accurately predicts hardness even though it fails to model dynamic plasticity. Complete dislocation mobility functions have been fitted for the edge and screw dislocations in both the thermally activated and phonon-drag regimes. The 300 K linear mobility of the edge dislocation using the Tseplyaev potential was found to be 817 Pa -1 ·s -1 , whereas that of the screw dislocation using the Kocevski potential was found to be 4546 Pa -1 ·s -1 . At intermediate stresses, we observed that the subsonic steady-state motion of the edge dislocation in UN is intermittently interrupted by velocity jumps, reaching the average sound velocity. Finally, the threshold Schmid stress is calculated as 179–197 MPa, which gives an upper-limit estimate of the uniaxial yield stress of polycrystalline UN of 548–603 MPa. These findings, including the fitted dislocation mobility function, provide essential input for future plasticity and dislocation dynamics models of nuclear fuels.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Integrated solver for fluid driven fracture and fragmentation

An integrated solver may be used for both a fluid phase and a solid phase. The integrated solver may use constitutive equations for the fluid phase and the solid phase, and calculate the stress in the fluid and solid using the same integration point. The integrated solver may also calculate other state variables using the same integration point.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

The Reynolds stress induced by weakly up-down asymmetric ion temperature gradient mode

A weakly up-down asymmetric ion temperature gradient (ITG) mode (formed in the up-down symmetric equilibrium), localized in the bad curvature region, may contribute dominantly to the nonlinear Reynolds stress because its growth rate is higher than the competing modes: the strongly up-down asymmetric modes, and other modes peaked (poloidally) away from the outboard mid-plane. Based on the fluid ITG model, the parameterized analytic wave function is constructed in the two dimensional (2D) Fourier-ballooning representation, and numerically verified by making use of the finite difference methods. Analytic expressions for the poloidal Reynolds stresses, defined on a rational surface, are given for both the weakly and the strongly up-down asymmetric ITG modes. It is shown here that the poloidal stresses of both modes can be interpreted as a superposition of a monopole and a dipole structure. For comparison with the experiment, the stresses are calculated for L-mode discharge parameters on five machines. Analytic expressions for the total parallel Reynolds stresses (contributed from all rational surfaces) are also derived for both modes: they are also exhibited, graphically, using a Gaussian wave intensity model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Stress field and interaction forces between dislocations and precipitate distributions

A computational method is developed for calculation of the stress field and interaction forces between dislocations and precipitates of arbitrary shape and distribution. The internal stress generated by precipitates due to coherency strain is implemented within the discrete dislocation dynamics (DDD) framework. The s-version finite element method (s-FEM), which models a precipitate of arbitrary shape using a local mesh is used to calculate coherency stress fields. The method facilitates meshing precipitate volumes of arbitrary geometry, and multiple local meshes can be superimposed at various positions of the global mesh. Accuracy and convergence conditions of the method are established. For a single precipiate, the method is shown to be 3.5 times faster than a standard FEM calculation for the same accuracy, and the gain in speed increases with the number of precipitates. The CRSS of spherical precipitates is found to be higher than disk-shaped precipitates of the same volume fraction.

42 ENGINEERING↗
DOE OSTI2026

The effect of dogbone sample size in the tensile testing of TATB-based plastic‑bonded explosive materials

Mechanical properties are of interest for many plastic‑bonded explosive (PBX) materials with tensile properties being of particular interest. Direct tensile measurements using dogbone-shaped samples are considered the gold standard, but they are fairly large, making testing more costly and less desirable from a safety perspective. We investigated whether the measured tensile strength depends on the dogbone specimen size, which to our knowledge, has not been reported in the literature for PBX materials. Understanding this should inform the feasibility of employing smaller samples and how sample size should be considered when comparing PBX dogbone values in the literature. The TATB-based PBX dogbone sample size was varied by (a) scaling all dimensions proportionally and (b) varying only the length of the samples. It was observed that the measured tensile peak stress (strength) was a function of the sample size, and was more dependent on the diameter (cross-sectional area) than the length of the samples. Since peak stress is calculated as peak force normalized to the diameter of the sample, one might not expect an explicit diameter dependence for the peak stress. Therefore, these results suggest there may be an additional strengthening effect as the sample diameter is increased.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Laboratory Evaluation of the Thermal Breakout Method for Maximum Horizontal Stress Measurement

Measuring in situ stress is essential for many problems in geomechanics, and the maximum horizontal stress is the most difficult to constrain. We are developing an extension of the breakout method to measure maximum horizontal stress in regions where natural breakouts do not occur. In the novel thermal breakout method, additional compression which leads to breakout development is induced by heating the borehole wall. In the present study, we validated the method experimentally in a true-triaxial apparatus on samples with predrilled boreholes. Two rocks were selected for laboratory testing: high-porosity Berea sandstone and low-porosity Niagaran dolomite. Prior to main true-triaxial tests, we carried out standard testing to characterize the strength, elasticity and thermal properties. The true-triaxial experiments consisted of: (1) room-temperature tests where samples were first loaded mechanically until the breakout formed, and (2) elevated-temperature tests where samples were loaded mechanically within the elastic range with additional compression induced thermally. Breakout initiation was monitored by acoustic emission sensors mounted on the pistons that applied horizontal stresses. Furthermore, the magnitude of induced thermal stress was calculated from temperature measurements around the borehole wall. In both rock types, we created thermally induced breakouts and examined analytical expressions to constrain maximum horizontal stress based on strength, elastic and thermal properties of the rocks.

58 GEOSCIENCES↗
DOE OSTI2019

Characterization and Modeling of Hot Deformation Behavior of a Copper-Bearing High-Strength Low-Carbon Steel Microalloyed with Nb

This study investigates hot deformation behavior of a newly developed Cu-bearing high-strength low-carbon steel microalloyed with Nb (Nb-HSLC). A computational method based on experimental data was employed to design the chemical composition of the alloy. Compression tests were carried out in the temperature range of 850-1100 °C as well as strain rates of 0.001-10 s{sup −1} using BAHR Dil 805 A/D thermo-analyzer equipment. The Arrhenius-type constitutive equations were used to model the hot working behavior of the designed steel. Effects of friction and temperature rise during deformation were corrected to obtain the actual stresses. The results showed that the peak flow stress was increased with increasing Zener–Hollomon parameter. The obtained flow curves at strain rates lower than 0.1 s{sup −1} and temperatures above 950 °C represented the typical dynamic recrystallization (DRX) behavior, while the flow curves at temperatures lower than 950 °C at all strain rates were associated with continuous strain hardening. This feature is in good agreement with the precipitation temperature range of Nb(C, N) particles, i.e., 800-1000 °C. Moreover, the flow curves showed the serrations during hot deformation at strain rates of 0.001 and 0.01 s{sup −1}, indicating that the dynamic strain aging (DSA) phenomenon occurred at low strain rates. The best fit between “peak stress” and “deformation conditions” was obtained by a hyperbolic sine-type equation (R{sup 2} = 0.993). Therefore, the average activation energy was determined as 348 kJ mol{sup −1}. The agreement between the achieved model and experimental flow data was verified using the results of additional tests at a strain rate of 5 s{sup −1}. The maximum difference between the measured and predicted “peak stresses” was calculated as 5 Mpa.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Pressure and tension effects on mechanical properties of ZrAl{sub 2}

Structural, elastic, thermodynamic of ZrAl{sub 2} under pressure, ideal strength and deformation mode under tension are investigated by the first-principles method. The calculated structural parameters at zero pressure are in consistent with experiments. Under pressure, elastic constants and their pressure dependence are obtained using the static finite strain technique. ZrAl{sub 2} exhibits lower elastic anisotropy. The linear thermal expansion coefficient shows greater effects of temperature at lower pressure. The ideal tensile have been investigated by stress–strain calculations. Finally, the microscopic mechanism that determines the structural stability is studied using the results of electronic structure calculations. We propose that the weakening of Zr-Zr leads to the significant change of stress–strain curve at strain ∼0.27, and the breaking of Zr{sub 2}-Zr{sub 3} leads to the structural instability of ZrAl{sub 2} under large tensile strains.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Geomechanical simulation of the stress tensor rotation caused by injection of cold water in a deep geothermal reservoir

We present a three-dimensional thermohydromechanical numerical study of the evolution and distribution of the stress tensor within the northwest part of The Geysers geothermal reservoir (in California), including a detailed study of the region around one injection well from 2003 to 2012. Initially, after imposing a normal faulting stress regime, we calculated local changes in the stress regime around injection wells. Our results were compared with previously published studies in which the stress state was inferred from inverting the focal plane mechanism of seismic events. Our main finding is that changes in stress tensor orientation are caused by injection-induced progressive cooling of the reservoir, as well as by the seasonal variations in injection rate. Because of the gravity flow and cooling around a liquid zone formed by the injection, the vertical stress reduction is larger and propagates far below the injection well. At the same time, the horizontal stress increases, mostly because of stress redistribution below and above the cooling area. These two phenomena cause the rotation of the stress tensor and the appearance of a strike-slip regime above, inside, and below the cooling area. The cooling and the associated rotation of the stress regime can play a significant role in the observed long-term deepening of the microseismicity below active injection wells.

58 GEOSCIENCES↗
DOE OSTI2022

Effects of edge-seal design on the mechanical and thermal performance of vacuum-insulated glazing

Although vacuum-insulated glazing (VIG) has been proposed as a promising solution towards developing energy-efficient buildings, VIGs have not become popular in the market due to several technical challenges including the complexity of the fabrication process. In particular, the edge-seal is a key component that significantly affects the thermal insulation and mechanical performance, and the development of edge-seal with adequate thermal insulation, mechanical strength, and reasonable processing cost is essential to overcome such technical issues in VIG. For this purpose, effects of edge-seal design parameters on the VIG performance should be identified. In this research, we analyzed the edge-seal for thermal transport as well as structural stresses to study the effects, and then identified and evaluated the material mixes for the edge-seal requirements. Here, the finite element simulations showed the significance of VIG corner calculation on overall thermal transmittance and the importance of seal conductivity below 1 W/m.K. The experiments with the flexible seals with different ratios of fine glass powder demonstrated that the measured shear strength values for the seal with less than 30% glass powder were more than 10 times larger than the calculated shear stress values. Based on these simulation and experimental results, a flexible sealant was developed using a proprietary mix of ceramic materials that meets the requirements of the designed VIG edge-seal, including structural as well as thermal stress resistance and a low conductivity. Moreover, the sealant is self-curing under atmospheric conditions, and thus it does not require costly inline process of laser curing or oven baking.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Complete Status Report Documenting the Development of Friction Stir Welding for Producing a Butt Joint in Thin Wall Tubing of ODS Alloys

A challenging attempt using friction stir welding to produce a butt joint between two 1.0 mm thick plates of 14YWT was nearly successful. The modified tool pin mostly penetrated the thickness of the adjoined thin plate samples resulting in a small blunt notch on the bottom surface of the butt joint. A bend occurred in the thin plate sample due to formation of stresses from clamping the thin plate sample during FSW. The peak displacement of the bend occurred along the center line of the butt joint. The weld zone showed good microstructure mixing with some flashing mostly on the retreating side that suggested excess heat input by friction between the traveling and spinning tool pin and thin plates of 14YWT. Neutron diffraction experiments were performed on the bent sample of 14YWT to measure the residual stresses in the butt joint on the HB-2B beam line at the Neutron Residual Stress Facility at HFIR. However, the neutron scattering to background noise ratio was very low during data acquisition. Measurements of lattice strains were obtained in the rolling (RD) and normal (ND) directions but not in the extrusion direction (ED) of the thin plate sample. Unfortunately, the quality of the data was not sufficient for determining the residual stresses in the butt joint. From the previous neutron scattering experiments conducted on four hot rolled plates (R1, R2, R3 and R4) of 14YWT, the residual stresses were calculated from lattice strain data and plotted in 2-D color contour maps for ED, RD and ND. The results showed variations in the magnitude and spatial position of residual stresses in the microstructure of the four plates. The magnitude in residual stresses in ED and RD increased gradually with increasing rolling deformation for R1 (50%) to R2 (76%), followed by larger increases in R3 (89%) and R4 (95%). The maximum peaks in residual stresses for R4 was ~710-730 MPa in ED and RD. The residual stresses in ND were below ~300 MPa for each of the 4 rolled plates, with a maximum residual stress in R4 (95%) of ~400 MPa. These results indicated that significant residual stresses representing ~57% of the yield stress of 14YWT (SM13) had accumulated in the 14YWT plate after hot rolling to the 1 mm thickness, which was used in the FSW experiments. The microstructure analysis of the 1.0 mm thick plate of 14YWT containing the butt joint produced by FSW revealed significant changes in the grain structure and texture between the unaffected zone (UZ) of the rolled plate and the butt joint weld. In the UZ, the grains were severely elongated after 95% deformation and mainly possessed the {101} texture component in the ED. Two predominant grain misorientation angles were observed: one consisting of small misorientations of <5º between highly elongated grains and the other consisting of large misorientations of >15º in pockets of small grains that might have formed by dynamic recrystallization processes during rolling of the plate at 1000ºC. The morphology and texture of grains in the heat affected zone (HAZ) were similar to those in UZ, except the grains had coarsened due to excess heat input during FSZ. A population of pores formed and the grains both coarsened and changed to isotropic morphology in the microstructure of the butt joint on both the advancing and retreating sides. The texture of grains in both regions were randomly orientated between the 3 poles {001}, {101} and {111}. The main difference in texture of grains between the advancing side and retreating side was with the misorientation angles of neighboring grains. Most of the grains in the advancing side had misorientation angles of >5º indicating that these grains experienced recrystallization. On the retreating side, areas containing larger grains showed high misorientation angles of >15º compared to areas appearing as bands of grains with misorientation angles of <5º and between 5º-15º. Numerous small pockets of grains with low misorientation angles of <5º were distributed in the microstructure among the grains showing high misorientation angles of >15º. It was determined that the grains associated with low misorientation angles contained sub-grain structures, which was consistent with both recrystallization and recovery mechanisms for grains located on the retreating side of the butt joint. These results demonstrate the complex microstructure that evolves in the microstructure of 14YWT resulting from the input of high temperatures and severe deformations during FSW.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Probing the properties and mechanisms of failure waves in soda-lime glass

Soda-lime glass (SLG) and other silica glasses exhibit the failure wave phenomenon under shock compression. The mechanism responsible for this peculiar behavior of glasses is still unresolved. In this study, a series of plate impact experiments was performed at three different impact stresses of 6.4, 8.3, and 10.8 GPa to better understand the mechanisms underlying the failure wave phenomenon. Specifically, spall experiments were conducted to probe the speed and existence of failure waves at different stresses in SLG. A layered glass target was used to probe the possibility of a “renucleation” of the wave at the SLG–SLG interface. When it existed, the failure wave was inferred to propagate at a speed of 1.3 km/s. However, it was observed that the failure wave phenomenon ceases to exist for impact stresses higher than 10 GPa. In experiments with a 6.4 GPa impact stress, the peak free surface velocity was significantly less than what is predicted by stress-Hugoniot calculations. This velocity deficit and other important features of the measured free surface velocity profiles were simulated using finite element analysis by incorporating an abrupt densification of SLG at a critical stress in the equation of state. Furthermore, this densification feature is similar to what would be expected of a phase transition. Although unable to unambiguously reveal the mechanism causing the failure wave phenomenon, the results of the present work clearly indicate that the failure wave causes a secondary compression and densification in SLG.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Prediction of continuous porosity gradients in ceramics using ZnO as a model material

Materials with gradient microstructures have a wide range of applications such as cutting tools, armor, and electronic devices. However, it is difficult to predict and control the gradient microstructure during processing. In the present work, a continuous porosity gradient was successfully achieved in ZnO material via spark plasma sintering with a large induced thermal gradient. The porosity is overestimated if isothermal prediction is applied. The current work proposed a more accurate prediction of the porosity by considering the stress-shielding effect caused by the thermal gradient. The shielding effect results from different stress states in the sample due to differential sintering: the hotter side of the specimen experiences a higher strain rate and more shrinkage while the colder side experiences a lower strain rate and less shrinkage simultaneously. Therefore, the axial strains are varied throughout the sample thickness. Using the constituent equations in advanced sintering analysis, the shield stress was calculated to be approximately 13 MPa for the viscoelastic assumption of sintering. To improve the accuracy of predicting porosity gradient, it is necessary to add a load to overcome the shield stress when the materials are sintered with a thermal gradient.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Stress Relaxation Related to Spontaneous Thin Film Buckling: Correlation between Finite Element Calculations and Micro Diffraction Analysis

Compressive residual stresses generated during thin film deposition may lead to undesirable film damage, such as delamination, buckling, and flaking, ultimately leading to the failure of the device employing the film. Understanding the residual stress generation and role in these damage mechanisms is necessary to preserve thin film integrity and optimize its functional properties. Thin shell theory has been used for decades to predict buckling but the results have not yet been correlated with experimental data since the techniques used to measure stress in metallic films were not able to do so at the required micron scale until recently. Micro scanning X-ray diffraction now enables the direct mapping of the local stress of metallic films. In this paper, finite element method based on thin shell theory and synchrotron X-ray micro diffraction have been used to determine stress maps of thin film buckling patterns. Calculations of the stress distribution in the metallic films have been performed taking into account the buckling geometry determined from optical measurements. Stress distributions over gold blisters and tungsten wrinkles obtained with the two techniques are in fair agreement and allow for the accurate determination of the stress relaxation profile from the bottom to the top of the buckling, validating the thin shell theory model.

36 MATERIALS SCIENCE↗