Convergence of a Boundary Integral Method for Water Waves

We prove nonlinear stability and convergence of certain boundary integral methods for time-dependent water waves in a two-dimensional, inviscid, irrotational, incompressible fluid, with or without surface tension. The methods are convergent as long as the underlying solution remains fairly regular (and a sign condition holds in the case without surface tension). Thus, numerical instabilities are ruled out even in a fully nonlinear regime. The analysis is based on delicate energy estimates, following a framework previously developed in the continuous case [Beale, Hou, and Lowengrub, Comm. Pure Appl. Math., 46 (1993), pp. 1269–1301]. No analyticity assumption is made for the physical solution. Our study indicates that the numerical methods must satisfy certain compatibility conditions in order to be stable. Violation of these conditions will lead to numerical instabilities. A breaking wave is calculated as an illustration

Wavelet transforms and their applications to MHD and plasma turbulence: a review

Wavelet analysis and compression tools are reviewed and different applications to study MHD and plasma turbulence are presented. We introduce the continuous and the orthogonal wavelet transform and detail several statistical diagnostics based on the wavelet coefficients. We then show how to extract coherent structures out of fully developed turbulent flows using wavelet-based denoising. Finally some multiscale numerical simulation schemes using wavelets are described. Several examples for analyzing, compressing and computing one, two and three dimensional turbulent MHD or plasma flows are presented.Comment: Journal of Plasma Physics, 201

Evolution of vortex-surface fields in viscous Taylor-Green and Kida-Pelz flows

In order to investigate continuous vortex dynamics based on a Lagrangian-like formulation, we develop a theoretical framework and a numerical method for computation of the evolution of a vortex-surface field (VSF) in viscous incompressible flows with simple topology and geometry. Equations describing the continuous, timewise evolution of a VSF from an existing VSF at an initial time are first reviewed. Non-uniqueness in this formulation is resolved by the introduction of a pseudo-time and a corresponding pseudo-evolution in which the evolved field is ‘advected’ by frozen vorticity onto a VSF. A weighted essentially non-oscillatory (WENO) method is used to solve the pseudo-evolution equations in pseudo-time, providing a dissipative-like regularization. Vortex surfaces are then extracted as iso-surfaces of the VSFs at different real physical times. The method is applied to two viscous flows with Taylor–Green and Kida–Pelz initial conditions respectively. Results show the collapse of vortex surfaces, vortex reconnection, the formation and roll-up of vortex tubes, vorticity intensification between anti-parallel vortex tubes, and vortex stretching and twisting. A possible scenario for understanding the transition from a smooth laminar flow to turbulent flow in terms of topology of vortex surfaces is discussed

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

The analytic structure of 2D Euler flow at short times

Using a very high precision spectral calculation applied to the incompressible and inviscid flow with initial condition $\psi_0(x_1, x_2) = \cos x_1+\cos 2x_2$, we find that the width $\delta(t)$ of its analyticity strip follows a $\ln(1/t)$ law at short times over eight decades. The asymptotic equation governing the structure of spatial complex-space singularities at short times (Frisch, Matsumoto and Bec 2003, J.Stat.Phys. 113, 761--781) is solved by a high-precision expansion method. Strong numerical evidence is obtained that singularities have infinite vorticity and lie on a complex manifold which is constructed explicitly as an envelope of analyticity disks.Comment: 19 pages, 14 figures, published versio

Singularity formation in three-dimensional vortex sheets

We study singularity formation of three-dimensional (3-D) vortex sheets without surface tension using a new approach. First, we derive a leading order approximation to the boundary integral equation governing the 3-D vortex sheet. This leading order equation captures the most singular contributions of the integral equation. By introducing an appropriate change of variables, we show that the leading order vortex sheet equation degenerates to a two-dimensional vortex sheet equation in the direction of the tangential velocity jump. This change of variables is guided by a careful analysis based on properties of certain singular integral operators, and is crucial in identifying the leading order singular behavior. Our result confirms that the tangential velocity jump is the physical driving force of the vortex sheet singularities. We also show that the singularity type of the three-dimensional problem is similar to that of the two-dimensional problem. Moreover, we introduce a model equation for 3-D vortex sheets. This model equation captures the leading order singularity structure of the full 3-D vortex sheet equation, and it can be computed efficiently using fast Fourier transform. This enables us to perform well-resolved calculations to study the generic type of 3-D vortex sheet singularities. We will provide detailed numerical results to support the analytic prediction, and to reveal the generic form of the vortex sheet singularity

Geometry of enstrophy and dissipation, grid resolution effects and proximity issues in turbulence

We perform a multi-scale non-local geometrical analysis of the structures extracted from the enstrophy and kinetic energy dissipation-rate, instantaneous fields of a numerical database of incompressible homogeneous isotropic turbulence decaying in time obtained by DNS in a periodic box. Three different resolutions are considered: 256^3, 512^3 and 1024^3 grid points, with k_(max)η(overbar) approximately 1, 2 and 4, respectively, the same initial conditions and Re_λ ≈ 77. This allows a comparison of the geometry of the structures obtained for different resolutions. For the highest resolution, structures of enstrophy and dissipation evolve in a continuous distribution from blob-like and moderately stretched tube-like shapes at the large scales to highly stretched sheet-like structures at the small scales. The intermediate scales show a predominance of tube-like structures for both fields, much more pronounced for the enstrophy field. The dissipation field shows a tendency towards structures with lower curvedness than those of the enstrophy, for intermediate and small scales. The 256^3 grid resolution case (k_(max)η(overbar) ≈ 1) was unable to detect the predominance of highly stretched sheet-like structures at the smaller scales in both fields. The same non-local methodology for the study of the geometry of structures, but without the multi-scale decomposition, is applied to two scalar fields used by existing local criteria for the eduction of tube- and sheet-like structures in turbulence, Q and [A_ij]_+, respectively, obtained from invariants of the velocity-gradient tensor and alike in the 1024^3 case. This adds the non-local geometrical characterization and classification to those local criteria, assessing their validity in educing particular geometries. Finally, we introduce a new methodology for the study of proximity issues among structures of different fields, based on geometrical considerations and non-local analysis, by taking into account the spatial extent of the structures. We apply it to the four fields previously studied. Tube-like structures of Q are predominantly surrounded by sheet-like structures of [A_ij]_+, which appear at closer distances. For the enstrophy, tube-like structures at an intermediate scale are primarily surrounded by sheets of smaller scales of the enstrophy and structures of dissipation at the same and smaller scales. A secondary contribution results from tubes of enstrophy at smaller scales appearing at farther distances. Different configurations of composite structures are presented

A Computational Analysis of Wind Turbine and Wind Farm Aerodynamics with a Focus on Dual Rotor Wind Turbines

This dissertation serves to summarize my research into wind farm and wind turbine aerodynamics. Included in this thesis is a summary of the methods I use as well as the four research problems that I investigated. Motivation is provided for my research as well as an overview of the computational methods that I use. These methods include analytical methods such as blade element momentum (BEM) theory and the vortex lattice method as well as computational fluid dynamic methods like the Reynolds averaged Navier-Stokes (RANS) equations and large eddy simulation (LES). These methods are used to investigate wind turbine and wind farm aerodynamics. In particular, I use these methods to confront the various forms of loss that wind turbines and wind farms experience. They include the losses that individual turbines experience due to swirl, induction, and viscosity as well as the loss that wind farms experience due to turbine-wake interaction. Horizontal axis wind turbines (HAWTs) suffer from aerodynamic ineffciencies in the blade root region (near the hub) due to several non-aerodynamic constraints. Aerodynamic interactions between turbines in a wind farm also lead to signifcant loss of wind farm efficiency. A new dual-rotor wind turbine (DRWT) concept is proposed that aims at mitigating these two losses. A DRWT is designed that uses an existing turbine rotor for the main rotor, while the secondary rotor is designed using a high lift-to-drag ratio airfoil. Reynolds Averaged Navier-Stokes computational fluid dynamics simulations are used to optimize the design. Large eddy simulations confirm the increase energy capture potential of the DRWT. Wake comparisons however do not show enhanced entrainment of axial momentum. I extend the prescribed wake vortex lattice method (VLM) to perform aerodynamic analysis and optimization of dual-rotor wind turbines. The additional vortex system introduced by the secondary rotor of a DRWT is modeled while taking into account the singularities that occur when the trailing vortices from the secondary (upstream) rotor interact with the bound vortices of the main (downstream) rotor. Pseduo-steady assumption is invoked and averaging over multiple relative rotor positions is performed to account for the primary and secondary rotors operating at different rotational velocities. This implementation of the VLM is first validated against experiments and blade element momentum theory results for a conventional, single rotor turbine. The solver is then verified against RANS CFD results for two DRWTs. Parametric sweeps are performed using the proposed VLM algorithm to optimize a DRWT design. The problem with the algorithm at high loading conditions is highlighted and a solution is proposed that uses RANS CFD results to calibrate the VLM model. In addition to wake losses, aerodynamic interaction between turbines in wind farms leads to surface flow convergence . This phenomenon has been observed in field tests with surface flux stations. A hypothesis is proposed to explain this surface flow convergence phenomenon - incomplete pressure recovery behind a turbine leading to successive pressure drops in tightly-spaced turbine arrays leads to drop in overall pressure deep inside a wind plant; this low-pressure acts as an attractor leading to flow convergence. Numerical investigations of the phenomenon of surface flow convergence are carried out that support this hypothesis. An actuator disk model to represent wind turbines in an LES CFD solver is used to simulate hypothetical wind plants. The flow convergence phenomenon reflects as change in flow velocity direction and is more prominent near the ground than at turbine hub height. Numerical simulations of wind plant aerodynamics are conducted with various approximations to investigate and explain the flow convergence phenomenon