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 181 records · Page 10

Population dynamics during cell proliferation and neuronogenesis in the developing murine neocortex

During the development of the neocortex, cell proliferation occurs in two specialized zones adjacent to the lateral ventricle. One of these zones, the ventricular zone, produces most of the neurons of the neocortex. The proliferating population that resides in the ventricular zone is a pseudostratified ventricular epithelium (PVE) that looks uniform in routine histological preparations, but is, in fact, an active and dynamically changing population. In the mouse, over the course of a 6-day period, the PVE produces approximately 95% of the neurons of the adult neocortex. During this time, the cell cycle of the PVE population lengthens from about 8 h to over 18 h and the progenitor population passes through a total of 11 cell cycles. This 6-day, 11-cell cycle period comprises the "neuronogenetic interval" (NI). At each passage through the cell cycle, the proportion of daughter cells that exit the cell cycle (Q cells) increases from 0 at the onset of the NI to 1 at the end of the NI. The proportion of daughter cells that re-enter the cell cycle (P cells) changes in a complementary fashion from 1 at the onset of the NI to 0 at the end of the NI. This set of systematic changes in the cell cycle and the output from the proliferative population of the PVE allows a quantitative and mathematical treatment of the expansion of the PVE and the growth of the cortical plate that nicely accounts for the observed expansion and growth of the developing neocortex. In addition, we show that the cells produced during a 2-h window of development during specific cell cycles reside in a specific set of laminae in the adult cortex, but that the distributions of the output from consecutive cell cycles overlap. These dynamic events occur in all areas of the PVE underlying the neocortex, but there is a gradient of maturation that begins in the rostrolateral neocortex near the striatotelencephalic junction and which spreads across the surface of the neocortex over a period of 24-36 h. The presence of the gradient across the hemisphere is a possible source of positional information that could be exploited during development to establish the areal borders that characterize the adult neocortex.

Non-NASA Center↗

Improved techniques for thermomechanical testing in support of deformation modeling

The feasibility of generating precise thermomechanical deformation data to support constitutive model development was investigated. Here, the requirement is for experimental data that is free from anomalies caused by less than ideal equipment and procedures. A series of exploratory tests conducted on Hastelloy X showed that generally accepted techniques for strain controlled tests were lacking in at least three areas. Specifically, problems were encountered with specimen stability, thermal strain compensation, and temperature/mechanical strain phasing. The source of these difficulties was identified and improved thermomechanical testing techniques to correct them were developed. These goals were achieved by developing improved procedures for measuring and controlling thermal gradients and by designing a specimen specifically for thermomechanical testing. In addition, innovative control strategies were developed to correctly proportion and phase the thermal and mechanical components of strain. Subsequently, the improved techniques were used to generate deformation data for Hastelloy X over the temperature range, 200 to 1000 C.

Castelli, Michael G.↗

Search for subgrid scale parameterization by projection pursuit regression

The dependence of subgrid-scale stresses on variables of the resolved field is studied using direct numerical simulations of isotropic turbulence, homogeneous shear flow, and channel flow. The projection pursuit algorithm, a promising new regression tool for high-dimensional data, is used to systematically search through a large collection of resolved variables, such as components of the strain rate, vorticity, velocity gradients at neighboring grid points, etc. For the case of isotropic turbulence, the search algorithm recovers the linear dependence on the rate of strain (which is necessary to transfer energy to subgrid scales) but is unable to determine any other more complex relationship. For shear flows, however, new systematic relations beyond eddy viscosity are found. For the homogeneous shear flow, the results suggest that products of the mean rotation rate tensor with both the fluctuating strain rate and fluctuating rotation rate tensors are important quantities in parameterizing the subgrid-scale stresses. A model incorporating these terms is proposed. When evaluated with direct numerical simulation data, this model significantly increases the correlation between the modeled and exact stresses, as compared with the Smagorinsky model. In the case of channel flow, the stresses are found to correlate with products of the fluctuating strain and rotation rate tensors. The mean rates of rotation or strain do not appear to be important in this case, and the model determined for homogeneous shear flow does not perform well when tested with channel flow data. Many questions remain about the physical mechanisms underlying these findings, about possible Reynolds number dependence, and, given the low level of correlations, about their impact on modeling. Nevertheless, demonstration of the existence of causal relations between sgs stresses and large-scale characteristics of turbulent shear flows, in addition to those necessary for energy transfer, provides important insight into the relation between scales in turbulent flows.

Meneveau, C.↗

Friction stir gradient alloying: A high-throughput method to explore the influence of V in enabling HCP to BCC transformation in a γ-FCC dominated high entropy alloy

The compositional possibilities in high entropy alloys (HEAs) is very vast and effective strategies are needed to establish potential alloy chemistries. In this study, friction stir gradient alloying (FSGA), a recently-introduced high-throughput (HT) technique incorporating compositional and microstructural gradients, was used to explore the possibility of introducing a bcc transformation domain in γ-fcc dominated TRIP HEA Fe 38.5 Mn 20 Cr 15 Co 20 Si 5 Cu 1.5 (at.%) by vanadium addition. FSGA has accelerated the process of attaining the composition-microstructure library and the results suggest that vanadium addition of ~1 at.% results in nucleation of α-bcc. Additionally, this is the first observation of ε-hcp to α-bcc transformation in a TRIP HEA and supports the Olson-Cohen model of martensitic transformation (γ-fcc → ε-hcp → α-bcc). ε-hcp is nucleated with a pyramidal orientation by the formation of planar faults on the {111} close-packed planes of γ-fcc, and α-bcc is nucleated in the vicinity of ε-hcp along the {111} γ-fcc trace. There are indications for the formation of basal orientated ε-hcp from intersecting ε-hcp planes orientated for the pyramidal slip. Thus, the study shows the nucleation and growth of the α-bcc phase attributed to the dual effects of chemistry and strain and potential co-existence of all the three phases.

36 MATERIALS SCIENCE↗

Edge delamination of composite laminates subject to combined tension and torsional loading

Delamination is a common failure mode of laminated composite materials. Edge delamination is important since it results in reduced stiffness and strength of the laminate. The tension/torsion load condition is of particular significance to the structural integrity of composite helicopter rotor systems. Material coupons can easily be tested under this type of loading in servo-hydraulic tension/torsion test stands using techniques very similar to those used for the Edge Delamination Tensile Test (EDT) delamination specimen. Edge delamination of specimens loaded in tension was successfully analyzed by several investigators using both classical laminate theory and quasi-three dimensional (Q3D) finite element techniques. The former analysis technique can be used to predict the total strain energy release rate, while the latter technique enables the calculation of the mixed-mode strain energy release rates. The Q3D analysis is very efficient since it produces a three-dimensional solution to a two-dimensional domain. A computer program was developed which generates PATRAN commands to generate the finite element model. PATRAN is a pre- and post-processor which is commonly used with a variety of finite element programs such as MCS/NASTRAN. The program creates a sufficiently dense mesh at the delamination crack tips to support a mixed-mode fracture mechanics analysis. The program creates a coarse mesh in those regions where the gradients in the stress field are low (away from the delamination regions). A transition mesh is defined between these regions. This program is capable of generating a mesh for an arbitrarily oriented matrix crack. This program significantly reduces the modeling time required to generate these finite element meshes, thus providing a realistic tool with which to investigate the tension torsion problem.

Hooper, Steven J.↗

A hybrid calibration approach to Hertz-type contact parameters for discrete element models

This study aims at providing a hybrid calibration framework to estimate Hertz-type contact parameters (particle-scale shear modulus and Poisson ratio) for both two-dimensional and three-dimensional discrete element modelling (DEM). On the basis of statistically isotropic granular packings, a set of analytical formulae between macroscopic material parameters (Young modulus and Poisson ratio) and particle-scale Hertz-type contact parameters for granular systems are derived under small-strain isotropic stress conditions. However, the derived analytical solutions are only estimated values for general models. By viewing each DEM modelling as an implicit mathematical function taking the particle-level parameters as independent variables and employing the derived analytical solutions as the initial input parameters, an automatic iterative scheme is proposed to obtain the calibrated parameters with higher accuracies. Considering highly nonlinear features and discontinuities of the macro-micro relationship in Hertz-based discrete element models, the adaptive moment estimation algorithm is adopted in this study because of its capacity of dealing with noise gradients of cost functions. Here, the proposed method is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. Noticeable improvements in terms of calibration efficiency and accuracy have been made.

Constitutive law↗

Elasto-plastic coupled temperature-displacement finite element analysis of two-dimensional rolling-sliding contact with a translating heat source

The present paper describes a transient translating elastoplastic thermomechanical finite element model to study two-dimensional frictional rolling contact. Frictional two-dimensional contact is simulated by repeatedly translating a nonuniform thermomechanical distribution across the surface of an elastoplastic half space. The half space is represented by a two-dimensional finite element mesh with appropriate boundaries. Calculations are for an elastic-perfectly plastic material and the selected thermophysical properties are assumed to be temperature independent. The paper presents temperature variations, stress and plastic strain distributions and deformations. Residual tensile stresses are observed. The magnitude and depth of these stresses depends on (1) the temperature gradients and (2) the magnitudes of the normal and tangential tractions.

Rubin, C. A.↗

On local approximations of the pressure-strain term in turbulence models

The results of numerical simulations of turbulent channel flows were used to examine the validity of the local approximation of the pressure-strain term in the Reynolds stress transport equation. Outside of the viscous sublayer the local approximation compares very well with the exact pressure strain. This agreement is due, at least in part, to the high correlation between the rapid pressure and its Laplacian, which suggests that only the near parts of the flow contribute to the rapid pressure at a point. In the viscous sublayer the distance over which the mean shear can be considered constant is comparable to the length scale in the normal direction of the correlations of velocity gradients, leading to failure of the local approximation.

Bradshaw, P.↗

Creep response of the lunar crust in mare regions from an analysis of crater deformation

The settling trends of 318 lunar mare craters are compared with predictions of numerical finite-element models in order to determine the creep response of the upper lunar mare crust. No settling is evident in craters smaller than 5 km in diameter. Settling rates of larger craters increase as function of crater size in a manner suggesting a nonlinear lunar creep response corresponding to the power law epsilon = 8.3 x 10 to the minus 34th sigma squared where epsilon is the strain rate and sigma is the differential stress. However, the observed nonlinearity is probably an apparent nonlinearity resulting from the temperature-induced viscosity decrease with depth due to a lunar crustal temperature gradient of 3 C/km and a creep activation energy of 20 kcal/mole. It is concluded that creep in the lunar medium is essentially Newtonian, and that the effective viscosity of the upper lunar mare is (1.6 plus or minus 0.3) x 10 to the 25th poise.

Kunze, A. W. G.↗

Detailed Studies on the Structure and Dynamics of Reacting Dusty Flows at Normal and Microgravity

Dusty reacting flows are of particular interest for a wide range of applications. Inert particles can alter the flammability and extinction limits of a combustible mixture. Reacting particles can release substantial amount of heat and can be used either for power generation or propulsion. Accumulation of combustible particles in air can result in explosions which, for example, can occur in grain elevators, during lumber milling and in mine galleries. Furthermore, inert particles are used as flow velocity markers in reacting flows, and their velocity is measured by non-intrusive laser diagnostic techniques. Despite their importance, dusty reacting flows have been less studied and understood compared to gas phase as well as sprays. The addition of solid particles in a flowing gas stream can lead to strong couplings between the two phases, which can be of dynamic, thermal, and chemical nature. The dynamic coupling between the two phases is caused by the inertia that causes the phases to move with different velocities. Furthermore, gravitational, thermophoretic, photophoretic, electrophoretic, diffusiophoretic, centrifugal, and magnetic forces can be exerted on the particles. In general, magnetic, electrophoretic, centrifugal, photophoretic, and diffusiophoretic can be neglected. On the other hand, thermophoretic forces, caused by steep temperature gradients, can be important. The gravitational forces are almost always present and can affect the dynamic response of large particles. Understanding and quantifying the chemical coupling between two phases is a challenging task. However, all reacting particles begin this process as inert particles, and they must be heated before they participate in the combustion process. Thus, one must first understand the interactions of inert particles in a combustion environment. The in-detail understanding of the dynamics and structure of dusty flows can be only advanced by considering simple flow geometries such as the opposed-jet, stagnation-type. In such configurations the imposed strain rate is well characterized, and the in-depth understanding of the details of the physico-chemical processes can be systematically obtained. A number of computational and experimental studies on spray and particle flows have been conducted in stagnation-type configurations. Numerically, the need for a hybrid Eulerian-Lagrangian approach has been identified by Continillo and Sirignano, and the use of such approach has allowed for the prediction of the phenomenon of droplet flow reversal. Gomez and Rosner have conducted a detailed study on the particle response in the opposed-jet configuration, and the particle thermophoretic diffusivities were determined experimentally. Sung, Law and co-workers have conducted numerical studies on the effect of strain rate and temperature gradients on the dynamics of inert particles, as a way of understanding potential errors in experimental LDV data that may arise from thermophoretic forces. This investigation is a combined experimental and numerical study on the details of reacting dusty flows. The specific tasks are: (1) Experimental determination of laminar flame speeds, and extinction strain rates of dusty flows at normal- and micro-gravity as functions of the particle type, particle initial diameter, particle initial number density, and gas phase chemical composition; (2) Detailed numerical simulation of the experiments. Results are compared with experiments and the adequacy of theoretical models is assessed; and (3) Provision of enhanced insight into the thermo-chemical coupling between the two phases.

Andac, M. Gurhan↗

Porous Media Approach for Modeling Closed Cell Foam

In order to minimize boil off of the liquid oxygen and liquid hydrogen and to prevent the formation of ice on its exterior surface, the Space Shuttle External Tank (ET) is insulated using various low-density, closed-cell polymeric foams. Improved analysis methods for these foam materials are needed to predict the foam structural response and to help identify the foam fracture behavior in order to help minimize foam shedding occurrences. This presentation describes a continuum based approach to modeling the foam thermo-mechanical behavior that accounts for the cellular nature of the material and explicitly addresses the effect of the internal cell gas pressure. A porous media approach is implemented in a finite element frame work to model the mechanical behavior of the closed cell foam. The ABAQUS general purpose finite element program is used to simulate the continuum behavior of the foam. The soil mechanics element is implemented to account for the cell internal pressure and its effect on the stress and strain fields. The pressure variation inside the closed cells is calculated using the ideal gas laws. The soil mechanics element is compatible with an orthotropic materials model to capture the different behavior between the rise and in-plane directions of the foam. The porous media approach is applied to model the foam thermal strain and calculate the foam effective coefficient of thermal expansion. The calculated foam coefficients of thermal expansion were able to simulate the measured thermal strain during heat up from cryogenic temperature to room temperature in vacuum. The porous media approach was applied to an insulated substrate with one inch foam and compared to a simple elastic solution without pore pressure. The porous media approach is also applied to model the foam mechanical behavior during subscale laboratory experiments. In this test, a foam layer sprayed on a metal substrate is subjected to a temperature variation while the metal substrate is stretched to simulate the structural response of the tank during operation. The thermal expansion mismatch between the foam and the metal substrate and the thermal gradient in the foam layer causes high tensile stresses near the metal/foam interface that can lead to delamination.

Ghosn, Louis J.↗

Reynolds stress development in pressure-driven three-dimensional turbulent boundary layers

The development of the Reynolds stress field was studied for flows in which an initially two-dimensional boundary layer was skewed sideways by a spanwise pressure gradient ahead of an upstream-facing wedge. Two different wedges were used, providing a variation in the boundary-layer skewing. Measurements of all components of the Reynolds stress tensor and all ten triple products were measured using a rotatable cross-wire anemometer. The results show the expected lag of the shear stress vector behind the strain rate. Comparison of the two present experiments with previous data suggests that the lag can be estimated if the radius of curvature of the free-stream streamline is known. The magnitude of the shear stress vector in the plane of the wall is seen to decrease rapidly as the boundary-layer skewing increases. The amount of decrease is apparently related to the skewing angle between the wall and the free stream. The triple products evolve rapidly and profiles in the three-dimensional boundary layer are considerably different than two-dimensional profiles, leaving little hope for gradient transport models for the Reynolds stresses. The simplified model presented by Rotta (1979) performs reasonably well providing that an appropriate value of the T-parameter is chosen.

Anderson, Shawn D.↗

Self-equilibrated backstresses induce compensation between hardening and softening: Micromechanical and microstructural features

For crystal plasticity formulation neglecting the strain gradients at the dislocation scale, backstresses are necessary to satisfy the mesoscopic stress equilibrium and strain compatibility requirements. Backstresses in a material point and those between neighboring points are treated separately. Physically backstresses can dynamically induce either hardening or softening effects to adjust the dislocation hardening of slip systems, whereas the sum of the stresses resulting in backstresses in the representative volume element (RVE) is zero. Backstress hardening suppresses certain slip activities, while backstress softening triggers excess dislocations. A two-phase brass is studied. Backstresses make the hard phase deform smaller and the soft phase deform more than those without employing backstresses in a model. In the former, more homogeneous lattice strain and higher yield stress are obtained. Backstress has limited effects on the evolution of the critical resolved shear stress and crystallographic texture. Furthermore, it is a major cause of the banded structures featuring sharp crystallographic contrasts which are not parallel to the slip planes. Tests indicate that backstress induces a serrated behavior at the early stages of deformation and induces the geometrically necessary bands (GNBs), which are parallel to the macroscopic axes.

36 MATERIALS SCIENCE↗

A gradient in the duration of the G1 phase in the murine neocortical proliferative epithelium

Neuronogenesis in the neocortical pseudostratified ventricular epithelium (PVE) is initiated rostrolaterally and progresses caudo-medially as development progresses. Here we have measured the cytokinetic parameters and the fractional neuronal output parameter, Q, of laterally located early-maturing regions over the principal embryonic days (E12-E15) of neocortical neuronogenesis in the mouse. These measures are compared with ones previously made of a medial, late-maturing portion of the PVE. Laterally, as medially, the duration of the neuronogenetic interval is 6 days and comprises 11 integer cell cycles. Also, in both lateral and medial areas the length of G1 phase (TG1) increases nearly 4-fold and is the only cell cycle parameter to change. Q progresses essentially identically laterally and medially with respect to the succession of integer cell cycles. Most importantly, from E12 to E13 there is a steeply declining lateral to medial gradient in TG1. The gradient is due both to the lateral to medial graded stage of neuronogenesis and to the stepwise increase in TG1 with each integer cycle during the neuronogenetic interval. To our knowledge this gradient in TG1 of the cerebral PVE is the first cell biological gradient to be demonstrated experimentally in such an extensive proliferative epithelial sheet. We suggest that this gradient in TG1 is the cellular mechanism for positionally encoding a protomap of the neocortex within the PVE.

Non-NASA Center↗

Validation of Hole-Drilling Residual Stress Measurements in Workpieces of Various Thickness

A recent revision to the ASTM E837 standard for near-surface residual stress measurement by the hole-drilling method describes a new thickness-dependent stress calculation procedure applicable to “thin” and “intermediate” workpieces for which strain versus depth response depends on workpiece thickness. This new calculation procedure differs from that of the prior standard, which applies only to thick workpieces with strain versus depth response independent of thickness. Herein we assess the new calculation procedures by performing hole-drilling residual stress measurements in samples with a range of thickness. Near-surface residual stress is measured in a thick aluminum plate containing near-surface residual stress from a uniform shot peening treatment, and in samples of different thickness removed from the plate at the peened surface. A finite element (FE) model is used to assess consistency between measured residual stress across the range of sample thickness. Measured residual stress varies with sample thickness, with thinner samples exhibiting smaller near-surface compressive stress and a larger gradient of subsurface stress. These trends are consistent with both observed bending (curvature) of the removed samples and the trend in FE-calculated expected residual stress. The measured and expected residual stresses are in good agreement for samples of intermediate thickness, but the agreement decreases with sample thickness. Measured residual stress is invariant with gage circle diameter. The new thickness-dependent stress calculation procedure for hole-drilling provides meaningful improvement compared to thick-workpiece calculations.

42 ENGINEERING↗

Measurement of the Turbulence Kinetic Energy Budget of a Turbulent Planar Wake Flow in Pressure Gradients

Turbulence kinetic energy (TKE) is a very important quantity for turbulence modeling and the budget of this quantity in its transport equation can provide insight into the flow physics. Turbulence kinetic energy budget measurements were conducted for a symmetric turbulent wake flow subjected to constant zero, favorable and adverse pressure gradients in year-three of research effort. The purpose of this study is to clarify the flow physics issues underlying the demonstrated influence of pressure gradient on wake development and provide experimental support for turbulence modeling. To ensure the reliability of these notoriously difficult measurements, the experimental procedure was carefully designed on the basis of an uncertainty analysis. Four different approaches, based on an isotropic turbulence assumption, a locally axisymmetric homogeneous turbulence assumption, a semi-isotropy assumption and a forced balance of the TKE equation, were applied for the estimate of the dissipation term. The pressure transport term is obtained from a forced balance of the turbulence kinetic energy equation. This report will present the results of the turbulence kinetic energy budget measurement and discuss their implication on the development of strained turbulent wakes.

Liu, Xiao-Feng↗

Effective drift velocity from turbulent transport by vorticity

Here, we highlight the differing roles of vorticity and strain in the transport of coarse-grained scalars at length scales larger than ℓ by smaller scale (subscale) turbulence. We use the first term in a multiscale gradient expansion due to which exhibits excellent correlation with the exact subscale physics when the partitioning length ℓ is any scale smaller than that of the spectral peak. We show that unlike subscale strain, which acts as an anisotropic diffusion/antidiffusion tensor, subscale vorticity's contribution is solely a conservative advection of coarse-grained quantities by an eddy-induced nondivergent velocity, v * , that is proportional to the curl of vorticity. Therefore, material (Lagrangian) advection of coarse-grained quantities is accomplished not by the coarse-grained flow velocity, u¯ ℓ , but by the effective velocity, u¯ ℓ + v * , the physics of which may improve commonly used LES models.

42 ENGINEERING↗

Use of photostress to characterize the mechanical behavior of weldments

Welded aluminum is an important part of many space structures. Knowledge of the properties and behavior of weld material and the material surrounding the weld is important for modeling and design of the structures. Photoelastic coatings (Photostress) and strain gages were used to determine behavior of heat treated and as welded joints made from 2219T87 parent material and 2319 weld material subjected to tensile loads. TIG welds of 1/8, 1/2, and 1.4 inches thickness were investigated. Discontinuous yielding was observed in all tests and highly non-uniform behavior through the weld thickness was observed in joints having welds 1.4 inches thick. Joints having welds 1/8 and 1/2 inches thick had only small differences in behavior through the thickness of the weld. Joints in the 1/2 inch thick material contained distinct zones of constant strain within the normal strain gradient extending outward from the weld centerline. These zones had different thickness and locations. Points at the weld centerline, and for a distance of nearly one inch from the centerline, exhibited very nonlinear behavior during the first loading but exhibited near perfect strain hardening during the second loading.

Gambrell, S. C., Jr.↗