Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CT Scanning”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

X-ray Computed Tomography Data of Dense Metallic Components

The data shared in here are X-ray computed tomography (XCT) scans of a hexagonal fuel nozzle in 3 sections with the Metrotom 800 system at the Manufacturing Demonstration Facility (MDF) at Oak Ridge National Laboratory. The data are used in the paper "Tomographic Sparse View Selection using the View Covariance Loss, by Lin et al. (doi:10.1109/TPAMI.2025.36000720), accepted to the international conference on computational imaging (ICCP 2025). Figures 4-7 in the paper describe the part/XCT scan. File name Descriptions: Bottom section: TCR- Single Channeled SRC L 2019-3-18 12-26-41.hdf5 Medium section: TCR- Single Channeled SRC M 2019-3-18 13-8-9.hdf5 Top section: TCR- Single Channeled SRC T 2019-3-18 13-45-39.hdf5 Each hdf5 file contains projection data, and all the relevant X-ray CT scan setting. The full list of included attributes: distance_unit: Units of all distances specified angle_unit : Units of the angles angles: Array of all angles used voxel_size_xy: Baseline recon (if any) has this voxel size in the in-plane direction voxel_size_z: Baseline recon (if any) has this voxel size in the cross-plane direction det_pixel_size_col: Size of the detector pixels in the column dimension det_pixel_size_row: Size of the detector pixels in the row dimension src_iso_dist: Source to iso-center distance iso_det_dist: Iso-center to detector distance det_angle: If the detector is rotated/tilted, this angle corresponds to that value det_row_offset: Center of rotation offset in the vertical direction det_col_offset: Center of rotation offset in the horizontal direction reconstruction: A baseline reconstruction stored as 3D array BHC params: Beam-hardening parameters - Van De Casteel Model - if it has been used to pre-process the projections We also provided a python script (hdf_io.py) that allows the user to read the relevant data from each hdf5 file.

Ziabari, Amir [Oak Ridge National Laboratory]↗
DOE OSTI2021

Reactivity of CO 2 with Utica, Marcellus, Barnett, and Eagle Ford Shales and Impact on Permeability

We report that in order to reduce greenhouse gas emissions while recovering hydrocarbons from unconventional shale formations, processes that make use of carbon dioxide to enhance oil recovery while storing carbon dioxide (CO 2 ) should be considered. Here, we examine samples from three shale basins across the United States (Utica and Marcellus Shales in the Appalachian Basin, Barnett Shale in the Bend Arch-Ft. Worth Basin, and Eagle Ford in the Western Gulf Basin) to address the following questions: (1) do changes from reaction with CO 2 and fluids at the micrometer and nanometer scale alter flow pathways and, in turn, impact hydrocarbon production, CO 2 storage, and seal integrity and (2) can CO 2 or fluid reactivity be predicted based on physical or chemical properties of shale formations? Experiments were conducted at 40 °C and 10.3 MPa to characterize the interaction between CO 2 and shale using X-ray diffraction (XRD), carbon and sulfur analysis, in situ Fourier transform infrared spectroscopy (FT-IR), feature relocation scanning electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), mercury (Hg) intrusion porosimetry, and Brunauer-Emmett-Teller (BET) surface area and pore size analysis coupled with density functional theory (DFT) methods. Changes in mechanical, physical, and flow properties of shale cores due to CO 2 exposure were addressed using a New England Research Autolab 1500 and Xenon X-ray computed tomography (CT) scanning. Results showed that CO 2 did not promote significant reactivity with the shale if water was not present; only shales with swelling clays or residual interstitial pore water reacted with dry CO 2 to promote reactivity in shale. When water was added as a reactant, CO 2 formed carbonic acid and reacted with the shale to dissolve carbonate pockets, etched and pitted the shale matrix surfaces, and increased the microporosity and decreased nanoporosity. Porosity and permeability increased appreciably in core shale samples after exposure to CO 2 saturated fluid due to dissolution of carbonate. Shale mechanical properties were not altered. Trends were not observed that could tie CO 2 or fluid reactivity to physical or chemical properties of the shale formations at the basin scale from the samples we examined. However, if the shale contained significant amounts of carbonate and water was available to react with the CO 2 , pore sizes were altered in the matrix and permeability and porosity increased.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2025

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗
DOE OSTI2022

Model-based reconstruction for enhanced x-ray CT of dense tri-structural isotropic particles

Tri-structural isotropic (TRISO) fuel particles are a key component of next generation nuclear fuels. Using x-ray computed tomography (CT) to characterize TRISO particles is challenging because of the strong attenuation of the x-ray beam by the uranium core, leading to severe photon starvation in a substantial fraction of the measurements. Furthermore, the overall acquisition time for a high-resolution CT scan can be very long when using conventional laboratory-based x-ray systems and reconstruction algorithms. Specifically, when analytic methods such as the Feldkamp–Davis–Kress (FDK) algorithm are used for reconstruction, it results in severe streak artifacts and noise in the corresponding 3D volume, which makes subsequent analysis of the particles challenging. In this paper, we develop and apply model-based image reconstruction (MBIR) algorithms to improve the quality of CT reconstructions for TRISO particles to facilitate better characterization. We demonstrate that the proposed MBIR algorithms can significantly suppress artifacts with minimal pre-processing compared to conventional approaches. We also demonstrate that the proposed MBIR approach can obtain high-quality reconstruction compared to the FDK approach even when using a fraction of the typically acquired measurements, thereby enabling dramatically faster measurement times for TRISO particles.

Venkatakrishnan, Singanallur↗
DOE OSTI2022

Hydro-mechanical behavior of heated bentonite buffer for geologic disposal of high-level radioactive waste: A bench-scale X-ray computed tomography investigation

The behavior of heated bentonite buffer is critical for the security and long-term performance of a geological repository for high-level radioactive waste (HLW). While laboratory column experiments have been conducted to investigate compacted bentonite and coupled THMC (thermal-hydro-mechanical and chemical) processes for a moderate temperature range of up to 100 °C, data for a higher temperature range are limited. Understanding bentonite behavior and coupled THMC processes under higher temperatures (e.g., up to 200 °C) could allow for a more economic repository design and would expand the data and knowledge base for more reliable modeling. In this study, a bench-scale experiment was conducted in a compacted bentonite column experiencing both heating up to 200 °C in the center and hydration from a sand-clay boundary surrounding the column. During the experiment run for 1.5 years, frequent X-ray computed tomography (CT) scanning of bentonite provided insights into the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. After the experiment, a comprehensive post-dismantling characterization of bentonite samples was conducted. Results showed that the bentonite hydration was axi-symmetrical despite the initial heterogeneity due to packing, confirming the ability of bentonite to seal fast flow/transport paths. Further, compared to a non-heated control experiment, the heated column showed greater CT density variations along the radial distance, indicating that homogenization of bentonite might be more difficult if a temperature gradient is maintained in the repository. Precipitation of an anhydrite layer occurred in the inner hot zone, pointing to potential concerns about salt precipitation causing canister corrosion. Ultimately, the experiments provided a high-resolution window into the strongly dynamic and coupled behavior of bentonite exposed to heating, hydration and swelling, which will be valuable for improving modeling of coupled processes, especially for the early state of a HLW repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1999

Quantifying Void Ratio Variation in Sand using Computed Tomography

A series of displacement-controlled, conventional, drained axisymmetric (triaxial) experiments were conducted on dry Ottawa sand specimens at very low effective confining stresses in a microgravity environment aboard the Space Shuttle during the NASA STS-89 mission. Post-flight analysis included studying the internal fabric and failure patterns of these specimens using Computed Tomography (CT). The CT scans of three specimens subjected to different compression levels (uncompressed specimen, a specimen compressed to 3.3% nominal axial strain (epsilon(sub a)), and a specimen compressed to 25% epsilon(sub a)) are presented to investigate the evolution of instability patterns and to quantify void ratio variation. The progress of failure is described and discussed. Also, specimens' densities were calibrated using standard ASTM procedures and void ratio spatial variation was calculated and represented by contour maps and histograms. The CT technique demonstrated good ability to detect specimen inhomogeneities, localization patterns, and quantifying void ratio variation within sand specimens.

Alshibli, Khalid A.↗
DOE OSTI2025

Distributed Stochastic Optimization of a Neural Representation Network for Time-Space Tomography Reconstruction

4D time-space reconstruction of dynamic events or deforming objects using X-ray computed tomography (CT) is an important inverse problem in non-destructive evaluation. Conventional back-projection based reconstruction methods assume that the object remains static for the duration of several tens or hundreds of X-ray projection measurement images (reconstruction of consecutive limited-angle CT scans). However, this is an unrealistic assumption for many in-situ experiments that causes spurious artifacts and inaccurate morphological reconstructions of the object. To solve this problem, we propose to perform a 4D time-space reconstruction using a distributed implicit neural representation (DINR) network that is trained using a novel distributed stochastic training algorithm. Our DINR network learns to reconstruct the object at its output by iterative optimization of its network parameters such that the measured projection images best match the output of the CT forward measurement model. Here, we use a forward measurement model that is a function of the DINR outputs at a sparsely sampled set of continuous valued 4D object coordinates. Unlike previous neural representation architectures that forward and back propagate through dense voxel grids that sample the object's entire time-space coordinates, we only propagate through the DINR at a small subset of object coordinates in each iteration resulting in an order-of-magnitude reduction in memory and compute for training. DINR leverages distributed computation across several compute nodes and GPUs to produce high-fidelity 4D time-space reconstructions. We use both simulated parallel-beam and experimental cone-beam X-ray CT datasets to demonstrate the superior performance of our approach.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

A Procedure to Quantify Interlaminar Damage in Refractory Composites by In-Situ Micro-X-Ray Computed Tomography

Refractory composites are relevant for hot-structure applications on re-entry vehicles. To develop optimized high-performance hot structures, accurate material property data is required. Material properties that often limit design in two-dimensional (2-D) composites are interlaminar strengths. A better understanding of the failure mechanisms for interlaminar shear (ILS), interlaminar tension (ILT), and curved-beam interlaminar tension (CB-ILT) tests can provide insight on improving the design and fabrication of refractory composites. In-situ tests were performed to capture the progressive failure in American Society for Testing and Materials (ASTM) sized ILS (ASTM D3846), ILT (ASTM D7291), and CB-ILT (ASTM D6415) specimens. Before measuring internal damage of refractory composites specimens, the methodology of the test setup warrants a complete and thorough discussion. Detail on the test procedure, test stand design, micro-x-ray computed tomography (CT) configuration, data processing, and CT system uncertainty is discussed. Through keeping the many parameters of the CT system consistent, the uncertainty of the CT system was minimized, optimizing the quality of measurements. The overall quality of a CT scan can be quantified by three metrics: voxel resolution, minimal voxel resolution feature detection, and directional variability. A better understanding of the progressive failure mechanisms of refractory composites is provided by the results of both the present discussion and the one found in the accompanying paper titled “Analysis of Interlaminar Damage in Refractory Composites by In-Situ Micro-X-Ray Computed Tomography”.

In-situ micro-x-ray CT↗
NASA NTRSDate not supplied

Characterization and Modeling of Cell Wall Imperfections in Aluminum Honeycomb Cores using X-ray CT Imaging

Honeycomb cell imperfections may significantly affect the failure of sandwich composites in core crushing and shear modes. Identification and quantification of the different types of imperfections and evaluation of their effects on the mechanical response of honeycomb cores in flatwise compression is presented. X-Ray computed tomography (CT) scans of honeycomb cores in co-cured composite sandwich panels were performed. Image analysis methods are applied to characterize and quantify core cell geometry and cell wall waviness. An approach is developed for extracting cell wall mid surface planes and waviness for constructing finite element models using shell elements. Analysis results from models with the as-manufactured imperfections are compared to results from models with idealized geometries to evaluate the effect of the as-manufactured imperfections on the initiation of cell wall buckling and subsequent core crushing.

aluminum honeycomb↗
NASA NTRS2021

Application of Radon Transform to Multi-Angle Measurements Made by the Research Scanning Polarimeter: A New Approach to Cloud Tomography. Part I: Theory and Tests on Simulated Data.

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP was a prototype for the Aerosol Polarimetry Sensor (APS) launched on-board the NASA Glory satellite. Currently the retrieval algorithms developed for the RSP are being adopted for the measurements of the space-borne polarimeters on the upcoming NASA’s Plankton, Aerosol, Cloud Ocean Ecosystem (PACE)satellite mission. The RSP’s uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid water cloud droplet sizes using the polarized rainbow structure. It also provides geometric constraints on the cumulus cloud’s 2D cross section yielding the cloud’s geometric shape estimates. In this study we further build on the latter technique to develop a new tomographic approach to retrieval of cloud internal structure from remote sensing measurements. While tomography in the strict definition is a technique based on active measurements yielding a tomogram (directional optical thickness as a function of angle and offset of the view ray), we developed a “semi-tomographic” approach in which tomogram of the cloud is estimated from passive observations instead of being measured directly. This tomogram is then converted into 2D spatial distribution of the extinction coefficient using inverse Radon transform (filtered back projection) which is the standard tomographic procedure used e.g., in medical CT scans. This algorithm is computationally inexpensive compared to techniques relying on highly-multi-dimensional least-square fitting; it does not require iterative 3D RT simulations. The resulting extinction distribution is defined up to an unknown constant factor, so we discuss the ways to calibrate it using additional independent measurements. In the next step we use the profile of the droplet size distribution parameters from the cloud’s side (derived by fitting the polarized rainbows) to convert the 2D extinction distribution into that of the droplet number concentration. We illustrate and validate the proposed technique using 3D-RT-simulatedRSP observations of a LES-generated Cu cloud. Quantitative comparisons between the retrieved and the original optical and microphysical parameters are presented.

clouds↗
DOE OSTI2021

Assessment of a tomography-informed polyhedral discrete element modelling approach for complex-shaped granular woody biomass in stress consolidation

The design of handling equipment for granular biomass primarily requires experiments. Discrete element models (DEM) can provide designers with detailed insight into the behaviour of granular materials. However, granular biomass comprising complex-shaped particles is difficult to model with DEM. A tomography-informed DEM approach and an exhaustive assessment of the ability of the approach to predict the bulk behaviour of milled pines using experimental data is presented. Nano-CT scan was conducted to obtain 3D particle surface geometries as the basis for particle shape approximation by a polyhedral model and a sphero-polyhedral model. These models were applied in the simulation of a compressibility test. Our parametric study showed that particle Young's modulus and restitution coefficient are the two main properties influencing the simulated bulk behaviour of the DEM particles and that the level of approximating DEM particle shape for real particles is critical for accurately replicating the bulk behaviour of milled pines. The polyhedral model demonstrated better suitability than its sphero-polyhedral counterpart for modelling pine particles. The polyhedral model calibrated in a compressibility test was then applied in the simulations of a friction test without additional parameter tuning. Interlocking was dominant in the shear of bulk pine particles. Remarkably, the polyhedral model predicted the frictional behaviour of the pine particles when compared to the experimental results. The limitations of the model, as well as possible ways for enhancement, are discussed. This work provided novel insights into the suitability of complex-shaped DEM models for granular woody biomass.

Biomass flow↗
DOE OSTI2023

Migration and Precipitation of Platinum in Anion‐Exchange Membrane Fuel Cells

Abstract Despite the recent progress in increasing the power generation of Anion‐exchange membrane fuel cells (AEMFCs), their durability is still far lower than that of Proton exchange membrane fuel cells (PEMFCs). Using the complementary techniques of X‐ray micro‐computed tomography (CT), Scanning Electron Microscopy (SEM) and Energy Dispersive X‐ray (EDX) spectroscopy, we have identified Pt ion migration as an important factor to explain the decay in performance of AEMFCs. In alkaline media Pt +2 ions are easily formed which then either undergo dissolution into the carbon support or migrate to the membrane. In contrast to PEMFCs, where hydrogen cross over reduces the ions forming a vertical “Pt line” within the membrane, the ions in the AEM are trapped by charged groups within the membrane, leading to disintegration of the membrane and failure. Diffusion of the metal components is still observed when the Pt/C of the cathode is substituted with a FeCo−N−C catalyst, but in this case the Fe and Co ions are not trapped within the membrane, but rather migrate into the anode, thereby increasing the stability of the membrane.

Raut, Aniket↗
DOE OSTI2023

Migration and Precipitation of Platinum in Anion–Exchange Membrane Fuel Cells

Despite the recent progress in increasing the power generation of Anion-exchange membrane fuel cells (AEMFCs), their durability is still far lower than that of Proton exchange membrane fuel cells (PEMFCs). Using the complementary techniques of X-ray micro-computed tomography (CT), Scanning Electron Microscopy (SEM) and Energy Dispersive X-ray (EDX) spectroscopy, we have identified Pt ion migration as an important factor to explain the decay in performance of AEMFCs. In alkaline media Pt+2 ions are easily formed which then either undergo dissolution into the carbon support or migrate to the membrane. In contrast to PEMFCs, where hydrogen cross over reduces the ions forming a vertical “Pt line” within the membrane, the ions in the AEM are trapped by charged groups within the membrane, leading to disintegration of the membrane and failure. Furthermore, diffusion of the metal components is still observed when the Pt/C of the cathode is substituted with a FeCo–N–C catalyst, but in this case the Fe and Co ions are not trapped within the membrane, but rather migrate into the anode, thereby increasing the stability of the membrane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Ability of x‐ray computed tomography to resolve critical flaw size in laser‐based, paste stereolithography ceramic printing of alumina

Abstract Complex alumina parts were printed using vat photopolymerization (VPP), which is a stereolithography‐based additive manufacturing (AM) technique used to shape ceramic preforms, or green parts. The critical flaw size was determined using classical fracture mechanics techniques. The strength and fracture toughness were measured and compared to flaws detected in x‐ray computed tomography (XCT or CT) distributions as well as the fracture surfaces. The strength was lower compared traditionally made alumina, and that is due to layering effects, slurry defects, and printing defects. The critical flaw size from fracture mechanics was 206 µm. XCT has high enough resolution to detect the critical flaw size and much smaller features, where the average flaw size observed in CT scans was around 80–100 µm. The fracture surfaces indicate that flaws causing failure are larger than that of the critical flaw size (∼300 µm), but fracture surfaces do not show definitive features compared to traditionally made ceramics. Since XCT can observe flaws smaller than the critical flaw size, this method can be used as a screening technique.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A deep learning framework for layer-wise porosity prediction in metal powder bed fusion using thermal signatures

Abstract Part quality manufactured by the laser powder bed fusion process is significantly affected by porosity. Existing works of process–property relationships for porosity prediction require many experiments or computationally expensive simulations without considering environmental variations. While efforts that adopt real-time monitoring sensors can only detect porosity after its occurrence rather than predicting it ahead of time. In this study, a novel porosity detection-prediction framework is proposed based on deep learning that predicts porosity in the next layer based on thermal signatures of the previous layers. The proposed framework is validated in terms of its ability to accurately predict lack of fusion porosity using computerized tomography (CT) scans, which achieves a F1-score of 0.75. The framework presented in this work can be effectively applied to quality control in additive manufacturing. As a function of the predicted porosity positions, laser process parameters in the next layer can be adjusted to avoid more part porosity in the future or the existing porosity could be filled. If the predicted part porosity is not acceptable regardless of laser parameters, the building process can be stopped to minimize the loss.

42 ENGINEERING↗
DOE OSTI2023

Path-Integrated X-Ray Images for Multi-Surface Digital Image Correlation (PI-DIC)

X-ray imaging offers unique possibilities for Digital Image Correlation (DIC), opening the door for full-field deformation measurements of a test article in complex environments where optical DIC suffers severe biases or is impossible. While X-ray DIC has been performed in the past with standard DIC codes designed for optical images, the path-integrated nature of X-ray images places constraints on the experimental setup, predominantly that only a single surface of interest moves/deforms. These requirements are difficult to realize for many practical situations and limit the amount of information that can be garnered in a single test. Other X-ray based diagnostics such as Digital Volume Correlation (DVC) and Projection DVC (P-DVC) overcome these obstacles, but DVC is limited to quasi-static tests, and both DVC and P-DVC necessitate high-resolution computed tomography (CT) scan(s) and often require a potentially invasive pattern throughout the volume of the specimen. Here this work presents a novel approach to measure time-resolved displacements and strains on multiple surfaces from a single series of 2D, path-integrated (PI) X-ray images, called PI-DIC. The principle of optical flow or conservation of intensity—the foundation of DIC—was reframed for path-integrated images, for an exemplar setup comprised of two plates moving and deforming independently. Synthetic images were generated for rigid translations, rigid rotations, and uniform stretches, where each plate underwent a unique motion/deformation. Experimental specimens were fabricated (either an aluminum plate with tantalum features or a plastic plate with steel features) and the two specimens were independently translated. PI-DIC was successfully demonstrated with the synthetic images and validated with the experimental images. Prescribed displacements were recovered for each plate from the single set of path-integrated, deformed images. Errors were approximately 0.02 px for the synthetic images with 1.5% image noise, and 0.05 px for the experimental images. These results provide the foundation for PI-DIC to measure motion and deformation of multiple, independent surfaces with subpixel accuracy from a single series of path-integrated X-ray images.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Multiscale porosity characterization in additively manufactured polymer nanocomposites using micro-computed tomography

Extrusion-based additive manufacturing (AM) of polymer composites exhibits complex thermally driven phenomena that introduce severe discontinuities in the internal structure across length scales, especially voids or porosity. This study utilizes a high-throughput porosity characterization technique to analyze large datasets from numerous micro-computed tomography (mu CT) scans to capture the influence of AM print parameters on the size, shape, and location of porosity across multiple print layers (up to a few cm) on fused granular fabrication (FGF) printers. The materials investigated include nanocomposite formulations based on commercially relevant nylon-12 and polyether ketone ketone (PEKK) materials comprising nano- or micro- sized fillers. The estimated global porosity follows an inverse linear correlation against the bulk density of the printed samples. Increasing the extrusion multiplier (EM) and the nozzle temperature while decreasing the print speeds reduces the global porosity. Outlier analyses (local porosity morphology) show that faster print speeds and higher extrusion rates result in long, slender inter-layer voids, while lower nozzle temperatures lead to large, symmetrical, inter-bead voids (at the bead junction). Lack of active chamber temperature increases inter-layer and intra-bead voids with a two-fold increase in global porosity. Overall, the micro filler-reinforced composites exhibit higher global porosity than nanofiller-reinforced composites, which is attributed to the increased mismatch in the thermal expansion coefficient between the filler and the polymers used in the study.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗