Vortex dislocations in wake-type flow induced by spanwise disturbances

Vortex dislocations in wake-type flow induced by three types of spanwise disturbances superimposed on an upstream velocity profile are investigated by direct numerical simulations. Three distinct modes of vortex dislocations and flow transitions have been found. A local spanwise exponential decay disturbance leads to the appearance of a twisted chainlike mode of vortex dislocation. A stepped spanwise disturbance causes a streamwise periodic spotlike mode of vortex dislocation. A spanwise sinusoidal wavy disturbance with a moderate waviness causes a strong unsteadiness of wake behavior. This unsteadiness starts with a systematic periodic mode of vortex dislocation in the spanwise direction followed by the spanwise vortex shedding suppressed completely with increased time and the near wake becoming a steady shear flow. Characteristics of these modes of vortex dislocation and complex vortex linkages over the dislocation, as well as the corresponding dynamic processes related to the appearance of dislocations, are described by examining the variations of vortex lines and vorticity distribution. The nature of the vortex dislocation is demonstrated by the substantial vorticity modification of the spanwise vortex from the original spanwise direction to streamwise and vertical directions, accompanied by the appearance of noticeable vortex branching and complex vortex linking, all of which are produced at the locations with the biggest phase difference or with a frequency discontinuity between shedding cells. The effect of vortex dislocation on flow transition, either to an unsteady irregular vortex flow or suppression of the Kaacutermaacuten vortex shedding making the wake flow steady state, is analyzed. Distinct similarities are found in the mechanism and main flow phenomena between the present numerical results obtained in wake-type flows and the experimental-numerical results of cylinder wakes reported in previous studies.

Two-dimensional Vortex-induced Vibrations of Submerged Floating Tunnel Tethers

A vortex-induced vibration (VIV) model is presented for predicting the nonlinear dynamic response of submerged floating tunnel (SFT) tethers which are subjected to wave, current and tunnel oscillatory displacements at their upper end in horizontal and vertical directions. A nonlinear fluid force formula is introduced in this model, and the effect of the nonlinearity of tether is investigated. First, the tunnel is stationary and the tether vibrates due to the vortices shedding. The calculated results show that the cross-flow amplitude of VIV decreases compared with the linear model. However the in-line amplitude of VIV increases. Next, the periodical oscillation of tunnel is considered. The oscillation caused by wave forces plays the roles of parametric exciter and forcing exciter to the VIV of tether. The time history of displacement of the tether mid-span is obtained by the proposed model. It is shown that the in-line amplitude increases obviously and the corresponding frequency is changed. The cross-flow amplitude exhibits a periodic behavior.

Topics in vortex motion

Six topics in incompressible, inviscid fluid flow involving vortex motion are presented. The stability of the unsteady flow field due to the vortex filament expanding under the influence of an axial compression is examined in the first chapter as a possible model of the vortex bursting observed in aircraft contrails. The filament with a stagnant core is found to be unstable to axisymmetric disturbances. For initial disturbances with the form of axisymmetric Kelvin waves, the filament with a uniformly rotating core is neutrally stable, but the compression causes the disturbance to undergo a rapid increase in amplitude. The time at which the increase occurs is, however, later than the observed bursting times, indicating the bursting phenomenon is not caused by this type of instability.

In the second and third chapters the stability of a steady vortex filament deformed by two-dimensional strain and shear flows, respectively, is examined. The steady deformations are in the plane of the vortex cross-section. Disturbances which deform the filament centerline into a wave which does not propagate along the filament are shown to be unstable and a method is described to calculate the wave number and corresponding growth rate of the amplified waves for a general distribution of vorticity in the vortex core.

In Chapter Four exact solutions are constructed for two-dimensional potential flow over a wing with a free ideal vortex standing over the wing. The loci of positions of the free vortex are found and the lift is calculated. It is found that the lift on the wing can be significantly increased by the free vortex.

The two-dimensional trajectories of an ideal vortex pair near an orifice are calculated in Chapter Five. Three geometries are examined, and the criteria for the vortices to travel away from the orifice are determined.

Finally, Chapter Six reproduces completely the paper, "Structure of a linear array of hollow vortices of finite cross-section," co-authored with G. R. Baker and P. G. Saffman. Free streamline theory is employed to construct an exact steady solution for a linear array of hollow, or stagnant cored vortices. If each vortex has area A and the separation is L, then there are two possible shapes if A^(1/2)/L is less than 0.38 and none if it is larger. The stability of the shapes to two-dimensional, periodic and symmetric disturbances is considered for hollow vortices. The more deformed of the two possible shapes is found to be unstable, while the less deformed shape is stable.

Studies in vortex motion

This thesis covers four different problems in the understanding of vortex sheets, and these are presented in four chapters.

In Chapter 1, free streamline theory is used to determine the steady solutions of an array of identical, hollow or stagnant core vortices in an inviscid, incompressible fluid. Assuming the array is symmetric to rotation through π radians about an axis through any vortex centre, there are two solutions or no solutions depending on whether A^(1/2)/L is less than or greater than 0.38 where A is the area of the vortex and L is the separation distance. Stability analysis shows that the more deformed shape is unstable to infinitesimal symmetric disturbances which leave the centres of the vortices undisplaced.

Chapter 2 is concerned with the roll-up of vortex sheets in homogeneous fluid. The flow over conventional and ring wings is used to test the method of Fink and Soh (1974). Despite modifications which improve the accuracy of the method, unphysical results occur. A possible explanation for this is that small scales are important and an alternate method based on "Cloud-in-Cell" techniques is introduced. The results show small scale growth and amalgamation into larger structures.

The motion of a buoyant pair of line vortices of opposite circulation is considered in Chapter 3. The density difference between the fluid carried by the vortices and the fluid outside is considered small, so that the Boussinesq approximation may be used. A macroscopic model is developed which shows the formation of a detrainment filament and this is included as a modification to the model. The results agree well with the numerical solution as developed by Hill (1975b) and show that after an initial slowdown, the vortices begin to accelerate downwards.

Chapter 4 reproduces completely a paper that has already been published (Baker, Barker, Bofah and Saffman (1974)) on the effect of "vortex wandering" on the measurement of velocity profiles of the trailing vortices behind a wing.

A dynamical systems analysis of vortex pinch-off

Vortex rings constitute the main structure in the wakes of a wide class of swimming and flying animals, as well as in cardiac flows and in the jets generated by some moss and fungi. However, there is a physical limit, determined by an energy maximization principle called the Kelvin-Benjamin principle, to the size that axisymmetric vortex rings can achieve. The existence of this limit is known to lead to the separation of a growing vortex ring from the shear layer feeding it, a process known as `vortex pinch-off', and characterized by the dimensionless vortex formation number. The goal of this thesis is to improve our understanding of vortex pinch-off as it relates to biological propulsion, and to provide future researchers with tools to assist in identifying and predicting pinch-off in biological flows.

To this end, we introduce a method for identifying pinch-off in starting jets using the Lagrangian coherent structures in the flow, and apply this criterion to an experimentally generated starting jet. Since most naturally occurring vortex rings are not circular, we extend the definition of the vortex formation number to include non-axisymmetric vortex rings, and find that the formation number for moderately non-axisymmetric vortices is similar to that of circular vortex rings. This suggests that naturally occurring vortex rings may be modeled as axisymmetric vortex rings. Therefore, we consider the perturbation response of the Norbury family of axisymmetric vortex rings. This family is chosen to model vortex rings of increasing thickness and circulation, and their response to prolate shape perturbations is simulated using contour dynamics. Finally, the response of more realistic models for vortex rings, constructed from experimental data using nested contours, to perturbations which resemble those encountered by forming vortices more closely, is simulated using contour dynamics. In both families of models, a change in response analogous to pinch-off is found as members of the family with progressively thicker cores are considered. We posit that this analogy may be exploited to understand and predict pinch-off in complex biological flows, where current methods are not applicable in practice, and criteria based on the properties of vortex rings alone are necessary.

Scanning tunneling spectroscopic studies on high-temperature superconductors and Dirac materials

This thesis details the investigations of the unconventional low-energy quasiparticle excitations in electron-type cuprate superconductors and electron-type ferrous superconductors as well as the electronic properties of Dirac fermions in graphene and three-dimensional strong topological insulators through experimental studies using spatially resolved scanning tunneling spectroscopy (STS) experiments.

Magnetic-field- and temperature-dependent evolution of the spatially resolved quasiparticle spectra in the electron-type cuprate La_0.1Sr_0.9CuO₂ (La-112) T_C = 43 K, are investigated experimentally. For temperature (T) less than the superconducting transition temperature (TC), and in zero field, the quasiparticle spectra of La-112 exhibits gapped behavior with two coherence peaks and no satellite features. For magnetic field measurements at T < TC, first ever observation of vortices in La-112 are reported. Moreover, pseudogap-like spectra are revealed inside the core of vortices, where superconductivity is suppressed. The intra-vortex pseudogap-like spectra are characterized by an energy gap of V_PG = 8.5 ± 0.6 meV, while the inter-vortex quasiparticle spectra shows larger peak-to-peak gap values characterized by Δ_pk-pk(H) >V_PG, and Δ_pk-pk (0)=12.2 ± 0.8 meV > Δ_pk-pk (H > 0). The quasiparticle spectra are found to be gapped at all locations up to the highest magnetic field examined (H = 6T) and reveal an apparent low-energy cutoff at the V_PG energy scale.

Magnetic-field- and temperature-dependent evolution of the spatially resolved quasiparticle spectra in the electron-type "122" iron-based Ba(Fe_1-xCo_x)₂A_s2 are investigated for multiple doping levels (x = 0.06, 0.08, 0.12 with T_C= 14 K, 24 K, and 20 K). For all doping levels and the T < T_C, two-gap superconductivity is observed. Both superconducting gaps decrease monotonically in size with increasing temperature and disappear for temperatures above the superconducting transition temperature, T_C. Magnetic resonant modes that follow the temperature dependence of the superconducting gaps have been identified in the tunneling quasiparticle spectra. Together with quasiparticle interference (QPI) analysis and magnetic field studies, this provides strong evidence for two-gap sign-changing s-wave superconductivity.

Additionally spatial scanning tunneling spectroscopic studies are performed on mechanically exfoliated graphene and chemical vapor deposition grown graphene. In all cases lattice strain exerts a strong influence on the electronic properties of the sample. In particular topological defects give rise to pseudomagnetic fields (B ~ 50 Tesla) and charging effects resulting in quantized conductance peaks associated with the integer and fractional Quantum Hall States.

Finally, spectroscopic studies on the 3D-STI, Bi₂Se₃ found evidence of impurity resonance in the surface state. The impurities are in the unitary limit and the spectral resonances are localized spatially to within ~ 0.2 nm of the impurity. The spectral weight of the impurity resonance diverges as the Fermi energy approaches the Dirac point and the rapid recovery of the surface state suggests robust topological protection against perturbations that preserve time reversal symmetry.

Identification of thermally-tagged coherent structures in the zero pressure gradient turbulent boundary layer

A zero pressure gradient boundary layer over a flat plate is subjected to step changes in thermal condition at the wall, causing the formation of internal, heated layers. The resulting temperature fluctuations and their corresponding density variations are associated with turbulent coherent structures. Aero-optical distortion occurs when light passes through the boundary layer, encountering the changing index of refraction resulting from the density variations. Instantaneous measurements of streamwise velocity, temperature and the optical deflection angle experienced by a laser traversing the boundary layer are made using hot and cold wires and a Malley probe, respectively. Correlations of the deflection angle with the temperature and velocity records suggest that the dominant contribution to the deflection angle comes from thermally-tagged structures in the outer boundary layer with a convective velocity of approximately 0.8U∞. An examination of instantaneous temperature and velocity and their temporal gradients conditionally averaged around significant optical deflections shows behavior consistent with the passage of a heated vortex. Strong deflections are associated with strong negative temperature gradients, and strong positive velocity gradients where the sign of the streamwise velocity fluctuation changes. The power density spectrum of the optical deflections reveals associated structure size to be on the order of the boundary layer thickness. A comparison to the temperature and velocity spectra suggests that the responsible structures are smaller vortices in the outer boundary layer as opposed to larger scale motions. Notable differences between the power density spectra of the optical deflections and the temperature remain unresolved due to the low frequency response of the cold wire.

General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

Dualidade na teoria de Landau-Ginzburg da supercondutividade

Neste trabalho abordamos a teoria de Ginzburg-Landau da supercondutividade (teoria GL). Apresentamos suas origens, características e resultados mais importantes. A idéia fundamental desta teoria e descrever a transição de fase que sofrem alguns metais de uma fase normal para uma fase supercondutora. Durante uma transição de fase em supercondutores do tipo II é característico o surgimento de linhas de fluxo magnético em determinadas regiões de tamanho finito chamadas comumente de vórtices. A dinâmica destas estruturas topológicas é de grande interesse na comunidade científica atual e impulsiona incontáveis núcleos de pesquisa na área da supercondutividade. Baseado nisto estudamos como essas estruturas topológicas influenciam em uma transição de fase em um modelo bidimensional conhecido como modelo XY. No modelo XY vemos que os principais responsáveis pela transição de fase são os vórtices (na verdade pares de vórtice-antivórtice). Villain, observando este fato, percebeu que poderia tornar explícita a contribuição desses defeitos topológicos na função de partição do modelo XY realizando uma transformação de dualidade. Este modelo serve como inspiração para a proposta deste trabalho. Apresentamos aqui um modelo baseado em considerações físicas sobre sistemas de matéria condensada e ao mesmo tempo utilizamos um formalismo desenvolvido recentemente na referência [29] que possibilita tornar explícita a contribuição dos defeitos topológicos na ação original proposta em nossa teoria. Após isso analisamos alguns limites clássicos e finalmente realizamos as flutuações quânticas visando obter a expressão completa da função correlação dos vórtices o que pode ser muito útil em teorias de vórtices interagentes (dinâmica de vórtices).

Behavior of spherical particles at low Reynolds numbers in a fluctuating translational flow

The behavior of spheres in non-steady translational flow has been studied experimentally for values of Reynolds number from 0.2 to 3000. The aim of the work was to improve our qualitative understanding of particle transport in turbulent gaseous media, a process of extreme importance in power plants and energy transfer mechanisms.

Particles, subjected to sinusoidal oscillations parallel to the direction of steady translation, were found to have changes in average drag coefficient depending upon their translational Reynolds number, the density ratio, and the dimensionless frequency and amplitude of the oscillations. When the Reynolds number based on sphere diameter was less than 200, the oscillation had negligible effect on the average particle drag.

For Reynolds numbers exceeding 300, the coefficient of the mean drag was increased significantly in a particular frequency range. For example, at a Reynolds number of 3000, a 25 per cent increase in drag coefficient can be produced with an amplitude of oscillation of only 2 per cent of the sphere diameter, providing the frequency is near the frequency at which vortices would be shed in a steady flow at the mean speed. Flow visualization shows that over a wide range of frequencies, the vortex shedding frequency locks in to the oscillation frequency. Maximum effect at the natural frequency and lock-in show that a non-linear interaction between wake vortex shedding and the oscillation is responsible for the increase in drag.

Distribuição do material particulado em suspensão na plataforma continental e talude superior entre Ubatuba (SP) e o cabo Búzios (RJ)

A pesquisa aqui apresentada amostrou a plataforma continental e talude superior entre Ubatuba (SP) e o cabo Búzios (RJ), em dois projetos de escala local e de meso-escala. A campanha entre Ubatuba e o cabo Frio foi realizada no âmbito do projeto Oceano-Rio: levantamentos oceanográficos integrados ao largo do Estado do Rio de Janeiro, realizado em colaboração com o Ministério de Ciência e Tecnologia (MCT) e com a cooperação da Marinha do Brasil (MB), a bordo do Navio Hidro-oceanográfico Cruzeiro do Sul. Dos dados coletados nesta campanha oceanográfica, foram analisados os dados de CTD, concentração de oxigênio dissolvido, fluorescência, transmissividade, atenuação, retro-espalhamento além de amostras de água para a determinação direta da concentração de material particulado em suspensão (MPS). No experimento de menor escala, entre o cabo Frio e o cabo Búzios, além destes dados, foi realizado um fundeio entre 02 e 30 de julho de 2011, que amostrou a intensidade do eco, as correntes e as ondas. As análises realizadas permitiram observar uma grande abrangência da massa dágua característica da Ressurgência sobre a plataforma continental entre Ubatuba (SP) e cabo Búzios (RJ), inclusive em regiões muito rasas. A ressurgência nem sempre ocorreu associada aos ventos NE, sugerindo que a ocorrência de vórtices tem grande influência sobre o Sistema de Ressurgência de Cabo Frio. O padrão de distribuição de MPS na região é muito influenciado pela ressuspensão dos sedimentos de fundo (aumentada durante as tempestades, quando a coluna dágua passa a apresentar maiores concentrações de MPS), produção fitoplanctônica e aporte das baías costeiras. Outros processos que parecem influenciar este padrão são os vórtices e as ondas internas no talude. A complexidade da interação destas diferentes fontes de MPS e processos associados à sua distribuição ficou evidenciada no processo de conversão dos sinais óticos do transmissômetro e nefelômetro em concentração de MPS. Apesar desta dificuldade associada à grande variabilidade das características do MPS, as análises dos testes de conversão do sinal ótico em concentração de MPS permitem concluir que as melhores conversões podem ser obtidas após a inspeção visual das concentrações de MPS medidos (filtração) e remoção das amostras que não seguem a tendência geral esperada de dispersão sinal ótico versus concentração de MPS.

Self-Similarity and helical symmetry in vortex generator flow simulations

According to experimental observations, the vortices generated by vortex generators have previously been observed to be self-similar for both the axial (u(z)) and azimuthal (u(circle minus)) velocity profiles. Further, the measured vortices have been observed to obey the criteria for helical symmetry. This is a powerful result, since it reduces the highly complex flow to merely four parameters. In the present work, corresponding computer simulations using Reynolds-Averaged Navier-Stokes equations have been carried out and compared to the experimental observations. The main objective of this study is to investigate how well the simulations can reproduce the physics of the flow and if the same analytical model can be applied. Using this model, parametric studies can be significantly reduced and, further, reliable simulations can substantially reduce the costs of the parametric studies themselves.

Produção de oscillons durante restaurações e quebras espontâneas de simetrias

Na natureza há vários fenômenos envolvendo transições de fase com quebra ou restauração de simetrias. Tipicamente, mudanças de fase, são associadas com uma quebra ou restauração de simetria, que acontecem quando um determinado parâmetro de controle é variado, como por exemplo temperatura, densidade, campos externos, ou de forma dinâmica. Essas mudanças que os sistemas sofrem podem levar a formação de defeitos topológicos, tais como paredes de domínios, vórtices ou monopolos magnéticos. Nesse trabalho estudamos particularmente mudanças de fase associadas com quebras ou restaurações dinâmicas de simetria que estão associadas com formação ou destruição de defeitos do tipo paredes de domínio em modelos de campos escalares com simetria discreta. Nesses processos dinâmicos com formação ou destruição de domínios, estudamos a possibilidade de formação de estruturas do tipo oscillons, que são soluções não homogêneas e instáveis de campo, mas que podem concentrar nelas uma quantidade apreciável de energia e terem uma vida (duração) suficientemente grande para serem de importância física. Estudamos a possibilidade de formação dessas soluções em modelos de dois campos escalares interagentes em que o sistema é preparado em diferentes situações, com a dinâmica resultante nesses sistemas estudada numa rede discreta.

Time-resolved volumetric measurement of fine-scale coherent structures in turbulence.

We present full volumetric (three-dimensional) time-resolved (+one-dimensional) measurements of the velocity field in a large water mixing tank, allowing us to assess spatial and temporal rotational energy (enstrophy) and turbulent energy dissipation intermittency. In agreement with previous studies, highly intermittent behavior is observed, with intense coherent flow structures clustering in the periphery of larger vortices. However, further to previous work the full volumetric measurements allow us to separate out the effects of advection from other effects, elucidating not only their topology but also the evolution of these intense events, through the local balance of stretching and diffusion. These findings contribute toward a better understanding of the intermittency phenomenon, which should pave the way for more accurate models of the small-scale motions based on an understanding of the underlying flow physics.