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.展开更多
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 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.展开更多
The average-passage equation system (APES) provides a rigorous mathematical framework for account- ing for the unsteady blade row interaction through multistage compressors in steady state environment by introducing...The average-passage equation system (APES) provides a rigorous mathematical framework for account- ing for the unsteady blade row interaction through multistage compressors in steady state environment by introducing de- terministic correlations (DC) that need to be modeled to close the equation system. The primary purpose of this study was to provide insight into the DC characteristics and the in- fluence of DC on the time-averaged flow field of the APES. In Part 2 of this two-part paper, the influence of DC on the time-averaged flow field was systematically studied; Several time-averaging computations boundary conditions and DC were conducted with various for the downstream stator in a transonic compressor stage, by employing the CFD solver developed in Part 1 of this two-part paper. These results were compared with the time-averaged unsteady flow field and the steady one. The study indicat;d that the circumferential- averaged DC can take into account major part of the unsteady effects on spanwise redistribution of flow fields in compres- sors. Furthermore, it demonstrated that both deterministic stresses and deterministic enthalpy fluxes are necessary to reproduce the time-averaged flow field.展开更多
The average-passage equation system (APES) provides a rigorous mathematical framework for account- ing for the unsteady blade row interaction through multi- stage compressors in steady state environment by introduc-...The average-passage equation system (APES) provides a rigorous mathematical framework for account- ing for the unsteady blade row interaction through multi- stage compressors in steady state environment by introduc- ing deterministic correlations (DC) that need to be modeled to close the equation system. The primary purpose of this study is to provide insight into the DC characteristics and the influence of DC on the time-averaged flow field of the APES. In Part 1 of this two-part paper, firstly a 3D viscous unsteady and time-averaging flow CFD solver is developed to investi- gate the APES technique. Then steady and unsteady simu- lations are conducted in a transonic compressor stage. The results from both simulations are compared to highlight the significance of the unsteady interactions. Furthermore, the distribution characteristics of DC are studied and the DC at the rotor/stator interface are compared with their spatial cor- relations (SC). Lastly, steady and time-averaging (employing APES with DC) simulations for the downstream stator alone are conducted employing DC derived from the unsteady re- suits. The results from steady and time-averaging simula- tions are compared with the time-averaged unsteady results. The comparisons demonstrate that the simulation employing APES with DC can reproduce the time-averaged field and the 3D viscous time-averaging flow solver is validated.展开更多
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.展开更多
The blades of large-scale wind turbines can obviously deform during operation,and such a deformation can affect the wind turbine’s output power to a certain extent.In order to shed some light on this phenomenon,for w...The blades of large-scale wind turbines can obviously deform during operation,and such a deformation can affect the wind turbine’s output power to a certain extent.In order to shed some light on this phenomenon,for which limited information is available in the literature,a bidirectional fluid-structure interaction(FSI)numerical model is employed in this work.In particular,a 5 MW large-scale wind turbine designed by the National Renewable Energy Laboratory(NREL)of the United States is considered as a testbed.The research results show that blades’deformation can increase the wind turbine’s output power by 135 kW at rated working conditions.Compared with the outcomes of the simulations conducted using the model with no blade deformation,the results obtained with the FSI model are closer to the experimental data.It is concluded that the bidirectional FSI model can replicate the working conditions of wind turbines with great fidelity,thereby providing an effective method for wind turbine design and optimization.展开更多
Rotor-to-rotor interaction among neighboring rotors of a multirotor has great significance for aerodynamically efficient multirotor design. Current research is conducted to analyze aerodynamic performance of different...Rotor-to-rotor interaction among neighboring rotors of a multirotor has great significance for aerodynamically efficient multirotor design. Current research is conducted to analyze aerodynamic performance of different octocopter configurations amid hover and forward flight. Conventional and coaxial configurations are studied and a hybrid configuration is also proposed to rectify the disadvantages associated with the earlier two. Comparison is carried out for the aforementioned configurations along with comparison of coaxial and hybrid octocopters with bigger diameter rotors in the same confined space for high thrust requirement missions. Vertical spacing of coaxial configuration is also studied. Virtual Blade Method (VBM) is considered herein due to its great computational efficiency. The results show that there are 11.89% and 14.22% loss in thrust for coaxial octocopter compared to conventional and hybrid configurations with normal size rotors and 15.61% loss compared to hybrid configuration with bigger rotors in hover, whereas coaxial square configuration performs the worst in forward flight with a lift loss of 9.1%, 14.77% and 18.8% compared to coaxial diamond, conventional and hybrid configurations with normal size rotors and 9.96% and 17.82% loss compared to coaxial diamond and hybrid configurations with bigger rotors. Combined FM shows that hybrid configuration outperforms other octocopter configurations in overall aerodynamic performance.展开更多
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 centrifugal pump is a prevalent power equipment widely used in different engineering patterns,and the impeller blade wrap angle significantly impacts its performance.A numerical investigation was conducted to anal...The centrifugal pump is a prevalent power equipment widely used in different engineering patterns,and the impeller blade wrap angle significantly impacts its performance.A numerical investigation was conducted to analyze the influence of the blade wrap angle on flow characteristics and energy distribution of a centrifugal pump evaluated as a low specific speed with a value of 69.This study investigates six impellermodels that possess varying blade wrap angles(95°,105°,115°,125°,135°,and 145°)that were created while maintaining the same volute and other geometrical characteristics.The investigation of energy loss was conducted to evaluate the values of total and entropy generation rates(TEG,EGR).The fluid-structure interaction was considered numerically using the software tools ANSYS Fluent and ANSYSWorkbench.The elastic structural dynamic equation was used to estimate the structural response,while the shear stress transport k–ωturbulence model was utilized for the fluid domain modeling.The findings suggest that the blade wrap angle has a significant influence on the efficiency of the pump.The impeller featuring a blade wrap angle of 145°exhibits higher efficiency,with a notable increase of 3.76%relative to the original model.Variations in the blade wrap angle impact the energy loss,shaft power,and pump head.The model with a 145°angle exhibited a maximum equivalent stress of 14.8MPa and a total deformation of 0.084 mm.The results provide valuable insights into the intricate flow mechanism of the centrifugal pump,particularly when considering various blade wrap angles.展开更多
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.展开更多
The effects of blade lean and vortex design on the aerodynamics of a turbine entry nozzle guide vane (NGV) are considered using computational fluid dynamics. The aim of the work is to address some of the uncertainti...The effects of blade lean and vortex design on the aerodynamics of a turbine entry nozzle guide vane (NGV) are considered using computational fluid dynamics. The aim of the work is to address some of the uncertainties which have arisen from previous studies where conflicting results have been reported for the effect on the NGV. The configuration was initially based on the energy efficient engine turbine which also served as the validation case for the computational method. A total of 17 NGV configurations were evaluated to study the effects of lean and vortex design on row efficiency and secondary kinetic energy. The distribution of mass flow ratio is introduced as an additional factor in the assessment of blade lean effects. The results show that in the turbine entry NGV, the secondary flow strength is not a dominant factor that determines NGV losses and therefore the changes of loading distribution due to blade lean and the associated loss mecha- nisms should be regarded as a key factor. Radial mass flow redistribution under different NGV lean and twist is demonstrated as an addition key factor influencin~ row efficiency.展开更多
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.展开更多
This study concerns the characterization of both the steady and unsteady flows and the analysis of stator/rotor interactions of a two-stage axial turbine. The predicted aerodynamic performances show noticeable differe...This study concerns the characterization of both the steady and unsteady flows and the analysis of stator/rotor interactions of a two-stage axial turbine. The predicted aerodynamic performances show noticeable differences when simulating the turbine stages simultaneously or separately. By considering the multi-blade per row and the scaling technique, the Computational fluid dynamics(CFD) produced better results concerning the effect of pitchwise positions between vanes and blades. The recorded pressure fluctuations exhibit a high unsteadiness characterized by a space–time periodicity described by a double Fourier decomposition. The Fast Fourier Transform FFT analysis of the static pressure fluctuations recorded at different interfaces reveals the existence of principal harmonics and their multiples, and each lobed structure of pressure wave corresponds to the number of vane/blade count. The potential effect is seen to propagate both upstream and downstream of each blade row and becomes accentuated at low mass flow rates. Between vanes and blades, the potential effect is seen to dominate the quasi totality of blade span, while downstream the blades this effect seems to dominate from hub to mid span. Near the shroud the prevailing effect is rather linked to the blade tip flow structure.展开更多
基金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.
基金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.
基金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 (51006006,51136003,50976010,50976009)the National Basic Research Program of China (2012CB72 0205)+2 种基金the Aeronautical Science Foundation of China (2010ZB51)the 111 Project (B08009)the National Science Special Foundation for Post-doctoral Scientists of China (201104049)
文摘The average-passage equation system (APES) provides a rigorous mathematical framework for account- ing for the unsteady blade row interaction through multistage compressors in steady state environment by introducing de- terministic correlations (DC) that need to be modeled to close the equation system. The primary purpose of this study was to provide insight into the DC characteristics and the in- fluence of DC on the time-averaged flow field of the APES. In Part 2 of this two-part paper, the influence of DC on the time-averaged flow field was systematically studied; Several time-averaging computations boundary conditions and DC were conducted with various for the downstream stator in a transonic compressor stage, by employing the CFD solver developed in Part 1 of this two-part paper. These results were compared with the time-averaged unsteady flow field and the steady one. The study indicat;d that the circumferential- averaged DC can take into account major part of the unsteady effects on spanwise redistribution of flow fields in compres- sors. Furthermore, it demonstrated that both deterministic stresses and deterministic enthalpy fluxes are necessary to reproduce the time-averaged flow field.
基金supported by the National Natural Science Foundation of China (51006006,51136003,50976010,50976009)the National Basic Research Program of China (2012CB720205)+2 种基金the Aeronautical Science Foundation of China (2010ZB51)the 111 Project (B08009)the National Science Special Foundation for Post-doctoral Scientists of China (201104049)
文摘The average-passage equation system (APES) provides a rigorous mathematical framework for account- ing for the unsteady blade row interaction through multi- stage compressors in steady state environment by introduc- ing deterministic correlations (DC) that need to be modeled to close the equation system. The primary purpose of this study is to provide insight into the DC characteristics and the influence of DC on the time-averaged flow field of the APES. In Part 1 of this two-part paper, firstly a 3D viscous unsteady and time-averaging flow CFD solver is developed to investi- gate the APES technique. Then steady and unsteady simu- lations are conducted in a transonic compressor stage. The results from both simulations are compared to highlight the significance of the unsteady interactions. Furthermore, the distribution characteristics of DC are studied and the DC at the rotor/stator interface are compared with their spatial cor- relations (SC). Lastly, steady and time-averaging (employing APES with DC) simulations for the downstream stator alone are conducted employing DC derived from the unsteady re- suits. The results from steady and time-averaging simula- tions are compared with the time-averaged unsteady results. The comparisons demonstrate that the simulation employing APES with DC can reproduce the time-averaged field and the 3D viscous time-averaging flow solver is validated.
基金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.
基金supported by the CHN Energy United Power Technology Co.,Ltd.,China(Contract No.2020-75).
文摘The blades of large-scale wind turbines can obviously deform during operation,and such a deformation can affect the wind turbine’s output power to a certain extent.In order to shed some light on this phenomenon,for which limited information is available in the literature,a bidirectional fluid-structure interaction(FSI)numerical model is employed in this work.In particular,a 5 MW large-scale wind turbine designed by the National Renewable Energy Laboratory(NREL)of the United States is considered as a testbed.The research results show that blades’deformation can increase the wind turbine’s output power by 135 kW at rated working conditions.Compared with the outcomes of the simulations conducted using the model with no blade deformation,the results obtained with the FSI model are closer to the experimental data.It is concluded that the bidirectional FSI model can replicate the working conditions of wind turbines with great fidelity,thereby providing an effective method for wind turbine design and optimization.
基金supported by the National Natural Science Foundation of China(No.11972190).
文摘Rotor-to-rotor interaction among neighboring rotors of a multirotor has great significance for aerodynamically efficient multirotor design. Current research is conducted to analyze aerodynamic performance of different octocopter configurations amid hover and forward flight. Conventional and coaxial configurations are studied and a hybrid configuration is also proposed to rectify the disadvantages associated with the earlier two. Comparison is carried out for the aforementioned configurations along with comparison of coaxial and hybrid octocopters with bigger diameter rotors in the same confined space for high thrust requirement missions. Vertical spacing of coaxial configuration is also studied. Virtual Blade Method (VBM) is considered herein due to its great computational efficiency. The results show that there are 11.89% and 14.22% loss in thrust for coaxial octocopter compared to conventional and hybrid configurations with normal size rotors and 15.61% loss compared to hybrid configuration with bigger rotors in hover, whereas coaxial square configuration performs the worst in forward flight with a lift loss of 9.1%, 14.77% and 18.8% compared to coaxial diamond, conventional and hybrid configurations with normal size rotors and 9.96% and 17.82% loss compared to coaxial diamond and hybrid configurations with bigger rotors. Combined FM shows that hybrid configuration outperforms other octocopter configurations in overall aerodynamic performance.
基金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.
文摘The centrifugal pump is a prevalent power equipment widely used in different engineering patterns,and the impeller blade wrap angle significantly impacts its performance.A numerical investigation was conducted to analyze the influence of the blade wrap angle on flow characteristics and energy distribution of a centrifugal pump evaluated as a low specific speed with a value of 69.This study investigates six impellermodels that possess varying blade wrap angles(95°,105°,115°,125°,135°,and 145°)that were created while maintaining the same volute and other geometrical characteristics.The investigation of energy loss was conducted to evaluate the values of total and entropy generation rates(TEG,EGR).The fluid-structure interaction was considered numerically using the software tools ANSYS Fluent and ANSYSWorkbench.The elastic structural dynamic equation was used to estimate the structural response,while the shear stress transport k–ωturbulence model was utilized for the fluid domain modeling.The findings suggest that the blade wrap angle has a significant influence on the efficiency of the pump.The impeller featuring a blade wrap angle of 145°exhibits higher efficiency,with a notable increase of 3.76%relative to the original model.Variations in the blade wrap angle impact the energy loss,shaft power,and pump head.The model with a 145°angle exhibited a maximum equivalent stress of 14.8MPa and a total deformation of 0.084 mm.The results provide valuable insights into the intricate flow mechanism of the centrifugal pump,particularly when considering various blade wrap angles.
基金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.
基金the support of Cranfield University and China Aviation Powerplant Research Institute of AVIC
文摘The effects of blade lean and vortex design on the aerodynamics of a turbine entry nozzle guide vane (NGV) are considered using computational fluid dynamics. The aim of the work is to address some of the uncertainties which have arisen from previous studies where conflicting results have been reported for the effect on the NGV. The configuration was initially based on the energy efficient engine turbine which also served as the validation case for the computational method. A total of 17 NGV configurations were evaluated to study the effects of lean and vortex design on row efficiency and secondary kinetic energy. The distribution of mass flow ratio is introduced as an additional factor in the assessment of blade lean effects. The results show that in the turbine entry NGV, the secondary flow strength is not a dominant factor that determines NGV losses and therefore the changes of loading distribution due to blade lean and the associated loss mecha- nisms should be regarded as a key factor. Radial mass flow redistribution under different NGV lean and twist is demonstrated as an addition key factor influencin~ row efficiency.
基金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 study concerns the characterization of both the steady and unsteady flows and the analysis of stator/rotor interactions of a two-stage axial turbine. The predicted aerodynamic performances show noticeable differences when simulating the turbine stages simultaneously or separately. By considering the multi-blade per row and the scaling technique, the Computational fluid dynamics(CFD) produced better results concerning the effect of pitchwise positions between vanes and blades. The recorded pressure fluctuations exhibit a high unsteadiness characterized by a space–time periodicity described by a double Fourier decomposition. The Fast Fourier Transform FFT analysis of the static pressure fluctuations recorded at different interfaces reveals the existence of principal harmonics and their multiples, and each lobed structure of pressure wave corresponds to the number of vane/blade count. The potential effect is seen to propagate both upstream and downstream of each blade row and becomes accentuated at low mass flow rates. Between vanes and blades, the potential effect is seen to dominate the quasi totality of blade span, while downstream the blades this effect seems to dominate from hub to mid span. Near the shroud the prevailing effect is rather linked to the blade tip flow structure.