考虑风剪切的漂浮式风力机气动特性研究

波浪与风载的共同作用下,海上风电机组平台存在多自由度运动,同时伴随着海上复杂的风况,其气动特性变化较为复杂。以NREL5MW风电机组为研究对象,在风剪切来流下,依据波浪和风载的作用规律,研究平台纵摇和纵荡运动对机组气动特性和绕流场细节的影响。结果表明:平台运动会造成风电机组气动性能的周期性波动,而风剪切作用使得风电机组平均发电量降低并加剧功率和推力的波动。风剪切会增大纵摇运动过程中展向截面的法向载荷幅值且波动加剧,但会降低与功率输出有关的切向受力;风剪切作用推迟了纵荡运动过程中展向截面法向和切向载荷峰值出现时刻,加剧载荷的波动。外叶展截面的法向和切向力系数曲线出现平台,加剧了叶片疲劳载荷,减小了功率输出。     

Hydrodynamic binary coalescence of droplets under air flow in a hydrophobic microchannel

Based on the volume of fluid(VOF) method, we conduct a numerical simulation to study the hydrodynamic binary coalescence of droplets under air flow in a hydrophobic rectangular microchannel. Two distinct regimes, coalescence followed by sliding motion and that followed by detaching motion, are identified and discussed. Additionally, the detailed hydrodynamic information behind the binary coalescence is provided, based on which a dynamic mechanical analysis is conducted to reveal the hydrodynamic mechanisms underlying these two regimes. The simulation results indicate that the sliding motion of droplets is driven by the drag force and restrained by the adhesion force induced by the interfacial tension along the main flow direction. The detachment(i.e., upward motion) of the droplet is driven by the lift force associated with an aerodynamic lifting pressure difference imposed on the coalescent droplet, and also restrained by the adhesion force perpendicular to the main flow direction. Especially, the lift force is mainly induced by an aerodynamic lifting pressure difference imposed on the coalescent droplet. Two typical regimes can be quantitatively recognized by a regime diagram depending on Re and We. The higher Re and We respectively lead to relatively larger lift forces and smaller adhesion forces acting on the droplet, both of which are helpful to detachment of the coalesced droplet.     

Experimental estimation of the electric force induced by a blade-plane actuator in dielectric liquids

A blade-plane actuator immersed in a dielectric liquid and connected to a high voltage power supply produces a Coulomb force which sets the liquid in motion from the blade tip to the plane. Jets of more than 1 m/s have already been observed. In this work, a method is proposed to estimate the electric force generated by a blade plane actuator from Particle Image Velocimetry (PIV) velocity field measurements. The originality of this paper comes from the fact that the volume force is not measured directly with an aerodynamic shielded balance but calculated from velocity fields obtained by Particle Image Velocimetry. In this article the global time-averaged electric force is computed in a fixed control volume by the use of the momentum equation in its integral form.     

Nonlinear analysis and enhancement of wing-based piezoaeroelastic energy harvesters

We investigate the level of harvested power from aeroelastic vibrations for an elastically mounted wing supported by nonlinear springs. The energy is harvested by attaching a piezoelectric transducer to the plunge degree of freedom. The considered wing has a low-aspect ratio and hence three dimensional aerodynamic effects cannot be neglected. To this end, the three dimensional unsteady vortex lattice method for the prediction of the unsteady aerodynamic loads is developed. A strong coupling scheme that is based on Hamming's fourth-order predictor–corrector method and accounts for the interaction between the aerodynamic loads and the motion of the wing is employed. The effects of the electrical load resistance, nonlinear torsional spring and eccentricity between the elastic axis and the gravity axis on the level of the harvested power, pitch and plunge amplitudes are investigated for a range of operating wind speeds. The results show that there is a specific wind speed beyond which the pitch motion does not pick any further energy from the incident flow. As such, the displacement in the plunge direction grows significantly and causes enhanced energy harvesting. The results also show that the nonlinear torsional spring plays an important role in enhancing the level of the harvested power. Furthermore, the harvested power can be increased by an order of magnitude by properly choosing the eccentricity and the load resistance. This analysis is helpful in designing piezoaeroelastic energy harvesters that can operate optimally at specific wind speeds.     

The quasiclassical Langevin equation and its application to the decay of a metastable state and to quantum fluctuations

We report on investigations on the consequences of the quasiclassical Langevin equation. This Langevin equation is an equation of motion of the classical type where, however, the stochastic Langevin force is correlated according to the quantum form of the dissipation-fluctuation theorem such that ultimately its power spectrum increases linearly with frequency. Most extensively, we have studied the decay of a metastable state driven by a stochastic force. For a particular type of potential well (piecewise parabolic), we have derived explicit expressions for the decay rate for an arbitrary power spectrum of the stochastic force. We have found that the quasiclassical Langevin equation leads to decay rates which are physically meaningful only within a very restricted range. We have also studied the influence of quantum fluctuations on a predominantly deterministic motion and we have found that there the predictions of the quasiclassical Langevin equations are correct.     

Edison revisited: Electro mechanical effects in wet porous materials

Edison discovered that the coefficient of friction between a metallic plate and a porous material moistened with a dilute electrolyte could be modulated by an electric field. In experiments on the same kind of contacts (clays or chalks on carbon or metals) but without continuous tangential relative motion we measure two electro-mechanical effects at frequencies of the order of 10 kHz. An alternating field induces an alternating normal force between the porous material and the conducting base. The force is lagging by versus the field. A forced normal relative motion induces through the contact a current nearly in phase with the motion. For an explanation we start from Helmholtz theory of stationary electrophoretic phenomena. We present a model in which liquid motions obey the Helmholtz laws. It explains decently the phase relations between causes and effects, and approximately the values of the effects. In optical experiments on contacts between a wet clay and the transparent conducting coating of a glass plate we measure in the frequency range 1-100 kHz a modulation of reflecting power induced by an alternating potential. The decrease of reflecting power is lagging by an angle close to behind the field. We think the modulation observed is induced by a modulation of the amount of liquid in the film present between glass and clay. In friction experiments this alternating liquid lubrication acting exclusively at the very places where friction occurs may have significant effect. Received: 22 July 1997 / Revised: 21 October 1997 / Accepted: 21 November 1997     

基于SPH方法的高速列车风阻制动板制动力研究

采用光滑粒子流体动力学（smoothed particle hydrodynamics,SPH）方法建立了高速列车风阻制动板流固耦合二维数值模型,实现了对风阻制动板动态开启过程的数值模拟,分析了风阻制动板的运动规律、受力分布和气动阻力变化规律.结果表明:当车速为300 km/h,最大开启角度为90°时,制动板整个开启过程的动作时间约为0.025 s,满足紧急制动时的快速响应要求;制动板开启初期,板的受力有较大突变且呈梯度分布,易弯曲变形;制动板的整个开启过程中,制动板整体受力水平逐渐提高且保持均匀分布,制动板直角边缘处受力始终较其他部位大,具有良好的增阻效果;最大制动力出现在最大开启角度为75°~85°时,而非垂直开启时,建议将最大开启角度设置为80°以获得最大制动力.      

Cutting process dynamics by nonlinear time series and wavelet analysis

We have modeled the dynamics of a cutting process by a two-degree-of-freedom mass-spring system with dry friction. Using nonlinear time series and wavelet analysis, we have investigated the vibrational instabilities of the system for different values of the cutting force. By constructing the phase portraits and calculating the Lyapunov exponents we have delineated the conditions for which a periodic or chaotic motion can occur. The results are verified by means of a time-scale representation of the wavelet power spectrum.     

Spin-dependent forces on trapped ions for phase-stable quantum gates and entangled states of spin and motion

Favored schemes for trapped-ion quantum logic gates use bichromatic laser fields to couple internal qubit states with external motion through a "spin-dependent force." We introduce a new degree of freedom in this coupling that reduces its sensitivity to phase decoherence. We demonstrate bichromatic spin-dependent forces on a single trapped 111Cd+ ion, and show that phase coherence of the resulting entangled states of spin and motion depends critically upon the spectral arrangement of the optical fields. This applies directly to the operation of entangling gates on multiple ions.     

Gravitomagnetism and Non-commutative Geometry

Similarity between the gravitoelectromagnetism and the electromagnetism is discussed. We show that the gravitomagnetic field (similar to the magnetic field) can be equivalent to the non-commutative effect of the momentum sector of the phase space when one maintains only the first order of the non-commutative parameters. This is performed through two approaches. In one approach, by employing the Feynman proof, the existence of a Lorentz-like force in the gravitoelectromagnetism is indicated. The appearance of such a force is subjected to the slow motion and the weak field approximations for stationary fields. The analogy between this Lorentz-like force and the motion equation of a test particle in a non-commutative space leads to the mentioned equivalency. In fact, this equivalency is achieved by the comparison of the two motion equations. In the other and quietly independent approach, we demonstrate that a gravitomagnetic background can be treated as a Dirac constraint. That is, the gravitoelectromagnetic field can be regarded as a constrained system from the sense of the Dirac theory. Indeed, the application of the Dirac formalism for the gravitoelectromagnetic field reveals that the phase space coordinates have non-commutative structure from the view of the Dirac bracket. Particularly, the gravitomagnetic field as a weak field induces the non-trivial Dirac bracket of the momentum sector which displays the non-commutativity.     

Analytical prediction of the non-linear response of a self-excited structure

The time-dependent amplitude response of the self-excited harmonium reed vibrating at finite amplitudes is investigated both analytically and experimentally. The analysis contained herein is based on the assumption that all of the significant non-linear forces that act on the reed are of aerodynamic origin and that these forces influence the reed behavior through the system damping. The analysis is carried out, without the aid of empirical techniques, by deriving the induced aerodynamic pressure force as a non-linear function of the amplitude of the reed motion, via an unsteady potential flow field analysis, and then applying the resulting forcing function to the equation of motion of the reed. An approximate solution of the equivalent lumped parameter equation of motion yields functional relationships that predict various observable phenomena, including the limit cycle amplitude. A comparison of the results with experimental data serves to substantiate the analysis. Also presented is a brief discussion of some of the phenomenological similarities that exist between the results of the harmonium reed analysis and the observed behavior of a flat-plate wing undergoing torsional stall flutter.     

Spin force and torque in non-relativistic Dirac oscillator on a sphere

The spin force operator on a non-relativistic Dirac oscillator (in the non-relativistic limit the Dirac oscillator is a spin one-half 3D harmonic oscillator with strong spin–orbit interaction) is derived using the Heisenberg equations of motion and is seen to be formally similar to the force by the electromagnetic field on a moving charged particle. When confined to a sphere of radius R, it is shown that the Hamiltonian of this non-relativistic oscillator can be expressed as a mere kinetic energy operator with an anomalous part. As a result, the power by the spin force and torque operators in this case are seen to vanish. The spin force operator on the sphere is calculated explicitly and its torque is shown to be equal to the rate of change of the kinetic orbital angular momentum operator, again with an anomalous part. This, along with the conservation of the total angular momentum, suggests that the spin force exerts a spin-dependent torque on the kinetic orbital angular momentum operator in order to conserve total angular momentum. The presence of an anomalous spin part in the kinetic orbital angular momentum operator gives rise to an oscillatory behavior similar to the Zitterbewegung. It is suggested that the underlying physics that gives rise to the spin force and the Zitterbewegung is one and the same in NRDO and in systems that manifest spin Hall effect.     

双三角翼非定常俯仰运动实验与数值模拟

文章针对双三角翼大振幅正弦俯仰运动过程中的非定常载荷和流动特性开展了实验与数值模拟研究,并与相同主翼后掠角的单三角翼进行了对比.实验研究在低速回流式水槽中开展,所采用的实验模型为边条后掠角为75°,主翼后掠角为50°的双三角翼全模,俯仰运动的旋转轴位于主翼弦长的2/3处,振幅为0~60°,运动的缩减频率k=0.03,0.06,0.12,0.24,0.48.实验Reynolds数以主翼弦长为参考Re=1.69×104.在水槽的测力实验中,发现非定常流动力的迟滞现象,并且随着非定常运动缩减频率的增大,流动的迟滞效应也随之增大.与相同主翼后掠角的单三角翼相比,双三角翼的迟滞环在低缩减频率下更小,但随着缩减频率的增大,这种差距逐渐减小.在数值模拟研究中,采用DDES湍流模型对俯仰双三角翼的流场进行了数值模拟.流场结果表明,在较低的缩减频率下,主翼吸力面的前缘涡是影响气动力的主要因素,非定常流动力的迟滞效应主要与前缘涡在上仰过程中的延迟破裂和下俯过程中的延迟恢复有关;在较高的缩减频率下,机翼前缘涡对气动力的影响减小,由机翼俯仰角速度而产生的环量力成为了气动力的主导因素,因此在较高缩减频率下,单三角翼与双三角翼的升力特性趋于一致.      

光学炮塔气动载荷的非稳态特性

采用耦合J-B模型的IDDES模型与双时间步LU-SGS方法开展炮塔非稳态气动载荷的数值仿真研究。炮塔流动会发生分离,天顶位置的分离角大于90°;当流动绕过炮塔时,形成马蹄涡、脱落涡街等非稳态流场结构,导致气动载荷也具有非稳态特性;炮塔顶点的脉动静压功率谱在1.6~40.0 kHz进入各向同性均匀湍流的惯性子区,基本满足Kolmogrov的-5/3定律;气动力以阻力为主,横向力的脉动幅值大,气动力矩则以俯仰力矩为主,滚转力矩的脉动幅值大,偏航力矩可以忽略不计;气动力和力矩的功率谱主要集中在1 kHz以下,存在多个尖峰频率,主频约为230 Hz (斯特劳哈尔数为0.15)。在ATP系统设计之初,需要考虑光学炮塔所受气动载荷的非稳态特性,并规避尖峰频率尤其是主频的谐振破坏问题。     

气温和颗粒密度对声场中颗粒动力学影响的数值模拟

综合考虑黏性夹带力、Basset力、虚拟质量力和压力梯度力,建立颗粒在声场中的动力学模型,利用变步长四阶RungeKutta算法和二阶隐式Adams插值算法对颗粒的受力和运动进行数值模拟。将模拟和实验得到的颗粒运动特性进行对比,验证数值模拟的正确性。在此基础上,研究气温和颗粒密度对颗粒动力学的影响规律。结果表明,黏性夹带力对颗粒运动起主导作用;气温升高,压力梯度力与黏性夹带力之间的相位差减小,Basset力、虚拟质量力与黏性夹带力之间的相位差增大。研究还发现,气温较低时,颗粒密度对颗粒运动有重要影响,夹带系数随着密度的增加而迅速下降;气温较高时,颗粒密度对颗粒运动的影响较小,颗粒位移振幅和夹带系数相对低温时明显增加。      

Magneto-elastic dynamics and bifurcation of rotating annular plate

In this paper, magneto-elastic dynamic behavior, bifurcation, and chaos of a rotating annular thin plate with various boundary conditions are investigated. Based on the thin plate theory and the Maxwell equations, the magneto-elastic dynamic equations of rotating annular plate are derived by means of Hamilton's principle. Bessel function as a mode shape function and the Galerkin method are used to achieve the transverse vibration differential equation of the rotating annular plate with different boundary conditions. By numerical analysis, the bifurcation diagrams with magnetic induction, amplitude and frequency of transverse excitation force as the control parameters are respectively plotted under different boundary conditions such as clamped supported sides, simply supported sides, and clamped-one-side combined with simply-anotherside. Poincare′ maps, time history charts, power spectrum charts, and phase diagrams are obtained under certain conditions,and the influence of the bifurcation parameters on the bifurcation and chaos of the system is discussed. The results show that the motion of the system is a complicated and repeated process from multi-periodic motion to quasi-period motion to chaotic motion, which is accompanied by intermittent chaos, when the bifurcation parameters change. If the amplitude of transverse excitation force is bigger or magnetic induction intensity is smaller or boundary constraints level is lower, the system can be more prone to chaos.     

Investigation of a hybrid approach combining experimental tests and numerical simulations to study vortex-induced vibration in a circular cylinder

In this study, a hybrid approach based on computational fluid dynamics (CFD) was used to investigate the aerodynamic forces associated with vortex-induced vibration (VIV) in a circular cylinder. The circular cylinder and the flow field were considered as two substructures of a system. Circular cylinder motion was produced in a wind tunnel test of the VIV prior to the numerical simulation; this motion was used as a known cylinder boundary condition and applied to the flow field. The flow field with the known moving boundary condition was then numerically simulated by the ANSYS CFX code. The transient aerodynamic coefficients of the circular cylinder with predetermined motion were obtained from the numerical simulation. To verify the feasibility and accuracy of the proposed hybrid approach and to calculate cylinder vibrations, the transient aerodynamic coefficients were applied to a single degree of freedom (SDOF) model of the circular cylinder. The oscillation responses of the circular cylinder from the calculated (SDOF model) and experimental results were compared, and the results indicate that the hybrid approach accurately simulated the transient aerodynamic coefficients of the circular cylinder. For further comparison, a nonlinear aerodynamic coefficient model based on a nonlinear least square technique was applied to the SDOF model. The nonlinear aerodynamic model can predict well the amplitude and lock-in region of the VIV of the circular cylinder model.     

Dynamical phase transition of a driven vortex lattice with disordered pinning

We have numerically solved the overdamped equation of vortex motion in a two-dimensional driven vortex lattice with disordered pinning, in which the driving Lorentz force, the pinning force due to point defects, the intervortex interacting force, and the thermal fluctuation force are taken into account. It is found that the vortex density and pinning strength are two important factors of affecting the melting transition of a vortex lattice. At low magnetic fields, there exist hysteresis loops of the average vortex velocity and the average pinning force vs. the driving force, from which the feature of a first-order melting transition of the vortex motion can be clearly seen. As the magnetic field is increased beyond a critical value, the hysteresis loops disappear and the melting transition is replaced by a second-order glass transition. We have also studied the influence of intervortex interactions on the vortex melting transition by comparing several forms of repulsive forces between the vortices.     

The flow induced by a jellyfish

The purpose of this study is to understand the propulsion mechanism of a jellyfish during its swimming. We observed the motion of a jellyfish (Aurelia aurita) by a motion-capture camera, and measured the vector field of flow around a jellyfish by using a PIV (Particle Image Velocimetry) measurement. A jellyfish is considered to be principally propelled by a jet at the contracting phase of its motion. If that is true, it is interesting that a jellyfish never stops traveling even at the expanding phase. We found that a vortex ring with the opposite vorticity to shed vortex ring was inside a jellyfish body in the expanding phase. We discussed a cause of an increase in thrust force and keeping constant speed in the expanding phase.     

Wall effect on fluid–structure interactions of a tethered bluff body

Wind tunnel experiments have shown an unexplained amplification of the free motion of a tethered bluff body in a small wind tunnel relative to that in a large wind tunnel. The influence of wall proximity on fluid–structure interaction is explored using a compound pendulum motion in the plane orthogonal to a steady freestream with a doublet model for aerodynamic forces. Wall proximity amplifies a purely symmetric single degree of freedom oscillation with the addition of an out-of-phase force. The success of this simple level of simulation enables progress to develop metrics for unsteady wall interference in dynamic testing of tethered bluff bodies.