Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRESS CALCULATIONS”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Validation of Hole-Drilling Residual Stress Measurements in Workpieces of Various Thickness

A recent revision to the ASTM E837 standard for near-surface residual stress measurement by the hole-drilling method describes a new thickness-dependent stress calculation procedure applicable to “thin” and “intermediate” workpieces for which strain versus depth response depends on workpiece thickness. This new calculation procedure differs from that of the prior standard, which applies only to thick workpieces with strain versus depth response independent of thickness. Herein we assess the new calculation procedures by performing hole-drilling residual stress measurements in samples with a range of thickness. Near-surface residual stress is measured in a thick aluminum plate containing near-surface residual stress from a uniform shot peening treatment, and in samples of different thickness removed from the plate at the peened surface. A finite element (FE) model is used to assess consistency between measured residual stress across the range of sample thickness. Measured residual stress varies with sample thickness, with thinner samples exhibiting smaller near-surface compressive stress and a larger gradient of subsurface stress. These trends are consistent with both observed bending (curvature) of the removed samples and the trend in FE-calculated expected residual stress. The measured and expected residual stresses are in good agreement for samples of intermediate thickness, but the agreement decreases with sample thickness. Measured residual stress is invariant with gage circle diameter. The new thickness-dependent stress calculation procedure for hole-drilling provides meaningful improvement compared to thick-workpiece calculations.

42 ENGINEERING↗
DOE OSTI2025

FY25 status report on the addition of candidate materials in Class B Code Case

This report provides the time-dependent allowable stress calculation strategy leveraging the limited creep rupture tests data generated to support the allowable stress for 100,000 hours in American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC), Section II, Part D. A variable confidence index procedure to extrapolate material properties to support 500,000 hours design life is discussed. Time-dependent allowable stresses for Class B component design and analysis are presented for Grade 1 and Grade 2 of Alloy 625. The presented data extrapolation and allowable stress calculation method will support new material addition using limited creep rupture data in the new ASME Boiler and BPVC, Section III, Division 5, Class B rules.

Part D↗
DOE OSTI2025

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗
DOE OSTI2026

Uncovering hidden bias in neutron diffraction residual strain measurements

When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation-dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e. noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (∼200 µm) of plastic strain, texture and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor and thirty-six measurement least-squares estimation. For the last two cases, the appropriateness of the uncertainty interval is statistically evaluated on the basis of a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that the interpretation of uncertainty in residual strain is strongly dependent on intrinsic sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pretest Modeling A Spent Nuclear Fuel Seismic Shake Test

The U.S. Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a dry canister storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental U.S. to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. This paper describes the pretest predictions made with LS-DYNA models of the mockup storage systems. The test will use two dry storage system configurations, a mockup vertical concrete cask (VCC) and a mockup horizontal storage module (HSM). The test will use a production-quality canister and basket. Within the canister will be four instrumented fuel assemblies with fuel rods containing surrogate mass and 28 instrumented dummy assemblies that are intended to match the mass and outer dimension of a fuel assembly. The finite element models include models of the VCC and HSM on the shake table to calculate the system level dynamic responses and separate single fuel assembly models to calculate stress and strain in fuel assembly components. Both types of models include nonlinear behavior like rod-to-rod contact and the ability for VCC’s to rock and slide. This paper presents the expected response of the VCC, HSM, and fuel assemblies to the shake table testing that is planned to start in April of 2024. The earthquake conditions represent seismic hazards in the 2,000-to-20,000-year return period range. The test data is expected to confirm the expectation that fuel rod cladding will remain intact, fuel assembly structural components like guide tubes will remain intact, and no significant VCC sliding or tipping will occur in the range of conditions to be tested with the shake table.

Klymyshyn, Nicholas A.↗
DOE OSTI2024

Stress field measurements using quantitative schlieren

Quantitative schlieren analysis is extended here to optically transparent solids in quasi-static and dynamic experiments to measure stress distributions. The quasi-static experiments in polymethyl methacrylate (PMMA) compared refraction angles and stress gradients calculated from schlieren images to the analytical Flamant solution of a line load on a half-space. The quantitative schlieren measurements of the stress field in the thin sample with a load compared well to the analytical solution. The analysis method was then extended to explosive induced shock waves in PMMA. The explosive induced response of PMMA was experimentally studied using high-speed schlieren to visualize the shock propagation in conjunction with Photon Doppler Velocimetry (PDV) to record surface velocity histories. The stress state estimated from the schlieren images was compared to the stress calculated from the PDV measurements. High-speed imaging limitations caused the shock wave to not be fully resolved in the images, but was resolved in the PDV measurement. The stress state behind the shock calculated from the high-speed images followed a similar trend to the stress calculated from the PDV measurements.

Physics↗
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 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 OSTI2025

CALPHAD-based ICME design of single-step aging to enhance mechanical strength of WAAM Haynes 282

To match the strength of wire-arc additive manufactured Haynes 282 to its wrought counterpart via a single-step aging heat treatment, the CALPHAD (Calculation of Phase Diagrams) method is integrated with physics-based process-structure-property models and experimental validation. The integrated computational materials engineering (ICME) framework simulates the effects of aging on γ′ and M 23 C 6 precipitation and the resulting yield strength. To improve simulation reliability, the interfacial energies between γ/γ′ and γ/M 23 C 6 carbides were estimated by comparison with precipitation kinetic modeling and measured precipitate sizes. γ′ and M23C6 were found to precipitate simultaneously between 640 and 860 °C, producing microstructures similar to those produced by two-step aging. The optimal γ′ size for peak yield stress was calculated to be 20–23 nm. WAAM Haynes 282 aged at 780 °C for 50 h exceeded the mechanical performance of its wrought counterpart subjected to two-step aging, though desired properties can also be achieved at 800 °C for 16 h or less. The error in yield strength is less than 20 MPa, demonstrating good agreement between the modeling framework and experiments. Creep studies showed that WAAM Haynes 282 exceeded the calculated rupture time, reaching 481 h. This proposed methodology can accelerate the design of aging heat treatments for any γ′-strengthened nickel-base alloy, minimizing the resources required for trial-and-error experiments.

CALPHAD↗
DOE OSTI2024

Interim Development of new Class B Code Case with Variable Design Lifetime and enhanced design rules to guard against cyclic structural failure mode

This report summarizes updates on the ongoing development of the new American Society of Mechanical Engineers Section III, Division 5, Class B rules to address the potential failure modes of Class B components. An initial design-by-analysis approach and supporting design rules are presented. The elastic-perfectly plastic (EPP) analysis approach has been adopted for primary stress limit and ratcheting check to prevent the component failure against the structural failure mode against primary load and strain accumulation due to cyclic loads. The creep-fatigue damage assessment uses explicitly defined elastic follow-up calculated from the stress concentration region of given component. Damage fraction calculation uses a novel coupled approach to capture the influence of the elastic follow-up and provide adequate conservatism for Class B components. The proposed Class B rules does not use the stress classification approach and uses the combined loads to assess component.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Study of hydride reorientation and ductility of SRA ZIRLO and PRXA ZIRLO using ring compression testing (RCT)

Before dry storage, spent nuclear fuel is dried at high temperatures, which induces hydride dissolution. Radial hydride particles precipitate when cooled under a stress above a threshold stress. This process can severely degrade the ductility of cladding tubes. Zirconium alloys with different microstructures, such as grain boundary density, were quantified and studied to examine their effect on hydride reorientation. The ring compression loading mode creates a range of stress states depending on the location within the ring with respect to the applied load, producing various hydride morphologies that can be used to study the effects of microstructures on hydride reorientation. Finite element modeling was performed to calculate the hoop stress, which is the dominant stress, during hydride reorientation. These results were coupled with hydride morphology to determine the threshold stress for hydride reorientation, depending on the microstructure. In conclusion, the samples were then mechanically tested using ring compression to evaluate the ductility of zirconium alloys as a function of different hydride morphologies.

Hydride embrittlement↗
DOE OSTI2025

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↗
DOE OSTI2026

Influence of Numerical Modeling Approaches on Damped Behavior of Flexible Beams: Preprint

Composites structures are widely used in aerospace and wind energy applications for their excellent stiffness and strength-to-weight properties. In these structures, structural damping is critical to predict vibration amplitudes, performance, and reliability. Structural damping is of particular interest for slender wings, rotorcraft blades, and wind turbine blades that can exhibit complex vibration phenomena and are frequently modeled with geometrically exact beam theory (GEBT). Standard approaches of stiffness proportional or modal damping merely assign user defined values and cannot predict damping behavior. This work compares stiffness proportional damping to two more advanced damping approaches: modal strain energy and Prony series. The modal strain energy approach uses a sectional analysis tool to calculate the beam stiffness and postprocess internal stresses from GEBT simulations. The internal stresses are then used to calculate modal damping factors. The Prony series is implemented within GEBT to directly model viscoelastic behavior of the composites. These approaches are compared by modeling the evolution of the damping factors of a realistic flexible wind turbine blade with varying rotational speed. Discrepancies between the approaches suggest areas for future modeling development, but differences in nonlinear damping values are less than current uncertainties about the magnitude of structural damping.

17 WIND ENERGY↗
DOE OSTI2025

Dynamic Behavior of Oval-Twisted Helical Tube Heat Exchanger: Numerical Study with RELAP5-3D

Convective heat transfer characteristics and theoretical thermal stress behaviors are numerically calculated using RELAP5-3D for the helical-coiled once-through steam generator (H-OTSG) and the novel heat exchanger design known as the oval-twisted helically coiled heat exchanger (OTHCHX) under (1) fluctuating wall temperature conditions, (2) square-wave pulsating flow conditions, and (3) the combined effects of fluctuating wall temperature and square-wave pulsating flow conditions. Heat transfer coefficient models for the H-OTSG and OTHCHX were developed based on existing data and implemented into RELAP5-3D, successfully capturing the N⁢uavg behavior within 8% to 10% of the reported data. Under fluctuating wall temperature conditions, the OTHCHX displayed higher N⁢u avg behavior than the H-OTSG. As 𝑓 increased, the $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ decreased. The $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ was higher for the OTHCHX than for the H-OTSG under fluctuating wall temperature conditions. Under pulsating flow conditions, the H-OTSG and OTHCHX displayed much higher 𝑁⁢𝑢 𝑎𝑣𝑔 than under constant flow conditions. The H-OTSG displayed a higher $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ over the OTHCHX. Under combined fluctuating wall temperature and pulsating flow conditions, the augmented heat transfer behavior from the pulsating flow was counteracted by the wall temperature fluctuations, producing slightly higher 𝑁⁢𝑢 𝑎𝑣𝑔 over constant wall temperature, constant flow conditions, but much lower than only constant pulsating flow under constant wall temperature conditions. The effects of simultaneous wall temperature fluctuations and square-wave pulsating flow caused higher $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ than that of only wall temperature fluctuations or pulsating flow. As the Reynolds number (Re) increased, $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ increased. However, when 𝑓=𝑓$_{\dot{m}}$, the $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ showed decreasing values as Re increased. In conclusion, the results indicate that thermal-fluid resonance can help mitigate thermal stresses.

Thermal stress↗
DOE OSTI2025

Report on the Integration of Experimental and Modeling Data for Initial Equivalence Study of Mechanical Performance in Irradiated LPBF 316 Stainless Steel

To ensure the rapid development, deployment, and use of advanced nuclear technologies, faster qualification approaches are needed. Typically, the primary pathway uses traditional data packages consistent with the American Society of Mechanical Engineers Boiler and Pressure Vessel Code, which does not consider the environmental effects the material will experience such as corrosion and radiation damage. Examining radiation effects requires a significant amount of space in US facilities at the Advanced Test Reactor and the High Flux Isotope Reactor (HFIR) and suffers from natural gradients in temperature and neutron flux profiles. Ion irradiation may enable rapid assessment of radiation-induced damage to a material and is proposed as part of an accelerated materials qualification framework through the Advanced Materials and Manufacturing Technologies program. To enhance the utility of ion irradiation as an examination tool, this report provides the initial assessment of engineering-relevant properties of microstructures produced from ion irradiation in the near-surface volume. Nanoindentation, Vickers hardness, and known tensile properties were brought together with simple mathematical models and experimental data for an initial equivalence study of the mechanical performance of irradiated laser powder bed fusion (LPBF) 316 stainless steels across length scales. Direct observation of the calculated ion irradiation yield stress and measured neutron irradiation yield stress at 2 dpa showed that both datasets exhibit the same trend with irradiation temperature and overlap within an acceptable band of stress values. Ion irradiations at 10 dpa serve as a prediction of properties to compare to postirradiation examination of HFIR-irradiated LPBF 316H further in the program. This work is a significant demonstration of the Licensing Approach with Ions and Neutrons, which uses ion irradiations to generate mechanical property information more rapidly than through neutron irradiations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗