Source author record

Michael Wilczek

Michael Wilczek appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

A minimal phase-coupling model for intermittency in turbulent systems

Turbulent systems exhibit a remarkable multi-scale complexity, in which spatial structures induce scale-dependent statistics with strong departures from Gaussianity. In Fourier space, this is reflected by pronounced phase synchronization. A quantitative relation between real-space structure, statistics, and phase synchronization is currently missing. Here, we address this problem in the framework of a minimal phase-coupling model, which enables a detailed investigation by means of dynamical systems theory and multi-scale high-resolution simulations. We identify the spectral power-law steepness, which controls the phase coupling, as the control parameter for tuning the non-Gaussian properties of the system. Whereas both very steep and very shallow spectra exhibit close-to-Gaussian statistics, the strongest departures are observed for intermediate slopes comparable to the ones in hydrodynamic and Burgers turbulence. We show that the non-Gaussian regime of the model coincides with a collapse of the dynamical system to a lower-dimensional attractor and the emergence of phase synchronization, thereby establishing a dynamical-systems perspective on turbulent intermittency.

preprint2022arXiv

An Efficient Particle Tracking Algorithm for Large-Scale Parallel Pseudo-Spectral Simulations of Turbulence

Particle tracking in large-scale numerical simulations of turbulent flows presents one of the major bottlenecks in parallel performance and scaling efficiency. Here, we describe a particle tracking algorithm for large-scale parallel pseudo-spectral simulations of turbulence which scales well up to billions of tracer particles on modern high-performance computing architectures. We summarize the standard parallel methods used to solve the fluid equations in our hybrid MPI/OpenMP implementation. As the main focus, we describe the implementation of the particle tracking algorithm and document its computational performance. To address the extensive inter-process communication required by particle tracking, we introduce a task-based approach to overlap point-to-point communications with computations, thereby enabling improved resource utilization. We characterize the computational cost as a function of the number of particles tracked and compare it with the flow field computation, showing that the cost of particle tracking is very small for typical applications.

preprint2022arXiv

Optimal Turbulent Transport in Microswimmer Suspensions

Microswimmer suspensions self-organize into complex spatio-temporal flow patterns, including vortex lattices and mesoscale turbulence. Here we explore the consequences for the motion of passive tracers, based on a continuum model for the microswimmer velocity field. We observe two qualitatively different regimes distinguished via the dimensionless Kubo number $K$. At advection strengths right above the transition to turbulence, the flow field evolves very slowly ($K \gg 1$) and the spatial vortex structures lead to dominant trapping effects. In contrast, deep in the turbulent state, much faster dynamics ($K \ll 1$) consistent with the so-called sweeping hypothesis leads to transport properties completely determined by the temporal correlations. In between ($K \approx 1$), we observe a regime of optimal transport, signaled by a maximum of the diffusion coefficient.

preprint2022arXiv

The statistical geometry of material loops in turbulence

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.

preprint2021arXiv

The role of advective inertia in active nematic turbulence

Suspensions of active agents with nematic interactions exhibit complex spatio-temporal dynamics such as mesoscale turbulence. Since the Reynolds number of microscopic flows is very small on the scale of individual agents, inertial effects are typically excluded in continuum theories of active nematic turbulence. Whether active stresses can collectively excite inertial flows is currently unclear. To address this question, we investigate a two-dimensional continuum theory for active nematic turbulence. In particular, we compare mesoscale turbulence with and without the effects of advective inertia. We find that inertial effects can influence the flow already close to the onset of the turbulent state and, moreover, give rise to large-scale fluid motion for strong active driving. A detailed analysis of the kinetic energy budget reveals an energy transfer to large scales mediated by inertial advection. While this transfer is small in comparison to energy injection and dissipation, its effects accumulate over time. The inclusion of friction, which is typically present in experiments, can compensate for this effect. The findings suggest that the inclusion of inertia and friction may be necessary for dynamically consistent theories of active nematic turbulence.

preprint2020arXiv

Resolved energy budget of superstructures in Rayleigh-Bénard convection

Turbulent superstructures, i.e. large-scale flow structures in turbulent flows, play a crucial role in many geo- and astrophysical settings. In turbulent Rayleigh-Bénard convection, for example, horizontally extended coherent large-scale convection rolls emerge. Currently, a detailed understanding of the interplay of small-scale turbulent fluctuations and large-scale coherent structures is missing. Here, we investigate the resolved kinetic energy and temperature variance budgets by applying a filtering approach to direct numerical simulations of Rayleigh-Bénard convection at high aspect ratio. In particular, we focus on the energy transfer rate between large-scale flow structures and small-scale fluctuations. We show that the small scales primarily act as a dissipation for the superstructures. However, we find that the height-dependent energy transfer rate has a complex structure with distinct bulk and boundary layer features. Additionally, we observe that the heat transfer between scales mainly occurs close to the thermal boundary layer. Our results clarify the interplay of superstructures and turbulent fluctuations and may help to guide the development of an effective description of large-scale flow features in terms of reduced-order models.

preprint2020arXiv

The challenges of containing SARS-CoV-2 via test-trace-and-isolate

Without a cure, vaccine, or proven long-term immunity against SARS-CoV-2, test-trace-and-isolate (TTI) strategies present a promising tool to contain its spread. For any TTI strategy, however, mitigation is challenged by pre- and asymptomatic transmission, TTI-avoiders, and undetected spreaders, who strongly contribute to hidden infection chains. Here, we studied a semi-analytical model and identified two tipping points between controlled and uncontrolled spread: (1) the behavior-driven reproduction number of the hidden chains becomes too large to be compensated by the TTI capabilities, and (2) the number of new infections exceeds the tracing capacity. Both trigger a self-accelerating spread. We investigated how these tipping points depend on challenges like limited cooperation, missing contacts, and imperfect isolation. Our model results suggest that TTI alone is insufficient to contain an otherwise unhindered spread of SARS-CoV-2, implying that complementary measures like social distancing and improved hygiene remain necessary.

preprint2015arXiv

Spatio-temporal spectra in the logarithmic layer of wall turbulence: large-eddy simulations and simple models

Motivated by the need to characterize the spatio-temporal structure of turbulence in wall-bounded flows, we study wavenumber-frequency spectra of the streamwise velocity component based on large-eddy simulation (LES) data. The LES data are used to measure spectra as a function of the two wall-parallel wavenumbers and the frequency in the equilibrium (logarithmic) layer. We then reformulate one of the simplest models that is able to reproduce the observations: the random sweeping model with a Gaussian large-scale fluctuating velocity and with additional mean flow. Comparison with LES data shows that the model captures the observed temporal decorrelation, which is related to the Doppler broadening of frequencies. We furthermore introduce a parameterization for the entire wavenumber-frequency spectrum $E_{11}(k_1,k_2,ω;z)$, where $k_1$, $k_2$ are the streamwise and spanwise wavenumbers, $ω$ is the frequency and $z$ is the distance to the wall. The results are found to be in good agreement with LES data.

preprint2015arXiv

Turbulent Rayleigh-Bénard convection described by projected dynamics in phase space

Rayleigh-Bénard convection, i.e. the flow of a fluid between two parallel plates that is driven by a temperature gradient, is an idealised setup to study thermal convection. Of special interest are the statistics of the turbulent temperature field, which we are investigating and comparing for three different geometries, namely convection with periodic horizontal boundary conditions in three and two dimensions as well as convection in a cylindrical vessel, in order to work out similarities and differences. To this end, we derive an exact evolution equation for the temperature probability density function (PDF). Unclosed terms are expressed as conditional averages of velocities and heat diffusion, which are estimated from direct numerical simulations. This framework lets us identify the average behaviour of a fluid particle by revealing the mean evolution of fluid of different temperatures in different parts of the convection cell. We connect the statistics to the dynamics of Rayleigh-Bénard convection, giving deeper insights into the temperature statistics and transport mechanisms. We find that the average behaviour is described by closed cycles in phase space that reconstruct the typical Rayleigh-Bénard cycle of fluid heating up at the bottom, rising up to the top plate, cooling down and falling down again. The detailed behaviour shows subtle differences between the three cases.

preprint2014arXiv

Large-eddy simulation study of the logarithmic law for second and higher-order moments in turbulent wall-bounded flow

The logarithmic law for the mean velocity in turbulent boundary layers has long provided a valuable and robust reference for comparison with theories, models, and large-eddy simulations (LES) of wall-bounded turbulence. More recently, analysis of high-Reynolds number experimental boundary layer data has shown that also the variance and higher-order moments of the streamwise velocity fluctuations $u'^{+}$ display logarithmic laws. Such experimental observations motivate the question whether LES can accurately reproduce the variance and the higher-order moments, in particular their logarithmic dependency on distance to the wall. In this study we perform LES of very high Reynolds number wall-modeled channel flow and focus on profiles of variance and higher-order moments of the streamwise velocity fluctuations. In agreement with the experimental data, we observe an approximately logarithmic law for the variance in the LES, with a `Townsend-Perry' constant of $A_1\approx 1.25$. The LES also yields approximate logarithmic laws for the higher-order moments of the streamwise velocity. Good agreement is found between $A_p$, the generalized `Townsend-Perry' constants for moments of order $2p$, from experiments and simulations. Both are indicative of sub-Gaussian behavior of the streamwise velocity fluctuations. The near-wall behavior of the variance, the ranges of validity of the logarithmic law and in particular possible dependencies on characteristic length scales such as the roughness scale $z_0$, the LES grid scale $Δ$, and sub-grid scale (SGS) mixing length $C_sΔ$ are examined. We also present LES results on moments of spanwise and wall-normal fluctuations of velocity.

preprint2014arXiv

Pressure Hessian and viscous contributions to velocity gradient statistics based on Gaussian random fields

Understanding the non-local pressure contributions and viscous effects on the small-scale statistics remains one of the central challenges in the study of homogeneous isotropic turbulence. Here we address this issue by studying the impact of the pressure Hessian as well as viscous diffusion on the statistics of the velocity gradient tensor in the framework of an exact statistical evolution equation. This evolution equation shares similarities with earlier phenomenological models for the Lagrangian velocity gradient tensor evolution, yet constitutes the starting point for a systematic study of the unclosed pressure Hessian and viscous diffusion terms. Based on the assumption of incompressible Gaussian velocity fields, closed expressions are obtained as the results of an evaluation of the characteristic functionals. The benefits and shortcomings of this Gaussian closure are discussed, and a generalization is proposed based on results from direct numerical simulations. This enhanced Gaussian closure yields, for example, insights on how the pressure Hessian prevents the finite-time singularity induced by the local self-amplification and how its interaction with viscous effects leads to the characteristic strain skewness phenomenon.

preprint2013arXiv

Nonuniversal power-law spectra in turbulent systems

Turbulence is generally associated with universal power-law spectra in scale ranges without significant drive or damping. Although many examples of turbulent systems do not exhibit such an inertial range, power-law spectra may still be observed. As a simple model for such situations, a modified version of the Kuramoto-Sivashinsky equation is studied. By means of semi-analytical and numerical studies, one finds power laws with nonuniversal exponents in the spectral range for which the ratio of nonlinear and linear time scales is (roughly) scale-independent.

preprint2012arXiv

Conditional vorticity budget of coherent and incoherent flow contributions in fully developed homogeneous isotropic turbulence

We investigate the conditional vorticity budget of fully developed three-dimensional homogeneous isotropic turbulence with respect to coherent and incoherent flow contributions. The Coherent Vorticity Extraction based on orthogonal wavelets allows to decompose the vorticity field into coherent and incoherent contributions, of which the latter are noise-like. The impact of the vortex structures observed in fully developed turbulence on statistical balance equations is quantified considering the conditional vorticity budget. The connection between the basic structures present in the flow and their statistical implications is thereby assessed. The results are compared to those obtained for large- and small-scale contributions using a Fourier decomposition, which reveals pronounced differences.

preprint2011arXiv

On the velocity distribution in homogeneous isotropic turbulence: correlations and deviations from Gaussianity

We investigate the single-point probability density function of the velocity in three-dimensional stationary and decaying homogeneous isotropic turbulence. To this end we apply the statistical framework of the Lundgren-Monin-Novikov hierarchy combined with conditional averaging, identifying the quantities that determine the shape of the probability density function. In this framework the conditional averages of the rate of energy dissipation, the velocity diffusion and the pressure gradient with respect to velocity play a key role. Direct numerical simulations of the Navier-Stokes equation are used to complement the theoretical results and assess deviations from Gaussianity.

preprint2011arXiv

Theory for the single-point velocity statistics of fully developed turbulence

We investigate the single-point velocity probability density function (PDF) in three-dimensional fully developed homogeneous isotropic turbulence within the framework of PDF equations focussing on deviations from Gaussianity. A joint analytical and numerical analysis shows that these deviations may be quantified studying correlations of dynamical quantities like pressure gradient, external forcing and energy dissipation with the velocity. A stationary solution for the PDF equation in terms of these quantities is presented, and the theory is validated with the help of direct numerical simulations indicating sub-Gaussian tails of the PDF.