Kinetic plasma turbulence cascade spans multiple scales ranging from macroscopic fluid flow to sub-electron scales. Mechanisms that dissipate large scale energy, terminate the inertial range cascade and convert kinetic energy into heat are hotly debated. Here we revisit these puzzles using fully kinetic simulation. By performing scale-dependent spatial filtering on the Vlasov equation, we extract information at prescribed scales and introduce several energy transfer functions. This approach allows highly inhomogeneous energy cascade to be quantified as it proceeds down to kinetic scales. The pressure work, − (P · ∇) · u, can trigger a channel of the energy conversion between fluid flow and random motions, which is a collision-free generalization of the viscous dissipation in collisional fluid. Both the energy transfer and the pressure work are strongly correlated with velocity gradients.

Many natural and engineering systems are simultaneously subjected to a driving force and a stabilizing force. The interplay between the two forces, especially for highly nonlinear systems such as fluid flow, often results in surprising features. Here we reveal such features in three different types of Rayleigh-Bénard (RB) convection, i.e. buoyancy-driven flow with the fluid density being affected by a scalar field. In the three cases different stabilizing forces are considered, namely (i) horizontal confinement, (ii) rotation around a vertical axis, and (iii) a second stabilizing scalar field. Despite the very different nature of the stabilizing forces and the corresponding equations of motion, at moderate strength we counterintuitively but consistently observe an enhancement in the flux, even though the flow motion is weaker than the original RB flow. The flux enhancement occurs in an intermediate regime in which the stabilizing force is strong enough to alter the flow structures in the bulk to a more organised morphology, yet not too strong to severely suppress the flow motions. Near the optimal transport enhancements all three systems exhibit a transition from a state in which the thermal boundary layer (BL) is nested inside the momentum BL to the one with the thermal BL being thicker than the momentum BL.

In this paper a numerical procedure to simulate low diffusivity scalar turbulence is presented. The method consists of using a grid for the advected scalar with a higher spatial resolutions than that of the momentum. The latter usually requires a less refined mesh and integrating both fields on a single grid tailored to the most demanding variable, produces an unnecessary computational overhead. A multiple resolution approach is used also in the time integration in order to maintain the stability of the scalars on the finer grid. The method is the more advantageous the less diffusive the scalar is with respect to momentum, therefore it is particularly well suited for large Prandtl or Schmidt number flows. However, even in the case of equal diffusivities the present procedure gives CPU time and memory occupation savings. The reason is that the absence of the pressure term in the scalar equation leads to much steeper gradients in the scalar field as compared to the velocity field.

The conservative cascade of kinetic energy is established using both Fourier analysis and a new exact physical-space flux relation in a simulated compressible turbulence. The subgrid scale (SGS) kinetic energy flux of the compressive mode is found to be significantly larger than that of the solenoidal mode in the inertial range, which is the main physical origin for the occurrence of Kolmogorov's -5/3 scaling of the energy spectrum in compressible turbulence. The perfect antiparallel alignment between the large-scale strain and the SGS stress leads to highly efficient kinetic energy transfer in shock regions, which is a distinctive feature of shock structures in comparison with vortex structures. The rescaled probability distribution functions of SGS kinetic energy flux collapse in the inertial range, indicating a statistical self-similarity of kinetic energy cascades.

The double diffusive convection between two parallel plates is numerically studied for a series of parameters. The flow is driven by the salinity difference and stabilised by the thermal field. Our simulations are directly compared to experiments by Hage and Tilgner (Phys. Fluids 22, 076603 (2010)) for several sets of parameters and reasonable agreement is found. This in particular holds for the salinity flux and its dependence on the salinity Rayleigh number. Salt fingers are present in all simulations and extend through the entire height. The thermal Rayleigh number seems to have minor influence on salinity flux but affects the Reynolds number and the morphology of the flow. Next to the numerical calculation, we apply the Grossmann-Lohse theory for Rayleigh-Bénard flow to the current problem without introducing any new coefficients. The theory successfully predicts the salinity flux both with respect to the scaling and even with respect to the absolute value for the numerical and experimental results.

In the traditional hybrid RANS/LES approaches for the simulation of wall-bounded fluid turbulence, such as detached-eddy simulation (DES), the whole flow domain is divided into an inner layer and an outer layer. Typically the Reynolds-averaged Navier–Stokes (RANS) equations are used for the inner layer, while large-eddy simulation (LES) is used for the outer layer. The transition from the inner-layer solution to the outer-layer solution is often problematic due to the lack of small-scale dynamics in the RANS region. In this paper, we propose to simulate the whole flow domain by large-eddy simulation while enforcing a Reynolds-stress constraint on the subgrid-scale (SGS) stress model in the inner layer. Both the algebraic eddy-viscosity model and the one-equation Spalart–Allmaras (SA) model have been used to constrain the Reynolds stress in the inner layer. In this way, we improve the LES methodology by allowing the mean flow of the inner layer to satisfy the RANS solution while small-scale dynamics is included. We validate the Reynolds-stress-constrained large-eddy simulation (RSC-LES) model by simulating three-dimensional turbulent channel flow and flow past a circular cylinder. Our model is able to predict mean velocity, turbulent stress and skin-friction coefficients more accurately in turbulent channel flow and to estimate the pressure coefficient after separation more precisely in flow past a circular cylinder compared with the pure dynamic Smagorinsky model (DSM) and DES using the same grid resolution. Furthermore, the computational cost of the RSC-LES is almost the same as that of DES.

The AFiD code, an open source solver for the incompressible Navier-Stokes equations (http://www.afid.eu), has been ported to GPU clusters to tackle large-scale wall-bounded turbulent flow simulations. The GPU porting has been carried out in CUDA Fortran with the extensive use of kernel loop directives (CUF kernels) in order to have a source code as close as possible to the original CPU version; just a few routines have been manually rewritten. A new transpose scheme, which is not limited to the GPU version only and can be generally applied to any CFD code that uses pencil distributed parallelization, has been devised to improve the scaling of the Poisson solver, the main bottleneck of incompressible solvers. The GPU version can reduce the wall clock time by an order of magnitude compared to the CPU version for large meshes. Due to the increased performance and efficient use of memory, the GPU version of AFiD can perform simulations in parameter ranges that are unprecedented in thermally-driven wall-bounded turbulence. To verify the accuracy of the code, turbulent Rayleigh-Bénard convection and plane Couette flow are simulated and the results are in good agreement with the experimental and computational data that published in previous literatures. PROGRAM SUMMARYProgram Title: AFiD-GPU Licensing provisions(please choose one): GPLv3 Programming language: Fortan 90, CUDA Fortan, MPI External routines: PGI, CUDA Toolkit, FFTW3, HDF5 Nature of problem(approx. 50-250 words): Solving the three-dimensional Navier-Stokes equations coupled with a scalar field in a cubic box bounded between two walls and other four periodic boundaries. Solution method(approx. 50-250 words): Second order finite difference method for spatial discretization, third order Runge-Kutta scheme and Crank-Nicolson method for time advancement, two dimensional pencil distributed MPI parallelization, GPU accelerated routines. Additional comments including Restrictions and Unusual features (approx. 50-250 words): The open-source code is supported and updated on http://www.afid.eu.

Direct numerical simulations are conducted for double diffusive convection (DDC) bounded by two parallel plates. The Prandtl numbers, i.e. the ratios between the viscosity and the molecular diffusivities of scalars, are similar to the values of seawater. The DDC flow is driven by an unstable salinity difference (here across the two plates) and stabilized at the same time by a temperature difference. For these conditions the flow can be in the finger regime. We develop scaling laws for three key response parameters of the system: The non-dimensional salinity flux Nu S mainly depends on the salinity Rayleigh number Ra S , which measures the strength of the salinity difference, and exhibits a very weak dependence on the density ratio Λ, which is the ratio of the buoyancy forces induced by two scalar differences. The non-dimensional flow velocity Re and the non-dimensional heat flux Nu T are dependent on both Ra S and Λ. However, the rescaled Reynolds number ReΛ α eff u and the rescaled convective heat flux (Nu T − 1)Λ α eff T depend only on Ra S . The two exponents are dependent on the fluid properties and are determined from the numerical results as α eff u = 0.25 ± 0.02 and α eff T = 0.75 ± 0.03. Moreover, the behaviors of Nu S and ReΛ α eff u agree with the predictions of the Grossmann-Lohse theory which was originally developed for the Rayleigh-Bénard flow. The non-dimensional salt-finger width and the thickness of the velocity boundary layers, after being rescaled by Λ α eff u /2 , collapse and obey a similar power-law scaling relation with Ra S . When Ra S is large enough, salt fingers do not extend from one plate to the other and horizontal zonal flows emerge in the bulk region. We then show that the current scaling strategy can be successfully applied to the experimental results of a heat-copper-ion system (Hage and Tilgner, Phys. Fluids, 22, 076603, 2010). The fluid has different properties and the exponent α eff u takes a different value 0.54 ± 0.10.