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 613 records · Page 34

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

Here, this study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 M w 4.49) and sixth (DPRK6, 7 September 2017 M w 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (⁠V P = 2.07 km/s, V S = 1.2 km/s⁠), a 250 m stratified volcanic deposit layer (⁠V P = 1.73 km/s, V S = 1.0 km/s⁠), a 700 m weathered granodiorite layer (⁠V P = 2.5 km/s, V S = 1.3 km/s⁠⁠), and a granodiorite half-space (⁠V P = 5.35 km/s, V S = 3.09 km/s⁠⁠). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

58 GEOSCIENCES↗

Learning Canonical Embeddings for Unsupervised Shape Correspondence With Locally Linear Transformations

We present a new approach to unsupervised shape correspondence learning between pairs of point clouds. We make the first attempt to adapt the classical locally linear embedding algorithm (LLE)-originally designed for nonlinear dimensionality reduction-for shape correspondence. The key idea is to find dense correspondences between shapes by first obtaining high-dimensional neighborhood-preserving embeddings of low-dimensional point clouds and subsequently aligning the source and target embeddings using locally linear transformations. We demonstrate that learning the embedding using a new LLE-inspired point cloud reconstruction objective results in accurate shape correspondences. More specifically, the approach comprises an end-to-end learnable framework of extracting high-dimensional neighborhood-preserving embeddings, estimating locally linear transformations in the embedding space, and reconstructing shapes via divergence measure-based alignment of probability density functions built over reconstructed and target shapes. Our approach enforces embeddings of shapes in correspondence to lie in the same universal/canonical embedding space, which eventually helps regularize the learning process and leads to a simple nearest neighbors approach between shape embeddings for finding reliable correspondences. Comprehensive experiments show that the new method makes noticeable improvements over state-of-the-art approaches on standard shape correspondence benchmark datasets covering both human and nonhuman shapes.

deformation↗

Wrinkle ridges on Mercury and the Moon within and outside of mascons

Found on all terrestrial planets, wrinkle ridges are anticlines formed by thrust faulting and folding resulting from crustal shortening. The MErcury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER)spacecraft's orbital phase returned high resolution images and topographic data of the previously unimaged northern high latitudes of Mercury where there are large expanses of smooth plains deformed by wrinkle ridges. Concurrently, the Lunar Reconnaissance Orbiter (LRO) is obtaining high resolution images and topographic data covering lunar mare wrinkle ridges. These data allow quantitative comparison of the scale of wrinkle ridges in smooth plains volcanic units on Mercury with mare wrinkle ridges. We evaluate the topographic relief of 300 wrinkle ridges within and outside of mascon basins on the Moon and Mercury. Measured wrinkle ridges range from ~112 to 776 min relief with a mean of ~350m (median=~340m, n=150) on Mercury and from ~47 to 678 min relief with a mean of ~198 m (median=~168m, n=150) on the Moon. Wrinkle ridges on Mercury thus are approximately twice as large in mean relief compared to their counterparts on the Moon. The larger scale of Mercury's wrinkle ridges suggests that their formation can be attributed, in part, to global contraction. As global contraction on the Moon is estimated to be an order of magnitude smaller than on Mercury, the smaller scale of lunar wrinkle ridges suggests they most likely form primarily by load induced subsidence of the mare basalt. Wrinkle ridges located in lunar mascon basins and in the Caloris mascon on Mercury are not statistically significantly different in relief than ridges in non-mascon regions, suggesting comparable levels of contractional strain. The fact that mascon basins do not host wrinkle ridges with greater structural relief relative to non-mascon units may indicate the critical role lithospheric thickness plays in controlling subsidence and contraction of thick volcanic sequences on the Moon and Mercury.

Lisa S. Schleicher↗

Data-driven surrogates for high dimensional models using Gaussian process regression on the Grassmann manifold

This paper introduces a surrogate modeling scheme based on Grassmannian manifold learning to be used for cost-efficient predictions of high-dimensional stochastic systems. The method exploits subspace-structured features of each solution by projecting it onto a Grassmann manifold. This point-wise linear dimensionality reduction harnesses the structural information to assess the similarity between solutions at different points in the input parameter space. The method utilizes a solution clustering approach in order to identify regions of the parameter space over which solutions are sufficiently similarly such that they can be interpolated on the Grassmannian. In this clustering, the reduced-order solutions are partitioned into disjoint clusters on the Grassmann manifold using the eigen-structure of properly defined Grassmannian kernels and, the Karcher mean of each cluster is estimated. Then, the points in each cluster are projected onto the tangent space with origin at the corresponding Karcher mean using the exponential mapping. For each cluster, a Gaussian process regression model is trained that maps the input parameters of the system to the reduced solution points of the corresponding cluster projected onto the tangent space. Using this Gaussian process model, the full-field solution can be efficiently predicted at any new point in the parameter space. In certain cases, the solution clusters will span disjoint regions of the parameter space. In such cases, for each of the solution clusters we utilize a second, density-based spatial clustering to group their corresponding input parameter points in the Euclidean space. The proposed method is applied to two numerical examples. Here, the first is a nonlinear stochastic ordinary differential equation with uncertain initial conditions where the surrogate is used to predict the time history solution. The second involves modeling of plastic deformation in a model amorphous solid using the Shear Transformation Zone theory of plasticity, where the proposed surrogate is used to predict the full strain field of a material specimen under large shear strains.

42 ENGINEERING↗

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗

Mechanical Loads on Spent Nuclear Fuel in the General 30 cm Package Drop Scenario

The U.S. Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) research program is guided by the high-level goal of closing prioritized knowledge gaps related to spent nuclear fuel (SNF) storage and transportation, which are summarized by Saltzstein et al. (2020). One of the high-priority knowledge gaps is the identification and quantification of mechanical loads that are expected to affect SNF during normal conditions of transportation and storage to inform the range of physical SNF test programs. This report uses modeling and analysis methods to estimate the mechanical loads on spent nuclear fuel (SNF) in the general 30 cm package drop scenario. The drop scenario assumes impact limiters are in place in the transportation configuration and the impact surface is perfectly rigid. The goal of this analysis is to consider the universe of potential mechanical loading conditions that can happen to SNF and present the results in a manner that is useful for materials testing, decision making, and regulatory rule making purposes. This study uses validated finite element models and methods to perform a broad parametric study of key variables that can affect the mechanical loads on SNF during a hypothetical package free drop scenario. Physical drop test data from a cask and fuel assembly drop test campaign is the basis for model validation. Additionally, the results of the parametric study are used to inform a damage model, which uses multiple nonlinear regression to estimate the relationships between input variables and output response. The parametric finite element analyses consider thousands of input variable combinations, while the damage model estimates millions of combinations. The breadth of this study provides confidence that the potential range of mechanical loads that SNF might experience during the general 30 cm package drop scenario are characterized well enough to consider this knowledge gap closed. While this report documents the overall peak values calculated in this study, the 95 th percentile values, the histograms, and the observed trends are equally important. This study covered a large range of SNF temperatures, room temperature to 300°C, and burnups, 10 GWd/MTU to 62 GWd/MTU. Each temperature and burnup combination has a different cladding yield strain, so it is more meaningful to summarize the calculated cladding strain response as its factor of safety, which is defined relative to the yield strain. The factor of safety is calculated as the yield strain divided by the peak cladding strain. A factor of safety greater than unity indicates that the cladding remains below yield, whereas a value less than unity is indicative of plastic deformation. In all cases of this study a safety factor over 1.0 was calculated, although in the most limiting case at 300°C the safety factor was only 1.01, which suggests that yielding could occur when additional loads like rod internal pressure are included. When the temperature is restricted to 200°C the limiting safety factor increases to 1.28, which has significant margin to accommodate internal pressure and potential local cladding defects that could cause a local stress concentration. An important trend in the calculated fuel rod mechanical loads is that the 2 nd highest loaded fuel rod in an assembly tends to be significantly lower than the highest loaded rod. The implication is that even if one rod in an assembly experiences a failure the loads would have to be significantly higher to cause two or more rods to fail.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Constraining ΛCDM cosmological parameters with Einstein Telescope mock data

ABSTRACT We investigate the capability of Einstein Telescope to constrain the cosmological parameters of the non-flat ΛCDM cosmological model. Two types of mock data sets are considered depending on whether or not a short gamma-ray burst is detected, and associated with the gravitational wave emitted by binary neutron stars merger, using the THESEUS satellite. Depending on the mock data set, two statistical estimators are applied: one assumes that the redshift is known, while the other marginalizes over it assuming a specific redshift prior distribution. We demonstrate that (i) using mock catalogues collecting gravitational wave signals emitted by binary neutron stars systems to which a short gamma-ray burst has been associated, Einstein Telescope may achieve an accuracy on the cosmological parameters of $\sigma _{H_0}\approx 0.40$ km s−1 Mpc−1, $\sigma _{\Omega _{k,0}}\approx 0.09$, and $\sigma _{\Omega _{\Lambda ,0}}\approx 0.07$; while (ii) using mock catalogues collecting all gravitational wave signals emitted by binary neutron stars systems for which an electromagnetic counterpart has not been detected, Einstein Telescope may achieve an accuracy on the cosmological parameters of $\sigma _{H_0}\approx 0.04$ km s−1 Mpc−1, $\sigma _{\Omega _{k,0}}\approx 0.01$, and $\sigma _{\Omega _{\Lambda ,0}}\approx 0.01$, once the redshift probability distribution of GW events is known from from population synthesis simulations and/or the measure of the tidal deformability parameter. These results show an improvement of a factor 2–75 with respect to earlier results using complementary data sets.

Califano, Matteo↗

Stratigraphy of the Caloris Basin, Mercury: Implications for Volcanic History and Basin Impact Melt

Caloris basin, Mercury's youngest large impact basin, is filled by volcanic plains that are spectrally distinct from surrounding material. Post-plains impact craters of a variety of sizes populate the basin interior, and the spectra of the material they have excavated enable the thickness of the volcanic fill to be estimated and reveal the nature of the subsurface. The thickness of the interior volcanic plains is consistently at least 2.5 km, reaching 3.5 km in places, with thinner fill toward the edge of the basin. No systematic variations in fill thickness are observed with long-wavelength topography or azimuth. The lack of correlation between plains thickness and variations in elevation at large horizontal scales within the basin indicates that plains emplacement must have predated most, if not all, of the changes in long-wavelength topography that affected the basin. There are no embayed or unambiguously buried (ghost) craters with diameters greater than 10 km in the Caloris interior plains. The absence of such ghost craters indicates that one or more of the following scenarios must hold: the plains are sufficiently thick to have buried all evidence of craters that formed between the Caloris impact event and the emplacement of the plains; the plains were emplaced soon after basin formation; or the complex tectonic deformation of the basin interior has disguised wrinkle-ridge rings localized by buried craters. That low-reflectance material (LRM) was exposed by every impact that penetrated through the surface volcanic plains provides a means to explore near-surface stratigraphy. If all occurrences of LRM are derived from a single layer, the subsurface LRM deposit is at least 7.5-8.5 km thick and its top likely once made up the Caloris basin floor. The Caloris-forming impact would have generated a layer of impact melt 3-15 km thick; such a layer could account for the entire thickness of LRM. This material would have been derived from a combination of lower crust and upper mantle.

Cratering↗

Determination of single crystal thermal expansion in Uranium-6wt%Niobium shape memory alloy using in-situ diffraction and modeling of textured polycrystalline samples

In-situ neutron diffraction and ex-situ dilatometry experiments were performed on a chemically banded, quenched uranium-6wt% niobium shape memory alloy to study the impact of deformation-(detwinning-) induced texture on the evolution of thermally induced strains (and associated stresses). Thermal heating and cooling cycles between 5K and 473K were performed in-situ with neutron diffraction, and the lattice strain evolution is reported for the observed monoclinic α" phase. Comparisons between the measured diffraction strains and macroscopic dilatometry experiments reveal relationships between micro- and macro-level thermal expansion. Softening of the texture during heating suggests that twin boundary motion can accommodate the large internal thermal strains which are approximately 2% greater than that typically observed during the 175K heating interval used to age more randomly oriented polycrystalline material. Assuming weak constraint of neighboring grains in randomly textured polycrystals, the single crystal thermal expansion tensor is extracted from measurements of lattice strains over the range 5K to 473K. Predictions of polycrystalline thermal expansion, based upon this single crystal thermal expansion tensor, are shown to compare favorably with bulk thermal expansion observations of the as quenched microstructure. However, such a lower-bound estimate is insufficient to explain all aspects of the behavior of the textured material, where the matrix is not isotropic. In conclusion, it is hypothesized that relaxation processes which occur within the quenched microstructure during heating are responsible for the distinct thermal expansion behavior observed during the initial heating cycle compared to cooling and subsequent cycling.

36 MATERIALS SCIENCE↗

Microwave Resonator Measurements of Atmospheric Absorption Coefficients: A Preliminary Design Study

A preliminary design study examined the feasibility of using microwave resonator measurements to improve the accuracy of atmospheric absorption coefficients and refractivity between 18 and 35 GHz. Increased accuracies would improve the capability of water vapor radiometers to correct for radio signal delays caused by Earth's atmosphere. Calibration of delays incurred by radio signals traversing the atmosphere has applications to both deep space tracking and planetary radio science experiments. Currently, the Cassini gravity wave search requires 0.8-1.0% absorption coefficient accuracy. This study examined current atmospheric absorption models and estimated that current model accuracy ranges from 5% to 7%. The refractivity of water vapor is known to 1% accuracy, while the refractivity of many dry gases (oxygen, nitrogen, etc.) are known to better than 0.1%. Improvements to the current generation of models will require that both the functional form and absolute absorption of the water vapor spectrum be calibrated and validated. Several laboratory techniques for measuring atmospheric absorption and refractivity were investigated, including absorption cells, single and multimode rectangular cavity resonators, and Fabry-Perot resonators. Semi-confocal Fabry-Perot resonators were shown to provide the most cost-effective and accurate method of measuring atmospheric gas refractivity. The need for accurate environmental measurement and control was also addressed. A preliminary design for the environmental control and measurement system was developed to aid in identifying significant design issues. The analysis indicated that overall measurement accuracy will be limited by measurement errors and imprecise control of the gas sample's thermodynamic state, thermal expansion and vibration- induced deformation of the resonator structure, and electronic measurement error. The central problem is to identify systematic errors because random errors can be reduced by averaging. Calibrating the resonator measurements by checking the refractivity of dry gases which are known to better than 0.1% provides a method of controlling the systematic errors to 0.1%. The primary source of error in absorptivity and refractivity measurements is thus the ability to measure the concentration of water vapor in the resonator path. Over the whole thermodynamic range of interest the accuracy of water vapor measurement is 1.5%. However, over the range responsible for most of the radio delay (i.e. conditions in the bottom two kilometers of the atmosphere) the accuracy of water vapor measurements ranges from 0.5% to 1.0%. Therefore the precision of the resonator measurements could be held to 0.3% and the overall absolute accuracy of resonator-based absorption and refractivity measurements will range from 0.6% to 1.

Walter, Steven J.↗

Using Digital Image Correlation for Material Testing at Cryogenic Temperatures

Recent efforts were made to parametrize an aluminum alloy at cryogenic temperatures to calibrate finite element models simulating material deformation. Digital image correlation (DIC) was used to measure material properties using front/back strains of specimens during tensile tests. Data acquired from DIC was compared to physical extensometer measurements. Two camera pairs (front and back) were mounted via aluminum extrusion to the exterior of an environmental chamber of a hydraulic test stand. Images were taken through single-pane optical glass windows of the environment chamber. Positive pressure bags encased the cameras, attaching to the perimeter of the glass to keep out moisture and prevent windows from frosting. Figure 1 shows the environmental setup and the camera setup during cryogenic testing. Some of the challenges during testing at cryogenic temperatures include the gaseous flow of liquid nitrogen in the optical path of the mounted cameras, frost buildup on the windows as seen in Figure 1b, and frost on the test coupon. The temperature gradient present during the introduction of liquid nitrogen into the chamber can increase the projection error of images captured of the specimen. Figure 2 shows the estimated projection error from DIC as a function of temperature during the cool down process. This presentation addresses some of the concerns and observations when using DIC for testing at cryogenic temperatures using methods of isolated camera system and techniques for maintaining the glass panel clear from frost.

Joseph Cochran↗

Robust deep learning framework for constitutive relations modeling

Modeling the full-range deformation behaviors of materials under complex loading and materials conditions is a significant challenge for constitutive relations (CRs) modeling. Here, we propose a general encoder-decoder deep learning framework that can model high-dimensional stress-strain data and complex loading histories with robustness and universal capability. The framework employs an encoder to project high-dimensional input information (e.g., loading history, loading conditions, and materials information) to a lower-dimensional hidden space and a decoder to map the hidden representation to the stress of interest. We evaluated various encoder architectures, including gated recurrent unit (GRU), GRU with attention, temporal convolutional network (TCN), and the Transformer encoder, on two complex stress-strain datasets that were designed to include a wide range of complex loading histories and loading conditions. All architectures achieved excellent test results with an root-mean-square error (RMSE) below 1 MPa. Additionally, we analyzed the capability of the different architectures to make predictions on out-of-domain applications, with an uncertainty estimation based on deep ensembles. The proposed approach provides a robust alternative to empirical/semi-empirical models for CRs modeling, offering the potential for more accurate and efficient materials design and optimization.

36 MATERIALS SCIENCE↗

NASA Tech Briefs, January 2011

The topics include: 1) Distributed Aerodynamic Sensing and Processing Toolbox; 2) Collaborative Supervised Learning for Sensor Networks; 3) Hazard Detection Software for Lunar Landing; 4) Onboard Nonlinear Engine Sensor and Component Fault Diagnosis and Isolation Scheme; 5) Network-Capable Application Process and Wireless Intelligent Sensors for ISHM; 6) Interface Supports Multiple Broadcast Transceivers for Flight Applications; 7) FPGA Sequencer for Radar Altimeter Applications; 8) Miniature Sapphire Acoustic Resonator - MSAR; 9) Process-Hardened, Multi-Analyte Sensor for Characterizing Rocket Plume Constituents; 10) SAD5 Stereo Correlation Line-Striping in an FPGA; 11) Hybrid Composite Cryogenic Tank Structure; 12) Nanoscale Deformable Optics; 13) Reliability-Based Design Optimization of a Composite Airframe Component; 14) Zinc Oxide Nanowire Interphase for Enhanced Lightweight Polymer Fiber Composites; 15) Plasma Igniter for Reliable Ignition of Combustion in Rocket Engines; 16) Wire Test Grip Fixture; 17) A Sub-Hertz, Low-Frequency Vibration Isolation Platform; 18) Carbon Nanofibers Synthesized on Selective Substrates for Nonvolatile Memory and 3D Electronics; 19) Nanoparticle/Polymer Nanocomposite Bond Coat or Coating; 20) High-Resolution Wind Measurements for Offshore Wind Energy Development; 21) Spring Tire; 22) Marsviewer 2008; 23) Mission Services Evolution Center Message Bus; 24) Major Constituents Analysis for the Vehicle Cabin Atmosphere Monitor; 25) Astronaut Health Participant Summary Application; 26) Adaption of the AMDIS Method to Flight Status on the VCAM Instrument; 27) Natural Language Interface for Safety Certification of Safety-Critical Software; 28) Cryogenic Caging for Science Instrumentation; 29) Wide-Range Neutron Detector for Space Nuclear Applications; 30) In Situ Guided Wave Structural Health Monitoring System; 31) Multiplexed Energy Coupler for Rotating Equipment; 32) Attitude Estimation in Fractionated Spacecraft Cluster Systems; 33) Full Piezoelectric Multilayer-Stacked Hybrid Actuation/Transduction Systems; 34) Active Flow Effectors for Noise and Separation Control; 35) Method and System for Temporal Filtering in Video Compression Systems; 36) Apparatus for Measuring Total Emissivity of Small, Low-Emissivity Samples; 37) Multiple-Zone Diffractive Optic Element for Laser Ranging Applications; 38) Simplified Architecture for Precise Aiming of a Deep-Space Communication Laser Transceiver; 39) Two-Photon-Absorption Scheme for Optical Beam Tracking; 40) High-Sensitivity, Broad-Range Vacuum Gauge Using Nanotubes for Micromachined Cavities; 41) Wide-Field Optic for Autonomous Acquisition of Laser Link; 42) Extracting Zero-Gravity Surface Figure of a Mirror; 43) Modeling Electromagnetic Scattering From Complex Inhomogeneous Objects; 44) Visual Object Recognition and Tracking of Tools; 45) Method for Implementing Optical Phase Adjustment; 46) Visual SLAM Using Variance Grid Maps; 47) Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration; 48) Efficient Kriging Algorithms; 49) Predicting Spacecraft Trajectories by the WeavEncke Method; 50) An Augmentation of G-Guidance Algorithms; 51) Comparison of Aircraft Icing Growth Assessment Software; 52) Silicon-Germanium Voltage-Controlled Oscillator at 105 GHz; 53) Estimation of Coriolis Force and Torque Acting on Ares-1; 54) Null Lens Assembly for X-Ray Mirror Segments; and 55) High-Precision Pulse Generator.

Source record↗

Hybrid Solar System (Final Scientific/Technical Report)

GTI Energy (GTI) teamed with the University of California at Merced (UCM) to scaleup the hybrid solar system (HSS) technology for demonstrating its performance at the US Gypsum (USG) plant in Plaster City, California. The technology integrates two-stage concentrating solar collector with matching particle thermal transport and storage (TSS) system to deliver cost-effective, and on-demand distributed high temperature industrial process heat up to 600°C with solar thermal, in this case to a gypsum kettle, to reduce its fuel use and carbon footprint. Current solar technologies, which reach these temperatures, are not distributable (towers) or cost-effective (dish). The research team developed a conceptual system design for host site retrofit, including preliminary heat balance, process flow diagram, particle to process heat exchanger and equipment placements at the site. Subsequently, parallel efforts were carried out at UCM to design, build and test a 12 m long commercial scale prototype concentrating thermal-only collector system and at GTI to design, build and test a matching 650°C capable particle TTS system. The nominal 50 kWth collector consists of a parabolic trough and three 4 m long two-stage receivers in series. Prior to on-sun testing, a 4 m long receiver was fabricated and successfully tested at 650 °C in a laboratory setting for 100 hrs of continuous operation showing less than 15% radiation loss. A 7 m wide x 17 m long parabolic trough was then installed at UCM for on-sun testing of the 12 m long receiver, and concurrently several 4 m long receivers were built. The optics of the parabolic trough were calibrated, and on-sun test were carried out on 12 m long receivers. During tests, the intense solar radiation (53x) caused the absorber tubes in the receivers to bend, reducing the overall optical efficiency. To address the bending issue, a self-consistent algorithm that includes ray tracing, thermal and deformation models was developed to perform thermal stress analysis on absorbers for parabolic solar collectors. Results obtained with this algorithm showed a dramatic rise in deformation as absorber tube length increases. A combined efficiency parameter that includes the occluded area for the mounts was developed to obtain an optimized tube length obtained. Based on the results, a length of 2.7 m for the absorber + 0.2 m for the coupler was chosen to minimize any bending and optimize optical efficiency while maintaining ease of mounting. The associated particle TTS system was designed, built and successfully tested at GTI. It includes storage, receiving and lock hoppers and piping that simulates the transfer of captured solar energy to an actual industrial furnace. Tests over 77 charge-discharge cycles demonstrated <2% particle degradation, with no problematic particle accumulations and no flow interruptions. The piping pressure drop was about 5 psi. The team also worked with Stanley Consultants (Stanley) to prepare conceptual and preliminary engineering packages to facilitate follow-on development and commercialization efforts. These include process and instrumentation diagram’s (P&ID’s), general arrangements, electrical one-line, project definitions document, equipment data sheets, schedule, and construction cost estimate for 2 MWth system. Updated HSS technology commercialization and customer engagement plans and detailed costs and evaluated market trade-offs and manufacturing.

03 NATURAL GAS↗

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

Observation of pressure ridges in SAR images of sea ice: Scattering theory and comparison with observations

Ridges and keels (hummocks and bummocks) in sea ice flows are important in sea ice research for both scientific and practical reasons. Sea ice movement and deformation is driven by internal and external stresses on the ice. Ridges and keels play important roles in both cases because they determine the external wind and current stresses via drag coefficients. For example, the drag coefficient over sea ice can vary by a factor of several depending on the fluid mechanical roughness length of the surface. This roughness length is thought to be strongly dependent on the ridge structures present. Thus, variations in ridge and keel structure can cause gradients in external stresses which must be balanced by internal stresses and possibly fracture of the ice. Ridging in sea ice is also a sign of fracture. In a practical sense, large ridges form the biggest impediment to surface travel over the ice or penetration through sea ice by ice-strengthened ships. Ridges also play an important role in the damage caused by sea ice to off-shore structures. Hence, observation and measurement of sea ice ridges is an important component of sea ice remote sensing. The research reported here builds on previous work, estimating the characteristics of ridges and leads in sea ice from SAR images. Our objective is to develop methods for quantitative measurement of sea ice ridges from SAR images. To make further progress, in particular, to estimate ridge height, a scattering model for ridges is needed. Our research approach for a ridge scattering model begins with a survey of the geometrical properties of ridges and a comparison with the characteristics of the surrounding ice. For this purpose we have used airborne optical laser (AOL) data collected during the 1987 Greenland Sea Experiment. These data were used to generate a spatial wavenumber spectrum for height variance for a typical ridge - the typical ridge is the average over 10 large ridges. Our first-order model radar scattering includes both the quasi-specular and Bragg resonant scatter mechanisms. This model is extended to include contributions from volume scatter and scatter from discrete objects. Geometrical characteristics from the AOL survey and model calculations imply that for radar wavelengths and observation geometries that are dominated by the quasi-specular scattering mechanism radar backscatter from a ridge is a measure of peak ridge height. We present scattering model results and compare them with ridges observed during the LEADEX experiment of March-April 1992 when both X, C, and L-band aircraft SAR and the ERS-1 satellite SAR observed a region in the Beaufort Sea near 86 deg N, 10 deg W. Data were also collected documenting ridge characteristics on the surface. The surface data are used to generate a SAR signature via the scattering model described above. The predicted SAR signatures compare well with the SAR observations.

Vesecky, J. F.↗

Oregon Roadway CSZ Liquefaction and Landslide Impact Screening Analysis

The roadway transportation routes and the surrounding environments in Oregon will be subject to ground acceleration and potential associated deformation due to the propagation of waves of energy during a Cascadia Subduction Zone (CSZ) earthquake. Linear features such as highways may be subject to partial or complete disruption of their function resulting from ground deformations at and along the roadway. A potential source of these disruptions associated with a CSZ event are landslides and liquefaction in the soils along and beneath the highways and local roadways. The construction of the Oregon roadway system and urban development have occurred largely along the banks of the rivers—those along the Pacific coast, and along the valleys throughout the Oregon Coast Range and Cascade Mountains. These locations also represent areas where landslides and liquefaction have occurred historically, or may potentially occur in the future, presenting a significant vulnerability for the roadway system and dependent post-CSZ earthquake emergency response and recovery effort. Many of the Oregon areas with landslide- and liquefaction-prone soils have the potential to impact roadway transportation routes after a major earthquake causing delays in establishing routes to transport emergency response personnel and resources to impacted areas and communities. This report presents methods used to project the impacts that CSZ earthquake-induced landslides and liquefaction may have on roadways across Oregon. First, this report presents an approach for estimating the risk of landslides to Oregon roadways, and determining approximate reopening times. Second, this report presents a similar approach for estimating the risk of liquefaction-induced ground displacements to roadway pavements, also resulting in approximate reopening times for roadway segments. The outcomes of these two ground failure analyses serve as inputs to a broader analysis of CSZ earthquake impacts to Oregon’s surface transportation system and its ability to support post-disaster response and recovery activities, as part of the Oregon Transportation Systems Regional Resiliency Assessment Program project (CISA 2021).

58 GEOSCIENCES↗