Regulation effect of the grille spacing of a funnel-type grating water–sediment separation structure on the debris flow performance

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.

Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties

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.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Temporal change of plankton size structure preserved by Lugol's solution:a FlowCAM study

Plankton size structure is crucial for understanding marine ecosystem dynamics and the associated biogeochemical processes.A fixation step by acid Lugol’s solution has been commonly employed to preserve plankton samples in the field.However,the acid Lugol’s solution can bias the estimation of size structure and the preserved plankton size structure can vary with time.Here,we explore the impact of sample storage time on the size-structure of the plankton community preserved by Lugol’s solution.Two short-term experiments and one long-term experiment were conducted to explore the change of plankton community size structure with the storage time:covering from a week to a month,and to nearly seven months based on particle-size data obtained by continuous Flow Cytometer and Microscope(FlowCAM)measurements.We found a linear change of plankton size with the storage time in short-term periods(less than 3 months)with a decrease of the slope but an increase of the intercept for the normalized biomass size spectrum(NBS S).However,there were opposite trends for NBSS with increasing slope but decreasing intercept after3 months.The potential causes of the distinct patterns of the NBSS parameters are addressed in terms of the interplay between particle aggregation and fragmentation.We found large changes in plankton biovolume and abundance among different size classes,which may indicate a distinct effect of acid Lugol’s solution on various plankton size classes.The mechanism driving temporal change in the size-structure of the Lugolfixed plankton community was further discussed in terms of particle aggregation and fragmentation.Finally,we emphasize that the effect of storage time should be taken into account when interpreting or comparing data of plankton community acquired from samples with various storage durations.

Flow Structure and Short-Term Riverbed Evolution in Curved Flumes

River bending is the major effect responsible for bed topography and bank changes.In this study,fluid velocity(measured by a three-dimensional Doppler advanced point current meter)and bed topographical data have been collected in 40 sections of an experimental model.The whole flume was composed of an organic glass bend,upstream and downstream water tanks,two transition straight sections,a circulation pump,and a connection pipeline.Each section has been found to be characterized by a primary circulation and a small reverse circulation,with some sections even presenting three more or more circulation structures.The minimum circulation intensity has been detected in proximity to the top of the curved channel,while a region with small longitudinal velocity has been observed near the concave bank of each bend,corresponding to the flat bed formed after a short period of scouring.The maximum sediment deposition and scour depth in the presence of a uniform distribution of living flexible vegetation within 10 cm of the flume wall have been found to be smaller than those observed in the tests conducted without vegetation.

Pore-scale modeling of pore structure properties and wettability effect on permeability of low-rank coal

Permeability is a key parameter for coalbed methane development.Although the absolute permeability of coal has been extensively studied,wettability and pore structure properties continue to challenge the microscopic description of water-gas flow in coal.For this purpose,we reconstructed the microstructures of low-rank coal using micro-computed tomography(micro-CT)images.Pore geometry and pore-throat parameters are introduced to establish a relationship with absolute permeability.A dual-porosity pore network model is developed to study water-gas displacement under different wetting and pore structure properties.Results show that absolute permeability is significantly affected by pore geometry and can be described using a binary quadratic function of porosity and fractal dimension.Water-gas relative permeability varies significantly and the residual gas saturation is lower;the crossover saturation first decreased and then increased with increasing porosity under hydrophobic conditions.While the water relative permeability is lower and a certain amount of gas is trapped in complex pore-throat networks;the crossover saturation is higher under hydrophilic conditions.Models with large percolating porosity and well-developed pore networks have high displacement efficiency due to low capillary resistance and avoidance of trapping.This work provides a systematic description of absolute permeability and water-gas relative permeability in coal microstructure for enhanced gas recovery.

Inlet Recirculation Influence to the Flow Structure of Centrifugal Impeller

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.

CFD Simulation of Flow Features and Vorticity Structures in Tuna-Like Swimming

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.

Instantaneous and time-averaged flow structures around a blunt double-cone with or without supersonic film cooling visualized via nano-tracer planar laser scattering

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.

Nonplanar flow-induced vibrations of a cantilevered PIP structure system concurrently subjected to internal and cross flows

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.

Power Flow Response Based Dynamic Topology Optimization of Bi-material Plate Structures

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.

Differential amplification method for flow structures analysis of centrifugal pump between design and off-design points

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.

Lagrangian-based investigation of the transient flow structures around a pitching hydrofoil

The objective of this paper is to address the transient flow structures around a pitching hydrofoil by com- bining physical and numerical studies. In order to predict the dynamic behavior of the flow structure effectively, the Lagrangian coherent structures （LCS） defined by the ridges of the finite-time Lyapunov exponent （FTLE） are utilized under the framework of Navier-Stokes flow computations. In the numerical simulations, the k-w shear stress trans- port （SST） turbulence model, coupled with a two-equation F-Reo transition model, is used for the turbulence closure. Results are presented for a NACA66 hydrofoil undergoing slowly and rapidly pitching motions from 0° to 15° then back to 0° at a moderate Reynolds number Re = 7.5 × 105. The results reveal that the transient flow structures can be observed by the LCS method. For the slowly pitching case, it consists of five stages： quasi-steady and laminar, transition from laminar to turbulent, vortex development, large-scale vortex shedding, and reverting to laminar. The observation of LCS and Lagrangian particle tracers elucidates that the trailing edge vortex is nearly attached and stable during the vortex development stage and the interaction between the leading and trailing edge vortex caused by the adverse pres- sure gradient forces the vortexes to shed downstream during the large-scale vortex shedding stage, which corresponds to obvious fluctuations of the hydrodynamic response. For the rapidly pitching case, the inflection is hardly to be observed and the stall is delayed. The vortex formation, interaction, and shedding occurred once instead of being repeated three times, which is responsible for just one fluctuation in the hydrody- namic characteristics. The numerical results also show that the FTLE field has the potential to identify the transient flows, and the LCS can represent the divergence extent of infinite neighboring particles and capture the interface of the vortex region.

A mass balanced model of trophic structure and energy flows of a semi-closed marine ecosystem

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.

SEDIMENTARY MACRO-STRUCTURES AND FORMING MECHANISM OF DEBRIS FLOW

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.

Vertical Structure of Interannual Variability in Cross-Equatorial Flows over the Maritime Continent and Indian Ocean in Boreal Summer

This study investigates the vertical structure of variability in the cross-equatorial flows(CEFs)over the Maritime Continent and Indian Ocean in boreal summer,based on three reanalysis datasets:ERA-Interim,JRA-55 and NCEP-2.The results show a high consistency in the interannual variability among the reanalysis datasets,especially between ERAInterim and JRA-55,while great uncertainty exists in the decadal or long-term changes of CEFs.Further analyses on the interannual variability in CEFs indicate that there is a significant negative relationship between the upper-and lower-level CEFs over the Maritime Continent—that is,the northerlies at the upper level and southerlies at the lower level are both enhanced or weakened.This seesaw pattern is also significantly related to the CEFs over the Indian Ocean at the upper level and lower level(i.e.,the Somali jet).This close relationship between the upper-and lower-level CEFs and between the Maritime Continent and Indian Ocean is manifested as the leading mode of equatorial meridional winds in the vertical-zonal section over the Maritime Continent and Indian Ocean.Finally,it is found that ENSO is closely related to the vertical leading mode,and to all the CEFs at the upper and lower levels over the Maritime Continent and Indian Ocean.

Fluid−Structure Interaction of Two-Phase Flow Passing Through 90° Pipe Bend Under Slug Pattern Conditions

Numerical simulations of evolution characteristics of slug flow across a 90°pipe bend have been carried out to study the fluid−structure interaction response induced by internal slug flow.The two-phase flow patterns and turbulence were modelled by using the volume of fluid(VOF)model and the Realizable k−εturbulence model respectively.Firstly,validation of the CFD model was carried out and the desirable results were obtained.The different flow patterns and the time-average mean void fraction was coincident with the reported experimental data.Simulations of different cases of slug flow have been carried out to show the effects of superficial gas and liquid velocity on the evolution characteristics of slug flow.Then,a one-way coupled fluid-structure interaction framework was established to investigate the slug flow interaction with a 90°pipe bend under various superficial liquid and gas velocities.It was found that the maximum total deformation and equivalent stress increased with the increasing superficial gas velocity,while decreased with the increasing superficial liquid velocity.In addition,the total deformation and equivalent stress has obvious periodic fluctuation.Furthermore,the distribution position of maximum deformation and stress was related to the evolution of slug flow.With the increasing superficial gas velocity,the maximum total deformation was mainly located at the 90°pipe bend.But as the superficial liquid velocity increases,the maximum total deformation was mainly located in the horizontal pipe section.Consequently,the slug flow with higher superficial gas velocity will induce more serious cyclical impact on the 90°pipe bend.

RESEARCH ON ABRASION OF DEBRIS FLOW TO HIGH-SPEED DRAINAGE STRUCTURE

As one weak topic in research of debris flow,abrasion of debris flow shortens obviously application life of control structure composed of concrete.High_speed drainage structure,one of the most effective techniques to control giant debris flow disaster,has shortened one_third application life due to abrasion by debris flow.Based on velocity calculation method founded by two_phase theory,research of abrasion mechanism of debris flow to high_speed drainage structure was made.The mechanism includes both abrasion mechanism of homogeneous sizing and shearing mechanism of particle of debris flow to high_speed drainage trough structure.Further abrasion equations of both sizing and particle were established by Newton movement theory of debris flow.And abrasion amount formula of the high_speed drainage trough structure is set up by dimensional analysis.Amount to calculating in the formula is consistent with testing data in_situ,which is valuable in design of high_speed drainage structure.

Three-dimensional spatial structure of the macro-pores and flow simulation in anthracite coal based on X-ray μ-CT scanning data

The three-dimensional(3 D) structures of pores directly affect the CH4 flow.Therefore,it is very important to analyze the3 D spatial structure of pores and to simulate the CH4 flow with the connected pores as the carrier.The result shows that the equivalent radius of pores and throats are 1-16 μm and 1.03-8.9 μm,respectively,and the throat length is 3.28-231.25 μm.The coordination number of pores concentrates around three,and the intersection point between the connectivity function and the X-axis is 3-4 μm,which indicate the macro-pores have good connectivity.During the single-channel flow,the pressure decreases along the direction of CH4 flow,and the flow velocity of CH4 decreases from the pore center to the wall.Under the dual-channel and the multi-channel flows,the pressure also decreases along the CH4 flow direction,while the velocity increases.The mean flow pressure gradually decreases with the increase of the distance from the inlet slice.The change of mean flow pressure is relatively stable in the direction horizontal to the bedding plane,while it is relatively large in the direction perpendicular to the bedding plane.The mean flow velocity in the direction horizontal to the bedding plane(Y-axis) is the largest,followed by that in the direction horizontal to the bedding plane(X-axis),and the mean flow velocity in the direction perpendicular to the bedding plane is the smallest.