SEARCH · Engineering Papers

Results for “fracture characterization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

DOE OSTI2023

The effect of inhomogeneities in coal-based mesophase pitch on filament formation and melt spinning for carbon fiber

Mesophase pitch-based carbon fibers excel as high specific strength and modulus materials, but challenges in stable, continuous melt spinning processing are a serious limiting factor. Inhomogeneities in the pitch may result in flow instabilities that can lead to filament breakage. This work aims to characterize and offer mitigation strategies for the negative effects that inhomogeneities such as isotropic binder, gels, and other inclusions impart on melt spinning of mesophase pitch. The effect of inhomogeneities on spinning stability are analyzed using capillary rheometry, borescope imaging of elongating filaments, and filtration. With no filter, immediate blockage of the spinning capillary occurred. Improvement in melt spinning duration was observed when a porous metal filter (20 μm nominal pore size) was utilized, but no significant further improvement in spinning duration was found upon decreasing the filter nominal pores size to 10, 5 and 2 μm. Moreover, flow instabilities that disrupted melt spinning stability were observed by capillary rheometry when the size of included isotropic pitch domains were larger than approximately 10 μm. Resultant carbon fiber properties and SEM imaging of fracture surfaces are reported.

01 COAL, LIGNITE, AND PEAT↗

NASA NTRS2002

Rapid Prototyping of Continuous Fiber Reinforced Ceramic Matrix Composites

For ceramics to be used as structural components in high temperature applications, their fracture toughness is improved by embedding continuous ceramic fibers. Ceramic matrix composite (CMC) materials allow increasing the overall operating temperature, raising the temperature safety margins, avoiding the need for cooling, and improving the damping capacity, while reducing the weight at the same time. They also need to be reliable and available in large quantities as well. In this paper, an innovative rapid prototyping technique to fabricate continuous fiber reinforced ceramic matrix composites is described. The process is simple, robust and will be widely applicable to a number of high temperature material systems. This technique was originally developed at the University of Delaware Center for Composite Materials (UD-CCM) for rapid fabrication of polymer matrix composites by a technique called automated tow placement or ATP. The results of mechanical properties and microstructural characterization are presented, together with examples of complex shapes and parts. It is believed that the process will be able to create complex shaped parts at an order of magnitude lower cost than current CVI and PIP processes.

Vaidyanathan, R.↗

NASA NTRS2003

Rapid Prototyping of Continuous Fiber Reinforced Ceramic Matrix Composites

Vaidyanathan, R.↗

NASA NTRS2024

Impact Response of Thermoplastic Composites – Experiments and Modeling

This study investigates the impact response of thermoplastic composite materials. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK reinforced with T700G (T700/LMPAEK ). Several High Energy Dynamic Impact (HEDI) tests were performed wherein an aluminum projectile impacted T700/LMPAEK panels at different velocities. In the impact experiment studied herein, the projectile impacted the T700/LMAPEK panel at 73.2 m/s and rebounded. The ensuing damage to the T700/LMPAEK panel was characterized via ultrasonic C-scans, visualizing damage modes including matrix cracking and delaminations. Since the projectile rebounded, the fiber breakage damage mode was not evident. The impact event was modeled using a combined continuum damage mechanics and cohesive zone analysis approach with the commercially available hydrocode LS-DYNA®. The dynamic deformation and damage of the continuum, i.e., the plies in the T700/LMPAEK laminate, was modeled using a Deformation Gradient Decomposition (DGD) material model, i.e., MAT299 available within LS-DYNA. The interlaminar regions were treated as zero-thickness cohesive zones, wherein the dynamic delamination mechanics was modeled using a traction-separation formulation. Mixed-mode delamination growth was captured using the Benzegaggh-Kenane (B-K) traction-separation interpolation law. Strain-rate sensitivity of interlaminar fracture toughness was included using approximations based on published literature. The outcomes of the modeling effort are discussed herein which demonstrate good correlations with test data regarding panel deflection, panel-projectile interactions, damage modes and damage extent.

Impact modeling↗

DOE OSTI2020

Preliminary Study of Potential Utilization of Triassic Rift Basins for Long-Term Carbon Sequestration

Triassic rift basins of the eastern United States present a potential option for long–term carbon dioxide (CO 2 ) sequestration. Exposed and buried basins are located near large point CO 2 sources. Because of their similar origins, comparable fill successions indicate a level of reproducibility that can be conveyed between basins. The thick rock successions of the exposed Culpeper and Gettysburg basins served as proxies for understanding of the basin sequences. Their study allowed recognition of five mappable assemblages of rock types, herein termed lithofacies associations. These lithofacies associations were formed by alluvial fan, braided and meandering streams and marginal and distal lake depositional processes.Stratigraphic architecture of the lithofacies associations within exposed basins suggests a vertical succession that consists of an initial coarse-grained fluvial succession that is progressively replaced upward and basinward by finer grained lacustrine deposits. Along the faulted margins of the basins alluvial fans transitioning to fluvial delta deposits reflect lateral progradation and contemporaneous lateral components of the basin filling.To test the model devised in exposed basins, facies associations were developed for the buried Taylorsville Basin. This basin preserves over 8,000 feet of Triassic rocks that are concealed beneath more than 2,000 feet of Cretaceous and Tertiary Coastal Plain sediments. Composition and stratigraphic architecture of these rocks are similar to those observed in exposed basins. Infilling of the basin was the result of vertical aggradation and lateral progradation of coarse-grained to fine-grained facies. Based upon study of outcrops in the Culpeper and Gettysburg basins, groups of recurring lithologic facies were identified as the fundamental constructs of the basin’s sediment infilling. These groups of facies, termed lithofacies associations, represent an amalgamation of lithologic components from broadly similar depositional systems. These lithologic associations were then extrapolated into the exposed and buried portions of the Taylorsville basins. This effort determined that Triassic rift basins provide several avenues for potential study of geologic sequestration of carbon. The characteristic rift basin succession in the basins presents possible reservoir targets within the marginal fluvial deposits. Furthermore, intrusive and extrusive mafic bodies present a potential source of fracture porosity that could serve as CO 2 reservoirs. Porosity values from thin section analysis of exposed basin samples are higher in alluvial fan and braided fluvial lithofacies. Porosity is highest in samples that show less compaction, due to a lack of ductile lithic fragments and/or higher stratigraphic position. Log data from the Taylorsville Basin indicates that porosity and permeability values in basin marginal fluvial strata are elevated. Near the center of the basin, porosity and permeability values are greatly reduced, primarily owing to the fine-grained character of the lacustrine deposits. Some Triassic basins also contain thick intervals of concordant extrusive and intrusive mafic igneous rocks. These igneous bodies are potential CO 2 reservoirs for several reasons. Firstly, extrusive lava flows provide potential storage in the layers of primary porosity that occur in the vesicules formed at the top of lava flows. Secondly, both lava flows and subsurface igneous sills exhibit extensive fracture porosity produced by the rapid cooling of the flows and intrusions. Thirdly, these igneous rocks are mafic in composition, and studies have shown that iron- and magnesium-rich mafic rocks provide sequestration opportunities through carbonate remineralization. Lastly, extrusive and intrusive igneous rocks are invariably preserved within fine-grained lake deposits. These lacustrine sediments can serve as a fine-grained confining layer that encases the igneous rocks both above and below.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Spark plasma sintering (SPS) densified U 3 Si 2 pellets: Microstructure control and enhanced mechanical and oxidation properties

Dense U 3 Si 2 pellets with controlled grain structure and enhanced thermal-mechanical and oxidation properties are synthesized with spark plasma sintering (SPS). Microstructure and phase composition of the SPS densified pellets are characterized systematically using SEM, EDS, and XRD. Thermal-mechanical properties and oxidation behavior of the sintered silicide fuel pellets are analyzed by laser flash, indentation, and dynamic thermogravimetric analysis. Dense U 3 Si 2 pellets are consolidated by combining high energy ball milling and rapid sintering by SPS, and the microstructure structures are controlled from micron-sized (~5.7 μm grain size) for conventional silicide to a nanocrystalline matrix with an average grain size of ~280 nm. A dominant phase of distorted U 3 Si 2 was identified with lattice expansion due to residual thermal stress upon SPS consolidation and rapid cooling processes. Both micron-sized and nano-sized pellets show exceptional thermal transport properties, consistent with monolithic silicides reported in literature. The SPS-densified pellets possess simultaneously high hardness and fracture toughness. The SPS-densified silicide pellets also demonstrate exceptional oxidation performance with extended onset oxidation temperature above 500 °C and reduced oxidation kinetics, particularly for nano-sized pellets. A strong strain effect was proposed in which compressive stress in nano-sized pellets enhances the oxidation resistance of silicide fuels, as evidenced by the degradation of oxidation performance upon strain relaxation by isothermal annealing. The correlation among the sintering process – microstructure control – physical properties and fuel behavior is established. A new concept of strain engineering is proposed further properties optimization, enabling the development of potential oxidation and corrosion-resistant silicides with extended performance, the key technological challenge of U 3 Si 2 as the leading concept of accident tolerant fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2000

Effects of Interface Coating and Nitride Enhancing Additive on Properties of Hi-Nicalon SiC Fiber Reinforced Reaction-Bonded Silicon Nitride Composites

Strong and tough Hi-Nicalon SiC fiber reinforced reaction-bonded silicon nitride matrix composites (SiC/ RBSN) have been fabricated by the fiber lay-up approach. Commercially available uncoated and PBN, PBN/Si-rich PBN, and BN/SiC coated SiC Hi-Nicalon fiber tows were used as reinforcement. The composites contained approximately 24 vol % of aligned 14 micron diameter SiC fibers in a porous RBSN matrix. Both one- and two-dimensional composites were characterized. The effects of interface coating composition, and the nitridation enhancing additive, NiO, on the room temperature physical, tensile, and interfacial shear strength properties of SiC/RBSN matrix composites were evaluated. Results indicate that for all three coated fibers, the thickness of the coatings decreased from the outer periphery to the interior of the tows, and that from 10 to 30 percent of the fibers were not covered with the interface coating. In the uncoated regions, chemical reaction between the NiO additive and the SiC fiber occurs causing degradation of tensile properties of the composites. Among the three interface coating combinations investigated, the BN/SiC coated Hi-Nicalon SiC fiber reinforced RBSN matrix composite showed the least amount of uncoated regions and reasonably uniform interface coating thickness. The matrix cracking stress in SiC/RBSN composites was predicted using a fracture mechanics based crack bridging model.

Bhatt, Ramakrishana T.↗

NASA NTRS2016

Pluto and Charon: Surface Colors and Compositions - A Hypothesis

The surface of Pluto displays an array of colors ranging from yellow to red to brown, while the surface of Charon is largely gray with a north polar zone of red color similar to regions on Pluto. Pluto's surface shows layers of intensely colored material in tilted and transported blocks, and fractured geo-graphical units. This arrangement suggests episodes of formation or deposition of that material interspersed with episodes of emplacement of ices having little or no color. The ices identified on the surfaces of these two bodies (N2, CH4, CO, C2H6, H2O on Pluto, and H2O and NH3 on Charon) are colorless, as are nearly all ices in a powdery state. The colors on Pluto probably arise from the in situ formation of a macro-molecular carbonaceous material generated by energetic processing of the ices on the surface. Laboratory experiments producing refractory tholins particularly relevant to Pluto explored the chemistry of both UV and low-energy electron bombardment of a mix of Pluto ices (N2:CH4:CO = 100:1:1). We can term this Pluto ice tholin PIT. Water ice in the crystalline state characterizes Charon's surface, and while most of Charon's surface is neutral in color, with geometric albedo approximately 0.38, the polar zone and a light cover of fainter but similar reddish color over some surface regions suggest a common origin with the colored material on Pluto. NH3 or NH3 x nH2O was identified from disk-integrated Earth-based spectra, and a few concentrated NH3 exposures have been found in the New Horizons spectral images.

Pluto↗

DOE OSTI2023

Degradation-Resistant Hypoxia Inducible Factor-2α in Murine Osteocytes Promotes a High Bone Mass Phenotype

Molecular oxygen levels vary during development and disease. Adaptations to decreased oxygen bioavailability (hypoxia) are mediated by hypoxia-inducible factor (HIF) transcription factors. HIFs are composed of an oxygen-dependent α subunit (HIF-α), of which there are two transcriptionally active isoforms (HIF-1α and HIF-2α), and a constitutively expressed β subunit (HIFβ). Under normoxic conditions, HIF-α is hydroxylated via prolyl hydroxylase domain (PHD) proteins and targeted for degradation via Von Hippel-Lindau (VHL). Under hypoxic conditions, hydroxylation via PHD is inhibited, allowing for HIF-α stabilization and induction of target transcriptional changes. Our previous studies showed that Vhl deletion in osteocytes (Dmp1-cre; Vhlf/f) resulted in HIF-α stabilization and generation of a high bone mass (HBM) phenotype. The skeletal impact of HIF-1α accumulation has been well characterized; however, the unique skeletal impacts of HIF-2α remain understudied. Because osteocytes orchestrate skeletal development and homeostasis, we investigated the role of osteocytic HIF-α isoforms in driving HBM phenotypes via osteocyte-specific loss-of-function and gain-of-function HIF-1α and HIF-2α mutations in C57BL/6 female mice. Deletion of Hif1a or Hif2a in osteocytes showed no effect on skeletal microarchitecture. Constitutively stable, degradation-resistant HIF-2α (HIF-2α cDR), but not HIF-1α cDR, generated dramatic increases in bone mass, enhanced osteoclast activity, and expansion of metaphyseal marrow stromal tissue at the expense of hematopoietic tissue. Our studies reveal a novel influence of osteocytic HIF-2α in driving HBM phenotypes that can potentially be harnessed pharmacologically to improve bone mass and reduce fracture risk.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Utah FORGE: Phase 3 Magnetotelluric (MT) Data

New high-quality tensor MT data at 122 sites, including the vertical magnetic field and utilizing ultra-remote referencing, have been acquired over the Utah FORGE project area. The results will be used to delineate the densities of faults and fractures in crystalline basement rocks, to define the heat sources, and to derive baseline 3D resistivity structure for later MT monitoring of temporal changes in resistivity structure following well stimulation in the EGS reservoir. There are three files here related to Utah FORGE magnetotelluric (MT) data acquisition and processing. The FORGE MT EDIs zip file contains the observed MT responses in industry-standard EDI format. For each site, there is an EDI response file that utilized a local independent reference for noise cancellation, and a file that utilized a distant reference cancelling noise associated with the DC transmission line of the Delta IPP passing down the west side of Milford Valley. These two site files could be merged as appropriate. The FORGE Model Cell Center file contains the model volume of the 3D Forge MT inversion for characterizing the resistivity structure in the project area. It was derived using finite element inversion methodology described in Wannamaker et al, in the attached FORGE Phase3 Geophysics paper, from the MT observation EDI files. It is ASCII format (.dat) and entries are defined at top of the file in a simple x-y-z-Rho listing in UTM coordinates. The element layers drape the topography so the Rho value layers are not purely horizontal slices. This greatly simplifies the listing.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

Characterization of Core from Sanford Underground Research Facility

Computed tomography and special core analysis data associated with the technical report series document Computed Tomography Scanning and Geophysical Measurements of the Enhanced Geothermal Systems (EGS) Collab SURF Core. Paronish, T.; Mackey, P.; Schmitt, R.; Crandall, D.; Moore, J.; Brown, S.; Roggenthen, W.; Schwering, P. C.; Dobson, P. F.; Kneafsey, T. Computed Tomography Scanning and Geophysical Measurements of the Enhanced Geothermal Systems (EGS) Collab SURF Core; DOE.NETL-2021.2866; NETL Technical Report Series; U.S. Department of Energy, National Energy Technology Laboratory: Morgantown, WV, 2022; p 76.

computed tomography↗

NASA NTRSDate not supplied

Engineering the Interface: Advanced Surface Technologies for Lunar Dust Management and Equipment Longevity

Through the Artemis program, NASA intends to develop a sustainable human foothold on the Moon, ultimately paving the way for crewed exploration of Mars. The Moon's hostile environment poses numerous obstacles, including exposure to radiation, temperature extremes, micrometeoroid threats, and particularly the persistent problem of lunar dust. Lunar dust impacts nearly every aspect of surface operations through adhesion and abrasion mechanisms, with contamination from anthropogenic activities (landing, rovers) far outweighing natural phenomena. Multiple adhesion pathways contribute to surface contamination in the lunar environment, including van der Waals forces, electrostatic forces, chemical reaction, and magnetic forces from elemental iron deposits. Sharp asperities from micrometeoroid bombardment and atmospheric absence increase interaction potential and enable mechanical interlocking. Low cohesion between dust particles exacerbates these challenges, as minimal interaction potential between dust and nearby surfaces overcomes particle cohesion, causing contamination. Lunar dust adhesion mitigation technologies can be categorized as either active, requiring external energy, or passive, relying on intrinsic material properties. Ultrasonic and electrodynamic technologies have been developed to the highest technology readiness level for active approaches. Passive strategies primarily focus on surface chemistry and topography modifications. At NASA Langley Research Center, approaches include surface migration agents to reduce surface energy, topographical modification using laser ablation patterning, and tailored surface conductivity to reduce intrinsic adhesion force. Performance has been evaluated using custom-built ultrasonic and centrifuge instruments. Plume-surface interactions from lunar landers can propel micrometer-sized particles at velocities up to 1000 m s-1.8 These particles pose risks to landers, habitats and infrastructure, leading to erosion, degradation, and reduced component lifespan. A panel recovered from Surveyor III was determined to have been severely abraded because of lunar dust displaced from the Apollo 12 lunar module that landed 160 m away. The performance of metallic surfaces has been evaluated via high velocity single particle impact using the laser-induced project impact test (LIPIT) facility at the University of Utah. Peridynamics modeling, a form of continuum mechanics that uses a nonlocal approach enabling greater simulation capabilities of crack initiation and fracture, has also been utilized to gain greater insight into material response during impact events. Lunar dust contamination challenges extend to power generation systems and moving equipment. Cables, rotation stages, and other mechanisms may experience limited range of motion and reduced lifetime due to dust infiltration. NASA Langley Research Center has evaluated traditional aerospace alloys, softgoods, wear resistant ceramics, and several polymer and polymer composite materials. Test methods have included traditional techniques like Taber abrasion testing, as well as designed test configurations developed in the DUSTE (dust, ultraviolet radiation, and space thermal environmental) chamber that reproduce mechanism functions in operational environment. Beyond laboratory experiments, several flight experiments have been conducted. Materials were exposed to the low Earth orbit environment on the Materials International Space Station Experiment (MISSE) and to the lunar surface environment through the Aegis Aerospace Regolith Adherence Characterization (RAC) payload and the Honeybee Robotics PlanetVac payload. Determining lunar dust's impact on surface exploration and habitation requires comprehensive experimental and computational capabilities combined with lessons learned from initial lunar activities. Identifying the greatest environmental challenges and developing mitigation technologies provides the clearest path toward successfully, expeditiously, and efficaciously completing NASA's mission. This presentation will discuss ongoing efforts at NASA Langley Research Center and collaborator contributions to these critical objectives.

Surface Engineering↗

NASA NTRS2001

Durability Characterization of Advanced Polymeric Composites at Cryogenic Temperatures

The next generation of reusable launch vehicles will require technology development in several key areas. Of these key areas, the development of polymeric composite cryogenic fuel tanks promises to present one of the most difficult technical challenges. It is envisioned that a polymer matrix composite (PMC) tank would be a large shell structure capable of containing cryogenic fuels and carrying a range of structural loads. The criteria that will be imposed on such a design include reduced weight, conformal geometry, and impermeability. It is this last criterion, impermeability, that will provide the focus of this paper. The essence of the impermeability criterion is that the tank remains leak free throughout its design lifetime. To address this criterion, one of the first steps is to conduct a complete durability assessment of the PMC materials. At Langley Research Center, a durability assessment of promising new polyimide-based PMCs is underway. This durability program has focused on designing a set of critical laboratory experiments that will determine fundamental material properties under combined thermal-mechanical loading at cryogenic temperatures. The test program provides measurements of lamina and laminate properties, including strength, stiffness, and fracture toughness. The performance of the PMC materials is monitored as a function of exposure conditions and aging time. Residual properties after exposure are measured at cryogenic temperatures and provide quantitative values of residual strength and stiffness. Primary degradation mechanisms and the associated damage modes are measured with both destructive and nondestructive techniques. In addition to mechanical properties, a range of physical properties, such as weight, glass transition, and crack density, are measured and correlated with the test conditions. This paper will report on the progress of this research program and present critical results and illustrative examples of current findings.

Gates, T. S.↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗

DOE OSTI2023

Computed Tomography Scanning and Petrophysical Measurements of the Lively Grove #1 Well Core

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core from the Lively Grove #1 Well (API 121892494700), drilled near Marissa, Washington County, Illinois. Core from the well was obtained as part of the Illinois Storage Corridor’s Carbon Storage Assurance Facility Enterprise (CarbonSAFE) project (DE FE0031892). The primary impetus of this work was to capture a detailed digital representation of the core from the Lively Grove #1 Well. The collaboration between the U.S. Department of Energy’s (DOE) NETL and the Illinois State Geological Survey (ISGS) enables other research entities to access information about this potential carbon storage location and Cambro-Ordovician Storage Complex formations of the Illinois Basin. The resultant datasets are presented in this report and can be accessed from NETL’s Energy Data eXchange (EDX) online system. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. Fractures, discontinuities, and millimeter-scale features were readily detectable with the medical CT scanner-acquired images. Imaging with the NETL medical CT scanner was performed on entire cores. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), gamma density, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for potential future analysis. Higher-resolution industrial and micro-CT images were acquired from selected zones along the depth of the core to visualize the structure in higher detail. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; with the resulting macro- and micro-descriptions relevant to many subsurface energy related examinations traditionally performed at NETL.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2020

Asteroid Cooling Rates Indicated by K-Feldspar Exsolution Textures in H4 Ordinary Chondrites

Undisturbed thermal metamorphism in ordinary chondrite (OC) asteroids, produced through the radioactive decay of 26Al, is expected to result in an onion-shell-like structure. In such a structure, the inner layers of the asteroid experience more extensive thermal metamorphism, as represented by higher petrologic type, than the exterior layers. Furthermore, cooling rates are expected to be slower for OCs of high petrologic type than those of low petrologic type. However, cooling rates determined using metallographic methods and pyroxene diffusion are inconsistent with onion-shell-style cooling and have resulted in new models. These models argue for the disruption of the asteroid after peak metamorphism followed by reaccretion into a rubble pile. Improved constraints on cooling rates would provide a better understanding of the timing and scale of disruptive events. Feldspar microtextures are another tool that can be used to determine asteroid cooling rates. In OC chondrules, plagioclase is present as either a primary phase, or a secondary phase forming from the crystallization of mesostasis glass through petrologic type 4, followed by chemical and textural equilibration. Potas-sium feldspar is observed in petrologic types 3.6-6, as either patches or lamellae exsolved from albite in a perthite texture, often near pores or fractures. Exsolution occurs most commonly, and most extensively, in petrologic type 4. Because the feldspar exsolution wavelength is related to the rate at which grains cooled from the solvus temperature, determined from the mineral’s bulk composition, the chondrite cooling rate can be measured from regions of exsolution. We have previously reported the perthite exsolution cooling rate of Avanhandava, an H4 chondrite, to be 1 °C per 1-4 months over a temperature interval of 765-670 °C. A peristerite exsolution texture was also present in the Na-rich lamellae for which we estimated a cooling rate of 1 °C in 103-104 years from 570-540 °C. Overall, the cooling rates determined from Avanhandava are consistent with pyroxene diffusion (fast cooling at high temperatures) and metallographic rates (slow cooling at low temperatures), hence with the rubble pile model of disruption and reaccretion. Here, we characterize feldspar microtextures in four additional H4 chondrites to test the consistency of feldspar cooling rates across a range of samples. We show that all H4s are similar and support rubble pile models.

Lewis, J. A.↗