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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Thermomechanical deformation testing and modeling in the presence of metallurgical instabilities

A number of viscoplastic constitutive models were developed to describe deformation behavior under complex combinations of thermal and mechanical loading. Questions remain, however, regarding the validity of procedures used to characterize these models for specific structural alloys. One area of concern is that the majority of experimental data available for this purpose are determined under isothermal conditions. This experimental study is aimed at determining whether viscoplastic constitutive theories characterized using an isothermal data base can adequately model material response under the complex thermomechanical loading conditions typical of power generation service. The approach adopted was to conduct a series of carefully controlled thermomechanical experiments on a nickel-based superalloy, Hastelloy Alloy X. Previous investigations had shown that this material experiences metallurgical instabilities leading to complex hardening behavior, termed dynamic strain aging. Investigating this phenomenon under full thermomechanical conditions leads to a number of challenging experimental difficulties which up to the present work were unresolved. To correct this situation, a number of advances were made in thermomechanical testing techniques. Advanced methods for dynamic temperature gradient control, phasing control and thermal strain compensation were developed and incorporated into real time test control software. These advances allowed the thermomechanical data to be analyzed with minimal experimental uncertainty. The thermomechanical results were evaluated on both a phenomenological and microstructural basis. Phenomenological results revealed that the thermomechanical hardening trends were not bounded by those displayed under isothermal conditions. For the case of Hastelloy Alloy X (and similar dynamic strain aging materials), this strongly suggests that some form of thermomechanical testing is necessary when characterizing a thermoviscoplastic deformation model. Transmission electron microscopy was used to study the microstructural physics, and analyze the unique phenomenological behavior.

Castelli, Michael G.↗
NASA NTRSDate not supplied

New Developments in the CompDam Fiber Kinking Model for the Interaction

This paper describes two new additions to the CompDam fiber kinking model: 1) a failure index for fiber kinking to facilitate interaction between kinking and matrix cracking and 2) a strain-based criterion for fiber rupture to account for fiber breakage before or after fiber kinking. The model formulation is described with the aid of the graphical Considère construction, which is used to determine the conditions for fiber kinking. Verification studies are conducted with single element models to show the behavior of the model in light of the new developments and to evaluate the effect of numerical parameters of the model. The model is used to analyze open hole compression (OHC) specimens with soft, quasi-isotropic, and hard layups. Analysis results are shown for a conventional continuum damage mechanics model using a trilinear constitutive law for fiber damage and for the proposed model to highlight the differences between the two. Comparison with test data highlights the capabilities and remaining gaps in the quantitative prediction of the damage and failure response. While both models can predict many features of the OHC response, accurate prediction of kinking onset in the OHC specimens remains elusive. The results suggest that it is necessary to predict the combined effects of compression and shear loading as well as the structural effects (e.g., ply thickness, strain gradient effect) accurately to match the experimental observations of fiber kinking onset in the three OHC laminates.

Andrew C. Bergan↗
NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗
NASA NTRS1992

Joints in Venusian lava flows

The bright polygonal patterns seen in some Magellan images of the volcanic plains are interpreted as joint networks. Several areas displaying such patterns are identified, noting their areal extent, the presence of volcanic features, and the regional topography. Intersection geometries of the lineations defining the polygonal patterns are typical of those found in tensile networks. In addition, the polygonal patterns exhibit no preferred orientation, implying that they are the result of horizontally isotropic stress fields. It is demonstrated that sufficient tensile stresses develop in a cooling Venusian flow to allow joints to form and that the thermal strain predicted is not incompatible with that observed. It is suggested that the tensioned networks are the result of a thermally stressed layer much thicker than a single lava flow. A second model is also considered in which thermal stresses develop as a result of regional elevation of the temperature gradient.

Johnson, Catherine L.↗
NASA NTRS2022

A Dynamic Nonlinear Subgrid-Scale Model for Large-Eddy Simulation of Complex Turbulent Flows

We present a new dynamic nonlinear subgrid-scale (SGS) model for large-eddy simulations (LES) and apply it to compute a flow involving pressure gradients, surface curvature and separation, for which data from a direct numerical simulation are available for comparison. The model, inspired by the triple model idea of Bardina et al. (“Improved Turbulence Models Based on Large Eddy Simulation of Homogeneous, Incompressible, Turbulent Flows,” Report No. TF-19, Thermosciences Division, Department of Mechanical Engineering, Stanford University, 1983), includes a Galilean-invariant term called the modified Leonard stress tensor, and two nonlinear terms comprised of the products of the strain-rate and rotation-rate tensors for an improved representation of the subgrid-scale dissipation, backscatter and anisotropy effects. The model does not employ any ad hoc averaging or clipping procedures, and does not require the specification of a characteristic length scale; hence, it naturally avoids the ambiguities associated with defining a proper length scale for anisotropic grids. Results from the wall-resolved LES of flow past a Gaussian bump using the new model demonstrate improved prediction of skin-friction, flow separation, mean flow profiles and turbulent quantities when compared to implicit LES as well as explicit LES using the Vreman SGS model on the same grid.

Turbulence↗
NASA NTRS1990

A study of turbulence on compression ramps with k-epsilon and Reynolds stress models

A theoretical study was conducted to determine the effects of adverse pressure gradient and compressibility in modeling turbulent compressible flows. The kinetic energy/dissipation and Reynolds stress model predictions are presented and compared with experimental data. The effects of compressibility, which include the mass averaged fluctuation term, the pressure dilatation term, and the dilatation dissipation, are important in modeling the turbulent compressible flows. The normal stresses and longitudinal strain rates also have an effect in the prediction of turbulent energy productions on the curved surfaces. A new compressible formulation of the pressure strain term, which includes the dilatation effects, in the Reynolds stress equation is presented.

Lee, J.↗
NASA NTRS1996

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. The extreme thermal expansion and contraction gradient of the fusion heat input across the groove tab thickness produces severe peaking, which induces bending under uniaxial loading. The filler strain-hardening decreased with increasing filler pass sequence, producing the weakest welds on the last pass side. Current welding schedules unknowingly compound these effects which reduce the weld strength. A depeaking index model was developed to select filler pass thicknesses, pass numbers, and sequences to improve depeaking in the welding process. The intent is to combine the strongest weld pass side with the peaking induced bending tension to provide a more uniform stress and stronger weld under axial tensile loading.

Verderaime, V.↗
NASA NTRS2016

A Continuum Damage Mechanics Model to Predict Kink-Band Propagation Using Deformation Gradient Tensor Decomposition

A new model is proposed that represents the kinematics of kink-band formation and propagation within the framework of a mesoscale continuum damage mechanics (CDM) model. The model uses the recently proposed deformation gradient decomposition approach to represent a kink band as a displacement jump via a cohesive interface that is embedded in an elastic bulk material. The model is capable of representing the combination of matrix failure in the frame of a misaligned fiber and instability due to shear nonlinearity. In contrast to conventional linear or bilinear strain softening laws used in most mesoscale CDM models for longitudinal compression, the constitutive response of the proposed model includes features predicted by detailed micromechanical models. These features include: 1) the rotational kinematics of the kink band, 2) an instability when the peak load is reached, and 3) a nonzero plateau stress under large strains.

Bergan, Andrew C.↗
NASA NTRS1988

On the origin of steep edges and filaments in vorticity and potential vorticity fields

High-resolution numerical calculations are shown which capture the fundamental process responsible for the intensification of vorticity gradients in an isolated vortex subject to externally imposed distrubances. Imposition of almost any weak strain or shear field in stripping away the relatively weak vorticity at the edge of the vortex and leaves it with gradients four to six orders of magnitude greater than in the initial state. Calculations displaying such enormous gradients have never been reported previously, because of the artificial eddy diffusivities that always limit such gradients in standard numerical models. The present calculations, which have no such limitations, have been made possible by the development of a novel and robust new numerical technique for vortex dynamics called contour surgery.

Dritschel, D. G.↗
NASA NTRS2024

Impact Response of Thermoplastic Composites – Experiments and Modeling

This study investigates the impact response of thermoplastic composite materials. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK reinforced with T700G (T700/LMPAEK ). Several High Energy Dynamic Impact (HEDI) tests were performed wherein an aluminum projectile impacted T700/LMPAEK panels at different velocities. In the impact experiment studied herein, the projectile impacted the T700/LMAPEK panel at 73.2 m/s and rebounded. The ensuing damage to the T700/LMPAEK panel was characterized via ultrasonic C-scans, visualizing damage modes including matrix cracking and delaminations. Since the projectile rebounded, the fiber breakage damage mode was not evident. The impact event was modeled using a combined continuum damage mechanics and cohesive zone analysis approach with the commercially available hydrocode LS-DYNA®. The dynamic deformation and damage of the continuum, i.e., the plies in the T700/LMPAEK laminate, was modeled using a Deformation Gradient Decomposition (DGD) material model, i.e., MAT299 available within LS-DYNA. The interlaminar regions were treated as zero-thickness cohesive zones, wherein the dynamic delamination mechanics was modeled using a traction-separation formulation. Mixed-mode delamination growth was captured using the Benzegaggh-Kenane (B-K) traction-separation interpolation law. Strain-rate sensitivity of interlaminar fracture toughness was included using approximations based on published literature. The outcomes of the modeling effort are discussed herein which demonstrate good correlations with test data regarding panel deflection, panel-projectile interactions, damage modes and damage extent.

Impact modeling↗
NASA NTRS1995

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. The extreme thermal expansion and contraction gradient of the fusion heat input across the groove tab thickness produces severe peaking which induces bending under uniaxial loading. The filler strain-hardening deceased with increasing filler pass sequence, producing the weakest welds on the last pass side. Current welding schedules unknowingly compound these effects which reduce the weld strength. A de-peaking index model was developed to select filler pass thicknesses, pass numbers, and sequences to improve de-peaking in the welding process. Intent is to combine the strongest weld pass side with the peaking induced bending tension to provide a more uniform stress and stronger weld under axial tensile loading.

Verderaime, V.↗
NASA NTRS1996

Turbulence Modeling and Computation of Turbine Aerodynamics and Heat Transfer

The objective of the present research is to develop improved turbulence models for the computation of complex flows through turbomachinery passages, including the effects of streamline curvature, heat transfer and secondary flows. Advanced turbulence models are crucial for accurate prediction of rocket engine flows, due to existance of very large extra strain rates, such as strong streamline curvature. Numerical simulation of the turbulent flows in strongly curved ducts, including two 180-deg ducts, one 90-deg duct and a strongly concave curved turbulent boundary layer have been carried out with Reynolds stress models (RSM) and algebraic Reynolds stress models (ARSM). An improved near-wall pressure-strain correlation has been developed for capturing the anisotropy of turbulence in the concave region. A comparative study of two modes of transition in gas turbine, the by-pass transition and the separation-induced transition, has been carried out with several representative low-Reynolds number (LRN) k-epsilon models. Effects of blade surface pressure gradient, freestream turbulence and Reynolds number on the blade boundary layer development, and particularly the inception of transition are examined in detail. The present study indicates that the turbine blade transition, in the presence of high freestream turbulence, is predicted well with LRN k-epsilon models employed. The three-dimensional Navier-Stokes procedure developed by the present authors has been used to compute the three-dimensional viscous flow through the turbine nozzle passage of a single stage turbine. A low Reynolds number k-epsilon model and a zonal k-epsilon/ARSM (algebraic Reynolds stress model) are utilized for turbulence closure. An assessment of the performance of the turbulence models has been carried out. The two models are found to provide similar predictions for the mean flow parameters, although slight improvement in the prediction of some secondary flow quantities has been obtained by the ARSM model. It's found that the wake profiles inside the endwall boundary layers are predicted better than those near the mid-span.

Lakshminarayana, B.↗
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.↗
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.↗
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.↗