Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “soil 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.

30 records · Page 2

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗

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↗

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗

Insights into the thermo-hydraulic properties of compacted MX80 bentonite during hydration under elevated temperature

In high-level radioactive waste geological repositories, compacted bentonite undergoes coupled thermo-hydraulic processes due to heat released from a central waste canister and groundwater imbibition from the surrounding host rock. An understanding of these processes is essential for long-term simulations of radionuclide migration and canister corrosion, which requires an understanding of temperature effects on the coupled thermo-hydraulic properties governing these processes. In this study, a tank-scale radial infiltration test was used to investigate water imbibition processes in compacted bentonite under a central heater temperature of 200 °C that simulates high thermal gradients in a repository. Interpretation of this test focuses on evaluation of the liquid water wetting front during hydration and the interpretation of the temperature-dependent transient soil water retention curve (SWRC), thermal conductivity function (TCF), and hydraulic conductivity function (HCF). The SWRC during imbibition follows a temperature-dependent wetting path. Temperature effects on the hydraulic conductivity of bentonite in saturated conditions had the greatest effect on the shape of the HCF, with minimal temperature effects at higher suctions. The transient thermal conductivity data matched well with a new TCF linked with the SWRC shape. Although the bentonite layer was restrained, local deformations during hydration may have affected the shapes of the TCF and HCF.

Engineering↗

The physics of desiccation cracks 1: Ductile fracturing and dependence on relative humidity

Ductile deformation is ubiquitously found in the shrinkage of geomaterials. The existence of ductility requires elastoplastic mechanics when analyzing the structure deformation under external stress. Here, in this work, we explore the physics of ductile fracturing based on the results from desiccation experiments and triaxial tests. By using the digital image correlation (DIC) method to generate the strain maps of samples undergoing desiccation cracking under different relative humidities, we obtain results showing the previously postulated Cnoidal Wave ductile failure patterns propagating under atmospheric condition-controlled crack velocities. We then correlate these observations with rate-dependent plasticity models calibrated through triaxial tests undergoing several unloading–reloading cycles and velocity stepping. This work demonstrates the necessity to consider ductility in soil cracking, indicating that the formation of crack patterns in soil desiccation is a slow and predictable process.

58 GEOSCIENCES↗

Simulation of axial tensile well deformation during reservoir compaction in offshore unconsolidated methane hydrate-bearing formation

Sand production encountered in the 2013 offshore field gas production tests at the Nankai Trough, Japan, could be attributed to well failure during reservoir compaction. In this study, well integrity under various reservoir compaction patterns for the Nankai Trough case is examined using a well-formation finite element model. Here, the modelling details include the inclusion of a cement sheath as well as the modelling of construction processes (such as cement shrinkage). Well elongation in the overburden layer becomes significant when the reservoir subsidence is localized near the wellbore under large depressurization. Results show that the maximum plastic deviatoric strain level in the cement could reach 0.7% when the maximum reservoir subsidence reaches 0.85 m and cement shrinkage is limited. When cement shrinkage rises to 0.75%, the maximum plastic deviatoric strain increases to 2.4% as the cement accumulates additional plastic strain during shrinkage due to its deformation being constrained by the casing. In order to prevent the cement from failure, it might be effective to hold the pressure drawdown at a low level (e.g., several MPa) until the hydrate dissociation front advances to a certain radius from the well (e.g., a couple of tens of metres).

58 GEOSCIENCES↗

Computational Analysis of Coupled Geoscience Processes in Fractured and Deformable Media

Prediction of flow, transport, and deformation in fractured and porous media is critical to improving our scientific understanding of coupled thermal-hydrological-mechanical processes related to subsurface energy storage and recovery, nonproliferation, and nuclear waste storage. Especially, earth rock response to changes in pressure and stress has remained a critically challenging task. In this work, we advance computational capabilities for coupled processes in fractured and porous media using Sandia Sierra Multiphysics software through verification and validation problems such as poro-elasticity, elasto-plasticity and thermo-poroelasticity. We apply Sierra software for geologic carbon storage, fluid injection/extraction, and enhanced geothermal systems. We also significantly improve machine learning approaches through latent space and self-supervised learning. Additionally, we develop new experimental technique for evaluating dynamics of compacted soils at an intermediate scale. Overall, this project will enable us to systematically measure and control the earth system response to changes in stress and pressure due to subsurface energy activities.

58 GEOSCIENCES↗

Assessing continuum plasticity postulates with grain stress and local strain measurements in triaxially compressed sand

Critical state and continuum plasticity theories have been used in research and engineering practice in soil and rock mechanics for decades. These theories rely on postulated relationships between material stresses and strains. Some classical postulates include coaxiality between stress and strain rates, stress–dilatancy relationships, and kinematic assumptions in shear bands. Although numerical and experimental data have quantified the strains and grain kinematics in such experiments, little data quantifying grain stresses are available. Here, we report the first-known grain stress and local strain measurements in triaxial compression tests on synthetic quartz sands using synchrotron X-ray tomography and 3D X-ray diffraction. We use these data to examine the micromechanics of shear banding, with a focus on coaxiality, stress-dilatancy, and kinematics within bands. Our results indicate the following: 1) elevated deviatoric stress, strain, and stress ratios in shear bands throughout experiments; 2) coaxial principal compressive stresses and strains throughout samples; 3) significant contraction along shear bands; 4) vanishing volumetric strain but nonvanishing stress fluctuations throughout samples at all stages of deformation. Our results provide some of the first-known in situ stress and strain measurements able to aid in critically evaluating postulates employed in continuum plasticity and strain localization theories for sands.

42 ENGINEERING↗

Clean Condensing Gas Furnace

Natural gas furnaces are the most common space heating equipment in the U.S. residential and commercial building markets. However, current residential natural gas condensing furnaces generate substantial acidic condensate as well as significant emissions of sulfur oxides (SOx), nitrogen oxides (NOx), carbon monoxide (CO), hydrocarbons (HC), and methane (CH 4 ) contributing to environmental degradation of air, water, and soil. This report describes a novel solution to reduce the environmental impact of natural gas condensing furnaces based on the technology of a monolithic acidic gas reduction (AGR) catalyst for SOx trapping, NOx redox to nitrogen, and oxidation of formic acid, CO, HC, and CH 4 . The AGR technology offers the following benefits: (1) a neutral furnace condensate with a pH of ~7, allowing its safe release into the sewer system thus eliminating a second drainage system; (2) trapping and removing nearly all SOx emissions; (3) NOx emissions nearly at nearly 1-2 ng/J, more than 95% lower than new emissions standards in California; (4) the use of a low-cost heat exchanger as a condensing heat exchanger (HX) since the condensate is not acidic, avoiding the need for expensive stainless steel alloys; and (5) unburnt fuel energy recovery to boost efficiency.The AGR component and AGR-enabled furnace performance were broadly tested to determine their effects on long-term reliability and durability, as well as SOx storage and regeneration activity. The AGR regeneration does not impair the performance in achieving neutral condensate and ultra-low NOx emissions, and the AGR catalyst subjected to regeneration activities continued to function well and achieved slightly better annual fuel utilization efficiency (AFUE). The 400-hour reliability and durability test of the retrofitted condensing furnace with the AGR component shows that the furnace unit still achieves a neutral furnace condensate with a pH of ~7 and enables 0~3 ng/J of NOx emissions. However, the 400-hour operation slightly degraded the AFUE because of soot particle accumulation caused by frequent incomplete combustion owing to inappropriate condensate drainage during testing. Thus, proper condensate drainage is critical for AGR-enabled furnaces. Furthermore, neutron computed tomography was employed to survey the aged AGR component and demonstrate high-resolution 2D and 3D representations for the nondestructive diagnosis of the AGR component. The tomography showed that the AGR component did not deform or suffer broken AGR channels. A new AGR catalyst with low precious metal loading was preliminarily explored to identify a pathway of optimizing AGR material loading and maximizing acidic gas reduction at low cost. The new AGR component can reduce precious metal loading by 38% and still achieve neutral condensate and ultralow NOx emissions. The furnace with the AGR component of low Pt/Rh loading enables a maximum AFUE of 97%, which is meaningfully higher than the original furnace. Long-duration testing for the furnace enabled with the low precious metal loading AGR component will be vital in future research. Although the current work demonstrates a proof of concept for the AGR-enabled furnace, the AGR assembly needs to be optimized and integrated into the design of new OEM furnace products. Furthermore, the AGR technology can be applied not only for residential gas furnaces, but also for commercial rooftop units, gas heat pumps, gas-fired water heaters, combustion boilers, and other systems.

03 NATURAL GAS↗

Benchmark study of a new simplified DFN model for shearing of intersecting fractures and faults

It is challenging to quantitatively predict shearing of intersecting fractures/faults because of dynamic frictional contacts accompanied by possible nonlinear rock deformation. To address such challenges, a new conceptual model—the simplified DFN model—was proposed and validated by Hu et al. 46 to use major paths (MPs) to represent complicated DFNs for calculation of shearing. In this work, we conducted a benchmark study for three examples that involve different levels of complexity of intersecting fractures, and correspondingly different numbers of MPs. The codes and software that were used in the benchmark cover a range of continuum, discontinuum and hybrid numerical methods: NMM (LBNL), FLAC3D (LBNL), GBDEM (KIGAM), FRACOD (DynaFrax), and CASRock (CAS). The general consistency between DFN and MP cases as predicted by all the codes/software demonstrates that major paths can be used to simplify the geometry of DFNs in a wide range of software. Disagreement in results made by some software and potential future improvements are discussed. We show that (1) shearing of one or multiple major fractures can be reduced if there are multiple smaller intersecting fractures in that area, which is a useful basis for understanding and controlling induced seismicity and merits further analysis, and (2) the agreement achieved in the benchmark examples provide confidence that the simplified DFN model is a promising conceptual model that can be used for different types of numerical approaches and software for simplifying the analysis of the shearing of intersecting fractures and faults.

58 GEOSCIENCES↗

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗

Hydraulic Conductivity Measurements, Utqiagvik (Barrow), Alaska, 2014

Six individual ice cores were collected from the Barrow Environmental Observatory in Barrow, Alaska, in May of 2013 as part of the Next Generation Ecosystem Experiment (NGEE). Each core was drilled at a different location to varying depths. After drilling, the cores were stored in coolers packed with dry ice and flown to Lawrence Berkeley National Laboratory (LBNL) in Berkeley, CA. 3-dimensional images of the cores were constructed using medical X-ray computed tomography (CT) scanner at 120kV. Hydraulic conductivity samples were extracted from these cores at LBNL Richmond Field Station in Richmond, CA, in February 2014 by cutting 5 to 8 inch segments using a chop saw. Samples were packed individually and stored at -20C freezing temperatures to minimize any changes in structure or loss of ice content prior to analysis. Hydraulic conductivity was determined through falling head tests using a permeameter [ELE International, Model #: K-770B] (Appendix A). Samples were placed in a latex membrane via a membrane stretcher while frozen. Use of a membrane stretcher made the membranes easier to secure and minimized contact with the sample. A clear polycarbonate sleeve, fabricated with a stainless steel ring at the bottom to keep the sleeve from floating, was placed around the sample inside the permeameter to minimize deformation during analysis. The permeameter was filled with water and 1.0 PSI of air was applied for confining pressure during sample defrost. Outflow valves were left open to allow for incremental thawing and samples were left to thaw for approximately 12 hours. After approximately 12 hours of thaw, initial falling head tests were performed. When the flow was significantly too fast or too slow, the analysis was stopped and the burette size was adjusted accordingly (i.e. a larger diameter burette was used for flows that were faster than desired or a smaller diameter burette was used for flows that were slower than desired). Two to four measurements were collected on each sample and collection stopped when the applied head load exceeded 25% change from the original load. Analyses were performed between 2 to 3 times for each sample. The final hydraulic conductivity calculations were computed using methodology of Das et al., 1985.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗