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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Rough-to-smooth transition of an equilibrium neutral constant stress layer

A model is proposed for low-level atmospheric flows over terrains of changing roughness length, such as those found at the windward end of landing strips adjoining rough terrain. The proposed model is used to develop a prediction technique for calculating transition wind and shear-stress profiles in the region following surface roughness discontinuity. The model for the transition region comprises two layers: a logarithmic layer and a buffer layer. The flow is assumed to be steady, two-dimensional, and incompressible, with neutral hydrostatic stability. A diagram is presented for a typical wind profile in the transition region, obtained from the logarithmic and velocity defect profiles using shear stress calculated by relevant equations.

Logan, E., Jr.↗
NASA NTRS1981

Evaluation of the Doraiswamy-Thompson winter wheat crop calendar model incorporating a modified spring restart sequence

The Robertson phenology was used to provide growth stage information to a wheat stress indicator mode. A stress indicator model demands two acurate predictions from a crop calendar: date of spring growth initiation; and crop calendar stage at growth initiation. Several approaches for restarting the Robertson phenology model at spring growth initiation were studied. Although best results were obtained with a solar thermal unit method, an alternate approach which indicates soil temperature as the controlling parameter for spring growth initiation was selected and tested. The modified model (Doraiswamy-Thompson) is compared to LACIE-Robertson model predictions.

Taylor, T. W.↗
DOE OSTI2023

Scaling Containment on a Large Centrifuge

Containment science concerns the trapping or leakage of radioactive cavity gases from underground nuclear explosions (UNEs). The most physically relevant data for validation of containment science comes from legacy UNE testing or large field-scale chemical explosions. Field-scale tests can inform on the degree of residual stress imparted to host geology and possible formation of an accompanying gas-containing “stress cage”, for example, as a function of chemical explosion yield. However, field testing can be expensive, difficult to conduct, and challenging to instrument with limitations on data coverage (i.e., the number of boreholes; difficult coring through damage zones). Small-scale laboratory tests are typically simpler to conduct and more thoroughly characterize, but may not include the right scale of containment processes for stress-cage or chimney formation. This report presents the theory and fundamentals of large centrifuge physical modeling, including a history of explosive or similar testing relevant to containment science investigations on centrifuges. Subscale models with embedded chemical explosives in the enhanced gravity of a large centrifuge can represent hundreds of meters of depth and large explosions not otherwise attainable in the laboratory as based on scaling of length, energy, and other processes by the g-factor of the centrifuge (i.e., the number of times larger the centrifugal force is than the gravitational force on the Earth’s surface). The centrifuge uniquely joins the simplicity and exhaustiveness of laboratory-scale characterization with the physics of field-scale processes. This report is part of Sandia National Laboratories’ (SNL’s) effort to prepare for physical modeling with approximately one-meter tall geologic models of stress cage and/or chimney formation on its 29-foot [8.84 m] radius centrifuge with a load capacity of 1.6 million 𝑔-pounds [726 𝑔-ton]. The scaling relationships herein will inform model design decisions and performance requirements for in-flight sensors.

42 ENGINEERING↗
NASA NTRS2003

Grid-Independent Large-Eddy Simulation in Turbulent Channel Flow using Three-Dimensional Explicit Filtering

In this paper, turbulence-closure models are evaluated using the 'true' LES approach in turbulent channel flow. The study is an extension of the work presented by Gullbrand (2001), where fourth-order commutative filter functions are applied in three dimensions in a fourth-order finite-difference code. The true LES solution is the grid-independent solution to the filtered governing equations. The solution is obtained by keeping the filter width constant while the computational grid is refined. As the grid is refined, the solution converges towards the true LES solution. The true LES solution will depend on the filter width used, but will be independent of the grid resolution. In traditional LES, because the filter is implicit and directly connected to the grid spacing, the solution converges towards a direct numerical simulation (DNS) as the grid is refined, and not towards the solution of the filtered Navier-Stokes equations. The effect of turbulence-closure models is therefore difficult to determine in traditional LES because, as the grid is refined, more turbulence length scales are resolved and less influence from the models is expected. In contrast, in the true LES formulation, the explicit filter eliminates all scales that are smaller than the filter cutoff, regardless of the grid resolution. This ensures that the resolved length-scales do not vary as the grid resolution is changed. In true LES, the cell size must be smaller than or equal to the cutoff length scale of the filter function. The turbulence-closure models investigated are the dynamic Smagorinsky model (DSM), the dynamic mixed model (DMM), and the dynamic reconstruction model (DRM). These turbulence models were previously studied using two-dimensional explicit filtering in turbulent channel flow by Gullbrand & Chow (2002). The DSM by Germano et al. (1991) is used as the USFS model in all the simulations. This enables evaluation of different reconstruction models for the RSFS stresses. The DMM consists of the scale-similarity model (SSM) by Bardina et al. (1983), which is an RSFS model, in linear combination with the DSM. In the DRM, the RSFS stresses are modeled by using an estimate of the unfiltered velocity in the unclosed term, while the USFS stresses are modeled by the DSM. The DSM and the DMM are two commonly used turbulence-closure models, while the DRM is a more recent model.

Gullbrand, Jessica↗
NASA NTRS1983

A new formulation of hybrid/mixed finite element

A new formulation of finite element method is accomplished by the Hellinger-Reissner principle for which the stress equilibrium conditions are not introduced initially but are brought-in through the use of additional internal displacement parameters. The method can lead to the same result as the assumed stress hybrid model. However, it is more general and more flexible. The use of natural coordinates for stress assumptions leads to elements which are less sensitive to the choice of reference coordinates. Numerical solutions by 3-D solid element indicate that more efficient elements can be constructed by assumed stresses which only partially satisfy the equilibrium conditions.

Pian, T. H. H.↗
DOE OSTI2023

Future Response of Coastal Wetlands to Environmental Stresses: Algorithm Comparison of Numerical Models

Coastal wetlands are a critical component of the earth system that strongly influence the global water and biogeochemical cycles. They are also likely important sentinel of climate change. Because eco-geomorphological processes have long been recognized to be important for coastal wetland survival under accelerated sea-level rise (SLR), many eco-geomorphology models have been developed to assess the impact of climate change on coastal wetlands. Although these models differ substantially in complexity and numerical methods, few studies have investigated the algorithm-level uncertainties in these models. In this study, we developed a multiple-algorithm model framework of coastal wetlands that represents coastal hydrodynamics (such as water level, significant wave height and bottom shear stress) and four eco-geomorphological processes: mineral accretion, organic matter accretion, storm surge erosion and landward migration. We validated the model at three representative coastal wetland sites (Venice Lagoon, Plum Island Estuary and Hunter Estuary) for hydrodynamics, mineral accretion and organic matter accretion. Through model-data comparison, we showed that the model can well capture the dynamics of hydrodynamical and eco-geomorphological conditions in the study sites. Importantly, analysis of the multiple-algorithm simulations suggests that differences in the process representation of mineral and organic matter accretion may contribute to the recent contradicting predictions of coastal wetland evolution under accelerated SLR.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Scalar flux transport models for self-similar turbulent mixing

A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.

Reynolds-averaged Navier Stokes↗
DOE OSTI2020

Stress-Drop and Source Scaling of the 2019 Ridgecrest, California, Earthquake Sequence

Stress drop, while difficult to measure reliably and at scale, is a key source parameter for understanding the earthquake rupture process and its relationship to strong ground motion. Here, we use a P-wave spectral decomposition approach, designed for large and densely sampled datasets, to measure earthquake stress drop in the region surrounding the 2019 Ridgecrest, California, earthquake sequence. With more than 11,000 measurements of earthquake stress drop in the 20-yr time period from 2000 through 2019, this dataset provides an opportunity to understand how coseismic stress changes and how other geophysical factors relate to the distribution of stress drop and its evolution in space and time. We observe a mild but persistent deviation from self-similar scaling, with larger events having systematically higher stress drops, though this trend depends on the assumption of an omega-square source spectral model. Earthquake stress drop increases with hypocentral depth in this study region, and the Ridgecrest aftershocks tend to have higher stress drops than the pre-event seismicity. This is in part due to their deeper hypocenters. Coherent spatial patterns of stress drop in the aftershock sequence correlate with the slip distribution of the M 7.1 mainshock, whose northwest rupture tip terminated in a long-lived zone of enervated stress drop. Although physical interpretation of these results is complicated by the trade-offs between the timing, depth, and location of these earthquakes, the observations provide new insight into the physics of the earthquake source in an area of renewed seismic activity in southern California.

58 GEOSCIENCES↗
NASA NTRS1990

Design of ceramic components with the NASA/CARES computer program

The ceramics analysis and reliability evaluation of structures (CARES) computer program is described. The primary function of the code is to calculate the fast-fracture reliability or failure probability of macro-scopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. CARES uses results from MSC/NASTRAN or ANSYS finite-element analysis programs to evaluate how inherent surface and/or volume type flaws component reliability. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effects of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or uniform uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for a single or multiple failure modes by using a least-squares analysis or a maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-to-fit-tests, 90 percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan 90 percent confidence band values are also provided. Examples are provided to illustrate the various features of CARES.

Nemeth, Noel N.↗
NASA NTRS1990

Ceramics Analysis and Reliability Evaluation of Structures (CARES). Users and programmers manual

This manual describes how to use the Ceramics Analysis and Reliability Evaluation of Structures (CARES) computer program. The primary function of the code is to calculate the fast fracture reliability or failure probability of macroscopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. The program uses results from MSC/NASTRAN or ANSYS finite element analysis programs to evaluate component reliability due to inherent surface and/or volume type flaws. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effect of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or unifrom uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for single or multiple failure modes by using the least-square analysis or the maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-of-fit tests, ninety percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan ninety percent confidence band values are also provided. The probabilistic fast-fracture theories used in CARES, along with the input and output for CARES, are described. Example problems to demonstrate various feature of the program are also included. This manual describes the MSC/NASTRAN version of the CARES program.

Nemeth, Noel N.↗
NASA NTRS1993

Turbulence modeling for non-equilibrium flows

Two projects are reported. The first is the development and testing of an eddy viscosity transport model. This project also is a starting point for our work on developing computational tools for solving turbulence models in complex geometries. The second project is a stochastic analysis of the realizability of Reynolds stress transport models.

Durbin, Paul A.↗
NASA NTRS1973

Numerical simulation of the world ocean circulation

A multi-level model, based on the primitive equations, is developed for simulating the temperature and velocity fields produced in the world ocean by differential heating and surface wind stress. The model ocean has constant depth, free slip at the lower boundary, and neglects momentum advection; so that there is no energy exchange between the barotropic and baroclinic components of the motion, although the former influences the latter through temperature advection. The ocean model was designed to be coupled to the UCLA atmospheric general circulation model, for the study of the dynamics of climate and climate changes. But here, the model is tested by prescribing the observed seasonally varying surface wind stress and the incident solar radiation, the surface air temperature and humidity, cloudiness and the surface wind speed, which, together with the predicted ocean surface temperature, determine the surface flux of radiant energy, sensible heat and latent heat.

Takano, K.↗
DOE OSTI2021

Assessment of a tomography-informed polyhedral discrete element modelling approach for complex-shaped granular woody biomass in stress consolidation

The design of handling equipment for granular biomass primarily requires experiments. Discrete element models (DEM) can provide designers with detailed insight into the behaviour of granular materials. However, granular biomass comprising complex-shaped particles is difficult to model with DEM. A tomography-informed DEM approach and an exhaustive assessment of the ability of the approach to predict the bulk behaviour of milled pines using experimental data is presented. Nano-CT scan was conducted to obtain 3D particle surface geometries as the basis for particle shape approximation by a polyhedral model and a sphero-polyhedral model. These models were applied in the simulation of a compressibility test. Our parametric study showed that particle Young's modulus and restitution coefficient are the two main properties influencing the simulated bulk behaviour of the DEM particles and that the level of approximating DEM particle shape for real particles is critical for accurately replicating the bulk behaviour of milled pines. The polyhedral model demonstrated better suitability than its sphero-polyhedral counterpart for modelling pine particles. The polyhedral model calibrated in a compressibility test was then applied in the simulations of a friction test without additional parameter tuning. Interlocking was dominant in the shear of bulk pine particles. Remarkably, the polyhedral model predicted the frictional behaviour of the pine particles when compared to the experimental results. The limitations of the model, as well as possible ways for enhancement, are discussed. This work provided novel insights into the suitability of complex-shaped DEM models for granular woody biomass.

Biomass flow↗
DOE OSTI2021

Characterizing Stress Orientations in Southern Kansas

We repot that Induced seismicity predominantly occurs along faults that are optimally oriented to the local principal compressive stress direction, and the characterization of these stress orientations is an important component of understanding seismic hazards. The seismicity rate in southern Kansas rapidly increased in 2013 primarily due to the disposal of large volumes of wastewater into the Arbuckle Group. Previously, local stress orientations in this area were poorly constrained, which limited our understanding of the complex faulting and diverse earthquake mechanisms in this region. We use shear-wave splitting and focal mechanism inversion techniques to create multiple, independent estimates of maximum horizontal stress directions (S Hmax ) and their spatial variation across the study area. We then create an integrated model of stress orientations for southern Kansas and northern Oklahoma using our local results in conjunction with previous, regional stress orientation estimates. We find that S Hmax in both southern Kansas and central Oklahoma exhibits an east-northeast (~N78°E) orientation, and these regions bound a northeast (~N59°E) rotation within a ~20km area in northern Oklahoma near the Nemaha ridge.

58 GEOSCIENCES↗
NASA NTRS2019

Structural Analysis and Performance-Based Validation of a Composite Wing Spar

Electric-motor powered aircraft possess the ability to operate with efficient energy delivery, but lack the operational range of internal combustion engine powered aircraft. This range limitation requires the use of high aspect ratio, thin-chord wings to minimize aerodynamic drag losses, which results in highly loaded composite spar structures. High aspect ratio wings are required to increase mission durations for a NASA-developed experimental multi-rotor electric powered aircraft denoted as the Scalable Convergent Electric Propulsion Technology and Operations Research (SCEPTOR) or “X-57”. This paper examines the structural performance of the composite main wing spars to validate spar strength using ply-based laminate finite element methods. Geometric scaling of a main spar test-section was initially proposed for proof-testing but sacrificed stability. Ply-based structures modeling with local structural features was implemented as a risk-reduction methodology. Ply-based modeling was selected to augment the conventional “building block” approach to reduce risk, and leverage a “performance-based” approval processes encouraged in Federal Aviation Administration (FAA) design guidance. Therefore, ply-based laminate modeling of the full-scale main spar and forward spar shear-web attachments were subsequently undertaken to determine load path complexity with predicted flight loads. Ply-based modeling included stress concentrations and interlaminate behavior at interface locations that can be obscured in traditional finite element sizing models. Analysis of the wing spar laminate ply-based models compared with bearing test coupon performance was used to reduce future wing assembly proof-testing burden and facilitate performance-based flight hardware safety for the X-57 experimental aircraft.

Moore, James B.↗