高超临界雷诺数区间内二维圆柱绕流的实测研究

实测强风工况下高度167 m的徐州彭城电厂冷却塔的表面风荷载,并归纳历史上其他研究人员给出的实测结果,以丰富高超临界雷诺数(Re)区间二维圆柱绕流的试验成果.在低湍流度均匀流场和高湍流度大气边界层流场中分别开展4种风速8类粗糙度条件下的冷却塔刚性模型测压风洞试验,通过对比低雷诺数(Re=2.1×10~5—4.19×10~5)条件下的风洞试验结果和高雷诺数(Re=5.4×10~7—1×10~8)条件下的现场实测结果研究各种静动态绕流特征随雷诺数的变化规律,重点考察雷诺数无关现象的产生条件.研究结果表明,对于物表相对粗糙度在0.01以上的圆柱绕流,雷诺数不相关现象存在于很宽的雷诺数范围(2×10~5Re1×10~8)内;增大来流湍流度亦能引起的雷诺数无关现象,但此时该现象可能仅存在于一个较窄的低雷诺数范围内.     

Investigation of turbulent flow past a flagged short finite circular cylinder

This study was conducted to investigate the flow structures of turbulent flow passing over a short finite cylinder in which a rigid flag is attached to the rear of the cylinder, in wake region. The length-to-diameter ratio of the cylinder was chosen to be L/D = 2, whereas the rigid flag had a width-to-diameter ratio of W/D = 1.5. Wall-adapted large-eddy simulation (LES-WALE) was used to resolve unsteady turbulent flow structures. The far field Reynolds number based on cylinder diameter was chosen to be 20,000. The results were compared with the regular case wherein no flag was attached to the cylinder. Results revealed that the flow pattern behind the cylinder with flag was totally different in comparison with the regular case one. However, top free end of the cylinder was not influenced by the flag in contrast with the wake region. At far downstream from the cylinder, most of the flow structures in both cases appeared the same. The horseshoe vortices in both cases appeared to be an unsteady phenomenon, with slightly different patterns. Moreover, in the case of flag attachment, the pressure coefficient was smaller than that of with no flag. Finally, it was shown that the main and secondary Strouhal numbers locations were different in both cases.     

Visualisation and analysis of large-scale vortex structures in three-dimensional turbulent lid-driven cavity flow

In this paper, large eddy simulation (LES) of a three-dimensional turbulent lid-driven cavity (LDC) flow at Re = 10,000 has been performed using the multiple relaxation time lattice Boltzmann method. A Smagorinsky eddy viscosity model was used to represent the sub-grid scale stresses with appropriate wall damping. The prediction for the flow field was first validated by comparing the velocity profiles with previous experimental and LES studies, and then subsequently used to investigate the large-scale three-dimensional vortical structures in the LDC flow. The instantaneous three-dimensional coherent structures inside the cavity were visualised using the second invariant (Q), Δ criterion, λ₂ criterion, swirling strength (λ_ci) and streamwise vorticity. The vortex structures obtained using the different criteria in general agree well with each other. However, a cleaner visualisation of the large vortex structures was achieved with the λ_ci criterion and also when the visualisation is based on the vortex identification criteria expressed in terms of the swirling strength parameters. A major objective of the study was to perform a three-dimensional proper orthogonal decomposition (POD) on the fluctuating velocity fields. The higher energy POD modes efficiently extracted the large-scale vortical structures within the flow which were then visualised with the swirling strength criterion. Reconstruction of the instantaneous fluctuating velocity field using a finite number of POD modes indicated that the large-scale vortex structures did effectively approximate the large-scale motion. However, such a reduced order reconstruction of the flow based on the large-scale vortical structures was clearly not as effective in predicting the small-scale details of the fluctuating velocity field which relate to the turbulent transport.     

Simulation of the inherent turbulence and wake interaction inside an infinitely long row of wind turbines

The turbulence in the interior of a wind farm is simulated using large eddy simulation and the actuator line technique implemented in the Navier–Stokes equations. The simulations are carried out for an infinitely long row of turbines simulated by applying cyclic boundary conditions at the inlet and outlet. The simulations investigate the turbulence inherent to the wind turbines as no ambient turbulence or shear is added to this idealised case. The simulated data give insight into the performance of the wind turbines operating in the wake of others as well as details on key turbulent quantities. One of the key features of wakes behind wind turbines is the dynamic wake meandering, which is shown to be related to the wind turbine spacing and the vortex shedding from the turbine as a bluff body. The flow is analysed and reconstructed by applying proper orthogonal decomposition.     

Control of a coupled map lattice model for vortex shedding in the wake of a cylinder

The flow behind a vibrating flexible cable at low Reynolds numbers can exhibit complex wake structures such as lace-like patterns, vortex dislocations and frequency cells. These structures have been observed in experiments and numerical simulations, and are predicted by a previously developed low-order coupled map lattice (CML). The discrete (in time and space) CML models consist of a series of diffusively coupled circle map oscillators along the cable span. Motivated by a desire to modify the complex wake patterns behind flexible vibrating cables we have studied the addition of control terms into the highly efficient CML models and explored the resulting dynamics. Proportional, adaptive proportional and discontinuous non-linear (DNL) control methods were used to derive the control laws. The first method employed occasional proportional feedback. The adaptive method used spatio-temporal feedback control. The DNL method used a discontinuous feedback linearization procedure, and the controller was designed for the resulting linearized system using eigenvalue assignment. These techniques were applied to a modeled vortex dislocation structure in the wake of a vibrating cable in uniform freestream flow. Parallel shedding patterns were achieved for a range of forcing frequency-forcing amplitude combinations studied to validate the control theory. The adaptive proportional and DNL methods were found to be more effective than the proportional control method due to the incorporation of a spatially varying feedback gain across the cylinder span. The DNL method was found to be the most efficient controller of the low-order CML model. The required control level across the cable span was correlated to the 1/1 lock-on behavior of the temporal circle map.     

Influences of initial velocity, diameter and Reynolds number on a circular turbulent air/air jet

This paper assesses the suitability of the inflow Reynolds number defined by Re_o ≡ U_oD/ν (here U_o and D are respectively the initial jet velocity and diameter while ν is kinematic viscosity) for a round air/air jet. Specifically an experimental investigation is performed for the influences of U_o, D and Re_o on the mean-velocity decay and spread coefficients (K_u, K_r) in the far field of a circular air jet into air from a smoothly contracting nozzle. Present measurements agree well with those previously obtained under similar inflow conditions. The relations K_u ∝ U_o and K_r ∝ 1/U_o for U_o < 5 m/s appear to work, while each coefficient approaches asymptotically to a constant for U_o > 6 m/s, regardless of the magnitudes of Re_o and D. It is revealed that Re_o may not be an appropriate dimensionless parameter to characterize the entire flow of a free air/air jet. This paper is the first paper that has challenged the suitability of Re_o for turbulent free jets.     

Scattering of an arbitrarily oriented dipole field by an infinite and a finite length PEMC circular cylinder

Electromagnetic scattering from an infinite and a finite length PEMC circular cylinder, illuminated by an arbitrarily oriented dipole, is investigated theoretically. An electric dipole as a source of excitation is considered first, and then a magnetic dipole as a source of excitation is treated. In contrast to the case of an axially directed dipole, it is shown that no additional terms are needed to incorporate the cross-polarized component of the field for the case of radial and circumferential dipoles. Numerical verifications are presented to verify the validity of derived results and numerical code by comparing results with the published literature.      

Vortex identification in the wake of a model wind turbine array

Vortex identification techniques are used to analyse the wake of a 4 × 3 array of model wind turbines. The Q-criterion, Δ-criterion, and λ₂-criterion are applied to particle image velocimetry data gathered fore and aft of the last row centerline turbine. The Q-criterion and λ₂-criterion provide a clear indication of regions where vortical activity exists, while the Δ-criterion is not successful. Galilean decomposition, Reynolds decomposition, vorticity, and swirling strength are used to further understand the location and behaviour of the vortices. The techniques identify and display the high-magnitude vortices in high-shear zones resulting from the blade tips. Using Galilean and Reynolds decomposition, swirling motions are shown encapsuling vortex regions in agreement with the identification criteria. The Galilean decompositions selected are 20% and 50% of a convective velocity of 7 m/s. As the vortices convect downstream, the strength of the vortices decreases in magnitude, particularly in the far wake of the array, to approximately 25% of those present in the near wake. A high level of vortex activity is visualised as a result of the top tip of the wind turbine blade -- the location where the highest vertical entrainment is present. Analysing the full frame set, the Q-criterion, λ₂-criterion, and swirling strength prove comparable, while the Δ-criterion under-performs in regions of high turbulence activity, namely in the back of the turbine. Entraining flow into the turbine canopy interacting with the turbine generates high-magnitude vortices concentrated at the blade tips. The count of vortices decreases when moving from the top tip down to the wall, as well as their strength for each Galilean technique when a non-zero threshold is applied. Vortex sizes in the near wake are found comparable to turbine blade, hub, and mast dimensions. In the far wake, the resulting size of the vortices is approximately 30% of those in the near wake. These vortices increase in velocity as they convect downstream, following the mean velocity behaviour. The lowest magnitude vortices reside at the hub height in the near-wake region, where they convect at nearly half the speed of those at the blade tips.     

The onset of chaos in the wake of an oscillating cylinder: Experiment and the dynamics of the circle map

This paper deals with a comparison between experimental observations in a low-Reynolds-number wake behind an oscillating cylinder and the universal properties of a sine circle map. When the limit cycle due to the natural vortex shedding in the wake is modulated at a second frequency by oscillating the cylinder transversely, one obtains in phase space a flow on a two torus. The nonlinear interaction between the two oscillators results in Arnol’d tongues due to phase locking, the devil’s staircase along the critical line, and a transition from order to chaosvia the quasiperiodic route. The sine circle map describes these features adequately. A comparison between the experiment and the theory is made in terms of multifractal formalism and trajectory scaling function.     

圆柱形人造原子的能级结构及电子几率密度分布

通过求解定态薛定格方程得到了园柱形人造原子的能级结构,进而求出了电子几率密度的分布。计算结果表明,电子被束缚在园柱内的条件与位阱深度Ｖ０及园柱的园柱的结构参数有关。当园柱的高ｈ＝４．００ｎｍ,半径α＝２．７５ｎｍ时,要求Ｖ０≥４．７４ｍｅＶ。在ｚ－轴方向上的电子几率密度分布是以ｈ／２为对称中心周期震荡,并且柱体内外的电子几率的大小是可比较的;而在ｘ－ｙ平面上的电子几率密度分布类似于氢原子,并且柱外     

Vortex shedding and vorticity fluxes in the wake of cylinders within a random array

This work investigates the impact of the background turbulence generated by randomly placed cylinders on the vortex shedding regime and the mechanisms associated to vorticity fluxes. The goals are achieved by exploring velocity databases acquired with a two-dimensional particle image velocimetry system in two types of turbulent flow experiments: flow around a single infinite cylinder and flow within random array of infinite cylinders. Formation lengths, power spectral density functions and vortex distributions are employed to discuss the vortex shedding regime. The effects of background turbulence and vorticity cancellation, due to opposite sign vorticity, on the vorticity fluxes are discussed. The results show that the background turbulence reduces the formation length and consequently increase the shedding frequency. The stronger decay of longitudinal vorticity flux in denser arrays is not accompanied by an increase of the lateral flux of vorticity. Furthermore, it was concluded that the decay of longitudinal vorticity flux is mainly caused by the vorticity cancellation due to the vorticity of opposite sign of close downstream cylinders.     

电磁力控制湍流边界层分离圆柱绕流场特性数值分析

在亚临界区高雷诺数Re=1.4×105下,采用脱体涡模拟结合湍流分离的方法对弱电解质中电磁力作用下湍流边界层分离圆柱绕流场及其升(阻)力特性进行了数值模拟和分析.结果表明,电磁力可以提高圆柱体湍流边界层内的流体动能,延缓圆柱体湍流边界层的流动分离,减弱圆柱体湍流绕流场中在流向和展向上大尺度漩涡的强度,减小圆柱体阻力时均值及其升力脉动幅值.当电磁力作用参数大于某个临界值后,湍流边界层流动分离消失,在圆柱体尾部产生射流现象,从而电磁力对圆柱体产生净推力作用,出现负阻力现象,而且升力脉动幅值接近于零,出现圆柱体升力消失现象.     

On the Scope of Lagrangian Vortex Methods for Two-Dimensional Flow Simulations and the POD Technique Application for Data Storing and Analyzing

The possibilities of applying the pure Lagrangian vortex methods of computational fluid dynamics to viscous incompressible flow simulations are considered in relation to various problem formulations. The modification of vortex methods—the Viscous Vortex Domain method—is used which is implemented in the VM2D code developed by the authors. Problems of flow simulation around airfoils with different shapes at various Reynolds numbers are considered: the Blasius problem, the flow around circular cylinders at different Reynolds numbers, the flow around a wing airfoil at the Reynolds numbers 104 and 105, the flow around two closely spaced circular cylinders and the flow around rectangular airfoils with a different chord to the thickness ratio. In addition, the problem of the internal flow modeling in the channel with a backward-facing step is considered. To store the results of the calculations, the POD technique is used, which, in addition, allows one to investigate the structure of the flow and obtain some additional information about the properties of flow regimes.     

Mean velocity behaviour in the near wake of long slender cylinder with a sharp trailing edge at a high Reynolds number

A simple eddy viscosity model is applied to the governing equations to establish the behaviour of the mean velocity in the turbulent axisymmetric near wake. The near wake develops from a long slender cylinder which is kept parallel to the flow and is developing under zero streamwise pressure gradient. The upstream turbulent boundary layer on the body of revolution is fully developed. In the inner layer of the flow downstream of the trailing edge, the turbulent inner layer of the upstream boundary layer grows into the initial logarithmic layer, and as a consequence, the centreline velocity in the near wake is shown to increase logarithmically with streamwise distances for large streamwise distances. The analysis further leads to two regions of the near wake flow (the inner near wake and the outer near wake), similar to that of a fully developed turbulent boundary layer, for which the governing equations have been derived. The matching between these two regions leads to a logarithmic variation in the normal direction. Also shown is the variation of the square of the wake width which varies logarithmically with streamwise distance in the near wake. These features are validated by comparison with available experimental data.     

Characterisation of the Eulerian and Lagrangian accelerations in the intermediate field of turbulent circular jets

A rich data-set of Lagrangian trajectories from 3D particle tracking velocimetry is used to study the structure of various acceleration components, vorticity, and strain in the intermediate field of a circular jet at Reynolds number Re = 6000. The total acceleration is decomposed into three distinctive sets: (1) streamwise–radial; (2) tangential–normal; and (3) local–convective components. Probability density function (PDF) and joint distributions of each set are characterised at various radial locations from the jet core within a streamwise band 16 ≤ x/d_h ≤ 17, where d_h is the diameter of the pipe. The PDF of the relative angle between the acceleration components and the velocity vector is also included to aid the characterisation. Results show that the acceleration components are described by two distinctive distributions: one of them exhibits symmetry and heavy tails, while the other is best fitted by a power-law type. The tails of acceleration PDFs are heavier with larger radial distance from the core. The increased departure from the Gaussian distribution with the distance from the core is a result of the increasing turbulence levels promoted by the mean shear. The variation of the third and fourth moments between the streamwise–tangential and the radial–normal accelerations indicate the anisotropy of the jet. Joint PDF of each acceleration decomposition exhibits distinctive distribution that appears to depend from the distance from the jet core. However, the vorticity and strain show similar PDF across radial distances. Finally, complementary analysis of a jet from a semicircular pipe shows the footprint of the nozzle geometry in the acceleration structure of jets.     

Particle‐in‐cell simulation of the effect of curved magnetic field on wall bombardment and erosion in a hall thruster

A curved, convex towards the channel bottom magnetic field is an important feature of an advanced Hall thruster that allows confining the plasma flow in the channel center, reducing the divergence angle of the ejected ion beam, and improving the discharge performance. In this article, the discharge behaviour of a Hall thruster in magnetic fields with different degrees of curvature is simulated with a particle‐in‐cell numerical method, and the effect of curved magnetic field on the ion bombardment and wall erosion and the associated mechanisms are studied and analysed. The results show that, as the curvature of the magnetic field increases, the propellant ionization becomes more confined at the channel center, the potential drop inside the channel decreases, and the acceleration region shifts outside the channel, which lead to the attenuation of the ion energy bombarding the wall and the deviation of the bombardment angle from the optimal sputtering angle. Conversely, the ion flux bombarding the wall near the channel exit increases. Nevertheless, the bombardment energy and angle are the dominant factors for the wall erosion, and the wall erosion rate clearly decreases with the increasing curvature of the magnetic field. These findings are closely related to the behaviour of electron conduction under a curved magnetic field; the relevant mechanisms are clarified in this article.     

Mechanism of the excitation of single pure mode L（0, 2） and its interaction with the defect in a hollow cylinder

Guided elastic waves have a great potential in pipe inspection as an efficient and low-cost nondestructive evaluation （NDE） technique, among which the wave of mode L（0, 2） receives a lot of attention because this mode is the fastest mode in a weakly dispersive region of frequency to minimize dispersion effects over a long distance and sensitive to the defects distributed circumferentially. Though many experimental and numerical researches have already been carried out about the excitation of L（0, 2） and its interaction with the defect in a hollow cylinder, its excitation mechanism has not been clarified yet. In this paper based on the transient response solution of the hollow cylinder, derived by the method of eigenfunction expansion, the theory about the exciting mechanism of mode L（0, 2） is advanced and the effects of the spatial distribution, vibration frequency and direction of the external force on the excitation are discussed. And the pure mode L（0, 2） is excited successfully under the parameters obtained through theoretical analysis. Furthermore, its interactions with some kinds of defects in hollow cylinders are simulated with the method of finite element analysis （FEA） and the results agree well with those obtained by other researchers.     

Theoretical and numerical calculation of the acoustic radiation force acting on a circular rigid cylinder near a flat wall in a standing wave excitation in an ideal fluid

The acoustic radiation force acting on a cylinder near a flat wall in a standing wave is calculated by analytical methods and numerical simulations. An exact theoretical solution is presented as well as an approximate solution. The approximate solution is in algebraic form and quite easy to compute. The numerical simulation is based on FVM (Finite Volume Method) on unstructured triangular meshes. The exact theoretical, approximate and numerical solutions are compared with each other and good agreements are obtained. Furthermore, the effects of the flat wall are investigated in detail by the three methods.     

First observation of natural circular dichroism spectra in the extreme ultraviolet region using a polarizing undulator‐based optical system and its polarization characteristics

Natural circular dichroism (CD) spectra in the extreme ultraviolet (EUV) region down to a wavelength of 80 nm have been observed for the first time, using an alanine thin film deposited on sodium salicylate coated glass as a sample. Calibrated EUV‐CD spectra of l ‐alanine exhibited a large negative peak at around 120 nm and a positive CD signal below 90 nm, which were roughly predicted by theoretical calculations. A CD measurement system with an Onuki‐type polarizing undulator was used to obtain the EUV‐CD spectra. This CD system, the development of which took five years, can be used to observe even weak natural CD spectra. The polarization characteristics of this system were also evaluated in order to calibrate the recorded CD spectra.