Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local stresses”

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 541 records · Page 30

Convolutional Neural Network–Aided Temperature Field Reconstruction: An Innovative Method for Advanced Reactor Monitoring

In this study, the capabilities of a physics-informed convolutional neural network (CNN) for reconstructing the temperature field from a limited set of measurements taken at the boundaries of internal flows are demonstrated. Such an approach enables the development of less invasive monitoring methods for real-time plant diagnostics. As a test case, a Molten Salt Fast Reactor (MSFR) design was selected. This circulating fuel reactor has received interest from both scientific and industrial communities due to its intrinsic safety and sustainability. Molten salt flows in such reactors, however, can present highly localized temperature peaks that can induce significant thermal stresses onto the vessel walls. At these local maxima, the salt temperature may exceed a thousand kelvins, which makes a direct measurement challenging or even unfeasible. The proposed CNN algorithm allows one to detect indirectly such discontinuities through an accurate, albeit indirect, temperature measurement method during reactor operation. The datasets employed to train and test the machine learning models in the present work were generated with Nek5000, a computational fluid dynamics (CFD) code developed at Argonne National Laboratory. The CNN algorithm is trained with CFD results that span a set of MSFR operational power and flow ranges. Here, to demonstrate the efficacy of the algorithm, predictions are made for test cases contained within the training range but for which the CFD data were not used when training. Results demonstrate that the proposed technique properly characterizes temperature peaks and distributions within the domain for a broad range of scenarios.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Overview of Boundary Layer Transition Research in Support of Orbiter Return To Flight

A predictive tool for estimating the onset of boundary layer transition resulting from damage to and/or repair of the thermal protection system was developed in support of Shuttle Return to Flight. The boundary layer transition tool is part of a suite of tools that analyze the aerothermodynamic environment to the local thermal protection system to allow informed disposition of damage for making recommendations to fly as is or to repair. Using mission specific trajectory information and details of each damage site or repair, the expected time (and thus Mach number) at transition onset is predicted to help define the aerothermodynamic environment to use in the subsequent thermal and stress analysis of the local thermal protection system and structure. The boundary layer transition criteria utilized for the tool was developed from ground-based measurements to account for the effect of both protuberances and cavities and has been calibrated against select flight data. Computed local boundary layer edge conditions were used to correlate the results, specifically the momentum thickness Reynolds number over the edge Mach number and the boundary layer thickness. For the initial Return to Flight mission, STS-114, empirical curve coefficients of 27, 100, and 900 were selected to predict transition onset for protuberances based on height, and cavities based on depth and length, respectively.

Berry, Scott A.↗

Implementation of thermal residual stresses in the analysis of fiber bridged matrix crack growth in titanium matrix composites

In this research, thermal residual stresses were incorporated in an analysis of fiber-bridged matrix cracks in unidirectional and cross-ply titanium matrix composites (TMC) containing center holes or center notches. Two TMC were investigated, namely, SCS-6/Timelal-21S laminates. Experimentally, matrix crack initiation and growth were monitored during tension-tension fatigue tests conducted at room temperature and at an elevated temperature of 200 C. Analytically, thermal residual stresses were included in a fiber bridging (FB) model. The local R-ratio and stress-intensity factor in the matrix due to thermal and mechanical loadings were calculated and used to evaluate the matrix crack growth behavior in the two materials studied. The frictional shear stress term, tau, assumed in this model was used as a curve-fitting parameter to matrix crack growth data. The scatter band in the values of tau used to fit the matrix crack growth data was significantly reduced when thermal residual stresses were included in the fiber bridging analysis. For a given material system, lay-up and temperature, a single value of tau was sufficient to analyze the crack growth data. It was revealed in this study that thermal residual stresses are an important factor overlooked in the original FB models.

Bakuckas, John G., Jr.↗

Approximate Analysis for Interlaminar Stresses in Composite Structures with Thickness Discontinuities

An efficient, approximate analysis for calculating complete three-dimensional stress fields near regions of geometric discontinuities in laminated composite structures is presented. An approximate three-dimensional local analysis is used to determine the detailed local response due to far-field stresses obtained from a global two-dimensional analysis. The stress results from the global analysis are used as traction boundary conditions for the local analysis. A generalized plane deformation assumption is made in the local analysis to reduce the solution domain to two dimensions. This assumption allows out-of-plane deformation to occur. The local analysis is based on the principle of minimum complementary energy and uses statically admissible stress functions that have an assumed through-the-thickness distribution. Examples are presented to illustrate the accuracy and computational efficiency of the local analysis. Comparisons of the results of the present local analysis with the corresponding results obtained from a finite element analysis and from an elasticity solution are presented. These results indicate that the present local analysis predicts the stress field accurately. Computer execution-times are also presented. The demonstrated accuracy and computational efficiency of the analysis make it well suited for parametric and design studies.

Rose, Cheryl A.↗

Load biaxiality and fracture - Synthesis and summary

The question of load biaxiality strongly affecting almost all of the characteristics of brittle fracture behavior of a cracked body is investigated. Results of research published piecemeal over the years by the authors have shown that the stress and displacement fields, the elastic strain energy density, and the maximum shear stress near the crack tip are all altered by loads applied parallel to the crack, as are the angle of initial crack extension, the strain energy of the entire body, the fracture load, and the rate of fatigue crack growth. A synthesis and summary of this research is presented showing the shortcomings inherent in both Griffith's global energy rate theory for crack instability and Irwin's local crack-tip stress intensity theory for fracture toughness. In both theories only the tensile load perpendicular to the crack influences fracture behavior of the body.

Eftis, J.↗

Wave-driven torques to drive current and rotation

In the classic Landau damping initial value problem, where a planar electrostatic wave transfers energy and momentum to resonant electrons, a recoil reaction occurs in the nonresonant particles to ensure momentum conservation. To explain how net current can be driven in spite of this conservation, the literature often appeals to mechanisms that transfer this nonresonant recoil momentum to ions, which carry negligible current. However, this explanation does not allow the transport of net charge across magnetic field lines, precluding E X B rotation drive. In this study, we show that in steady state, this picture of current drive is incomplete. Using a simple Fresnel model of the plasma, we show that for lower hybrid waves, the electromagnetic energy flux (Poynting vector) and momentum flux (Maxwell stress tensor) associated with the evanescent vacuum wave become the Minkowski energy flux and momentum flux in the plasma and are ultimately transferred to resonant particles. Thus, the torque delivered to the resonant particles is ultimately supplied by the electromagnetic torque from the antenna, allowing the nonresonant recoil response to vanish and rotation to be driven. We present a warm fluid model that explains how this momentum conservation works out locally, via a Reynolds stress that does not appear in the one-dimensional initial value problem. This model is the simplest that can capture both the nonresonant recoil reaction in the initial-value problem, and the absence of a nonresonant recoil in the steady-state boundary value problem, thus forbidding rotation drive in the former while allowing it in the latter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A scheme for computing surface layer turbulent fluxes from mean flow surface observations

A physical model and computational scheme are developed for generating turbulent surface stress, sensible heat flux and humidity flux from mean velocity, temperature and humidity at some fixed height in the atmospheric surface layer, where conditions at this reference level are presumed known from observations or the evolving state of a numerical atmospheric circulation model. The method is based on coupling the Monin-Obukov surface layer similarity profiles which include buoyant stability effects on mean velocity, temperature and humidity to a force-restore formulation for the evolution of surface soil temperature to yield the local values of shear stress, heat flux and surface temperature. A self-contained formulation is presented including parameterizations for solar and infrared radiant fluxes at the surface. Additional parameters needed to implement the scheme are the thermal heat capacity of the soil per unit surface area, surface aerodynamic roughness, latitude, solar declination, surface albedo, surface emissivity and atmospheric transmissivity to solar radiation.

Hoffert, M. I.↗

A comparison of three algebraic stress closures for combustor flow calculations

A comparison is made of the performance of two locally nonequilibrium and one equilibrium algebraic stress closures in calculating combustor flows. Effects of four different pressure-strain models on these closure models are also analyzed. The results show that the pressure-strain models have a much greater influence on the calculated mean velocity and turbulence field than the algebraic stress closures, and that the best mean strain model for the pressure-strain terms is that proposed by Launder, Reece and Rodi (1975). However, the equilibrium algebraic stress closure with the Rotta return-to-isotropy model (1951) for the pressure-strain terms gives as good a correlation with measurements as when the Launder et al. mean strain model is included in the pressure-strain model. Finally, comparison of the calculations with the standard k-epsilon closure results show that the algebraic stress closures are better suited for simple turbulent flow calculations.

Nikjooy, M.↗

Comparison of three state-of-the-art crystal plasticity based deformation twinning models for magnesium alloys

In magnesium alloys, deformation twinning and its interactions with dislocation slip are responsible for a sigmoidal shape stress–strain behavior and an asymmetrical tension–compression yield strength in magnesium alloys. The sensitivity of twinning to the underlying microstructure renders the crystal plasticity method the most commonly adopted modeling approach for magnesium-twinning. This paper compares three state-of-the-art crystal plasticity-based twinning models from the literature, namely the elastic-viscoplastic self-consistent twinning-detwinning (EVPSC-TDT) model, crystal plasticity finite element model based on enhanced predominate twin reorientation approach (CPFE-ePTR), and the crystal plasticity finite element model based on “discrete twinning” approach (CPFE-DT). A polycrystalline microstructure is simulated with all three methods to compare the resulting stress–strain curves and lattice strains to those from the experimentally measured data. All three methods showed the capability of capturing the experimental results with different levels of accuracy. Additionally, the EVPSC-TDT method avoids solving the finite element matrices and showed the highest computational efficiency. The CPFE-ePTR model shows a higher accuracy in capturing the lattice strain. The CPFE-DT relies on high-resolution finite element mesh and is much slower than the other two methods, but it captured the local deformation concentration and stress reversal phenomena near the twin band, which was not possible with the other two methods. Based on the comparison, guidance for the selection of the appropriate model based on the specific modeling target is provided in this paper.

36 MATERIALS SCIENCE↗

Synergistic effect of microstructure and defects on the initiation of fatigue cracks in additively manufactured Inconel 718

Fatigue cracks in additively manufactured (AMed) Inconel 718 (IN-718) in machined surface condition often initiate from persistent slip bands (PSBs) unlike other popular AM alloys such as 17-4 PH stainless steel or Ti-6Al-4V, where fatigue crack initiation is exclusively from volumetric defects; therefore, a competition between PSB- vs. defect- mediated crack initiation clearly exists. To shed light on the factors governing the competition, this study investigates the characteristics of cyclic strain localization, PSB formation, and crack initiation via crystal plasticity (CP) modeling of cyclic loading on polycrystalline aggregates which are then validated by experiments. In this work, a physics-based, free slip distance (FSD) dependent slip strength evolution law is proposed, which is shown to enable the CP model to simulate the heterogeneous strain distribution in IN-718. Implementing a crack initiation criterion based on strain contrast, the locations and lives for crack initiation can also be calculated. It is shown that both FSD and resolve shear stress influence the strain localization and crack initiation behaviors. The distribution of a localization parameter calculated based on the multiplication of FSD and Schmid factor within a grain is found to correlate well with the locations of PSBs. The maximum values of the localization parameters within a microstructure are shown to correlate well with the experimentally obtained crack initiation lives. The presence of volumetric defects in IN-718 generally do not impact strain localization behavior unless their size is large compared to the grain size.

42 ENGINEERING↗

Theoretical studies on the mechanical behavior of granular materials under very low intergranular stresses

The salient aspects of the theoretical modeling of a conventional triaxial test (CTC) of a cohesionless granular medium with stress and strain rate loading are described. Included are a controllable gravitational body force and provision for low confining pressure and/or very low intergranular stress. The modeling includes rational, analytic, and numerical phases, all in various stages of development. The numerical evolutions of theoretical models will be used in final design stages and in the analysis of the experimental data. In this the experimental design stage, it is of special interest to include in the candidate considerations every anomaly found in preliminary terrestrial experimentation. Most of the anomalies will be eliminated by design or enhanced for measurement as the project progresses. The main aspect of design being not the physical apparatus but the type and trajectories of loading elected. The major considerations that have been treated are: appearance and growth of local surface aberrations, stress-power coefficients, strain types, optical strain, radial bead migration, and measures of rotation for the proper stress flux.

French, Kenneth W., Jr.↗

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗

In situ high-temperature 3D imaging of the damage evolution in a SiC nuclear fuel cladding material

Silicon carbide (SiC)-based nuclear fission fuel rod cladding has been considered as one of the possible designs for accident tolerant fuels. It is in the form of a SiC fibre reinforced SiC matrix composite tube (SiC f -SiC m ) with monolithic SiC outer and/or inner coating layers. This study focuses on the deformation and fracture processes in this material using in situ X-ray micro-computed tomography (XCT) at room temperature (RT) and 1200 °C in an inert gas environment in a C-ring compression loading configuration. Prior to testing, local properties and residual stresses were characterised using nanoindentation and Raman spectroscopy since they can impact the mechanical behaviour of the material. The 3D strain distribution, crack formation and propagation processes including the toughening mechanisms (e.g., crack deflection, micro-cracking, crack bridging and bifurcation) are investigated in the coating and underlying composites at RT and 1200 °C. There is no particular sequence which toughening mechanism occurs first – this is very different from the conventional toughening theory in ceramic-matrix composites under uniaxial tension loading. Indeed, no evidence of fibre pull-out or fibre fracture was observed in this SiC f -SiC m nuclear cladding material in the current C-ring compression configuration. The correlation between local measurements and bulk mechanical behaviour are discussed.

3D imaging↗

Two-Dimensional Nonlinear Finite Element Analysis of CMC Microstructures

A research program has been developed to quantify the effects of the microstructure of a woven ceramic matrix composite and its variability on the effective properties and response of the material. In order to characterize and quantify the variations in the microstructure of a five harness satin weave, chemical vapor infiltrated (CVI) SiC/SiC composite material, specimens were serially sectioned and polished to capture images that detailed the fiber tows, matrix, and porosity. Open source quantitative image analysis tools were then used to isolate the constituents, from which two dimensional finite element models were generated which approximated the actual specimen section geometry. A simplified elastic-plastic model, wherein all stress above yield is redistributed to lower stress regions, is used to approximate the progressive damage behavior for each of the composite constituents. Finite element analyses under in-plane tensile loading were performed to examine how the variability in the local microstructure affected the macroscopic stress-strain response of the material as well as the local initiation and progression of damage. The macroscopic stress-strain response appeared to be minimally affected by the variation in local microstructure, but the locations where damage initiated and propagated appeared to be linked to specific aspects of the local microstructure.

Mital, Subodh K.↗

Determining Stress Orientation in Rock Valley, Nevada, Using Ambient Seismic Noise

The stress field and the mechanical properties of rocks are important to consider for nuclear explosion monitoring due to their effect on seismic wave radiation from earthquakes and explosions. At the Rock Valley Direct Comparison site, the regional orientation of the maximum horizontal compressive stress (SH max ) is well constrained, but it is unknown whether there are local heterogeneities. Here, I show that stress-induced anisotropy in nonlinear elasticity can be used to estimate the orientation of SH max . Rocks have compliant internal contacts, such as fractures and mineral grain boundaries, that respond to applied strains more strongly than individual mineral crystals. This strain response is asymmetric between compressive and dilatational strains and is affected by anisotropy in the ambient stress field. Traditional seismic velocity measurements are of strain-averaged velocity, which is less sensitive than nonlinear elasticity to the behavior of compliant internal contacts and therefore the stress field and fracture behavior. My results show that the orientation of SH max measured using ambient seismic noise aligns with regional estimates, with some heterogeneity that can also be explained by the limitations of the seismic array. I demonstrate the potential of a passive technique for monitoring the stress field in places that currently lack local measurements.

58 GEOSCIENCES↗

Tuned Reactivity at the Lithium Metal–Argyrodite Solid State Electrolyte Interphase

Thin intermetallic Li 2 Te–LiTe 3 bilayer (0.75 µm) derived from 2D tellurene stabilizes the solid electrolyte interphase (SEI) of lithium metal and argyrodite (LPSCl, Li 6 PS 5 Cl) solid-state electrolyte (SSE). Tellurene is loaded onto a standard battery separator and reacted with lithium through single-pass mechanical rolling or transferred directly to SSE surface by pressing. State-of-the-art electrochemical performance is achieved, e.g., symmetric cell stable for 300 cycles (1800 h) at 1 mA cm -2 and 3 mAh cm -2 (25% DOD, 60 µm foil). Cryo-stage focused ion beam (Cryo-FIB) sectioning and Raman mapping demonstrate that the Li 2 Te–LiTe 3 bilayer impedes SSE decomposition. The unmodified Li–LPSCl interphase is electrochemically unstable with a geometrically heterogeneous reduction decomposition reaction front that extends deep into the SSE. Decomposition drives voiding in Li metal due to its high flux to the reaction front, as well as voiding in the SSE due to the associated volume changes. Analysis of cycled SSE found no evidence for pristine (unreacted) lithium metal filaments/dendrites, implying failure driven by decomposition phases with sufficient electrical conductivity that span electrolyte thickness. In conclusion, DFT calculations clarify thermodynamic stability, interfacial adhesion, and electronic transport properties of interphases, while mesoscale modeling examines interrelations between reaction front heterogeneity (SEI heterogeneity), current distribution, and localized chemo-mechanical stresses.

25 ENERGY STORAGE↗

Multiphysics simulation of recent experiments on alkali‐silica reaction expansion in reinforced concrete members

Alkali‐silica reaction (ASR) is an important degradation process that causes volumetric expansion and damage in concrete, and is affected significantly by the local temperature, moisture and stress conditions that often vary across the regions of a structure. Numerical simulation is essential to predict the progression and effects of ASR on the performance of structures. Because of the interactions between thermal and moisture transport and mechanical deformation, it is important for numerical models to represent all these physical phenomena and the coupling between them. Simulations of ASR in reinforced concrete (RC) structures are further complicated by the need to capture interactions between concrete and embedded reinforcing bars. Here, this paper describes the implementation of a scalable, coupled‐physics ASR model for simulating RC structures and assesses the ability of that model to predict ASR‐induced expansion in recent laboratory tests on RC block and beam specimens. These laboratory tests and the simulation approach were selected because of their applicability to RC structural‐scale simulations. This validation study helps builds confidence the ability of this approach to model ASR expansion in large, complex RC structures, which is a current high‐priority need.

36 MATERIALS SCIENCE↗