Direct numerical simulations of capillary wave turbulence

This work presents Direct Numerical Simulations of capillary wave turbulence solving the full 3D Navier Stokes equations of a two-phase flow. When the interface is locally forced at large scales, a statistical stationary state appears after few forcing periods. Smaller wave scales are generated by nonlinear interactions, and the wave height spectrum is found to obey a power law in both wave number and frequency in good agreement with weak turbulence theory. By estimating the mean energy flux from the dissipated power, the Kolmogorov-Zakharov constant is evaluated and found to be compatible with the exact theoretical value. The time scale separation between linear, nonlinear interaction and dissipative times is also observed. These numerical results confirm the validity of weak turbulence approach to quantify out-of equilibrium wave statistics.Comment: Physical Review Letters (2014) in pres

Optimal Taylor-Couette flow: direct numerical simulations

We numerically simulate turbulent Taylor-Couette flow for independently rotating inner and outer cylinders, focusing on the analogy with turbulent Rayleigh-B\'enard flow. Reynolds numbers of $Re_i=8\cdot10^3$ and $Re_o=\pm4\cdot10^3$ of the inner and outer cylinders, respectively, are reached, corresponding to Taylor numbers Ta up to $10^8$. Effective scaling laws for the torque and other system responses are found. Recent experiments with the Twente turbulent Taylor-Couette ($T^3C$) setup and with a similar facility in Maryland at very high Reynolds numbers have revealed an optimum transport at a certain non-zero rotation rate ratio $a = -\omega_o / \omega_i$ of about $a_{opt}=0.33-0.35$. For large enough $Ta$ in the numerically accessible range we also find such an optimum transport at non-zero counter-rotation. The position of this maximum is found to shift with the driving, reaching a maximum of $a_{opt}=0.15$ for $Ta=2.5\cdot10^7$. An explanation for this shift is elucidated, consistent with the experimental result that $a_{opt}$ becomes approximately independent of the driving strength for large enough Reynolds numbers. We furthermore numerically calculate the angular velocity profiles and visualize the different flow structures for the various regimes. By writing the equations in a frame co-rotating with the outer cylinder a link is found between the local angular velocity profiles and the global transport quantities.Comment: Under consideration for publication in JFM, 31 pages, 25 figure

Direct numerical simulations of aeolian sand ripples

Aeolian sand beds exhibit regular patterns of ripples resulting from the interaction between topography and sediment transport. Their characteristics have been so far related to reptation transport caused by the impacts on the ground of grains entrained by the wind into saltation. By means of direct numerical simulations of grains interacting with a wind flow, we show that the instability turns out to be driven by resonant grain trajectories, whose length is close to a ripple wavelength and whose splash leads to a mass displacement towards the ripple crests. The pattern selection results from a compromise between this destabilizing mechanism and a diffusive downslope transport which stabilizes small wavelengths. The initial wavelength is set by the ratio of the sediment flux and the erosion/deposition rate, a ratio which increases linearly with the wind velocity. We show that this scaling law, in agreement with experiments, originates from an interfacial layer separating the saltation zone from the static sand bed, where momentum transfers are dominated by mid-air collisions. Finally, we provide quantitative support for the use the propagation of these ripples as a proxy for remote measurements of sediment transport.Comment: 21 pages, 12 figure

Direct Numerical Simulations of Electrophoresis of Charged Colloids

We propose a numerical method to simulate electrohydrodynamic phenomena in charged colloidal dispersions. This method enables us to compute the time evolutions of colloidal particles, ions, and host fluids simultaneously by solving Newton, advection-diffusion, and Navier--Stokes equations so that the electrohydrodynamic couplings can be fully taken into account. The electrophoretic mobilities of charged spherical particles are calculated in several situations. The comparisons with approximation theories show quantitative agreements for dilute dispersions without any empirical parameters, however, our simulation predicts notable deviations in the case of dense dispersions.Comment: 4pages, 3figures, to appear in Phys. Rev. Let

Direct numerical simulations of vortex rings at ReΓ = 7500

We present direct numerical simulations of the turbulent decay of vortex rings with ReΓ = 7500. We analyse the vortex dynamics during the nonlinear stage of the instability along with the structure of the vortex wake during the turbulent stage. These simulations enable the quantification of vorticity dynamics and their correlation with structures from dye visualization and the observations of circulation decay that have been reported in related experimental works. Movies are available with the online version of the paper

A weakly nonlinear theory for wave-vortex interactions in curved channel flow

A weakly nonlinear theory is developed to study the interaction of Tollmien-Schlichting (TS) waves and Dean vortices in curved channel flow. The predictions obtained from the theory agree well with results obtained from direct numerical simulations of curved channel flow, especially for low amplitude disturbances. Some discrepancies in the results of a previous theory with direct numerical simulations are resolved

Advanced in turbulence physics and modeling by direct numerical simulations

The advent of direct numerical simulations of turbulence has opened avenues for research on turbulence physics and turbulence modeling. Direct numerical simulation provides values for anything that the scientist or modeler would like to know about the flow. An overview of some recent advances in the physical understanding of turbulence and in turbulence modeling obtained through such simulations is presented

High performance Python for direct numerical simulations of turbulent flows

Direct Numerical Simulations (DNS) of the Navier Stokes equations is an invaluable research tool in fluid dynamics. Still, there are few publicly available research codes and, due to the heavy number crunching implied, available codes are usually written in low-level languages such as C/C++ or Fortran. In this paper we describe a pure scientific Python pseudo-spectral DNS code that nearly matches the performance of C++ for thousands of processors and billions of unknowns. We also describe a version optimized through Cython, that is found to match the speed of C++. The solvers are written from scratch in Python, both the mesh, the MPI domain decomposition, and the temporal integrators. The solvers have been verified and benchmarked on the Shaheen supercomputer at the KAUST supercomputing laboratory, and we are able to show very good scaling up to several thousand cores. A very important part of the implementation is the mesh decomposition (we implement both slab and pencil decompositions) and 3D parallel Fast Fourier Transforms (FFT). The mesh decomposition and FFT routines have been implemented in Python using serial FFT routines (either NumPy, pyFFTW or any other serial FFT module), NumPy array manipulations and with MPI communications handled by MPI for Python (mpi4py). We show how we are able to execute a 3D parallel FFT in Python for a slab mesh decomposition using 4 lines of compact Python code, for which the parallel performance on Shaheen is found to be slightly better than similar routines provided through the FFTW library. For a pencil mesh decomposition 7 lines of code is required to execute a transform