Search NASA⌕ Search

SEARCH · Search NASA

Results for “composite porous media”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multi-stage preconditioners for thermal–compositional–reactive flow in porous media

We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.

42 ENGINEERING↗
DOE OSTI2020

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗
DOE OSTI2023

A state-of-the-art review of experimental and computational studies of granular materials: Properties, advances, challenges, and future directions

Modeling of heterogeneous materials and media is a problem of fundamental importance to a wide class of phenomena and systems, ranging from condensed matter physics, soft materials, and composite media to porous media, biological systems, geosystems, ceramic engineering, pharmaceutical science and even in space discoveries. Among the most important materials are granular systems, which have received intense interest from the engineering, physics, and mathematics communities. In this review paper, the recent developments and new advances in experimental, and computational methods on a variety of scales and physics that extend understanding to a wide range of materials and phenomena are reviewed. Experimental advances include computed neutron and nanometer-scale tomography, magnetic resonance imaging, refractive index matching, digital image correlation, acoustic emission analysis, and the most recent 4D techniques. Furthermore, a tremendous shift has occurred from the continuum scale to micro-scale and developing multiscale approaches. As such, various computational methods, including, constitutive modeling, discrete modeling, and multiscale approaches, have been developed. In conclusion, aside from all these evolutions, more complicated modeling called coupled, or multiphysics, systems representing a simultaneous presence of heat, fluid, chemical variation, and mechanical effect are also explored.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Co-Transport of PFAS and Cr(VI) in porous media

PFAS and Cr are present at some sites as co-contaminants. Here, the objective of this research was to investigate the co-transport behavior of per- and polyfluoroalkyl substances (PFAS) and hexavalent chromium (Cr(VI)) in porous media. Miscible-displacement experiments were conducted using two soils and an aquifer sediment with different geochemical properties. Perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) were employed as model PFAS. The retardation of PFOS was decreased in the presence of Cr(VI). Conversely, the transport and retardation of PFOA was not affected by the presence of Cr(VI). The reduction of PFOS retardation caused by Cr(VI) is likely due to sorption competition for both organic-carbon and inorganic (metal-oxides and clay minerals) domains. The relative contributions of the three soil constituents to PFOS sorption and the potential for competition between PFOS and Cr(VI) is a function of the geochemical composition of the porous media (i.e., organic carbon, metal-oxides and clay minerals). The PFAS had minimal impact on the retention and transport of Cr(VI). To our knowledge, the results presented herein represent the first reported data for PFOS and Cr(VI) co-transport in porous media. The results of this study indicate that the presence of Cr(VI) has the potential to increase the migration potential of PFOS in soil and groundwater, which should be considered when characterizing electroplating facilities, leather tanning facilities, and other co-contaminated sites.

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

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗
DOE OSTI2020

Estimating diffusion coefficients of shale oil, gas, and condensate with nano-confinement effect

Molecular diffusion plays an important role in oil and gas migration in tight shale formations. However, there are insufficient reference data in the literature to specify the diffusion coefficients within porous media. This study aims at calculating diffusion coefficients of shale gas, shale condensate, and shale oil at reservoir conditions with CO 2 injection for EOR/EGR. The large nano-confinement effects including large gas-oil capillary pressure and critical property shifts could alter the phase behaviors. This paper estimates the diffusivities of shale fluids in nanometer-scale shale rock from two perspectives: 1) examining the shift of diffusivity caused by nanopore confinement effects from phase change (phase composition and fluid property) perspective, and 2) calculating the effective diffusion coefficient in porous media by incorporating rock intrinsic properties (porosity and tortuosity factor). The tortuosity is obtained by using tortuosity-porosity relations as well as the measured tortuosity of shale from 3D imaging techniques. The results indicated that nano-confinement effects could affect the diffusion coefficient through altering the phase properties, such as phase compositions and densities. Compared to bulk phase diffusivity, the effective diffusion coefficient in porous shale rock is reduced by 10 2 to 10 4 times as porosity decreases from 0.1 to 0.03.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

An experimental investigation on the CO 2 storage capacity of the composite confining system

Assuring secure containment of stored CO 2 is of paramount importance—for climate change mitigation, for permitting, and for reassuring the public. Regional seals, such as those sealing petroleum accumulations, have proven to be effective for securing CO 2 . However, the goal of permanent sequestration can also be satisfied with “composite confining systems” consisting of multiple, possibly discontinuous flow barriers that, in aggregate, create a system with very high permeability anisotropy and effectively retard the vertical migration of CO 2 . This study focuses on investigating the barrier characteristics necessary for effective containment of CO 2 , using both physical flow experiments and modified invasion percolation simulations. The simulations are calibrated to the physical experiments and are used to further extend the analysis. Results show that for a composite confining system, a) even barriers with low capillary entry pressure contrast to the underlying flow unit can divert rising CO 2 ; b) curved or anticlinal barrier topography can enhance CO 2 trapping; c) fining-upward gradations make little difference to barrier effectiveness or CO 2 retention; d) longer barriers retain more CO 2 regardless of barrier topography. Finally, field-scale simulations have demonstrated the importance of barrier length for increasing CO 2 storage capacity. Furthermore, the results presented can be used to inform the development of new screening criteria for characterization and effectiveness of composite confining systems.

58 GEOSCIENCES↗
DOE OSTI2021

Modulating Interactions between Molten Polystyrene and Porous Solids Using Atomic Layer Deposition

Understanding and modulating the interactions between molten polymers and porous solids is important for numerous processes and phenomena including catalytic conversion of polymers and fabrication of nanocomposites and nanostructured materials. Although changing the surface composition of pores would enable modulation of interactions between polymer and nanoporous solids, it is challenging to achieve such a control without inducing significant changes to the size and structure of nanopores. In this work, we demonstrate that the interactions between molten polystyrene (PS) and disordered packings of SiO 2 nanoparticles (NPs) can be modulated by changing the surface composition of the NPs using atomic layer deposition (ALD). A disordered packing of silica NPs is modified with varying surface coverages of TiO 2 , WO 3 and CaCO 3 , with coverages estimated by the mass gain and the refractive index change of NP packings. Based on the time required to fully infiltrate these ALD-modified NP packings via capillarity, the contact angles for PS on different surfaces prepared via ALD are determined. The contact angle gradually changes from that of pure SiO 2 to that of the fully covered surfaces. The contact angles for PS on SiO 2 , TiO 2 , WO 3 and CaCO 3 are found to be 20°, 62°, 70° and 10°, respectively. Interestingly, the contact angles and interfacial energies between PS and the ALD-modified surfaces do not correlate strongly with the water contact angle of these surfaces; thus, caution must be exercised in predicting how a polymer would wet or interact with porous solids solely based on their hydrophilicity. Furthermore, the method presented in this work can be extended to study the interactions between a wide range of polymers and surfaces in porous media, which will have important implications for designing new catalytic materials for polymer upcycling reactions as well as novel NP-polymer composite films and membranes with enhanced mechanical and transport properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The hard x-ray nanotomography microscope at the advanced light source

Beamline 11.3.1 at the Advanced Light Source is a tender/hard (6–17 keV) x-ray bend magnet beamline recently re-purposed with a new full-field, nanoscale transmission x-ray microscope. The microscope is designed to image composite and porous materials possessing a submicrometer structure and compositional heterogeneity that determine materials’ performance and geologic behavior. The theoretical and achieved resolutions are 55 and <100 nm, respectively. The microscope is used in tandem with a <25 nm eccentricity rotation stage for high-resolution volume imaging using nanoscale computed tomography. The system also features a novel bipolar illumination condenser for the illumination of an ~100 μm spot of interest on the sample, followed by a phase-type zone plate magnifying objective of ~52 µm field of view and a phase detection ring. The zone plate serves as the system objective and magnifies the sample with projection onto an indirect x-ray detection system, consisting of a polished single crystal CsI(Tl) scintillator and a range of high-quality Plan Fluorite visible light objectives. The objectives project the final visible light image onto a water-cooled CMOS 2048 × 2048-pixel 2 detector. Here, in this paper, we will discuss the salient features of this instrument and describe early results from imaging the internal three-dimensional microstructure and nanostructure of target materials, including fiber-reinforced composites and geomaterials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Effects of temperature inversion on densification in chemical vapor infiltration

In a classical chemical vapor infiltration (CVI) process, the competing effects of chemical kinetics and reagent gas transport lead to non–uniform depositions such that outer layers of a preform densify faster leaving the core highly porous. Currently, CVI must be performed at a sufficiently low temperature to achieve good densification quality which leads to high processing time and cost. Volumetric heating of the preform, especially through microwaves, can create temperature inversion such that the core is hotter than the outer surface and potentially, overcome the challenges associated with isothermal CVI. Direct numerical simulations (DNS) of densification under various such temperature distributions indicate that microwave heating in CVI processing can lead to better (uniform) densification of porous preforms. Here the role of key parameters describing the temperature distributions on the densification behavior is investigated. Strategic temporal control of the temperature distribution shows that processing times can be reduced by almost half while maintaining a good densification quality similar to that of low–temperature isothermal processing. Inside–out densification due to the inverted temperature profile is a key distinguishing characteristic of microwave assisted CVI.

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

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗
DOE OSTI2021

Alumina-based filters made via binder jet 3D printing of alumina powder, colloidal silica infiltration, and sintering

Alumina-based, porous filter media was made via a binder jet 3D printing process consisting of an alumina powder printing step with subsequent heating, colloidal silica infiltration, drying, and sintering to consolidate particles yet retain a net open porous microstructure. The composites made were alumina-silica or alumina-mullite, where the silica sintering aid was used to densify and join the alumina particles. The resulting composite structures had open porosities in the 25–31 vol% range as measured by Archimedes density. Pressure drops were measured across the filter media at constant flow rates to compare disc shapes and complex, 3D printed filters based on the N95 design requirements. Complex, 3D-printed alumina composites were produced with acceptable pressure drops for N95 implementation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Synthesis of PHI-type zeolite harmotome using barium as an inorganic structure-directing agent inspired from natural mineral compositions

Zeolites are nanoporous aluminosilicates widely used as commercial adsorbents, heterogeneous catalysts, and ion-exchange materials due to their unique porosity, acidity, and thermal stability. The vast majority of synthetically-realized zeolite structures are prepared under hydrothermal conditions in alkaline media, and often in the presence of an organic structure-directing agent that facilitates the crystallization of diverse porous networks. Due to economic and environmental disadvantages of organic-based syntheses, it is often desirable to produce zeolites in organic-free media using alkali metal ions as the most commonly employed inorganic structure-directing agents (ISDAs). Here, in this study, we were inspired by the elemental compositions of natural zeolite minerals to examine the use of an alkaline earth metal ion as the ISDA to prepare zeolites with the PHI framework type. Two natural minerals exhibiting this crystal structure are phillipsite and harmotome. Both structures contain sodium and potassium as extra-framework cations, but only harmotome contains barium cations. Here we examined how the combination of alkali and alkaline earth ISDAs promote the formation of PHI-type zeolites. Our findings reveal that each isostructure is prepared under different synthesis conditions where simply adding barium to a growth mixture does not always lead to harmotome formation. Harmotome is formed in more dilute and alkaline growth media; however, time-resolved analyses of zeolite crystallization reveal that harmotome can be synthesized at low temperatures (ca. 65 °C) with rates that are 5-times faster than syntheses of phillipsite. This seemingly indicates that barium reduces the energetic barrier for PHI-type zeolite nucleation in ways that are not fully understood. Moreover, harmotome crystals exhibit an average particle size of ca. 60 nm, which is nearly 40-fold smaller than phillipsite crystals prepared under similar conditions. Collectively, these findings identify routes to control crystallization of PHI-type zeolites and may provide a broader strategy for using the composition of natural zeolite minerals to optimize the preparation of synthetic analogues.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulating the non-monotonic strain response of nanoporous multiferroic composites under electric field control

In this work, we simulate and analyze the mechanical response of a class of multiferroic materials consisting of a templated porous nanostructure made out of cobalt ferrite (CFO) partially filled by atomic layer deposition (ALD) with a ferroelectric phase of lead zirconate titanate (PZT). The strain in the device is measured when an electric field is applied for varying ALD thicknesses, displaying a non-monotonic dependence with a maximum strain achieved for a coating thickness of 3 nm. To understand this behavior, we apply finite element modeling to the smallest repeatable unit of the nanoporous template and simulate the mechanical response as a function of PZT coating thickness. We find that this non-monotonic response is caused by the interplay between two driving forces opposing one another. First, increased porosity works toward increasing the strain due to a reduced system stiffness. Second, decreased porosity involves a larger mass fraction of PZT, which drives the electro-mechanical response of the structure, thus leading to a larger strain. The balance between these two driving forces is controlled by the shear coupling at the CFO/PZT interface and the effective PZT cross section along the direction of the applied electric field. Here, our numerical results show that considering a nonlinear piezoelectric response for PZT leads to an improved agreement with the experimental data, consistent with ex situ poling of the nanostructure prior to magnetic measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Model-based optimization strategies for direct hydrogenation of carbon dioxide to dimethyl ether

Here, this study discusses model-based optimization strategies for CO 2 hydrogenation to dimethyl ether (DME) over aCuZnZr (CZZ) and ferrite (FER) mixed catalyst system in a packed-bed reactor configuration. A two-dimensional axisymmetric, nonisothermal packed-bed reactor model was developed using COMSOL Multiphysics 6.2 software. The model solves two-dimensional (radial and axial) heat and mass transport equations in the packed-bed and integrates intraparticle diffusion and heat transfer in a 1D approach. This powerful feature differs from a traditional porous media approach and takes into account any heat and mass transfer limitations that may exist. Analysis shows that the heat transfer limitations are negligible, but strong internal mass transfer limitations were observed at 10 ≤ WHSV ≤ 90 h -1 on the FER catalyst and at 240 °C. The optimum catalyst composition (i.e., mixing ratio) varies depending on the operating regime. The FER catalyst weight in the mixture can be as low as 5 wt.%, but the ideal composition de-pends on the internal mass transfer limitation and its relationship with the operating regime (i.e., weight hourly space velocity, temperature). A catalyst composition of 80 wt.% CZZ and 20 wt.% FER was suggested; this composition can provide high CO 2 conversion and DME production rates at a wide range of temperatures and flow rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗