Search NASA⌕ Search

SEARCH · Search NASA

Results for “Deformation Gradient”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Damaged stiffened shell research at NASA. Langley Research Center

The structural behavior of stiffened fuselage shells with local cracks is influenced by the local deformations and stress gradients near the cracks, and this behavior must be understood to predict fuselage structural integrity and residual strength. As a crack grows, the stiffness and internal load distributions in a stiffened shell change, affecting the local deformations and stress gradients near the crack. Transport fuselage shells are designed to support internal pressure and mechanical loads which cause a geometrically nonlinear response. Research is being conducted at NASA LaRC to develop and demonstrate a verified nonlinear stiffened shell analysis methodology for accurately predicting the global and local behavior of nonlinear stiffened fuselage shells with local cracks and with combined internal pressure and mechanical loads. The nonlinear structural analysis methodology being developed and planned at NASA LaRC for damaged stiffened fuselage shells and the experiments being planned to verify this analysis methodology are summarized in the present paper.

Aging Aircraft↗
DOE OSTI2022

Atomistic simulations of nanoindentation on nanoglasses: Effects of grain size and gradient microstructure on the mechanical properties

The use of nanoglass (NG) microstructures has been shown as an effective strategy to improve the ductility of rather brittle metallic glass alloys. To evaluate the effects of grain size and gradient design on the mechanical properties of NGs, we perform molecular dynamics simulations of nanoindentation on Cu 64 Zr 36 NGs. We consider samples with uniform 3 and 7 nm grain sizes as well as gradient microstructures with gradient grain sizes varying from 3 to 7 nm. Here, the results show that the deformation mechanism in NGs with small grain sizes is dominated by the activation and evolution of multiple shear transformation zones. Increasing grain sizes enhances the elastic modulus and hardness at the cost of reduced abrasion resistance. Though the average grain size at the indentation surface plays a crucial role in the deformation behavior of gradient NG models, their strain localization and plastic deformation states are affected by the grain sizes far from the indenter. Desired combinations of mechanical properties can be realized by different NG microstructure designs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Eulerian thermo-mechanical simulations of heterogeneous solid propellants using an approximate projection method

We present a mathematical model that describes thermo-mechanical deformations and thermal gradients on the unsteady burning of a heterogeneous solid propellant. A scaling study shows that the deformations in solid at combustion timescales can be treated as quasi-static. The resulting thermo-mechanical formulation is formulated on a Cartesian grid and makes use of a weak form of Chorin-type projection method to deal with large difference in shear modulus of constituent materials. Further, a one-dimensional verification study is carried out by comparing numerical simulations with those of an analytical model. In addition, convergence studies for a two-dimensional propellant sandwich configuration are presented for the stress, velocity, and reference map components. Finally, simulations are carried out for a two-dimensional random propellant pack and the time-averaged burn rate and strains are reported.

42 ENGINEERING↗
DOE OSTI2020

On the use of adjoints in the inversion of observed quasi-static deformation

An adjoint-based conjugate gradient algorithm provides an efficient means for imaging sources of deformation within the Earth, such as volume stresses associated with fluid flow in aquifers and reservoirs. For time intervals over which the overburden deforms elastically, one can calculate the gradient elements for a single model update using just two numerical simulations. The first is a forward run that is used to compute the residuals associated with the given iteration. The second simulation is to evaluate the application of the adjoint operator to the residuals. In this adjoint calculation, the residual displacements are applied as sources at the measurement locations, driving the deformation in the simulation. The volume stress on the source grid blocks, in response to the residual displacements, provide the gradient components. We apply this technique to satellite-based interferometric synthetic aperture radar (InSAR) line-of-sight displacements that were observed over an oil reservoir in California’s Central Valley. Furthermore, we find that the adjoint-based gradient estimates, requiring 18 CPU seconds, agree with conventional numerical calculations that take over 3700 CPU seconds to compute. Conjugate gradient algorithms utilizing the conventional approach and adjoint-based gradient computations give roughly the same reductions in misfit and similar final estimates of reservoir volume change.

58 GEOSCIENCES↗
DOE OSTI2022

Sand Migration Simulation during Gas Production from Gas Hydrate Reservoir at Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska

Uncontrolled sand production impedes continuous gas production from a hydrate reservoir as observed in the past field-scale gas production tests. Sand mobilization is strongly linked with sediment deformation and high pressure gradient. Sand production occurs when the mobilized sands reach the well. Throughout gas production from a hydrate reservoir, because of hydrate dissociation, deformation and pressure gradient change both in time and space. Therefore, it is necessary to consider the entire process to identify where sand is likely mobilized and how much mobilized sand could reach the well. Here, this study utilizes a coupled thermal, hydrological, chemical, and mechanical (thermo-hydro-chemo-mechanical) sand migration model to simulate a year long gas production from the hydrate reservoir at the Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska. It is found that sand would mainly come from the lower portion of the production zone where faster hydrate dissociation occurs. The relatively faster hydrate dissociation coupled with the fact that hydrate-bearing sediments and hydrate-free sediments have similar stiffness but different strengths causes complex stress transfer, which results in excessive shear deformation when the upper portion of the production zone starts to dissociate. This is evident not only near the well but also away from the well, leading to a large amount of sand mobilization. Another focus of the modeling study is evaluation of the effect of sand migration on gas production. The comparison of two extreme cases–with sandscreen and without sandscreen–suggests that migrated and settled solids around a sand control device can prevent pressure drop across the near-wellbore zone and lead to reduction of the gas production rate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Deformation accommodating periodic computational domain for a uniform velocity gradient

Many multiscale methods for granular materials use periodic computational domains to consider particle scale interactions and then calculate the stress to drive the continuum scale calculations. For problems involving large material deformations, the computation domain often needs to be reinitialized because of the distortion, causing the loss of the history information of the system. This work introduces an algorithm to accommodate a large deformation of the material while maintaining the computational domain cuboid to avoid the domain reinitialization during the computation. The algorithm uses a rotating frame of reference, in which the velocity gradient can be represented by an upper triangular matrix. The deformation caused by the upper triangular matrix is treated by the image system implied in the periodicity to maintain the computational domain cuboid. The effect from the rotation of the reference frame is considered using the inertial forces. Finally, simulations of simple and pure shear motions are carried out to illustrate the algorithm.

42 ENGINEERING↗
NASA NTRS1986

The development of an environmental disturbance model for large space structures after the onset of thermal shock

Expressions for the temperature response across a thin large space structure during and immediately following the onset of thermal shock are obtained as a function of the properties of the material and the solar incidence angle. The thermal gradients induced due to solar radiation heating result in thermal deformation. From the temperature gradient information, expressions are developed for the time history of the thermally induced deflections based on the thermoelastic relationships. The effect of the variation of the solar incidence angle and the variation of the emissivity of the surface is considered. From the deflection time history, a model of the disturbance moments is developed and evaluated for different solar incidence angles.

Hamsath, N.↗
NASA NTRS2002

A geometric deformable model for echocardiographic image segmentation

Gradient vector flow (GVF), an elegant external force for parametric deformable models, can capture object boundaries from both sides. A new geometric deformable model is proposed that combines GVF and the geodesic active contour model. The level set method is used as the numerical method of this model. The model is applied for echocardiographic image segmentation.

NASA Program Biomedical Research and Countermeasur↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
NASA NTRS1992

Progress in integrated analysis with adaptive unstructured meshing

Design of lightweight structures and thermal protection systems for hypersonic vehicles depend on accurate prediction of aerothermal loads, structural temperatures and their gradients, and structural deformations and stresses. Concentration is on an alternative meshing technique which generates an entirely new adaptive unstructured mesh based on the solution obtained from the earlier mesh. The technique combined with the finite element method has been shown to significantly improve the efficiency and accuracy of the fluid, thermal, and structural analyses. Current capability of the adaptive unstructured meshing technique for the integrated fluid-thermal-structural analysis is described first. The technique was extended to transient thermal analysis of structures with time-dependent adaptive meshing to capture the detailed temperature response with a minimum number of unknowns and computational cost. Both linear and higher-order finite elements are implemented to demonstrate the generality of the technique and to investigate their solution accuracy. Currently, the adaptive meshing technique is being developed for plane structures that can be modeled with membrane elements and built-up structures modeled with membrane and bending elements. The capability of the technique to these different disciplinary problems is demonstrated by several examples.

Dechaumphai, Pramote↗
DOE OSTI2021

Constraint and size effects in confined layer plasticity

Plane strain finite element analyses are used to model the experiments for a thin metal layer confined between elastic solids. Here, the thin metal layer undergoes elastic-plastic deformations, here modelled by strain gradient plasticity, while the rest of the solid deforms only elastically. Plane strain is assumed and finite strains are accounted for. The microscopic boundary condition at the band interface is initially taken to be zero plastic strain increment, but plastic straining at the interface is allowed when the plastic strain gradient at the interface has exceeded a critical value. When the narrow metal band is inclined 45° relative to the direction of compression of the solid simple shear develops in the band and predictions are rather similar to those of a one-dimensional analysis for an infinite band. For a band at 90° the elastic parts give significant constraint on plastic flow, the stress fields are non-uniform along the band and the average nominal compressive stress shows strong sensitivity to the layer thickness and to the width of the specimen. Remarkably good qualitative, and in some cases quantitative, agreement is found with the experimental observations. Considering also other band angles, both the average nominal shear and compressive stress show variations with the angle.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
DOE OSTI2023

Gradient nanostructured steel with superior tensile plasticity

Nanostructured metallic materials with abundant high-angle grain boundaries exhibit high strength and good radiation resistance. While the nanoscale grains induce high strength, they also degrade tensile ductility. We show that a gradient nanostructured ferritic steel exhibits simultaneous improvement in yield strength by 36% and uniform elongation by 50% compared to the homogenously structured counterpart. In situ tension studies coupled with electron backscattered diffraction analyses reveal intricate coordinated deformation mechanisms in the gradient structures. The outermost nanolaminate grains sustain a substantial plastic strain via a profound deformation mechanism involving prominent grain reorientation. This synergistic plastic co-deformation process alters the rupture mode in the post-necking regime, thus delaying the onset of fracture. The present discovery highlights the intrinsic plasticity of nanolaminate grains and their significance in simultaneous improvement of strength and tensile ductility of structural metallic materials.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Stresses and deformations in composite tubes due to a circumferential temperature gradient

A linear elasticity solution for determining the response of composite tubes subjected to a circumferential temperature gradient is presented. Numerical examples are used to show that, in a single layer tube, fiber orientation strongly influences response. When the fibers are aligned axially, all stress components in the tube are small. When the fibers are aligned circumferentially, the hoop stress becomes large. This difference in behavior is due to the large difference between the radial and circumferential coefficients of thermal expansion when the fibers are oriented circumferentially. In multilayer tubes, stresses are quite high and just two constants characterize the overall bending and axial deformations of the tubes.

Hyer, M. W.↗
NASA NTRS2000

A Dynamical Model for The Jupiter Thermocline

A new conceptual model for the observed banded wind structure of the giant outer planets is proposed, assuming a smoothly distributed potential vorticity (PV) atop a layer of deeply seated, latitudinally variable stratification. With a vanishing horizontal entropy contrast presumably enforced at the bottom of the flow layer by the underlying convection, the thermal-wind balance of observed cloud-top motions implies a mapping of constant potential temperature surfaces, nearly vertical at upper tropospheric levels, down to a deeper, flatter, but variably sloped 'thermocline.' For a generally colder-poleward distribution of isentropes consistent with the strong equatorial jets of both Jupiter and Saturn, this kind of mapping would also imply a poleward decreasing static stability. The proposed temperature-stability distribution is just the reverse of the Earth's troposphere, where the strongest latitudinal potential temperature gradients are at the bottom instead of the top, with static stability generally increasing toward the pole. The warmer-stabler correlation for Jupiter would be consistent, however, with a nearly monotonic distribution of PV from low to high latitudes, over planetary scales for which the planetary vorticity dominates the relative vorticity of the jets. In this way efficient PV mixing for the isentropically bounded thermocline can account for the dynamical maintenance of the cyclonic flanks of the equatorial jet. Over latitudinal intervals comparable to the internal deformation scale, however, the gradient of the absolute vorticity is dominated by the flow curvature, with the local latitudinal gradient of the static stability proportional to the local potential temperature gradient and therefore the local variation in geostrophic velocity. The resulting correlation of velocities and vorticity gradients implies an alternation of the flow over the internal deformation radius set by the vertical entropy contrast.

Allison, Michael↗
NASA NTRS2005

Tools for Designing and Analyzing Structures

Structural Design and Analysis Toolset is a collection of approximately 26 Microsoft Excel spreadsheet programs, each of which performs calculations within a different subdiscipline of structural design and analysis. These programs present input and output data in user-friendly, menu-driven formats. Although these programs cannot solve complex cases like those treated by larger finite element codes, these programs do yield quick solutions to numerous common problems more rapidly than the finite element codes, thereby making it possible to quickly perform multiple preliminary analyses - e.g., to establish approximate limits prior to detailed analyses by the larger finite element codes. These programs perform different types of calculations, as follows: 1. determination of geometric properties for a variety of standard structural components; 2. analysis of static, vibrational, and thermal- gradient loads and deflections in certain structures (mostly beams and, in the case of thermal-gradients, mirrors); 3. kinetic energies of fans; 4. detailed analysis of stress and buckling in beams, plates, columns, and a variety of shell structures; and 5. temperature dependent properties of materials, including figures of merit that characterize strength, stiffness, and deformation response to thermal gradients

Luz, Paul L.↗
DOE OSTI2021

Global and local texture development during initial plastic deformation of cold-pilgered Zircaloy-4 tubing

In this work we present a thorough characterization of the texture development during the initial step of a cold pilgering process of a Zircaloy-4 TREX tube, on which the axial cross section is reduced by 80% over a 40cm length with an overall Q-ratio of 2.34, followed by a recrystallization thermal treatment. Texture across a tapered specimen representing this working zone was characterized by neutron diffraction and high energy synchrotron x-ray diffraction at 28 positions (4 locations along its length and 7 locations across its thickness with 0.3 mm spatial resolution), and on the deformed end after the thermal treatment. Results show that the starting material (TREX tube) presents a texture gradient through its wall thickness, with the c poles aligned to the hoop direction at the outer surface (with radial and hoop Kearns factors: f R =0.36, f H = 0.50), and gradually rotating around the axial direction to become a quasi-radial orientation (~35 o out of the radial direction) at the inner surface ( f R =0.43, f H =0.35). The overall texture becomes more radial due to the cold pilgering, but the texture gradient remains during the whole deformation process, being always more radial at the inner surface, ( f R =0.53, f H =0.37 at 80%) than at the outer surface ( f R =0.41, f H =0.50at 80%). The most remarkable changes in texture produced by the rolling process is a 30 o rotation of the grains around their c-axis, producing a marked increase of the <11 2 ¯ 0> // axial fiber texture at expenses of the original <10 1 ¯ 0> // axial fiber texture observed in the TREX tube. The observed changes in texture during deformation were correctly explained by visco-plastic self-consistent (VPSC) modelling. This was achieved after the inclusion of the basal <a> slip system, which plays a central role to explain the 30 o rotation around the c-axis during deformation. The texture evolution was correctly described after the inclusion of a drag effect to grain rotation due to the interaction with neighbouring grains. The through thickness texture differences are greatly reduced by the recrystallization thermal treatment, with similar Kearns factors at the inner surface ( f R =0.50, f H =0.28) and the outer surface ( f R =0.48, f H =0.39). The predominance of the <10 1 ¯ 0> // axial fiber texture is recovered after recrystallization.

36 MATERIALS SCIENCE↗