Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress Modeling”

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

Initial Developments in Modeling Graphite Behavior

Graphite is used in High Temperature Reactor designs as a moderator, reflector and core material responsible for protecting the fuel and maintaining structural stability. From a safety perspective, it is imperative to be able to predict how graphite will behave under reactor operating conditions which may compromise the core’s structural integrity or fuel’s safety performance. The following report summarizes recent modeling work undertaken at Idaho National Laboratory which focuses on graphite degradation behavior. One phenomenon which is explored in this report and can cause significant degradation to graphite properties is oxidation. Oxidation behavior in graphite is highly temperature dependent. At low temperatures, oxygen can fully penetrate a component and cause a homogeneous damage profile. At high temperatures, only the graphite near the surface of a component is affected. Understanding and modeling this temperature dependence is essential to predicting graphite behavior. In this report, three graphite models are discussed. The first model is used to investigate strength loss after low temperature oxidation. It does this by generating an approximate graphite microstructure then determining the required load to cause crack propagation. It has been shown that the model reproduces the quasi brittle stress versus strain relationship observed in graphite. Strength loss results from the model are shown to match experimental values. The second model investigates strength loss under high temperature oxidation conditions. This model is more applicable to full scale graphite components where the oxidation damage is often inhomogeneous throughout a component. The third model computes stresses in graphite under reactor conditions. Specifically, the model incorporates the effect of creep, irradiation dose, thermal strains, and oxidation in order to predict internal stresses in a component. All of these models are implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open source, parallel finite element framework.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2024

Out of Plane Stresses Using a Beam-Based Elastic Foundation Bonded Joint Model

Efficient, semi-analytical models of bonded joints are utilized for initial joint design and sizing due to their convenience and speed. Many of these models are based on beam and plate theories, which are most applicable when in-plane stresses/strains dominate. While this can be a sensible assumption for obtaining the deformations of the joint, out of plane (transverse) stresses are often of critical importance for failure in bonded joints composed of fiber reinforced composite materials. Since through-thickness strains are neglected in these semi-analytical theories, the present study shows a method to back-calculate the out of plane stresses based on stress balance equations, which use the in-plane stresses and displacements that have already been obtained. The method was implemented within the Joint Element Designer (JED) software, and verification and validation examples have been shown. As is typical for these types of models, corner stress singularities are not captured by the JED model, but the general trends of the stress and deformation fields agreed with a higher fidelity finite element model. As such, this method for obtaining out of plane stresses, when coupled with the JED code, is very useful for the complete stress analysis of bonded composite joints.

composites↗
NASA NTRS1981

An experimental and computational investigation of the transport of Reynolds stress in an axisymmetric swirling boundary layer

A coordinated experimental and computational investigation of the transport and interaction of the individual Reynolds stresses within a turbulent axisymmetric swirling boundary layer flowing over a stationary cylinder is described. The cylinder is instrumented with new directional surface fence skin-friction gages developed for the experiment. The longitudinal and transverse mean velocity profiles are measured with a miniature directional pressure probe. These data provided a standard on which to compare computed results employing a variety of turbulence models. The turbulence models tested range from a simple mixing length model, through a two-equation model, on to the Reynolds stress equation models. The full Reynolds stress equation models are tested with various sets of modeling coefficients for the pressure rate-of-strain correlation. It is shown that while even the simplest model shows the general features of the data, agreement with the experimental data is enhanced through the use of the more complex models. For a better prediction of the data, however, further model improvement is still required.

Higuchi, H.↗
DOE OSTI2020

Evaluation of a Geotechnical Compaction Model to Estimate Sediment Deformation

When surface ruts from vehicle activity occur above archeological deposits, there is a potential for an impact to the archeological record to occur. While impacts to the soil surface are relatively straightforward to describe, sub-surface disturbances are more difficult to quantify. In particular, disturbances to buried archeological materials can be difficult to characterize. Materials can be damaged by surface disturbances, but they can also be impacted even if they stay intact. Interpretation of intact materials relocated within the sediment column can be problematic. Traditional characterization methods (e.g. excavation) can also disturb archaeological materials. Here, we present a geotechnical compaction model that can be used to assist with sub-surface characterization under vehicle ruts formed during off-road travel. The model is based on a Boussinesq analysis of stress increase below a surface load. The change in bulk density resulting from compaction is assumed to be proportional to the stress increase. The model results were compared against pilot scale tests conducted with homogeneous sediment and field scale in-situ tests. The results indicate that the model generally agreed with measurements to within 20%, although only a single soil type was evaluated in these studies, with a narrow range of initial bulk densities and rut depths. While application of the model across a broader matrix of soil types, bulk densities and rut depths would increase confidence in the model, it appears that the model may be suitable for use in archeological field investigations. Potential uses include both compliance and research based field applications.

Fritz, Brad G.↗
NASA NTRS1981

Forced shelf circulation by an alongshore wind band

The forced circulation over a continental shelf generated by the alongshore wind stress is studied within the frictional regime. The model alongshore wind stress has a finite extent in the alongshore direction and oscillates monochromatically, resembling a series of anticyclones traveling across the coastline. Both the shelf-wave response and the localized non-wavelike response exist within the wind band. The numerical results show that the interaction of shelf waves with locally wind-forced response generates a large increase in the phase speed within the wind band. Given an alongshore variation in the alongshore wind stress, a right-bounded phase propagation is possible even in the absence of continental shelf waves. Given a reasonable friction currently accepted for the east coast of the United States, the resonance mechanism at the cutoff frequency may not be important, and lower frequency wind events generate larger amplitude continental shelf circulation. By reducing the friction, energy at cutoff fequencies leaks out of the wind band in both directions effectively. It also is shown that non-wavelike response depends only on the local wind stress and is not affected significantly by friction.

Chao, S.-Y.↗
NASA NTRS1986

Direct simulation of a turbulent boundary layer up to R(sub)(theta)= 1410

The turbulent boundary layer on a flat plate, with zero pressure gradient, is simulated numerically at four stations between R sub theta = 225 and R sub theta = 1410. The three-dimensional time-dependent Navier-Stokes equations are solved using a spectra method with up to about 10 to the 7th power grid points. Periodic spanwise and stream-wise conditions are applied, and a multiple-scale procedure is applied to approximate the slow streamwise growth of the boundary layer. The flow is studied, primarily, from a statistical point of view. The solutions are compared with experimental results. The scaling of the mean and turbulent quantities with Reynolds number is examined and compared with accepted laws, and the significant deviations are documented. The turbulence at the highest Reynolds number is studied in detail. The spectra are compared with various theoretical models. Reynolds-stress budget data are provided for turbulence-model testing.

Spalart, Philippe R.↗
NASA NTRS1988

A comparison of experimental and calculated thin-shell leading-edge buckling due to thermal stresses

High-temperature thin-shell leading-edge buckling test data are analyzed using NASA structural analysis (NASTRAN) as a finite element tool for predicting thermal buckling characteristics. Buckling points are predicted for several combinations of edge boundary conditions. The problem of relating the appropriate plate area to the edge stress distribution and the stress gradient is addressed in terms of analysis assumptions. Local plasticity was found to occur on the specimen analyzed, and this tended to simplify the basic problem since it effectively equalized the stress gradient from loaded edge to loaded edge. The initial loading was found to be difficult to select for the buckling analysis because of the transient nature of thermal stress. Multiple initial model loadings are likely required for complicated thermal stress time histories before a pertinent finite element buckling analysis can be achieved. The basic mode shapes determined from experimentation were correctly identified from computation.

Jenkins, Jerald M.↗
NASA NTRS1988

Direct simulation of a turbulent boundary layer up to R sub theta = 1410

The turbulent boundary layer on a flat plate, with zero pressure gradient, is simulated numerically at four stations between R sub theta = 225 and R sub theta = 1410. The three-dimensional time-dependent Navier-Stokes equations are solved using a spectral method with up to about 10 to the 7th grid points. Periodic spanwise and streamwise conditions are applied, and a multiple-scale procedure is applied to approximate the slow streamwise growth of the boundary layer. The flow is studied, primarily, from a statistical point of view. The solutions are compared with experimental results. The scaling of the mean and turbulent quantities with Reynolds number is compared with accepted laws, and the significant deviations are documented. The turbulence at the highest Reynolds number is studied in detail. The spectra are compared with various theoretical models. Reynolds-stress budget data are provided for turbulence-model testing.

Spalart, Philippe R.↗
NASA NTRS1990

A non-isotropic multiple-scale turbulence model

A newly developed non-isotropic multiple scale turbulence model (MS/ASM) is described for complex flow calculations. This model focuses on the direct modeling of Reynolds stresses and utilizes split-spectrum concepts for modeling multiple scale effects in turbulence. Validation studies on free shear flows, rotating flows and recirculating flows show that the current model perform significantly better than the single scale k-epsilon model. The present model is relatively inexpensive in terms of CPU time which makes it suitable for broad engineering flow applications.

Chen, C. P.↗
NASA NTRS1995

Lagrangian mixed layer modeling of the western equatorial Pacific

Processes that control the upper ocean thermohaline structure in the western equatorial Pacific are examined using a Lagrangian mixed layer model. The one-dimensional bulk mixed layer model of Garwood (1977) is integrated along the trajectories derived from a nonlinear 1 1/2 layer reduced gravity model forced with actual wind fields. The Global Precipitation Climatology Project (GPCP) data are used to estimate surface freshwater fluxes for the mixed layer model. The wind stress data which forced the 1 1/2 layer model are used for the mixed layer model. The model was run for the period 1987-1988. This simple model is able to simulate the isothermal layer below the mixed layer in the western Pacific warm pool and its variation. The subduction mechanism hypothesized by Lukas and Lindstrom (1991) is evident in the model results. During periods of strong South Equatorial Current, the warm and salty mixed layer waters in the central Pacific are subducted below the fresh shallow mixed layer in the western Pacific. However, this subduction mechanism is not evident when upwelling Rossby waves reach the western equatorial Pacific or when a prominent deepening of the mixed layer occurs in the western equatorial Pacific or when a prominent deepening of the mixed layer occurs in the western equatorial Pacific due to episodes of strong wind and light precipitation associated with the El Nino-Southern Oscillation. Comparison of the results between the Lagrangian mixed layer model and a locally forced Eulerian mixed layer model indicated that horizontal advection of salty waters from the central Pacific strongly affects the upper ocean salinity variation in the western Pacific, and that this advection is necessary to maintain the upper ocean thermohaline structure in this region.

Shinoda, Toshiaki↗
NASA NTRS1999

Turbulence Model Comparisons for a High-Speed Aircraft

Four turbulence models are described and evaluated for transonic flows over the High-Speed Research/industry baseline configuration known as Reference H by using the thin-layer, upwind, Navier-Stokes solver known as CFL3D. The turbulence models studied are the equilibrium model of Baldwin-Lomax (B-L) with the Degani-Schiff (D-S) modifications, the one-equation Baldwin-Barth (B-B) model, the one-equation Spalart-Allmaras (S-A) model, and Menter's two-equation Shear Stress Transport (SST) model. The flow conditions, which correspond to tests performed in the National Transonic Facility (NTF) at Langley Research Center, are a Mach number of 0.90 and a Reynolds number of 30 x 10 (exp. 6) based on mean aerodynamic chord for angles of attack of 1 deg., 5 deg., and 10 deg. The effects of grid topology and the representation of the actual wind tunnel model geometry are also investigated. Computed forces and surface pressures compare reasonably well with the experimental data for all four turbulence models.

Rivers, Melissa B.↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
NASA NTRS1985

Lunar Mascon Loading: an Example of the Importance of Fault Initiation Depth on Fault Type Predictions from Stresses on Planetary Surfaces

Geophysical models that predict fault type from stresses calculated at the planetary surface (zero depth) can be misleading because faults generally initiate at depth under a different stress regime and propagate to the surface. The generally accepted and most commonly employed criteria for fault prediction from calculated stresses on planetary surfaces points out problems in neglecting the effects of non-isotropic overburden stresses and the use of realistic material constants and shows a well known example (mascon loading on the moon) where neglecting this effect predicts a different type of fault than is actually observed.

Golombek, M. P.↗
DOE OSTI2024

pvcracks

SAND2024-00922O This software uses electroluminescence images to predict power loss due to cell cracks in photovoltaic modules. The software will incorporate trained variational autoencoder(s) to parameterize cell cracks detected in electroluminescence images of photovoltaic modules and relate cracks to electrical properties; reduced-order models of finite element simulations of electrical behavior of photovoltaic modules under thermomechanical stresses; reduced-order models of stress distributions; inside photovoltaic modules resulting from x-ray topography experiments; and image segmentation and splining methods to analyze x-ray topography measurements of cracked photovoltaic cells. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Hartley, James↗
NASA NTRS2001

Constituent Effects on the Stress-Strain Behavior of Woven Melt-Infiltrated SiC Composites

The stress-strain behavior of 2D woven SiC fiber reinforced, melt-infiltrated SiC matrix composites with BN interphases were studied for composites fabricated with different fiber tow ends per unit length, different composite thickness, and different numbers of plies. In general, the stress-strain behavior, i.e., the 'knee' in the curve and the final slope of the stress-strain curve, was controlled by the volume fraction of fibers. Some of the composites exhibited debonding and sliding in between the interphase and the matrix rather than the more common debonding and sliding interface between the fiber and the interphase. Composites that exhibited this 'outside debonding' interface, in general, had lower elastic moduli and higher ultimate strains as well as longer pull-out lengths compared to the 'inside debonding' interface composites. Stress-strain curves were modeled where matrix crack formation as a function of stress was approximated from the acoustic emission activity and the measured crack density from the failed specimens. Interfacial shear strength measurements from individual fiber push-in tests were in good agreement with the interfacial shear strength values used to model the stress-strain curves.

Morscher, Gregory N.↗
DOE OSTI2024

Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress

Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.

58 GEOSCIENCES↗
DOE OSTI2023

Evolution of eddy viscosity in the wake of a wind turbine

The eddy viscosity hypothesis is a popular method in wind turbine wake modeling for estimating turbulent Reynolds stresses. We document the downstream evolution of eddy viscosity in the wake of a wind turbine from experimental and large-eddy-simulation data. Wake eddy viscosity is isolated from its surroundings by subtracting the inflow profile, and the driving forces are identified in each wake region. Eddy viscosity varies in response to changes in turbine geometry and nacelle misalignment with larger turbines generating stronger velocity gradients and shear stresses. We propose a model for eddy viscosity based on a Rayleigh distribution. Model parameters are obtained from scaling the eddy viscosity hypothesis and demonstrate satisfactory agreement with the reference data. The model is implemented in the curled wake formulation in the FLOw Redirection and Induction in Steady State (FLORIS) framework and assessed through comparisons with the previous formulation. Our approach produced more accurate flow field estimates with lower total error for the majority of cases.

17 WIND ENERGY↗