The size of pores or the grille spacing of water–sediment separation structures directly affects their regulation effect on the debris flow performance.A suitable pore size or grille spacing can effectively improve t...The size of pores or the grille spacing of water–sediment separation structures directly affects their regulation effect on the debris flow performance.A suitable pore size or grille spacing can effectively improve the water–sediment separation ability of the structure.The new funnel-type grating water–sediment separation structure(FGWSS)combines vertical and horizontal structures and provides a satisfactory water–sediment separation effect.However,the regulation effect of the grille spacing of the structure on the debris flow performance has not been studied.The regulation effect of the structure grille spacing on the debris flow performance is studied through a flume test,and the optimal structure grille spacing is obtained.An empirical equation of the relationship between the relative grille spacing of the structure and the sediment separation rate is established.Finally,the influence of the water–sediment separation structure on the regulation effect of debris flows is examined from two aspects:external factors(properties of debris flows)and internal factors(structural factors).The experimental results show that the gradation characteristics of solid particles in debris flows constitute a key factor affecting the regulation effect of the structure on the debris flow performance.The optimum grille spacing of the FGWSS matches the particle size corresponding to the material distribution curves d85~d90 of the debris flow.The total separation rate of debris flow particles is related to the grille spacing of the structure and the content of coarse and fine particles in the debris flow.展开更多
Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the ef...Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.展开更多
An experimental system for forming a rotating paraboloid shaped shallow water with a free surface was conducted to study coherent vortex structures in a differentially rotating quasi two dimensional zonal flow.Flow...An experimental system for forming a rotating paraboloid shaped shallow water with a free surface was conducted to study coherent vortex structures in a differentially rotating quasi two dimensional zonal flow.Flow visualization and laser light scattering techniques were used to obtain the information of spatial flow patterns.Experimental results show that the coexistence of Coriolis effect and strong shear in latitudinal zones may lead to formation of coherent vortices.Power spectra analysis and photographs which were taken in a reference frame rotating with the observed vortices also justified the emergence,drift and evolution of persistent vortices on the large scale.Locked vortex state manifests the cyclone and anticyclone asymmetry.展开更多
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 theoretical research on the propulsive principle of aquatic animal becomes more important and attracted more researchers to make efforts on it. In the present study, a computational fluid dynamic (CFD) simulatio...The theoretical research on the propulsive principle of aquatic animal becomes more important and attracted more researchers to make efforts on it. In the present study, a computational fluid dynamic (CFD) simulation of a three-dimensional traveling-wave undulations body of tuna has been developed to investigate the fluid flow features and vorticity structures around this body when moving in a straight line. The undulation only takes place in the posterior half of the fish, and the tuna-tail is considered as a lunate fin oscillating with the mode combined swaying with yawing. A Reynolds-averaged Navier-Stokes (RANS) equation is developed, employing a control-volume method and a k-omega SST turbulent model; meanwhile an unstructured tetrahedral grid, which is generated for the three-dimensional geometry, is used based on the deformation of the hind parts of the body and corresponding movement of the tail. We calculated the hydrodynamic performance of tuna-like body when a tuna swims in a uniform velocity, and compared the input power coefficient, output power coefficient and propulsive efficiency of the oscillating tuna-tail with or without body vortex shedding. Additionally, the load distribution on the body, flow features and vorticity structures around the body were demonstrated. The effect of interaction between the body-generated vortices and the tail-generated vorticity on the hydrodynamic performance can be obtained.展开更多
According to the analysis of the turbulent intensity level around the high-speed train, the maximum turbulent intensity ranges from 0.2 to 0.5 which belongs to high turbulent flow. The flow field distribution law was ...According to the analysis of the turbulent intensity level around the high-speed train, the maximum turbulent intensity ranges from 0.2 to 0.5 which belongs to high turbulent flow. The flow field distribution law was studied and eight types of flow regions were proposed. They are high pressure with air stagnant region, pressure decreasing with air accelerating region, low pressure with high air flow velocity region I, turbulent region, steady flow region, low pressure with high air flow velocity region II,pressure increasing with air decelerating region and wake region. The analysis of the vortex structure around the train shows that the vortex is mainly induced by structures with complex mutation and large curvature change. The head and rear of train, the underbody structure, the carriage connection section and the wake region are the main vortex generating sources while the train body with even cross-section has rare vortexes. The wake structure development law studied lays foundation for the train drag reduction.展开更多
Work on dynamic topology optimization of engineering structures for vibration suppression has mainly addressed the maximization of eigenfrequencies and gaps between consecutive eigenfrequencies of free vibration, mini...Work on dynamic topology optimization of engineering structures for vibration suppression has mainly addressed the maximization of eigenfrequencies and gaps between consecutive eigenfrequencies of free vibration, minimization of the dynamic compliance subject to forced vibration, and minimization of the structural frequency response. A dynamic topology optimization method of bi-material plate structures is presented based on power flow analysis. Topology optimization problems formulated directly with the design objective of minimizing the power flow response are dealt with. In comparison to the displacement or velocity response, the power flow response takes not only the amplitude of force and velocity into account, but also the phase relationship of the two vector quantities. The complex expression of power flow response is derived based on time-harmonic external mechanical loading and Rayleigh damping. The mathematical formulation of topology optimization is established based on power flow response and bi-material solid isotropic material with penalization(SIMP) model. Computational optimization procedure is developed by using adjoint design sensitivity analysis and the method of moving asymptotes(MMA). Several numerical examples are presented for bi-material plate structures with different loading frequencies, which verify the feasibility and effectiveness of this method. Additionally, optimum results between topological design of minimum power flow response and minimum dynamic compliance are compared, showing that the present method has strong adaptability for structural dynamic topology optimization problems. The proposed research provides a more accurate and effective approach for dynamic topology optimization of vibrating structures.展开更多
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.展开更多
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.展开更多
Pipe-in-pipe(PIP)structures are widely used in offshore oil and gas pipelines to settle thermal insulation issues.A PIP structure system usually consists of two concentric pipes and one softer layer for thermal insula...Pipe-in-pipe(PIP)structures are widely used in offshore oil and gas pipelines to settle thermal insulation issues.A PIP structure system usually consists of two concentric pipes and one softer layer for thermal insulation consideration.The total response of the system is related to the dynamics of both pipes and the interactions between these two concentric pipes.In the current work,a theoretical model for flow-induced vibrations of a PIP structure system is proposed and analyzed in the presence of an internal axial flow and an external cross flow.The interactions between the two pipes are modeled by a linear distributed damper,a linear distributed spring and a nonlinear distributed spring along the pipe length.The unsteady hydrodynamic forces due to cross flow are modeled by two distributed van der Pol wake oscillators.The nonlinear partial differential equations for the two pipes and the wake are further discretized by the aid of Galerkin’s technique,resulting in a set of ordinary differential equations.These ordinary differential equations are further numeri cally solved by using a fourth-order Runge-Kutta integration algorithm.Phase portraits,bifurcation diagrams,an Argand diagram and oscillation shape diagrams are plotted,showing the existence of a lock-in phenomenon and figure-of-eight trajectory.The PIP system subjected to cross flow displays some interesting dynamical behaviors different from that of a single-pipe structure.展开更多
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.展开更多
The three-dimensional internal flow field of centrifugal pump is complex and variable with design parameters and operation conditions. The post-processing technique named differential amplification method was proposed...The three-dimensional internal flow field of centrifugal pump is complex and variable with design parameters and operation conditions. The post-processing technique named differential amplification method was proposed for the comparison study of different flow structures. The full steady flow fields of an industrial centrifugal pump working on-design and off-design points were numerically investigated by solving Reynolds average Navier-Stokes equations together with a shear-stress transport(SST) k-? turbulence model. And the numerically predicted performance curves of the studied pump agree well with test measurement results. Compared with the flow flied on design point under the help of differential amplification method, the disturbance caused by interaction between blade and volute tongue is very obvious and it extends to the diffuser pipe on the working point with 0.8 times rated flux. While on the point with 1.2 times rated flux, the flow distribution in impeller region is roughly even and it flows more to the bottom section of the diffuser pipe. The above method was proved to be good at displaying the subtle secondary flow structure changes with a higher resolution effect relative to single isolated case observation, which helps the optimization decision-making from multiple design cases.展开更多
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.展开更多
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.展开更多
Characterizing the complex two-phase hydrodynamics in structured packed columns requires a power- ful modeling tool. The traditional two-dimensional model exhibits limitations when one attempts to model the de- tailed...Characterizing the complex two-phase hydrodynamics in structured packed columns requires a power- ful modeling tool. The traditional two-dimensional model exhibits limitations when one attempts to model the de- tailed two-phase flow inside the columns. The present paper presents a three-dimensional computational fluid dy- namics (CFD) model to simulate the two-phase flow in a representative unit of the column. The unit consists of an CFD calculations on column packed with Flexipak 1Y were implemented within the volume of fluid (VOF) mathe- matical framework. The CFD model was validated by comparing the calculated thickness of liquid film with the available experimental data. Special attention was given to quantitative analysis of the effects of gravity on the hy- drodynamics. Fluctuations in the liquid mass flow rate and the calculated pressure drop loss were found to be quali- tatively in agreement with the experimental observations.展开更多
Based on patent cooperation data,this study used a range of city network analysis approaches in order to explore the structure of the Chinese city network which is driven by technological knowledge flows.The results r...Based on patent cooperation data,this study used a range of city network analysis approaches in order to explore the structure of the Chinese city network which is driven by technological knowledge flows.The results revealed the spatial structure,composition structure,hierarchical structure,group structure,and control structure of Chinese city network,as well as its dynamic factors.The major findings are:1) the spatial pattern presents a diamond structure,in which Wuhan is the central city;2) although the invention patent knowledge network is the main part of the broader inter-city innovative cooperation network,it is weaker than the utility model patent;3) as the senior level cities,Beijing,Shanghai and the cities in the Zhujiang(Pearl) River Delta Region show a strong capability of both spreading and controlling technological knowledge;4) whilst a national technology alliance has preliminarily formed,regional alliances have not been adequately established;5) even though the cooperation level amongst weak connection cities is not high,such cities still play an important role in the network as a result of their location within ′structural holes′ in the network;and 6) the major driving forces facilitating inter-city technological cooperation are geographical proximity,hierarchical proximity and technological proximity.展开更多
A new type of drainage channel with an energy dissipation structure has been proposed based on previous engineering experiences and practical requirements for hazard mitigation in earthquakeaffected areas.Experimental...A new type of drainage channel with an energy dissipation structure has been proposed based on previous engineering experiences and practical requirements for hazard mitigation in earthquakeaffected areas.Experimental studies were performed to determine the characteristics of viscous debris flow in a drainage channel of this type with a slope of 15%.The velocity and depth of the viscous debris flow were measured,processed,and subsequently used to characterize the viscous debris flow in the drainage channel.Observations of this experiment showed that the surface of the viscous debris flow in a smooth drainage channel was smoother than that of a similar debris flow passing through the energy dissipation section in a channel of the new type studied here.However,the flow patterns in the two types of channels were similar at other points.These experimental results show that the depth of the viscous debris flow downstream of the energy dissipation structure increased gradually with the length of the energy dissipation structure.In addition,in the smooth channel,the viscous debris-flow velocity downstream of the energy dissipation structure decreased gradually with the length of the energy dissipation structure.Furthermore,theviscous debris-flow depth and velocity were slightly affected by variations in the width of the energy dissipation structure when the channel slope was 15%.Finally,the energy dissipation ratio increased gradually as the length and width of the energy dissipation structure increased;the maximum energy dissipation ratio observed was 62.9%(where B = 0.6m and L/w = 6.0).展开更多
The discussion on sedimentary macro-structures and their forming mechanics of debris flow is based on thedata of present processes and ancient deposits of Dongchuan in Yunnan Province, Wudu in Gansu Province and Fuxin...The discussion on sedimentary macro-structures and their forming mechanics of debris flow is based on thedata of present processes and ancient deposits of Dongchuan in Yunnan Province, Wudu in Gansu Province and Fuxin inLiaoning Province. Non-cohesive debris flow, which is 1. 3 -1. 7 t/m3 in density, follows hydraulic fluid and flow model. In the hydraulic fluid of flood, electrolytic water combines clay into pulp to transport solid debris. The sedimentary structures show fluid processes as stone-line structure, imbricated structure and stone supporting-sustructure. A part of non-cohesive debris flow with a density of 1. 7 -1. 9t/m3 follows granular model. Debris is in the action of friction, collision and dispersion which forces as debris moving upward to form reverse graded load. The sedimentary structures of granular flow show reverse graded bedding, reverse-normal graded bedding, imbricated verticalstructure and circling linear structure.Cohesive debris flow, which is 1. 9 - 2. 3 t/m3 in density, follows structural two phase flow (viscoplastic) model. The sedimentary structures of cohesive debris flow show reverse graded-chaotic structure,bottom mud-chaotic structure and outwedging structure.展开更多
Field geological work, field engineering monitoring, laboratory experiments and numerical simulation were used to study the development characteristics of pore-fracture system and hydraulic fracture of No.3 coal reser...Field geological work, field engineering monitoring, laboratory experiments and numerical simulation were used to study the development characteristics of pore-fracture system and hydraulic fracture of No.3 coal reservoir in Southern Qinshui Basin. Flow patterns of methane and water in pore-fracture system and hydraulic fracture were discussed by using limit method and average method. Based on the structure model and flow pattern of post-fracturing high-rank coal reservoir, flow patterns of methane and water were established. Results show that seepage pattern of methane in pore-fracture system is linked with pore diameter, fracture width, coal bed pressure and flow velocity. While in hydraulic fracture, it is controlled by fracture height, pressure and flow velocity. Seepage pattern of water in pore-fracture system is linked with pore diameter, fracture width and flow velocity. While in hydraulic fracture, it is controlled by fracture height and flow velocity. Pores and fractures in different sizes are linked up by ultramicroscopic fissures, micro-fissures and hydraulic fracture. In post-fracturing high-rank coal reservoir, methane has level-three flow and gets through triple medium to the wellbore; and water passes mainly through double medium to the wellbore which is level-two flow.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.42027806 and 42041006)。
文摘The size of pores or the grille spacing of water–sediment separation structures directly affects their regulation effect on the debris flow performance.A suitable pore size or grille spacing can effectively improve the water–sediment separation ability of the structure.The new funnel-type grating water–sediment separation structure(FGWSS)combines vertical and horizontal structures and provides a satisfactory water–sediment separation effect.However,the regulation effect of the grille spacing of the structure on the debris flow performance has not been studied.The regulation effect of the structure grille spacing on the debris flow performance is studied through a flume test,and the optimal structure grille spacing is obtained.An empirical equation of the relationship between the relative grille spacing of the structure and the sediment separation rate is established.Finally,the influence of the water–sediment separation structure on the regulation effect of debris flows is examined from two aspects:external factors(properties of debris flows)and internal factors(structural factors).The experimental results show that the gradation characteristics of solid particles in debris flows constitute a key factor affecting the regulation effect of the structure on the debris flow performance.The optimum grille spacing of the FGWSS matches the particle size corresponding to the material distribution curves d85~d90 of the debris flow.The total separation rate of debris flow particles is related to the grille spacing of the structure and the content of coarse and fine particles in the debris flow.
基金supported by National Natural Science Foundation of China(Grant No.42172159)Science Foundation of China University of Petroleum,Beijing(Grant No.2462023XKBH002).
文摘Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.
文摘An experimental system for forming a rotating paraboloid shaped shallow water with a free surface was conducted to study coherent vortex structures in a differentially rotating quasi two dimensional zonal flow.Flow visualization and laser light scattering techniques were used to obtain the information of spatial flow patterns.Experimental results show that the coexistence of Coriolis effect and strong shear in latitudinal zones may lead to formation of coherent vortices.Power spectra analysis and photographs which were taken in a reference frame rotating with the observed vortices also justified the emergence,drift and evolution of persistent vortices on the large scale.Locked vortex state manifests the cyclone and anticyclone asymmetry.
文摘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(Grant No. 50579007)the Doctoral Program of Higher Education of China(Grant No. 200802170010)
文摘The theoretical research on the propulsive principle of aquatic animal becomes more important and attracted more researchers to make efforts on it. In the present study, a computational fluid dynamic (CFD) simulation of a three-dimensional traveling-wave undulations body of tuna has been developed to investigate the fluid flow features and vorticity structures around this body when moving in a straight line. The undulation only takes place in the posterior half of the fish, and the tuna-tail is considered as a lunate fin oscillating with the mode combined swaying with yawing. A Reynolds-averaged Navier-Stokes (RANS) equation is developed, employing a control-volume method and a k-omega SST turbulent model; meanwhile an unstructured tetrahedral grid, which is generated for the three-dimensional geometry, is used based on the deformation of the hind parts of the body and corresponding movement of the tail. We calculated the hydrodynamic performance of tuna-like body when a tuna swims in a uniform velocity, and compared the input power coefficient, output power coefficient and propulsive efficiency of the oscillating tuna-tail with or without body vortex shedding. Additionally, the load distribution on the body, flow features and vorticity structures around the body were demonstrated. The effect of interaction between the body-generated vortices and the tail-generated vorticity on the hydrodynamic performance can be obtained.
基金Project(U1134203)supported by the National Natural Science Foundation of China
文摘According to the analysis of the turbulent intensity level around the high-speed train, the maximum turbulent intensity ranges from 0.2 to 0.5 which belongs to high turbulent flow. The flow field distribution law was studied and eight types of flow regions were proposed. They are high pressure with air stagnant region, pressure decreasing with air accelerating region, low pressure with high air flow velocity region I, turbulent region, steady flow region, low pressure with high air flow velocity region II,pressure increasing with air decelerating region and wake region. The analysis of the vortex structure around the train shows that the vortex is mainly induced by structures with complex mutation and large curvature change. The head and rear of train, the underbody structure, the carriage connection section and the wake region are the main vortex generating sources while the train body with even cross-section has rare vortexes. The wake structure development law studied lays foundation for the train drag reduction.
基金supported by China Armament Pre-research Foundation(Grant No. 51318010402)UK Engineering and Physical Science Research Council (EPSRC), and China Scholarship Council (Grant No.2010611054)
文摘Work on dynamic topology optimization of engineering structures for vibration suppression has mainly addressed the maximization of eigenfrequencies and gaps between consecutive eigenfrequencies of free vibration, minimization of the dynamic compliance subject to forced vibration, and minimization of the structural frequency response. A dynamic topology optimization method of bi-material plate structures is presented based on power flow analysis. Topology optimization problems formulated directly with the design objective of minimizing the power flow response are dealt with. In comparison to the displacement or velocity response, the power flow response takes not only the amplitude of force and velocity into account, but also the phase relationship of the two vector quantities. The complex expression of power flow response is derived based on time-harmonic external mechanical loading and Rayleigh damping. The mathematical formulation of topology optimization is established based on power flow response and bi-material solid isotropic material with penalization(SIMP) model. Computational optimization procedure is developed by using adjoint design sensitivity analysis and the method of moving asymptotes(MMA). Several numerical examples are presented for bi-material plate structures with different loading frequencies, which verify the feasibility and effectiveness of this method. Additionally, optimum results between topological design of minimum power flow response and minimum dynamic compliance are compared, showing that the present method has strong adaptability for structural dynamic topology optimization problems. The proposed research provides a more accurate and effective approach for dynamic topology optimization of vibrating structures.
基金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.
基金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.
基金The work was supported by the National Natural Science Foundation of China(Grant 11622216).
文摘Pipe-in-pipe(PIP)structures are widely used in offshore oil and gas pipelines to settle thermal insulation issues.A PIP structure system usually consists of two concentric pipes and one softer layer for thermal insulation consideration.The total response of the system is related to the dynamics of both pipes and the interactions between these two concentric pipes.In the current work,a theoretical model for flow-induced vibrations of a PIP structure system is proposed and analyzed in the presence of an internal axial flow and an external cross flow.The interactions between the two pipes are modeled by a linear distributed damper,a linear distributed spring and a nonlinear distributed spring along the pipe length.The unsteady hydrodynamic forces due to cross flow are modeled by two distributed van der Pol wake oscillators.The nonlinear partial differential equations for the two pipes and the wake are further discretized by the aid of Galerkin’s technique,resulting in a set of ordinary differential equations.These ordinary differential equations are further numeri cally solved by using a fourth-order Runge-Kutta integration algorithm.Phase portraits,bifurcation diagrams,an Argand diagram and oscillation shape diagrams are plotted,showing the existence of a lock-in phenomenon and figure-of-eight trajectory.The PIP system subjected to cross flow displays some interesting dynamical behaviors different from that of a single-pipe structure.
基金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.
基金Project(2014GK3150)supported by Science and Technology Plan of Hunan Province,China
文摘The three-dimensional internal flow field of centrifugal pump is complex and variable with design parameters and operation conditions. The post-processing technique named differential amplification method was proposed for the comparison study of different flow structures. The full steady flow fields of an industrial centrifugal pump working on-design and off-design points were numerically investigated by solving Reynolds average Navier-Stokes equations together with a shear-stress transport(SST) k-? turbulence model. And the numerically predicted performance curves of the studied pump agree well with test measurement results. Compared with the flow flied on design point under the help of differential amplification method, the disturbance caused by interaction between blade and volute tongue is very obvious and it extends to the diffuser pipe on the working point with 0.8 times rated flux. While on the point with 1.2 times rated flux, the flow distribution in impeller region is roughly even and it flows more to the bottom section of the diffuser pipe. The above method was proved to be good at displaying the subtle secondary flow structure changes with a higher resolution effect relative to single isolated case observation, which helps the optimization decision-making from multiple design cases.
基金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.
基金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 Major State Basic Research Development Program of China(2011CB706501)the National Natural Science Foundation of China(51276157)
文摘Characterizing the complex two-phase hydrodynamics in structured packed columns requires a power- ful modeling tool. The traditional two-dimensional model exhibits limitations when one attempts to model the de- tailed two-phase flow inside the columns. The present paper presents a three-dimensional computational fluid dy- namics (CFD) model to simulate the two-phase flow in a representative unit of the column. The unit consists of an CFD calculations on column packed with Flexipak 1Y were implemented within the volume of fluid (VOF) mathe- matical framework. The CFD model was validated by comparing the calculated thickness of liquid film with the available experimental data. Special attention was given to quantitative analysis of the effects of gravity on the hy- drodynamics. Fluctuations in the liquid mass flow rate and the calculated pressure drop loss were found to be quali- tatively in agreement with the experimental observations.
基金Under the auspices of Major Project of National Social Science Foundation of China(No.13&ZD027)National Natural Science Foundation of China(No.41201128,71433008)
文摘Based on patent cooperation data,this study used a range of city network analysis approaches in order to explore the structure of the Chinese city network which is driven by technological knowledge flows.The results revealed the spatial structure,composition structure,hierarchical structure,group structure,and control structure of Chinese city network,as well as its dynamic factors.The major findings are:1) the spatial pattern presents a diamond structure,in which Wuhan is the central city;2) although the invention patent knowledge network is the main part of the broader inter-city innovative cooperation network,it is weaker than the utility model patent;3) as the senior level cities,Beijing,Shanghai and the cities in the Zhujiang(Pearl) River Delta Region show a strong capability of both spreading and controlling technological knowledge;4) whilst a national technology alliance has preliminarily formed,regional alliances have not been adequately established;5) even though the cooperation level amongst weak connection cities is not high,such cities still play an important role in the network as a result of their location within ′structural holes′ in the network;and 6) the major driving forces facilitating inter-city technological cooperation are geographical proximity,hierarchical proximity and technological proximity.
基金supported by the Key Deployment Project of Chinese Academy of Sciences (Grant No.KZZD-EW-05-01)the National Natural Science Foundation of China (Grant No.41302283)the West Light Foundation of Chinese Academy of Sciences
文摘A new type of drainage channel with an energy dissipation structure has been proposed based on previous engineering experiences and practical requirements for hazard mitigation in earthquakeaffected areas.Experimental studies were performed to determine the characteristics of viscous debris flow in a drainage channel of this type with a slope of 15%.The velocity and depth of the viscous debris flow were measured,processed,and subsequently used to characterize the viscous debris flow in the drainage channel.Observations of this experiment showed that the surface of the viscous debris flow in a smooth drainage channel was smoother than that of a similar debris flow passing through the energy dissipation section in a channel of the new type studied here.However,the flow patterns in the two types of channels were similar at other points.These experimental results show that the depth of the viscous debris flow downstream of the energy dissipation structure increased gradually with the length of the energy dissipation structure.In addition,in the smooth channel,the viscous debris-flow velocity downstream of the energy dissipation structure decreased gradually with the length of the energy dissipation structure.Furthermore,theviscous debris-flow depth and velocity were slightly affected by variations in the width of the energy dissipation structure when the channel slope was 15%.Finally,the energy dissipation ratio increased gradually as the length and width of the energy dissipation structure increased;the maximum energy dissipation ratio observed was 62.9%(where B = 0.6m and L/w = 6.0).
文摘The discussion on sedimentary macro-structures and their forming mechanics of debris flow is based on thedata of present processes and ancient deposits of Dongchuan in Yunnan Province, Wudu in Gansu Province and Fuxin inLiaoning Province. Non-cohesive debris flow, which is 1. 3 -1. 7 t/m3 in density, follows hydraulic fluid and flow model. In the hydraulic fluid of flood, electrolytic water combines clay into pulp to transport solid debris. The sedimentary structures show fluid processes as stone-line structure, imbricated structure and stone supporting-sustructure. A part of non-cohesive debris flow with a density of 1. 7 -1. 9t/m3 follows granular model. Debris is in the action of friction, collision and dispersion which forces as debris moving upward to form reverse graded load. The sedimentary structures of granular flow show reverse graded bedding, reverse-normal graded bedding, imbricated verticalstructure and circling linear structure.Cohesive debris flow, which is 1. 9 - 2. 3 t/m3 in density, follows structural two phase flow (viscoplastic) model. The sedimentary structures of cohesive debris flow show reverse graded-chaotic structure,bottom mud-chaotic structure and outwedging structure.
基金Projects(41330638,41272154)supported by the National Natural Science Foundation of ChinaProject supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD),ChinaProject(2014M551705)supported by the China Postdoctoral Science Foundation
文摘Field geological work, field engineering monitoring, laboratory experiments and numerical simulation were used to study the development characteristics of pore-fracture system and hydraulic fracture of No.3 coal reservoir in Southern Qinshui Basin. Flow patterns of methane and water in pore-fracture system and hydraulic fracture were discussed by using limit method and average method. Based on the structure model and flow pattern of post-fracturing high-rank coal reservoir, flow patterns of methane and water were established. Results show that seepage pattern of methane in pore-fracture system is linked with pore diameter, fracture width, coal bed pressure and flow velocity. While in hydraulic fracture, it is controlled by fracture height, pressure and flow velocity. Seepage pattern of water in pore-fracture system is linked with pore diameter, fracture width and flow velocity. While in hydraulic fracture, it is controlled by fracture height and flow velocity. Pores and fractures in different sizes are linked up by ultramicroscopic fissures, micro-fissures and hydraulic fracture. In post-fracturing high-rank coal reservoir, methane has level-three flow and gets through triple medium to the wellbore; and water passes mainly through double medium to the wellbore which is level-two flow.