Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “capillary flows”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An Experimental Investigation of Sintered Particle Effect on Heat Transfer Performance in an “Annular Flow” Evaporation Tube

Wicking structures have been widely used within passive heat transfer devices with high heat fluxes, such as heat pipes, to enhance their thermal performance. While wicking structures promote capillary pumping of the working fluid and thin film evaporation, they can result in capillary evaporation and further enhance the evaporation heat transfer. In this study, a 0.5 mm thick layer of 105 µm sintered copper particles was added to the inner wall of a copper tube, aiming to form an “annular flow” and enhance the heat transfer characteristics by taking advantage of thin film and capillary evaporation. Acetone was chosen as the working fluid, and the performance of an evaporation tube was tested for power inputs of 10, 30, 50, and 70 W. For each power input, trials were run at inclination angles varying from -90 deg to 90 deg to investigate the capillary effects. The temperature measurements showed that the temperature distribution along the evaporation tube is always downward sloping, meaning the temperature at the fluid inlet is larger than the outlet. In conclusion, results show that an “annular flow” formed by a thin layer of sintered particles can promote thin film and capillary evaporation and, therefore, boost the evaporation heat transfer coefficient.

42 ENGINEERING↗
DOE OSTI2024

Quantifying the Effect of Pore‐Size Dependent Wettability on Relative Permeability Using Capillary Bundle Model

Abstract Relative permeability is a key parameter for characterizing the multiphase flow dynamics in porous media at macroscopic scale while it can be significantly impacted by wettability. Recently, it has been reported in microfluidic experiments that wettability is dependent on the pore size (Van Rooijen et al., 2022). To investigate the effect of pore‐size‐dependent wettability on relative permeability, we propose a theoretical framework informed by digital core samples to quantify the deviation of relative permeability curves due to wettability change. We find that the significance of impact is highly dependent on two factors: (i) the function between contact angle and pore size (ii) overall pore size distribution. Under linear function, this impact can be significant for tight porous media with a maximum deviation of 1,000%.

Yu, Siqin↗
DOE OSTI2023

Capillary barriers for staged loading of microfluidic devices

Various aspects of the present disclosure are directed toward methods and apparatuses for interacting a first liquid and a second liquid in one or more fluidic channels of a capillary structure. The methods and apparatuses can include providing at least one capillary barrier that positions a meniscus of the first liquid at a fluid-interface region using capillary forces within the capillary structure. Additionally, a path is provided along one of the channels for the second liquid to flow toward the fluid-interface region. Additionally, gas pressure is released, via a gas-outflow port, from the fluid-interface region while flow of the first liquid is arrested. Further, the first liquid and the second liquid contact in the fluid-interface region with the capillary barrier holding the first liquid at the fluid-interface region.

Santiago, Juan G.↗
DOE OSTI2023

Balanced Fast-SpaciMS capillary configurations provide practically noninvasive channel-average measurements in catalytic monoliths

Spatially resolved capillary-inlet mass spectrometry (SpaciMS) provides a detailed picture of the spatiotemporal evolution of reaction network in catalytic monoliths. In the present work, we combine the SpaciMS experiments with a newly developed non-isothermal 3D CFD model for heterogeneously catalyzed reactive flows, including diffusion and permeation through the coated catalyst and channel wall. We explore how the capillary size and sampling rate can be balanced to minimize the impact on probed-channel conversions and provide species concentration measurements representative of the free channel average. In all studied configurations, the balanced sampling-rate fraction is noticeably higher than the corresponding capillary occlusion fraction. For a typical 350-micron capillary that occupies 11% of the channel cross- section, the balanced sampling rate represents 43% of the channel flow. Furthermore, the present work shows that together, a balanced SpaciMS configuration and numerical simulation provide accurate channel-averaged information of a monolithic catalyst under realistic operation conditions.

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

Liquid–Liquid Microfluidic Flows for Ultrafast 5-Hydroxymethyl Furfural Extraction

We study the effects of liquid–liquid flow patterns on the extraction of 5-hydroxymethyl furfural from the aqueous phase into an organic phase (methyl isobutyl ketone, ethyl acetate, and 2-pentanol) at ultrashort residence times of 0.2–20 s in a coiled capillary microchannel. Slug, droplet, slug-droplet, dispersed, parallel, annular, and irregular two-phase flow patterns are generated. We evaluate the extraction efficiency and the volumetric mass transfer coefficient via experiments and computational fluid dynamics (CFD). The CFD predictions are in quantitative agreement with the experimental data. The convective and diffusive contributions to the mass transfer of different flow patterns are analyzed. For all solvent pairs, the extraction efficiency varies irregularly with increasing flow rate due to changes in the mass transfer mechanism caused by a transition between flow regimes. Overall, we demonstrate >90% extraction efficiencies with high volumetric mass transfer coefficients ranging from 0.006 to 2.17 s–1. Furthermore, a merit index is defined to assess the trade-off between the fast mass transfer and the pressure drop in biphasic microchannels. In general, the mass transfer rate increases significantly at high flow rates at the expense of a high-pressure drop. Furthermore, the ethyl acetate/water irregular flow pattern is optimal due to a large throughput, extraction efficiency, and mass transfer coefficient, with a minimal pressure drop.

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

Experimental Investigation of Buoyant Flow in Realistic Bedforms With Heterogeneous Wettability

Submeter-scale geologic heterogeneity greatly affects CO 2 plume migration and retention. In this work, we present meter-scale laboratory experiments that can capture the impact of realistic submeter-scale geologic heterogeneity on multiphase flow and trapping. We produce realistic sedimentary formations consisting of ripple deposits with varying grain size contrast and wettability in a meter-scale slab chamber. Then, we conduct multiphase flow experiments with analog fluids through these structures and measure the saturation patterns, capillary heterogeneity trapping (CHT), and overall trapping performance. When we alter the ripple bedform architecture, variations in trapped saturation and CHT (10–20%) increment are exhibited. Similar growth in trapping performance is also observed when grain size contrast increases. Finally, wettability changes (water- to oil-wet) can increase nonwetting saturation and CHT up to 5% and 10–20%, respectively. These results emphasize the importance of correctly characterizing the impact of small-scale heterogeneities and wettability changes. We believe this is the first time that multiphase flow experiments were conducted in meter-scale domains with realistic ripple bedforms and heterogeneous wettability to investigate plume migration and trapping.

58 GEOSCIENCES↗
DOE OSTI2021

Heterogeneous multiphase flow properties of volcanic rocks and implications for noble gas transport from underground nuclear explosions

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

When ancient numerical demons meet physics-informed machine learning: adjoint-based gradients for implicit differentiable modeling

Recent advances in differentiable modeling, a genre of physics-informed machine learning that trains neural networks (NNs) together with process-based equations, have shown promise in enhancing hydrological models' accuracy, interpretability, and knowledge-discovery potential. Current differentiable models are efficient for NN-based parameter regionalization, but the simple explicit numerical schemes paired with sequential calculations (operator splitting) can incur numerical errors whose impacts on models' representation power and learned parameters are not clear. Implicit schemes, however, cannot rely on automatic differentiation to calculate gradients due to potential issues of gradient vanishing and memory demand. Here we propose a “discretize-then-optimize” adjoint method to enable differentiable implicit numerical schemes for the first time for large-scale hydrological modeling. The adjoint model demonstrates comprehensively improved performance, with Kling–Gupta efficiency coefficients, peak-flow and low-flow metrics, and evapotranspiration that moderately surpass the already-competitive explicit model. Therefore, the previous sequential-calculation approach had a detrimental impact on the model's ability to represent hydrological dynamics. Furthermore, with a structural update that describes capillary rise, the adjoint model can better describe baseflow in arid regions and also produce low flows that outperform even pure machine learning methods such as long short-term memory networks. The adjoint model rectified some parameter distortions but did not alter spatial parameter distributions, demonstrating the robustness of regionalized parameterization. Despite higher computational expenses and modest improvements, the adjoint model's success removes the barrier for complex implicit schemes to enrich differentiable modeling in hydrology.

58 GEOSCIENCES↗
DOE OSTI2025

Three-Dimensional Bubble Fluidics in Architected Porous Media

Gas bubble flows in porous media often exhibit complex and seemingly unpredictable behaviors that are difficult to control. This lack of control limits the ability to design effective devices which manage multiphase flows. Here, we show how the design of 3D printed pores can deterministically control the flow path of an injected gas stream. Open cell structures can be designed to shape the gas/liquid interface with fidelity to control how the two phases are distributed throughout a porous material. The distributed gas volume is free to interact physically and chemically with the surrounding liquid phase, an effect we exploit to create a logical control gate to redirect flows within a lattice. This also allows us to design architectures for reactive capture and aerating bioreactors, resulting in patterned boundaries which can make more effective use of the liquid and gas reagents.

3D microfluidics↗
DOE OSTI2021

Cellular fluidics

The natural world provides many examples of multiphase transport and reaction processes that have been optimized by evolution. These phenomena take place at multiple length and time scales and typically include gas–liquid–solid interfaces and capillary phenomena in porous media. Many biological and living systems have evolved to optimize fluidic transport. However, living things are exceptionally complex and very difficult to replicate, and human-made microfluidic devices (which are typically planar and enclosed) are highly limited for multiphase process engineering. In this paper, we introduce the concept of cellular fluidics: a platform of unit-cell-based, three-dimensional structures—enabled by emerging 3D printing methods—for the deterministic control of multiphase flow, transport and reaction processes. We show that flow in these structures can be ‘programmed’ through architected design of cell type, size and relative density. We demonstrate gas–liquid transport processes such as transpiration and absorption, using evaporative cooling and CO 2 capture as examples. We design and demonstrate preferential liquid and gas transport pathways in three-dimensional cellular fluidic devices with capillary-driven and actively pumped liquid flow, and present examples of selective metallization of pre-programmed patterns. Our results show that the design and fabrication of architected cellular materials, coupled with analytical and numerical predictions of steady-state and dynamic behaviour of multiphase interfaces, provide deterministic control of fluidic transport in three dimensions. Cellular fluidics may transform the design space for spatial and temporal control of multiphase transport and reaction processes.

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

The effect of inhomogeneities in coal-based mesophase pitch on filament formation and melt spinning for carbon fiber

Mesophase pitch-based carbon fibers excel as high specific strength and modulus materials, but challenges in stable, continuous melt spinning processing are a serious limiting factor. Inhomogeneities in the pitch may result in flow instabilities that can lead to filament breakage. This work aims to characterize and offer mitigation strategies for the negative effects that inhomogeneities such as isotropic binder, gels, and other inclusions impart on melt spinning of mesophase pitch. The effect of inhomogeneities on spinning stability are analyzed using capillary rheometry, borescope imaging of elongating filaments, and filtration. With no filter, immediate blockage of the spinning capillary occurred. Improvement in melt spinning duration was observed when a porous metal filter (20 μm nominal pore size) was utilized, but no significant further improvement in spinning duration was found upon decreasing the filter nominal pores size to 10, 5 and 2 μm. Moreover, flow instabilities that disrupted melt spinning stability were observed by capillary rheometry when the size of included isotropic pitch domains were larger than approximately 10 μm. Resultant carbon fiber properties and SEM imaging of fracture surfaces are reported.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Effect of topology changes on the breakup of a periodic liquid jet

Here, the breakup of a periodic jet is examined computationally, using a front-tracking/finite-volume method, where the interface is represented by connected marker points moving with the fluid, while the governing equations are solved on a fixed grid. Tracking the interface allows control of whether topology changes take place or not. The Reynolds and Capillary numbers are kept relatively low ($Re = 150$ and $Ca = 2$) so most of the flow is well resolved. The effect of topology changes is examined by following the jet until it has mostly disintegrated, for different “coalescence criterion,” based on the thickness of thin films and threads. The evolution of both two-dimensional and fully three-dimensional flows is examined. It is found that although there is a significant difference between the evolution when no breakup takes place and when it does, once breakup takes place the evolution is relatively insensitive to exactly how it is triggered for a range of coalescence criterion, and any differences are mostly confined to the smallest scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Topology-Dependent Polymer Stretching and Scission in Solution at Extreme Shear Rates

There has been significant interest in the engineering of polymer topology to control rheology and mechanical stability in dilute solutions for applications involving extreme shear rate flows. However, methods to experimentally probe properties at relevant shear rates (≈ 10 4 –10 6 s –1 ) have remained ex situ, obscuring access to measures of polymer deformation and rheology that would otherwise provide mechanistic insight into the topology-dependent properties that control their behavior in extreme shear flows. In this study, we used novel in situ small angle neutron scattering measurements in a capillary rheometer (capillary rheo-SANS) to simultaneously measure solution viscosities and polymer deformations in high shear on a series of chemically homologous topology-defined polymers including linear, randomly branched, and star-shaped molecules. We demonstrate that differences in the onset of chain stretching and shear thinning of these polymers in dilute solution are controlled primarily by differences in their molecular relaxation time. These differences correlate with differences in chain scission inferred from ex situ measurements at more extreme shear rates. Together, the results demonstrate a direct coupling between chain deformation and scission, and suggest that the dominant effect of branching as a means to impart resilience against mechanical degradation is through differences in relaxation dynamics due to branching. We anticipate that these results will provide key insights to engineer topology-controlled polymers for rheological modification, mechanical stability, and controlled mechano-chemistry.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Using Direct Numerical Simulation of Pore-Level Events to Improve Pore-Network Models for Prediction of Residual Trapping of CO2

Direct numerical simulation and pore-network modeling are common approaches to study the physics of two-phase flow through natural rocks. For assessment of the long-term performance of geological sequestration of CO 2 , it is important to model the full drainage-imbibition cycle to provide an accurate estimate of the trapped CO 2 . While direct numerical simulation using pore geometry from micro-CT rock images accurately models two-phase flow physics, it is computationally prohibitive for large rock volumes. On the other hand, pore-network modeling on networks extracted from micro-CT rock images is computationally efficient but utilizes simplified physics in idealized geometric pore elements. This study uses the lattice-Boltzmann method for direct numerical simulation of CO 2 -brine flow in idealized pore elements to develop a new set of pore-level flow models for the pore-body filling and snap-off events in pore-network modeling of imbibition. Lattice-Boltzmann simulations are conducted on typical idealized pore-network configurations, and the interface evolution and local capillary pressure are evaluated to develop modified equations of local threshold capillary pressure of pore elements as a function of shape factor and other geometrical parameters. The modified equations are then incorporated into a quasi-static pore-network flow solver. The modified model is applied on extracted pore-network of sandstone samples, and saturation of residual trapped CO 2 is computed for a drainage-imbibition cycle. The modified model yields different statistics of pore-level events compared with the original model; in particular, the occurrence of snap-off in pore-throats is reduced resulting in a more frontal displacement pattern along the main injection direction. Compared to the original model, the modified model is in closer agreement with the residual trapped CO 2 obtained from core flow experiments and direct numerical simulation.

Kohanpur, Amir H.↗
DOE OSTI2023

Effects of grain size and small-scale bedform architecture on CO 2 saturation from buoyancy-driven flow

Small-scale (mm-dm scale) heterogeneity has been shown to significantly impact CO 2 migration and trapping. To investigate how and why different aspects of small-scale heterogeneity affect the amount of capillary trapping during buoyancy-driven upward migration of CO 2 , we conducted modified invasion percolation simulations on heterogeneous domains. Realistic simulation domains are constructed by varying two important aspects of small-scale geologic heterogeneity: sedimentary bedform architecture and grain size contrast between the matrix and the laminae facies. Buoyancy-driven flow simulation runs cover 59 bedform architecture and 40 grain size contrast cases. Simulation results show that the domain effective CO 2 saturation is strongly affected by both grain size and bedform architecture. At high grain size contrasts, bedforms with continuous ripple lamination at the cm scale tend to retain higher CO 2 saturation than bedforms with discontinuous or cross lamination. In addition, the “extremely well sorted” grain sorting cases tend to have lower CO 2 saturation than expected for cross-laminated domains. Finally, both a denser CO 2 phase and greater interfacial tension increase CO 2 saturation. Again, variation in fluid properties seems to have a greater effect on CO 2 saturation for cross-laminated domains. This result suggests that differences in bedform architecture can impact how CO 2 saturation values respond to other variables such as grain sorting and fluid properties.

58 GEOSCIENCES↗
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↗
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 OSTI2024

Capillary-Enhanced Two-Phase Micro-Cooler Using Copper-Inverse-Opal Wick with Silicon Microchannel Manifold for High-Heat-Flux Cooling Application

In this work, we demonstrate a two-phase capillary-fed boiling micro-cooler that consists of a ~ 25-..mu..m-thick copper inverse opal (CIO) porous wicking structure for high-heat-flux boiling and a silicon 3D-manifold for distributed liquid delivery and vapor extraction across a 0.5 cm x 0.5 cm heated area. At low inlet water mass flow rates of 1.5 to 1.9 g(min)-1, the micro-cooler displays nearly two-phase boiling with exit vapor quality ~ 1 and a high critical heat flux (CHF) of 253 to 320 W cm-2 with low superheat of ~ 10 degrees C resulting in a thermal resistance of boiling ~ 0.025 cm2 degrees C W-1 or heat transfer coefficient of 0.4 MW m-2 degrees C-1. For higher flow rates of 5, 10, and 15 g(min)-1, the micro-cooler exhibits a hybrid single-phase and two-phase cooling regime where the contribution of the sensible heat (single-phase) cooling is linearly added to that of the two-phase cooling. For the highest flow rate of 15 g(min)-1, the CHF is increased to ~ 500 W cm-2 resulting in an overall thermal resistance of ~ 0.18 cm2 degrees C W-1. However, the two-phase heat transfer effectiveness, which estimates the utilization level of the inlet mass flow rate for two-phase boiling, is reduced to ~ 0.11. To achieve the best cooling system performances, the micro-cooler must operate entirely within the two-phase boiling regime (exit vapor quality or two-phase heat transfer effectiveness ~ 1). Ideally, the "coolant" should be delivered near its saturation temperature (~ 100 degrees C for water), which provides significant advantages for the energy-efficient operation of data centers and power electronics. We present detailed analysis with Infrared and high speed camera images at various inlet flow rates and heat fluxes to understand complex heat transfer in the micro-cooler. Furthermore, a conjugate thermofluidic simulation model, which incorporates the physics of capillary-fed boiling in a porous copper layer, agrees well with the experimental data.

capillary-enhanced boiling↗