Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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

At least 361 records · Page 20

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗

Computational Micromechanics Model for the Analysis of Fiber Kinking in Unidirectional FRP

A computational micromechanics (CMM) model is developed to analyze fiber kinking, which is a failure mechanism that takes place in fiber-reinforced composites when they are loaded under longitudinal compression. The CMM model consists of a single AS4 carbon fiber with an initial misalignment embedded in an 8552 polymer matrix. The deformation of the model is governed by periodic boundary conditions (PBC). The relatively simple CMM model enables the evaluation of the role played by initial misalignment of the fiber, shear yielding of the matrix and fiber-matrix debonding. A novel microscale experimental technique devoted to the characterization of the longitudinal compressive strength of the fibers, X(sup f)(sub c) , was developed. By exercising the model and comparing it with several existing models in the literature, the nonlinear shear response of the composite lamina was shown to play a fundamental role not only in the prediction of the compressive strength X(sub c), but also during the post-peak regime in terms of residual stress σ(sub r) and fiber rotation φ. Finally, the influence of the fiber-matrix interface damage (not considered in most other fiber kinking models) on the fiber kinking phenomenon was assessed through a parametric study.

Polymer-matrix composites (PMC)↗

The "Puck" Energetic Charged Particle Detector: Design, Heritage, and Advancements

Energetic charged particle detectors characterize a portion of the plasma distribution function that plays critical roles in some physical processes, from carrying the currents in planetary ring currents to weathering the surfaces of planetary objects. For several low-resource missions in the past, the need was recognized for a low-resource but highly capable, mass-species-discriminating energetic particle sensor that could also obtain angular distributions without motors or mechanical articulation. This need led to the development of a compact Energetic Particle Detector (EPD), known as the "Puck" EPD (short for hockey puck), that is capable of determining the flux, angular distribution, and composition of incident ions between an energy range of approximately 10 kiloelectronvolts to several megaelectronvolts. This sensor makes simultaneous angular measurements of electron fluxes from the tens of kiloelectronvolts to about 1 megaelectronvolt. The same measurements can be extended down to approximately 1 kiloelectronvolt per nucleon,with some composition ambiguity. These sensors have a proven flight heritage record that includes missions such as MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) and New Horizons, with multiple sensors on each of Juno, Van Allen Probes, and Magnetospheric Multiscale. In this review paper we discuss the Puck EPD design, its heritage, unexpected results from these past missions and future advancements. We also discuss high-voltage anomalies that are thought to be associated with the use of curved foils, which is a new foil manufacturing processes utilized on recent Puck EPD designs. Finally, we discuss the important role Puck EPDs can potentially play in upcoming missions.

Clark, G.↗

Effects of surface vacancy-adatom pair formation on PFC tungsten surface morphological response

Abstract We report a simulation study of the effect of He-irradiation-induced surface vacancy-adatom pair formation on the surface morphological evolution of plasma-facing component (PFC) tungsten and examine a number of factors that impact such evolution. Our analysis is based on self-consistent dynamical simulations according to an atomistically-informed, continuum-scale surface evolution model that has been developed following a hierarchical multiscale modeling strategy and can access the spatiotemporal scales of relevance to fuzz formation. The model accounts for the flux of surface adatoms generated as a result of the surface vacancy-adatom pair formation effect upon He implantation, which contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. We have carried out atomic-scale computations of optimal diffusion pathways along and across island step edges on the W(110) surface and calculated Ehrlich–Schwoebel (ES) barriers in adatom diffusion along and across such step edges. This aspect of surface adatom diffusion contributes to anisotropic surface atomic fluxes, terrace and step diffusive currents, and has been incorporated into our PFC surface evolution model, which predicts the formation of preferentially aligned nanoridge stripe patterns on the PFC surface. We establish that these anisotropic diffusive currents accelerate nanotendril growth on the PFC surface and the onset of surface nanostructure pattern formation. We also explore systematically the dependence of the PFC surface morphological response on the surface temperature and He ion incident flux, characterize in detail the resulting surface topographies and growth kinetics, and compare the predicted surface morphologies with experimental observations. Our simulation predictions for the emerging surface nanostructure patterns under certain plasma exposure conditions are consistent with experimental findings in the literature.

36 MATERIALS SCIENCE↗

Using In Situ High-Energy X-ray Diffraction to Quantify Electrode Behavior of Li-Ion Batteries from Extreme Fast Charging

Extreme fast charging (XFC, =15 min charging time) of Li-ion batteries (LIBs) has been proposed as an immediate target to increase the commercial appeal of electric vehicles. However, XFC of LIBs is associated with the degradation of battery performance and safety concerns. Quantitative and simultaneous characterization of various components during cell degradation represents a major experimental challenge. In this work, we outline a methodology for the use of spatially resolved, high-energy X-ray diffraction as a quantitative, in situ method of mapping the degradation of LIBs. We use this approach to study the battery cell capacity loss, both locally (mm scale) and globally over the entire cell (cm scale). Specifically, our workflow allows us to quantify the total amount of plated Li on the anode, as well as its spatial correlation to the structural properties of the anode and cathode. The method complements existing optical methods to resolve the spatial heterogeneity of local degradation mechanisms such as Li plating and provides simultaneous insights into concomitant anode state-of-charge variability. We apply it to commercially relevant single-layer pouch cells with the graphite anode and the LiNi 0.5 Mn 0.3 Co 0.2 O 2 cathode. Our results show that Li plating occurs heterogeneously on the graphite anode and that it is spatially correlated to the extent of anode lithiation. In this work, we anticipate that the described workflow will allow for understanding multiscale degradation in energy-storage devices beyond LIBs, where quantitative analysis at a local and global length scale can be performed without the necessity to tear down the device, due to the applicability of high-energy X-rays to probe in situ degradation.

25 ENERGY STORAGE↗

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗

Multiphase, Multiscale Chemomechanics at Extreme Low Temperatures: Battery Electrodes for Operation in a Wide Temperature Range

Understanding the behavior of lithium-ion batteries (LIBs) under extreme conditions, for example, low temperature, is key to broad adoption of LIBs in various application scenarios. LIBs, poor performance at low temperatures is often attributed to the inferior lithium-ion transport in the electrolyte, which has motivated new electrolyte development as well as the battery preheating approach that is popular in electric vehicles. A significant irrevocable capacity loss, however, is not resolved by these measures nor well understood. Herein, multiphase, multiscale chemomechanical behaviors in composite LiNi x Mn y Co z O 2 (NMC, x + y + z = 1) cathodes at extremely low temperatures are systematically elucidated. The low-temperature storage of LIBs can result in irreversible structural damage in active electrodes, which can negatively impact the subsequent battery cycling performance at ambient temperature. Furthermore, beside developing electrolytes that have stable performance, designing batteries for use in a wide temperature range also calls for the development of electrode components that are structurally and morphologically robust when the cell is switched between different temperatures.

25 ENERGY STORAGE↗

Spatial Growth Modeling and High Resolution Remote Sensing Data Coupled with Air Quality Modeling to Assess the Impact of Atlanta, Georgia on the Local and Regional Environment

The growth of cities, both in population and areal extent, appears as an inexorable process. Urbanization continues at a rapid rate, and it is estimated that by the year 2025, 60 percent of the world s population will live in cities. Urban expansion has profound impacts on a host of biophysical, environmental, and atmospheric processes within an urban ecosystems perspective. A reduction in air quality over cities is a major result of these impacts. Because of its complexity, the urban landscape is not adequately captured in air quality models such as the Community Multiscale Air Quality (CMAQ) model that is used to assess whether urban areas are in attainment of EPA air quality standards, primarily for ground level ozone. This inadequacy of the CMAQ model to sufficiently respond to the heterogeneous nature of the urban landscape can impact how well the model predicts ozone levels over metropolitan areas and ultimately, whether cities exceed EPA ozone air quality standards. We are exploring the utility of high-resolution remote sensing data and urban spatial growth modeling (SGM) projections as improved inputs to a meteorological/air quality modeling system focusing on the Atlanta, Georgia metropolitan area as a case study. These growth projections include business as usual and smart growth scenarios out to 2030. The growth projections illustrate the effects of employing urban heat island mitigation strategies, such as increasing tree canopy and albedo across the Atlanta metro area, which in turn, are used to model how air temperature can potentially be moderated as impacts on elevating ground-level ozone, as opposed to not utilizing heat island mitigation strategies. The National Land Cover Dataset at 30m resolution is being used as the land use/land cover input and aggregated to the 4km scale for the MM5 mesoscale meteorological model and the CMAQ modeling schemes. Use of these data has been found to better characterize low density/suburban development as compared with USGS lkm land use/land cover data that have traditionally been used in modeling. Air quality prediction for future scenarios to 2030 is being facilitated by land use projections using a spatial growth model. Land use projections were developed using the 2030 Regional Transportation Plan developed by the Atlanta Regional Commission, the regional planning agency for the area. This allows the Georgia Environmental Protection Division to evaluate how these transportation plans will affect future air quality. The coupled SGM and air quality modeling approach provides insight on what the impacts of Atlanta s growth will be on the local and regional environment and exists as a mechanism that can be used by policy makers to make rational decisions on urban growth and sustainability for the metropolitan area in the future.

Quattrochi, Dale A.↗

Remote Sensing and Spatial Growth Modeling Coupled with Air Quality Modeling to Assess the Impact of Atlanta, Georgia on the Local and Regional Environment

The growth of cities, both in population and areal extent, appears as an inexorable process. Urbanization continues at a rapid rate, and it is estimated that by the year 2025, 80 percent of the world s population will live in cities. Directly aligned with the expansion of cities is urban sprawl. Urban expansion has profound impacts on a host of biophysical, environmental, and atmospheric processes. A reduction in air quality over cities is a major result of these impacts. Strategies that can be directly or indirectly implemented to help remediate air quality problems in cities and that can be accepted by political decision makers and the general public are now being explored to help bring down air pollutants and improve air quality. The urban landscape is inherently complex and this complexity is not adequately captured in air quality models, particularly the Community Multiscale Air Quality (CMAQ) model that is used to assess whether urban areas are in attainment of EPA air quality standards, primarily for ground level ozone. This inadequacy of the CMAQ model to sufficiently respond to the heterogeneous nature of the urban landscape can impact how well the model predicts ozone pollutant levels over metropolitan areas and ultimately, whether cities exceed EPA ozone air quality standards. We are exploring the utility of high-resolution remote sensing data and urban spatial growth modeling (SGM) projections as improved inputs to the meteorology component of the CMAQ model focusing on the Atlanta, Georgia metropolitan area as a case study. These growth projections include "business as usual" and "smart growth" scenarios out to 2030. The growth projections illustrate the effects of employing urban heat island mitigation strategies, such as increasing tree canopy and albedo across the Atlanta metro area, which in turn, are used to model how ozone and air temperature can potentially be moderated as impacts on elevating ground-level ozone, as opposed to not utilizing heat island mitigation strategies. The National Land Cover Dataset at 30m resolution is being used as the land use/land cover input and aggregated to the 4km scale for the MM5 mesoscale meteorological model and the (CMAQ) modeling schemes. Use of these data have been found to better characterize low density/suburban development as compared with USGS 1km land use/land cover data that have traditionally been used in modeling. Air quality prediction for future scenarios to 2030 is being facilitated by land use projections using a spatial growth model. Land use projections were developed using the 2030 Regional Transportation Plan developed by the Atlanta Regional Commission, the regional planning agency for the area. This allows the State Environmental Protection agency to evaluate how these transportation plans will affect future air quality. The coupled SGM and air quality modeling approach provides insight on what the impacts of Atlanta s growth will be on the local and regional environment and exists as a mechanism that can be used by policy makers to make rationale decisions on urban growth and sustainability for the metropolitan area in the future.

Quattrochi, Dale A.↗

Going mobile to address emerging climate equity needs in the heterogeneous urban environment

We report the Brookhaven National Laboratory Center for Multiscale Applied Sensing (CMAS) aims to address environmental equity needs in the context of a changing climate. As a first step towards this goal, the center developed a one-of-a-kind observatory tailored to the study of highly heterogeneous urban environments. This article describes the features of the mobile observatory that enable its rapid deployment either on or off the power grid, as well as its instrument payload. Beyond its unique design, the observatory optimizes data collection within the obstacle-laden urban environment using a new smart sampling paradigm. This setup facilitated the collection of previously poorly documented environmental properties including wind profiles throughout the atmospheric column. The mobile observatory captured unique observations during its first few intensive observation periods (IOPs).. Vertical air motion and infrared temperature measurements collected along the faces of the supertall One Vanderbilt skyscraper in Manhattan, NY reveal how solar and anthropogenic heating affect wind flow and thus the venting of heat, pollution, and contaminants in urban street canyons. Also, air temperature measurements collected during travel along a 150-km transect between Upton and Manhattan, NY offer a high-resolution view of the urban heat island and reveal that temperature disparities also exist within the city across different neighborhoods. Ultimately, the datasets collected by CMAS are poised to help guide equitable urban planning by highlighting existing disparities and characterizing the impact of urban features on the urban microclimate with the goal of improving human comfort.

54 ENVIRONMENTAL SCIENCES↗

Driving macro-scale transformations in three-dimensional-printed biopolymers through controlled induction of molecular anisotropy at the nanoscale

Motivated by the need to harness the properties of renewable and biodegradable polymers for the design and manufacturing of multi-scale structures with complex geometries, we have employed our additive manufacturing platform that leverages molecular self-assembly for the production of metre-scale structures characterized by complex geometries and heterogeneous material composition. As a precursor material, we used chitosan, a chemically modified form of chitin, an abundant and sustainable structural polysaccharide. We demonstrate the ability to control concentration-dependent crystallization as well as the induction of the preferred orientation of the polymer chains through the combination of extrusion-based robotic fabrication and directional toolpathing. Anisotropy is demonstrated and assessed through high-resolution micro-X-ray diffraction in conjunction with finite element simulations. Using this approach, we can leverage controlled and user-defined small-scale propagation of residual stresses to induce large-scale folding of the resulting structures.

36 MATERIALS SCIENCE↗

The Colorado East River Community Observatory Data Collection

Abstract The U.S. Department of Energy's (DOE) Colorado East River Community Observatory (ER) in the Upper Colorado River Basin was established in 2015 as a representative mountainous, snow‐dominated watershed to study hydrobiogeochemical responses to hydrological perturbations in headwater systems. The ER is characterized by steep elevation, geologic, hydrologic and vegetation gradients along floodplain, montane, subalpine, and alpine life zones, which makes it an ideal location for researchers to understand how different mountain subsystems contribute to overall watershed behaviour. The ER has both long‐term and spatially‐extensive observations and experimental campaigns carried out by the Watershed Function Scientific Focus Area (SFA), led by Lawrence Berkeley National Laboratory, and researchers from over 30 organizations who conduct cross‐disciplinary process‐based investigations and modelling of watershed behaviour. The heterogeneous data generated at the ER include hydrological, genomic, biogeochemical, climate, vegetation, geological, and remote sensing data, which combined with model inputs and outputs comprise a collection of datasets and value‐added products within a mountainous watershed that span multiple spatiotemporal scales, compartments, and life zones. Within 5 years of collection, these datasets have revealed insights into numerous aspects of watershed function such as factors influencing snow accumulation and melt timing, water balance partitioning, and impacts of floodplain biogeochemistry and hillslope ecohydrology on riverine geochemical exports. Data generated by the SFA are managed and curated through its Data Management Framework. The SFA has an open data policy, and over 70 ER datasets are publicly available through relevant data repositories. A public interactive map of data collection sites run by the SFA is available to inform the broader community about SFA field activities. Here, we describe the ER and the SFA measurement network, present the public data collection generated by the SFA and partner institutions, and highlight the value of collecting multidisciplinary multiscale measurements in representative catchment observatories.

54 ENVIRONMENTAL SCIENCES↗

Improved Detection System Description and New Method for Accurate Calibration of Micro-Channel Plate Based Instruments and Its Use in the Fast Plasma Investigation on NASA's Magnetospheric MultiScale Mission

The Fast Plasma Investigation (FPI) on NASAs Magnetospheric MultiScale (MMS) mission employs 16 Dual Electron Spectrometers (DESs) and 16 Dual Ion Spectrometers (DISs) with 4 of each type on each of 4 spacecraft to enable fast (30 ms for electrons; 150 ms for ions) and spatially differentiated measurements of the full 3D particle velocity distributions. This approach presents a new and challenging aspect to the calibration and operation of these instruments on ground and in flight. The response uniformity, the reliability of their calibration and the approach to handling any temporal evolution of these calibrated characteristics all assume enhanced importance in this application, where we attempt to understand the meaning of particle distributions within the ion and electron diffusion regions of magnetically reconnecting plasmas. Traditionally, the micro-channel plate (MCP) based detection systems for electrostatic particle spectrometers have been calibrated using the plateau curve technique. In this, a fixed detection threshold is set. The detection system count rate is then measured as a function of MCP voltage to determine the MCP voltage that ensures the count rate has reached a constant value independent of further variation in the MCP voltage. This is achieved when most of the MCP pulse height distribution (PHD) is located at higher values (larger pulses) than the detection system discrimination threshold. This method is adequate in single-channel detection systems and in multi-channel detection systems with very low crosstalk between channels. However, in dense multi-channel systems, it can be inadequate. Furthermore, it fails to fully describe the behavior of the detection system and individually characterize each of its fundamental parameters. To improve this situation, we have developed a detailed phenomenological description of the detection system, its behavior and its signal, crosstalk and noise sources. Based on this, we have devised a new detection system calibration method that enables accurate and repeatable measurement and calibration of MCP gain, MCP efficiency, signal loss due to variation in gain and efficiency, crosstalk from effects both above and below the MCP, noise margin, and stability margin in one single measurement. More precise calibration is highly desirable as the instruments will produce higher quality raw data that will require less post-acquisition data correction using results from in-flight pitch angle distribution measurements and ground calibration measurements. The detection system description and the fundamental concepts of this new calibration method, named threshold scan, will be presented. It will be shown how to derive all the individual detection system parameters and how to choose the optimum detection system operating point. This new method has been successfully applied to achieve a highly accurate calibration of the DESs and DISs of the MMS mission. The practical application of the method will be presented together with the achieved calibration results and their significance. Finally, it will be shown that, with further detailed modeling, this method can be extended for use in flight to achieve and maintain a highly accurate detection system calibration across a large number of instruments during the mission.

detection system↗

Toughening of Boron Carbide Composites by Hierarchical Microstructuring

Due to a unique combination of properties including high hardness, low density, chemical and thermal stability, semi-conductivity, and high neutron absorption, boron carbide (B4C) is a potential candidate for various applications involving extreme environment. However, B4C’s current application is limited because of its low fracture toughness. In this study, a hierarchical microstructure design with features including TiB2 grains and graphite platelets was used to toughen B4C by simultaneously utilizing multiple toughening mechanisms including crack deflection, bridging, and micro-crack toughening. Using field-assisted sintering technology (FAST), B4C composites with dense and hierarchical microstructure were fabricated. Previously, the fracture toughness of fabricated B4C composites was measured at micro-scale using micro-indentation to have up to 56% improvement. In this work, the B4C composites’ fracture toughness was characterized at macro-scale using four-point bending methods and compared with previous results obtained at micro-scale. Micromechanics modeling of fracture behaviors for B4C-TiB2 composites was also performed to evaluate the contributions from experimentally observed toughening mechanisms. From four-point bending tests, B4C composites reinforced with both TiB2 grains (~15 vol%) and graphite platelets (~8.7 vol%) exhibited the highest fracture toughness enhancement from 2.38 to 3.65 MPa∙m^(1/2). The measured values were lower than those obtained using micro-indentation but maintained the general trends. The discrepancy between the indentation and four-point bending test results originated from the complex deformation behaviors triggered by the high contact load during indentation tests. Through micromechanics modeling, introduced thermal residual stress due to thermal expansion mismatch between B4C and TiB2, and weak interphases at B4C-TiB2 boundaries were identified as the main causes for experimentally observed toughness enhancement. These results proved the effectiveness of hierarchical microstructure designs for B4C toughening and can provide reference for the future design of B4C composites with optimized microstructures for further fracture toughness enhancement.

composites↗

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

ExaAM: Metal additive manufacturing simulation at the fidelity of the microstructure

Additive manufacturing (AM), or 3D printing, of metals is transforming the fabrication of components, in part by dramatically expanding the design space, allowing optimization of shape and topology. However, although the physical processes involved in AM are similar to those of welding, a field with decades of experimental, modeling, simulation, and characterization experience, qualification of AM parts remains a challenge. The availability of exascale computational systems, particularly when combined with data-driven approaches such as machine learning, enables topology and shape optimization as well as accelerated qualification by providing process-aware, locally accurate microstructure and mechanical property models. We describe the physics components comprising the Exascale Additive Manufacturing simulation environment and report progress using highly resolved melt pool simulations to inform part-scale finite element thermomechanics simulations, drive microstructure evolution, and determine constitutive mechanical property relationships based on those microstructures using polycrystal plasticity. We report on implementation of these components for exascale computing architectures, as well as the multi-stage simulation workflow that provides a unique high-fidelity model of process–structure–property relationships for AM parts. In addition, we discuss verification and validation through collaboration with efforts such as AM-Bench, a set of benchmark test problems under development by a team led by the National Institute of Standards and Technology.

3D printing↗

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗