In this article,the method of volume averaging of flow in porous media and averaging theorem was employed to the Navier-Stokes equations.The total drag force per unit volume was considered as the bulk damping resistan...In this article,the method of volume averaging of flow in porous media and averaging theorem was employed to the Navier-Stokes equations.The total drag force per unit volume was considered as the bulk damping resistance due to the porous structure（i.e.,the Darcy resistance）and the resistance due to the inertia force（i.e.,the Forchheimer resistance）,then the Brinkman-Forchheimer extended Darcy model was obtained from average momentum equation.The structural flow in a pipe containing porous medium saturated with power law fluid was applied in the present study.A new theoretical analysis of fully developed non-Darcy flow in a pipe containing a fibrous medium saturated with power law fluid was conducted.Using the integration method,the boundary layer solutions were obtained for flow core velocity and its radius.These theoretical solutions were used to analyze the effects of the Darcy number and inertia parameter on the axial velocity profile in the porous medium pipe.展开更多
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.展开更多
The modern definition of the wave concept,which is based on the functional connection between the parameters of the spatial structure of an instantaneous flow pattern and the characteristics of the temporal variabilit...The modern definition of the wave concept,which is based on the functional connection between the parameters of the spatial structure of an instantaneous flow pattern and the characteristics of the temporal variability at a given point,is discussed.The dispersion relation for 2D plane periodic perturbations on the surface of viscous stratified fluid is selected as the characteristic function defining the wave motion.Using the theory of singular perturbations,a method for calculating complete solutions to the dispersion relations of periodic flows,including regular wave and singular ligament solutions is presented.Properties of the complete exact solution of the dispersion relation containing regular and singular functions are compared with asymptotic solutions.In limiting cases,obtained dispersion relations are matched with well⁃known expressions for waves in homogeneous viscous and ideal liquids.展开更多
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.展开更多
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 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.展开更多
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.展开更多
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.展开更多
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.展开更多
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,展开更多
The marine ecosystem of the Jiaozhou Bay has degraded significantly in fisheries productivity and its ecological roles as spawning and nursery ground for many species of commercial importance has been declining in rec...The marine ecosystem of the Jiaozhou Bay has degraded significantly in fisheries productivity and its ecological roles as spawning and nursery ground for many species of commercial importance has been declining in recent years. A mass-balanced trophic model was developed using Ecopath with Ecosim to evaluate the trophic structure of the Jiaozhou Bay for improving ecosystem management. The model were parameterized based on the fisheries survey data in the Jiaozhou Bay in 2011, including 23 species groups and one detritus group according to their ecological roles. The trophic levels of these ecological groups ranged from 1（primary producers and detritus） to4.3（large demersal fishes）. The estimated total system throughput was 12 917.10 t/（km^2·a）, with 74.59% and25.41% contribution of the total energy flows from phytoplankton and detritus, respectively. Network analyses showed that the overall transfer efficiency of the ecosystem was 14.4%, and the mean transfer efficiency was 14.5%for grazing food chain and 13.9% for detritus food chain. The system omnivory index（SOI）, Finn＇s cycled index（FCI） and connectance index（CI） were relatively low in this area while the total primary production/total respiration（TPP/TR） was high, indicating an immature and unstable status of the Jiaozhou Bay ecosystem. Mixed trophic impact analysis revealed that the cultured shellfish had substantial negative impacts on most functional groups. This study contributed to ecosystem-level evaluation and management planning of the Jiaozhou Bay ecosystem.展开更多
Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,bound...Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,boundary layer thickness and Stanton number demonstrate that the flow changes dramatically in the shock wave/turbulent boundary layer interaction area.It is found that the steady jet has no effect on suppressing flow separation unexpectedly,but increases its spatial scale instead.Instantaneous flow structures show that the turbulence amplification can be observed after the application of flow control,and abundant Görtler-like vorticities appear,but the strength of the main shock decreases.Analyzing the wall fluctuating pressure signals using weighted power spectral density,we found an interesting thing.That is,although the low-frequency oscillation phenomenon induced by separation shock is suppressed by the steady jet,wall fluctuating pressure beneath the jet shock is oscillating at a frequency lower than 0.1u∞/δref.Results of coherent and intermittency factor reveal that it is related to the backand-forth movement of the jet shock itself.展开更多
Experimental studies which focus on flow visualization and the velocity field of a supersonic laminar/turbulent flow over a compression ramp were carried out in a Mach 3.0 wind tunnel. Fine flow structures and velocit...Experimental studies which focus on flow visualization and the velocity field of a supersonic laminar/turbulent flow over a compression ramp were carried out in a Mach 3.0 wind tunnel. Fine flow structures and velocity field structures were obtained via NPLS （nanoparticle-tracer planar laser scattering） and PIV （particle image velocimetry） techniques, time- averaged flow structures were researched, and spatiotemporal evolutions of transient flow structures were analyzed. The flow visualization results indicated that when the ramp angles were 25~, a typical separation occurred in the laminar flow, some typical flow structures such as shock induced by the boundary layer, separation shock, reversed flow and reattachment shock were visible clearly. While a certain extent separation occurred in turbulent flow, the separation region was much smaller. When the ramp angles were 28~, laminar flow separated further, and the separation region expanded evidently, flow structures in the separation region were complex. While a typical separation occurred in turbulent flow, reversed flow structures were significant, flow structures in the separation region were relatively simple. The experimental results of velocity field were corresponding to flow visualization, and the velocity field structures of both compression ramp flows agreed with the flow structures well. There were three layered structures in the U component velocity, and the V component velocity appeared like an oblique ＂v＂. Some differences between these two compression ramp flows can be observed in the velocity profiles of the shear layer and the shearing intensity.展开更多
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.展开更多
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.展开更多
Unsteady-state operation has been widely applied in chemical engineering, such as optimizing a process, increasing yield and saving energy, etc. But the knowledge of the flow characteristics in bubble column reactors(...Unsteady-state operation has been widely applied in chemical engineering, such as optimizing a process, increasing yield and saving energy, etc. But the knowledge of the flow characteristics in bubble column reactors(BCRs) under unsteady state control is far from enough. In order to study the flow structures in this operation, the volume of fluid (VOF) model and the standard k-ε model to simulate the evolution of gas-liquid flow in BCRs under the start-up state are combined. For both the symmetry and asymmetry flow, the layout of the gas-inlets, the gas-in velocity, the liquid viscosity and the aspect ratio of the BCR all have effects on the liquid velocity distribution. The simulation results could provide some information for the design and scale-up of the BCRs.展开更多
Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate ...Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.50474021,51074102)
文摘In this article,the method of volume averaging of flow in porous media and averaging theorem was employed to the Navier-Stokes equations.The total drag force per unit volume was considered as the bulk damping resistance due to the porous structure（i.e.,the Darcy resistance）and the resistance due to the inertia force（i.e.,the Forchheimer resistance）,then the Brinkman-Forchheimer extended Darcy model was obtained from average momentum equation.The structural flow in a pipe containing porous medium saturated with power law fluid was applied in the present study.A new theoretical analysis of fully developed non-Darcy flow in a pipe containing a fibrous medium saturated with power law fluid was conducted.Using the integration method,the boundary layer solutions were obtained for flow core velocity and its radius.These theoretical solutions were used to analyze the effects of the Darcy number and inertia parameter on the axial velocity profile in the porous medium pipe.
基金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.
基金Sponsored by Ministry of Science and Higher Education within the Framework of Russian State Assignment(Grant No.124012500442⁃3).
文摘The modern definition of the wave concept,which is based on the functional connection between the parameters of the spatial structure of an instantaneous flow pattern and the characteristics of the temporal variability at a given point,is discussed.The dispersion relation for 2D plane periodic perturbations on the surface of viscous stratified fluid is selected as the characteristic function defining the wave motion.Using the theory of singular perturbations,a method for calculating complete solutions to the dispersion relations of periodic flows,including regular wave and singular ligament solutions is presented.Properties of the complete exact solution of the dispersion relation containing regular and singular functions are compared with asymptotic solutions.In limiting cases,obtained dispersion relations are matched with well⁃known expressions for waves in homogeneous viscous and ideal liquids.
文摘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.
基金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.
基金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.
基金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.
基金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.
基金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.
文摘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,
基金The Public Science and Technology Research Funds Projects of Ministry of Agriculture under contract No.201303050-02the Scientific and Technological Innovation Project for the Qingdao National Laboratory for Marine Science and Technology under contract No.2015ASKJ02+2 种基金the Fundamental Research Funds for the Central Universities under contract No.201262004the National Natural Science Foundation of China under contract No.41006083the Shandong Provincial Natural Science Foundation,China under contract No.ZR2010DQ026
文摘The marine ecosystem of the Jiaozhou Bay has degraded significantly in fisheries productivity and its ecological roles as spawning and nursery ground for many species of commercial importance has been declining in recent years. A mass-balanced trophic model was developed using Ecopath with Ecosim to evaluate the trophic structure of the Jiaozhou Bay for improving ecosystem management. The model were parameterized based on the fisheries survey data in the Jiaozhou Bay in 2011, including 23 species groups and one detritus group according to their ecological roles. The trophic levels of these ecological groups ranged from 1（primary producers and detritus） to4.3（large demersal fishes）. The estimated total system throughput was 12 917.10 t/（km^2·a）, with 74.59% and25.41% contribution of the total energy flows from phytoplankton and detritus, respectively. Network analyses showed that the overall transfer efficiency of the ecosystem was 14.4%, and the mean transfer efficiency was 14.5%for grazing food chain and 13.9% for detritus food chain. The system omnivory index（SOI）, Finn＇s cycled index（FCI） and connectance index（CI） were relatively low in this area while the total primary production/total respiration（TPP/TR） was high, indicating an immature and unstable status of the Jiaozhou Bay ecosystem. Mixed trophic impact analysis revealed that the cultured shellfish had substantial negative impacts on most functional groups. This study contributed to ecosystem-level evaluation and management planning of the Jiaozhou Bay ecosystem.
基金supported by the National Natural Science Foundation of China(Grant Nos.12202488 and 12002377)the Natural Science Program of National University of Defense Technology(Grant No.ZK22-30)Independent Cultivation Project for Young Talents in College of Aerospace Science and Engineering.
文摘Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,boundary layer thickness and Stanton number demonstrate that the flow changes dramatically in the shock wave/turbulent boundary layer interaction area.It is found that the steady jet has no effect on suppressing flow separation unexpectedly,but increases its spatial scale instead.Instantaneous flow structures show that the turbulence amplification can be observed after the application of flow control,and abundant Görtler-like vorticities appear,but the strength of the main shock decreases.Analyzing the wall fluctuating pressure signals using weighted power spectral density,we found an interesting thing.That is,although the low-frequency oscillation phenomenon induced by separation shock is suppressed by the steady jet,wall fluctuating pressure beneath the jet shock is oscillating at a frequency lower than 0.1u∞/δref.Results of coherent and intermittency factor reveal that it is related to the backand-forth movement of the jet shock itself.
基金supported by the National Natural Science Foundation of China(Grant Nos.11172326 and 11302256)
文摘Experimental studies which focus on flow visualization and the velocity field of a supersonic laminar/turbulent flow over a compression ramp were carried out in a Mach 3.0 wind tunnel. Fine flow structures and velocity field structures were obtained via NPLS （nanoparticle-tracer planar laser scattering） and PIV （particle image velocimetry） techniques, time- averaged flow structures were researched, and spatiotemporal evolutions of transient flow structures were analyzed. The flow visualization results indicated that when the ramp angles were 25~, a typical separation occurred in the laminar flow, some typical flow structures such as shock induced by the boundary layer, separation shock, reversed flow and reattachment shock were visible clearly. While a certain extent separation occurred in turbulent flow, the separation region was much smaller. When the ramp angles were 28~, laminar flow separated further, and the separation region expanded evidently, flow structures in the separation region were complex. While a typical separation occurred in turbulent flow, reversed flow structures were significant, flow structures in the separation region were relatively simple. The experimental results of velocity field were corresponding to flow visualization, and the velocity field structures of both compression ramp flows agreed with the flow structures well. There were three layered structures in the U component velocity, and the V component velocity appeared like an oblique ＂v＂. Some differences between these two compression ramp flows can be observed in the velocity profiles of the shear layer and the shearing intensity.
基金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.
基金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.
文摘Unsteady-state operation has been widely applied in chemical engineering, such as optimizing a process, increasing yield and saving energy, etc. But the knowledge of the flow characteristics in bubble column reactors(BCRs) under unsteady state control is far from enough. In order to study the flow structures in this operation, the volume of fluid (VOF) model and the standard k-ε model to simulate the evolution of gas-liquid flow in BCRs under the start-up state are combined. For both the symmetry and asymmetry flow, the layout of the gas-inlets, the gas-in velocity, the liquid viscosity and the aspect ratio of the BCR all have effects on the liquid velocity distribution. The simulation results could provide some information for the design and scale-up of the BCRs.
基金supported by the National Natural Science Foundation of China (10872096)the Open Fund of State Key Laboratory of Explosion Science and Technology, Beijing University of Science and Technology (KFJJ09-13)
文摘Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.
