1
|
Mukherjee S, Paul MR. Velocity and geometry of propagating fronts in complex convective flow fields. Phys Rev E 2019; 99:012213. [PMID: 30780300 DOI: 10.1103/physreve.99.012213] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2018] [Indexed: 11/07/2022]
Abstract
We numerically study the propagation of reacting fronts through three-dimensional flow fields composed of convection rolls that include time-independent cellular flow, spatiotemporally chaotic flow, and weakly turbulent flow. We quantify the asymptotic front velocity and determine its scaling with system parameters including the local angle of the convection rolls relative to the direction of front propagation. For cellular flow fields, the orientation of the convection rolls has a significant effect upon the front velocity and the front geometry remains relatively smooth. However, for chaotic and weakly turbulent flow fields, the front velocity depends upon the geometric complexity of the wrinkled front interface and does not depend significantly upon the local orientation of the convection rolls. Using the box counting dimension we find that the front interface is fractal for chaotic and weakly turbulent flows with a dimension that increases with flow complexity.
Collapse
Affiliation(s)
- S Mukherjee
- Department of Biomedical Engineering and Mechanics, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - M R Paul
- Department of Mechanical Engineering, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
2
|
Doan M, Simons JJ, Lilienthal K, Solomon T, Mitchell KA. Barriers to front propagation in laminar, three-dimensional fluid flows. Phys Rev E 2018; 97:033111. [PMID: 29776060 DOI: 10.1103/physreve.97.033111] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2017] [Indexed: 11/07/2022]
Abstract
We present experiments on one-way barriers that block reaction fronts in a fully three-dimensional (3D) fluid flow. Fluorescent Belousov-Zhabotinsky reaction fronts are imaged with laser-scanning in a laminar, overlapping vortex flow. The barriers are analyzed with a 3D extension to burning invariant manifold (BIM) theory that was previously applied to two-dimensional advection-reaction-diffusion processes. We discover tube and sheet barriers that guide the front evolution. The experimentally determined barriers are explained by BIMs calculated from a model of the flow.
Collapse
Affiliation(s)
- Minh Doan
- Department of Physics and Astronomy, Bucknell University, Lewisburg, Pennsylvania 17837, USA
| | - J J Simons
- Department of Physics and Astronomy, Bucknell University, Lewisburg, Pennsylvania 17837, USA
| | - Katherine Lilienthal
- Department of Physics and Astronomy, Bucknell University, Lewisburg, Pennsylvania 17837, USA
| | - Tom Solomon
- Department of Physics and Astronomy, Bucknell University, Lewisburg, Pennsylvania 17837, USA
| | - Kevin A Mitchell
- School of Natural Sciences, University of California, Merced, California 95344, USA
| |
Collapse
|
3
|
Nevins TD, Kelley DH. Front tracking velocimetry in advection-reaction-diffusion systems. CHAOS (WOODBURY, N.Y.) 2018; 28:043122. [PMID: 31906630 DOI: 10.1063/1.5020055] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
In advection-reaction-diffusion systems, the spreading of a reactive scalar can be significantly influenced by the flow field in which it grows. In systems with sharp boundaries between reacted and unreacted regions, motion of the reaction fronts that lie at those boundaries can quantify spreading. Here, we present an algorithm for measuring the velocity of reaction fronts in the presence of flow, expanding previous work on tracking reaction fronts without flow. The algorithm provides localized measurements of front speed and can distinguish its two components: one from chemical dynamics and another from the underlying flow. We validate that the algorithm returns the expected front velocity components in two simulations and then show that in complex experimental flows, the measured front velocity maps fronts from one time step to the next self-consistently. Finally, we observe a variation of the chemical speed with flow speed in a variety of experiments with different time scales and length scales.
Collapse
Affiliation(s)
- Thomas D Nevins
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Douglas H Kelley
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
4
|
Wang J, Tithof J, Nevins TD, Colón RO, Kelley DH. Optimal stretching in the reacting wake of a bluff body. CHAOS (WOODBURY, N.Y.) 2017; 27:123109. [PMID: 29289053 DOI: 10.1063/1.5004649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We experimentally study spreading of the Belousov-Zhabotinsky reaction behind a bluff body in a laminar flow. Locations of reacted regions (i.e., regions with high product concentration) correlate with a moderate range of Lagrangian stretching and that range is close to the range of optimal stretching previously observed in topologically different flows [T. D. Nevins and D. H. Kelley, Phys. Rev. Lett. 117, 164502 (2016)]. The previous work found optimal stretching in a closed, vortex dominated flow, but this article uses an open flow and only a small area of appreciable vorticity. We hypothesize that optimal stretching is common in advection-reaction-diffusion systems with an excitation threshold, including excitable and bistable systems, and that the optimal range depends on reaction chemistry and not on flow shape or characteristic speed. Our results may also give insight into plankton blooms behind islands in ocean currents.
Collapse
Affiliation(s)
- Jinge Wang
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| | - Jeffrey Tithof
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| | - Thomas D Nevins
- Department of Physics & Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Rony O Colón
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| | - Douglas H Kelley
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
5
|
Nevins TD, Kelley DH. Front tracking for quantifying advection-reaction-diffusion. CHAOS (WOODBURY, N.Y.) 2017; 27:043105. [PMID: 28456164 DOI: 10.1063/1.4979668] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We present an algorithm for measuring the speed and thickness of reaction fronts, and from those quantities, the diffusivity and the reaction rate of the active chemical species. This front-tracking algorithm provides local measurements suitable for statistics and requires only a sequence of concentration fields. Though our eventual goal is front tracking in advection-reaction-diffusion, here we demonstrate the algorithm in reaction-diffusion. We test the algorithm with validation data in which front speed and thickness are prescribed, as well as simulation results in which diffusivity and reaction rate are prescribed. In all tests, measurements closely match true values. We apply the algorithm to laboratory experiments using the Belousov-Zhabotinsky reaction, producing speed, diffusivity, and reaction rate measurements that are statistically more robust than in prior studies. Finally, we use thickness measurements to quantify the concentration profile of chemical waves in the reaction.
Collapse
Affiliation(s)
- Thomas D Nevins
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Douglas H Kelley
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
6
|
Nevins TD, Kelley DH. Optimal Stretching in Advection-Reaction-Diffusion Systems. PHYSICAL REVIEW LETTERS 2016; 117:164502. [PMID: 27792376 DOI: 10.1103/physrevlett.117.164502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/27/2016] [Indexed: 06/06/2023]
Abstract
We investigate growth of the excitable Belousov-Zhabotinsky reaction in chaotic, time-varying flows. In slow flows, reacted regions tend to lie near vortex edges, whereas fast flows restrict reacted regions to vortex cores. We show that reacted regions travel toward vortex centers faster as flow speed increases, but nonreactive scalars do not. For either slow or fast flows, reaction is promoted by the same optimal range of the local advective stretching, but stronger stretching causes reaction blowout and can hinder reaction from spreading. We hypothesize that optimal stretching and blowout occur in many advection-diffusion-reaction systems, perhaps creating ecological niches for phytoplankton in the ocean.
Collapse
Affiliation(s)
- Thomas D Nevins
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Douglas H Kelley
- Department of Mechanical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
7
|
Peacock T, Froyland G, Haller G. Introduction to Focus Issue: Objective Detection of Coherent Structures. CHAOS (WOODBURY, N.Y.) 2015; 25:087201. [PMID: 26328571 DOI: 10.1063/1.4928894] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Affiliation(s)
- T Peacock
- Mechanical Engineering Department, Massachusetts Institute of Technology, Cambridge, Massachusetts 20139, USA
| | - G Froyland
- School of Mathematics and Statistics, University of New South Wales, Sydney, Sydney NSW 2052, Australia
| | - G Haller
- Institute for Mechanical Systems, ETH Zürich, Zurich 8092, Switzerland
| |
Collapse
|
8
|
Mahoney JR, Mitchell KA. Finite-time barriers to front propagation in two-dimensional fluid flows. CHAOS (WOODBURY, N.Y.) 2015; 25:087404. [PMID: 26328575 DOI: 10.1063/1.4922026] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Recent theoretical and experimental investigations have demonstrated the role of certain invariant manifolds, termed burning invariant manifolds (BIMs), as one-way dynamical barriers to reaction fronts propagating within a flowing fluid. These barriers form one-dimensional curves in a two-dimensional fluid flow. In prior studies, the fluid velocity field was required to be either time-independent or time-periodic. In the present study, we develop an approach to identify prominent one-way barriers based only on fluid velocity data over a finite time interval, which may have arbitrary time-dependence. We call such a barrier a burning Lagrangian coherent structure (bLCS) in analogy to Lagrangian coherent structures (LCSs) commonly used in passive advection. Our approach is based on the variational formulation of LCSs using curves of stationary "Lagrangian shear," introduced by Farazmand et al. [Physica D 278-279, 44 (2014)] in the context of passive advection. We numerically validate our technique by demonstrating that the bLCS closely tracks the BIM for a time-independent, double-vortex channel flow with an opposing "wind."
Collapse
Affiliation(s)
- John R Mahoney
- School of Natural Sciences, University of California, Merced, California 95343
| | - Kevin A Mitchell
- School of Natural Sciences, University of California, Merced, California 95343
| |
Collapse
|