Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain gradient 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

An Experimental and Analytical Evaluation of a Biaxial Test for Determining Shear Properties of Composite Materials

The results of an experimental and analytical investigation of a biaxial tension/compression test for determining shear properties of composite materials are reported. Using finite element models of isotropic and orthotropic laminates, a specimen geometry was optimized. A kinematic fixture was designed to introduce an equal and opposite pair of forces into a specimen with a one inch square test section. Aluminum and several composite laminates with the optimized geometry and a configuration with large stress gradients were tested in the fixture. The specimens were instrumented with strain gages in the center of the test section for shear stiffness measurements. Pure shear strain was measured. The results from the experiments correlated well with finite element results. Failure of the specimens occurred through the center of the test section and appeared to have initiated at the high stress points. The results lead to the conclusion that the optimized specimen is suitable for determining shear modulus for composite materials. Further revisions to the specimen geometry are necessary if the method is to give shear strength data.

Kennedy, John M.↗
NASA NTRS1990

Simulation of rotation effect in turbulence dissipation rate equation

It was observed that the Reynolds stress transport equations using the usual isotropic dissipation rate transport equations do not accurately predict the correct evolution of the kinetic energy in both homogeneous shear and isotropic decay. This same effect carries over to more practical flows such as those with strong irrotational strains and/or adverse pressure gradients. The deficiencies in these predictions were attributed to the need for modifications in the isotropic dissipation rate equation. Pope and Hanjalic and Launder have proposed modifications to this equation. These modifications were tested in the cases of homogeneous shear and iso-decay and their performance was analyzed. The results are compared to the rotating flow models of Bardina and Raj.

Raj, Rishi S.↗
NASA NTRS2004

Thermoelastic Formulation of Stiffened, Unsymmetric Composite Panels for Finite Element Analysis of High Speed Aircraft

An emerging technology need for capturing 3-D panel thermoelastic response with 2-D planar finite element models (FEMs) is aided with an equivalent plate stiffness and thermal coefficient formulation. The formulation is general and applies to all panel concepts. Included with the formulation is the ability to provide membrane-bending coupling of unsymmetric sections and calculation of all thermal expansion and bending responses from in-plane and through-the-thickness temperature gradients. Thermal residual strains for both the laminates and plies are included. The general formulation is defined and then applied to a hat-shaped, corrugated stiffened panel. Additional formulations are presented where required to include all of the hat's unique characteristics. Each formulation is validated independently with 3-D FEA.

Collier, Craig S.↗
DOE OSTI2024

An Engineering Approach for Weld Creep Lifetime Assessment Based on Local Property Measurement

The localized creep failure in the heat-affected zone (HAZ) of Grade 91 steel weldments has been identified as one of the most important factors causing significantly shortened service lifetime and structural integrity issues of welded components in advanced fossil and nuclear power plants. To conduct a reliable creep lifetime assessment, a new engineering assessment approach has been developed by incorporating the experimentally determined local properties of the heterogeneous HAZ. By creep testing a purposely simulated HAZ specimen with in situ digital image correlation (DIC) technique, the highly gradient creep properties across the HAZ of Grade 91 steel was quantitatively measured. A physical creep cavitation constitutive model was proposed to investigate the local creep deformation and damage accumulation within the heterogeneous HAZ, which takes into account the nucleation of creep cavities and their growth by both grain boundary diffusion and creep deformation. The relationship among the local material property, creep strain accumulation, and evolution characteristic of creep cavities was established. The approach was then utilized to investigate the creep response and subsequent life for an ex-service 9% Cr steel weldment by incorporating the effects of pre-existing damages which developed and accumulated during long-term services. The predicted results exhibited quantitative agreement with the DIC measurement in terms of both nominal/local creep deformation as well as the subsequent life under the test conditions at 650 and 80 MPa.

Zhang, Wei↗
NASA NTRS1988

MHOST: An efficient finite element program for inelastic analysis of solids and structures

An efficient finite element program for 3-D inelastic analysis of gas turbine hot section components was constructed and validated. A novel mixed iterative solution strategy is derived from the augmented Hu-Washizu variational principle in order to nodally interpolate coordinates, displacements, deformation, strains, stresses and material properties. A series of increasingly sophisticated material models incorporated in MHOST include elasticity, secant plasticity, infinitesimal and finite deformation plasticity, creep and unified viscoplastic constitutive model proposed by Walker. A library of high performance elements is built into this computer program utilizing the concepts of selective reduced integrations and independent strain interpolations. A family of efficient solution algorithms is implemented in MHOST for linear and nonlinear equation solution including the classical Newton-Raphson, modified, quasi and secant Newton methods with optional line search and the conjugate gradient method.

Nakazawa, S.↗
NASA NTRS1983

HOST structural analysis program overview

Hot section components of aircraft gas turbine engines are subjected to severe thermal structural loading conditions, especially during the start up and take off portions of the engine cycle. The most severe and damaging stresses and strains are those induced by the steep thermal gradients induced during the start up transient. These transient stresses and strains are also the most difficult to predict, in part because of the temperature gradients and distributions are not well known or readily predictable, and also because the cyclic elastic viscoplastic behavior of the materials at these extremes of temperature and strain are not well known or readily predictable. A broad spectrum of structures related technology programs is underway to address these deficiencies. One element of the structures program is developing improved time varying thermal mechanical load models for the entire engine mission cycle from start up to shutdown. Another major part of the program is the development of new and improved nonlinear 3-D finite elements and associated structural analysis programs, including the development of temporal elements with time dependent properties to account for creep effects in the materials and components.

Johns, R. H.↗
DOE OSTI2022

Report Summarizing Progress in Digital Image Correlation Analysis of Burst Phenomenon

This report is an update on the progress made towards utilizing in-situ methods to observe and quantify Zry-4 deformation behavior during burst testing. A heavily modified burst test was developed that allows both digital image correlation and infrared thermography to be utilized during testing to quantify strain and axial thermal gradient data during burst. A variety of tests ranging from 6.2-15.9 MPa internal pressure and utilizing both constant and transient pressure conditions were conducted and presently reported on. The test modifications, which are detailed extensively in a previous report, are effectively found to have no discernible impact on burst temperatures when compared to historical Zircaloy burst data. To support future modeling studies, the IR thermography data was used to approximate the thermal gradient down the axial length at burst for all the claddings. Then, DIC is first correlated to thermocouple temperature for a 10.3 MPa Transient pressure test, and then DIC correlated to IR thermal gradient data down the axial length of a 8.2 MPa Transient pressure test at several points.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unraveling kinking: A plasticity enhancing failure mode in high strength nano metallic laminates

Kinking is an important and plasticity-enhancing deformation/failure mode in numerous mechanically anisotropic materials including high-strength nano metallic laminates (NMLs). However, our current limited understanding of the mechanics of kinking and its dependence on microstructural attributes is insufficient for thoroughly comprehending and eventually being able to control failure behaviors of materials. In this study, we investigate kinking dependencies on microstructural attributes in NMLs via in situ micropillar compression, multiscale microstructure characterization, dislocation dynamic simulations, and crystal plasticity modeling. Additionally, by examining several NML systems (Cu/Fe, Ag/Fe, Al-4Mg/Fe), we demonstrate that the development of internal stresses during loading activates local layer-parallel glide triggering kinking in NMLs. Furthermore, this work reveals the effect of key microstructural features including layer thickness, layer waviness, interface barrier strength, and work hardening capacity on kink band formation in NMLs. More broadly, our efforts represent a generically applicable approach for probing large-strain deformation behavior of complex materials via synergetic modeling and experimental efforts.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Thrust chamber life prediction

The reusable life of the Space Shuttle main engine (SSME) is influenced by the cyclic life of the regeneratively liquid cooled main combustion chamber (MCC). During an operational duty cycle the MCC liner is subjected to a large transient thermal gradient that imparts a high thermal cyclic strain to the liner hot gas wall. Life predictions of such chambers have usually been based on low cycle fatigue (LCF) evaluations. Hot-fire testing, however, has shown significant mid-channel wall deformation and thinning during accrued cyclic testing. This phenomenon is termed cyclic creep and appears to be significantly accelerated at elevated temperatures. An analytical method that models the cyclic creep phenomenon and its application to thrust chamber life prediction is presented. The chamber finite element geometry is updated periodically to account for accrued wall thinning and distortion. Failure is based on the tensile instability failure criterion. Cyclic life results for several chamber life enhancing coolant channel designs are compared to the typically used LCF analysis that neglects cyclic creep. The results show that the usable cyclic creep life is approximately 30 to 50% of the commonly used LCF life.

Kasper, H. J.↗
NASA NTRS1997

U-Groove Aluminum Weld Strength Improvement

Though butt-welds are among the most preferred joining methods in aerostructures, their strength dependence on inelastic mechanics is generally the least understood. This study investigated experimental strain distributions across a thick aluminum U-grooved weld and identified two weld process considerations for improving the multipass weld strength. One is the source of peaking in which the extreme thermal expansion and contraction gradient of the fusion heat input across the groove tab thickness produces severe angular distortion that induces bending under uniaxial loading. The other is the filler strain hardening decreasing with increasing filler pass sequences, producing the weakest welds on the last weld pass side. Both phenomena are governed by weld pass sequences. Many industrial welding schedules unknowingly compound these effects, which reduce the weld strength. A depeaking index model was developed to select filler pass thickness, pass numbers, and sequences to improve depeaking in the welding process. The result was to select the number and sequence of weld passes to reverse the peaking angle such as to combine the strongest weld pass side with the peaking induced bending tension component side to provide a more uniform stress and stronger weld under axial tensile loading.

Verderaime, V.↗
DOE OSTI2023

Reaction–Diffusion Coupling Facilitates the Sequential Precipitation of Metal Ions from Battery Feedstock Solutions

Here, the development of new technologies for chemical separations is urgently needed to meet the surging demand for critical materials that has strained resources and caused environmental challenges. Inspired by the classic Liesegang experiment, we demonstrated the separation of critical metal ions based on the coupling of ion diffusion and precipitation kinetics. For this purpose, a model feedstock solution simulating dissolved battery electrodes was placed on top of a hydrogel loaded with a precipitating agent, namely sodium hydroxide. As the lithium, manganese, cobalt, and nickel ions diffused into the gel, a gradient of precipitates formed along the length of the reactor. Elemental analysis of the spatially distributed precipitates showed the enrichment of nickel near the gel-solution interface, followed by the formation of an almost pure (>96%) manganese product further along the reactor. Optimization experiments revealed that a sodium hydroxide concentration of 10 mM and a gel/solution volume ratio of 2:1 favored efficient separations. The robustness of the method was demonstrated in four out of five feedstock compositions of typically used battery cathodes. Our proof-of-concept experiments present a paradigm for critical materials separations that does not require specialty chemicals, binding agents, membranes, or toxic solvents.

25 ENERGY STORAGE↗
DOE OSTI2025

Frontal Maintenance in Submesoscale Flows

Classic deformation theory includes parameters—divergence, total strain, and vorticity—that are invariant to changes in the coordinate system. However, these parameters are sometimes ambiguous with respect to characterizing how fronts are formed and maintained because the presence of a front imposes a reference coordinate system. To help remedy this ambiguity, we propose a framework in frontal coordinates based on along- and cross-front velocity gradients to better characterize frontal maintenance, which can also be used to define divergence and normal strain in frontal coordinates. The framework with these four parameters (along-, cross-front velocity gradients, divergence, and normal strain in frontal coordinate) defines eight characteristic flow types at a front, providing a complete characterization of the flow that strengthens or weakens a front. This framework highlights the importance of the “strain efficiency” concept, which unambiguously defines the contribution of total strain to frontogenesis. Two examples, one based on a realistic simulation of submesoscales in the northern Gulf of Mexico and the other based on an idealized model with similar flow characteristics, are provided to demonstrate how this framework can be used to enhance our understanding of frontal dynamics in submesoscale flows.

frontogenesis↗
NASA NTRS2008

Performance Evaluation and Modeling of Erosion Resistant Turbine Engine Thermal Barrier Coatings

The erosion resistant turbine thermal barrier coating system is critical to the rotorcraft engine performance and durability. The objective of this work was to determine erosion resistance of advanced thermal barrier coating systems under simulated engine erosion and thermal gradient environments, thus validating a new thermal barrier coating turbine blade technology for future rotorcraft applications. A high velocity burner rig based erosion test approach was established and a new series of rare earth oxide- and TiO2/Ta2O5- alloyed, ZrO2-based low conductivity thermal barrier coatings were designed and processed. The low conductivity thermal barrier coating systems demonstrated significant improvements in the erosion resistance. A comprehensive model based on accumulated strain damage low cycle fatigue is formulated for blade erosion life prediction. The work is currently aiming at the simulated engine erosion testing of advanced thermal barrier coated turbine blades to establish and validate the coating life prediction models.

Miller, Robert A.↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
NASA NTRS1985

Dynamic Creep Buckling: Analysis of Shell Structures Subjected to Time-dependent Mechanical and Thermal Loading

The objective of the present research is to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermomechanical loads. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) can be anticipated and must be considered in developing the mathematical model. A complete, true ab-initio rate theory of kinematics and kinetics for continuum and curved thin structures, without any restriction on the magnitude of the strains or the deformations, was formulated. The time dependence and large strain behavior are incorporated through the introduction of the time rates of metric and curvature in two coordinate systems: fixed (spatial) and convected (material). The relations between the time derivative and the covariant derivative (gradient) were developed for curved space and motion, so the velocity components supply the connection between the equations of motion and the time rates of change of the metric and curvature tensors.

Simitses, G. J.↗
NASA NTRS1997

Enhancement of Mixing and Reaction in High-speed Combustor Flowfields

Research has been underway for a number of years, both in the United States and abroad, to develop advanced aerospace propulsion systems for use late in this century and beyond. One program is now underway at the NASA Langley Research Research has been underway for a number of years, both in the United States and abroad, to develop advanced aerospace propulsion systems for use late in this century and beyond. One program is now underway at the NASA Langley Research Center to develop a hydrogen-fueled supersonic combustion ramjet (scramjet) that is capable of propelling a vehicle at hypersonic speeds in the atmosphere. A part of that research has been directed toward the optimization of the scramjet combustor and, in particular, the efficiency of fuel-air mixing and reaction taking place in the engine. In the high-speed vehicle configurations currently being considered, achieving a high combustor efficiency becomes particularly difficult. With increasing combustor Mach number, the degree of fuel-air mixing that can be achieved through natural convective and diffusive processes is reduced leading to an overall decrease in combustion efficiency and thrust. Compressible shear layers and jets provide a good model for studying the physical processes occurring in high speed mixing and combustion in a scramjet. Mixing layers are characterized by large scale eddies that form due to the high shear that is present between the fuel and the air streams. These eddies entrain the fuel and air into the mixing region between the fluids leading to increased surface area and locally steep concentration gradients. Molecular diffusion then occurs across the strained interfaces. In an early study of high-speed mixing, Brown and Roshko show that the spreading rate of a supersonic mixing layer decreases with increasing Mach number, exhibiting a factor of three decrease in spread rate as compared with an incompressible mixing layer with the same density ratio. They conclude that the reduced spread rate is primarily due to compressibility. Papamoschou and Roshko also observe that the spreading rate of compressible mixing layers is significantly reduced over that of incompressible layers. To characterize the structure of the flow quantitatively, they define a convective Mach number. The reduction in mixing layer spreading rate is shown in these experiments to correlate well with increasing convective Mach number.

Drummond, J. Philip↗
NASA NTRS1999

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert, spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport, For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗
NASA NTRS1999

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport. For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗