Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “capillary pressure”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fundamentals of Using Cracked Film Lithography to Pattern Transparent Conductive Metal Grids for Photovoltaics

The fundamentals of using cracked film lithography (CFL) to fabricate metal grids for transparent contacts in solar cells were studied. The underlying physics of drying-induced cracks were well-predicted by an empirical correlation relating crack spacing to capillary pressure. CFL is primarily controlled by varying the crack template thickness, which establishes a three-way tradeoff between the areal density of cracks, crack width, and spacing between cracks, which in turn determine final grid transmittance, grid sheet resistance, and the semiconductor resistance for a given solar cell. Since CFL uses a lift-off process, an additional constraint is that the metal thickness must be less than 1/3 of the crack template thickness. Furthermore, the transmittance/grid sheet resistance/wire spacing tradeoffs measured in this work were used to calculate solar cell performance: CFL-patterned grids should outperform screen-printed grids for narrow cells (0.5–2 cm wide) and/or cells with high semiconductor sheet resistance (=100 O/sq), making CFL attractive for monolithically integrated thin-film photovoltaic modules.

14 SOLAR ENERGY↗
DOE OSTI2020

Scaling of Imbibition Front Dynamics in Heterogeneous Porous Media

Here we report a quantitative study on the dynamics of countercurrent imbibition in fractured rock samples with the application of CO 2 geosequestration, as the rock matrix is saturated with CO 2 and brine exists in the adjacent microfracture. Direct pore-scale simulations of spontaneous and forced imbibition in a three-dimensional heterogeneous porous matrix enable us to capture the effect of microscopic properties on the imbibition dynamics as these pore-level descriptions are ignored in the existing coarse-grained models. We have established new scaling classes for the propagation of the center of mass of the imbibition front and fluid-fluid interface broadening. It has become evident that the wettability is a major factor which leads to anomalies in these scaling laws. While the governing capillarity under the water-wet condition induces fast roughening dynamics of the interface, the adverse effect of capillary pressure under the intermediate-wet condition leads to a damping effect on the imbibition front fluctuations.

58 GEOSCIENCES↗
DOE OSTI2021

Prototype equation of state for phase transition of confined fluids based on the generalized van der Waals partition function

A simple self-consistent prototype equation of state (EOS) based on the Generalized van der Waals (vdW) partition function has been demonstrated to describe the phase transition of simple fluids in nanopores with uniform size. Different from those commonly presented in the literature, the new EOS does not need an auxiliary equation that is conventionally applied to provide the capillary pressure derived from surface tension. Lastly, the encouraging performance of the EOS calls for further extension to applications with more complex fluids and porous media.

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

Understanding Carbon Dioxide Transfer in Direct Methanol Fuel Cells Using a Pore-Scale Model

Abstract The gas flow of carbon dioxide from the catalyst layer (CL) through the microporous layer (MPL) and gas diffusion layer (GDL) has great impacts on the water and fuel management in direct methanol fuel cells (DMFCs). This work has developed a liquid–vapor two-phase model considering the counter flow of carbon dioxide gas, methanol, and water liquid solution in porous electrodes of DMFC. The model simulation includes the capillary pressure as well as the pressure drop due to flow resistance through the fuel cell components. The pressure drop of carbon dioxide flow is found to be about two to three orders of magnitude higher than the pressure drop of the liquid flow. The big difference between liquid and gas pressure drops can be explained by two reasons: volume flowrate of gas is three orders of magnitude higher than that of liquid; only a small fraction of pores (<5%) in hydrophilic fuel cell components are available for gas flow. Model results indicate that the gas pressure and the mass transfer resistance of liquid and gas are more sensitive to the pore size distribution than the thickness of porous components. To buildup high gas pressure and high mass transfer resistance of liquid, the MPL and CL should avoid micro-cracks during manufacture. Distributions of pore size and wettability of the GDL and MPL have been designed to reduce the methanol crossover and improve fuel efficiency. The model results provide design guidance to obtain superior DMFC performance using highly concentrated methanol solutions or even pure methanol.

Electrochemistry↗
DOE OSTI2020

Editors’ Choice—Flooded by Success: On the Role of Electrode Wettability in CO 2 Electrolyzers that Generate Liquid Products

Economic operation of carbon dioxide (CO 2 ) electrolyzers generating liquid products will likely require high reactant conversions and product concentrations, conditions anticipated to challenge existing gas diffusion electrodes (GDEs). Notably, electrode wettability will increase as lower surface tension products (e.g., formic acid, alcohols) are introduced into electrolyte streams, potentially leading to flooding. To understand the hydraulically stable operating envelopes in mixed aqueous-organic liquid domains, we connect intrinsic electrode wettability descriptors to operating parameters such as electrolyte flow rate and current. We first measure contact angles of water-organic product dilutions on polytetrafluoroethylene (PTFE) and graphite surfaces as planar analogues for GDE components. We then use material balances around the reactive gas-liquid-solid interface to calculate product mass fractions as functions of water sweep rate and current. Product composition maps visualize the extent to which changes in cell performance influence capillary pressure, a determinant of GDE saturation. Analyses suggest that formic acid mixtures pose little risk for GDE flooding across a wide range of conditions, but effluents containing <30% alcohol by mass may cause flooding. This study reveals opportunities to integrate microstructural features and oleophobic surface treatments into GDEs to repel aqueous-organic mixtures and expand the window of stable operating conditions.

25 ENERGY STORAGE↗
DOE OSTI2022

Effects of Carbonic Acid-Rock Interactions on CO2/Brine Multiphase Flow Properties in the Upper Minnelusa Sandstones

Summary Carbon dioxide (CO2) injection into a deep saline aquifer can dissolve into formation brine and generate carbonic acid. The resulting acid can drive fluid-rock geochemical reactions. The impact of these fluid-rock geochemical reactions on porosity, permeability, and multiphase flow responses is relevant to the determination of CO2 storage capacity of deep saline aquifers. In this research, carbonic acid flooding experiments were performed on core samples consisting of poorly sorted, quartz-rich sand with laminated bedding from a possible CO2 storage target in northwest Wyoming. Complementary pre- and post-injection porosity and permeability, thin-section, Brunauer-Emmett-Teller (BET) surface area, mercury intrusion capillary pressure (MICP), and time-domain nuclear magnetic resonance (TD-NMR) measurements were conducted. Overall, both core porosity and permeability increased after a 7-day carbonic acid injection, from 6.2 to 8.4% and 1.6 to 3.7 md, respectively. We attributed these changes to carbonate mineral dissolution, which was evidenced by the effluent brine geochemistry, pore-throat size distribution (PTSD), and BET surface area. To be more specific, within the more-permeable section of core samples containing larger pore size, the permeability increment is apparent due to dolomite mineral grains and cements dissolution. However, for the lower-permeability section corresponding to the smaller pore size, mineral precipitation possibly lessened dissolution effects, leading to insignificant petrophysical properties changes. Consequently, the observed heterogeneous carbonic acid-rock interactions resulted in alterations of CO2/brine relative permeability (i.e., the initial CO2 saturation decreased and the CO2 flow capacity was enhanced). This research provides a fundamental understanding regarding effects of fluid-rock reactions on changes in static and multiphase flow properties of eolian sandstones, which lays the foundation for more accurate prediction/simulation of CO2 injection into deep saline aquifers.

Engineering↗
DOE OSTI2020

Crystalline Disposal R&D at LBNL: FY20 Progress Report

Within the Spent Fuel & Waste Science and Technology (SFWST) Program, research work continues further with the goal to better understand long-term performance of disposal systems in three main geologic rock types: clay/shale, salt, and crystalline rock. This report documents LBNL’s research activities related to investigations of crystalline host rock according to the scope of two work packages: SF-20LB01030207 “Crystalline International Collaborations – LBNL,” and SF-20LB01030203 “Crystalline Disposal R&D – LBNL.” These research activities correspond are related to key Features, Events, and Processes (FEPs), ranked in importance from medium to high, as listed in Table 7 of the Used Fuel Disposition Campaign Disposal Research and Development Roadmap (FCR&D-USED-2011-000065 REV0) (Nutt, 2011). Specifically, these research activities address FEP 2.2.01, Excavation Disturbed Zone (EZD). The results of these research activities provide important insights into understanding and predicting flow and transport processes that could occur in low-permeability crystalline rocks, in which fractures might serve as main conduits for fluid flow and radionuclide transport. The evolution of the EDZ during the excavation of the tunnel as well as the evolution of microcrack growth within the EDZ after the emplacement of backfill are critical for predicting the long- term behavior of the EDZ. A number of factors including stress, temperature, water activity, capillary pressure, chemistry, and mineralogy can affect the rock evolution, which require advanced experimental tools to study it.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Calculation of Friction Factors in Heat Pipes Using CFD in Support of the Sockeye Code

A heat pipe is two-phase heat transfer device which relies on surface tension and capillary pressure to provide a very efficient heat transfer mechanism. Currently, there is interest in designing micro nuclear power reactors using high temperature heat pipes to provide totally passive cooling. To support this, Sockeye is being developed as a heat pipe simulation capability. To support the development of Sockeye, individual phenomena must be tested. This work focuses on applying a CFD model to evaluate friction factors for use in Sockeye. A CFD model has been developed in Nek5000 which simulates each component fluid, i.e. sodium liquid and vapor, independently and captures the effects of the wick on friction factor. The model tests the effect of the wick permeability and the thickness of the wick occupied by each fluid component. The thickness of the wick occupied by each component can be directly related to the local vapor volume fraction from the Sockeye heat pipe model. Results from the CFD model indicate that friction factor in the liquid flow decreases with increasing wick permeability, while the friction factor in the vapor flow increases with increasing wick permeability. In the liquid, this is likely a direct result of the increased cross sectional area available for the flow. For the vapor, this is explained as the wick acting like a rough wall, while the increased area has only a negligible impact. Results for the tests of the thickness of the wick occupied by each fluid showed a thresholding behavior for both liquid and vapor. When the thickness of the wick occupied by each fluid was increased, it increased the effect of the wick, but only up to a certain point. Once a certain thickness was reached, the friction factors remained constant. When the fluid occupied only a very small fraction of the wick, friction factors returned to nominal values.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Ordovician-Cambrian Units: Hierarchical Evaluation of Geologic Carbon Storage Resource Estimates

The Indiana Geological and Water Survey (IGWS) led subtask 1.1 to assess the regional distribution and estimate the storage capacity of Ordovician-Cambrian stratigraphic units located within the partnership region. A comprehensive data set of wireline logs and petrophysical information was used to generate these interpretations. These data include core analysis for porosity and permeability, mercury injection capillary pressure (MICP), and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for three potential storage reservoirs(limestone and dolostone from the Upper Ordovician Trenton Limestone/Black River Group and equivalent units, the Middle Ordovician St. Peter Sandstone, and primary target reservoir rocks of the Lower Ordovician and Upper Cambrian Knox Supergroup and equivalent units) calculated using six methodologies: (1) a fixed value of porosity of 10 percent in all units evaluated; (2) a unique average porosity (per well) from wireline-derived porosity (neutron, sonic, and/or density porosity for each unit); (3) porosity values from core analysis; (4) a depth-dependent porosity model (Knox Supergroup only); (5) porosity based on a model based on petrophysical facies; and (6) SREs using National Energy Technology Laboratory’s CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN beta V2). All methods used the same values for thickness for each unit. However, the areal extent of each assessment was limited by the data available for each method. Estimated volumes were calculated in 1-by-1 kilometer grid cells and summarized as county and total stratigraphic unit volumes. The resultant SREs mass are displayed using boxplots, which allow for comparing data statistics (mean values and variability) between methods. Differences observed in SRE results from the six methods are mainly attributable to differences in the data and conceptual models used to interpret or estimate porosity in each method. Based on this systematic variability between methods, it is inferred that methods 1, 4, and 6 are best used for regional-scale reconnaissance estimates of storage capacity while methods 2, 3, and 5 are more appropriate for local scales where more data is required. All estimates are data-density dependent and different methods require different amounts of data for reasonable assessments. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Spontaneous Imbibition Tests and Parameter Estimation in Volcanic Tuff

We present a dynamic laboratory spontaneous imbibition test and interpretation method, demonstrated on volcanic tuff samples from the Nevada National Security Site. The method includes numerical inverse modeling to quantify uncertainty of estimated two-phase fluid flow properties. As opposed to other approaches requiring multiple different laboratory instruments, the dynamic imbibition method simultaneously estimates capillary pressure and relative permeability from one test apparatus.

36 MATERIALS SCIENCE↗
DOE OSTI2021

GDSA PFLOTRAN Development (FY2021)

The Spent Fuel & Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is to develop a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2021 advances of the PFLOTRAN Development group of the SFWST Campaign. The mission of this group is to develop a geologic disposal system modeling capability for nuclear waste that can be used to probabilistically assess the performance of generic disposal concepts. In FY 2021, development proceeded along three main thrusts: software infrastructure, code performance, and process model advancement. Software infrastructure improvements included implementing an Agile software development framework and making improvements to the QA Test Suite. Code performance improvements included development of advanced linear and nonlinear solvers as well as design of flexible smoothing algorithms for capillary pressure functions. Process modeling advancements included the addition of flexible thermal conductivity function definitions and refinement of multi-continuum reactive transport to support Sandia’s participation in DECOVALEX

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Simple Heat Pipe Model

This is a simple model designed to run fast but still maintain the key physics and feedback mechanisms of a heat pipe. First, the capillary pressure is a function of the liquid working fluid volume fraction. Second, the boiling and condensation are based on the saturation temperature that is based on the heat pipe pressure. When the pressure goes up, the saturation temperature goes up and the vapor rains on the wick. When the pressure goes down, the saturation temperature goes down and the liquid in the entire wick boils. This is how the heat pipe adjusts to stay robust under different temperatures and heat fluxes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

FY2022 Progress on Imbibition Testing in Containment Science

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as for vapor-phase contaminant transport and/or containment of noble gases in the subsurface. In this second progress report of FY22, we present two ongoing activities related to imbibition testing on volcanic rock samples. We present the development of a new analytical solution predicting the temperature response observed during imbibition into dry samples, as discussed in our previous first progress report for FY22. We also illustrate the use of a multi-modal capillary pressure distribution to simulate both early- and late-time imbibition data collected on tuff core that can exhibit multiple pore types. These FY22 imbibition tests were conducted for an extended period (i.e., far beyond the time required for the wetting front to reach the top of the sample), which is necessary for parameter estimation and characterization of two different pore types within the samples.

42 ENGINEERING↗
DOE OSTI2022

The Cypress Sandstone Seal System

The Cypress Sandstone is the youngest and shallowest unit in the Illinois Basin that was featured in the United States Carbon Utilization and Storage Atlas IV as a target for saline carbon storage with an estimated 0.2 to 2.3 GT of storage potential. Additional research on a residual oil zone (ROZ) developed within the Cypress Sandstone has delineated 27 prospects with approximately 290.8 million m3 (1.8 billion barrels) of oil in place. 21 to 31 million m 3 (144 to 196 million barrels) of oil is estimated to be recoverable using carbon dioxide enhanced oil recovery (CO 2 -EOR). Storage of CO 2 associated with EOR in these ROZ prospects alone, not accounting for associated main pay zones (MPZs), underlying brine formation, or intervals adjacent to or between prospects, is estimated to be up to 10.4 billion tonnes. The Cypress Sandstone is thus well understood to be a CO 2 injection target, both for EOR and associated storage. However, the seal system overlying the Cypress Sandstone is poorly understood. Unlike deeper formations such as the Mt. Simon Sandstone or the St. Peter Sandstone which are either in use as a CO 2 sink or being characterized for prospective storage, respectively, the Cypress is not overlain by hundreds of feet of impermeable shale. Rather, the Cypress is overlain by a lithologically variable interval that is composed generally of shales and limestones with some sandstone in the part of the Basin where the Cypress is deep enough to facilitate CO 2 storage. Also, due to its status as one of the shallowest and most prolific oil reservoirs in the Basin, the seal system overlying the Cypress Sandstone has a relatively high number of legacy well penetrations. The purpose of this report is to characterize the Cypress Sandstone seal system using well logs and available core. Gross thickness, lithology (facies), and mineralogy of seals is described and mapped across the Basin. The column height of CO 2 that can be held is calculated using capillary pressure data from a representative core.

02 PETROLEUM↗
DOE OSTI2025

CalTestBed - Twelve Benefit Co - Understanding Degradation of Water Management Properties in Gas Diffusion Layers (CRADA Final Report)

The research between Lawrence Berkeley National Lab and Twelve Benefit Co. enhanced Twelve’s understanding of how chemical changes to the gas diffusion layer lead to failure of their CO 2 electrolyzers. During the CRADA period, LBNL team tested Twelve’s gas diffusion layer samples (GDL), some with microporous layer (MPL) and some without MPL using the water-air capillary pressure setup. LBNL team also tested microporous layer samples using the fuel cell test stands and compared those results with commercial GDL and MPL. Understanding the hydrophobicity of the material, along with electrochemical performance is the key to determining Twelve’s ability to scale their platform. Further improvements to Twelve’s CO 2 electrolyzer can enhance the net zero emissions with their technology.

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

PFLOTRAN Development FY2021

The Spent Fuel & Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is to develop a disposal system modeling and analysis capability for valuating disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2021 advances of the PFLOTRAN Development group of the SFWST Campaign. The mission of this group is to develop a geologic disposal system modeling capability for nuclear waste that can be used to probabilistically assess the performance of generic disposal concepts. In FY 2021, development proceeded along three main thrusts: software infrastructure, code performance, and process model advancement. Software infrastructure improvements included implementing an Agile software development framework and making improvements to the QA Test Suite. Code performance improvements included development of advanced linear and nonlinear solvers as well as design of flexible smoothing algorithms for capillary pressure functions. Process modeling advancements included the addition of flexible thermal conductivity function definitions and refinement of multi-continuum reactive transport to support Sandia’s participation in DECOVALEX. This report fulfills the GDSA PFLOTRAN Development Work Package Level 3 Milestone – PFLOTRAN Development, FY2021, M3SF-21SN010304072.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗