Rare-event sampling applied to the simulation of extreme mechanical eorts exerted by a turbulent ow on a blu body

This study evaluates the relevance of rare-event sampling techniques to accelerate the simulation of extreme mechanical eorts exerted by a turbulent ow impinging onto a blu body. The main idea is to replace a long simulation by a set of much shorter ones, running in parallel, with dynamics that are replicated or pruned in order to sample large-amplitude events more frequently. Such techniques have been shown to be ecient for a wide range of problems in statistical physics, computer science, biochemistry, enabling the simulation of rare events otherwise out of reach by direct sampling. This work is the rst application to uid-structure interaction problems. The drag experienced by a squared obstacle placed in a turbulent ow (in two dimensions) is taken as a representative case study to investigate the performance of two major rare-event sampling algorithms, namely the Adaptive Multilevel Splitting (AMS) and the Giardina-Kurchan-Tailleur-Lecomte (GKTL) algorithms. Practical evidence is given that the fast sweeping-time of uid structures past the obstacle has a drastic inuence on the eciency of these two algorithms. While it is shown that the AMS algorithm does not yield signicant run-time savings, the GKTL algorithm appears to be ecient to sample extreme uctuations of the time-averaged drag and estimate related statistics such as return times. Beyond the study of applicability of rare-event sampling techniques to a uid-mechanical problem, this work also includes a detailed phenomenological description of extreme-drag events of a turbulent ow on a blu body.

[1]  Luca Biferale,et al.  Lattice Monte Carlo methods for systems far from equilibrium , 2013, 1311.4386.

[2]  T. Nemoto,et al.  Finite-time and finite-size scalings in the evaluation of large-deviation functions: Analytical study using a birth-death process. , 2016, Physical review. E.

[3]  F. Bouchet,et al.  Computation of extreme heat waves in climate models using a large deviation algorithm , 2017, Proceedings of the National Academy of Sciences.

[4]  Alessandro De Rosis,et al.  Non-orthogonal central moments relaxing to a discrete equilibrium: A D2Q9 lattice Boltzmann model , 2016 .

[5]  Katepalli R Sreenivasan,et al.  Extreme events in computational turbulence , 2015, Proceedings of the National Academy of Sciences.

[6]  P. Bolhuis Kinetic Pathways of β-Hairpin (Un)folding in Explicit Solvent , 2005 .

[7]  Tony Lelièvre,et al.  ADAPTIVE MULTILEVEL SPLITTING IN MOLECULAR DYNAMICS SIMULATIONS. , 2015, ESAIM. Proceedings and surveys.

[8]  G. Huber,et al.  Weighted-ensemble Brownian dynamics simulations for protein association reactions. , 1996, Biophysical journal.

[9]  J P Garrahan,et al.  Dynamical first-order phase transition in kinetically constrained models of glasses. , 2007, Physical review letters.

[10]  R. Kraichnan,et al.  Is there a statistical mechanics of turbulence , 1988 .

[11]  R. Grauer,et al.  The instanton method and its numerical implementation in fluid mechanics , 2015, 1506.08745.

[12]  Migdal,et al.  Instantons in the Burgers equation. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[13]  Michael C. Sukop,et al.  Lattice Boltzmann Modeling: An Introduction for Geoscientists and Engineers , 2005 .

[14]  Eric Vanden-Eijnden,et al.  Experimental Evidence of Hydrodynamic Instantons: The Universal Route to Rogue Waves , 2019, Physical Review X.

[15]  Vivien Lecomte,et al.  A numerical approach to large deviations in continuous time , 2007 .

[16]  Elaborating transition interface sampling methods , 2005, cond-mat/0405116.

[17]  Freddy Bouchet,et al.  Langevin Dynamics, Large Deviations and Instantons for the Quasi-Geostrophic Model and Two-Dimensional Euler Equations , 2014, 1403.0216.

[18]  Paul Glasserman,et al.  Multilevel Splitting for Estimating Rare Event Probabilities , 1999, Oper. Res..

[19]  Rainer Grauer,et al.  Instanton filtering for the stochastic Burgers equation , 2012, 1209.0905.

[20]  Jorge Kurchan,et al.  Probing rare physical trajectories with Lyapunov weighted dynamics , 2007 .

[21]  Diego Donzis,et al.  Short-term forecasts and scaling of intense events in turbulence , 2008, Journal of Fluid Mechanics.

[22]  H. Touchette The large deviation approach to statistical mechanics , 2008, 0804.0327.

[23]  P. Grassberger,et al.  PERM: A Monte Carlo Strategy for Simulating Polymers and Other Things , 1998 .

[24]  Freddy Bouchet,et al.  Computing Transition Rates for the 1-D Stochastic Ginzburg–Landau–Allen–Cahn Equation for Finite-Amplitude Noise with a Rare Event Algorithm , 2015, 1507.05577.

[25]  Eric Simonnet,et al.  Statistical behaviour of adaptive multilevel splitting algorithms in simple models , 2014, J. Comput. Phys..

[26]  Arnaud Guyader,et al.  Adaptive multilevel splitting: Historical perspective and recent results. , 2019, Chaos.

[27]  Klaus Schulten,et al.  Adaptive Multilevel Splitting Method for Molecular Dynamics Calculation of Benzamidine-Trypsin Dissociation Time. , 2016, Journal of chemical theory and computation.

[28]  J. Rolland Extremely rare collapse and build-up of turbulence in stochastic models of transitional wall flows. , 2018, Physical review. E.

[29]  Freddy Bouchet,et al.  Computation of rare transitions in the barotropic quasi-geostrophic equations , 2014, 1409.3219.

[30]  Robert L. Jack,et al.  Efficient characterisation of large deviations using population dynamics , 2017, 1711.05095.

[31]  P. Moral Feynman-Kac Formulae: Genealogical and Interacting Particle Systems with Applications , 2004 .

[32]  Daniel M Zuckerman,et al.  Weighted Ensemble Simulation: Review of Methodology, Applications, and Software. , 2017, Annual review of biophysics.

[33]  Instantons in a Lagrangian model of turbulence , 2016, 1608.07332.

[34]  P. Glasserman,et al.  A Look At Multilevel Splitting , 1998 .

[35]  A I Chernykh,et al.  Large negative velocity gradients in Burgers turbulence. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[37]  P. Grassberger Go with the Winners: a General Monte Carlo Strategy , 2002, cond-mat/0201313.

[38]  Marcel Lesieur,et al.  Turbulence in Fluids: Stochastic and Numerical Modelling , 1987 .

[39]  Vivien Lecomte,et al.  Simulating Rare Events in Dynamical Processes , 2011, 1106.4929.

[40]  Orestis Malaspinas,et al.  Straight velocity boundaries in the lattice Boltzmann method. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[41]  J. Araque,et al.  Transition path sampling and forward flux sampling. Applications to biological systems , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[42]  F. Bouchet,et al.  Rare Event Algorithm Links Transitions in Turbulent Flows with Activated Nucleations. , 2018, Physical review letters.

[43]  Eric Vanden-Eijnden,et al.  Rogue waves and large deviations in deep sea , 2017, Proceedings of the National Academy of Sciences.

[44]  Thibault Lestang,et al.  Computing return times or return periods with rare event algorithms , 2017, 1711.08428.

[45]  D. Mesterházy,et al.  Anomalous scaling in the random-force-driven Burgers' equation: a Monte Carlo study , 2011, 1111.2304.

[46]  J. Kurchan,et al.  Large deviations of Lyapunov exponents , 2013, 1302.6254.

[47]  José Villén-Altamirano,et al.  RESTART: a straightforward method for fast simulation of rare events , 1994, Proceedings of Winter Simulation Conference.

[48]  F. Cérou,et al.  Adaptive Multilevel Splitting for Rare Event Analysis , 2007 .

[49]  Umesh V. Vazirani,et al.  "Go with the winners" algorithms , 1994, Proceedings 35th Annual Symposium on Foundations of Computer Science.

[50]  Freddy Bouchet,et al.  Rare event computation in deterministic chaotic systems using genealogical particle analysis , 2015, 1511.02703.