Phase and vortex dynamics in Josephson junction arrays with percolative disorder

The AC conductance G of site-diluted Josephson junction arrays close to the percolation threshold was studied over a broad frequency range. As a function of frequency, G exhibits a crossover from a low-frequency Euclidean to a high-frequency fractal regime. At low-frequency the response is dominated by vortices and its temperature dependence is well-described by dynamical extensions of the Kosterlitz–Thouless theory for two-dimensional systems.     

PET/PTT双组分弹性长丝的结晶取向结构和卷曲性能

为研制军官礼服用PET/PTT双组分弹性长丝,在纺丝加工工艺研究的基础上,通过声速法、WAXD、DSC、Instron5566对典型工艺下的弹性长丝进行了结晶和取向结构及卷曲性能的测试分析.在可纺的前提下,PET/PTT两组分复合纺丝中,PET组分优先结晶,具有高于其单组分纤维的拉伸诱导取向和结晶;而PTT组分只有形变,其结晶度和晶区取向均低于其对应的单组分纤维.在实验条件范围内,两组分粘度差异越大,纤维的卷曲伸长率和收缩率越大、声速取向因子增加、各单组分结晶度增加;两组分质量比为50/50时,纤维有最大的卷曲伸长率和收缩率,且各单组分结晶度随该两组分含量差异的增加而减少,而声速取向变化相反;随牵伸比的增加,纤维的整体取向、各组分结晶度均有所增加,卷曲伸长和收缩率也增加.牵伸温度和定型温度对双组分纤维的结构和卷曲性能影响较小.     

IMPLEMENTATION OF VORTEX FILAMENT METHODS ON PARALLEL MACHINES WITH DISTRIBUTED ADAPTIVE DATA STRUCTURE

This paper addressed the implementation of vortex filament methods on parallel machines with distributed memory to simulate a three-dimensionally evolving jet. Vortical structure developments due to Kelvin–Helmholtz instability of the axially perturbed jet are also examined. The implementation is conducted in a single-programme multiple-data (SPMD) environment and the parallelism is focused on issues of data distribution, efficient support of parallel I/O and overlapping of communications with computations. In addition, since the number of segment markers in a filament is dynamically growing according to the requirement of numerical accuracy, a novel packet-oriented data structure is proposed not only to partition filament segment markers among distributed processors but also to support dynamical load balancing at run time. This work is the first to apply packet-oriented structures to implement a parallel vortex filament method. Experimental results indicate performance improvement from 1·5 to 2·6 times over static schemes on nCUBE2, DEC Alpha and IBM SP2 by incorporating the proposed scheme with packet-oriented structures. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids 24: 939–951, 1997.     

Numerical simulation of unsteady interaction of centrifugal impeller with its diffuser using Lagrangian discrete vortex method

In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.     

Exotic vortex structures of the dipolar Bose–Einstein condensates trapped in harmonic-like and toroidal potential

Based on the tunable intensity and waist of Gaussian laser, harmonic-like and toroidal potentials can be achieved and the ground-state properties of the dipolar Bose–Einstein condensate (BEC) trapped in such potentials are investigated. It is found that, in the harmonic-like potential, the singly and doubly quantized vortices can exist in the scale condensate and translate respectively into vortex pairs and triangular vortex lattice with increasing dipole–dipole interaction (DDI). Especially, the sandwich-like structure can be observed in the ground-state density profiles by tuning the direction and strength of DDI for some rotating frequency. In the toroidal potential, the competition between the inter-component interaction and DDI can induce the transition between immiscible and miscible states, and results in the structures of a doubly quantized vortex surrounded by a vortex ring. It is worth emphasizing that, with the increasing of DDI, the doubly quantized vortex in the harmonic-like potential becomes two singly quantized vortices, while in the toroidal potential it is no happen due to the presence of Gaussian barrier.     

A three-dimensional numerical study on instability of sinusoidal flame induced by multiple shock waves

The instabilities of a three-dimensional sinusoidally premixed flame induced by an incident shock wave with Mach = 1.7 and its reshock waves were studied by using the Navier–Stokes(NS) equations with a single-step chemical reaction and a high resolution, 9th-order weighted essentially non-oscillatory scheme. The computational results were validated by the grid independence test and the experimental results in the literature. The computational results show that after the passage of incident shock wave the flame interface develops in symmetric structure accompanied by large-scale transverse vortex structures. After the interactions by successive reshock waves, the flame interface is gradually destabilized and broken up, and the large-scale vortex structures are gradually transformed into small-scale vortex structures. The small-scale vortices tend to be isotropic later.The results also reveal that the evolution of the flame interface is affected by both mixing process and chemical reaction. In order to identify the relationship between the mixing and the chemical reaction, a dimensionless parameter, η, that is defined as the ratio of mixing time scale to chemical reaction time scale, is introduced. It is found that at each interaction stage the effect of chemical reaction is enhanced with time.The enhanced effect of chemical reaction at the interaction stage by incident shock wave is greater than that at the interaction stages by reshock waves. The result suggests that the parameter η can reasonably character the features of flame interface development induced by the multiple shock waves.     

Vortex shedding around a near-wall circular cylinder induced by a solitary wave

This study developed a two-dimensional generalized vortex method to analyze the shedding of vortices and the hydrodynamic forces resulting from a solitary wave passing over a submerged circular cylinder placed near a flat seabed. Numerical results for validation are compared with other numerical and experimental results, and satisfactory agreement is found. A series of simulations were performed to study the effects of gap-to-diameter ratio and incident wave height on vorticity pattern as well as the forces exerted on the cylinder. The range of the heights of incident waves is from 0.3h to 0.7h, where h is the still water depth. The range of the gap-to-diameter ratios is from 0.1 to 0.8. The results indicate that the flow pattern and the pressure distribution change significantly because of the close proximity of the seabed where the vorticity flux on the seabed-side surface of the cylinder is suppressed. Placing the cylinder nearer the seabed increases the drag and the positive lift on the cylinder. When the gap-to-diameter ratio increases, the pattern of vortices changes because of the interaction between the main recirculation zone and the shear layers separated from the gap. The maxima of drag, lift and total force increase linearly with the height of the incident wave.     

Numerical investigations of the vortex interactions for a flow over a pitching foil at different stages

The fluid–structure interaction is investigated numerically for a two-dimensional flow (Re=2.5·10⁶) over a sinusoid-pitching foil by the SST (Shear Stress Transport) k–ω model. Although discrepancies in the downstroke phase, which are also documented in other numerical studies, are observed by comparing with experimental results, our current numerical results are sufficient to predict the mean features and qualitative tendencies of the dynamic stall phenomenon. These discrepancies are evaluated carefully from the numerical and experimental viewpoints.In this study, we have utilized Λ, which is the normalized second invariant of the velocity gradient tensor, to present the evolution of the Leading Edge Vortex (LEV) and Trailing Edge Vortex (TEV). The convective, pressure, and diffusion terms during the dynamic stall process are discussed based on the transport equation of Λ. It is found that the pressure term dominates the rate of the change of the rotation strength inside the LEV. This trend can hardly be observed directly by using the vorticity transport equation due to the zero baroclinic term for the incompressible flow.The mechanisms to delay the stall are categorized based on the formation of the LEV. At the first stage before the formation of the LEV in the upper surface, the pitching foil provides extra momentum into the fluid flows to resist the flow separation, and hence the stall is delayed. At the second stage, a low-pressure area travels with the evolution of the LEV such that the lift still can be maintained. Three short periods at the second stage corresponds to different flow patterns during the dynamic stall, and these short periods can be distinguished according to the trend of the pressure variation inside the LEV. The lift stall occurs when a reverse flow from the lower surface is triggered during the shedding of the LEV. For a reduced frequency k_f=0.15, the formation of the TEV happens right after the lift stall, and the lift can drop dramatically. With a faster reduced frequency k_f=0.25, the shedding of the LEV is postponed into the downstroke, and the interaction between the LEV and TEV becomes weaker correspondingly. Thus, the lift drops more gently after the stall. In order to acquire more reliable numerical results within the downstroke phase, the Large Eddy Simulation (LES), which is capable of better predictions for the laminar-to-turbulent transition and flow reattachment process, will be considered as the future work.     

Direct numerical simulation of free and forced square jets

Direct numerical simulation (DNS) of incompressible, spatially developing square jets in the Reynolds number range of 500–2000 is reported. The three-dimensional unsteady Navier–Stokes equations are solved using high order spatial and temporal discretization. The objective of the present work is to understand the evolution of free and forced square jets by examining the formation of large-scale structures. Coherent structures and related interactions of free jets suggest control strategies that can be used to achieve enhanced spreading and mixing of the jet with the surrounding fluid. The critical Reynolds number for the onset on unsteadiness in an unperturbed free square jet is found to be 875–900 while it reduces to the range 500–525 in the presence of small-scale perturbations. Disturbances applied at the flow inlet cause saturation of KH-instability and early transition to turbulence. Forced jet calculations have been carried out using varicose perturbation with amplitude of 15%, while frequency is independently varied. Simulations show that the initial development of the square jet is influenced by the four corners leading to the appearance hairpin structures along with the formation of vortex rings. Farther downstream, adjacent vortices strongly interact leading to their rapid breakup. Excitation frequencies in the range 0.4–0.6 cause axis-switching of the jet cross-section. Results show that square jets achieve greater spreading but are less controllable in comparison to the circular ones.