Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deformation estimation”

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 361 records · Page 20

Accurate Realization of GPS Vertical Global Reference Frame

The goal of this project is to improve our current understanding of GPS error sources associated with estimates of radial velocities at global scales. An improvement in the accuracy of radial global velocities would have a very positive impact on a large number of geophysical studies of current general interest such as global sea-level and climate change, coastal hazards, glacial isostatic adjustment, atmospheric and oceanic loading, glaciology and ice mass variability, tectonic deformation and volcanic inflation, and geoid variability. A set of GPS error sources relevant to this project are those related to the combination of the positions and velocities of a set of globally distributed stations as determined &om the analysis of GPS data, including possible methods of combining and defining terrestrial reference frames. This is were our research activities during this reporting period have concentrated. During this reporting period, we have researched two topics: (1) The effect of errors on the GPS satellite antenna models (or lack thereof) on global GPS vertical position and velocity estimates; (2) The effect of reference W e definition and practice on estimates of the geocenter variations.

Elosegui, Pedro↗

Relative motions of the Australian, Pacific and Antarctic plates estimated by the Global Positioning System

Global Positioning System (GPS) measurements spanning approximately 3 years have been used to determine velocities for 7 sites on the Australian, Pacific and Antarctic plates. The site velocities agree with both plate model predictions and other space geodetic techniques. We find no evidence for internal deformation of the interior of the Australian plate. Wellington, New Zealand, located in the Australian-Pacific plate boundary zone, moves 20 +/- 5 mm/yr west-southwest relative to the Australian plate. Its velocity lies midway between the predicted velocities of the two plates. Relative Euler vectors for the Australia-Antarctica and Pacific-Antarctica plates agree within one standard deviation with the NUVEL-1A predictions.

Larson, Kristine M.↗

Adjoint-Based Inversion of Geodetic Data for Sources of Deformation and Strain

An adjoint-based formulation leads to a particularly efficient approach for inverting geodetic measurements for the source of the deformation. Specifically, the quantities necessary to iteratively improve the fit to the observations can be computed with just three forward calculations, one to obtain the current residuals, another to solve the adjoint problem, and a third to compute the step length. An inversion algorithm utilizing the adjoint-based gradient is applied to a set of Interferometric Synthetic Aperture Radar (InSAR) data gathered between 2016 and 2018 over the Tulare Basin in California's Central Valley. Because the measured deformation is due to groundwater withdrawal, a penalty function is included in the inversion to avoid placing aquifer volume change in locations that are far from any documented wells. The solution of the inverse problem provides estimates of aquifer compaction that provide a match to the observed range changes while honoring the well data. The solution indicates an average aquifer volume loss of 2.17 km 3 /year over the two year period from January 2016 to January 2018, encompassing one drought year (2016) and one wet year (2017). Finally, this magnitude of lost volume is compatible with the 3.1 km 3 /year decrease in water volume for the entire Central Valley, estimated from GRACE satellite gravity data.

58 GEOSCIENCES↗

Tuscaloosa Marine Shale Laboratory

The Tuscaloosa Marine Shale (TMS) in Louisiana and Mississippi is an Upper Cretaceous source rock formation sandwiched between the sands of the upper and lower Tuscaloosa sections. The TMS is believed to be the source rock for underlying prolific Tuscaloosa sand formation. The TMS has an unproven estimate of 7,000,000,000 bbls of recoverable oil while its current total average production is about 3,000 bbls of oil per day in 2017. In 2013 and 2014, more than 80 wells were drilled horizontally into the TMS that were fractured using multi-stage fracturing technology. The results from this have been mixed, but recent production for several wells show an appealing initial oil production rate of more than 1000 bbl/day. The preliminary core analysis by industry partners and a few literature studies shows that the TMS is one of the most clay-rich and sensitive shales to water. Due to these and other technical problems, there is high risk for the economic development of TMS compared to other shale plays. The experiences of major industrial players in the TMS show the necessity of open and collaborative efforts to better understand the critical gaps in the development of this challenging and potentially highly economic shale play to enable more cost-efficient and environmentally-sound recovery from this unconventional liquid-rich shale play. The overall objective of this project is to form a consortium of science and industry partners to address the following six major objectives using scientific and technical approaches: 1. To improve wellbore integrity by better understanding the sources of the wellbore instability issues, proposing innovative mud and cement design for the TMS. 2. To improve formation evaluation using laboratory techniques for the evaluation of mineralogical composition, organic content, and produced-water chemistry as well as well log and geophysical analysis. 3. To determine the role of geologic discontinuities on fracture growth and shale creep behavior using digital image correlation technique. 4. To investigate the application of stable CO 2 foam and super-hydrophobic proppants for improved reservoir stimulation. 5. To better understand the nature of water/hydrocarbon/CO 2 flow in clay and organic-rich formation and the role of water/fluid interaction on recovery. 6. To prepare better socio-economic environment for TMS development by community engagement. Subsequently, the TMS virtual laboratory conducted testing and analysis of various properties of rock and formation fluids from the TMS, including but not limited to the following: Analyzing reports and logs to better understand the source of wellbore instability in TMS wells; Experiments to design a customized cement based on TMS requirements; Experiments to obtain the mineralogical and geochemical composition of TMS samples; Seismic analysis of TMS geophysical data to better predict total organic carbon (TOC) content and brittleness in TMS; Well log analysis to better estimate the TOC and geo-mechanical properties of TMS; Experiments on formation water to understand the chemistry of produced water; Experiments to determine the role of lamination and natural fractures on fracture propagation or rock deformation using digital image correlation technique in in-direct tensile tests, semi-circular bend test and creep tests Experiments to determine the stability and rheological properties of nanoparticle-stabilized CO 2 foam in TMS rock samples; Experiments to determine fluid dynamics in un-propped TMS fractures and the role of nano-coating of proppants on fluid dynamics in fractures with proppants; Micro-fluidics experiments to enhance the understanding of fluid dynamics in tight liquidrich pores with high clay content; Socio-economic studies to better engage communities in TMS development.

58 GEOSCIENCES↗

Diamond Smoothing Tools

Diamond smoothing tools have been proposed for use in conjunction with diamond cutting tools that are used in many finish-machining operations. Diamond machining (including finishing) is often used, for example, in fabrication of precise metal mirrors. A diamond smoothing tool according to the proposal would have a smooth spherical surface. For a given finish machining operation, the smoothing tool would be mounted next to the cutting tool. The smoothing tool would slide on the machined surface left behind by the cutting tool, plastically deforming the surface material and thereby reducing the roughness of the surface, closing microcracks and otherwise generally reducing or eliminating microscopic surface and subsurface defects, and increasing the microhardness of the surface layer. It has been estimated that if smoothing tools of this type were used in conjunction with cutting tools on sufficiently precise lathes, it would be possible to reduce the roughness of machined surfaces to as little as 3 nm. A tool according to the proposal would consist of a smoothing insert in a metal holder. The smoothing insert would be made from a diamond/metal functionally graded composite rod preform, which, in turn, would be made by sintering together a bulk single-crystal or polycrystalline diamond, a diamond powder, and a metallic alloy at high pressure. To form the spherical smoothing tip, the diamond end of the preform would be subjected to flat grinding, conical grinding, spherical grinding using diamond wheels, and finally spherical polishing and/or buffing using diamond powders. If the diamond were a single crystal, then it would be crystallographically oriented, relative to the machining motion, to minimize its wear and maximize its hardness. Spherically polished diamonds could also be useful for purposes other than smoothing in finish machining: They would likely also be suitable for use as heat-resistant, wear-resistant, unlubricated sliding-fit bearing inserts.

Voronov, Oleg↗

Estimate Low and High Order Wavefront Using P1640 Calibrator Measurements

P1640 high contrast imaging system on the Palomar 200 inch Telescope consists of an apodized-pupil Lyot coronagraph, the PALM-3000 adaptive optics (P3K-AO), and P1640 Calibrator (CAL). Science images are recorded by an integral field spectrograph covering J-H bands for detecting and characterizing stellar companions. With aberrations from atmosphere corrected by the P3K-AO, instrument performance is limited mainly by the quasi-static speckles due to noncommon path wavefront aberrations for the light to propagate to the P3K-AO wavefront sensor and to the coronagraph mask. The non-common path wavefront aberrations are sensed by CAL, which measures the post-coronagraph E-field using interferometry, and can be effectively corrected by offsetting the P3K-AO deformable mirror target position accordingly. Previously, we have demonstrated using CAL measurements to correct high order wavefront aberrations, which is directly connected to the static speckles in the image plane. Low order wavefront, on the other hand, usually of larger amplitudes, causes light to leak through the coronagraph making the whole image plane brighter. Knowledge error in low order wavefront aberrations can also affect the estimation of the high order wavefront. Even though, CAL is designed to sense efficiently high order wavefront aberrations, the low order wavefront front can be inferred with less sensitivity. Here, we describe our method for estimating both low and high order wavefront aberrations using CAL measurements by propagating the post-coronagraph E-field to a pupil before the coronagraph. We present the results from applying this method to both simulated and experiment data.

aberrations↗

Diffusion and Deformation Mechanism Maps in Li Metal from Atomistic Simulations

The rate of Li transport in the Li metal anode is important for the operation of Li metal-solid state batteries. Here, transport rates due to diffusion and creep are predicted in Li using atomistic simulations. First, molecular dynamics is used to estimate the rate of Li diffusion along dislocations and in grain boundary triple junctions. By combining this data with that from a prior study of grain boundary diffusion the dominant mechanisms and rates of self-diffusion in Li polycrystals are predicted as a function of grain size, grain shape, dislocation density, and temperature. Second, the dominant creep mechanisms are predicted and used to estimate critical current densities and void annihilation times. Grain boundary sliding and coble creep are the dominant mechanisms for micron-sized grains. Lastly, a continuum model for interfacial contact loss reveals that high dislocation densities of ∼10 12 /cm 2 enable achieving battery performance targets for Li grain sizes of ∼10 μm.

Batteries↗

Fracture mechanics validity limits

Fracture behavior is characteristics of a dramatic loss of strength compared to elastic deformation behavior. Fracture parameters have been developed and exhibit a range within which each is valid for predicting growth. Each is limited by the assumptions made in its development: all are defined within a specific context. For example, the stress intensity parameters, K, and the crack driving force, G, are derived using an assumption of linear elasticity. To use K or G, the zone of plasticity must be small as compared to the physical dimensions of the object being loaded. This insures an elastic response, and in this context, K and G will work well. Rice's J-integral has been used beyond the limits imposed on K and G. J requires an assumption of nonlinear elasticity, which is not characteristic of real material behavior, but is thought to be a reasonable approximation if unloading is kept to a minimum. As well, the constraint cannot change dramatically (typically, the crack extension is limited to ten-percent of the initial remaining ligament length). Rice, et al investigated the properties required of J-type parameters, J(sub x), and showed that the time rate, dJ(sub x)/dt, must not be a function of the crack extension rate, da/dt. Ernst devised the modified-J parameter, J(sub M), that meets this criterion. J(sub M) correlates fracture data to much higher crack growth than does J. Ultimately, a limit of the validity of J(sub M) is anticipated, and this has been estimated to be at a crack extension of about 40-percent of the initial remaining ligament length. None of the various parameters can be expected to describe fracture in an environment of gross plasticity, in which case the process is better described by deformation parameters, e.g., stress and strain. In the current study, various schemes to identify the onset of the plasticity-dominated behavior, i.e., the end of fracture mechanics validity, are presented. Each validity limit parameter is developed in detail, and then data is presented and the various schemes for establishing a limit of the validity are compared. The selected limiting parameter is applied to a set of fracture data showing the improvement of correlation gained.

Lambert, Dennis M.↗

The interpretation of crustal dynamics data in terms of plate motions and regional deformation near plate boundaries

During our participation in the NASA Crustal Dynamics Project under NASA contract NAS-27339 and grant NAG5-814 for the period 1982-1991, we published or submitted for publication 30 research papers and 52 abstracts of presentations at scientific meetings. In addition, five M.I.T. Ph.D. students (Eric Bergman, Steven Bratt, Dan Davis, Jeanne Sauber, Anne Sheehan) were supported wholly or in part by this project during their thesis research. Highlights of our research progress during this period include the following: application of geodetic data to determine rates of strain in the Mojave block and in central California and to clarify the relation of such strain to the San Andreas fault and Pacific-North American plate motions; application of geodetic data to infer post seismic deformation associated with large earthquakes in the Imperial Valley, Hebgen Lake, Argentina, and Chile; determination of the state of stress in oceanic lithosphere from a systematic study of the centroid depths and source mechanisms of oceanic intraplate earthquakes; development of models for the state of stress in young oceanic regions arising from the differential cooling of the lithosphere; determination of the depth extent and rupture characteristics of oceanic transform earthquakes; improved determination of earthquake slip vectors in the Gulf of California, an important data set for the estimation of Pacific-North American plate motions; development of models for the state of stress and mechanics of fold-and-thrust belts and accretionary wedges; development of procedures to invert geoid height, residual bathymetry, and differential body wave travel time residuals for lateral variations in the characteristic temperature and bulk composition of the oceanic upper mantle; and initial GPS measurements of crustal deformation associated with the Imperial-Cerro Prieto fault system in southern California and northern Mexico. Full descriptions of the research conducted on these topics may be found in the Semi-Annual status Reports submitted regularly to NASA over the course of this project and in the publications listed.

Solomon, Sean C.↗

An elastoplastic analysis of a uniaxially loaded sheet with an interference-fit bolt

The stresses and strains in a uniaxially loaded sheet with an unloaded interference-fit bolt were calculated by an elastoplastic finite-element analysis. The material properties represented a 7075-T6 aluminum alloy sheet and a steel bolt. The analysis considered the two ideal cases of no slip and no friction at the bolt-sheet interface for a single combination of bolt diameter, interference level, and cyclic loading. When the bolt was inserted, the sheet deformed plastically near the hole; the first tensile load cycle produced additional yielding, but subsequent cycles to the same level caused only elastic cyclic stresses. These stresses together with fatigue data for unnotched specimens were used to estimate crack initiation periods and initiation sites. The cases analyzed with interference-fit bolts were predicted to have crack initiation periods which were about 50 times that for a clearance-fit bolt. Crack initiation was predicted to occur on the transverse axis at a distance of about one radius from the hole.

Crews, J. H., Jr.↗

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗

Plate motions and deformations from geologic and geodetic data

A satellite laser ranging experiment conducted by NASA since 1972 has measured the relative motion between the North America and Pacific plates in California. Based on these measurements, the 896-km distance between San Diego and Quincy, California, is shortening at 62 + or - 9 mm/yr. This geodetic estimate is consistent with the rate of motion between the two plates, calculated from geological data to be 53 + or - 3 mm/yr averaged over the past few million years.

Jordan, T. H.↗

Elastic thermobarometry on metapelites across the crustal section of the Famatinian Arc, Argentina

We applied elastic thermobarometry on garnet-bearing migmatites along two transects through the crustal section at Sierra Valle Fértil-La Huerta, Argentina. We performed quartz-in-garnet barometry and zircon-in-garnet thermometry on metapelites from different paleo-depths across the crustal section. Our work recovers entrapment pressures ranging from 240 to 1330 MPa and entrapment temperatures between 691 and 1574 °C. The entrapment conditions are broadly consistent with anticipated pressures and temperatures along the crustal section derived previously using conventional, thermodynamic thermobarometers. The quartz-in-garnet barometer reproduces those conventionally established entrapment conditions when samples only experienced conditions within the alpha-quartz stability field. Raman-derived pressures for samples that experienced beta-quartz reference conditions are commonly much higher than those established by conventional barometry. Samples that preserve compressive (positive) residual pressures best reproduce reference entrapment pressures. Entrapment temperatures show high variability and overestimation of temperature conditions compared to conventional results. These results indicate elastic thermobarometry over- or under-estimates crystallization conditions in rocks crystallized at high temperatures, as is common in the Famatinian Arc deep-crust. We suggest that modeling quartz behavior across the alpha–beta transition may present challenges, as does shape maturation, viscous deformation, and radiation damage in zircon.

58 GEOSCIENCES↗

Impact of Upheaval Dome, Utah

Reexamination of Upheaval Dome in the Canyonlands National Park, Utah, shows that the structure of this remarkable feature conforms with that expected for a deeply eroded astrobleme. The structure is definitely not compatible with an origin due simply to plastic flowage of salt and other rocks in the underlying Paradox Formation. The most strongly deformed rocks are bounded by a series of circumferential listric faults. The convergent displacement of the rocks corresponds to the deformation that results from collapse of a transient cavity produced by high speed impact. From considerations of the probable depth of exposure of the impact structure and upward extrapolation of the listric faults, the final collapsed crater is estimated to be about 8 to 10 km in diameter; the impacting body was on the order of 0.5 km in diameter.

Shoemaker, E. M.↗

Seismicity, fault plane solutions, depth of faulting, and active tectonics of the Andes of Peru, Ecuador, and southern Colombia

The long-period P waveforms observed for 17 earthquakes in the Peruvian Andes during 1963-1976 are compared with synthetic waveforms to obtain fault-plane solutions and focal depths. The morphological units of the Peruvian Andes are characterized: coastal plains, Cordillera Occidental, altiplano and central high plateau, Cordillera Oriental, and sub-Andes. The data base and analysis methodology are discussed, and the results are presented in tables, diagrams, graphs, maps, and photographs illustrating typical formations. Most of the earthquakes are shown to occur in the transition zone from the sub-Andes to the Cordillera Oriental under formations of about 1 km elevation at focal depths of 10-38 km. It is suggested that the sub-Andean earthquakes reflect hinterland deformation of a detached fold and thrust belt, perhaps like that which occurred in parts of the Canadian Rockies. From the total crustal shortening evident in Andean morphology and the shortening rate of the recent earthquakes it is estimated that the topography and crustal root of the Andes have been formed during the last 90-135 Myr.

Suarez, G.↗

In situ study of microstructure evolution and α → ω phase transition in annealed and pre-deformed Zr under hydrostatic loading

The detailed study of the effect of the initial microstructure on its evolution under hydrostatic compression before, during, and after the irreversible α → ω phase transformation and during pressure release in Zr using in situ x-ray diffraction is presented. Two samples were studied: one is plastically pre-deformed Zr with saturated hardness and the other is annealed. Phase transformation α → ω initiates at lower pressure for a pre-deformed sample but for a volume fraction of ω Zr, c > 0.7, a larger volume fraction is observed for the annealed sample. This implies that the proportionality between the athermal resistance to the transformation and the yield strength in the continuum phase transformation theory is invalid; an advanced version of the theory is outlined. Phenomenological plasticity theory under hydrostatic loading is outlined in terms of microstructural parameters, and plastic strain is estimated. During transformation, the first rule is suggested, i.e., the average domain size, microstrain, and dislocation density in ω Zr for c < 0.8 are functions of the volume fraction, c of ω Zr only, which are independent of the plastic strain tensor prior to transformation and pressure. The microstructure is not inherited during phase transformation. Surprisingly, for the annealed sample, the final dislocation density and the average microstrain after pressure release in the ω phase are larger than for the severely pre-deformed sample. The results suggest that an extended experimental basis is required for the predictive models for the combined pressure-induced phase transformations and microstructure evolutions.

36 MATERIALS SCIENCE↗

Measurement of crack length in width tapered beam experiments

The width tapered beam method for measuring fracture toughness has been contemplated for use in measuring the adhesion of photovoltaic materials, is promoted by it being a viable method for use directly on modules rather than to engineered test specimens. However, precise determination of crack length is often difficult as failure often involves tendril formation, cavitation voids, and other plastic deformation at the crack front in addition to obfuscation by opaque materials. Because fracture energy varies as the square of crack length, imprecise crack tip length results in large measurement uncertainties. In this work, we develop a method to estimate the crack tip length using the relationship between the applied force and the displacement of the load frame. The crack tip is determined for every data point and can be defined from a calibration curve related to the amount of energy put into the beam, or from a theoretically idealized estimation. Once a calibration curve is set up and the equations programmed, the fracture energy as a function of position on the newly formed crack interface can be accurately estimated. The test additionally provides much easier analysis and insight into the joint failure through more precise knowledge of fracture toughness at each de-bonded area.

42 ENGINEERING↗

Free vibration of arches flexible in shear.

An analysis reported by Veletsos et al. (1972) concerning the free vibrational characteristics of circular arches vibrating in their own planes is considered. The analysis was based on a theory which neglects the effects of rotatory inertia and shearing deformation. A supplementary investigation is conducted to assess the effects of the previously neglected factors and to identify the conditions under which these effects are of practical significance or may be neglected. A simple approximate procedure is developed for estimating the natural frequencies of arches, giving due consideration to the effects of the previously neglected factors.

Austin, W. J.↗