In a Mach 3.8 wind tunnel, both instantaneous and time-averaged flow structures of different scales around a blunt double-cone with or without supersonic film cooling were visualized via nano-tracer planar laser scatt...In a Mach 3.8 wind tunnel, both instantaneous and time-averaged flow structures of different scales around a blunt double-cone with or without supersonic film cooling were visualized via nano-tracer planar laser scattering (NPLS), which has a high spatiotemporal resolution. Three experimental cases with different injection mass flux rates were carried out. Many typical flow structures were clearly shown, such as shock waves, expansion fans, shear layers, mixing layers, and turbulent boundary layers. The analysis of two NPLS images with an interval of 5 us revealed the temporal evolution characteristics of flow structures. With matched pressures, the laminar length of the mixing layer was longer than that in the case with a larger mass flux rate, but the full covered region was shorter. Structures like K-H (Kelvin-Helmholtz) vortices were clearly seen in both flows. Without injection, the flow was similar to the supersonic flow over a backward- facing step, and the structures were relatively simpler, and there was a longer laminar region. Large scale structures such as hairpin vortices were visualized. In addition, the results were compared in part with the schlieren images captured by others under similar conditions.展开更多
The objective of this paper is to address the transient flow structures around a pitching hydrofoil by com- bining physical and numerical studies. In order to predict the dynamic behavior of the flow structure effecti...The objective of this paper is to address the transient flow structures around a pitching hydrofoil by com- bining physical and numerical studies. In order to predict the dynamic behavior of the flow structure effectively, the Lagrangian coherent structures (LCS) defined by the ridges of the finite-time Lyapunov exponent (FTLE) are utilized under the framework of Navier-Stokes flow computations. In the numerical simulations, the k-w shear stress trans- port (SST) turbulence model, coupled with a two-equation F-Reo transition model, is used for the turbulence closure. Results are presented for a NACA66 hydrofoil undergoing slowly and rapidly pitching motions from 0° to 15° then back to 0° at a moderate Reynolds number Re = 7.5 × 105. The results reveal that the transient flow structures can be observed by the LCS method. For the slowly pitching case, it consists of five stages: quasi-steady and laminar, transition from laminar to turbulent, vortex development, large-scale vortex shedding, and reverting to laminar. The observation of LCS and Lagrangian particle tracers elucidates that the trailing edge vortex is nearly attached and stable during the vortex development stage and the interaction between the leading and trailing edge vortex caused by the adverse pres- sure gradient forces the vortexes to shed downstream during the large-scale vortex shedding stage, which corresponds to obvious fluctuations of the hydrodynamic response. For the rapidly pitching case, the inflection is hardly to be observed and the stall is delayed. The vortex formation, interaction, and shedding occurred once instead of being repeated three times, which is responsible for just one fluctuation in the hydrody- namic characteristics. The numerical results also show that the FTLE field has the potential to identify the transient flows, and the LCS can represent the divergence extent of infinite neighboring particles and capture the interface of the vortex region.展开更多
A bionic experimental platform was designed for the purpose of investigating time accurate three-dimensional flow field, using digital particle image velocimetry (DSPIV). The wake behind the flapping trail of a robo...A bionic experimental platform was designed for the purpose of investigating time accurate three-dimensional flow field, using digital particle image velocimetry (DSPIV). The wake behind the flapping trail of a robotic fish model was studied at high spatial resolution. The study was performed in a water channel. A robot fish model was designed and built. The model was fixed onto a rigid support frame- work using a cable-supporting method, with twelve stretched wires. The entire tail of the model can perform prescribed motions in two degrees of freedom, mainly in carangiform mode, by driving its afterbody and lunate caudal fin respectively. The DSPIV system was set up to operate in a trans- lational manner, measuring velocity field in a series of parallel slices. Phase locked measurements were repeated for a number of runs, allowing reconstruction of phase average flow field. Vortex structures with phase history of the wake were obtained. The study reveals some new and complex three-dimensional flow structures in the wake of the fish, including "reverse hairpin vortex" and "reverse Karman S-H vortex rings", allowing insight into physics of this complex flow.展开更多
In this work,the dielectric-barrier-discharge plasma actuator was employed to study the flow structures induced by the plasma actuator over a flat plate and a wall-mounted hump.A phenomenological dielectric-barrier-di...In this work,the dielectric-barrier-discharge plasma actuator was employed to study the flow structures induced by the plasma actuator over a flat plate and a wall-mounted hump.A phenomenological dielectric-barrier-discharge plasma model which regarded the plasma effect as the body force was implemented into the Navier–Stokes equations solved by the method of large eddy simulations.The results show that a series of vortex pairs,which indicated dipole formation and periodicity distribution were generated in the boundary layer when the plasma was applied to the flow over a flat plane.They would enhance the energy exchanged between the near wall region and the free stream.Besides,their spatial trajectories are deeply affected by the actuation strength.When the actuator was engaged in the flow over a wall-mounted hump,the vortex pairs were also produced,which was able to delay flow separation as well as to promote flow reattachment and reduce the generation of a vortex,achieving the goal of reducing dissipation and decreasing flow resistance.展开更多
The evolution of wake structures and variation of the forces on a flat plate in harmonic oscillatory and in-line combined flows are obtained numerically by improved discrete vortex method. For the oscillatory oncoming...The evolution of wake structures and variation of the forces on a flat plate in harmonic oscillatory and in-line combined flows are obtained numerically by improved discrete vortex method. For the oscillatory oncoming flow cases, when K_c number varies from 2 to 40, the vortex pattern changes from a 'harmonic wave' shaped (in a range of small K_c numbers) to a slight inclined 'harmonic wave' shaped (in a range of moderate K_c numbers), then to inclined vortex clusters with an angle of 50 ° to the oncoming flow direction (at K_c = 20), at last, as K_c number becomes large, the vortex pattern is like a normal Karman vortex street. The well predicted drag and inertia force coefficients are obtained, which are more close to the results of Keulegan & Carpenter's experiment as compared with previous vortex simulation by other au- thors. The existence of minimum point of inertia force coefficient C_m near K_c = 20 is also well predicted and this phenomenon can be interpreted according to the vortex structure. For steady-oscillatory in-line combined flow cases, the vortex modes behave like a vortex street, exhibit a 'longitudinal wave' structure, and a vor- tex cluster shape corresponding to the ratios of U_m to U_0 which are of O (10^(-1)), O(1)and O (10), respectively. The effect on the prediction of forces on the flat plate from the disturbance component in a combined flow has been demon- strated qualitatively. In addition to this, the lock-in phenomenon of vortex shedding has been checked.展开更多
The stall is an unsteady flow phenomenon that always causes instabilities and low efficiency for pumps. This paper focuses on the unsteady flow structures and evolutions under two types of stall conditions in centrifu...The stall is an unsteady flow phenomenon that always causes instabilities and low efficiency for pumps. This paper focuses on the unsteady flow structures and evolutions under two types of stall conditions in centrifugal pump impellers. Two centrifugal pump impellers, one with 6 and another with 5 blades, are considered and a developed large-eddy simulation method is adopted. The results show that the alternative stall occurs in the impeller with 6 blades, while, the rotating stall is observed in that with 5 blades. The flow structure and the pressure fluctuation characteristics are further analyzed. For the alternative stall, the stall cells are fixed relative to the impeller, but a large vortex in the stalled passage is always swaying. The outlet vortex is generated from it, and then develops and sheds periodically. For the rotating stall, the stall cells first occur in the suction side of the blade. With the growth of the stall cells, the block area gradually increases until the inlet region is almost blocked, then moves to the pressure side with a continuous decay. When the rotating stall occurs, the amplitude of the pressure fluctuation is much larger than that under the alternative stall condition. The propagation of the stall cells has a significant effect on the pressure fluctuations in the impeller.展开更多
This paper investigates the unsteady structures and the hydrodynamics of cavitating flows. Experimental results are presented for a Clark-Y hydrofoil fixed at ?= 0?, 5? and 8?, for various cavitation numbers, from...This paper investigates the unsteady structures and the hydrodynamics of cavitating flows. Experimental results are presented for a Clark-Y hydrofoil fixed at ?= 0?, 5? and 8?, for various cavitation numbers, from subcavitating flow to supercavitation. The high-speed video camera and the particle image velocimetry(PIV) are applied to observe the transient flow structures. Statistics of the cavity lengths, the velocity and vorticity distributions, as well as the turbulent intensities are presented to quantify the unsteady process. Meanwhile, the dynamic measurement system is used to record the dynamic characteristics. The experimental results show that the flow structures and the hydrodynamics of the cavitation vary considerably with various combinations of angles of attack and cavitation numbers. Under various conditions, the cavitation can be generally grouped as the inception cavitation, the sheet cavitation, the cloud cavitation and the supercavitation. The cloud cavitation exhibits noticeable unsteady characteristics. Experimental evidence indicates that the hydrodynamics are clearly affected by the cavitating flow structures, the amplitude of the load fluctuation is much higher in the cloud cavitating cases.展开更多
A series of numerical analyses have been performed to investigate the flow structures in a narrow confined channel with 12 staggered circular impingement holes and one bigger exit hole. The flow enters the channel thr...A series of numerical analyses have been performed to investigate the flow structures in a narrow confined channel with 12 staggered circular impingement holes and one bigger exit hole. The flow enters the channel through the impingement holes and exits through the far end outlet. The flow fields corresponding to two jet Reynolds numbers (25000 and 65000) and three channel con- figurations with different ratios of the channel height to the impingement hole diameter (Zr 1, 3, 5) are analyzed by solving the Reynolds averaged Navier-Stokes equations with the realizable k-e turbulence model. Detailed flow field information including the secondary flow, the interaction between the jets and the cross flow, and flow distribution along the channel has been obtained. Comparisons between the numerical and experimental results of the flow fields at the four planes along the channel are performed to validate the numerical method. The calculated impingement pattern, high velocity flow distribution, low velocity separation region and vortices are in good agreement with the experimental data, implying the validity and effectiveness of the employed numerical approach for analyzing relevant flow field.展开更多
Studies on the high-lift mechanisms of butterfly gliding flights shed light on the design of the micro air vehicles.The flow field around a simplified Danaus plexippus model is investigated using the hydrogen bubble v...Studies on the high-lift mechanisms of butterfly gliding flights shed light on the design of the micro air vehicles.The flow field around a simplified Danaus plexippus model is investigated using the hydrogen bubble visualization and the Particle Image Velocimetry(PIV)techniques.There are three near-wall topological patterns with different Angles of Attack(AoAs):the separation bubble,the Leading-Edge Vortex(LEV)and the high Ao As flow.For the separation bubble pattern,two saddles and two foci form in the middle of the model.The features of the LEV pattern are the leading-edge separation lines.The topological characteristics of the separation lines are changed by the interaction between the LEV and the Wing-Tip Vortex(WTV).For the high Ao As flow pattern,four unstable foci are found at the forewing and the hindwing respectively.The angle between the trajectory of the WTV and the model increases with increasing Ao A even though the slope of the WTV angle versus Ao A curve declines at the moderate AoAs.展开更多
By means of ink trace visualization of the flows in conventional straight, positively curved and negatively curved cascades with tip clearance, and measurement of the aerodynamic parameters in transverse section, and...By means of ink trace visualization of the flows in conventional straight, positively curved and negatively curved cascades with tip clearance, and measurement of the aerodynamic parameters in transverse section, and by appling topology theory, the topological structures and vortex structure in the transverse section of a blade cascade were analyzed. Compared with conventional straight cascade, blade positive curving eliminates the separation line of the upper passage vortex, and leads the secondary vortex to change from close separation to open separation, while blade negative curving effects merely the positions of singular points and the intensities and scales of vortex.展开更多
River bending is the major effect responsible for bed topography and bank changes.In this study,fluid velocity(measured by a three-dimensional Doppler advanced point current meter)and bed topographical data have been ...River bending is the major effect responsible for bed topography and bank changes.In this study,fluid velocity(measured by a three-dimensional Doppler advanced point current meter)and bed topographical data have been collected in 40 sections of an experimental model.The whole flume was composed of an organic glass bend,upstream and downstream water tanks,two transition straight sections,a circulation pump,and a connection pipeline.Each section has been found to be characterized by a primary circulation and a small reverse circulation,with some sections even presenting three more or more circulation structures.The minimum circulation intensity has been detected in proximity to the top of the curved channel,while a region with small longitudinal velocity has been observed near the concave bank of each bend,corresponding to the flat bed formed after a short period of scouring.The maximum sediment deposition and scour depth in the presence of a uniform distribution of living flexible vegetation within 10 cm of the flume wall have been found to be smaller than those observed in the tests conducted without vegetation.展开更多
In practical problems, the variation of the free surface around a circular pier is severe. For the Fluent CFD, the multiphase flow models cannot be used together with the LES model to calculate the free surface. Prese...In practical problems, the variation of the free surface around a circular pier is severe. For the Fluent CFD, the multiphase flow models cannot be used together with the LES model to calculate the free surface. Present paper provides a two-step method which joints the 2D compressible ideal-gas equations and the LES model to calculate the 3D flow field with free surface around the pier. The effects of the free surface on the flow structures are studied in detail.展开更多
Inlet recirculation is proved as an effective way for centrifugal compressor surge margin extension,and is successively used in some engineering applications.Unfortunately its working mechanism is still not being well...Inlet recirculation is proved as an effective way for centrifugal compressor surge margin extension,and is successively used in some engineering applications.Unfortunately its working mechanism is still not being well understood,which leads to redesigning of inlet recirculation mostly by experience.Also,most study about inlet recirculation is steady to date.It is necessary to study surge margin extension mechanism about inlet recirculation.To expose the mechanism in detail,steady and unsteady numerical simulations were performed on a centrifugal compressor with and without inlet recirculation.The results showed that,with inlet recirculation,the inlet axial velocity is augmented,relative Mach number around blade tip leading edge area is significantly reduced and so is the flow angle.As the flow angle decreased,the incidence angle reduced which greatly improves the flow field inside the impeller.Moreover,inlet recirculation changes the blade loading around blade tip and restrains the flow separation on the blade suction side at the leading edge area.The unsteady results of static pressure around blade surface,entropy at inlet crossflow section and vorticity distributions at near tip span surface indicated that,at near stall condition,strong fluctuation exists in the vicinity of tip area due to the interaction between tip leakage flow and core flow.By inlet recirculation these strong flow fluctuations are eliminated so the flow stability is greatly enhanced.All these improvements mentioned above are the reason for inlet recirculation delays compressor stall.This research reveals the surge margin extension reason of inlet recirculation from an unsteady flow viewpoint and provides important reference for inlet recirculation structure design.展开更多
Wall-bounded turbulent flow involves the development of multi-scale turbulent eddies, as well as a sharply varying boundary layer. Its theoretical descriptions are yet phenomenological. We present here a new framework...Wall-bounded turbulent flow involves the development of multi-scale turbulent eddies, as well as a sharply varying boundary layer. Its theoretical descriptions are yet phenomenological. We present here a new framework called structural ensemble dynamics (SED), which aims at using systematically all relevant statistical properties of turbulent structures for a quantitative description of ensemble means. A new set of closure equations based on the SED approach for a turbulent channel flow is presented. SED order functions are defined, and numerically determined from data of direct numerical simulations (DNS). Computational results show that the new closure model reproduces accurately the solution of the original Navier-Stokes simulation, including the mean velocity profile, the kinetic energy of the streamwise velocity component, and every term in the energy budget equation. It is suggested that the SED-based studies of turbulent structure builds a bridge between the studies of physical mechanisms of turbulence and the development of accurate model equations for engineering predictions.展开更多
This study investigates the heterogeneous structure and its influence on drag coefficient for concurrent-up gas-solid flow. The energy-minimization multi-scale (EMMS) model is modified to simulate the variation of str...This study investigates the heterogeneous structure and its influence on drag coefficient for concurrent-up gas-solid flow. The energy-minimization multi-scale (EMMS) model is modified to simulate the variation of structure parameters with solids concentration, showing the tendency for particles to aggregate to form clusters and for fluid to pass around clusters. The global drag coefficient is resolved into that for the dense phase, for the dilute phase and for the so-called inter-phase, all of which can be obtained from their respective phase-specific structure parameters. The computational results show that the drag coefficients of the different phases are quite different, and the global drag coefficient calculated from the EMMS approach is much lower than that from the correlation of Wen and Yu. The simulation results demonstrate that the EMMS approach can well describe the heterogeneous flow structure, and is very promising for incorporation into the two-fluid model or the discrete particle model as the closure law for drag coefficient.展开更多
Gas-liquid two-phase flow widely exits in production and transportation of petroleum industry.Characterizing gas-liquid flow and measuring flow parameters represent challenges of great importance,which contribute to t...Gas-liquid two-phase flow widely exits in production and transportation of petroleum industry.Characterizing gas-liquid flow and measuring flow parameters represent challenges of great importance,which contribute to the recognition of flow regime and the optimal design of industrial equipment.In this paper,we propose a novel complex network-based deep learning method for characterizing gas-liquid flow.Firstly,we map the multichannel measurements to multiple limited penetrable visibility graphs(LPVGs)and obtain their degree sequences as the graph representation.Based on the degree distribution,we analyze the complicated flow behavior under different flow structures.Then,we design a dual-input convolutional neural network to fuse the raw signals and the graph representation of LPVGs for the classification of flow structures and measurement of gas void fraction.We implement the model with two parallel branches with the same structure,each corresponding to one input.Each branch consists of a channel-projection convolutional part,a spatial-temporal convolutional part,a dense block and an attention module.The outputs of the two branches are concatenated and fed into several full connected layers for the classification and measurement.At last,our method achieves an accuracy of 95.3%for the classification of flow structures,and a mean squared error of 0.0038 and a mean absolute percent error of 6.3%for the measurement of gas void fraction.Our method provides a promising solution for characterizing gas-liquid flow and measuring flow parameters.展开更多
In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge mod...In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge model are employed, Numerical wind tunnel technology based on computational fluid dynamics(CFD) is used, and the CFD models are set as stationary models. The Reynolds number of the flow, based on the inflow velocity and the height of the vehicle, is 1.9×10~6. The computations are conducted under three cases, train on the windward track on the bridge(WWC), train on the leeward track on the bridge(LWC) and train on the flat ground(FGC). Commercial software FLUENT is used and the mesh sensitivity research is carried out by three different grids: coarse, medium and fine. Results show that compared with FGC case, the side force coefficients of the head cars for the WWC and LWC cases increases by 14% and 29%, respectively; the coefficients of middle cars for the WWC and LWC increase by 32% and 10%, respectively; and that of the tail car increases by 45% for the WWC whereas decreases by 2% for the LWC case. The most notable thing is that the side force and the rolling moment of the head car are greater for the LWC, while the side force and the rolling moment of the middle car and the tail car are greater for the WWC. Comparing the velocity profiles at different locations, the flow is significantly influenced by the bridge-train system when the air is close to it. For the three cases(WWC, LWC and FGC), the pressure on the windward side of train is mostly positive while that of the leeward side is negative. The discrepancy of train's aerodynamic force is due to the different surface area of positive pressure and negative pressure zone. Many vortices are born on the leeward edge of the roofs. Theses vortices develop downstream, detach and dissipate into the wake region. The eddies develop irregularly, leading to a noticeably turbulent flow at leeward side of train.展开更多
The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Prov...The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Province, which is one of the most typical coastal aeolian distribution regions in China and famous for the tall and typical coastal transverse ridges. The measurement results show that, on the conditions of approximate wind velocities and same surface materials and environments, some changes happen to the structure of wind-sand flow with the increase of total sand transport rate on the crest of coastal transverse ridge. First, the sand transport rates of layers at different heights in the wind-sand flow increase, with the maximum increase at the height layer of 4-8cm. Second, the ratios of sand trans-port rates of layers at different heights to total sand transport rate decrease at the low height layer (0-4cm), but increase at the high height layer (4-60cm). Third, the distribution of the sand transport rate in the wind-sand flow can be expressed by an exponential function at the height layer of 0-40cm, but it changes from power function model to ex-ponential function model in the whole height layer (0-60cm) and changes into polynomial function model at the height layer of 40-60cm with the increase of total sand transport rate. Those changes have a close relationship with the limit of sand grain size of wind flow transporting and composition of sand grain size in the wind-sand flow.展开更多
This paper presents the results of an experimental study on the influences of floodplain impermeable groynes on flow structure, velocity, and water depth around the groyne(s). A wooden symmetrical compound channel w...This paper presents the results of an experimental study on the influences of floodplain impermeable groynes on flow structure, velocity, and water depth around the groyne(s). A wooden symmetrical compound channel was used. Groyne models with three different groyne relative lengths, 0.5, 0.75, and 1.0, were used on one floodplain with single and series arrangements. Analysis of the experimental results using the measured flow velocity and water depth values showed that flow structure, velocity, and water depth mainly depend on groyne relative length and the relative distance between series groynes. The flow velocity at the main channel centerline increased by about 40%, 60%, and 85%, and in other parts on the horizontal plane at the floodplain mid-water, depth by about 75%, 125%, and 175% of its original value in eases of one-side floodplain groyne(s) with relative lengths of 0.5, 0.75, and 1.0, respectively. The effective distance between two groynes in series arrangement ranges from 3 to 4 times the groyne length. Using an impermeable groyne with a large relative length in river floodplains increases the generation of eddy and roller zones downstream of the groyne, leading to more scouring and deposition. To avoid that, the groyne relative length must be kept below half the floodplain width,展开更多
基金Project supported by the National Basic Research Program of China (Grant No. 2009 CB724100)the National Natural Science Foundation of China (Grant No. 11172326)
文摘In a Mach 3.8 wind tunnel, both instantaneous and time-averaged flow structures of different scales around a blunt double-cone with or without supersonic film cooling were visualized via nano-tracer planar laser scattering (NPLS), which has a high spatiotemporal resolution. Three experimental cases with different injection mass flux rates were carried out. Many typical flow structures were clearly shown, such as shock waves, expansion fans, shear layers, mixing layers, and turbulent boundary layers. The analysis of two NPLS images with an interval of 5 us revealed the temporal evolution characteristics of flow structures. With matched pressures, the laminar length of the mixing layer was longer than that in the case with a larger mass flux rate, but the full covered region was shorter. Structures like K-H (Kelvin-Helmholtz) vortices were clearly seen in both flows. Without injection, the flow was similar to the supersonic flow over a backward- facing step, and the structures were relatively simpler, and there was a longer laminar region. Large scale structures such as hairpin vortices were visualized. In addition, the results were compared in part with the schlieren images captured by others under similar conditions.
基金project was supported by the National Natural Science Foundation of China (Grants 51306020, 11172040)the Natural Science Foundation of Beijing (Grant 3144034)the Excellent Young Scholars Research Fund of Beijing Institute of Technology
文摘The objective of this paper is to address the transient flow structures around a pitching hydrofoil by com- bining physical and numerical studies. In order to predict the dynamic behavior of the flow structure effectively, the Lagrangian coherent structures (LCS) defined by the ridges of the finite-time Lyapunov exponent (FTLE) are utilized under the framework of Navier-Stokes flow computations. In the numerical simulations, the k-w shear stress trans- port (SST) turbulence model, coupled with a two-equation F-Reo transition model, is used for the turbulence closure. Results are presented for a NACA66 hydrofoil undergoing slowly and rapidly pitching motions from 0° to 15° then back to 0° at a moderate Reynolds number Re = 7.5 × 105. The results reveal that the transient flow structures can be observed by the LCS method. For the slowly pitching case, it consists of five stages: quasi-steady and laminar, transition from laminar to turbulent, vortex development, large-scale vortex shedding, and reverting to laminar. The observation of LCS and Lagrangian particle tracers elucidates that the trailing edge vortex is nearly attached and stable during the vortex development stage and the interaction between the leading and trailing edge vortex caused by the adverse pres- sure gradient forces the vortexes to shed downstream during the large-scale vortex shedding stage, which corresponds to obvious fluctuations of the hydrodynamic response. For the rapidly pitching case, the inflection is hardly to be observed and the stall is delayed. The vortex formation, interaction, and shedding occurred once instead of being repeated three times, which is responsible for just one fluctuation in the hydrody- namic characteristics. The numerical results also show that the FTLE field has the potential to identify the transient flows, and the LCS can represent the divergence extent of infinite neighboring particles and capture the interface of the vortex region.
基金supported by the National Natural Science Foundation of China (10772017 and 10472011)BUAA-985 Foundation
文摘A bionic experimental platform was designed for the purpose of investigating time accurate three-dimensional flow field, using digital particle image velocimetry (DSPIV). The wake behind the flapping trail of a robotic fish model was studied at high spatial resolution. The study was performed in a water channel. A robot fish model was designed and built. The model was fixed onto a rigid support frame- work using a cable-supporting method, with twelve stretched wires. The entire tail of the model can perform prescribed motions in two degrees of freedom, mainly in carangiform mode, by driving its afterbody and lunate caudal fin respectively. The DSPIV system was set up to operate in a trans- lational manner, measuring velocity field in a series of parallel slices. Phase locked measurements were repeated for a number of runs, allowing reconstruction of phase average flow field. Vortex structures with phase history of the wake were obtained. The study reveals some new and complex three-dimensional flow structures in the wake of the fish, including "reverse hairpin vortex" and "reverse Karman S-H vortex rings", allowing insight into physics of this complex flow.
基金supported by the Scientific Research Project Fund of Middle East Technical University,under project # BAP-08-11-2016-044
文摘In this work,the dielectric-barrier-discharge plasma actuator was employed to study the flow structures induced by the plasma actuator over a flat plate and a wall-mounted hump.A phenomenological dielectric-barrier-discharge plasma model which regarded the plasma effect as the body force was implemented into the Navier–Stokes equations solved by the method of large eddy simulations.The results show that a series of vortex pairs,which indicated dipole formation and periodicity distribution were generated in the boundary layer when the plasma was applied to the flow over a flat plane.They would enhance the energy exchanged between the near wall region and the free stream.Besides,their spatial trajectories are deeply affected by the actuation strength.When the actuator was engaged in the flow over a wall-mounted hump,the vortex pairs were also produced,which was able to delay flow separation as well as to promote flow reattachment and reduce the generation of a vortex,achieving the goal of reducing dissipation and decreasing flow resistance.
基金The project supported by National Natural Science Foundation of China & LNM, Institute of Mechanics, CAS
文摘The evolution of wake structures and variation of the forces on a flat plate in harmonic oscillatory and in-line combined flows are obtained numerically by improved discrete vortex method. For the oscillatory oncoming flow cases, when K_c number varies from 2 to 40, the vortex pattern changes from a 'harmonic wave' shaped (in a range of small K_c numbers) to a slight inclined 'harmonic wave' shaped (in a range of moderate K_c numbers), then to inclined vortex clusters with an angle of 50 ° to the oncoming flow direction (at K_c = 20), at last, as K_c number becomes large, the vortex pattern is like a normal Karman vortex street. The well predicted drag and inertia force coefficients are obtained, which are more close to the results of Keulegan & Carpenter's experiment as compared with previous vortex simulation by other au- thors. The existence of minimum point of inertia force coefficient C_m near K_c = 20 is also well predicted and this phenomenon can be interpreted according to the vortex structure. For steady-oscillatory in-line combined flow cases, the vortex modes behave like a vortex street, exhibit a 'longitudinal wave' structure, and a vor- tex cluster shape corresponding to the ratios of U_m to U_0 which are of O (10^(-1)), O(1)and O (10), respectively. The effect on the prediction of forces on the flat plate from the disturbance component in a combined flow has been demon- strated qualitatively. In addition to this, the lock-in phenomenon of vortex shedding has been checked.
基金Project supported by the National Nature Science Foundation of China(Grant Nos.51709234,51779226)the Zhejiang Provincial Natural Science Foundation of China(Grant No.LQ17E090005)
文摘The stall is an unsteady flow phenomenon that always causes instabilities and low efficiency for pumps. This paper focuses on the unsteady flow structures and evolutions under two types of stall conditions in centrifugal pump impellers. Two centrifugal pump impellers, one with 6 and another with 5 blades, are considered and a developed large-eddy simulation method is adopted. The results show that the alternative stall occurs in the impeller with 6 blades, while, the rotating stall is observed in that with 5 blades. The flow structure and the pressure fluctuation characteristics are further analyzed. For the alternative stall, the stall cells are fixed relative to the impeller, but a large vortex in the stalled passage is always swaying. The outlet vortex is generated from it, and then develops and sheds periodically. For the rotating stall, the stall cells first occur in the suction side of the blade. With the growth of the stall cells, the block area gradually increases until the inlet region is almost blocked, then moves to the pressure side with a continuous decay. When the rotating stall occurs, the amplitude of the pressure fluctuation is much larger than that under the alternative stall condition. The propagation of the stall cells has a significant effect on the pressure fluctuations in the impeller.
基金Project supported by the National Natural Science Foundation of China(Grant No.51306020)the Natural Science Foundation of Beijing(Grant No.3144034)
文摘This paper investigates the unsteady structures and the hydrodynamics of cavitating flows. Experimental results are presented for a Clark-Y hydrofoil fixed at ?= 0?, 5? and 8?, for various cavitation numbers, from subcavitating flow to supercavitation. The high-speed video camera and the particle image velocimetry(PIV) are applied to observe the transient flow structures. Statistics of the cavity lengths, the velocity and vorticity distributions, as well as the turbulent intensities are presented to quantify the unsteady process. Meanwhile, the dynamic measurement system is used to record the dynamic characteristics. The experimental results show that the flow structures and the hydrodynamics of the cavitation vary considerably with various combinations of angles of attack and cavitation numbers. Under various conditions, the cavitation can be generally grouped as the inception cavitation, the sheet cavitation, the cloud cavitation and the supercavitation. The cloud cavitation exhibits noticeable unsteady characteristics. Experimental evidence indicates that the hydrodynamics are clearly affected by the cavitating flow structures, the amplitude of the load fluctuation is much higher in the cloud cavitating cases.
基金supported by the National Natural Science Foundation of China(No.51206180)the Natural Science Basic Research Plan in Shaanxi Province of China(No.2014JQ7276)
文摘A series of numerical analyses have been performed to investigate the flow structures in a narrow confined channel with 12 staggered circular impingement holes and one bigger exit hole. The flow enters the channel through the impingement holes and exits through the far end outlet. The flow fields corresponding to two jet Reynolds numbers (25000 and 65000) and three channel con- figurations with different ratios of the channel height to the impingement hole diameter (Zr 1, 3, 5) are analyzed by solving the Reynolds averaged Navier-Stokes equations with the realizable k-e turbulence model. Detailed flow field information including the secondary flow, the interaction between the jets and the cross flow, and flow distribution along the channel has been obtained. Comparisons between the numerical and experimental results of the flow fields at the four planes along the channel are performed to validate the numerical method. The calculated impingement pattern, high velocity flow distribution, low velocity separation region and vortices are in good agreement with the experimental data, implying the validity and effectiveness of the employed numerical approach for analyzing relevant flow field.
基金funded by the National Natural Science Foundation of China (No. 11721202)。
文摘Studies on the high-lift mechanisms of butterfly gliding flights shed light on the design of the micro air vehicles.The flow field around a simplified Danaus plexippus model is investigated using the hydrogen bubble visualization and the Particle Image Velocimetry(PIV)techniques.There are three near-wall topological patterns with different Angles of Attack(AoAs):the separation bubble,the Leading-Edge Vortex(LEV)and the high Ao As flow.For the separation bubble pattern,two saddles and two foci form in the middle of the model.The features of the LEV pattern are the leading-edge separation lines.The topological characteristics of the separation lines are changed by the interaction between the LEV and the Wing-Tip Vortex(WTV).For the high Ao As flow pattern,four unstable foci are found at the forewing and the hindwing respectively.The angle between the trajectory of the WTV and the model increases with increasing Ao A even though the slope of the WTV angle versus Ao A curve declines at the moderate AoAs.
文摘By means of ink trace visualization of the flows in conventional straight, positively curved and negatively curved cascades with tip clearance, and measurement of the aerodynamic parameters in transverse section, and by appling topology theory, the topological structures and vortex structure in the transverse section of a blade cascade were analyzed. Compared with conventional straight cascade, blade positive curving eliminates the separation line of the upper passage vortex, and leads the secondary vortex to change from close separation to open separation, while blade negative curving effects merely the positions of singular points and the intensities and scales of vortex.
基金supported in part by the Special Fund for Basic Scientific Research Business Expenses of Central Public Welfare Scientific Research Institutes under Grant TKS20210103the Open Fund of Key Laboratory of Ocean Observation Technology,Ministry of Natural Resources of China(2021klootA06).
文摘River bending is the major effect responsible for bed topography and bank changes.In this study,fluid velocity(measured by a three-dimensional Doppler advanced point current meter)and bed topographical data have been collected in 40 sections of an experimental model.The whole flume was composed of an organic glass bend,upstream and downstream water tanks,two transition straight sections,a circulation pump,and a connection pipeline.Each section has been found to be characterized by a primary circulation and a small reverse circulation,with some sections even presenting three more or more circulation structures.The minimum circulation intensity has been detected in proximity to the top of the curved channel,while a region with small longitudinal velocity has been observed near the concave bank of each bend,corresponding to the flat bed formed after a short period of scouring.The maximum sediment deposition and scour depth in the presence of a uniform distribution of living flexible vegetation within 10 cm of the flume wall have been found to be smaller than those observed in the tests conducted without vegetation.
文摘In practical problems, the variation of the free surface around a circular pier is severe. For the Fluent CFD, the multiphase flow models cannot be used together with the LES model to calculate the free surface. Present paper provides a two-step method which joints the 2D compressible ideal-gas equations and the LES model to calculate the 3D flow field with free surface around the pier. The effects of the free surface on the flow structures are studied in detail.
文摘Inlet recirculation is proved as an effective way for centrifugal compressor surge margin extension,and is successively used in some engineering applications.Unfortunately its working mechanism is still not being well understood,which leads to redesigning of inlet recirculation mostly by experience.Also,most study about inlet recirculation is steady to date.It is necessary to study surge margin extension mechanism about inlet recirculation.To expose the mechanism in detail,steady and unsteady numerical simulations were performed on a centrifugal compressor with and without inlet recirculation.The results showed that,with inlet recirculation,the inlet axial velocity is augmented,relative Mach number around blade tip leading edge area is significantly reduced and so is the flow angle.As the flow angle decreased,the incidence angle reduced which greatly improves the flow field inside the impeller.Moreover,inlet recirculation changes the blade loading around blade tip and restrains the flow separation on the blade suction side at the leading edge area.The unsteady results of static pressure around blade surface,entropy at inlet crossflow section and vorticity distributions at near tip span surface indicated that,at near stall condition,strong fluctuation exists in the vicinity of tip area due to the interaction between tip leakage flow and core flow.By inlet recirculation these strong flow fluctuations are eliminated so the flow stability is greatly enhanced.All these improvements mentioned above are the reason for inlet recirculation delays compressor stall.This research reveals the surge margin extension reason of inlet recirculation from an unsteady flow viewpoint and provides important reference for inlet recirculation structure design.
基金supported by the National Natural Science Foundation of China (90716008)the MOST under 973 project (2009CB724100)
文摘Wall-bounded turbulent flow involves the development of multi-scale turbulent eddies, as well as a sharply varying boundary layer. Its theoretical descriptions are yet phenomenological. We present here a new framework called structural ensemble dynamics (SED), which aims at using systematically all relevant statistical properties of turbulent structures for a quantitative description of ensemble means. A new set of closure equations based on the SED approach for a turbulent channel flow is presented. SED order functions are defined, and numerically determined from data of direct numerical simulations (DNS). Computational results show that the new closure model reproduces accurately the solution of the original Navier-Stokes simulation, including the mean velocity profile, the kinetic energy of the streamwise velocity component, and every term in the energy budget equation. It is suggested that the SED-based studies of turbulent structure builds a bridge between the studies of physical mechanisms of turbulence and the development of accurate model equations for engineering predictions.
基金Supported by the National Key Program for Developing Basic Sciences of China (No. G1999022103) and the National Natural Science Foundation of China (No. 20176059).
文摘This study investigates the heterogeneous structure and its influence on drag coefficient for concurrent-up gas-solid flow. The energy-minimization multi-scale (EMMS) model is modified to simulate the variation of structure parameters with solids concentration, showing the tendency for particles to aggregate to form clusters and for fluid to pass around clusters. The global drag coefficient is resolved into that for the dense phase, for the dilute phase and for the so-called inter-phase, all of which can be obtained from their respective phase-specific structure parameters. The computational results show that the drag coefficients of the different phases are quite different, and the global drag coefficient calculated from the EMMS approach is much lower than that from the correlation of Wen and Yu. The simulation results demonstrate that the EMMS approach can well describe the heterogeneous flow structure, and is very promising for incorporation into the two-fluid model or the discrete particle model as the closure law for drag coefficient.
基金supported by the National Natural Science Foundation of China under Grants 61922062 and 61873181。
文摘Gas-liquid two-phase flow widely exits in production and transportation of petroleum industry.Characterizing gas-liquid flow and measuring flow parameters represent challenges of great importance,which contribute to the recognition of flow regime and the optimal design of industrial equipment.In this paper,we propose a novel complex network-based deep learning method for characterizing gas-liquid flow.Firstly,we map the multichannel measurements to multiple limited penetrable visibility graphs(LPVGs)and obtain their degree sequences as the graph representation.Based on the degree distribution,we analyze the complicated flow behavior under different flow structures.Then,we design a dual-input convolutional neural network to fuse the raw signals and the graph representation of LPVGs for the classification of flow structures and measurement of gas void fraction.We implement the model with two parallel branches with the same structure,each corresponding to one input.Each branch consists of a channel-projection convolutional part,a spatial-temporal convolutional part,a dense block and an attention module.The outputs of the two branches are concatenated and fed into several full connected layers for the classification and measurement.At last,our method achieves an accuracy of 95.3%for the classification of flow structures,and a mean squared error of 0.0038 and a mean absolute percent error of 6.3%for the measurement of gas void fraction.Our method provides a promising solution for characterizing gas-liquid flow and measuring flow parameters.
基金Project(U1534210)supported by the National Natural Science Foundation of ChinaProject(14JJ1003)supported by the Natural Science Foundation of Hunan Province,China+2 种基金Project(2015CX003)supported by the Project of Innovation-driven Plan in Central South University,ChinaProject(14JC1003)supported by the Natural Science Foundation of Hunan Province,ChinaProject(2015T002-A)supported by the Technological Research and Development program of China Railways Cooperation
文摘In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge model are employed, Numerical wind tunnel technology based on computational fluid dynamics(CFD) is used, and the CFD models are set as stationary models. The Reynolds number of the flow, based on the inflow velocity and the height of the vehicle, is 1.9×10~6. The computations are conducted under three cases, train on the windward track on the bridge(WWC), train on the leeward track on the bridge(LWC) and train on the flat ground(FGC). Commercial software FLUENT is used and the mesh sensitivity research is carried out by three different grids: coarse, medium and fine. Results show that compared with FGC case, the side force coefficients of the head cars for the WWC and LWC cases increases by 14% and 29%, respectively; the coefficients of middle cars for the WWC and LWC increase by 32% and 10%, respectively; and that of the tail car increases by 45% for the WWC whereas decreases by 2% for the LWC case. The most notable thing is that the side force and the rolling moment of the head car are greater for the LWC, while the side force and the rolling moment of the middle car and the tail car are greater for the WWC. Comparing the velocity profiles at different locations, the flow is significantly influenced by the bridge-train system when the air is close to it. For the three cases(WWC, LWC and FGC), the pressure on the windward side of train is mostly positive while that of the leeward side is negative. The discrepancy of train's aerodynamic force is due to the different surface area of positive pressure and negative pressure zone. Many vortices are born on the leeward edge of the roofs. Theses vortices develop downstream, detach and dissipate into the wake region. The eddies develop irregularly, leading to a noticeably turbulent flow at leeward side of train.
基金Under the auspices of National Natural Science Foundation of China (No 40571019)
文摘The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Province, which is one of the most typical coastal aeolian distribution regions in China and famous for the tall and typical coastal transverse ridges. The measurement results show that, on the conditions of approximate wind velocities and same surface materials and environments, some changes happen to the structure of wind-sand flow with the increase of total sand transport rate on the crest of coastal transverse ridge. First, the sand transport rates of layers at different heights in the wind-sand flow increase, with the maximum increase at the height layer of 4-8cm. Second, the ratios of sand trans-port rates of layers at different heights to total sand transport rate decrease at the low height layer (0-4cm), but increase at the high height layer (4-60cm). Third, the distribution of the sand transport rate in the wind-sand flow can be expressed by an exponential function at the height layer of 0-40cm, but it changes from power function model to ex-ponential function model in the whole height layer (0-60cm) and changes into polynomial function model at the height layer of 40-60cm with the increase of total sand transport rate. Those changes have a close relationship with the limit of sand grain size of wind flow transporting and composition of sand grain size in the wind-sand flow.
文摘This paper presents the results of an experimental study on the influences of floodplain impermeable groynes on flow structure, velocity, and water depth around the groyne(s). A wooden symmetrical compound channel was used. Groyne models with three different groyne relative lengths, 0.5, 0.75, and 1.0, were used on one floodplain with single and series arrangements. Analysis of the experimental results using the measured flow velocity and water depth values showed that flow structure, velocity, and water depth mainly depend on groyne relative length and the relative distance between series groynes. The flow velocity at the main channel centerline increased by about 40%, 60%, and 85%, and in other parts on the horizontal plane at the floodplain mid-water, depth by about 75%, 125%, and 175% of its original value in eases of one-side floodplain groyne(s) with relative lengths of 0.5, 0.75, and 1.0, respectively. The effective distance between two groynes in series arrangement ranges from 3 to 4 times the groyne length. Using an impermeable groyne with a large relative length in river floodplains increases the generation of eddy and roller zones downstream of the groyne, leading to more scouring and deposition. To avoid that, the groyne relative length must be kept below half the floodplain width,