Modelling plane mixing layers using vortex points and sheets

We present here a critical assessment of two vortex approaches (both two-dimensional) to the modelling of turbulent mixing layers. In the first approach the flow is represented by point vortices, and in the second it is simulated as the evolution of a continuous vortex sheet composed of short linear elements or ''panels''. The comparison is based on fresh simulations using approximately the same number of elements in either model, paying due attention in both to the boundary conditions far downstream as well as those on the splitter plate from which the mixing layer issues. The comparisons show that, while both models satisfy the well-known invariants of vortex dynamics approximately to the same accuracy, the vortex panel model, although ultimately not convergent, leads to smoother roll-up and values of stresses and moments that are in closer agreement with the experiment, and has a higher computational efficiency for a given degree of convergence on moments. The point vortex model, while faster for a given number of elements, produces an unsatisfactory roll-up which (for the number of elements used) is rendered worse by the incorporation of the Van der Vooren correction for sheet curvature.

Vortex migration caused by travelling magnetic wave in a 2 in. diameter YBa2Cu3O7-δ thin film

We provide experimental evidence for a vortex migration phenomenon in YBa2Cu3O7-δ (YBCO) thin film caused by travelling magnetic wave. The experiment is carried out on a 2 in. diameter YBCO thin film with a circular-type magnetic flux pump. We found that the travelling wave helps the vortices migrate into the centre of the sample: after the zero-field cooling process, the increase of the flux density in the centre is four times larger than the amplitude of the travelling wave. The reason for this massive vortex migration is probably due to the magnetic stress variation caused by the travelling wave: the magnetic stress increases locally in the crest region while decreases locally in the trough region, which could help the vortices to move locally. A comparison shows that the magnetization by standing wave can be easily predicted by Bean's model while travelling wave causes vortex migration generally much larger than the prediction of Bean's model. It is possible that travelling magnetic wave can be an effective way to magnetize a type II superconductor in considering this unusual vortex dynamics. © 2013 AIP Publishing LLC.

Critical currents and melting temperature of a two-dimensional vortex lattice with periodic pinning

The critical current and melting temperature of a vortex system are analyzed. Calculations are made for a two-dimensional film at finite temperature with two kinds of periodic pinning: hexagonal and Kagomé. A transport current parallel and perpendicular to the main axis of the pinning arrays is applied and molecular dynamics simulations are used to calculate the vortex velocities to obtain the critical currents. The structure factor and displacements of vortices at zero transport current are used to obtain the melting temperature for both pinning arrays. The critical currents are higher for the hexagonal pinning lattice and anisotropic for both pinning arrays. This anisotropy is stronger with temperature for the hexagonal array. For the Kagomé pinning lattice, our analysis shows a multi stage phase melting; that is, as we increase the temperature, each different dynamic phase melts before reaching the melting temperature. Both the melting temperature and critical currents are larger for the hexagonal lattice, indicating the role for the interstitial vortices in decreasing the pinning strength. © 2012 Springer Science+Business Media New York.

Dynamics of the corotating vortices in dipolar Bose-Einstein condensates in the presence of dissipation

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Influence of thermal gradient in vortex states of mesoscopic superconductors

In general, the studies of finite size effects in mesoscopic superconductors have been carried out in such a way that the temperature parameter is constant in the entire system. However, we could have situations where a real sample is near a heater source, as an example. In such situations, gradients of temperature are present. On the other hand, mesoscopic superconductors are interesting systems due to the fact that they present confinement effects which influence all the vortex dynamics. Thus, in this work we studied the influence of thermal gradients on the vortex dynamics in mesoscopic superconductors. For this purposes, we used the time dependent Ginzburg-Landau equations. The thermal gradients produce an asymmetric distribution of the currents around the system which, in turn, yield interesting vortex configurations and difficult the formation of giant vortices.

The impacts of volcanic aerosol on stratospheric ozone and the Northern Hemisphere polar vortex: separating radiative-dynamical changes from direct effects due to enhanced aerosol heterogeneous chemistry

After major volcanic eruptions the enhanced aerosol causes ozone changes due to greater heterogeneous chemistry on the particle surfaces (HET-AER) and from dynamical effects related to the radiative heating of the lower stratosphere (RAD-DYN). We carry out a series of experiments with an atmosphere–ocean–chemistry–climate model to assess how these two processes change stratospheric ozone and Northern Hemispheric (NH) polar vortex dynamics. Ensemble simulations are performed under present day and preindustrial conditions, and with aerosol forcings representative of different eruption strength, to investigate changes in the response behaviour. We show that the halogen component of the HET-AER effect dominates under present-day conditions with a global reduction of ozone (−21 DU for the strongest eruption) particularly at high latitudes, whereas the HET-AER effect increases stratospheric ozone due to N2O5 hydrolysis in a preindustrial atmosphere (maximum anomalies +4 DU). The halogen-induced ozone changes in the present-day atmosphere offset part of the strengthening of the NH polar vortex during mid-winter (reduction of up to −16 m s-1 in January) and slightly amplify the dynamical changes in the polar stratosphere in late winter (+11 m s-1 in March). The RAD-DYN mechanism leads to positive column ozone anomalies which are reduced in a present-day atmosphere by amplified polar ozone depletion (maximum anomalies +12 and +18 DU for present day and preindustrial, respectively). For preindustrial conditions, the ozone response is consequently dominated by RAD-DYN processes, while under present-day conditions, HET-AER effects dominate. The dynamical response of the stratosphere is dominated by the RAD-DYN mechanism showing an intensification of the NH polar vortex in winter (up to +10 m s-1 in January). Ozone changes due to the RAD-DYN mechanism slightly reduce the response of the polar vortex after the eruption under present-day conditions.

Nonlinear evolution of harmonically forced perturbations on a wingtip vortex

Wingtip vortices are created by flying airplanes due to lift generation. The vortex interaction with the trailing aircraft has sparked researchers’ interest to develop an efficient technique to destroy these vortices. Different models have been used to describe the vortex dynamics and they all show that, under real flight conditions, the most unstable modes produce a very weak amplification. Another linear instability mechanism that can produce high energy gains in short times is due to the non-normality of the system. Recently, it has been shown that these non-normal perturbations also produce this energy growth when they are excited with harmonic forcing functions. In this study, we analyze numerically the nonlinear evolution of a spatially, pointwise and temporally forced perturbation, generated by a synthetic jet at a given radial distance from the vortex core. This type of perturbation is able to produce high energy gains in the perturbed base flow (10^3), and is also a suitable candidate for use in engineering applications. The flow field is solved for using fully nonlinear three-dimensional direct numerical simulation with a spectral multidomain penalty method model. Our novel results show that the nonlinear effects are able to produce locally small bursts of instability that reduce the intensity of the primary vortex.

Vortex evolution in the wake of accelerating low-aspect-ratio plates and three-dimensional bodies of revolution

To find examples of effecient locomotion and manoeuvrability, one need only turn to the elegant solutions natural flyers and swimmers have converged upon. This dissertation is specifically motivated by processes of evolutionary convergence, which have led to the propulsors and body shapes in nature that exhibit strong geometric collapse over diverse scales. These body features are abstracted in the studies presented herein using low-aspect-ratio at plates and a three-dimensional body of revolution (a sphere). The highly-separated vortical wakes that develop during accelerations are systematically characterized as a function of planform shape, aspect ratio, Reynolds number, and initial boundary conditions. To this end, force measurements and time-resolved (planar) particle image velocimetry have been used throughout to quantify the instantaneous forces and vortex evolution in the wake of the bluff bodies. During rectilinear motions, the wake development for the flat plates is primarily dependent on plate aspect ratio, with edge discontinuities and curvature playing only a secondary role. Furthermore, the axisymmetric case, i.e. the circular plate, shows strong sensitivity to Reynolds number, while this sensitivity quickly diminishes with increasing aspect ratio. For rotational motions, global insensitivity to plate aspect ratio has been observed. For the sphere, it has been shown that accelerations play an important role in the mitigation of flow separation. These results - expounded upon in this dissertation - have begun to shed light on the specific vortex dynamics that may be coopted by flying and swimming species of all shapes and sizes towards efficient locomotion.

Effect of granularity on the insulator-superconductor transition in ultrathin Bi films

We have studied the insulator-superconductor transition (IST) by tuning the thickness in quench-condensed Bi films. The resistive transitions of the superconducting films are smooth and can be considered to represent ''homogeneous'' films. The observation of an IST very close to the quantum resistance for pairs R-square(N) similar to h/4e(2) on several substrates supports this idea. The relevant length scales here are the localization length, and the coherence length. However, at the transition, the localization length is much higher than the superconducting coherence length, contrary to expectation for a "homogeneous" transition. This suggests the invalidity of a purely fermionic model for the transition. Furthermore, the current-voltage characteristics of the superconducting films are hysteretic, and show the films to be granular. The relevant energy scales here are the Josephson coupling energy and the charging energy. However, Josephson coupling energies (E-J) and the charging energies (E-c) at the IST, they are found to obey the relation E-j < E-c. This is again contrary to expectation, for the IST in a granular or inhomogeneous system. Hence, a purely bosonic picture of the transition is also inconsistent with our observations. We conclude that the IST observed in our experiments may be either an intermediate case between the fermioinc and bosonic mechanisms, or in a regime of charge and vortex dynamics for which a quantitative analysis has not yet been done.

Flexibility in flapping foil suppresses meandering of induced jet in absence of free stream

Thrust-generating flapping foils are known to produce jets inclined to the free stream at high Strouhal numbers St = fA/U-infinity, where f is the frequency and A is the amplitude of flapping and U-infinity is the free-stream velocity. Our experiments, in the limiting case of St —> infinity (zero free-stream speed), show that a purely oscillatory pitching motion of a chordwise flexible foil produces a coherent jet composed of a reverse Benard-Karman vortex street along the centreline, albeit over a specific range of effective flap stiffnesses. We obtain flexibility by attaching a thin flap to the trailing edge of a rigid NACA0015 foil; length of flap is 0.79 c where c is rigid foil chord length. It is the time-varying deflections of the flexible flap that suppress the meandering found in the jets produced by a pitching rigid foil for zero free-stream condition. Recent experiments (Marais et al., J. Fluid Mech., vol. 710, 2012, p. 659) have also shown that the flexibility increases the St at which non-deflected jets are obtained. Analysing the near-wake vortex dynamics from flow visualization and particle image velocimetry (PIV) measurements, we identify the mechanisms by which flexibility suppresses jet deflection and meandering. A convenient characterization of flap deformation, caused by fluid-flap interaction, is through a non-dimensional effective stiffness', EI* = 8 EI/(rho V-TEmax(2) s(f) c(f)(3)/2), representing the inverse of the flap deflection due to the fluid-dynamic loading; here, EI is the bending stiffness of flap, rho is fluid density, V-TEmax is the maximum velocity of rigid foil trailing edge, s(f) is span and c(f) is chord length of the flexible flap. By varying the amplitude and frequency of pitching, we obtain a variation in EI* over nearly two orders of magnitude and show that only moderate EI*. (0.1 less than or similar to EI * less than or similar to 1 generates a sustained, coherent, orderly jet. Relatively `stiff' flaps (EI* greater than or similar to 1), including the extreme case of no flap, produce meandering jets, whereas highly `flexible' flaps (EI* less than or similar to 0.1) produce spread-out jets. Obtained from the measured mean velocity fields, we present values of thrust coefficients for the cases for which orderly jets are observed.

A Note On The Numerical Simulations Of Flow Past A Wavy Square-Section Cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Karman vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode - "dual vortex pair mode" downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f(m) in addition to the vortex shedding frequency f(s) is identified. Frequency coupling between f(m) and f(s) occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

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).

An exact local conservation theorem for finite-amplitude disturbances to non-parallel shear flows, with remarks on Hamiltonian structure and on Arnol'd's stability theorems

Disturbances of arbitrary amplitude are superposed on a basic flow which is assumed to be steady and either (a) two-dimensional, homogeneous, and incompressible (rotating or non-rotating) or (b) stably stratified and quasi-geostrophic. Flow over shallow topography is allowed in either case. The basic flow, as well as the disturbance, is assumed to be subject neither to external forcing nor to dissipative processes like viscosity. An exact, local ‘wave-activity conservation theorem’ is derived in which the density A and flux F are second-order ‘wave properties’ or ‘disturbance properties’, meaning that they are O(a2) in magnitude as disturbance amplitude a [rightward arrow] 0, and that they are evaluable correct to O(a2) from linear theory, to O(a3) from second-order theory, and so on to higher orders in a. For a disturbance in the form of a single, slowly varying, non-stationary Rossby wavetrain, $\overline{F}/\overline{A}$ reduces approximately to the Rossby-wave group velocity, where (${}^{-}$) is an appropriate averaging operator. F and A have the formal appearance of Eulerian quantities, but generally involve a multivalued function the correct branch of which requires a certain amount of Lagrangian information for its determination. It is shown that, in a certain sense, the construction of conservable, quasi-Eulerian wave properties like A is unique and that the multivaluedness is inescapable in general. The connection with the concepts of pseudoenergy (quasi-energy), pseudomomentum (quasi-momentum), and ‘Eliassen-Palm wave activity’ is noted. The relationship of this and similar conservation theorems to dynamical fundamentals and to Arnol'd's nonlinear stability theorems is discussed in the light of recent advances in Hamiltonian dynamics. These show where such conservation theorems come from and how to construct them in other cases. An elementary proof of the Hamiltonian structure of two-dimensional Eulerian vortex dynamics is put on record, with explicit attention to the boundary conditions. The connection between Arnol'd's second stability theorem and the suppression of shear and self-tuning resonant instabilities by boundary constraints is discussed, and a finite-amplitude counterpart to Rayleigh's inflection-point theorem noted