Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Deformation”

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 91 records · Page 5

The influence of kinematics of deformation on polycrystalline halite dynamic recrystallization: Full-field simulation of simple shear versus pure shear

Rock salt, composed mainly of halite, flows viscoplastically over a wide range of geological conditions, strongly impacting the dynamic evolution of sedimentary basins and orogens. Understanding how dislocation creep, which involves dislocation glide, intracrystalline recovery and dynamic recrystallization, influences the microstructure and rheology of halite under various deformation kinematics and temperatures is crucial for enhancing knowledge of salt flow dynamics. Here, this study employs a full-field numerical simulation method to compare the viscoplastic deformation of polycrystalline halite under simple shear and pure shear conditions up to a natural strain of ε = 1.5 at temperatures ranging from 100 °C to 300 °C. The results are presented in terms of crystallographic preferred orientation (CPO), grain shape preferred orientation (SPO), subgrain boundary direction, grain size and strain rate distribution. The results indicate that the crystallographic anisotropy of individual halite crystals is transferred to the polycrystalline scale, resulting in strain localization, particularly in simple shear simulations. The kinematics of deformation affect the evolution and distribution of high strain-rate bands, determining the direction of intragranular substructures and the morphology of strain-induced grain boundaries, with minimal impact on grain size. The intensity of grain boundary migration increases with temperature, significantly influencing grain morphology and size, thereby obscuring strain localization, while having little effect on CPOs. At low strain (ε < 1.0), CPOs relative to both the maximum shortening direction and the grain SPO are similar regardless of the deformation kinematics. At high strain (ε > 1.0), simple shear CPOs exhibit three stronger {100} maxima with a monoclinic symmetry relative to the grain SPO compared to the six {100} maxima with an orthotropic symmetry relative to the grain SPO generated under pure shear. Therefore, microstructures and CPOs can serve as indicators of the strain path in polycrystalline halite under various conditions, aiding in determining the shear sense and elucidating the deformation kinematics of salt structures.

58 GEOSCIENCES↗

Ductile Glassy Polymer Networks Capable of Large Plastic Deformation and Heat-Induced Elastic Recovery

Many thermoplastic polymers are ductile by combining strength and large deformations. These deformations are irreversible - known as plastic deformation. Elastomers can deform reversibly but have low strength. To this end, we developed glassy and ductile polyamide networks capable of large plastic deformation (>200% strain) and high strength (~50 MPa tensile strength and ~1500 MPa Young’s modulus), similar to those of polyolefins and Nylon-66. We discovered that hydrogen bonding between meta-phthalamide groups was essential to the ductility. Since these polyamide networks are covalently bonded, we demonstrated their unique durability by repeatable elastic recovery at elevated temperatures, exhibiting indifferent tensile properties in each cycle. Furthermore, when we fixed the strain during the elastic recovery, these polyamide networks actuated stresses of 9–18 MPa, among the highest reported in shape-memory polymer actuators. Here, we envision these ductile, glassy polymer networks as promising alternatives to ductile thermoplastics, given the combined benefit of ductility and durability.

36 MATERIALS SCIENCE↗

Surface deformation coupled with self-organized pattern on a liquid anode of 1 atm DC glow discharge

Intricate, self-organized plasma structures observed above the surface of a liquid anode of atmospheric DC glow discharge were found to give rise to coherent, organized surface deformation and mechanical wave formation at the plasma–liquid interface. This new phenomenon indicates that the liquid is closely coupled to the plasma by the anode sheath’s electrohydrodynamic (EHD) force. A scientific question then arises: Do surface perturbations, coupled with the nonuniform surface charge distribution, enhance the electric field and induce self-organization? Using the reflective background-oriented schlieren technique, the liquid surface profile under the plasma pattern was measured for the first time. The results show that surface distortions are driven by the repulsive Coulomb force of nonuniform net-negative surface charge acted by the anode sheath field. The impacts of various operating parameters on the patterns and surface waves were examined, revealing the significance of gas heating and liquid charge relaxation time in the pattern formation mechanism. Time-resolved dynamics of a pulsed DC discharge indicated that the surface deformation only appeared after the establishment of plasma patterns. Statistically, the surface wave under the plasma has high wave numbers (8000–16000 m −1 ) and small amplitudes ($<$ 10 µm), generally found in the capillary wave regime. Yet the motion of surface deformations is in tandem with the plasma pattern and exhibits a nondispersive nature of constant phase velocity (0.1–0.4 m s −1 ), suggesting the dominant role of EHD force over the surface tension in the observed surface wave. These results indicate that the deformed liquid surface is driven by the EHD force of the plasma sheath. Although they share a similar geometry, the deformation and wave dynamics of the liquid surface do not stimulate a plasma pattern. Importantly, the complex EHD coupling in the plasma–liquid system raises awareness and new challenges for plasma control engineering, and the quantitative characteristics of the nondispersive surface wave are informative for advancing relevant theory and modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coulomb Branch Amplitudes from a Deformed Amplituhedron Geometry

The amplituhedron provides, via geometric means, the all-loop integrand of scattering amplitudes in maximally supersymmetric Yang-Mills theory. Unfortunately, dimensional regularization, used conventionally for integration, breaks the beautiful geometric picture. This motivates us to propose a “deformed” amplituhedron. Focusing on the four-particle amplitude, we introduce two deformation parameters, which can be interpreted as particle masses. We provide evidence that the mass pattern corresponds to a specific choice of vacuum expectation values on the Coulomb branch. The deformed amplitude is infrared finite, making the answer well defined in four dimensions. Leveraging four-dimensional integration techniques based on differential equations, we compute the amplitude up to two loops. In the limit where the deformation parameters are taken to zero, we recover the known Bern-Dixon-Smirnov amplitude. In the limit where only one deformation parameter is taken to zero, we find a connection to the angle-dependent cusp anomalous dimension. Published by the American Physical Society 2024

Physics↗

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗

Effect of propellant deformation on ignition and combustion processes in solid propellant cracks

A comprehensive theoretical model was formulated to study the development of convective burning in a solid propellant crack which continually deforms due to burning and pressure loading. In the theoretical model, the effect of interrelated structural deformation and combustion processes was taken into account by considering (1) transient, one dimensional mass, momentum, and energy conservation equations in the gas phase; (2) a transient, one dimensional heat conduction equation in the solid phase; and (3) quasi-static deformation of the two dimensional, linear viscoelastic propellant crack caused by pressure loading. Partial closures may generate substantial local pressure peaks along the crack, implying a strong coupling between chamber pressurization, crack combustion, and propellant deformation, especially when the cracks are narrow and the chamber pressurization rates high. The maximum pressure in the crack cavity is generally higher than that in the chamber. The initial flame-spreading process is not affected by propellant deformation.

Kumar, M.↗

Crustal deformation, the earthquake cycle, and models of viscoelastic flow in the asthenosphere

The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.

Cohen, S. C.↗

Finite elastic-plastic deformation of polycrystalline metals

Applying Hill's self-consistent method to finite elastic-plastic deformations, the overall moduli of polycrystalline solids are estimated. The model predicts a Bauschinger effect, hardening, and formation of vertex or corner on the yield surface for both microscopically non-hardening and hardening crystals. The changes in the instantaneous moduli with deformation are examined, and their asymptotic behavior, especially in relation to possible localization of deformations, is discussed. An interesting conclusion is that small second-order quantities, such as shape changes of grains and residual stresses (measured relative to the crystal elastic moduli), have a first-order effect on the overall response, as they lead to a loss of the overall stability by localized deformation. The predicted incipience of localization for a uniaxial deformation in two dimensions depends on the initial yield strain, but the orientation of localization is slightly less than 45 deg with respect to the tensile direction, although the numerical instability makes it very difficult to estimate this direction accurately.

Iwakuma, T.↗

Crustal deformation, the earthquake cycle, and models of viscoelastic flow in the asthenosphere

The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.

Cohen, S. C.↗

Geometry of miocene extensional deformation, lower Colorado River Region, Southeastern California and Southwestern Arizona: Evidence for the presence of a regional low-angle normal fault

The geometry of Miocene extensional deformation, which changes along a 120 km-long, northeast-trending transect from the southestern Chocolate Mountains, southeastern California, to the Trigo and southern Dome Rock Mountains, southwestern Arizona is discussed. Based upon regional differences in the structural response to extension and estimated extensional strain, the transet can be divided into three northwesterly-trending structural domains. From southwest to northeast, these domains are: (1) southestern Chocolate-southernmost Trigo Mountains; (2) central to northern Trigo Mountains; and (3) Trigo Peaks-southern Dome Rock Mountains. All structures formed during the deformation are brittle in style; fault rocks are composed of gouge, cohesive gouge, and local microbreccia. In each structural domain, exposed lithologic units are composed of Mesozoic crystalline rocks unconformably overlain by Oligocene to Early Miocene volcanic and minor interbedded sedimentary rocks. Breccia, conglomerate, and sandstone deposited synchronously with regional extension locally overlie the volcanic rocks. Extensional deformation largely postdated the main phase of volcanic activity, but rare rhyolitic tuff and flows interbedded with the syndeformational clastic rocks suggest that deformation began during the waning stages of valcanism. K-Ar isotopic ages indicate that deformation occurred in Miocene time, between about 22 and m.y. ago.

Tosdal, R. M.↗

Structural deformation of the Leoville chondrite

A petrofabric study of the Leoville chondrite was performed in order to gain an understanding of its deformation process. Strain values were determined from the chondrules' axial ratios or by the all-object-separation method, and the compaction required to deform Leoville chondrules was calculated from the decrease of radius in the circle-to-ellipse transformation. Structural and strain studies indicate the following order of events experienced by Leoville: (1) formation as an accretionary breccia, (2) compaction and development of foliation, (3) fracturing, possibly due to impact on earth, and (4) terrestrial weathering. Similarities between the axial ratios of chondrules (1.9-2.0) and of the whole meteorite (2.0) indicate either that Leoville deformed homogeneously or that it deformed as a framework of touching chondrules. This amount of strain corresponds to approximately 33 percent uniaxial shortening, assuming constant volume, and is likely to have been caused by compaction due to overburden from progressive accretion on the chondrite body. Diffusional flow is believed to be the mechanism by which olivines in these chondrules deformed.

Cain, P. M.↗

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗

Effects of transverse shear deformation on buckling of laminated cylinders as a function of thickness and ply orientation

Buckling loads of thick-walled, orthotropic, simply-supported right circular cylinders are predicted using a new higher-order transverse shear deformation theory. The higher-order theory shows that, by more accurately accounting for transverse shear deformation effects, the predicted buckling load may be reduced by as much as 80 percent compared to predictions based on conventional transverse shear deformation theory. A parametric study of the effect of ply orientation on the buckling load of axially compressed cylinders indicates that laminates containing 0 deg plies are the most sensitive to transverse shear deformation effects. Interaction curves for buckling of cylinders with axial compressive and external pressure loadings indicate that buckling loads due to external pressure loadings are much less sensitive to transverse shear deformation effects than those due to axial compressive loadings.

Jegley, Dawn C.↗

Vector constraints on western U.S. deformation from space geodesy, neotectonics, and plate motions

The rate-of-slip vector on the San Andreas fault in central California estimated from geodetic and Holocene geological data is inconsistent with the prediction of rigid plate models such as RM2. This well-known 'San Andreas discrepancy' is diagnostic of plate deformation distributed both east of the fault in the Basin and Range and west of the fault along the California continental margin. Constraints on the integrated deformation rates across these two regions consistent with: (1) the kinematical boundary conditions imposed by the rigid plate model; (2) neotectonic and paleoseismic estimates of deformation rates; (3) ground-based geodetic measurements; and (4) rates of change observed by very long baseline interferometry along seven baselines to western U.S. sites are constructed. The space-geodetic data on Basin and Range extension taken over a 4-year interval are compatible with geological observations averaged over the Holocene; the best estimate of its integrated deformation rate, provided by the joint inversion of both data types, is 9.7 + or - 2.1 mm/yr, N 56 deg W + or - 10 deg, too small and in the wrong direction to account entirely for the San Andreas discrepancy. The integral of this deformation, estimated by subtracting the Basin and Range contribution from the discrepancy vector, requires significant right-lateral shear parallel to the San Andreas (13 + or - 5 mm/yr) and some compression perpendicular to it (9 + or - 3 mm/yr).

Minster, J. Bernard↗

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗

An algorithmic approach to crustal deformation analysis

In recent years the analysis of crustal deformation measurements has become important as a result of current improvements in geodetic methods and an increasing amount of theoretical and observational data provided by several earth sciences. A first-generation data analysis algorithm which combines a priori information with current geodetic measurements was proposed. Relevant methods which can be used in the algorithm were discussed. Prior information is the unifying feature of this algorithm. Some of the problems which may arise through the use of a priori information in the analysis were indicated and preventive measures were demonstrated. The first step in the algorithm is the optimal design of deformation networks. The second step in the algorithm identifies the descriptive model of the deformation field. The final step in the algorithm is the improved estimation of deformation parameters. Although deformation parameters are estimated in the process of model discrimination, they can further be improved by the use of a priori information about them. According to the proposed algorithm this information must first be tested against the estimates calculated using the sample data only. Null-hypothesis testing procedures were developed for this purpose. Six different estimators which employ a priori information were examined. Emphasis was put on the case when the prior information is wrong and analytical expressions for possible improvements under incompatible prior information were derived.

Iz, Huseyin Baki↗

Gravity-driven deformation of the crust on Venus

Intense horizontal deformation features are observed on the surface of Venus despite the apparent absence of terrestrial-style plate tectonics. High surface temperature, low erosion rate, and regions of elevated terrain suggest that gravity spreading may be responsible for some of these features. Gravity spreading on the Earth is known to cause extensive deformation in certain thin-skinned tectonic regions, where the uppermost layers become detached along a decollement. While the high pore pressures and soft sedimentary rocks frequently observed along decollements on the Earth are not plausible for Venus, temperature induced regions of ductile deformation may exist. Simple viscous models of gravity spreading above a ductile region of the crust are developed to examine the effects of rheology, thermal gradient, surface slope and crustal thickness on deformation. For a reasonable range of parameters, the models show that surface-slope-driven ductile deformation at relatively shallow depths in the crust is possible and may be capable of detaching the uppermost rigid portion of the crust, resulting in horizontal translation.

Smrekar, Suzanne↗