skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Bayesian conditional diffusion models for versatile spatiotemporal turbulence generation
Turbulent flows, characterized by their chaotic and stochastic nature, have historically presented formidable challenges to predictive computational modeling. Traditional eddy-resolved numerical simulations often require vast computational resources, making them impractical or infeasible for numerous engineering applications. As an alternative, deep learning-based surrogate models have emerged, offering data-drive solutions. However, these are typically constructed within deterministic settings, leading to shortfall in capturing the innate chaotic and stochastic behaviors of turbulent dynamics. In this study, we introduce a novel generative framework grounded in probabilistic diffusion models for versatile generation of spatiotemporal turbulence under various conditions. Our method unifies both unconditional and conditional sampling strategies within a Bayesian framework, which can accommodate diverse conditioning scenarios, including those with a direct differentiable link between specified conditions and generated unsteady flow outcomes, as well as scenarios lacking such explicit correlations. A notable feature of our approach is the method proposed for long-span flow sequence generation, which is based on autoregressive gradient-based conditional sampling, eliminating the need for cumbersome retraining processes. We evaluate and showcase the versatile turbulence generation capability of our framework through a suite of numerical experiments, including: (1) the synthesis of Large Eddy Simulations (LES) simulated instantaneous flow sequences from unsteady Reynolds-Averaged Navier–Stokes (URANS) inputs; (2) holistic generation of inhomogeneous, anisotropic wall-bounded turbulence, whether from given initial conditions, prescribed turbulence statistics, or entirely from scratch; (3) super-resolved generation of high-speed turbulent boundary layer flows from low-resolution data across a range of input resolutions. Collectively, our numerical experiments highlight the merit and transformative potential of the proposed methods, making a significant advance in the field of turbulence generation.  more » « less
Award ID(s):
2047127
PAR ID:
10519148
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Computer Methods in Applied Mechanics and Engineering
Volume:
427
Issue:
C
ISSN:
0045-7825
Page Range / eLocation ID:
117023
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The atmospheric boundary layer (ABL) is a highly turbulent geophysical flow, which has chaotic and often too complex dynamics to unravel from limited data. Characterizing coherent turbulence structures in complex ABL flows under various atmospheric regimes is not systematically well established yet. This study aims to bridge this gap using large eddy simulations (LESs), Koopman theory, and unsupervised classification techniques. To this end, eight LESs of different convective, neutral, and unsteady ABLs are conducted. As the ratio of buoyancy to shear production increases, the turbulence structures change from roll vortices to convective cells. The quadrant analysis indicated that as this ratio increases, the sweep and ejection events decrease, and inward/outward interactions increase. The Koopman mode decomposition (KMD) is then used to characterize their turbulence structures. Our results showed that KMD can reveal non-trivial modes of highly turbulent ABL flows (e.g., transverse to the mean flow direction) and can reconstruct the primary dynamics of ABLs even under unsteady conditions with only ∼5% of the modes. We attributed the detected modes to the imposed pressure gradient (shear), Coriolis (inertial oscillations), and buoyancy (convection) forces by conducting novel timescale and quadrant analyses. We then applied the convolutional neural network combined with the K-means clustering to group the Koopman modes. This approach is displacement and rotation invariant, which allows efficiently reducing the number of modes that describe the overall ABL dynamics. Our results provide new insights into the dynamics of ABLs and present a systematic data-driven method to characterize their complex spatiotemporal patterns. 
    more » « less
  2. Abstract Material elements – which are lines, surfaces, or volumes behaving as passive, non-diffusive markers – provide an inherently geometric window into the intricate dynamics of chaotic flows. Their stretching and folding dynamics has immediate implications for mixing in the oceans or the atmosphere, as well as the emergence of self-sustained dynamos in astrophysical settings. Here, we uncover robust statistical properties of an ensemble of material loops in a turbulent environment. Our approach combines high-resolution direct numerical simulations of Navier-Stokes turbulence, stochastic models, and dynamical systems techniques to reveal predictable, universal features of these complex objects. We show that the loop curvature statistics become stationary through a dynamical formation process of high-curvature folds, leading to distributions with power-law tails whose exponents are determined by the large-deviations statistics of finite-time Lyapunov exponents of the flow. This prediction applies to advected material lines in a broad range of chaotic flows. To complement this dynamical picture, we confirm our theory in the analytically tractable Kraichnan model with an exact Fokker-Planck approach. 
    more » « less
  3. Low-fidelity engineering wake models are often combined with linear superposition laws to predict wake velocities across wind farms under steady atmospheric conditions. While convenient for wind farm planning and long-term performance evaluation, such models are unable to capture the time-varying nature of the waked velocity field, as they are agnostic to the complex aerodynamic interactions among wind turbines and the effects of atmospheric boundary layer turbulence. To account for such effects while remaining amenable to conventional system-theoretic tools for flow estimation and control, we propose a new class of data-enhanced physics-based models for the dynamics of wind farm flow fluctuations. Our approach relies on the predictive capability of the stochastically forced linearized Navier–Stokes equations around static base flow profiles provided by conventional engineering wake models. We identify the stochastic forcing into the linearized dynamics via convex optimization to ensure statistical consistency with higher-fidelity models or experimental measurements while preserving model parsimony. We demonstrate the utility of our approach in completing the statistical signature of wake turbulence in accordance with large-eddy simulations of turbulent flow over a cascade of yawed wind turbines. Our numerical experiments provide insight into the significance of spatially distributed field measurements in recovering the statistical signature of wind farm turbulence and training stochastic linear models for short-term wind forecasting. 
    more » « less
  4. Accurate simulation of turbulent flows is of crucial importance in many branches of science and engineering. Direct numerical simulation (DNS) provides the highest fidelity means of capturing all intricate physics of turbulent transport. However, the method is computationally expensive because of the wide range of turbulence scales that must be accounted for in such simulations. Large eddy simulation (LES) provides an alternative. In such simulations, the large scales of the flow are resolved, and the effects of small scales are modelled. Reconstruction of the DNS field from the low-resolution LES is needed for a wide variety of applications. Thus the construction of super-resolution methodologies that can provide this reconstruction has become an area of active research. In this work, a new physics-guided neural network is developed for such a reconstruction. The method leverages the partial differential equation that underlies the flow dynamics in the design of spatio-temporal model architecture. A degradation-based refinement method is also developed to enforce physical constraints and to further reduce the accumulated reconstruction errors over long periods. Detailed DNS data on two turbulent flow configurations are used to assess the performance of the model. 
    more » « less
  5. The deformation of elementary fluid volumes by velocity gradients is a key process for scalar mixing, chemical reactions and biological processes in flows. Whilst fluid deformation in unsteady, turbulent flow has gained much attention over the past half-century, deformation in steady random flows with complex structure – such as flow through heterogeneous porous media – has received significantly less attention. In contrast to turbulent flow, the steady nature of these flows constrains fluid deformation to be anisotropic with respect to the fluid velocity, with significant implications for e.g. longitudinal and transverse mixing and dispersion. In this study we derive an ab initio coupled continuous-time random walk (CTRW) model of fluid deformation in random steady three-dimensional flow that is based upon a streamline coordinate transform which renders the velocity gradient and fluid deformation tensors upper triangular. We apply this coupled CTRW model to several model flows and find that these exhibit a remarkably simple deformation structure in the streamline coordinate frame, facilitating solution of the stochastic deformation tensor components. These results show that the evolution of longitudinal and transverse fluid deformation for chaotic flows is governed by both the Lyapunov exponent and power-law exponent of the velocity probability distribution function at small velocities, whereas algebraic deformation in non-chaotic flows arises from the intermittency of shear events following similar dynamics as that for steady two-dimensional flow. 
    more » « less