Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid flow regimes”

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

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2023

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
DOE OSTI2024

Prandtl number effects on extreme mixing events in forced stratified turbulence

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ( $U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$ , where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ( $Fr=0.5, 2$ ) and Prandtl ( $Pr=1, 7$ ) numbers forced so that $Re_{b}=50$ , with resolutions up to $30\,240 \times 30\,240 \times 3780$ . We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$ ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$ .

Mechanics↗
DOE OSTI2020

Viscometric flow of dense granular materials under controlled pressure and shear stress

This study examines the flow of dense granular materials under external shear stress and pressure using discrete element method simulations. In this method, the material is allowed to strain along all periodic directions and adapt its solid volume fraction in response to an imbalance between the internal state of stress and the external applied stress. By systematically varying the external shear stress and pressure, the steady rheological response is simulated for: (1) rate-independent quasi-static flow; and (2) rate-dependent inertial flow. The simulated flow is viscometric with non-negligible first and second normal stress differences. While both normal stress differences are negative in inertial flows, the first normal stress difference switches from negative to slightly positive, and second normal stress difference tends to zero in quasi-static flows. The first normal stress difference emerges from a lack of coaxiality between a second-rank contact fabric tensor and strain rate tensor in the flow plane, while the second normal stress difference is linked to an excess of contacts in the shear plane compared with the vorticity direction. A general rheological model of second order (in terms of strain rate tensor) is proposed to describe the two types of flow, and the model is calibrated for various values of interparticle friction from simulations on nearly monodisperse spheres. The model incorporates normal stress differences in both regimes of flow and provides a complete viscometric description of steady dense granular flows.

42 ENGINEERING↗
DOE OSTI2020

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗
DOE OSTI2021

Influence of Heave Plate Topology on Reaction Force

In this study, we explore the role of heave plate topology on fluid reaction forces using three, representative shapes: a hexagonal flat plate, a hexagonal conic with an open top, and the same with a closed top that encloses a fluid mass. We force each test article sinusoidally in a quiescent tank and decompose the reaction force using forms of the Morison equation for phase-invariant and phase-dependent parameterizations. We find that a flat plate generates 5.3 % more fluid reaction force than the open conic topology, and 21.4% more than the enclosed conic. Asymmetric topologies generate asymmetric reaction forces, but the magnitude of asymmetry is limited by nearly symmetric fluid inertia forces, which dominate over drag for these test articles. Additionally, we observe asymmetric vortex dynamics for the flat plate when the Keulegan-Carpenter number is between 1 and 2, accompanied by a shift in the phase of the peak force by of the oscillation period. As a consequence of this shift, the hydrodynamic coefficients estimated from the phase-dependent Morison equation decomposition are asymmetric, suggesting that phase-dependent representations may not provide physical insight in some hydrodynamic regimes.

added mass↗
DOE OSTI2023

Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows

Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.

17 WIND ENERGY↗
DOE OSTI2025

Lithium Divertor Targets and Walls for the ASTER Liquid Stellarator Reactor, Distributed Divertor

Stellarators may have advantages for certain liquid metal options as Plasma Facing Components (PFC) for divertor targets and first walls due to the wide range of possible magnetic configurations, which additionally are free of disruptions and fast field variations. In a previous work (V. Queral et al., IEEE Trans. Plasma Sci. 52, 2024), a concept of stellarator reactor (ASTER-CP) based on swirling Li-molten salts and liquid lithium floating on the molten salt as PFC was presented. The divertor matters were not studied then and, thus, they are being studied and experimentally tested now. The ASTER-CP reactor concept, the initial liquid metal experiments and potential concepts for the ASTER-CP divertor and first wall are reported. Concerning the experiments, several small scale experiments of galinstan in a small rotating cylinder under magnetic field have been produced, including one experiment with high viscosity galinstan-mixture for increased thickness of layer. An experiment of floating lithium on the molten salt LiCl-PbCl2 gave fast volatilization/decomposition of the molten salt. Particularly for divertors, the traditional free-flow, Capillary Porous Systems and ‘divertorlets’ have been studied for application to ASTER-CP. Surface waves (hot spots), lack of enough surface fluid turbulence and excessive fluid speed are the main issues found in fast free-flow. The perhaps original concept of Distributed Divertor and Equi-power Surface is tentatively proposed and studied, taking advantage of stellarator fields and low recycling regime.

distributed divertor↗
DOE OSTI2025

Dimensionless Regime Mapping and Chaos Analysis of Solids Volume Fraction Data in a Circulating Fluidized Bed Riser

Presentation of experimental results of solids volume fraction data measurements in a 12-inch diameter cold flow circulating fluidized bed in relation to a dimensionless fluidization regime map. Statistical and chaotic system analysis results are shown to illustrate differences in hydrodynamic behavior at different operating conditions. Presentation was made at the ASTFE 10th Thermal and Fluids Engineering Conference, March 09-12, 2025.

chaos↗
DOE OSTI2026

Gas-driven short disconnection mitigates thermal runaway in Li-ion batteries under mechanical abuse

Mechanical abuse poses a critical safety risk to Li-ion batteries by inducing internal short circuits that initiate thermal runaway. Remarkably, voltage recovery frequently emerges during thermal runaway initiation, a phenomenon that conventional theories fail to explain. Here, our study develops a new multiphysics mechanism to explain voltage recovery, termed gas-driven short disconnection, whereby internal gas pressure mechanically disengages short-circuit contacts and causes the cell voltage to rebound. This mechanism incorporates gas generation and its structural impact on the short circuit. Real-time optical and thermal imaging and X-ray computed tomography reveal fluid-structure interaction between internal gas flow and adjacent shorting contacts. We establish a mechanistic framework linking gas-driven short disconnection to cell-level voltage and temperature responses, elucidating the extension–truncation pattern of voltage recovery. Furthermore, thermal regime maps show that a voltage recovery duration exceeding 5 s correlates with limited temperature rise below 150 °C, indicating that sustained short-circuit disconnection suppresses Joule heating. Additionally, a dimensionless criterion is deduced from scaling analysis for physical plausibility of gas-driven short disconnection in mechanically abused cells. This finding inspires smart venting control, which regulates gas release to maintain the internal pressure while dissipating gas enthalpy, thereby providing a device-level strategy for thermal runaway mitigation.

Gas-driven short disconnection↗
DOE OSTI2023

Dynamic structural and microstructural responses of a metal–organic framework type material to carbon dioxide under dual gas flow and supercritical conditions

The structural and microstructural responses of a model metal–organic framework material, Ni(3-methyl-4,4'-bipyridine)[Ni(CN) 4 ] (Ni-BpyMe or PICNIC-21), to CO 2 adsorption and desorption are reported for in situ small-angle X-ray scattering and X-ray diffraction measurements under different gas pressure conditions for two technologically important cases. These conditions are single or dual gas flow (CO 2 with N 2 , CH 4 or H 2 at sub-critical CO 2 partial pressures and ambient temperatures) and supercritical CO 2 (with static pressures and temperatures adjusted to explore the gas, liquid and supercritical fluid regimes on the CO 2 phase diagram). The experimental results are compared with density functional theory calculations that seek to predict where CO 2 and other gas molecules are accommodated within the sorbent structure as a function of gas pressure conditions, and hence the degree of swelling and contraction in the associated structure spacings and void spaces. Furthermore, these predictions illustrate the insights that can be gained concerning how such sorbents can be designed or modified to optimize the desired gas sorption properties relevant to enhanced gas recovery or to addressing carbon dioxide reduction through carbon mitigation, or even direct air capture of CO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Standing shock prevents propagation of sparks in supersonic explosive flows

Abstract Volcanic jet flows in explosive eruptions emit radio frequency signatures, indicative of their fluid dynamic and electrostatic conditions. The emissions originate from sparks supported by an electric field built up by the ejected charged volcanic particles. When shock-defined, low-pressure regions confine the sparks, the signatures may be limited to high-frequency content corresponding to the early components of the avalanche-streamer-leader hierarchy. Here, we image sparks and a standing shock together in a transient supersonic jet of micro-diamonds entrained in argon. Fluid dynamic and kinetic simulations of the experiment demonstrate that the observed sparks originate upstream of the standing shock. The sparks are initiated in the rarefaction region, and cut off at the shock, which would limit their radio frequency emissions to a tell-tale high-frequency regime. We show that sparks transmit an impression of the explosive flow, and open the way for novel instrumentation to diagnose currently inaccessible explosive phenomena.

58 GEOSCIENCES↗
DOE OSTI2026

The Kinetic-Energy–Momentum–Mass 5-Flux of a Baryon Fluid in Bargmann Spacetimes

A Bargmann spacetime is a constrained five-dimensional setting that, while introducing no new physical degrees of freedom beyond those of ordinary four-dimensional spacetime, permits Galilei physics to be expressed with a tensor formalism that respects the distinction between mass and energy while affording the conceptual and technical advantages of a spacetime metric. This framework offers a route to a strong-field ‘Galilei general relativity’ approximating the usual Poincaré general relativity introduced by Einstein. In preparation for modeling core-collapse supernovae, where such an approximation would be useful, this work generalizes the kinetic-energy–momentum–mass 5-flux 𝒯 and its associated spacetime tensor law from a simple fluid of constant particle mass to a baryon fluid whose multiple nuclear species can interconvert rest mass and internal energy. The spacetime tensor law on Bargmann–Galilei spacetime 𝐵𝒢 and its decompositions relative to comoving (‘Lagrangian’) and fiducial (‘Eulerian’) observers are derived in detail. The formalism is rendered more suitable for core-collapse supernova modeling by an extension from strict 𝐵𝒢 to a regime that might be denoted as 𝐵𝒢+: microscopically Poincaré yet macroscopically Galilei. This extension accommodates energy generation by nuclear composition changes and allows comoving energy density and pressure to contribute relative to mass density, while preserving the simplifications of Galilei bulk fluid flow and the streamlined geometry governed by the Bargmann–Galilei spacetime metric.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗
DOE OSTI2020

Investigation of countercurrent flow profile and liquid holdup in random packed column with local CFD data

Liquid holdup and mass transfer area are critical parameters used in packed column design and CO 2 capture efficiency prediction. Compared to experimental methods, 3D computational fluid dynamics (CFD) can provide high-resolution data for studying countercurrent flow hydrodynamics in a packed column. In this paper, a framework was established for modeling the liquid–gas countercurrent flow in a random packed column with pall rings. Besides the column-averaged information, the radial pall ring distribution, velocity, and liquid holdup profiles are obtained in the simulations. The radial two-phase velocity and liquid holdup profiles show significant fluctuations caused by the non-uniform pall ring distribution. The entrance effect and wall influence are also observed in the packed column. By using local countercurrent flow information, the packing specific area $a_p$ and liquid velocity $u_L$ range are significantly expanded for liquid holdup hL prediction. This enables us to develop the liquid holdup correlation with one packing size and six CFD runs only. The proposed power law of the liquid holdup correlation $h_L~u_L^{0.44}$ indicates the random packed column falls in a viscous to turbulent transition regime with a Reynolds Number range of [6.7–40.2]. The derived liquid holdup correlation is in good agreement with existing correlations that are developed using the column-averaged experimental data.

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

Effects of variable deceleration periods on Rayleigh-Taylor instability with acceleration reversals

The dynamics of an interfacial flow that is initially Rayleigh-Taylor unstable but becomes statically stable for some intermediate period due to the reversal of the externally imposed acceleration field is studied. We discuss scenarios that consider both single and double-acceleration reversals. The accel-decel (AD) case consists of a single reversal imposed at an instant after the constant acceleration instability has entered a self-similar regime. The layer of mixed fluid ceases to grow upon acceleration reversal, and the dominant mechanics are due to internal wave oscillations. Variation of mass flux and the Reynolds stress anisotropy is observed due to the action of the internal waves. Here, a second reversal of the AD case that is termed as accel-decel-accel, ADA is then explored; the response of the mixing layer is shown to depend strongly on the duration and the periodicity of the Reynolds stress anisotropy of the mixing layer during the deceleration period. We explore the effect of this variable deceleration period after the second acceleration reversal where the flow once again becomes Rayleigh-Taylor unstable based on metrics that include the integral mixing-layer width, bubble and spike amplitudes, mass flux, Reynolds stress anisotropy tensor, and the molecular mixing parameter.

42 ENGINEERING↗
DOE OSTI2020

Enhancement of chemical oscillations by self-generated convective flows

Chemical feedback loops in fluids can produce not only chemical oscillations, but also density variations that generate solutal buoyancy forces, which in turn initiate fluid flow. Using analytical and computational models, we herein examine how the reaction-induced flows alter the chemical oscillations in a fluid-filled chamber whose top and bottom walls are coated with different enzymes. Due to this chemo-fluidic coupling, the systems form oscillating flow patterns, which combine the characteristic size of the buoyancy-driven convection rolls with the frequency of the chemical oscillations. With changes in the distance between the enzyme-coated walls, the convective flows not only enhance or suppress the chemical oscillations, but also substantially increase the amplitude and frequency of the oscillations and extend the regime of the oscillatory behavior. These design principles can facilitate the development of artificial biochemical networks that act as chemical clocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗