Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two phase 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Evaporator Temperature Transient Testing of a High-Performance Sodium Filled Heat Pipe

Heat pipes are two-phase heat transfer devices that enable passive removal of heat from the reactor core to the power conversion system in heat pipe-cooled microreactor designs. Experimental investigations of heat pipe transients are needed for technology demonstration, verification and validation of numerical codes, and the establishment of regulatory requirements. The Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL) serves as a platform for evaluating the dynamic response of high-temperature heat pipes under a variety of operating conditions. The present work details the experimental investigation of a high-performance, defined as over 2 kW sodium heat pipe subjected to rapid input power fluctuations induced by sudden changes in the evaporator temperature setpoint. In addition, the heat pipe was subjected to an asymmetrical heat load where a subset of heaters operated at 30% and 70% below their nominal power. These experimental conditions were chosen to simulate thermal and operational stresses expected to be encountered in microreactors to provide data on heat pipe behavior during such important transient events. Key data and performance metrics, including time series of temperatures and strains, axial temperature profiles, thermal response times, and heat transfer capabilities, the thermal output over thermal input, were reported and discussed. The results highlight the resilience of heat pipes, revealing their potential to maintain thermal stability and efficiency under varying power loads. Lastly, the paper concludes with a discussion on the significance of the results and their implications for future research.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

FY24 Report on Water NSTF Testing: Lower Tank Inlet Piping Configuration

The following report serves as a summary of the accomplishments and testing results by the Natural convection Shutdown heat removal Test Facility (NSTF) experimental program over the past 12-month period. A major activity included reconfiguration of the chimney piping geometry which altered the discharge position into the tank from the original 50% tank height to a new lower position at a 10.9% tank height. This modification increases the volume of available coolant thus extending long-term operating capacity, however also results in a decreased liquid driving head which has the potential to reduce the natural circulation efficiency and decrease overall performance. Examination of the tradeoffs for these two configurations is an area of interest for RCCS designers and drove planned test activities this year. Ten matrix tests were performed in FY24, totaling 195 hours of heated operations and 12.9 MWh of electrical heating throughout the year, with eight classified as Accepted per NQA-1 and two classified as Trending. Testing prior to the facility reconfiguration examined the facility response to throttling at the tank inlet, demonstrating increased sensitivity of flow instabilities to throttling within the two-phase region when compared to throttling at the lower sensitivity single-phase inlet region. Testing after the chimney reconfiguration began with a baseline test at conditions of 80% inventory fill and prototypic thermal power input of 2.1 MW t . In addition to establishing a reference for nominal system behavior and performance, repeat testing with multiple subsequent runs demonstrated strong repeatability between tests performed at the same conditions in this new configuration. These tests also examined if influences would occur to system behavior with installation of new higher-resolution instrumentation within the upper chimney. In this critical region where boiling and flashing phenomena dominate, the test results provided confidence that the new instrumentation does not uncharacteristically influence the observed behaviors. An inventory parametric series was then initiated to examine system behavior with the new chimney configuration at six varying initial inventory levels, ranging from high 80% to low 20% fill. Two-phase oscillations with similar peak and mean flow rates were observed when comparing to the previous mid tank configuration. Generally, similar system response trends with inventory were also observed, such as two phase oscillations suppressing as inventories were lowered. However, in one absolute inventory comparison at the highest fill of 80%, flow oscillations saw a gradual growth over the 4 hours of two-phase operation in the new lower tank inlet configuration, opposite of the gradual dampening observed in the mid tank inlet configuration. This can be attributed to, in part, a greater hydrostatic head pressure above the two-phase discharge region where the boiling front is developed. Furthermore, the change resulted in greater sensitivity to liquid degassing phenomena during single-phase heating, causing loop instabilities to form which trigged moderate flow degradation during the period approaching saturation and boiling conditions. This behavior had been observed previously under some conditions but was a common occurrence in recent testing with the newer lower tank inlet configuration. Initial observations from these first data sets suggest an overall larger window of stability for the mid tank inlet configuration compared to the lower tank inlet. Lastly, the lowest inventory fill test was repeated over an extended testing window to examine depletion behavior. Natural circulation flow and effective heat removal performance were observed during most of the testing period; only after the tank became fully drained (0% fill) did flow stagnate and violent geysering events occur. This early observation confirms one relative advantage over the mid tank inlet configuration, which stagnated under comparable conditions at ~20% inventory remaining in the tank.

42 ENGINEERING↗
DOE OSTI2024

FY-24 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the computational modeling of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY24. This year’s modeling campaign focuses on the continual testing of the RELAP5 model against experimental data. Several fault cases were simulated with the RELAP5 model in FY23 to study the accuracy of the model under complex flow conditions. Similarly, in FY24, a case is simulated where the chimney outlet is throttled with a valve while the flow is undergoing two-phase oscillations. The purpose of this case is to study the effects of increasing pressure drop in the two-phase region on the system behavior. Even though the model is able to capture the experimental data qualitatively, it overpredicts the pressure drops experienced by the flow in the two-phase regime. Nevertheless, the model and the experiment show that oscillations are stabilized with increased pressure drop in the two-phase region.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Boiling-Water Reactor Testing Capability in the Advanced Test Reactor

I-Loop is an irradiation facility that is currently being installed at the Advanced Test Reactor. It is a two-loop test facility capable of performing Light Water Reactor (LWR) irradiations in prototypic coolant conditions. The two loops are being installed to be capable of both Boiling Water Reactor (BWR) and Pressurized Water Reactor (PWR) pressure, temperature, and chemistry environments. Each loop is nominally dedicated as a BWR or PWR for simplicity of operations. In-reactor water loop testing that an I-Loop provides is key to the deployment of new accident tolerant fuel technologies and other advanced LWR fuel concepts. Currently, pressurized water loops are the only testing facilities available to test BWR fuel concepts. Their test environments are non-prototypic at higher pressure/temperature and at single-phase fluid flow conditions. This void in the LWR test bed capabilities is one that the I-Loop is uniquely situated to provide. This report discusses the mechanical design, thermal hydraulic calculations, and neutronic calculations of a proposed standard experiment of accident tolerant BWR fuel concepts. Mechanical design examines the geometry and features of the main components. Thermal hydraulic calculations examine the modeling and results of the two-phase flow options available. Lastly, neutronic calculations examine the Monte Carlo analysis of enrichment, heat rates, and flux spectrum.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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↗
DOE OSTI2025

Validation of Pronghorn for Natural-Circulation Molten Salt Loops

This paper presents the development and validation of a high-fidelity thermal-hydraulic model of a molten salt natural circulation flow loop, designed for integration within a digital twin framework. The study evaluates the performance of Idaho National Laboratory’s Pronghorn against experimental data from Texas A&M University Molten Salt Flow Loop (MSFL) four Hitec-salt test benchmark data. Natural circulation of high-Prandtl-number fluids exhibits complex, counter-intuitive flow patterns that make pointwise thermocouple readings unreliable. Experimental work at TAMU’s MSFL provides benchmark data, including flow visualization at a test-section and centerline steady-state temperature measurements along the loop. Validation includes four single-phase natural circulation test cases with Hitec salt. Key metrics include flow profile agreement and steady-state temperature accuracy. Pronghorn results for two-dimensional single-phase agree qualitatively with the experimental flow profile. This paper illustrates the importance of Computational Fluid Dynamics (CFD) in elucidating the behavior of high-Prandtl-number thermal-hydraulics, along with how misleading centerline temperature measurements can be. Pronghorn reproduces the axial and radial stratification that makes single thermocouple readings unreliable. Future research will focus on reduced-order modeling techniques to enable rapid simulation suitable for real-time digital twin applications. The validated cases provide a basis for developing reduced-order surrogates aimed at real-time digital-twin applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Flow maldistribution in plate heat exchangers – Impact, analysis, and solutions

Plate heat exchangers are commonly used in various industrial applications, such as refrigeration, air conditioning, heat pumps, powerplants, and chemical industries. Plate heat exchangers are well known for their superior heat transfer performance, compactness, and low refrigerant charge. Despite offering several advantages, they suffer from flow maldistribution issues. The flow maldistribution can deteriorate both the heat transfer and pressure drop performance, ultimately resulting in a lower system efficiency where plate heat exchangers are deployed. The flow maldistribution issues become more pronounced when the heat exchanger size is relatively large and the number of plates is higher, limiting the deployment of plate heat exchangers in larger industrial systems. Consequently, analyzing and understanding the flow maldistribution behavior in plate heat exchangers and finding ways to mitigate flow maldistribution related issues become essential topics of interest. Further, this review aims to address the effect of flow maldistribution on plate heat exchanger characteristics. First, the experimental and numerical works on flow maldistribution under single-phase and two-phase conditions are detailed. Subsequently, the end-channel and end-plate effects are discussed. Then, the methods to mitigate flow maldistribution in plate heat exchangers are outlined. Finally, based on a thorough literature survey and industrial requirements, future research directions are recommended.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

zmodel-kokkos

The code implements the so-called ZModel on an unstructured mesh. This models two-phase (or feasibly multi-phase) incompressible flow by meshing and evolving only the interface between materials, and never directly discretizing the bulk flow.

May, Ian↗
DOE OSTI2025

Unraveling the Determinant Mechanisms in Flow-Mediated Crystal Growth and Phase Behaviors

To uncover the critical mechanisms responsible for mesoscopic level development during flow-mediated crystal growth, we develop a semi-two-way hydrodynamic coupled structural phase-field crystal formalism (HXPFC-s2). The new formalism, inspired by previous attempts at coupling hydrodynamic and phase-field crystal (PFC) equations, allows for studying mesoscopic flow-mediated crystallization at diffusive timescales pertinent to industrial applications. Unlike previous efforts, the devised coupling to the structural PFC (XPFC) equations allows generalization to more complex crystal structures through explicit parameterization of the direct correlation function (DCF). Utilizing the HXPFC-s2 formalism, we seek to uncover the determinant physical mechanisms in crystallization under simple shear flows by comparing temperature-driven crystallization to flow-mediated crystallization under varying flow-strengths. Parallels and deviations of under-cooling and flow-strength effects on crystal growth are drawn using the crystal cluster-size and system ordering time evolutions. In doing so, we identify scaling behaviors with a Peclet-like number, Pe∼, a critical Peclet-like number, Pe∼*, and flow-field-crystal plane-dependent interactions. Our findings may be relevant for controlling crystal growth and phase behaviors in flow applications.

Willis, L. Connor (ORCID:0009000961321848)↗
DOE OSTI2025

Experimental Validation of a 40kW, 480V Point-to-Point DC Interlinks for Controller-Agnostic, Interoperable Networked Microgrids

This paper presents the experimental validation of point-to-point dc-interlinks for interconnecting two solar-based, laboratory-scale AC microgrids. DC interlinks provide a solution to numerous technical and operational challenges encountered in networked microgrids, including precise power flow control, stable and fast synchronization, enhanced stability, and improved voltage and frequency regulation. By decoupling microgrids through power converters, dc-interlinks enable power exchange among microgrids that can be owned by different entities (such as communities, utilities, or universities) and managed by diverse microgrid controller vendors. This characteristic makes this dc-interlinks a promising solution for networking real-world microgrids. The presented point-to-point dc interlink utilizes two four-quadrant converters: a 40kW 3-phase ac/dc regulating the dc-link voltage to 800V, and a 3-phase dc/ac controlling power flow. These converters, connected via a 20-foot dc cable, interconnect two ac microgrids operating at 480V, each featuring energy storage, photovoltaic generation, and load emulation. Experimental validation employs commercially available off-the-shelf (COTS) converters and real-world data from solar-powered microgrids in Adjuntas, Puerto Rico. To the authors’ knowledge, this work provides the first at-scale experimental validation of dc-interlinks for networked ac microgrids using COTS inverters, demonstrating their practicality and effectiveness in addressing real-world operational challenges.

Ferrari Maglia, Max [ORNL]↗
DOE OSTI2024

Flexible dynamic boundary microgrid operation considering network and load unbalances

Flexible microgrids with dynamic boundaries have recently been introduced in the literature. With the ability to reconfigure the topology of the microgrids dynamically through remotely controlled switches, flexible microgrids with dynamic boundaries can further improve the resiliency and energy efficiency of microgrids with distributed energy resources (DERs). This paper focuses on the optimal operation considering one of the predominant characteristics of microgrids and distribution systems – unbalanced networks and loads. In existing literature, balanced modeling of microgrids is more common due to its attractive simplicity. The three-phase power unbalance has not been considered as a constraint on the generation units in a microgrid. Further, negative sequence constraints have also been neglected. In this article, we propose a set of constraints that is specifically related to the capabilities of inverter interfaced resources to supply unbalanced current/power when the microgrid is islanded from the main distribution grid. We incorporate the new set of constraints into two optimization formulations leveraging two convex relaxations of the three-phase power flow equations: mixed-integer linear programming (MILP) and mixed-integer semidefinite programming (MISDP) that optimize the dispatch of controllable switches and DERs in the microgrid. The algorithms are then extended to networked microgrids with grid-forming sources. We test the algorithms on a realistic community microgrid model in Puerto Rico as well as standardized IEEE distribution test feeders. The testing results demonstrate the performance of the proposed algorithms. The MILP is fast and scalable, and the MISDP enforces the negative sequence voltage constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Microstructural, Oxidation, and Mechanical Behavior of NbTi-Based Refractory Alloys with 5 to 10 Pct Co, Cr, and Ni Additions

NbTi-based refractory alloys with additions of Co, Cr, and Ni represent an interesting medium-entropy alloy system with potential for protective oxide film formation, high strength, and ductility. This study investigates the microstructural evolution, oxidation behavior, and mechanical properties of NbTi-based alloys containing 5 to 10 at. pct Co, Cr, and Ni. CALPHAD predictions suggest that this composition range can be heat treated to obtain a predominantly body-centered cubic matrix phase. Mechanical properties, including microhardness, yield strength, maximum strength, and specific strength are evaluated through isothermal compression tests conducted between room temperature and 800 °C. The oxidation kinetics of these alloys are assessed through discontinuous oxidation tests. Parabolic oxidation kinetics were observed for NbTi–10Ni and NbTi–5Co, while linear oxidation kinetics were found for NbTi–10Cr and NbTi–10(CoCrNi). Microstructures and oxide layers are characterized using X-ray diffraction, electron backscatter diffraction, energy-dispersive X-ray spectroscopy, and scanning electron microscopy. All alloys exhibit significant mechanical softening between room temperature and 800 °C, with elastic-perfectly plastic flow observed at 800 °C. The addition of 10 pct Cr to NbTi resulted in two BCC phases up to 1050 °C, conflicting with CALPHAD predictions of a single-phase solid solution at this temperature, and resulting in higher flow stress at 800 °C. NbTi–10(CoCrNi) exhibited the lowest flow stress at 800 °C despite having more ‘cocktail effect’ potential and insufficient molar fractions of Co, Cr, or Ni to form a desirable protective oxide film.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GalaxyFlow: upsampling hydrodynamical simulations for realistic mock stellar catalogues

ABSTRACT Cosmological N-body simulations of galaxies operate at the level of ‘star particles’ with a mass resolution on the scale of thousands of solar masses. Turning these simulations into stellar mock catalogues requires ‘upsampling’ the star particles into individual stars following the same phase-space density. In this paper, we introduce two new upsampling methods. First, we describe GalaxyFlow, a sophisticated upsampling method that utilizes normalizing flows to both estimate the stellar phase-space density and sample from it. Secondly, we improve on existing upsamplers based on adaptive kernel density estimation (KDE), using maximum likelihood estimation to fine-tune the bandwidth for such algorithms in a way that improves both the density estimation accuracy and upsampling results. We demonstrate our upsampling techniques on a neighbourhood of the Solar location in two simulated galaxies: Auriga 6 and h277. Both yield smooth stellar distributions that closely resemble the stellar densities seen in the Gaia DR3 catalogue. Furthermore, we introduce a novel multimodel classifier test to compare the accuracy of different upsampling methods quantitatively. This test confirms that GalaxyFlow more accurately estimates the density of the underlying star particles than methods based on KDE, at the cost of being more computationally intensive.

Lim, Sung Hak (ORCID:0000000330981092)↗
DOE OSTI2026

Experimental Investigation of Uranium and Iron Condensation from High-Temperature Plasma Conditions

We used a plasma flow reactor (PFR) to generate synthetic fallout nanoparticles from vapor-phase condensation of two different input concentrations of uranium and iron analytes (U/Fe = 1:1 and 1:2). Synthetic fallout from complex chemical matrices (e.g., mixtures of U and Fe) has not been generated in this setup before and allows for the observation of relative condensation and fractionation of nuclear debris. Experiments were conducted under two different temperature histories with variations in particle flow patterns along the PFR. Transmission electron microscopy (TEM) observation and analysis of the nanoparticles revealed variations in speciation of uranium oxides (UO 2 and α-UO 3 ) depending on the competition between flow mixing and oxygen sequestration by iron. A ternary metal oxide, UFeO 4 , was observed in addition to iron oxides (e.g., FeO and Fe 3 O 4 ), which suggests that fallout models should account for chemical speciation of ternary metal oxides (i.e., UFeO 4 ) and their relative condensation behaviors in addition to those of singular metal oxides (e.g., FeO and UO 2 ). X-ray Energy Dispersive Spectroscopy (EDS) elemental maps showed that some particles had Fe-rich cores surrounded by U-rich regions. This suggests that either condensed U oxides coagulate onto molten Fe oxides or that the increase in iron analyte concentration might drive the system toward a higher degree of supersaturation, leading to earlier formation of iron oxide particles and providing an energetically favored pathway for nucleation of uranium oxides around iron oxide particles.

and nuclear chemistry↗
DOE OSTI2025

Benchmark study of the DTU OWC chamber with both two-way and one-way absorption

This paper reports on a benchmark study based on small-scale (1:50) measurements of a single, oscillating water column chamber mounted sideways in a long flume. The geometry of the OWC chamber is extracted from a barge-like, attenuator-type floating concept “KNSwing” with 40 chambers targeted for deployment in the Danish part of the North Sea. In addition to traditional two-way energy extraction we also consider one-way energy extraction with passive venting and compare chamber response, pressures and total absorbed energy between the two methods. A blind study was established for the numerical modeling, with participants applying several implementations of weakly nonlinear potential flow theory and commercial Navier–Stokes solvers (CFD). Both compressible and incompressible models were used for the air phase. Potential flow calculations predict more energy absorption near the chamber resonance for one-way absorption than for two-way absorption, but the opposite is found from the experimental measurements. This outcome is mainly attributed to energy losses in the experimental passive valve system, but this conclusion must be confirmed by better experimental measurements. Modeling the one-way valve in CFD proved to be very challenging and only one team was able to provide results which were generally closer to the experiments. The study illustrates the challenges associated with both numerical and experimental analysis of OWC chambers. Air compressibility effects were not found to be important at this scale, even with the large volume of additional air used for the one-way case.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗
DOE OSTI2024

Correlation Between Microscopic Current Fluctuations Observed at Ultra-Microelectrodes and Macroscopic Bulk Electrolysis Performance in Redox-Active Microemulsions

Microemulsions (μEs) have been proposed as redox flow battery (RFB) electrolytes that maximize ionic conductivity and charge capacity by synergizing two immiscible phases. However, charge transfer during electrolysis in μEs is poorly understood. Here, we show that ultramicroelectrode electrolysis of ferrocene-loaded μEs –20%, 60%, and 90% water - reveals stochastic current fluctuations. These are differentiated in the scanning electrochemical microscopy (SECM) geometry, where power spectral density analysis showed distinct changes in the frequency contributions. SECM in the substrate generation-tip collection mode showed that fluctuations arise under mass-transfer control. Significant differences in the diffusion coefficient of ferrocene species were deducted from SECM approach curves, suggesting phase transfer behavior. Using bulk electrolysis, we calculated the charge accessibility and cycling behavior in the μEs. A decrease in the stochastic behavior of the μEs seems to correlate to a higher accessibility and cycling performance, with the 90% water μE displaying the best reversibility and the 60% the lowest. Altogether, these results suggest that Marangoni-type convection driven by concentration gradients and/or μE restructuring during charge transfer play a role in the electrochemical performance of μEs. This presents opportunities for screening and diagnosing the performance of these emerging RFB electrolytes.

25 ENERGY STORAGE↗
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↗