In bio-inspired flapping-wing flight,lift generation and flexible deformation are intrinsically coupled.Thus,an experimental study is conducted to reveal the fluid-structure interaction mechanism for a flexible plate ...In bio-inspired flapping-wing flight,lift generation and flexible deformation are intrinsically coupled.Thus,an experimental study is conducted to reveal the fluid-structure interaction mechanism for a flexible plate undergoing pitching and plunging motion,and theoretical models are proposed to predict either lift or deformation based on a series of simplifying assumptions.It is indicated that flexible plates can effectively reduce the amplitude of lift and pitching moment coefficients during dynamic stall,with increased flexibility leading to higher load reduction.To investigate the effect of plate flexibility on lift,a definition of effective angle of attack is proposed,incorporating the pitching and plunging motion and chordwise deformation of the flexible plate,which can reduce the hysteresis effect of the lift coefficient during dynamic stall.As a consequence,a theoretical model is developed to predict lift based on observed motion and deformation.On the other hand,another theoretical model is developed to predict flexible-plate deformation utilizing aerodynamic forces,revealing the effect of leading-edge vortex evolution on passive deformation.The influence of kinematic parameters,including the maximum effective angle of attack,reduced frequency,and Strouhal number,on the aerodynamic forces is further studied.Compared with the rigid plate,flexible plates exhibit lower sensitivity of aerodynamic forces to changes in kinematic parameters due to their inherent compliance and resulting deformation.The proposed theoretical models can serve as a reference for aerodynamic and deformation prediction in bio-inspired flexible structures.展开更多
While spin-orbit interaction has been extensively studied,few investigations have reported on the interaction between orbital angular momenta(OAMs).In this work,we study a new type of orbit-orbit coupling between the ...While spin-orbit interaction has been extensively studied,few investigations have reported on the interaction between orbital angular momenta(OAMs).In this work,we study a new type of orbit-orbit coupling between the longitudinal OAM and the transverse OAM carried by a three-dimensional(3D)spatiotemporal optical vortex(STOV)in the process of tight focusing.The 3D STOV possesses orthogonal OAMs in the x-y,t-x,and y-t planes,and is preconditioned to overcome the spatiotemporal astigmatism effect.x,y,and t are the axes in the spatiotemporal domain.The corresponding focused wavepacket is calculated by employing the Debye diffraction theory,showing that a phase singularity ring is generated by the interactions among the transverse and longitudinal vortices in the highly confined STOV.The Fourier-transform decomposition of the Debye integral is employed to analyze the mechanism of the orbit-orbit interaction.This is the first revelation of coupling between the longitudinal OAM and the transverse OAM,paving the way for potential applications in optical trapping,laser machining,nonlinear light-matter interactions,and more.展开更多
The vortex interference mechanism on low Reynolds number between the canard and main wing of the canard-forward sweep wing （Canard-FSW） configurations is simulated numerically by employing the numerical wind tunnel ...The vortex interference mechanism on low Reynolds number between the canard and main wing of the canard-forward sweep wing （Canard-FSW） configurations is simulated numerically by employing the numerical wind tunnel method. The variations of aerodynamic characteristics of Canard-FSW configurations with different positions of the canard are investigated,finding that the aerodynamic interference and mutual coupling effect between the canard and main wing have made great contributions to the lift and stability characteristics of the whole aircraft. Canard can radically improve the surface flow pattern of the main wing. And its own vortex can have a favorable interference on the main wing and can effectively control the airflow boundary layer separation. At small angles of attack,the aerodynamic characteristics are sensitive to the positions of the canard and the main wing,but at high angles of attack,the aerodynamic performances of the configuration are not only related to the shape of the canard （forward or backward）,but also with the size of control force as well as the features of the vortices generated above the main wing and the canard. The different configurations and vortices are illustrated using the velocity vector,streamlines and pressure contours.展开更多
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. I...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.展开更多
Observations are presented from experiments and calculations where a laminar spherical CH4/air flame is perturbed successively by incident and reflected shock waves. The experiments are performed in a standard shock t...Observations are presented from experiments and calculations where a laminar spherical CH4/air flame is perturbed successively by incident and reflected shock waves. The experiments are performed in a standard shock tube arrangement, in which a high-speed shadowgraph imaging system is used to record evolutions of the flame. Numerical simulations are conducted by using second-order wave propagation algorithms, based on two-dimensional axisymmetric Navier-Stokes equations with detailed chemical reactions. Qualitative agreements are obtained between the experimental and numerical results. Under actions of incident shock waves, Richtmyer-Meshkov instability responsible for the flame deformation is induced in the flame, and the distoned flame takes a barrel shape. Then, under subsequent actions of the shock wave reflected from a planar wall, the flame takes an inclined non-symmetrical kidney shape in a symmetric cross section, which means a mushroom-like shape of the flame comes finally into being. The vorticity direction in the ring cap has been altered by the reflected shock＇s action, which makes the head of the mushroom-like flame extend quickly to the side wall.展开更多
In order to understand the interaction between large-scale vortex structure and particles, a two-way coupling temporal mixing layer laden with particles at a Stokes number of 5 with different mass loading planted init...In order to understand the interaction between large-scale vortex structure and particles, a two-way coupling temporal mixing layer laden with particles at a Stokes number of 5 with different mass loading planted initially in the upper half region is numerically studied. The pseudospectral method is used for the flow fluid and the Lagrangian approach is employed to trace particles. The momentum coupling effect introduced by a particle is approximated to a point force. The simulation results show that the coherent structures are still dominant in the mixing layer, but the large-scale vortex structure and particle dispersion are modulated. The length of large-scale vortex structure is shortened and the pairing is delayed. At the same time, the particles are distributed more evenly in the whole flow field as the mass loading is increased, but the particle dispersion along the transverse direction differs from that along the spanwise direction, which indicates that the effect by the addition of particle on the spanwise large-scale vortex structure is different from the streamwise counterpart.展开更多
Numerical simulation methods for unsteady vortex field of helicopter rotor with high resolution and low dissipation TENO8-AA primitive variables reconstruction schemes are established based on moving-embedded grid and...Numerical simulation methods for unsteady vortex field of helicopter rotor with high resolution and low dissipation TENO8-AA primitive variables reconstruction schemes are established based on moving-embedded grid and Navier-Stokes equations.Firstly,the Targeted Essentially Non-Oscillatory(TENO)scheme are developed by employing ENO-like candidate stencil selection strategy,and the candidate stencil is adopted with optimal weight in smooth region while it is discarded completely in discontinuous region,which reduces the dissipation and dispersion errors and approaches better spectral properties.Then,the aerodynamic characteristics of Helishape-7A model rotor in Blade Vortex Interaction(BVI)state and the flowfield of Lynx rotor in hover are simulated,which validates that the blade tip vortex trajectory with larger wake age and more details of vortex can be captured by TENO8-AA scheme with only a quarter of grid points and half time comparing to WENO-JS scheme.Moreover,the simulation accuracy of thrust coefficient is improved by up to 36%.Finally,the analyses for BVI and aeroacoustic characteristics of Operational Loads Survey(OLS)rotor are conducted,and the different forms of interaction mechanism are explored,such as oblique and parallel interactions.The results indicate that TENO scheme not only ensures the resolution of simulation in discontinuous region,but also minimizes the numerical dissipation in smooth region dominated by blade tip vortex.Therefore,the acoustic pressure peak prediction error of rotor in BVI state is significantly reduced to 5.6%and 0.8%at two microphone locations,respectively.展开更多
The effects of inlet gas parameters and sloping sidewall angle on the flame structure and combustion limit with and without sidewall were experimentally investigated.Flame height and impact angle were obtained by chem...The effects of inlet gas parameters and sloping sidewall angle on the flame structure and combustion limit with and without sidewall were experimentally investigated.Flame height and impact angle were obtained by chemiluminescence intensity analysis of CH*distribution.First,the combustion characteristics of flame with and without sidewall at different equivalence ratios were explored;then,the influence of Reynolds number and inlet gas temperature on flame structure and combustion limit of v-shaped flame with sidewall were analyzed,and the results with sidewall were compared with those without sidewall.Finally,the variation trend of flame parameters with different sloping sidewall angles was analyzed.The experimental results show that the existence of sidewall makes flame shape change from“M-shaped”to“inverted N-shaped”,and conical shape to trapezoidal shape.The inhibition effect of sidewall on flame stretching downstream is strengthened with the increase in Reynolds number;but as the temperature of the inlet gas increases,the inhibitory effect is obviously weakened.When sloping sidewall angle decreases from 90°to 55°at 5°intervals,flame height and impact angle of v-shaped flame reach the extreme value whenβ=80°.Compared with the case without sidewall,the range of v-shaped flame with sidewall has no obvious trend of broadening or shrinking when inlet gas temperature is increased;however,as sloping sidewall angle decreases,the range of the v-shaped flame shrinks obviously and flammability limit increases significantly.展开更多
Based on 6-hourly sensible heat flux and latent heat flux from the NCEP Climate Forecast System Reanalysis(CFSR) and circulation data from the Japanese 25-year Reanalysis(JRA-25),the initial developing process of trop...Based on 6-hourly sensible heat flux and latent heat flux from the NCEP Climate Forecast System Reanalysis(CFSR) and circulation data from the Japanese 25-year Reanalysis(JRA-25),the initial developing process of tropical cyclone Mindulle(1005) in 2010 has been diagnosed to reveal the impact of air-sea interaction over the South China Sea(SCS) on the genesis of its incipient vortex.The results show that the incipient vortex first occurred east of the Luzon Island on 0000 UTC 20 August,suggesting that the topographic forcing of the Luzon Island for easterly winds over the western Pacific might be one of the factors responsible for the formation of the incipient vortex.During the formation stage of the incipient vortex,strong southeasterlies over the SCS caused warm water of the middle and eastern SCS to flow toward the Luzon Island due to Ekman transport resulting from wind stress,leading to an increase of the sea surface temperature and sensible heat flux into the atmosphere.Although the anomalous sensible heating favored surface pressure to reduce,it was not conducive to the increase of local vorticity associated with the vortex above the heating area because,according to the atmospheric thermal adaptation theory,the anticyclonic vorticity would be created in the lower troposphere due to the decreased vertical gradient of the sensible heating.However,the ascending motions occurred over the eastern area of the anomalous sensible heating due to the augmentation of the vorticity advection with increasing height,causing water vapor to condense in the middle and upper troposphere.In turn,cyclonic vorticity was generated in the lower troposphere due to the increased vertical gradient of the condensation latent heating,resulting in the formation and further growth of the incipient vortex.Therefore,the vorticity creation due to the condensation heating played a dominant role during the subsequent enhancing stage of the incipient vortex.展开更多
When a Bose-Einstein condensate is set to rotate, superfluid vortices will be formed, which finally condense into a vortex lattice as the rotation frequency further increases. We show that the dipole-dipole interactio...When a Bose-Einstein condensate is set to rotate, superfluid vortices will be formed, which finally condense into a vortex lattice as the rotation frequency further increases. We show that the dipole-dipole interactions renormalize the short-range interaction strength and result in a distinction between interactions of parallel-polarized atoms and interactions of antiparallel-polarized atoms. This effect may lead to a spontaneous breakdown of the rapidly rotating Bose condensate into a novel anti-ferromagnetic-like vortex lattice. The upward-polarized Bose condensate forms a vortex lattice, which is staggered against a downward-polarized vortex lattice. A phase diagram related to the coupling strength is obtained.展开更多
The underlying effect of vortex interaction characterized by the merging and non-merging on mixing enhancement is of fundamental significance to understand the flow dynamics of strut injectors in scramjets.Starting fr...The underlying effect of vortex interaction characterized by the merging and non-merging on mixing enhancement is of fundamental significance to understand the flow dynamics of strut injectors in scramjets.Starting from a simplified configuration of a vortex generator,this study focuses on the influence of geometric parameters on vortex structures and fluid mixing through compressible Navier-Stokes(NS)simulations.By adjusting the induction of outer vortices,the inner co-rotating vortex pair exhibits two modes of interaction(merging/separation regime)reflected by closer/farther vortex centers.Defined by the zero variation rate of the inner vortex spacing,the critical state of equilibrium is determined.The critical condition is well predicted by a theoretical model based on the Biot-Savart law.Through the introduction of mixedness and mixing time,the intrinsic impact of interaction modes on fluid mixing is revealed.Compared with the vortex dynamics in the merging regime,the one in the separation regime is more effective for passive scalar mixing augmentation.With efficient material stretching characterized by the higher interface stretching factor and averaging finite-time Lyapunov exponent(FTLE),the mixing time is shortened by as much as 2.5 times in the separation regime.The implication of the present two interaction regimes in mixing enhancement physically reflected by the averaging FTLE has the potential to improve the combustion performance and shorten the combustor chamber.展开更多
The dynamic properties of interacting vortex-antivortex pairs in thin film are studied by analytical calculations. An- alytical expressions for the magnetization vector distribution of vortex-antivortex pairs and the ...The dynamic properties of interacting vortex-antivortex pairs in thin film are studied by analytical calculations. An- alytical expressions for the magnetization vector distribution of vortex-antivortex pairs and the trivortex states are given. The magnetic states of the vortices are treated as having rigid structures, i.e., the vortex maintains its spin distribution when moving. The trajectories of the vortex cores are calculated by the Thiele＇s equation. It is found that the vortex-antivortex pair rotates around each other when they have opposite polarities, however, vortex and antivortex cores move along straight lines when they have the same polarity. The frequency of the rotation decreases with increasing the distance between the two cores of vortex-antivortex pair, and it has a lower value when a third vortex is introduced.展开更多
The influences due to several AVIs （airfoil-vortex interactions） are studied by using a two-dimensional CFD （computational fluid dynamics） method. The primary goal of this effort is to assess the variation of vort...The influences due to several AVIs （airfoil-vortex interactions） are studied by using a two-dimensional CFD （computational fluid dynamics） method. The primary goal of this effort is to assess the variation of vortex center location and vortex circulation associated with sequential AVI toward an improvement of the hybrid method of CFD and prescribed wake model, which closely relates to predicting the BVI （blade-vortex interaction） noise radiated from a helicopter rotor. The representative of sequential AVI is performed by single vortex and two airfoils. Investigations with respect to vortex center location and vortex circulation after AVIs have been made by varying the miss-distance, which is the vertical distance between the airfoil leading edge and the vortex center. Correlations between miss-distance and vorticity field show that there exists complicated vortex wake flow with several vortices newly induced in 1st AVI. The pressure fluctuation amplitude clarifies that the intensity in 2nd AV1 is significantly low compared to the intensity in 1st AVI due to the influence of vortex dissipation. Simulations turned out to modify the vortex center location represented by the hybrid method using an offset value for a streamwise direction and to dissipate the vortex circulation for improving the accuracy of BVI noise prediction.展开更多
Large eddy simulation is performed to study three-dimensional wave-current interaction with a square cylinder at different Reynolds numbers, ranging from 1,000 to 600,000. The Keulegan-Carpenter number is relevantly a...Large eddy simulation is performed to study three-dimensional wave-current interaction with a square cylinder at different Reynolds numbers, ranging from 1,000 to 600,000. The Keulegan-Carpenter number is relevantly a constant of 0.6 for all cases. The Strouhal number, the mean and the RMS values of the effective drag coefficient in the streamwise and transverse directions are computed for various Reynolds numbers, and the velocity of a rep- resentative point in the turbulent zone is simulated to find the turbulent feature. It is found that the wave-current interaction should be considered as three-dimensional flow when the Reynolds number is high; under wave-current effect, there exists a critical Reynolds number, and when the Reynolds number is smaller than the critical one, current effect on wave can be nearly neglected; conversely, with the Reynolds number increasing, wave-currentstructure interaction is sensitive to the Reynolds number.展开更多
This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier...This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier-Stokes equation is discretized spatially with collocated finite volume method and Eulerian implicit method in time domain. The hybrid method that combines immersed boundary method （IBM） and volume of fluid （VOF） method is used to deal with rigid body motion in fluid domain. The details of movement of immersed boundary （IB） and calculation of VOF are also described. This method can be easily applied to any existing finite-volume-based computational fluid dynamics （CFD） solver without complex operation, with which fluid flow interaction of arbitrarily complex geometry can be realized on a fixed mesh. The method is verified by low Reynolds number flows passing both stationary and oscillating cylinders. The drag and lift coefficients acquired by the study well accord with other published results, which indicate the reasonability of the proposed method.展开更多
This paper studies in detail the interaction of two edge dislocations nested in a Gaussian beam propagating in free space. It shows that in free-space propagation the edge dislocations are unstable and vanish, and two...This paper studies in detail the interaction of two edge dislocations nested in a Gaussian beam propagating in free space. It shows that in free-space propagation the edge dislocations are unstable and vanish, and two noncanonical vortices with opposite topological charge take place when off-axis distances cl and c2 of two edge dislocations are nonzero, and the condition k2w08＋ 32c1c2（w02- 2C1C2）Z2 〉 0 is fulfilled （k-wave number, w0-waist width）. A noncanonical vortex appears when one off-axis distance is zero. However, one edge dislocation is stable when two edge dislocations are perpendicular and one off-axis distance is zero. Two perpendicular edge dislocations both with zero off-axis distance are also stable. The analytical results are illustrated by numerical examples.展开更多
Effects of Reynolds(Re)number and Schmidt(Sc)number on the flow structures and variable density mixing are numerically investigated through the canonical shock cylindrical bubble interaction.By determining the viscosi...Effects of Reynolds(Re)number and Schmidt(Sc)number on the flow structures and variable density mixing are numerically investigated through the canonical shock cylindrical bubble interaction.By determining the viscosity and diffusivity within a wide range,the controlling parameters,total vortex circulation,and compression rate,are conservative under a broad range of Re and Sc numbers(Re≈10^(3)-10^(5)and Sc≈0.1-5)in the same shock Mach(Ma)number condition(Ma=2.4).As for the Re number effect,the circulation of secondary baroclinic vorticity(SBV),induced by the main vortex centripetal acceleration,is observed to be higher in high Re number and vice versa.Based on the vorticity transport equation decomposition,a growth-inhibition vorticity dynamics balance mechanism is revealed:the vorticity viscous term grows synchronously with baroclinic production to inhibit SBV production in low Re number.By contrast,the viscous term terminates the baroclinic term with a time lag in high Re number,leading to the SBV production.Since the SBV reflects the local stretching enhancement based on the advection-diffusion equation,mixing is influenced by the Sc number in a different behavior if different Re numbers are considered.The time-averaged variable density mixing rate emerges a scaling law with Sc number asχ^(∗)=β·Sc^(−α),where the coefficientβ∼Re−0.2 and the scaling exponentα∼Re−0.385.The understanding of Re number and Sc number effect on variable density mixing provides an opportunity for mixing enhancement from the perspective of designing the viscosity and diffusivity of the fluid mixture.展开更多
Interactions of shock waves and rows of vortices are studied by solving the two-dimensional,compressible Euler equations with a fifth-order weighted essentially non-oscillatory finite difference scheme.For a compressi...Interactions of shock waves and rows of vortices are studied by solving the two-dimensional,compressible Euler equations with a fifth-order weighted essentially non-oscillatory finite difference scheme.For a compressible flow the Mach number of the shock wave and vortex equals to 1.05 and 0.25,respectively.The resulting flow field contains complex shock structures,such as multiple shock focusing and reflecting regions.At the meantime,sound waves are generated,interrupted and reformed when they touch the main and reflected shock waves.展开更多
Manipulating nonlinear excitations,including solitons and vortices,is an essential topic in quantum many-body physics.A new progress in this direction is a protocol proposed in[Phys.Rev.Res.2043256(2020)]to produce da...Manipulating nonlinear excitations,including solitons and vortices,is an essential topic in quantum many-body physics.A new progress in this direction is a protocol proposed in[Phys.Rev.Res.2043256(2020)]to produce dark solitons in a one-dimensional atomic Bose–Einstein condensate(BEC)by quenching inter-atomic interaction.Motivated by this work,we generalize the protocol to a two-dimensional BEC and investigate the generic scenario of its post-quench dynamics.For an isotropic disk trap with a hard-wall boundary,we find that successive inward-moving ring dark solitons(RDSs)can be induced from the edge,and the number of RDSs can be controlled by tuning the ratio of the after-and before-quench interaction strength across different critical values.The role of the quench played on the profiles of the density,phase,and sound velocity is also investigated.Due to the snake instability,the RDSs then become vortex–antivortex pairs with peculiar dynamics managed by the initial density and the after-quench interaction.By tuning the geometry of the box traps,demonstrated as polygonal ones,more subtle dynamics of solitons and vortices are enabled.Our proposed protocol and the discovered rich dynamical effects on nonlinear excitations can be realized in near future cold-atom experiments.展开更多
Even today, with the great progress that has been made in the scientific, technological and computational fields, we are still stunned by the devastating effects brought about by atmospheric phenomena. This paper aims...Even today, with the great progress that has been made in the scientific, technological and computational fields, we are still stunned by the devastating effects brought about by atmospheric phenomena. This paper aims to propose new hypotheses in the field of dynamics to enhance our understanding of the behaviour of atmospheric disturbances caused by rotating winds. I believe that the criteria of classical dynamics that are applied to vortex systems in the atmosphere should be rigorously reviewed. I propose to establish new hypotheses in the field of dynamics, in order to better interpret rotation in nature. These hypotheses have been structured into a new theory that has been tested experimentally by both ourselves and third parties, with positive results. I propose to use the Theory of Dynamic Interactions (TDI) to interpret the behaviour of systems undergoing successive rotations around different axes—which we will refer to as non-coaxial rotations. I hold that this theory applies to air masses and groups of particles in suspension that are accelerated by rotations. Accordingly, it should be used to interpret the behaviour of tornadoes, cyclones and hurricanes. I believe that this proposal could enhance our understanding of these atmospheric phenomena and improve predictions about them.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.12472279)。
文摘In bio-inspired flapping-wing flight,lift generation and flexible deformation are intrinsically coupled.Thus,an experimental study is conducted to reveal the fluid-structure interaction mechanism for a flexible plate undergoing pitching and plunging motion,and theoretical models are proposed to predict either lift or deformation based on a series of simplifying assumptions.It is indicated that flexible plates can effectively reduce the amplitude of lift and pitching moment coefficients during dynamic stall,with increased flexibility leading to higher load reduction.To investigate the effect of plate flexibility on lift,a definition of effective angle of attack is proposed,incorporating the pitching and plunging motion and chordwise deformation of the flexible plate,which can reduce the hysteresis effect of the lift coefficient during dynamic stall.As a consequence,a theoretical model is developed to predict lift based on observed motion and deformation.On the other hand,another theoretical model is developed to predict flexible-plate deformation utilizing aerodynamic forces,revealing the effect of leading-edge vortex evolution on passive deformation.The influence of kinematic parameters,including the maximum effective angle of attack,reduced frequency,and Strouhal number,on the aerodynamic forces is further studied.Compared with the rigid plate,flexible plates exhibit lower sensitivity of aerodynamic forces to changes in kinematic parameters due to their inherent compliance and resulting deformation.The proposed theoretical models can serve as a reference for aerodynamic and deformation prediction in bio-inspired flexible structures.
基金supported by the National Natural Science Foun-dation of China(12274299 and 92050202)the Shanghai Science and Technology Committee(22QA1406600).
文摘While spin-orbit interaction has been extensively studied,few investigations have reported on the interaction between orbital angular momenta(OAMs).In this work,we study a new type of orbit-orbit coupling between the longitudinal OAM and the transverse OAM carried by a three-dimensional(3D)spatiotemporal optical vortex(STOV)in the process of tight focusing.The 3D STOV possesses orthogonal OAMs in the x-y,t-x,and y-t planes,and is preconditioned to overcome the spatiotemporal astigmatism effect.x,y,and t are the axes in the spatiotemporal domain.The corresponding focused wavepacket is calculated by employing the Debye diffraction theory,showing that a phase singularity ring is generated by the interactions among the transverse and longitudinal vortices in the highly confined STOV.The Fourier-transform decomposition of the Debye integral is employed to analyze the mechanism of the orbit-orbit interaction.This is the first revelation of coupling between the longitudinal OAM and the transverse OAM,paving the way for potential applications in optical trapping,laser machining,nonlinear light-matter interactions,and more.
基金National Defense Basic Research Program (B222006060)
文摘The vortex interference mechanism on low Reynolds number between the canard and main wing of the canard-forward sweep wing （Canard-FSW） configurations is simulated numerically by employing the numerical wind tunnel method. The variations of aerodynamic characteristics of Canard-FSW configurations with different positions of the canard are investigated,finding that the aerodynamic interference and mutual coupling effect between the canard and main wing have made great contributions to the lift and stability characteristics of the whole aircraft. Canard can radically improve the surface flow pattern of the main wing. And its own vortex can have a favorable interference on the main wing and can effectively control the airflow boundary layer separation. At small angles of attack,the aerodynamic characteristics are sensitive to the positions of the canard and the main wing,but at high angles of attack,the aerodynamic performances of the configuration are not only related to the shape of the canard （forward or backward）,but also with the size of control force as well as the features of the vortices generated above the main wing and the canard. The different configurations and vortices are illustrated using the velocity vector,streamlines and pressure contours.
文摘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.
基金supported by the National Natural Science Foundation of China (10472047)the Open Fund of State Key Laboratory of Explosion Science Technology, Beijing University of Science and Technology (KFJJ06-3)
文摘Observations are presented from experiments and calculations where a laminar spherical CH4/air flame is perturbed successively by incident and reflected shock waves. The experiments are performed in a standard shock tube arrangement, in which a high-speed shadowgraph imaging system is used to record evolutions of the flame. Numerical simulations are conducted by using second-order wave propagation algorithms, based on two-dimensional axisymmetric Navier-Stokes equations with detailed chemical reactions. Qualitative agreements are obtained between the experimental and numerical results. Under actions of incident shock waves, Richtmyer-Meshkov instability responsible for the flame deformation is induced in the flame, and the distoned flame takes a barrel shape. Then, under subsequent actions of the shock wave reflected from a planar wall, the flame takes an inclined non-symmetrical kidney shape in a symmetric cross section, which means a mushroom-like shape of the flame comes finally into being. The vorticity direction in the ring cap has been altered by the reflected shock＇s action, which makes the head of the mushroom-like flame extend quickly to the side wall.
基金Supported by the National Natural Science Foundation of China (No. 50236030, No. 50076038) and the Major State Basic Research Development Program of China (No. G19990222).
文摘In order to understand the interaction between large-scale vortex structure and particles, a two-way coupling temporal mixing layer laden with particles at a Stokes number of 5 with different mass loading planted initially in the upper half region is numerically studied. The pseudospectral method is used for the flow fluid and the Lagrangian approach is employed to trace particles. The momentum coupling effect introduced by a particle is approximated to a point force. The simulation results show that the coherent structures are still dominant in the mixing layer, but the large-scale vortex structure and particle dispersion are modulated. The length of large-scale vortex structure is shortened and the pairing is delayed. At the same time, the particles are distributed more evenly in the whole flow field as the mass loading is increased, but the particle dispersion along the transverse direction differs from that along the spanwise direction, which indicates that the effect by the addition of particle on the spanwise large-scale vortex structure is different from the streamwise counterpart.
基金the National Natural Science Foundation of China(Nos.12032012,12072156)the National Key Laboratory Foundation of China(No.61422202103)the Priority Academic Program Development of Jiangsu Higher Education Institutions,China(PAPD).
文摘Numerical simulation methods for unsteady vortex field of helicopter rotor with high resolution and low dissipation TENO8-AA primitive variables reconstruction schemes are established based on moving-embedded grid and Navier-Stokes equations.Firstly,the Targeted Essentially Non-Oscillatory(TENO)scheme are developed by employing ENO-like candidate stencil selection strategy,and the candidate stencil is adopted with optimal weight in smooth region while it is discarded completely in discontinuous region,which reduces the dissipation and dispersion errors and approaches better spectral properties.Then,the aerodynamic characteristics of Helishape-7A model rotor in Blade Vortex Interaction(BVI)state and the flowfield of Lynx rotor in hover are simulated,which validates that the blade tip vortex trajectory with larger wake age and more details of vortex can be captured by TENO8-AA scheme with only a quarter of grid points and half time comparing to WENO-JS scheme.Moreover,the simulation accuracy of thrust coefficient is improved by up to 36%.Finally,the analyses for BVI and aeroacoustic characteristics of Operational Loads Survey(OLS)rotor are conducted,and the different forms of interaction mechanism are explored,such as oblique and parallel interactions.The results indicate that TENO scheme not only ensures the resolution of simulation in discontinuous region,but also minimizes the numerical dissipation in smooth region dominated by blade tip vortex.Therefore,the acoustic pressure peak prediction error of rotor in BVI state is significantly reduced to 5.6%and 0.8%at two microphone locations,respectively.
基金This work was supported by the National Natural Science Foundation of China(Grant No.51976082)Qing Lan Project.
文摘The effects of inlet gas parameters and sloping sidewall angle on the flame structure and combustion limit with and without sidewall were experimentally investigated.Flame height and impact angle were obtained by chemiluminescence intensity analysis of CH*distribution.First,the combustion characteristics of flame with and without sidewall at different equivalence ratios were explored;then,the influence of Reynolds number and inlet gas temperature on flame structure and combustion limit of v-shaped flame with sidewall were analyzed,and the results with sidewall were compared with those without sidewall.Finally,the variation trend of flame parameters with different sloping sidewall angles was analyzed.The experimental results show that the existence of sidewall makes flame shape change from“M-shaped”to“inverted N-shaped”,and conical shape to trapezoidal shape.The inhibition effect of sidewall on flame stretching downstream is strengthened with the increase in Reynolds number;but as the temperature of the inlet gas increases,the inhibitory effect is obviously weakened.When sloping sidewall angle decreases from 90°to 55°at 5°intervals,flame height and impact angle of v-shaped flame reach the extreme value whenβ=80°.Compared with the case without sidewall,the range of v-shaped flame with sidewall has no obvious trend of broadening or shrinking when inlet gas temperature is increased;however,as sloping sidewall angle decreases,the range of the v-shaped flame shrinks obviously and flammability limit increases significantly.
基金National Basic Research Program of China(2011CB403505,2010CB950402)National Natural Science Foundation of China(40975052,41175059)
文摘Based on 6-hourly sensible heat flux and latent heat flux from the NCEP Climate Forecast System Reanalysis(CFSR) and circulation data from the Japanese 25-year Reanalysis(JRA-25),the initial developing process of tropical cyclone Mindulle(1005) in 2010 has been diagnosed to reveal the impact of air-sea interaction over the South China Sea(SCS) on the genesis of its incipient vortex.The results show that the incipient vortex first occurred east of the Luzon Island on 0000 UTC 20 August,suggesting that the topographic forcing of the Luzon Island for easterly winds over the western Pacific might be one of the factors responsible for the formation of the incipient vortex.During the formation stage of the incipient vortex,strong southeasterlies over the SCS caused warm water of the middle and eastern SCS to flow toward the Luzon Island due to Ekman transport resulting from wind stress,leading to an increase of the sea surface temperature and sensible heat flux into the atmosphere.Although the anomalous sensible heating favored surface pressure to reduce,it was not conducive to the increase of local vorticity associated with the vortex above the heating area because,according to the atmospheric thermal adaptation theory,the anticyclonic vorticity would be created in the lower troposphere due to the decreased vertical gradient of the sensible heating.However,the ascending motions occurred over the eastern area of the anomalous sensible heating due to the augmentation of the vorticity advection with increasing height,causing water vapor to condense in the middle and upper troposphere.In turn,cyclonic vorticity was generated in the lower troposphere due to the increased vertical gradient of the condensation latent heating,resulting in the formation and further growth of the incipient vortex.Therefore,the vorticity creation due to the condensation heating played a dominant role during the subsequent enhancing stage of the incipient vortex.
基金The project supported by National Natural Science Foundation of China under Grant No. 10574012
文摘When a Bose-Einstein condensate is set to rotate, superfluid vortices will be formed, which finally condense into a vortex lattice as the rotation frequency further increases. We show that the dipole-dipole interactions renormalize the short-range interaction strength and result in a distinction between interactions of parallel-polarized atoms and interactions of antiparallel-polarized atoms. This effect may lead to a spontaneous breakdown of the rapidly rotating Bose condensate into a novel anti-ferromagnetic-like vortex lattice. The upward-polarized Bose condensate forms a vortex lattice, which is staggered against a downward-polarized vortex lattice. A phase diagram related to the coupling strength is obtained.
基金Project supported by the National Natural Science Foundation of China(Nos.91741113,91841303,and 91941301)。
文摘The underlying effect of vortex interaction characterized by the merging and non-merging on mixing enhancement is of fundamental significance to understand the flow dynamics of strut injectors in scramjets.Starting from a simplified configuration of a vortex generator,this study focuses on the influence of geometric parameters on vortex structures and fluid mixing through compressible Navier-Stokes(NS)simulations.By adjusting the induction of outer vortices,the inner co-rotating vortex pair exhibits two modes of interaction(merging/separation regime)reflected by closer/farther vortex centers.Defined by the zero variation rate of the inner vortex spacing,the critical state of equilibrium is determined.The critical condition is well predicted by a theoretical model based on the Biot-Savart law.Through the introduction of mixedness and mixing time,the intrinsic impact of interaction modes on fluid mixing is revealed.Compared with the vortex dynamics in the merging regime,the one in the separation regime is more effective for passive scalar mixing augmentation.With efficient material stretching characterized by the higher interface stretching factor and averaging finite-time Lyapunov exponent(FTLE),the mixing time is shortened by as much as 2.5 times in the separation regime.The implication of the present two interaction regimes in mixing enhancement physically reflected by the averaging FTLE has the potential to improve the combustion performance and shorten the combustor chamber.
基金Project supported by the National Natural Science Foundation of China(Grant No.11204026)the Fundamental Research Funds for Central Universities of the Ministry of Education of China(Grant No.n130405011)
文摘The dynamic properties of interacting vortex-antivortex pairs in thin film are studied by analytical calculations. An- alytical expressions for the magnetization vector distribution of vortex-antivortex pairs and the trivortex states are given. The magnetic states of the vortices are treated as having rigid structures, i.e., the vortex maintains its spin distribution when moving. The trajectories of the vortex cores are calculated by the Thiele＇s equation. It is found that the vortex-antivortex pair rotates around each other when they have opposite polarities, however, vortex and antivortex cores move along straight lines when they have the same polarity. The frequency of the rotation decreases with increasing the distance between the two cores of vortex-antivortex pair, and it has a lower value when a third vortex is introduced.
文摘The influences due to several AVIs （airfoil-vortex interactions） are studied by using a two-dimensional CFD （computational fluid dynamics） method. The primary goal of this effort is to assess the variation of vortex center location and vortex circulation associated with sequential AVI toward an improvement of the hybrid method of CFD and prescribed wake model, which closely relates to predicting the BVI （blade-vortex interaction） noise radiated from a helicopter rotor. The representative of sequential AVI is performed by single vortex and two airfoils. Investigations with respect to vortex center location and vortex circulation after AVIs have been made by varying the miss-distance, which is the vertical distance between the airfoil leading edge and the vortex center. Correlations between miss-distance and vorticity field show that there exists complicated vortex wake flow with several vortices newly induced in 1st AVI. The pressure fluctuation amplitude clarifies that the intensity in 2nd AV1 is significantly low compared to the intensity in 1st AVI due to the influence of vortex dissipation. Simulations turned out to modify the vortex center location represented by the hybrid method using an offset value for a streamwise direction and to dissipate the vortex circulation for improving the accuracy of BVI noise prediction.
基金supported by the National Natural Science Foundation of China (No. 51178397)Technological Research and Development Programs of the Ministry of Railways (No. 2010G004-L)
文摘Large eddy simulation is performed to study three-dimensional wave-current interaction with a square cylinder at different Reynolds numbers, ranging from 1,000 to 600,000. The Keulegan-Carpenter number is relevantly a constant of 0.6 for all cases. The Strouhal number, the mean and the RMS values of the effective drag coefficient in the streamwise and transverse directions are computed for various Reynolds numbers, and the velocity of a rep- resentative point in the turbulent zone is simulated to find the turbulent feature. It is found that the wave-current interaction should be considered as three-dimensional flow when the Reynolds number is high; under wave-current effect, there exists a critical Reynolds number, and when the Reynolds number is smaller than the critical one, current effect on wave can be nearly neglected; conversely, with the Reynolds number increasing, wave-currentstructure interaction is sensitive to the Reynolds number.
文摘This article introduces a numerical scheme on the basis of semi-implicit method for pressure-linked equations （SIMPLE） algorithm to simulate incompressible unsteady flows with fluid-structure interaction. The Navier-Stokes equation is discretized spatially with collocated finite volume method and Eulerian implicit method in time domain. The hybrid method that combines immersed boundary method （IBM） and volume of fluid （VOF） method is used to deal with rigid body motion in fluid domain. The details of movement of immersed boundary （IB） and calculation of VOF are also described. This method can be easily applied to any existing finite-volume-based computational fluid dynamics （CFD） solver without complex operation, with which fluid flow interaction of arbitrarily complex geometry can be realized on a fixed mesh. The method is verified by low Reynolds number flows passing both stationary and oscillating cylinders. The drag and lift coefficients acquired by the study well accord with other published results, which indicate the reasonability of the proposed method.
基金Project supported by the National Natural Science Foundation of China(Grant No.10874125)
文摘This paper studies in detail the interaction of two edge dislocations nested in a Gaussian beam propagating in free space. It shows that in free-space propagation the edge dislocations are unstable and vanish, and two noncanonical vortices with opposite topological charge take place when off-axis distances cl and c2 of two edge dislocations are nonzero, and the condition k2w08＋ 32c1c2（w02- 2C1C2）Z2 〉 0 is fulfilled （k-wave number, w0-waist width）. A noncanonical vortex appears when one off-axis distance is zero. However, one edge dislocation is stable when two edge dislocations are perpendicular and one off-axis distance is zero. Two perpendicular edge dislocations both with zero off-axis distance are also stable. The analytical results are illustrated by numerical examples.
基金This work was supported by the National Natural Science Foundation of China(NSFC)(Grant No.91941301)the Key Research and Development Project of Sichuan Province(Grant No.2019ZYZF0002)。
文摘Effects of Reynolds(Re)number and Schmidt(Sc)number on the flow structures and variable density mixing are numerically investigated through the canonical shock cylindrical bubble interaction.By determining the viscosity and diffusivity within a wide range,the controlling parameters,total vortex circulation,and compression rate,are conservative under a broad range of Re and Sc numbers(Re≈10^(3)-10^(5)and Sc≈0.1-5)in the same shock Mach(Ma)number condition(Ma=2.4).As for the Re number effect,the circulation of secondary baroclinic vorticity(SBV),induced by the main vortex centripetal acceleration,is observed to be higher in high Re number and vice versa.Based on the vorticity transport equation decomposition,a growth-inhibition vorticity dynamics balance mechanism is revealed:the vorticity viscous term grows synchronously with baroclinic production to inhibit SBV production in low Re number.By contrast,the viscous term terminates the baroclinic term with a time lag in high Re number,leading to the SBV production.Since the SBV reflects the local stretching enhancement based on the advection-diffusion equation,mixing is influenced by the Sc number in a different behavior if different Re numbers are considered.The time-averaged variable density mixing rate emerges a scaling law with Sc number asχ^(∗)=β·Sc^(−α),where the coefficientβ∼Re−0.2 and the scaling exponentα∼Re−0.385.The understanding of Re number and Sc number effect on variable density mixing provides an opportunity for mixing enhancement from the perspective of designing the viscosity and diffusivity of the fluid mixture.
基金Supported by the National Natural Science Foundation of China(11072053)
文摘Interactions of shock waves and rows of vortices are studied by solving the two-dimensional,compressible Euler equations with a fifth-order weighted essentially non-oscillatory finite difference scheme.For a compressible flow the Mach number of the shock wave and vortex equals to 1.05 and 0.25,respectively.The resulting flow field contains complex shock structures,such as multiple shock focusing and reflecting regions.At the meantime,sound waves are generated,interrupted and reformed when they touch the main and reflected shock waves.
基金Project supported by the Natural Science Foundation of Zhejiang Province of China(Grant Nos.LQ22A040006,LY21A040004,LR22A040001,and LZ21A040001)the National Natural Science Foundation of China(Grant Nos.11835011 and 12074342).
文摘Manipulating nonlinear excitations,including solitons and vortices,is an essential topic in quantum many-body physics.A new progress in this direction is a protocol proposed in[Phys.Rev.Res.2043256(2020)]to produce dark solitons in a one-dimensional atomic Bose–Einstein condensate(BEC)by quenching inter-atomic interaction.Motivated by this work,we generalize the protocol to a two-dimensional BEC and investigate the generic scenario of its post-quench dynamics.For an isotropic disk trap with a hard-wall boundary,we find that successive inward-moving ring dark solitons(RDSs)can be induced from the edge,and the number of RDSs can be controlled by tuning the ratio of the after-and before-quench interaction strength across different critical values.The role of the quench played on the profiles of the density,phase,and sound velocity is also investigated.Due to the snake instability,the RDSs then become vortex–antivortex pairs with peculiar dynamics managed by the initial density and the after-quench interaction.By tuning the geometry of the box traps,demonstrated as polygonal ones,more subtle dynamics of solitons and vortices are enabled.Our proposed protocol and the discovered rich dynamical effects on nonlinear excitations can be realized in near future cold-atom experiments.
文摘Even today, with the great progress that has been made in the scientific, technological and computational fields, we are still stunned by the devastating effects brought about by atmospheric phenomena. This paper aims to propose new hypotheses in the field of dynamics to enhance our understanding of the behaviour of atmospheric disturbances caused by rotating winds. I believe that the criteria of classical dynamics that are applied to vortex systems in the atmosphere should be rigorously reviewed. I propose to establish new hypotheses in the field of dynamics, in order to better interpret rotation in nature. These hypotheses have been structured into a new theory that has been tested experimentally by both ourselves and third parties, with positive results. I propose to use the Theory of Dynamic Interactions (TDI) to interpret the behaviour of systems undergoing successive rotations around different axes—which we will refer to as non-coaxial rotations. I hold that this theory applies to air masses and groups of particles in suspension that are accelerated by rotations. Accordingly, it should be used to interpret the behaviour of tornadoes, cyclones and hurricanes. I believe that this proposal could enhance our understanding of these atmospheric phenomena and improve predictions about them.
