Methodology for tidal turbine representation in ocean circulation model

The present method proposes the use and adaptation of ocean circulation models as an assessment tool framework for tidal current turbine (TCT) array layout optimization. By adapting both momentum and turbulence transport equations of an existing model, the present TCT representation method is proposed to extend the actuator disc concept to 3-D large-scale ocean circulation models. Through the reproduction of experimental flume tests and grid dependency tests, this method has shown its numerical coherence as well as its ability to simulate accurately both momentum and turbulent turbine-induced perturbations in both near and far wakes in a relatively short period of computation time. Consequently the present TCT representation method is a very promising basis for the development of a TCT array layout optimization tool

Geometric study of Lagrangian and Eulerian structures in turbulent channel flow

We report the detailed multi-scale and multi-directional geometric study of both evolving Lagrangian and instantaneous Eulerian structures in turbulent channel flow at low and moderate Reynolds numbers. The Lagrangian structures (material surfaces) are obtained by tracking the Lagrangian scalar field, and Eulerian structures are extracted from the swirling strength field at a time instant. The multi-scale and multi-directional geometric analysis, based on the mirror-extended curvelet transform, is developed to quantify the geometry, including the averaged inclination and sweep angles, of both structures at up to eight scales ranging from the half-height δ of the channel to several viscous length scales δ_ν. Here, the inclination angle is on the plane of the streamwise and wall-normal directions, and the sweep angle is on the plane of streamwise and spanwise directions. The results show that coherent quasi-streamwise structures in the near-wall region are composed of inclined objects with averaged inclination angle 35°–45°, averaged sweep angle 30°–40° and characteristic scale 20δ_ν, and 'curved legs' with averaged inclination angle 20°–30°, averaged sweep angle 15°–30° and length scale 5δ_ν–10δ_ν. The temporal evolution of Lagrangian structures shows increasing inclination and sweep angles with time, which may correspond to the lifting process of near-wall quasi-streamwise vortices. The large-scale structures that appear to be composed of a number of individual small-scale objects are detected using cross-correlations between Eulerian structures with large and small scales. These packets are located at the near-wall region with the typical height 0.25δ and may extend over 10δ in the streamwise direction in moderate-Reynolds-number, long channel flows. In addition, the effects of the Reynolds number and comparisons between Lagrangian and Eulerian structures are discussed

Multi-scale geometric analysis of Lagrangian structures in isotropic turbulence

We report the multi-scale geometric analysis of Lagrangian structures in forced isotropic turbulence and also with a frozen turbulent field. A particle backward-tracking method, which is stable and topology preserving, was applied to obtain the Lagrangian scalar field φ governed by the pure advection equation in the Eulerian form ∂_tφ + u · ∇φ = 0. The temporal evolution of Lagrangian structures was first obtained by extracting iso-surfaces of φ with resolution 1024^3 at different times, from t = 0 to t = T_e, where T_e is the eddy turnover time. The surface area growth rate of the Lagrangian structure was quantified and the formation of stretched and rolled-up structures was observed in straining regions and stretched vortex tubes, respectively. The multi-scale geometric analysis of Bermejo-Moreno & Pullin (J. Fluid Mech., vol. 603, 2008, p. 101) has been applied to the evolution of φ to extract structures at different length scales and to characterize their non-local geometry in a space of reduced geometrical parameters. In this multi-scale sense, we observe, for the evolving turbulent velocity field, an evolutionary breakdown of initially large-scale Lagrangian structures that first distort and then either themselves are broken down or stretched laterally into sheets. Moreover, after a finite time, this progression appears to be insensible to the form of the initially smooth Lagrangian field. In comparison with the statistical geometry of instantaneous passive scalar and enstrophy fields in turbulence obtained by Bermejo-Moreno & Pullin (2008) and Bermejo-Moreno et al. (J. Fluid Mech., vol. 620, 2009, p. 121), Lagrangian structures tend to exhibit more prevalent sheet-like shapes at intermediate and small scales. For the frozen flow, the Lagrangian field appears to be attracted onto a stream-surface field and it develops less complex multi-scale geometry than found for the turbulent velocity field. In the latter case, there appears to be a tendency for the Lagrangian field to move towards a vortex-surface field of the evolving turbulent flow but this is mitigated by cumulative viscous effects

The interplay between helicity and rotation in turbulence: implications for scaling laws and small-scale dynamics

Invariance properties of physical systems govern their behavior: energy conservation in turbulence drives a wide distribution of energy among modes, observed in geophysical or astrophysical flows. In ideal hydrodynamics, the role of helicity conservation (correlation between velocity and its curl, measuring departures from mirror symmetry) remains unclear since it does not alter the energy spectrum. However, with solid body rotation, significant differences emerge between helical and non-helical flows. We first outline several results, like the energy and helicity spectral distribution and the breaking of strict universality for the individual spectra. Using massive numerical simulations, we then show that small-scale structures and their intermittency properties differ according to whether helicity is present or not, in particular with respect to the emergence of Beltrami-core vortices (BCV) that are laminar helical vertical updrafts. These results point to the discovery of a small parameter besides the Rossby number; this could relate the problem of rotating helical turbulence to that of critical phenomena, through renormalization group and weak turbulence theory. This parameter can be associated with the adimensionalized ratio of the energy to helicity flux to small scales, the three-dimensional energy cascade being weak and self-similar

On the non-local geometry of turbulence

A multi-scale methodology for the study of the non-local geometry of eddy structures in turbulence is developed. Starting from a given three-dimensional field, this consists of three main steps: extraction, characterization and classification of structures. The extraction step is done in two stages. First, a multi-scale decomposition based on the curvelet transform is applied to the full three-dimensional field, resulting in a finite set of component three-dimensional fields, one per scale. Second, by iso-contouring each component field at one or more iso-contour levels, a set of closed iso-surfaces is obtained that represents the structures at that scale. The characterization stage is based on the joint probability density function (p.d.f.), in terms of area coverage on each individual iso-surface, of two differential-geometry properties, the shape index and curvedness, plus the stretching parameter, a dimensionless global invariant of the surface. Taken together, this defines the geometrical signature of the iso-surface. The classification step is based on the construction of a finite set of parameters, obtained from algebraic functions of moments of the joint p.d.f. of each structure, that specify its location as a point in a multi-dimensional ‘feature space’. At each scale the set of points in feature space represents all structures at that scale, for the specified iso-contour value. This then allows the application, to the set, of clustering techniques that search for groups of structures with a common geometry. Results are presented of a first application of this technique to a passive scalar field obtained from 5123 direct numerical simulation of scalar mixing by forced, isotropic turbulence (Reλ = 265). These show transition, with decreasing scale, from blob-like structures in the larger scales to blob- and tube-like structures with small or moderate stretching in the inertial range of scales, and then toward tube and, predominantly, sheet-like structures with high level of stretching in the dissipation range of scales. Implications of these results for the dynamical behaviour of passive scalar stirring and mixing by turbulence are discussed

An efficient surrogate model for emulation and physics extraction of large eddy simulations

In the quest for advanced propulsion and power-generation systems, high-fidelity simulations are too computationally expensive to survey the desired design space, and a new design methodology is needed that combines engineering physics, computer simulations and statistical modeling. In this paper, we propose a new surrogate model that provides efficient prediction and uncertainty quantification of turbulent flows in swirl injectors with varying geometries, devices commonly used in many engineering applications. The novelty of the proposed method lies in the incorporation of known physical properties of the fluid flow as {simplifying assumptions} for the statistical model. In view of the massive simulation data at hand, which is on the order of hundreds of gigabytes, these assumptions allow for accurate flow predictions in around an hour of computation time. To contrast, existing flow emulators which forgo such simplications may require more computation time for training and prediction than is needed for conducting the simulation itself. Moreover, by accounting for coupling mechanisms between flow variables, the proposed model can jointly reduce prediction uncertainty and extract useful flow physics, which can then be used to guide further investigations.Comment: Submitted to JASA A&C

A Multiresolution Census Algorithm for Calculating Vortex Statistics in Turbulent Flows

The fundamental equations that model turbulent flow do not provide much insight into the size and shape of observed turbulent structures. We investigate the efficient and accurate representation of structures in two-dimensional turbulence by applying statistical models directly to the simulated vorticity field. Rather than extract the coherent portion of the image from the background variation, as in the classical signal-plus-noise model, we present a model for individual vortices using the non-decimated discrete wavelet transform. A template image, supplied by the user, provides the features to be extracted from the vorticity field. By transforming the vortex template into the wavelet domain, specific characteristics present in the template, such as size and symmetry, are broken down into components associated with spatial frequencies. Multivariate multiple linear regression is used to fit the vortex template to the vorticity field in the wavelet domain. Since all levels of the template decomposition may be used to model each level in the field decomposition, the resulting model need not be identical to the template. Application to a vortex census algorithm that records quantities of interest (such as size, peak amplitude, circulation, etc.) as the vorticity field evolves is given. The multiresolution census algorithm extracts coherent structures of all shapes and sizes in simulated vorticity fields and is able to reproduce known physical scaling laws when processing a set of voriticity fields that evolve over time