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,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗
NASA NTRS1992

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗
NASA NTRS1992

Simplified micromechanical equations for thermal residual stress analysis of coated fiber composites

The fabrication of metal matrix composites poses unique problems to the materials engineer. The large thermal expansion coefficient mismatch between the fiber and matrix leads to high tensile residual stresses at the fiber/matrix interface which could lead to premature matrix cracking during cooldown. Fiber coating could be used to reduce thermal residual stresses. A simple closed-form analysis, based on a three-phase composite cylinder model, was developed to calculate thermal residual stresses in a fiber/interface/matrix system. Guidelines, in the form of simple equations, for the selection of appropriate material properties of the fiber coating, were also derived to minimize thermal residual stresses in the matrix during fabrication.

Naik, Rajiv A.↗
NASA NTRSDate not supplied

Application of the MAT 213 Composite Impact Model to NASA Problems of Interest

A material model has been developed which incorporates several key capabilities which have been identified as lacking in currently available composite impact models. The material model utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength). The material model has been implemented into LS-DYNA® as MAT 213. The model can simulate the nonlinear deformation, damage and failure that takes place in a composite under dynamic loading conditions. The deformation model utilizes an orthotropic plasticity formulation. For the damage model, the nonlinear unloading response that is observed prior to the point where the peak stress is reached can be simulated, as well as the stress degradation response that occurs after the peak stress is reached. A variety of failure models, including a generalized tabulated failure model which facilitates the utilization of general failure surfaces, have been implemented into MAT 213. Recent studies at NASA have concentrated on using MAT 213 to analyze both the impact and crush response of a variety of laminated and textile architectures. Several of these studies will be discussed in the presented paper. For example, a woven carbon/Kevlar composite is currently being examined for use in an energy absorber system for rotorcraft structures. MAT 213 analyses have been conducted to examine the ability of the model to accurately simulate the dynamic crush response of this material. Studies are also being conducted to examine the ability of MAT 213 to simulate the ballistic impact response of representative laminated and woven thermoplastic and thermoset matrix composites. NASA efforts are also concentrated on developing methods and processes for improving characterization methods and developing “best practices” for using MAT 213, a summary of which will be discussed in the paper.

Polymer Matrix Composites↗
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗
DOE OSTI2025

Lattice Structured Lightweight Structural Materials

The development of lightweight structural materials is crucial for enhancing the performance and deployment feasibility of fission batteries. This study aims to produce lightweight structural materials whose strength-to-weight ratios exceed those of current widely used structural materials. To achieve this, advanced modeling and simulation tools were employed to design lattice structures with different lattice parameters and different lattice types. A process was successfully developed for transforming lattice-structured models into Multiphysics Object Oriented Simulation Environment (MOOSE) inputs. Finite element modeling (FEM) was used to simulate the uniaxial tensile testing of the lattice-structured parts to investigate the stress distribution at a given displacement. The modeling results showed that the lattice-structured sample displayed a lower Young’s modulus in comparison to the solid material; the increase in solid shell thickness and blend radius enhances the mechanical performance; and the effect of unit cell size on macro scale stress is minimal. Tensile testing was conducted on the solid and lattice-structured materials fabricated by laser powder bed fusion (LPBF) additive manufacturing. The experimental results agreed well with the model prediction. The approach of using modeling as a guiding tool for preliminary material design can significantly save time and cost for new material development.

lightweight material↗
DOE OSTI2026

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗
NASA NTRS1993

Accuracy of the Generalized Self-Consistent Method in Modelling the Elastic Behaviour of Periodic Composites

Local stress and strain fields in the unit cell of an infinite, two-dimensional, periodic fibrous lattice have been determined by an integral equation approach. The effect of the fibres is assimilated to an infinite two-dimensional array of fictitious body forces in the matrix constituent phase of the unit cell. By subtracting a volume averaged strain polarization term from the integral equation we effectively embed a finite number of unit cells in a homogenized medium in which the overall stress and strain correspond to the volume averaged stress and strain of the constrained unit cell. This paper demonstrates that the zeroth term in the governing integral equation expansion, which embeds one unit cell in the homogenized medium, corresponds to the generalized self-consistent approximation. By comparing the zeroth term approximation with higher order approximations to the integral equation summation, both the accuracy of the generalized self-consistent composite model and the rate of convergence of the integral summation can be assessed. Two example composites are studied. For a tungsten/copper elastic fibrous composite the generalized self-consistent model is shown to provide accurate, effective, elastic moduli and local field representations. The local elastic transverse stress field within the representative volume element of the generalized self-consistent method is shown to be in error by much larger amounts for a composite with periodically distributed voids, but homogenization leads to a cancelling of errors, and the effective transverse Young's modulus of the voided composite is shown to be in error by only 23% at a void volume fraction of 75%.

Walker, Kevin P.↗
NASA NTRS1993

Fabrication and properties of functionally graded NiAl/Al2O3 composites

A modified sedimentation process was used in the production of a functionally gradient material (FGM), NiAl/Al2O3. A simple finite element model was used to guide our design and fabrication efforts by estimating residual stress states as a function of composite structure. This approach could lead to tailored designs that enhance or avoid specific residual stress states. Thermal cycling tests were factored into the model to predict time dependent or steady-state internal temperature and stress profiles. Four-point bend tests were conducted to establish the mechanical load-displacement behavior of a single interlayer FGM at room temperature, 800 and 1000 K. Room temperature bend strength of the FGM was 3-4 times that of the base NiAl. At elevated temperatures, composite fracture occurred in a gradual, noncatastrophic mode involving NiAl retardation of a succession of cracks originating in the alumina face.

Miller, D. P.↗
DOE OSTI2025

Quantifying precursors to void nucleation and coalescence in aluminum

Ductile rupture is a common failure mode for engineering alloys. It is generally comprised of three mechanisms: void nucleation at secondary particles, void growth, and coalescence of voids. Conventional models for ductile rupture have limited precision for predicting macroscopic failure. This is partially because they have been corroborated using ex-situ observations of these three mechanisms and calibrated by macroscopic strain or stress metrics. In this study, in-situ high-energy X-ray characterization was conducted to correlate sites of void nucleation via particle cracking and sites of void growth with grain-scale metrics. Here it is shown that particle cracking is not predicated by elevated stress metrics, which disagrees with classical models. Instead, particle cracking tended to occur for the largest, least spherical particles. This trend persisted on the aggregate-scale and within individual neighborhoods around particles. Furthermore, an Eshelby analysis illuminated a statistically significant increase in maximum principal stress within a cracked particle, compared to an uncracked particle. In-situ observations also revealed two separate examples of intragranular void growth comprised of flat, crack-like features that grew showing alignment with slip systems of the highest resolved shear stress. Post-mortem fractography revealed a secondary population of dimples on these crack-like features, implying that void sheeting may be occurring during this crack-like void growth. Furthermore, these in-situ observations imply that classical models for void growth, that assume a homogenous medium, should be extended to account for the anisotropic grain-level behavior to correctly capture distinct features of void nucleation and growth.

Al-2219↗
NASA NTRS1992

Finite Element Modeling of the Bulk Magnitization of Railroad Wheels to Improve Test Conditions for Magnetoacoustic Residual Stress Measurements

The magnetoacoustic measurement technique has been used successfully for residual stress measurements in laboratory samples. However, when used to field test samples with complex geometries, such as railroad wheels, the sensitivity of the method declines dramatically. It has been suggested that the decrease in performance may be due, in part, to an insufficient or nonuniform magnetic induction in the test sample. The purpose of this paper is to optimize the test conditions by using finite element modeling to predict the distribution of the induced bulk magnetization of railroad wheels. The results suggest that it is possible to obtain a sufficiently large and uniform bulk magnetization by altering the shape of the electromagnet used in the tests. Consequently, problems associated with bulk magnetization can be overcome, and should not prohibit the magnetoacoustic technique from being used to make residual stress measurements in railroad wheels. We begin by giving a brief overview of the magnetoacoustic technique as it applies to residual stress measurements of railroad wheels. We then define the finite element model used to predict the behavior of the current test configuration along with the nonlinear constitutive relations which we obtained experimentally through measurements on materials typically used to construct both railroad wheels and electromagnets. Finally, we show that by modifying the pole of the electromagnet it is possible to obtain a significantly more uniform bulk magnetization in the region of interest.

Fulton, J. P.↗
DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗
DOE OSTI2025

Bayesian calibration of irradiated graphite property models under high temperatures

Graphite under high temperatures and irradiation is central to advanced reactors. We develop a Bayesian calibration framework for graphite property models that explicitly represents model-data mismatch via a Gaussian-process discrepancy. The approach propagates uncertainty from parameters, experimental noise, and model form, with a hierarchical variance structure to capture group and cross-group noise. Using two predictive models across five grades (IG-110, NBG-18, PCEA, NBG-17, 2114) and four properties-irradiation-induced dimension change, creep, Young’s modulus change ratio, and coefficient of thermal expansion change ratio-we obtain average predictive-error reductions of 54%, 65%, 17%, and 17% when discrepancy is included. We illustrate engineering impact with a multiphysics model of a very-high-temperature reactor prismatic reflector brick, analyzing stresses under high fluence and temperature. Accounting for model discrepancy markedly improves predictive accuracy and provides a robust basis for reliable graphite component design in advanced reactors.

36 - MATERIALS SCIENCE↗
NASA NTRS1981

Ore deposits in Africa and their relation to the underlying mantle

African magmatism is largely related to the tensional stress regimes of the crust which are induced by the hotter upwelling mantle rocks. These mantle rocks may provide emanating forces and thermal energy for the upward movements of primary ore bodies with fluid inclusions in the tensional stress regimes of the crust. In this paper, the Goddard Earth Gravity Model is used to calculate a detailed subcrustal stress system exerted by mantle convection under Africa. The resulting system is found to be correlated with the African metallogenic provinces. Recognition of the full spectrum of ore deposits in Africa that may be associated with the hotter upwelling mantle rocks has provided an independent evidence to support the hypothesis of mantle-derived heat source for ore deposits.

Liu, H.-S.↗

Refine your search to explore more results.