Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porous flow”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Extending Acoustic Liner Bandwidth with Simple Embedded Septa

An acoustic liner design incorporating single-hole perforate septa was evaluated by the NASA Langley Liner Physics Team. The design, termed the ‘Simplified Septa’ concept, incorporates embedded septa as found in traditional, multidegree-of-freedom liners. Generally, such septa are composed of multihole perforates or porous mesh, whereas this concept uses septa with one hole per cell. An optimization scheme was employed to determine septa placement and hole diameters based on a target impedance cost function. Versions of the concept are tested in the NASA Langley Normal Incidence Tube (NIT) to determine impedance and absorption spectra for swept tonal excitation at 120 and 140 dB. Experimental results are compared to predictions from a liner model based on the Zwikker-Kosten Transmission Line (ZKTL) code and show good agreement for no-flow conditions.

acoustics↗
NASA NTRSDate not supplied

Spectral Proper Orthogonal Decomposition of uPSP Measurements in Recent NASA Ames Wind Tunnel Test

This paper discusses Spectral Proper Orthogonal Decomposition (SPOD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements in recent NASA Ames wind tunnel test. The uPSP measurements were collected using Innovative Scientific Solutions, Inc. (ISSI) porous, fast-response pressure-sensitive paint, 40 ISSI four-inch air-cooled Light-Emitting Diodes, and 8 Phantom v2512 high-speed cameras at 10,000 frames per second in the uPSP Launch Vehicle Demonstration Test (LVDT) of the Space Launch System (SLS) vehicle in the 11-by 11-foot transonic test section of the Unitary Plan Wind Tunnel at NASA Ames Research Center in April 2024. SPOD is derived from a space-time proper orthogonal decomposition problem for statistically stationary flows. SPOD modes are determined in the frequency domain. Each SPOD mode oscillates at a single frequency. SPOD can be viewed as an extension of the Discrete Fourier Transform composition and the Dynamic Mode Decomposition. In this paper, the outputs of SPOD of the uPSP measurements in the uPSP LVDT are presented and the effectiveness of SPOD in the identification, diagnosis and analysis of the aerodynamic and aeroacoustic phenomena is demonstrated. The unsteady and dynamic property of the pressure field on the surface of the SLS Block 1B crew vehicle is presented with the visualization of the SPOD modes of the uPSP measurements in the tests of a Mach sweep run of the uPSP LVDT. The SPOD outputs were generated with the execution in parallel of a code in Python, with the library of Message Passing Interface for parallel processing, on the NASA Pleiades supercomputer. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities Portfolio Office.

Aeroacoustics↗
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
NASA NTRS1997

Smoldering News From STS-77 Endeavour

The Microgravity Smoldering Combustion (MSC) experiment lifted off aboard the Space Shuttle Endeavour for its second flight in May 1996, as part of the STS-77 mission. This experiment is part of a series of studies focused on the smolder characteristics of porous combustible materials in a microgravity environment. Smoldering is a nonflaming form of combustion that takes place in the interior of combustible materials. Common examples of smoldering are nonflaming embers, charcoal briquettes, and cigarettes. The objective of this study is to provide a better understanding of the controlling mechanisms of smoldering in microgravity and normal Earth gravity (1g). As with other forms of combustion, gravity affects the availability of air and transport of heat, and therefore, the rate of combustion. The results of the microgravity experiments will be compared with identical ones carried out in 1g. In addition, they will be used to verify present theories of smolder combustion and will provide new insights into the process of smoldering combustion, enhancing our fundamental understanding of this frequently encountered combustion process and guiding improvements in fire safety practices. Two smoldering combustion tests with polyurethane foam were successfully accomplished during the STS-77 mission. The tests investigated smoldering combustion in a quiescent (no-flow) enriched oxygen environment, and in an air environment with a 2-mm/sec airflow through the fuel sample. The primary data from the tests are the ignition characteristics, spread rate, smolder reaction temperature, and products of combustion (solid and gas). On both the first mission on STS-69 and the second mission on STS-77, a smolder front propagated the length of the forced-flow samples, with the spread rate between the corresponding upward and downward 1g smolder rates. Neither of the quiescent cases propagated combustion (the first case was due in part to a problem with the experiment electronics). These results show a dramatic difference from the normal gravity results, where smolder propagation is very rapid and complete for both of these conditions. The experiment was conceived by Prof. A. Carlos Fernandez-Pello at the University of California-Berkeley. The MSC hardware was designed and built at the NASA Lewis Research Center by a team of civil servants and contractors from NYMA, Inc., and Aerospace Design & Fabrication, Inc. (ADF). The hardware consists of two sealed aluminum combustion chambers (each being a half a cylinder). The chambers hold the MSC test section, data acquisition electronics, power distribution electronics, and instrumentation. The hardware is fitted into a 5-ft(sup 3) Get-Away-Special (GAS) canister that is mounted in the shuttle cargo bay. The test section consists of a quartz cylinder that contains the polyurethane foam sample and an igniter. This igniter, which is an electrically heated wire sandwiched between two porous ceramic disks, is mounted in contact with the end of the foam sample. An array of 12 thermocouples placed axially and radially along the foam sample provide temperature histories, which are used to determine the rate of smolder propagation, and the characteristics of the reaction.

Koudelka, John M.↗
DOE OSTI2020

Water-Soluble Flexible Organic Frameworks That Include and Deliver Proteins

Four water-soluble hydrazone-based three-dimensional (3D) flexible organic frameworks FOF-1-4 have been synthesized from a semirigid tetracationic tetraaldehyde and four flexible dihydrazides. 1 H NMR spectroscopy indicated the quantitative formation of FOF-1-4 in D 2 O, while dynamic light scattering experiments revealed that, depending on the concentration, these porous frameworks display hydrodynamic diameters ranging from 50 to 120 nm. The porosity of the frameworks is confirmed by ethanol vapor adsorption experiments of the solid samples as well as the high loading capacity for a 2.3 nm porphyrin guest in water. The new water-soluble frameworks exhibit low cytotoxicity and form inherent pores with diameters of 5.3 or 6.7 nm, allowing rapid inclusion of proteins such as bovine serum albumin and green and orange fluorescent proteins, and efficient delivery of the proteins into normal and cancer cells. Flow cytometric analysis reveals percentages of the delivered cells up to 99.8%.

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

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2023

Pore-Scale Analysis of the Permeability Damage and Recovery during Cyclic Freshwater and Brine Injection in Porous Media Containing Non-Swelling Clays

Permeability damage in subsurface porous media caused by clay mobilization is encountered in many engineering applications, such as geothermal energy, water disposal, oil recovery, and underground CO2 storage. During the freshwater injection into rocks containing brine, the sudden decrease in salinity causes native clay fines to detach and clog pore throats, leading to a significant decline in permeability. The clay fines detach due to weakened net-attractive forces binding them to each other and the grain. Past experiments link this permeability damage on the immediate history of the salinity and the direction of flow. To better understand this phenomenon, we conducted pore-scale simulations of cyclic injection of freshwater and brine into sandstone containing Kaolinite clay. Our simulations establish a link between the clay-fine trajectory and the permeability trend observed by Khilar and Fogler (1983). For a uniform clay size of 3 microns, we observe a permeability decline by two orders of magnitude during freshwater injection with respect to brine injection. Increasing salinity and simultaneously reversing flow direction restores the permeability. The permeability restoration upon reversing the brine flow direction is attributed to the unblocking of pore throats in the reverse direction by the movement of the clay particles along the grain surfaces by the hydrodynamic force and the strong net-attractive force under high salinity.

Bhuvankar, Pramod (ORCID:000000027739973X)↗
NASA NTRS1998

Pump, and earth-testable spacecraft capillary heat transport loop using augmentation pump and check valves

A spacecraft includes heat-generating payload equipment, and a heat transport system with a cold plate thermally coupled to the equipment and a capillary-wick evaporator, for evaporating coolant liquid to cool the equipment. The coolant vapor is coupled to a condenser and in a loop back to the evaporator. A heated coolant reservoir is coupled to the loop for pressure control. If the wick is not wetted, heat transfer will not begin or continue. A pair of check valves are coupled in the loop, and the heater is cycled for augmentation pumping of coolant to and from the reservoir. This augmentation pumping, in conjunction with the check valves, wets the wick. The wick liquid storage capacity allows the augmentation pump to provide continuous pulsed liquid flow to assure continuous vapor transport and a continuously operating heat transport system. The check valves are of the ball type to assure maximum reliability. However, any type of check valve can be used, including designs which are preloaded in the closed position. The check valve may use any ball or poppet material which resists corrosion. For optimum performance during testing on Earth, the ball or poppet would have neutral buoyancy or be configured in a closed position when the heat transport system is not operating. The ball may be porous to allow passage of coolant vapor.

Baker, David↗
NASA NTRS1973

Surface brightness temperatures at the Apollo 17 heat flow site - Thermal conductivity of the upper 15 cm of regolith

Lunar surface brightness temperatures derived as part of the Apollo 17 heat flow experiment are reported. Nighttime surface temperatures, calculated from the data provided by two thermocouples suspended about 15 cm above the lunar surface, are used to determine the conductivity profile of the upper 15 cm of regolith at the ALSEP site. The surface reaches a maximum temperature of 384 (plus or minus 6) K at lunar noon and cools to a minimum temperature of 102 (plus or minus 1.5 K) at the end of the lunar night. Conductivities of the order of .000015 W/cm-deg K are postulated for a 2-cm porous surface layer overlying more compact regolith material with conductivities in the range of .0001 to .00015 W/cm-deg K between 2 and 15 cm. A mean surface temperature of 216 (plus or minus 5) K is deduced from the thermocouple data. The 256 K temperature measured by the probe sensors at 130 cm thus indicates that a large mean temperature gradient exists at the Apollo 17 site.

Keihm, S. J.↗
DOE OSTI2020

Thermal Analysis Evaluations Using the Dry Cask Simulator

Different thermal analysis models were developed to simulate the dry cask simulator (DCS). The DCS is an experiment designed to simulate dry storage of a single boiling water reactor fuel assembly under a variety of heat loads and internal pressures. The DCS was set up and tested in both a vertical and horizontal configuration to determine cladding temperatures in vertical and horizontal dry cask storage systems. The models included a detailed STAR-CCM+ model with the fuel assembly geometry explicitly modeled, a porous STAR-CCM+ model with the fuel assembly geometry modeled as porous media region with calculated effective properties, and a COBRA-SFS model. COBRA-SFS is a thermal-hydraulic code developed for steady-state and transient analysis of multi assembly spent-fuel storage and transportation systems. STAR-CCM+ is a commercial computation fluid dynamics (CFD) code. Both a detailed and porous STARCCM+ model were developed to look at the effective thermal conductivity (keff) approach to modeling a fuel assembly. A keff fuel model is typically model in CFD thermal analyses due to its significantly lower computational costs. The models were run for a combination of low and high canister pressures (100 kPa and 800 kPa) and low and high internal heat loads (0.5 kW and 5 kW). Results from all three models were compared against experimental data taken from the DCS for the peak cladding temperature (PCT) and inlet air mass flow.

Suffield, Sarah R.↗
DOE OSTI2020

Predicting Saturated Hydraulic Conductivity over Time for Degrading Saltstone Vault Concrete - 20376

At the Savannah River Site (SRS), low-level liquid tank waste is mixed with cementitious materials to produce a solid, monolithic waste form referred to as saltstone. The liquid saltstone mixture is pumped into massive concrete structures referred to as saltstone disposal units (SDUs) for curing and long-term storage. SDUs are vaults made of high performance concrete forming a barrier between saltstone and the surrounding environment. Predicting the long-term performance of SDU vault concrete, as well as the saltstone material itself, is critical to understanding the potential releases of radionuclides to the surrounding environment as described in the Saltstone Disposal Facility (SDF) performance assessment (PA). The hydraulic conductivity profile for SDU vault concrete was conservatively assumed to degrade linearly from the initial saturated hydraulic conductivity of the intact concrete over a long period of time. The initial hydraulic properties of saturated SDU vault concrete have been measured, and the final hydraulic properties of degraded vault concrete are assumed to be the same as the porous backfill material that will surround the vault after closure. The assumed final condition of the vault concrete is also conservative in that it provides a greater estimate of water penetration through the vault than would generally be expected. Therefore, insights that provide more realistic but conservative estimates as to how degradation processes may change the hydraulic properties of the SDU vault concrete, the time dependent relationship of the vault concrete hydraulic properties to degradation (and thus time), and the final hydraulic properties of the vault concrete provides useful additions to the SDF PA. This paper focuses on the assumption that the saltstone vault concrete saturated hydraulic conductivity would decrease linearly as a function of degradation time; this assumption resulted from a worst-case, weighted (by relative layer thickness), arithmetic average of intact and backfill material hydraulic conductivities representing flow parallel to a layered system. This research indicates that the effective saturated hydraulic conductivity of the degraded SDU vault concrete should be based on the relative thickness-weighted geometric mean of the intact and degraded layers hydraulic conductivities. This alternative, which is supported extensively in the relevant literature and evaluated in this paper, leads to a more realistic and defensible estimate of the effective saturated hydraulic conductivity in degraded vault concrete while remaining conservative. The impact of using the more realistic estimate of the effective saturated hydraulic conductivity for the degraded concrete decreases the effective saturated hydraulic conductivity by up to several orders of magnitude for early times (i.e., when the SDU vault concrete is assumed degraded via sulfate attack and carbonation). This revised estimate of the effective saturated hydraulic conductivity has been implemented in the recently revised SDF PA. (authors)

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

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗
DOE OSTI2022

Perspectives on Microfluidics for the Study of Asphaltenes in Upstream Hydrocarbon Production: A Minireview

The utilization of microfluidics has generated deep insights into asphaltene precipitation mechanisms and oil–water emulsion stabilization. Agglomeration and precipitation of asphaltenes can cause flow assurance problems during the extraction and transportation of crude oil. Change in temperature, pressure, reservoir conditions, and solvents can change the local environment, leading to asphaltene precipitation. Understanding asphaltene properties and precipitation pathways becomes critical in devising mitigation methods, demulsifiers, and suitable conditions during hydrocarbon processing. Microfluidics has helped in high throughput measurement studies, understanding critical processing conditions, fast demulsifier screening, and the effect of solvent concentration on deposition, generating useful information for utilization at the point of resource extraction facilitating improved resource management. It has become possible to capture the porous, complex nature of reservoir formations and the interaction of chemicals during precipitation through integrated analytics and visualization studies available only through microfluidics. The use of droplet microfluidics, with optical microscopy and high-speed imaging to study the oil–water interface, has resulted in greater understanding of the role of asphaltenes in interfacial properties and emulsion stabilization. Here, this minireview highlights the crucial aspects of microfluidics that have been used to understand physicochemical behavior and dynamics of asphaltene deposition. Some of the unique devices have been presented focusing on the key elements of microfluidics design, fabrication, and analysis, as the insight obtained from microfluidics strongly depends on the device design and the controllability of the experimental parameters. Successful implementation of microfluidics for efficient and controlled experiments, short analysis time scales and rapid screening, and generation of high-quality, reliable data that convey asphaltene deposition issues and interface behavior in emulsions shows the importance of microsystems for advancing knowledge in hydrocarbon production and processing.

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

Comparison of microwave and conventional heating for CO 2 desorption from zeolite 13X

In this study, we investigate microwave irradiation as an alternative to conventional heating for temperature swing adsorption processes. The performance of microwave and conventional heating during sorbent regeneration was evaluated by measuring CO 2 desorption from zeolite 13X at different temperatures. Experimentally, a fixed bed of zeolite 13X was saturated by a 150 sccm flow of 15 % CO 2 at room temperature followed by sorbent regeneration under nitrogen at 55, 100, or 150 °C by applying either microwave irradiation or conventional heating. Microwaves reduced sorbent regeneration times by at least half compared to regeneration by conventional heating. Under conventional regeneration, desorption curves were resolved into two peaks representing physisorbed CO 2 (mass diffusion limited) at low temperature and bicoordinated CO 2 (thermally limited) with increasing temperature. Under microwave regeneration, only one desorption peak was observed suggesting that CO 2 desorption was limited by mass diffusion through the porous zeolite 13X structure, rather than by temperature. Depending on microwave power, apparent activation energy of the microwave-assisted regeneration was 15.8–18.1 kJ/mol, compared to 41.5 kJ/mol for conventional regeneration. The reduction in apparent activation energy is mainly attributed to selective microwave heating of CO 2 adsorption sites (Na + sites) resulting in greater steady state temperatures of Na + relative to framework atoms, suggesting greater heating efficiency due to microwaves compared to conventional heat transfer. Due to rapid cycling and efficient heat transfer to CO 2 sites on zeolite 13X, microwave regeneration is found to increase adsorption/desorption cycling productivity and potentially reduce the energy penalty of temperature swing capture.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

Argon/UF6 plasma exhaust gas reconstitution experiments using preheated fluorine and on-line diagnostics

The feasibility of employing a flowing, high-temperature, pure fluorine/UF6 regeneration system to efficiently convert a large fraction of the effluent plasma exhaust back to pure UF6 was demonstrated. The custom built T.O.F. mass spectrometer sampling system permitted on-line measurements of the UF6 concentration at different locations in the exhaust system. Negligible amounts ( 100 ppm) of UF6 were detected in the axial bypass exhaust duct and the exhaust ducts downstream of the cryogenic trap system used to collect the UF6, thus verifying the overall system efficiency over a range of operating conditions. Use of a porous Monel duct as part of the exhaust duct system, including provision for injection of pure fluorine, provided a viable technique to eliminate uranium compound residue on the inside surface of the exhaust ducts. Typical uranium compound mass deposition per unit area of duct was 2 micron g/sq cm. This porous duct technique is directly applicable to future uranium compound transfer exhaust systems. Throughout these experiments, additional basic data on the corrosion aspects of hot, pressurized UF6/fluorine were also accumulated.

Roman, W. C.↗
NASA NTRS2003

Modeling Candle Flame Behavior In Variable Gravity

The burning of a candle, as typical non-propagating diffusion flame, has been used by a number of researchers to study the effects of electric fields on flame, spontaneous flame oscillation and flickering phenomena, and flame extinction. In normal gravity, the heat released from combustion creates buoyant convection that draws oxygen into the flame. The strength of the buoyant flow depends on the gravitational level and it is expected that the flame shape, size and candle burning rate will vary with gravity. Experimentally, there exist studies of candle burning in enhanced gravity (i.e. higher than normal earth gravity, g(sub e)), and in microgravity in drop towers and space-based facilities. There are, however, no reported experimental data on candle burning in partial gravity (g < g(sub e)). In a previous numerical model of the candle flame, buoyant forces were neglected. The treatment of momentum equation was simplified using a potential flow approximation. Although the predicted flame characteristics agreed well with the experimental results, the model cannot be extended to cases with buoyant flows. In addition, because of the use of potential flow, no-slip boundary condition is not satisfied on the wick surface. So there is some uncertainty on the accuracy of the predicted flow field. In the present modeling effort, the full Navier-Stokes momentum equations with body force term is included. This enables us to study the effect of gravity on candle flames (with zero gravity as the limiting case). In addition, we consider radiation effects in more detail by solving the radiation transfer equation. In the previous study, flame radiation is treated as a simple loss term in the energy equation. Emphasis of the present model is on the gas-phase processes. Therefore, the detailed heat and mass transfer phenomena inside the porous wick are not treated. Instead, it is assumed that a thin layer of liquid fuel coated the entire wick surface during the burning process. This is the limiting case that the mass transfer process in the wick is much faster than the evaporation process at the wick surface.

Alsairafi, A.↗
DOE OSTI2026

Geochemical evolution in Cacapon member: Fluid-rock interaction experiments and model insights for Appalachian Basin geothermal development

Here, this study combines recirculated flow-through experimental results conducted for 17 days at 90C and 200 PSI with reactive transport modeling to estimate fluid-rock interactions occurring in a sandy mudstone using an interbedded sandstone-shale sample from the Cacapon Member of the upper Tuscarora Sandstone/lower Rose Hill Formation for the purpose of geothermal exploration. Results suggest that the fluid and rock are likely to be in or near partial-equilibrium after approximately one year. In addition, after >400 h of continuous injection at 0.05 mL/min (5*10 −8 m 3 /min) the reactive front is restricted to the first ∼13% of the 4 cm*1.6 cm experimental rock length, whereas after >9100 h (∼1 year), the reactive front extends to 30% total length. The rate of changes in dissolution or precipitation are however, very minimal, with all major rock-forming minerals having rates <10 −11 mol/L porous media/s. Reservoir rock in the presence of dilute brine as may be the case during operation of an enhanced geothermal system would experience little alteration during the shut-in phase, and possibly up to one year. These results have utility in geothermal exploration for reservoirs at similar temperatures as well as general fluid-sandy mudstone rock interaction geochemistry.

Appalachian Basin↗
DOE OSTI2025

Quantifying the Effect of Pore‐Size Dependent Wettability on Relative Permeabilities Using Lattice Boltzmann Simulation

Relative permeability is a crucial two-phase property in porous media that can be significantly impacted by wettability conditions. While traditional research has predominantly examined homogeneous wettability, this work explores the less studied pore-size dependent (PSD) wettability, featured by a pore-size dependent wettability distribution. Leveraging high-fidelity Lattice Boltzmann simulations on CT-scanned porous samples, we demonstrate how PSD wettability would impact relative permeability at the pore scale. Our findings reveal that the deviation of relative permeability from the homogeneous wettabilitity induced by PSD wettability can be 5%–20%. The deviation of relative permeability curves increases as the spanning range of the contact angle increases. We also find that this impact is less pronounced as the capillary number increases. By adopting a pore-size-dependent contact angle relationship, our approach provides a more accurate and nuanced understanding of how PSD wettability would impact two-phase flow. These findings discovered at the pore scale may also provide valuable insights on relative permeability at the core to reservoir scales.

25 ENERGY STORAGE↗