1
|
Hallatschek O, Datta SS, Drescher K, Dunkel J, Elgeti J, Waclaw B, Wingreen NS. Proliferating active matter. NATURE REVIEWS. PHYSICS 2023; 5:1-13. [PMID: 37360681 PMCID: PMC10230499 DOI: 10.1038/s42254-023-00593-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Accepted: 05/02/2023] [Indexed: 06/28/2023]
Abstract
The fascinating patterns of collective motion created by autonomously driven particles have fuelled active-matter research for over two decades. So far, theoretical active-matter research has often focused on systems with a fixed number of particles. This constraint imposes strict limitations on what behaviours can and cannot emerge. However, a hallmark of life is the breaking of local cell number conservation by replication and death. Birth and death processes must be taken into account, for example, to predict the growth and evolution of a microbial biofilm, the expansion of a tumour, or the development from a fertilized egg into an embryo and beyond. In this Perspective, we argue that unique features emerge in these systems because proliferation represents a distinct form of activity: not only do the proliferating entities consume and dissipate energy, they also inject biomass and degrees of freedom capable of further self-proliferation, leading to myriad dynamic scenarios. Despite this complexity, a growing number of studies document common collective phenomena in various proliferating soft-matter systems. This generality leads us to propose proliferation as another direction of active-matter physics, worthy of a dedicated search for new dynamical universality classes. Conceptual challenges abound, from identifying control parameters and understanding large fluctuations and nonlinear feedback mechanisms to exploring the dynamics and limits of information flow in self-replicating systems. We believe that, by extending the rich conceptual framework developed for conventional active matter to proliferating active matter, researchers can have a profound impact on quantitative biology and reveal fascinating emergent physics along the way.
Collapse
Affiliation(s)
- Oskar Hallatschek
- Departments of Physics and Integrative Biology, University of California, Berkeley, CA US
- Peter Debye Institute for Soft Matter Physics, Leipzig University, Leipzig, Germany
| | - Sujit S. Datta
- Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ USA
| | | | - Jörn Dunkel
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA USA
| | - Jens Elgeti
- Theoretical Physics of Living Matter, Institute of Biological Information Processing, Forschungszentrum Jülich, Jülich, Germany
| | - Bartek Waclaw
- Dioscuri Centre for Physics and Chemistry of Bacteria, Institute of Physical Chemistry PAN, Warsaw, Poland
- School of Physics and Astronomy, The University of Edinburgh, JCMB, Edinburgh, UK
| | - Ned S. Wingreen
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ USA
- Department of Molecular Biology, Princeton University, Princeton, NJ USA
| |
Collapse
|
2
|
Cao X, Le Doussal P, Rosso A. Clusters in an Epidemic Model with Long-Range Dispersal. PHYSICAL REVIEW LETTERS 2022; 129:108301. [PMID: 36112459 DOI: 10.1103/physrevlett.129.108301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Revised: 05/27/2022] [Accepted: 08/17/2022] [Indexed: 06/15/2023]
Abstract
In the presence of long-range dispersal, epidemics spread in spatially disconnected regions known as clusters. Here, we characterize exactly their statistical properties in a solvable model, in both the supercritical (outbreak) and critical regimes. We identify two diverging length scales, corresponding to the bulk and the outskirt of the epidemic. We reveal a nontrivial critical exponent that governs the cluster number and the distribution of their sizes and of the distances between them. We also discuss applications to depinning avalanches with long-range elasticity.
Collapse
Affiliation(s)
- Xiangyu Cao
- Laboratoire de Physique de l'École normale supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, F-75005 Paris, France
| | - Pierre Le Doussal
- Laboratoire de Physique de l'École normale supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, F-75005 Paris, France
| | - Alberto Rosso
- Université Paris-Saclay, CNRS, LPTMS, 91405 Orsay, France
| |
Collapse
|
3
|
Benzi R, Nelson DR, Shankar S, Toschi F, Zhu X. Spatial population genetics with fluid flow. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2022; 85:096601. [PMID: 35853344 DOI: 10.1088/1361-6633/ac8231] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/02/2022] [Accepted: 07/19/2022] [Indexed: 06/15/2023]
Abstract
The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2D), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1D and 2D flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how thesegeneticdrops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well.
Collapse
Affiliation(s)
- Roberto Benzi
- Department of Physics and INFN, University of Rome Tor Vergata, I-00133 Rome, Italy
| | - David R Nelson
- Department of Physics, Harvard University, Cambridge, MA 02138, United States of America
| | - Suraj Shankar
- Department of Physics, Harvard University, Cambridge, MA 02138, United States of America
| | - Federico Toschi
- Department of Applied Physics, Department of Mathematics and Computer Science, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
- CNR-IAC, I-00185 Rome, Italy
| | - Xiaojue Zhu
- Solar and Stellar Interiors, Max Planck Institute for Solar System Research, Justus-von-Liebig-Weg 3, Göttingen 37077, Germany
- Center of Mathematical Sciences and Applications, and School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, United States of America
| |
Collapse
|
4
|
Labavić D, Loverdo C, Bitbol AF. Hydrodynamic flow and concentration gradients in the gut enhance neutral bacterial diversity. Proc Natl Acad Sci U S A 2022; 119:e2108671119. [PMID: 34969835 PMCID: PMC8740595 DOI: 10.1073/pnas.2108671119] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/08/2021] [Indexed: 01/23/2023] Open
Abstract
The gut microbiota features important genetic diversity, and the specific spatial features of the gut may shape evolution within this environment. We investigate the fixation probability of neutral bacterial mutants within a minimal model of the gut that includes hydrodynamic flow and resulting gradients of food and bacterial concentrations. We find that this fixation probability is substantially increased, compared with an equivalent well-mixed system, in the regime where the profiles of food and bacterial concentration are strongly spatially dependent. Fixation probability then becomes independent of total population size. We show that our results can be rationalized by introducing an active population, which consists of those bacteria that are actively consuming food and dividing. The active population size yields an effective population size for neutral mutant fixation probability in the gut.
Collapse
Affiliation(s)
- Darka Labavić
- CNRS, Institut de Biologie Paris-Seine, Laboratoire Jean Perrin (UMR 8237), Sorbonne Université, F-75005 Paris, France
| | - Claude Loverdo
- CNRS, Institut de Biologie Paris-Seine, Laboratoire Jean Perrin (UMR 8237), Sorbonne Université, F-75005 Paris, France;
| | - Anne-Florence Bitbol
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland;
- SIB Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| |
Collapse
|
5
|
Turbulent coherent structures and early life below the Kolmogorov scale. Nat Commun 2020; 11:2192. [PMID: 32366844 PMCID: PMC7198613 DOI: 10.1038/s41467-020-15780-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2019] [Accepted: 03/27/2020] [Indexed: 01/28/2023] Open
Abstract
Major evolutionary transitions, including the emergence of life, likely occurred in aqueous environments. While the role of water’s chemistry in early life is well studied, the effects of water’s ability to manipulate population structure are less clear. Population structure is known to be critical, as effective replicators must be insulated from parasites. Here, we propose that turbulent coherent structures, long-lasting flow patterns which trap particles, may serve many of the properties associated with compartments — collocalization, division, and merging — which are commonly thought to play a key role in the origins of life and other evolutionary transitions. We substantiate this idea by simulating multiple proposed metabolisms for early life in a simple model of a turbulent flow, and find that balancing the turnover times of biological particles and coherent structures can indeed enhance the likelihood of these metabolisms overcoming extinction either via parasitism or via a lack of metabolic support. Our results suggest that group selection models may be applicable with fewer physical and chemical constraints than previously thought, and apply much more widely in aqueous environments. Models of the origin of life generally require a mechanism to structure emerging populations. Here, Krieger et al. develop spatial models showing that coherent structures arising in turbulent flows in aquatic environments could have provided compartments that facilitated the origin of life.
Collapse
|
6
|
Guccione G, Benzi R, Plummer A, Toschi F. Discrete Eulerian model for population genetics and dynamics under flow. Phys Rev E 2020; 100:062105. [PMID: 31962443 DOI: 10.1103/physreve.100.062105] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2019] [Indexed: 11/07/2022]
Abstract
Marine species reproduce and compete while being advected by turbulent flows. It is largely unknown, both theoretically and experimentally, how population dynamics and genetics are changed by the presence of fluid flows. Discrete agent-based simulations in continuous space allow for accurate treatment of advection and number fluctuations, but can be computationally expensive for even modest organism densities. In this report, we propose an algorithm to overcome some of these challenges. We first provide a thorough validation of the algorithm in one and two dimensions without flow. Next, we focus on the case of weakly compressible flows in two dimensions. This models organisms such as phytoplankton living at a specific depth in the three-dimensional, incompressible ocean experiencing upwelling and/or downwelling events. We show that organisms born at sources in a two-dimensional time-independent flow experience an increase in fixation probability.
Collapse
Affiliation(s)
- Giorgia Guccione
- Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands and Department of Physics and INFN, University of Tor Vergata, Via della Ricerca Scientifica 1, I-00133 Rome, Italy
| | - Roberto Benzi
- Department of Physics and INFN, University of Tor Vergata, Via della Ricerca Scientifica 1, I-00133 Rome, Italy
| | - Abigail Plummer
- Department of Physics, Harvard University, 17 Oxford Street, Cambridge, Massachusetts 02138, USA
| | - Federico Toschi
- Department of Applied Physics, Department of Mathematics and Computer Science, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands and CNR-IAC, Via dei Taurini 19, I-00185 Rome, Italy
| |
Collapse
|
7
|
Abstract
Competition between biological species in marine environments is affected by the motion of the surrounding fluid. An effective 2D compressibility can arise, for example, from the convergence and divergence of water masses at the depth at which passively traveling photosynthetic organisms are restricted to live. In this report, we seek to quantitatively study genetics under flow. To this end, we couple an off-lattice agent-based simulation of two populations in 1D to a weakly compressible velocity field-first a sine wave and then a shell model of turbulence. We find for both cases that even in a regime where the overall population structure is approximately unaltered, the flow can significantly diminish the effect of a selective advantage on fixation probabilities. We understand this effect in terms of the enhanced survival of organisms born at sources in the flow and the influence of Fisher genetic waves.
Collapse
|
8
|
Uppal G, Vural DC. Shearing in flow environment promotes evolution of social behavior in microbial populations. eLife 2018; 7:34862. [PMID: 29785930 PMCID: PMC6002248 DOI: 10.7554/elife.34862] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2018] [Accepted: 05/10/2018] [Indexed: 11/23/2022] Open
Abstract
How producers of public goods persist in microbial communities is a major question in evolutionary biology. Cooperation is evolutionarily unstable, since cheating strains can reproduce quicker and take over. Spatial structure has been shown to be a robust mechanism for the evolution of cooperation. Here we study how spatial assortment might emerge from native dynamics and show that fluid flow shear promotes cooperative behavior. Social structures arise naturally from our advection-diffusion-reaction model as self-reproducing Turing patterns. We computationally study the effects of fluid advection on these patterns as a mechanism to enable or enhance social behavior. Our central finding is that flow shear enables and promotes social behavior in microbes by increasing the group fragmentation rate and thereby limiting the spread of cheating strains. Regions of the flow domain with higher shear admit high cooperativity and large population density, whereas low shear regions are devoid of life due to opportunistic mutations. According to the principle of the ‘survival of the fittest’, selfish individuals should be better off compared to peers that cooperate with each other. Indeed, even though a population of organisms benefits from working together, selfish members can exploit the cooperative behavior of others without doing their part. These ‘cheaters’ then use their advantage to reproduce faster and take over the population. Yet, social cooperation is widespread in the natural world, and occurs in creatures as diverse as bacteria and whales. How can it arise and persist then? One idea is that when individuals form distinct groups, the ones with cheaters will perish. Even though a selfish individual will fare better than the rest of its team, overall, cooperating groups will survive more and reproduce faster; ultimately, they will be favored by evolution. This is called group selection. Here, Uppal and Vural examine how the physical properties of the environment can influence the evolution of social interactions between bacteria. To this end, mathematical models are used to simulate how bacteria grow, evolve and drift in a flowing fluid. These are based on equations worked out from the behavior of real-life populations. The results show that flow patterns in a fluid habitat govern the social behavior of bacteria. When different regions of the fluid are moving at different speeds, ‘shear forces’ are created that cause bacterial colonies to distort and occasionally break apart to form two groups. As such, cooperative groups will rapidly form new cooperating colonies, whereas groups with cheaters will reproduce slower or perish. Furthermore, results show that when different areas of the fluid have different shear forces, social cooperation will only prevail in certain places. This makes it possible to use flow patterns to fine tune social evolution so that cooperating bacteria will be confined in a certain region. Outside of this area, these bacteria would be taken over by cheaters and go extinct. Bacteria are both useful and dangerous to humans: for example, certain species can break down pollutants in the water, when others cause deadly infections. These results show it could be possible to control the activity of these microorganisms to our advantage by changing the flow of the fluids in which they live. More broadly, the simulations developed by Uppal and Vural can be applied to a variety of ecosystems where microscopic organisms inhabit fluids, such as plankton flowing in oceanic currents.
Collapse
Affiliation(s)
- Gurdip Uppal
- Department of Physics, University of Notre Dame, Notre Dame, United States
| | - Dervis Can Vural
- Department of Physics, University of Notre Dame, Notre Dame, United States
| |
Collapse
|
9
|
Krieger MS, McAvoy A, Nowak MA. Effects of motion in structured populations. J R Soc Interface 2017; 14:rsif.2017.0509. [PMID: 28978749 DOI: 10.1098/rsif.2017.0509] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Accepted: 09/05/2017] [Indexed: 11/12/2022] Open
Abstract
In evolutionary processes, population structure has a substantial effect on natural selection. Here, we analyse how motion of individuals affects constant selection in structured populations. Motion is relevant because it leads to changes in the distribution of types as mutations march towards fixation or extinction. We describe motion as the swapping of individuals on graphs, and more generally as the shuffling of individuals between reproductive updates. Beginning with a one-dimensional graph, the cycle, we prove that motion suppresses natural selection for death-birth (DB) updating or for any process that combines birth-death (BD) and DB updating. If the rule is purely BD updating, no change in fixation probability appears in the presence of motion. We further investigate how motion affects evolution on the square lattice and weighted graphs. In the case of weighted graphs, we find that motion can be either an amplifier or a suppressor of natural selection. In some cases, whether it is one or the other can be a function of the relative reproductive rate, indicating that motion is a subtle and complex attribute of evolving populations. As a first step towards understanding less restricted types of motion in evolutionary graph theory, we consider a similar rule on dynamic graphs induced by a spatial flow and find qualitatively similar results, indicating that continuous motion also suppresses natural selection.
Collapse
Affiliation(s)
- Madison S Krieger
- Program for Evolutionary Dynamics, Harvard University, One Brattle Square, Suite 6, Cambridge, MA 02138, USA
| | - Alex McAvoy
- Program for Evolutionary Dynamics, Harvard University, One Brattle Square, Suite 6, Cambridge, MA 02138, USA
| | - Martin A Nowak
- Program for Evolutionary Dynamics, Harvard University, One Brattle Square, Suite 6, Cambridge, MA 02138, USA
| |
Collapse
|
10
|
Finite-size effects on bacterial population expansion under controlled flow conditions. Sci Rep 2017; 7:43903. [PMID: 28262769 PMCID: PMC5338255 DOI: 10.1038/srep43903] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2016] [Accepted: 01/26/2017] [Indexed: 11/09/2022] Open
Abstract
The expansion of biological species in natural environments is usually described as the combined effect of individual spatial dispersal and growth. In the case of aquatic ecosystems flow transport can also be extremely relevant as an extra, advection induced, dispersal factor. We designed and assembled a dedicated microfluidic device to control and quantify the expansion of populations of E. coli bacteria under both co-flowing and counter-flowing conditions, measuring the front speed at varying intensity of the imposed flow. At variance with respect to the case of classic advective-reactive-diffusive chemical fronts, we measure that almost irrespective of the counter-flow velocity, the front speed remains finite at a constant positive value. A simple model incorporating growth, dispersion and drift on finite-size hard beads allows to explain this finding as due to a finite volume effect of the bacteria. This indicates that models based on the Fisher-Kolmogorov-Petrovsky-Piscounov equation (FKPP) that ignore the finite size of organisms may be inaccurate to describe the physics of spatial growth dynamics of bacteria.
Collapse
|
11
|
Chatterjee R, Joshi AA, Perlekar P. Front structure and dynamics in dense colonies of motile bacteria: Role of active turbulence. Phys Rev E 2016; 94:022406. [PMID: 27627334 DOI: 10.1103/physreve.94.022406] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2016] [Indexed: 11/07/2022]
Abstract
We study the spreading of a bacterial colony undergoing turbulentlike collective motion. We present two minimalistic models to investigate the interplay between population growth and coherent structures arising from turbulence. Using direct numerical simulation of the proposed models we find that turbulence has two prominent effects on the spatial growth of the colony: (a) the front speed is enhanced, and (b) the front gets crumpled. Both these effects, which we highlight by using statistical tools, are markedly different in our two models. We also show that the crumpled front structure and the passive scalar fronts in random flows are related in certain regimes.
Collapse
Affiliation(s)
- Rayan Chatterjee
- TIFR Centre for Interdisciplinary Sciences, 21 Brundavan Colony, Narsingi, Hyderabad
| | - Abhijeet A Joshi
- TIFR Centre for Interdisciplinary Sciences, 21 Brundavan Colony, Narsingi, Hyderabad
| | - Prasad Perlekar
- TIFR Centre for Interdisciplinary Sciences, 21 Brundavan Colony, Narsingi, Hyderabad
| |
Collapse
|
12
|
Lai H, Xu A, Zhang G, Gan Y, Ying Y, Succi S. Nonequilibrium thermohydrodynamic effects on the Rayleigh-Taylor instability in compressible flows. Phys Rev E 2016; 94:023106. [PMID: 27627391 DOI: 10.1103/physreve.94.023106] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2015] [Indexed: 06/06/2023]
Abstract
The effects of compressibility on Rayleigh-Taylor instability (RTI) are investigated by inspecting the interplay between thermodynamic and hydrodynamic nonequilibrium phenomena (TNE, HNE, respectively) via a discrete Boltzmann model. Two effective approaches are presented, one tracking the evolution of the local TNE effects and the other focusing on the evolution of the mean temperature of the fluid, to track the complex interfaces separating the bubble and the spike regions of the flow. It is found that both the compressibility effects and the global TNE intensity show opposite trends in the initial and the later stages of the RTI. Compressibility delays the initial stage of RTI and accelerates the later stage. Meanwhile, the TNE characteristics are generally enhanced by the compressibility, especially in the later stage. The global or mean thermodynamic nonequilibrium indicators provide physical criteria to discriminate between the two stages of the RTI.
Collapse
Affiliation(s)
- Huilin Lai
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P. O. Box 8009-26, Beijing 100088, China
- School of Mathematics and Computer Science, Fujian Normal University, Fuzhou 350007, China
| | - Aiguo Xu
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P. O. Box 8009-26, Beijing 100088, China
- Center for Applied Physics and Technology, MOE Key Center for High Energy Density Physics Simulations, College of Engineering, Peking University, Beijing 100871, China
| | - Guangcai Zhang
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P. O. Box 8009-26, Beijing 100088, China
| | - Yanbiao Gan
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P. O. Box 8009-26, Beijing 100088, China
- North China Institute of Aerospace Engineering, Langfang 065000, China
| | - Yangjun Ying
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P. O. Box 8009-26, Beijing 100088, China
| | - Sauro Succi
- Istituto Applicazioni Calcolo, CNR, Viale del Policlinico 137, 00161 Roma, Italy
| |
Collapse
|
13
|
Chotibut T, Nelson DR. Evolutionary dynamics with fluctuating population sizes and strong mutualism. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:022718. [PMID: 26382443 DOI: 10.1103/physreve.92.022718] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/18/2014] [Indexed: 06/05/2023]
Abstract
Game theory ideas provide a useful framework for studying evolutionary dynamics in a well-mixed environment. This approach, however, typically enforces a strictly fixed overall population size, deemphasizing natural growth processes. We study a competitive Lotka-Volterra model, with number fluctuations, that accounts for natural population growth and encompasses interaction scenarios typical of evolutionary games. We show that, in an appropriate limit, the model describes standard evolutionary games with both genetic drift and overall population size fluctuations. However, there are also regimes where a varying population size can strongly influence the evolutionary dynamics. We focus on the strong mutualism scenario and demonstrate that standard evolutionary game theory fails to describe our simulation results. We then analytically and numerically determine fixation probabilities as well as mean fixation times using matched asymptotic expansions, taking into account the population size degree of freedom. These results elucidate the interplay between population dynamics and evolutionary dynamics in well-mixed systems.
Collapse
Affiliation(s)
- Thiparat Chotibut
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
| | - David R Nelson
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
| |
Collapse
|
14
|
De Pietro M, van Hinsberg MAT, Biferale L, Clercx HJH, Perlekar P, Toschi F. Clustering of vertically constrained passive particles in homogeneous isotropic turbulence. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:053002. [PMID: 26066244 DOI: 10.1103/physreve.91.053002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/07/2014] [Indexed: 06/04/2023]
Abstract
We analyze the dynamics of small particles vertically confined, by means of a linear restoring force, to move within a horizontal fluid slab in a three-dimensional (3D) homogeneous isotropic turbulent velocity field. The model that we introduce and study is possibly the simplest description for the dynamics of small aquatic organisms that, due to swimming, active regulation of their buoyancy, or any other mechanism, maintain themselves in a shallow horizontal layer below the free surface of oceans or lakes. By varying the strength of the restoring force, we are able to control the thickness of the fluid slab in which the particles can move. This allows us to analyze the statistical features of the system over a wide range of conditions going from a fully 3D incompressible flow (corresponding to the case of no confinement) to the extremely confined case corresponding to a two-dimensional slice. The background 3D turbulent velocity field is evolved by means of fully resolved direct numerical simulations. Whenever some level of vertical confinement is present, the particle trajectories deviate from that of fluid tracers and the particles experience an effectively compressible velocity field. Here, we have quantified the compressibility, the preferential concentration of the particles, and the correlation dimension by changing the strength of the restoring force. The main result is that there exists a particular value of the force constant, corresponding to a mean slab depth approximately equal to a few times the Kolmogorov length scale η, that maximizes the clustering of the particles.
Collapse
Affiliation(s)
- Massimo De Pietro
- Dipartimento di Fisica and Istituto Nazionale di Fisica Nucleare, Università "Tor Vergata," Via della Ricerca Scientifica 1, I-00133 Roma, Italy
| | - Michel A T van Hinsberg
- Department of Applied Physics, J. M. Burgerscentrum, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
| | - Luca Biferale
- Dipartimento di Fisica and Istituto Nazionale di Fisica Nucleare, Università "Tor Vergata," Via della Ricerca Scientifica 1, I-00133 Roma, Italy
| | - Herman J H Clercx
- Department of Applied Physics, J. M. Burgerscentrum, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
| | - Prasad Perlekar
- TIFR Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, 21 Brundavan Colony, Narsingi, Hyderabad 500075, India
| | - Federico Toschi
- Department of Applied Physics and Department of Mathematics and Computer Science, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands and IAC, Consiglio Nazionale delle Ricerche, Via dei Taurini 19, I-00185 Roma, Italy
| |
Collapse
|
15
|
Grošelj D, Jenko F, Frey E. How turbulence regulates biodiversity in systems with cyclic competition. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:033009. [PMID: 25871204 DOI: 10.1103/physreve.91.033009] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/14/2014] [Indexed: 06/04/2023]
Abstract
Cyclic, nonhierarchical interactions among biological species represent a general mechanism by which ecosystems are able to maintain high levels of biodiversity. However, species coexistence is often possible only in spatially extended systems with a limited range of dispersal, whereas in well-mixed environments models for cyclic competition often lead to a loss of biodiversity. Here we consider the dispersal of biological species in a fluid environment, where mixing is achieved by a combination of advection and diffusion. In particular, we perform a detailed numerical analysis of a model composed of turbulent advection, diffusive transport, and cyclic interactions among biological species in two spatial dimensions and discuss the circumstances under which biodiversity is maintained when external environmental conditions, such as resource supply, are uniform in space. Cyclic interactions are represented by a model with three competitors, resembling the children's game of rock-paper-scissors, whereas the flow field is obtained from a direct numerical simulation of two-dimensional turbulence with hyperviscosity. It is shown that the space-averaged dynamics undergoes bifurcations as the relative strengths of advection and diffusion compared to biological interactions are varied.
Collapse
Affiliation(s)
- Daniel Grošelj
- Max-Planck-Institut für Plasmaphysik, Boltzmannstraße 2, D-85748 Garching, Germany
| | - Frank Jenko
- Max-Planck-Institut für Plasmaphysik, Boltzmannstraße 2, D-85748 Garching, Germany
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095-1547, USA
| | - Erwin Frey
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstraße 37, D-80333 München, Germany
| |
Collapse
|
16
|
Abstract
The spreading of evolutionary novelties across populations is the central element of adaptation. Unless populations are well mixed (like bacteria in a shaken test tube), the spreading dynamics depend not only on fitness differences but also on the dispersal behavior of the species. Spreading at a constant speed is generally predicted when dispersal is sufficiently short ranged, specifically when the dispersal kernel falls off exponentially or faster. However, the case of long-range dispersal is unresolved: Although it is clear that even rare long-range jumps can lead to a drastic speedup--as air-traffic-mediated epidemics show--it has been difficult to quantify the ensuing stochastic dynamical process. However, such knowledge is indispensable for a predictive understanding of many spreading processes in natural populations. We present a simple iterative scaling approximation supported by simulations and rigorous bounds that accurately predicts evolutionary spread, which is determined by a trade-off between frequency and potential effectiveness of long-distance jumps. In contrast to the exponential laws predicted by deterministic "mean-field" approximations, we show that the asymptotic spatial growth is according to either a power law or a stretched exponential, depending on the tails of the dispersal kernel. More importantly, we provide a full time-dependent description of the convergence to the asymptotic behavior, which can be anomalously slow and is relevant even for long times. Our results also apply to spreading dynamics on networks with a spectrum of long-range links under certain conditions on the probabilities of long-distance travel: These are relevant for the spread of epidemics.
Collapse
|
17
|
Elperin T, Kleeorin N, Liberman M, Rogachevskii I. Turbulent diffusion of chemically reacting gaseous admixtures. Phys Rev E 2014; 90:053001. [PMID: 25493875 DOI: 10.1103/physreve.90.053001] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2014] [Indexed: 11/07/2022]
Abstract
We study turbulent diffusion of chemically reacting gaseous admixtures in a developed turbulence. In our previous study [Phys. Rev. Lett. 80, 69 (1998)PRLTAO0031-900710.1103/PhysRevLett.80.69] using a path-integral approach for a delta-correlated in a time random velocity field, we demonstrated a strong modification of turbulent transport in fluid flows with chemical reactions or phase transitions. In the present study we use the spectral τ approximation that is valid for large Reynolds and Peclet numbers and show that turbulent diffusion of the reacting species can be strongly depleted by a large factor that is the ratio of turbulent and chemical times (turbulent Damköhler number). We have demonstrated that the derived theoretical dependence of a turbulent diffusion coefficient versus the turbulent Damköhler number is in good agreement with that obtained previously in the numerical modeling of a reactive front propagating in a turbulent flow and described by the Kolmogorov-Petrovskii-Piskunov-Fisher equation. We have found that turbulent cross-effects, e.g., turbulent mutual diffusion of gaseous admixtures and turbulent Dufour effect of the chemically reacting gaseous admixtures, are less sensitive to the values of stoichiometric coefficients. The mechanisms of the turbulent cross-effects differ from the molecular cross-effects known in irreversible thermodynamics. In a fully developed turbulence and at large Peclet numbers the turbulent cross-effects are much larger than the molecular ones. The obtained results are applicable also to heterogeneous phase transitions.
Collapse
Affiliation(s)
- T Elperin
- The Pearlstone Center for Aeronautical Engineering Studies, Department of Mechanical Engineering, Ben-Gurion University of the Negev, P. O. Box 653, Beer-Sheva 84105, Israel
| | - N Kleeorin
- The Pearlstone Center for Aeronautical Engineering Studies, Department of Mechanical Engineering, Ben-Gurion University of the Negev, P. O. Box 653, Beer-Sheva 84105, Israel
| | - M Liberman
- Nordita, KTH Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, 10691 Stockholm, Sweden and Moscow Institute of Physics and Technology, Dolgoprudnyi, 141700, Russia
| | - I Rogachevskii
- The Pearlstone Center for Aeronautical Engineering Studies, Department of Mechanical Engineering, Ben-Gurion University of the Negev, P. O. Box 653, Beer-Sheva 84105, Israel
| |
Collapse
|
18
|
Kussell E, Vucelja M. Non-equilibrium physics and evolution--adaptation, extinction, and ecology: a key issues review. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2014; 77:102602. [PMID: 25303141 DOI: 10.1088/0034-4885/77/10/102602] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Evolutionary dynamics in nature constitute an immensely complex non-equilibrium process. We review the application of physical models of evolution, by focusing on adaptation, extinction, and ecology. In each case, we examine key concepts by working through examples. Adaptation is discussed in the context of bacterial evolution, with a view toward the relationship between growth rates, mutation rates, selection strength, and environmental changes. Extinction dynamics for an isolated population are reviewed, with emphasis on the relation between timescales of extinction, population size, and temporally correlated noise. Ecological models are discussed by focusing on the effect of spatial interspecies interactions on diversity. Connections between physical processes--such as diffusion, turbulence, and localization--and evolutionary phenomena are highlighted.
Collapse
Affiliation(s)
- E Kussell
- Department of Biology and Center for Genomics and Systems Biology, New York University, 12 Waverly Place, New York, NY 10003, USA. Department of Physics, New York University, New York, NY 10003, USA
| | | |
Collapse
|
19
|
Volk R, Mauger C, Bourgoin M, Cottin-Bizonne C, Ybert C, Raynal F. Chaotic mixing in effective compressible flows. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:013027. [PMID: 25122389 DOI: 10.1103/physreve.90.013027] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/28/2014] [Indexed: 06/03/2023]
Abstract
We study numerically joint mixing of salt and colloids by chaotic advection and how salt inhomogeneities accelerate or delay colloid mixing by inducing a velocity drift V(dp) between colloids and fluid particles as proposed in recent experiments [J. Deseigne et al., Soft Matter 10, 4795 (2014)]. We demonstrate that because the drift velocity is no longer divergence free, small variations to the total velocity field drastically affect the evolution of colloid variance σ(2) = 〈C(2)〉-〈C〉(2). A consequence is that mixing strongly depends on the mutual coherence between colloid and salt concentration fields, the short time evolution of scalar variance being governed by a new variance production term P = -〈C(2)∇ · V(dp)〉/2 when scalar gradients are not developed yet so that dissipation is weak. Depending on initial conditions, mixing is then delayed or enhanced, and it is possible to find examples for which the two regimes (fast mixing followed by slow mixing) are observed consecutively when the variance source term reverses its sign. This is indeed the case for localized patches modeled as Gaussian concentration profiles.
Collapse
Affiliation(s)
- R Volk
- Laboratoire de Physique de l'ENS de Lyon, CNRS UMR5672 and Université de Lyon, Lyon, France
| | - C Mauger
- LMFA, CNRS UMR5509, École Centrale Lyon, INSA de Lyon and Université de Lyon, Lyon, France
| | - M Bourgoin
- LEGI, CNRS UMR5519, Université Joseph Fourier, Université de Grenoble, Grenoble INP, France
| | - C Cottin-Bizonne
- Institut Lumière Matière, CNRS UMR5306, Université Claude Bernard Lyon 1, Université de Lyon, Lyon, France
| | - C Ybert
- Institut Lumière Matière, CNRS UMR5306, Université Claude Bernard Lyon 1, Université de Lyon, Lyon, France
| | - F Raynal
- LMFA, CNRS UMR5509, École Centrale Lyon, INSA de Lyon and Université de Lyon, Lyon, France
| |
Collapse
|
20
|
Bianco F, Chibbaro S, Vergni D, Vulpiani A. Front speed in reactive compressible stirred media. Phys Rev E 2013; 87:042924. [PMID: 23679506 DOI: 10.1103/physreve.87.042924] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2013] [Indexed: 11/07/2022]
Abstract
We investigated a nonlinear advection-diffusion-reaction equation for a passive scalar field. The purpose is to understand how the compressibility can affect the front dynamics and the bulk burning rate. We study two classes of flows: periodic shear flow and cellular flow, analyzing the system by varying the extent of compressibility and the reaction rate. We find that the bulk burning rate v(f) in a shear flow increases with compressibility intensity ε, following the relation Δv(f)~ε(2). Furthermore, the faster the reaction is, the more important the difference is with respect to the laminar case. The effect has been quantitatively measured, and it turns out to be generally small. For the cellular flow, two extreme cases have been investigated, with the whole perturbation situated either in the center of the vortex or in the periphery. The dependence in this case does not show a monotonic scaling with different behavior in the two cases. The enhancing remains modest and is always less than 20%.
Collapse
Affiliation(s)
- Federico Bianco
- Institut D'Alembert, University Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 5, France
| | | | | | | |
Collapse
|
21
|
Growth, competition and cooperation in spatial population genetics. Theor Popul Biol 2013; 84:72-86. [PMID: 23298763 DOI: 10.1016/j.tpb.2012.12.002] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2012] [Revised: 12/18/2012] [Accepted: 12/20/2012] [Indexed: 11/21/2022]
|
22
|
De Monte S, Cotté C, d'Ovidio F, Lévy M, Le Corre M, Weimerskirch H. Frigatebird behaviour at the ocean-atmosphere interface: integrating animal behaviour with multi-satellite data. J R Soc Interface 2012; 9:3351-8. [PMID: 22951344 DOI: 10.1098/rsif.2012.0509] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Marine top predators such as seabirds are useful indicators of the integrated response of the marine ecosystem to environmental variability at different scales. Large-scale physical gradients constrain seabird habitat. Birds however respond behaviourally to physical heterogeneity at much smaller scales. Here, we use, for the first time, three-dimensional GPS tracking of a seabird, the great frigatebird (Fregata minor), in the Mozambique Channel. These data, which provide at the same time high-resolution vertical and horizontal positions, allow us to relate the behaviour of frigatebirds to the physical environment at the (sub-)mesoscale (10-100 km, days-weeks). Behavioural patterns are classified based on the birds' vertical displacement (e.g. fast/slow ascents and descents), and are overlaid on maps of physical properties of the ocean-atmosphere interface, obtained by a nonlinear analysis of multi-satellite data. We find that frigatebirds modify their behaviours concurrently to transport and thermal fronts. Our results suggest that the birds' co-occurrence with these structures is a consequence of their search not only for food (preferentially searched over thermal fronts) but also for upward vertical wind. This is also supported by their relationship with mesoscale patterns of wind divergence. Our multi-disciplinary method can be applied to forthcoming high-resolution animal tracking data, and aims to provide a mechanistic understanding of animals' habitat choice and of marine ecosystem responses to environmental change.
Collapse
Affiliation(s)
- Silvia De Monte
- Centre National de Recherche Scientifique, Unité Mixte de Recherche 7625, Écologie et Évolution, 46 rue d'Ulm, 75005 Paris, France.
| | | | | | | | | | | |
Collapse
|
23
|
Abstract
Three areas where time-independent disorder plays a key role in biological dynamics far from equilibrium are reviewed. We first discuss the anomalous localization dynamics that arises when a single species spreads in space and time via diffusion and fluid advection in the presence of frozen heterogeneities in the growth rate. Next we treat the unzipping of double-stranded DNA as a function of force and temperature, a challenge that must be surmounted every time a cell divides. Heterogeneity in the DNA sequence dominates the physics of single-molecule force-extension curves for a broad range of forces upon approaching a sharp unzipping transition. The dynamics of the unzipping fork exhibits anomalous drift and diffusion in a similar range above this transition, with energy barriers that scale as the square root of the genome size. Finally, we describe how activated peptidoglycan strand extension sites, called dislocations in materials science, can mediate the growth of bacterial cell walls. Enzymatically driven circumferential motions of a few dozen of these defects are sufficient to describe the exponential elongation rates observed in experiments on Escherichia coli in a nutrient-rich environment. However, long-range elastic forces transmitted by the peptidoglycan meshwork cause the moving dislocations to interact not only with each other, but also with a disorderly array of frozen, inactivated strand ends.
Collapse
Affiliation(s)
- David R. Nelson
- Department of Physics, Department of Molecular and Cellular Biology, and School of Engineering and Applied Science, Harvard University, Cambridge, Massachusetts 02138
| |
Collapse
|
24
|
Pigolotti S, Benzi R, Jensen MH, Nelson DR. Population genetics in compressible flows. PHYSICAL REVIEW LETTERS 2012; 108:128102. [PMID: 22540626 DOI: 10.1103/physrevlett.108.128102] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/17/2011] [Indexed: 05/31/2023]
Abstract
We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We investigate three examples of compressible advecting fields: a shell model of turbulence, a sinusoidal velocity field and a linear velocity sink. In all cases, advection leads to a striking drop in the fixation time, as well as a large reduction in the global carrying capacity. We find localization on convergence zones, and very rapid extinction compared to well-mixed populations. For a linear velocity sink, one finds a bimodal distribution of fixation times. The long-lived states in this case are demixed configurations with a single interface, whose location depends on the fitness advantage.
Collapse
Affiliation(s)
- Simone Pigolotti
- The Niels Bohr Institut, Blegdamsvej 17, DK-2100 Copenhagen, Denmark
| | | | | | | |
Collapse
|
25
|
Perlekar P, Benzi R, Pigolotti S, Toschi F. Particle algorithms for population dynamics in flows. ACTA ACUST UNITED AC 2011. [DOI: 10.1088/1742-6596/333/1/012013] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
26
|
Burgener T, Kadau D, Herrmann HJ. Simulation of particle mixing in turbulent channel flow due to intrinsic fluid velocity fluctuation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:066301. [PMID: 21797471 DOI: 10.1103/physreve.83.066301] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/25/2010] [Indexed: 05/31/2023]
Abstract
We combine a discrete-element-method simulation with a stochastic process to model the movement of spherical particles in a turbulent channel flow. With this model we investigate the mixing properties of two species of particles flowing through the channel. We find a linear increase of the mixing zone with the length of the pipe. Flows at different Reynolds number are studied. Below a critical Reynolds number at the Taylor microscale of around Rc ≈ 300 the mixing rate is strongly dependent on the Reynolds number. Above Rc the mixing rate stays nearly constant.
Collapse
Affiliation(s)
- Thomas Burgener
- Computational Physics, IfB, ETH-Hönggerberg, Schafmattstrasse 6, CH-8093 Zürich, Switzerland.
| | | | | |
Collapse
|
27
|
Straube AV. Small-scale particle advection, manipulation and mixing: beyond the hydrodynamic scale. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2011; 23:184122. [PMID: 21508483 DOI: 10.1088/0953-8984/23/18/184122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
In this paper we discuss the problems of particle advection, manipulation and mixing at small scales. We start by considering reaction-advection-diffusion systems with the focus on mixing. We show how mixing advection affects the processes of reaction-diffusion and discuss mixing-induced instabilities. Further, we consider the problem of particle manipulation and discuss collective effects in systems comprising solid and compressible particles. We particularly discuss mechanisms of particle entrapment, the role of compressibility in the dynamics of bubbly liquids and nonequilibrium colloidal explosion. Finally, we address two issues related to the problem of wetting. First, we study the role of contact line motion for a sessile droplet (or a bubble) on an oscillating substrate. Second, we discuss an instability of a thin film leading to the formation of a fractal structure of droplets.
Collapse
Affiliation(s)
- Arthur V Straube
- Department of Physics, Humboldt University of Berlin, Berlin, Germany.
| |
Collapse
|
28
|
Brandenburg A, Haugen NEL, Babkovskaia N. Turbulent front speed in the Fisher equation: dependence on Damköhler number. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:016304. [PMID: 21405771 DOI: 10.1103/physreve.83.016304] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2010] [Indexed: 05/30/2023]
Abstract
Direct numerical simulations and mean-field theory are used to model reactive front propagation in a turbulent medium. In the mean-field approach, memory effects of turbulent diffusion are taken into account to estimate the front speed in cases in which the Damköhler number is large. This effect is found to saturate the front speed to values comparable with the speed of the turbulent motions. By comparing with direct numerical simulations, it is found that the effective correlation time is much shorter than for nonreacting flows. The nonlinearity of the reaction term is found to make the front speed slightly faster.
Collapse
Affiliation(s)
- Axel Brandenburg
- NORDITA, AlbaNova University Center, Roslagstullsbacken 23, SE-10691 Stockholm, Sweden.
| | | | | |
Collapse
|