1
|
Bandak D, Goldenfeld N, Mailybaev AA, Eyink G. Dissipation-range fluid turbulence and thermal noise. Phys Rev E 2022; 105:065113. [PMID: 35854607 DOI: 10.1103/physreve.105.065113] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2021] [Accepted: 04/26/2022] [Indexed: 11/07/2022]
Abstract
We revisit the issue of whether thermal fluctuations are relevant for incompressible fluid turbulence and estimate the scale at which they become important. As anticipated by Betchov in a prescient series of works more than six decades ago, this scale is about equal to the Kolmogorov length, even though that is several orders of magnitude above the mean free path. This result implies that the deterministic version of the incompressible Navier-Stokes equation is inadequate to describe the dissipation range of turbulence in molecular fluids. Within this range, the fluctuating hydrodynamics equation of Landau and Lifschitz is more appropriate. In particular, our analysis implies that both the exponentially decaying energy spectrum and the far-dissipation-range intermittency predicted by Kraichnan for deterministic Navier-Stokes will be generally replaced by Gaussian thermal equipartition at scales just below the Kolmogorov length. Stochastic shell model simulations at high Reynolds numbers verify our theoretical predictions and reveal furthermore that inertial-range intermittency can propagate deep into the dissipation range, leading to large fluctuations in the equipartition length scale. We explain the failure of previous scaling arguments for the validity of deterministic Navier-Stokes equations at any Reynolds number and we provide a mathematical interpretation and physical justification of the fluctuating Navier-Stokes equation as an "effective field theory" valid below some high-wave-number cutoff Λ, rather than as a continuum stochastic partial differential equation. At Reynolds number around a million, comparable to that in Earth's atmospheric boundary layer, the strongest turbulent excitations observed in our simulation penetrate down to a length scale of about eight microns, still two orders of magnitude greater than the mean free path of air. However, for longer observation times or for higher Reynolds numbers, more extreme turbulent events could lead to a local breakdown of fluctuating hydrodynamics.
Collapse
Affiliation(s)
- Dmytro Bandak
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Nigel Goldenfeld
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Alexei A Mailybaev
- Instituto de Matemática Pura e Aplicada-IMPA, Rio de Janeiro, 22460-320, Brazil
| | - Gregory Eyink
- Department of Applied Mathematics & Statistics, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
2
|
Ebener L, Margazoglou G, Friedrich J, Biferale L, Grauer R. Instanton based importance sampling for rare events in stochastic PDEs. CHAOS (WOODBURY, N.Y.) 2019; 29:063102. [PMID: 31266309 DOI: 10.1063/1.5085119] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2018] [Accepted: 05/10/2019] [Indexed: 06/09/2023]
Abstract
We present a new method for sampling rare and large fluctuations in a nonequilibrium system governed by a stochastic partial differential equation (SPDE) with additive forcing. To this end, we deploy the so-called instanton formalism that corresponds to a saddle-point approximation of the action in the path integral formulation of the underlying SPDE. The crucial step in our approach is the formulation of an alternative SPDE that incorporates knowledge of the instanton solution such that we are able to constrain the dynamical evolutions around extreme flow configurations only. Finally, a reweighting procedure based on the Girsanov theorem is applied to recover the full distribution function of the original system. The entire procedure is demonstrated on the example of the one-dimensional Burgers equation. Furthermore, we compare our method to conventional direct numerical simulations as well as to Hybrid Monte Carlo methods. It will be shown that the instanton-based sampling method outperforms both approaches and allows for an accurate quantification of the whole probability density function of velocity gradients from the core to the very far tails.
Collapse
Affiliation(s)
- Lasse Ebener
- Institut für Theoretische Physik I, Ruhr-Universität Bochum, Universitätsstraße 150, 44780 Bochum, Germany
| | - Georgios Margazoglou
- Department of Physics and INFN, University of Rome "Tor Vergata," Via della Ricerca Scientifica 1, I-00133 Roma, Italy
| | - Jan Friedrich
- Institut für Theoretische Physik I, Ruhr-Universität Bochum, Universitätsstraße 150, 44780 Bochum, Germany
| | - Luca Biferale
- Department of Physics and INFN, University of Rome "Tor Vergata," Via della Ricerca Scientifica 1, I-00133 Roma, Italy
| | - Rainer Grauer
- Institut für Theoretische Physik I, Ruhr-Universität Bochum, Universitätsstraße 150, 44780 Bochum, Germany
| |
Collapse
|
3
|
Debue P, Shukla V, Kuzzay D, Faranda D, Saw EW, Daviaud F, Dubrulle B. Dissipation, intermittency, and singularities in incompressible turbulent flows. Phys Rev E 2018; 97:053101. [PMID: 29906866 DOI: 10.1103/physreve.97.053101] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Indexed: 11/07/2022]
Abstract
We examine the connection between the singularities or quasisingularities in the solutions of the incompressible Navier-Stokes equation (INSE) and the local energy transfer and dissipation, in order to explore in detail how the former contributes to the phenomenon of intermittency. We do so by analyzing the velocity fields (a) measured in the experiments on the turbulent von Kármán swirling flow at high Reynolds numbers and (b) obtained from the direct numerical simulations of the INSE at a moderate resolution. To compute the local interscale energy transfer and viscous dissipation in experimental and supporting numerical data, we use the weak solution formulation generalization of the Kármán-Howarth-Monin equation. In the presence of a singularity in the velocity field, this formulation yields a nonzero dissipation (inertial dissipation) in the limit of an infinite resolution. Moreover, at finite resolutions, it provides an expression for local interscale energy transfers down to the scale where the energy is dissipated by viscosity. In the presence of a quasisingularity that is regularized by viscosity, the formulation provides the contribution to the viscous dissipation due to the presence of the quasisingularity. Therefore, our formulation provides a concrete support to the general multifractal description of the intermittency. We present the maps and statistics of the interscale energy transfer and show that the extreme events of this transfer govern the intermittency corrections and are compatible with a refined similarity hypothesis based on this transfer. We characterize the probability distribution functions of these extreme events via generalized Pareto distribution analysis and find that the widths of the tails are compatible with a similarity of the second kind. Finally, we make a connection between the topological and the statistical properties of the extreme events of the interscale energy transfer field and its multifractal properties.
Collapse
Affiliation(s)
- P Debue
- DSM/IRAMIS/SPEC, CNRS UMR 3680, CEA, Université Paris-Saclay, 91190 Gif sur Yvette, France
| | - V Shukla
- DSM/IRAMIS/SPEC, CNRS UMR 3680, CEA, Université Paris-Saclay, 91190 Gif sur Yvette, France
| | - D Kuzzay
- DSM/IRAMIS/SPEC, CNRS UMR 3680, CEA, Université Paris-Saclay, 91190 Gif sur Yvette, France.,LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92195 Meudon, France
| | - D Faranda
- DSM/LSCE, CNRS UMR 8212, CEA, Université Paris-Saclay, 91190 Gif sur Yvette, France.,London Mathematical Laboratory, 14 Buckingham Street, London WC2N 6DF, United Kingdom
| | - E-W Saw
- DSM/IRAMIS/SPEC, CNRS UMR 3680, CEA, Université Paris-Saclay, 91190 Gif sur Yvette, France
| | - F Daviaud
- DSM/IRAMIS/SPEC, CNRS UMR 3680, CEA, Université Paris-Saclay, 91190 Gif sur Yvette, France
| | - B Dubrulle
- DSM/IRAMIS/SPEC, CNRS UMR 3680, CEA, Université Paris-Saclay, 91190 Gif sur Yvette, France
| |
Collapse
|
4
|
Gürcan ÖD. Nested polyhedra model of isotropic magnetohydrodynamic turbulence. Phys Rev E 2018; 97:063111. [PMID: 30011494 DOI: 10.1103/physreve.97.063111] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2017] [Indexed: 06/08/2023]
Abstract
A nested polyhedra model has been developed for magnetohydrodynamic turbulence. Driving only the velocity field at large scales with random, divergence-free forcing results in a clear, stationary k^{-5/3} spectrum for both kinetic and magnetic energies. Since the model naturally effaces disparate scale interactions, does not have a guide field, and avoids injecting any sign of helicity by random forcing, the resulting three-dimensional k spectrum is statistically isotropic. The strengths and weaknesses of the model are demonstrated by considering large or small magnetic Prandtl numbers. It was also observed that the timescale for the equipartition offset with those of the smallest scales shows a k^{-1/2} scaling.
Collapse
Affiliation(s)
- Ö D Gürcan
- CNRS, Laboratoire de Physique des Plasmas, Ecole Polytechnique, 91128 Palaiseau, France and Sorbonne Universités, UPMC Univ Paris 06, Paris, France
| |
Collapse
|
5
|
De Pietro M, Biferale L, Mailybaev AA. Inverse energy cascade in nonlocal helical shell models of turbulence. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:043021. [PMID: 26565346 DOI: 10.1103/physreve.92.043021] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2015] [Indexed: 06/05/2023]
Abstract
Following the exact decomposition in eigenstates of helicity for the Navier-Stokes equations in Fourier space [F. Waleffe, Phys. Fluids A 4, 350 (1992)], we introduce a modified version of helical shell models for turbulence with nonlocal triadic interactions. By using both an analytical argument and numerical simulation, we show that there exists a class of models, with a specific helical structure, that exhibits a statistically stable inverse energy cascade, in close analogy with that predicted for the Navier-Stokes equations restricted to the same helical interactions. We further support the idea that turbulent energy transfer is the result of a strong entanglement among triads possessing different transfer properties.
Collapse
Affiliation(s)
- Massimo De Pietro
- Dip. di Fisica and INFN, Università "Tor Vergata," Via della Ricerca Scientifica 1, I-00133 Roma, Italy
| | - Luca Biferale
- Dip. di Fisica and INFN, Università "Tor Vergata," Via della Ricerca Scientifica 1, I-00133 Roma, Italy
| | - Alexei A Mailybaev
- Instituto Nacional de Matemática Pura e Aplicada-IMPA, Est. Dona Castorina 110, Rio de Janeiro 22460-320 Brazil
| |
Collapse
|