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 1,315 records · Page 73

Effects of moisture, residual thermal curing stresses, and mechanical load on the damage development in quasi-isotropic laminates

This investigation demonstrates how moisture absorbed in (0/+ or - 45/90)s and (0/90/+ or - 45)s graphite epoxy laminates significantly alters the stress state and chronology of damage development along the laminate edge during static tension and tension-tension cyclic loading. Emphasis is placed on using reasonable approximations for wet and dry elastic properties, including out-of-plane properties (nu sub 23 and G sub 23), since these properties are required by finite element and shear lag models to predict the stress state at the laminate edge. Moisture was observed to alter the dry edge stress state in the 90-deg plies of the (0/+ or - 45/90)s laminate such that delaminations occurred at a lower load and transverse cracks occurred at a higher load. A model was developed which predicted the differences in loads required to initiate damage in the 90-deg plies of the two laminates in the wet and dry conditions. Although moisture can alter the chronology of damage development, the damage state in each laminate observed prior to fracture appeared to be independent of moisture content.

Kriz, R. D.↗

Summary of Methodology for Mitigating Risks Associated with Licensing and Qualifying AM Nuclear Materials

The US Department of Energy’s Advanced Materials and Manufacturing Technologies (AMMT) program focuses on accelerating the development, qualification, demonstration, and deployment of advanced materials and manufacturing technologies to enable reliable and economical nuclear energy. Laser powder bed fusion (LPBF) is one of the most popular additive manufacturing (AM) processes for fabricating components with intrinsically complex geometries. LPBF was extensively explored for nuclear applications under the previous Transformational Challenge Reactor program. Additionally, Oak Ridge National Laboratory developed and licensed the Peregrine software and larger digital platform that couples machine learning and in situ data collection during AM to detect anomalies and any evolved defects. The digital platform will be critical to (1) the qualification of AM components for nuclear applications that link location-specific data to macroscopic properties and (2) predict final component performance. Current in situ process monitoring tools are valuable for observing the formation of stochastic flaws, but additional data are needed to predict the resulting microstructures and associated material performance. Rapid cooling rates and large thermal gradients have caused large heterogeneities in the microstructure, which cause anisotropy in mechanical performance. The AMMT program is evaluating the best approaches for addressing these heterogeneities and their effect on component performance using a combination of multiscale modeling, enhanced in situ process monitoring, and high throughput experimental testing. This report summarizes strategies for mitigating the risks associated with qualifying AM components, including developing new sensing capabilities for in situ process monitoring and characterizing melt pool solidification and residual stresses to inform multiscale modeling efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Wave effects on ocean-ice interaction in the marginal ice zone

The effects of wave train on ice-ocean interaction in the marginal ice zone are studied through numerical modeling. A coupled two-dimensional ice-ocean model has been developed to include wave effects and wind stress for the predictions of ice edge dynamics. The sea ice model is coupled to the reduced-gravity ocean model through interfacial stresses. The main dynamic balance in the ice momentum is between water-ice stress, wind stress, and wave radiation stresses. By considering the exchange of momentum between waves and ice pack through radiation stress for decaying waves, a parametric study of the effects of wave stress and wind stress on ice edge dynamics has been performed. The numerical results show significant effects from wave action. The ice edge is sharper, and ice edge meanders form in the marginal ice zone owing to forcing by wave action and refraction of swell system after a couple of days. Upwelling at the ice edge and eddy formation can be enhanced by the nonlinear effects of wave action; wave action sharpens the ice edge and can produce ice meandering, which enhances local Ekman pumping and pycnocline anomalies. The resulting ice concentration, pycnocline changes, and flow velocity field are shown to be consistent with previous observations.

Liu, Antony K.↗

Temperature-time distribution and thermal stresses on the RTG fins and shell during water cooling

Radioisotope thermoelectric generator (RTG) packages designed for space missions generally do not require active cooling. However, the heat they generate cannot remain inside of the launch vehicle bay and requires active removal. Therefore, before the Shuttle bay door is closed, the RTG coolant tubes attached to the heat rejection fins must be filled with water, which will circulate and remove most of the heat from the cargo bay. There is concern that charging a system at initial temperature around 200 C with water at 24 C can cause unacceptable thermal stresses in the RTG shell and fins. A computer model is developed to estimate the transient temperature distribution resulting from such charging. The thermal stresses resulting from the temperature gradients do not exceed the elastic deformation limit for the material. Since the simplified mathematical model for thermal stresses tends to overestimate stresses, it is concluded that the RTG can be cooled by introducing water at 24 C to the initially hot fin coolant tubes while the RTG is in the Shuttle cargo bay.

Turner, R. H.↗

Fracture behavior of unidirectional boron/aluminum composite laminates

An experiment was conducted to verify the results of mathematical models which predict the stresses and displacements of fibers and the amount of damage growth in a center-notched lamina as a function of the applied remote stress and the matrix and fiber material properties. A brittle lacquer coating was used to detect the yielding in the matrix while X-ray techniques were used to determine the number of broken fibers in the laminate. The notched strengths and the amounts of damage found in the specimens agree well with those predicted by the mathematical model. It is shown that the amount of damage and the crack opening displacement does not depend strongly on the number of plies in the laminate for a given notch width. By heat-treating certain laminates to increase the yield stress of the alumina matrix, the effect of different matrix properties on the fracture behavior was investigated. The stronger matrix is shown to weaken the notched laminate by decreasing the amount of matrix damage, thereby making the laminate more notch sensitive.

Goree, J. G.↗

A New Simplified Discrete Fracture Model for Shearing of Intersecting Fractures and Faults

Abstract Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model , to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.

58 GEOSCIENCES↗

Analyzing Mechanical integrity of Microwave resonant Cavity flowmeter

The microwave resonant cavity flow sensor has been proposed to address the difficulties of measuring high-temperature molten salt flow in reactor vessels. The flow sensor consists of a metallic cylindrical cavity coupled to a microwave waveguide. The cylindrical cavity has one of the flat surfaces thin enough to undergo microscopic deflection due to dynamic pressure of the exerted by the flowing fluid. Deflection of the membrane causes a shift in the microwave resonant frequency. Fluid velocity can be estimated from the calibration curve. The objective of this report is to investigate mechanical integrity of the flowmeter’s membrane. The membrane is a thin stainless-steel circular plate with a clamped edge, subjected to uniform load. By calculating the stress on the plate due to deflection and, comparing the stress to the material ultimate tensile strength and yield strength, it can be estimated if the plate will fail. If the maximum stress is greater than the yield strength, the plate will undergo plastic deformation, and the sensor will be disabled. A stress greater than the ultimate tensile strength will cause the plate to fracture. The stress on the plate was calculated with an analytic closed form solution model, and with COMSOL Structural Mechanics Module which does not involve any approximations. Stresses in both the analytic model and COMSOL numerical model were calculated for stainless steel 316 thin circular plates under room temperature conditions, where the fluid is water. The diameter of the circular plate is approximately one inch, and thickness is 10mil. Fluid flow velocity is in the range from 0.5m/s to 1m/s. The results from the two models are similar. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Evaluation of Inelastic Constitutive Models for Nonlinear Structural Analysis

The influence of inelastic material models on computed stress-strain states, and therefore predicted lives, was studied for thermomechanically loaded structures. Nonlinear structural analyses were performed on a fatigue specimen which was subjected to thermal cycling in fluidized beds and on a mechanically load cycled benchmark notch specimen. Four incremental plasticity creep models (isotropic, kinematic, combined isotropic-kinematic, combined plus transient creep) were exercised. Of the plasticity models, kinematic hardening gave results most consistent with experimental observations. Life predictions using the computed strain histories at the critical location with a Strainrange Partitioning approach considerably overpredicted the crack initiation life of the thermal fatigue specimen.

Kaufman, A.↗

Estimation of total above-ground phytomass production using remotely sensed data

Remote sensing potentially offers a quick and nondestructive method for monitoring plant canopy condition and development. In this study, multispectral reflectance and thermal emittance data were used in conjunction with micrometeorological data in a simple model to estimate above-ground total dry phytomass production of several spring wheat canopies. The fraction of absorbed photosynthetic radiation (PAR) by plants was estimated from measurements of visible and near-infrared canopy reflectance. Canopy radiation temperature was used as a crop stress indicator in the model. Estimated above-ground phytomass values based on this model were strongly correlated with the measured phytomass values for a wide range of climate and plant-canopy conditions.

Asrar, G.↗

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗

First Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Sensor Performance Model of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. Feasibility of the sensor was initially investigated with analytical derivations and with COMSOL RF Module computer simulations of resonant frequency spectral shift due to uniform load. We also investigated the mechanical integrity of the flowmeter’s membrane through analytical modelling and COMSOL Structural Mechanics Module computer simulations. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the material’s yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure. Using results from models, we have developed an initial design for a microwave K-band sensor. A cylindrical resonator prototype was fabricated from brass for the initial tests. The external dimensions of the cavity are matched to the flange of a standard WR-42 waveguide. Microwave field is coupled into the resonant cavity through a subwavelength-size aperture. A test article was developed consisting of a piping Tee with a bulkhead WR-42 microwave waveguide installed in a leak-proof assembly. A microwave waveguide circulator was installed in the setup to suppress the effect of reflections at the cavity entrance by increasing the isolation between the input and the output port. Preliminary spectral characterization of cavity spectral response was performed with a portable PXIe chassis microwave VNA with a custom GUI. Preliminary dry tests of the transducer response were conducted with a set of calibrated weights. Transducer frequency shift was shown to be monotonically increasing with increasing pressure. The next steps will involve investigation of the transducer performance for water flow sensing.

42 ENGINEERING↗

Evaluation of inelastic constitutive models for nonlinear structural analysis

The influence of inelastic material models on computed stress-strain states, and therefore predicted lives, was studied for thermomechanically loaded structures. Nonlinear structural analyses were performed on a fatigue specimen which had been subjected to thermal cycling in fluidized beds and on a mechanically load cycled benchmark notch specimen. Four incremental plasticity creep models (isotropic, kinematic, combined isotropic kinematic, combined plus transient creep) were exercised using the MARC program. Of the plasticity models, kinematic hardening gave results most consistent with experimental observations. Life predictions using the computed strain histories at the critical location with a strainrange partitioning approach considerably overpredicted the crack initiation life of the thermal fatigue specimen.

Kaufman, A.↗

Reynolds-Averaged Navier-Stokes Simulation of a 2D Circulation Control Wind Tunnel Experiment

Numerical simulations are performed using a Reynolds-averaged Navier-Stokes (RANS) flow solver for a circulation control airfoil. 2D and 3D simulation results are compared to a circulation control wind tunnel test conducted at the NASA Langley Basic Aerodynamics Research Tunnel (BART). The RANS simulations are compared to a low blowing case with a jet momentum coefficient, C(sub u), of 0:047 and a higher blowing case of 0.115. Three dimensional simulations of the model and tunnel walls show wall effects on the lift and airfoil surface pressures. These wall effects include a 4% decrease of the midspan sectional lift for the C(sub u) 0.115 blowing condition. Simulations comparing the performance of the Spalart Allmaras (SA) and Shear Stress Transport (SST) turbulence models are also made, showing the SST model compares best to the experimental data. A Rotational/Curvature Correction (RCC) to the turbulence model is also evaluated demonstrating an improvement in the CFD predictions.

Allan, Brian G.↗

Residual stresses in welded plates

The purpose of this project was to develop a simple model which could be used to study residual stress. The mechanism that results in residual stresses in the welding process starts with the deposition of molten weld metal which heats the immediately adjacent material. After solidification of weld material, normal thermal shrinkage is resisted by the adjacent, cooler material. When the thermal strain exceeds the elastic strain corresponding to the yield point stress, the stress level is limited by this value, which decreases with increasing temperature. Cooling then causes elastic unloading which is restrained by the adjoining material. Permanent plastic strain occurs, and tension is caused in the region immediately adjacent to the weld material. Compression arises in the metal farther from the weld in order to maintain overall static equilibrium. Subsequent repair welds may add to the level of residual stresses. The level of residual stress is related to the onset of fracture during welding. Thus, it is of great importance to be able to predict the level of residual stresses remaining after a weld procedure, and to determine the factors, such as weld speed, temperature, direction, and number of passes, which may affect the magnitude of remaining residual stress. It was hoped to use traditional analytical modeling techniques so that it would be easier to comprehend the effect of these variables on the resulting stress. This approach was chosen in place of finite element methods so as to facilitate the understanding of the physical processes. The accuracy of the results was checked with some existing experimental studies giving residual stress levels found from x-ray diffraction measurements.

Bernstein, Edward L.↗

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗

Negative Stress Margins - Are They Real?

Advances in modeling and simulation, new finite element software, modeling engines and powerful computers are providing opportunities to interrogate designs in a very different manner and in a more detailed approach than ever before. Margins of safety are also often evaluated using local stresses for various design concepts and design parameters quickly once analysis models are defined and developed. This paper suggests that not all the negative margins of safety evaluated are real. The structural areas where negative margins are frequently encountered are often near stress concentrations, point loads and load discontinuities, near locations of stress singularities, in areas having large gradients but with insufficient mesh density, in areas with modeling issues and modeling errors, and in areas with connections and interfaces, in two-dimensional (2D) and three-dimensional (3D) transitions, bolts and bolt modeling, and boundary conditions. Now, more than ever, structural analysts need to examine and interrogate their analysis results and perform basic sanity checks to determine if these negative margins are real.

Raju, Ivatury S.↗

Negative stress margins – are they real?

Advances in modeling and simulation, new finite element software, modeling engines and powerful computers are providing opportunities to interrogate designs in a very different manner and in a more detailed approach than ever before. Margins of safety are also often evaluated using local stresses for various design concepts and design parameters quickly once analysis models are defined and developed. This paper suggests that not all the negative margins of safety evaluated are real. The structural areas where negative margins are frequently encountered are often near stress concentrations, point loads and load discontinuities, near locations of stress singularities, in areas having large gradients but with insufficient mesh density, in areas with modeling issues and modeling errors, and in areas with connections and interfaces, in two-dimensional (2D) and three-dimensional (3D) transitions, bolts and bolt modeling, and boundary conditions. Now, more than ever, structural analysts need to examine and interrogate their analysis results and perform basic sanity checks to determine if these negative margins are real.

Raju, Ivatury S.↗

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗