118 research outputs found
Waves and vortices in the inverse cascade regime of stratified turbulence with or without rotation
We study the partition of energy between waves and vortices in stratified
turbulence, with or without rotation, for a variety of parameters, focusing on
the behavior of the waves and vortices in the inverse cascade of energy towards
the large scales. To this end, we use direct numerical simulations in a cubic
box at a Reynolds number Re=1000, with the ratio between the
Brunt-V\"ais\"al\"a frequency N and the inertial frequency f varying from 1/4
to 20, together with a purely stratified run. The Froude number, measuring the
strength of the stratification, varies within the range 0.02 < Fr < 0.32. We
find that the inverse cascade is dominated by the slow quasi-geostrophic modes.
Their energy spectra and fluxes exhibit characteristics of an inverse cascade,
even though their energy is not conserved. Surprisingly, the slow vortices
still dominate when the ratio N/f increases, also in the stratified case,
although less and less so. However, when N/f increases, the inverse cascade of
the slow modes becomes weaker and weaker, and it vanishes in the purely
stratified case. We discuss how the disappearance of the inverse cascade of
energy with increasing N/f can be interpreted in terms of the waves and
vortices, and identify three major effects that can explain this transition
based on inviscid invariants arguments
On the emergence of helicity in rotating stratified turbulence
We perform numerical simulations of decaying rotating stratified turbulence
and show, in the Boussinesq framework, that helicity (velocity-vorticity
correlation), as observed in super-cell storms and hurricanes, is spontaneously
created due to an interplay between buoyancy and rotation common to large-scale
atmospheric and oceanic flows. Helicity emerges from the joint action of eddies
and of inertia-gravity waves (with inertia and gravity with respective
associated frequencies and ), and it occurs when the waves are
sufficiently strong. For the amount of helicity produced is correctly
predicted by a quasi-linear balance equation. Outside this regime, and up to
the highest Reynolds number obtained in this study, namely ,
helicity production is found to be persistent for as large as , and for and respectively as large as and
.Comment: 10 pages, 5 figure
A hybrid MPI-OpenMP scheme for scalable parallel pseudospectral computations for fluid turbulence
A hybrid scheme that utilizes MPI for distributed memory parallelism and
OpenMP for shared memory parallelism is presented. The work is motivated by the
desire to achieve exceptionally high Reynolds numbers in pseudospectral
computations of fluid turbulence on emerging petascale, high core-count,
massively parallel processing systems. The hybrid implementation derives from
and augments a well-tested scalable MPI-parallelized pseudospectral code. The
hybrid paradigm leads to a new picture for the domain decomposition of the
pseudospectral grids, which is helpful in understanding, among other things,
the 3D transpose of the global data that is necessary for the parallel fast
Fourier transforms that are the central component of the numerical
discretizations. Details of the hybrid implementation are provided, and
performance tests illustrate the utility of the method. It is shown that the
hybrid scheme achieves near ideal scalability up to ~20000 compute cores with a
maximum mean efficiency of 83%. Data are presented that demonstrate how to
choose the optimal number of MPI processes and OpenMP threads in order to
optimize code performance on two different platforms.Comment: Submitted to Parallel Computin
GPU parallelization of a hybrid pseudospectral geophysical turbulence framework using CUDA
An existing hybrid MPI-OpenMP scheme is augmented with a CUDA-based fine grain parallelization approach for multidimensional distributed Fourier transforms, in a well-characterized pseudospectral fluid turbulence code. Basics of the hybrid scheme are reviewed, and heuristics provided to show a potential benefit of the CUDA implementation. The method draws heavily on the CUDA runtime library to handle memory management and on the cuFFT library for computing local FFTs. The manner in which the interfaces to these libraries are constructed, and ISO bindings utilized to facilitate platform portability, are discussed. CUDA streams are implemented to overlap data transfer with cuFFT computation. Testing with a baseline solver demonstrated significant aggregate speed-up over the hybrid MPI-OpenMP solver by offloading to GPUs on an NVLink-based test system. While the batch streamed approach provided little benefit with NVLink, we saw a performance gain of 30% when tuned for the optimal number of streams on a PCIe-based system. It was found that strong GPU scaling is nearly ideal, in all cases. Profiling of the CUDA kernels shows that the transform computation achieves 15% of the attainable peak FlOp-rate based on a roofline model for the system. In addition to speed-up measurements for the fiducial solver, we also considered several other solvers with different numbers of transform operations and found that aggregate speed-ups are nearly constant for all solvers.Fil: Rosenberg, Duane. State University of Colorado - Fort Collins; Estados UnidosFil: Mininni, Pablo Daniel. Consejo Nacional de Investigaciones CientÃficas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de FÃsica de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de FÃsica de Buenos Aires; ArgentinaFil: Reddy, Raghu. Environmental Modeling Center; Estados UnidosFil: Pouquet, Annick. State University of Colorado at Boulder; Estados Unidos. National Center for Atmospheric Research; Estados Unido
Anisotropy and non-universality in scaling laws of the large scale energy spectrum in rotating turbulence
Rapidly rotating turbulent flow is characterized by the emergence of columnar
structures that are representative of quasi-two dimensional behavior of the
flow. It is known that when energy is injected into the fluid at an
intermediate scale , it cascades towards smaller as well as larger scales.
In this paper we analyze the flow in the \textit{inverse cascade} range at a
small but fixed Rossby number, {}. Several
{numerical simulations with} helical and non-helical forcing functions are
considered in periodic boxes with unit aspect ratio. In order to resolve the
inverse cascade range with {reasonably} large Reynolds number, the analysis is
based on large eddy simulations which include the effect of helicity on eddy
viscosity and eddy noise. Thus, we model the small scales and resolve
explicitly the large scales. We show that the large-scale energy spectrum has
at least two solutions: one that is consistent with
Kolmogorov-Kraichnan-Batchelor-Leith phenomenology for the inverse cascade of
energy in two-dimensional (2D) turbulence with a {}
scaling, and the other that corresponds to a steeper {}
spectrum in which the three-dimensional (3D) modes release a substantial
fraction of their energy per unit time to 2D modes. {The spectrum that} emerges
{depends on} the anisotropy of the forcing function{,} the former solution
prevailing for forcings in which more energy is injected into 2D modes while
the latter prevails for isotropic forcing. {In the case of anisotropic forcing,
whence the energy} goes from the 2D to the 3D modes at low wavenumbers,
large-scale shear is created resulting in another time scale ,
associated with shear, {thereby producing} a spectrum for the
{total energy} with the 2D modes still following a {}
scaling
- …