This paper presents a numerical study of the flow topologies of three-dimensional （3D） flows in a high pressure compressor stator blade row without and with boundary layer aspiration on the hub wall. The stator blad...This paper presents a numerical study of the flow topologies of three-dimensional （3D） flows in a high pressure compressor stator blade row without and with boundary layer aspiration on the hub wall. The stator blade is representative of the first stage operating under transonic inlet conditions and the blade design encourages development of highly complex 3D flows. The blade has a small tip clearance. The computational fluid dynamics （CFD） studies show progressive increase of hub corner stall with the increase in incidence. Aspiration is implemented on the hub wall via a slot in the comer between the hub wall and the suction surface. The CFD studies show aspiration to be sensitive to the suction flow rate; lower rate leads to very complex flow struc- tures and increased level of losses whereas higher rate renders aspiration effective for control of hub comer separation. The flow topologies are studied by trace of skin friction lines on the walls. The nature of flow can be explained by the topological rules of closed separation. Furthermore, a deeper analysis is done for a particular case with advanced criterion to test the non-degeneracy of critical points in the flow field.展开更多
The increasing performance demands of modern aero engines necessitate the integrated design of compressor transition ducts with upstream components to reduce the axial length of the engine.However,this design approach...The increasing performance demands of modern aero engines necessitate the integrated design of compressor transition ducts with upstream components to reduce the axial length of the engine.However,this design approach narrows the spacing between the stator and the strut,making traditional research on transition ducts only with struts unsuitable.The numerical results and experimental oil flow visualization results were utilized to reconstruct the three-dimensional flow structures in the stator passages under various operating conditions.Additionally,numerical methods were employed to analyze the mechanisms of the strut's effect on the upstream stator in an aggressive transition duct.The results show that the strut potential field increases the load on the upstream stator,leading to severe blade surface separation and corner separation/stall,and redistributes the inflow angle of the upstream stators circumferentially,resulting in significant differences in the flow structures within the stator passages on both sides.The separation flows within the stator passages mainly manifest in five types:pressure surface separation vortex,suction surface concentrated shedding vortex,suction surface separation vortex,suction surface-corner stall separation vortex,and suction surface separation vortex pair.Under different operating conditions,the separation flows within the stator passages are always composed of a part of these five types or a transitional state between two of them.展开更多
The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet ...The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet dynamics,are demonstrated in the current research.Hereby,the 2nd Stokes problem assumes that the surface,with a thin shared layer of the fluid on it,oscillates in a harmonic manner along the x-axis of the rivulet flow,which coincides with the main flow direction streaming down the underlying surface.We obtain the exact extension of the rivulet flow family,clarifying the structure of the pressure field,which fully absorbs the arising perturbation.The profile of the velocity field is assumed to be Gaussian-type with a non-zero level of plasticity.Hence,the absolutely non-Newtonian case of the viscoplastic flow solution,which satisfies the motion and continuity equations,is considered(with particular cases of exact solutions for pressure).The perturbed governing equations of motion for rivulet flows then result in the Riccati-type ordinary differential equation(ODE),describing the dynamics of the coordinate x(t).The approximated schematic dynamics are presented in graphical plots.展开更多
Particle-and droplet-laden flows are central to many problems in mechanics and transport.They occur in sedimentladen boundary layers,gas-solid and gas-liquid dispersions,and surface-water films driven by external forc...Particle-and droplet-laden flows are central to many problems in mechanics and transport.They occur in sedimentladen boundary layers,gas-solid and gas-liquid dispersions,and surface-water films driven by external forcing.They also underpin practical applications ranging from environmental transport to high-speed and aerothermal systems.Despite decades of progress,prediction remains difficult.The physics spans a wide range of scales and often couples turbulence,interphase momentum exchange,collisions,and interfacial transport.Reliable computation therefore requires both robust numerical methodology and careful physical interpretation.展开更多
According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing ...According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.展开更多
Large eddy simulations were used to model the three-dimensional flows around a non-submerged spur dike. Based on the rigid lid assumption, all these vortex flows around the spur dike, which probably affected the whole...Large eddy simulations were used to model the three-dimensional flows around a non-submerged spur dike. Based on the rigid lid assumption, all these vortex flows around the spur dike, which probably affected the whole flow field, were numerically simulated and analyzed. The largest circulating flow region caused by the spur dike behind it was weak, which would naturally lead to sedimentation for silt-laden twophase flows. Based on the Smagorinsky model, the finite volume method was used to discretize the NavierStokes equations, and the SIMPLEC algorithm was used to solve them. Meanwhile, these flows were investigated experimentally in a 0.5-m wide flume with a dike placed as a barrier. The computational results are in a fairly good agreement with the experimental data.展开更多
A comprehensive assessment of grain supply,demand,and ecosystem service flows is essential for identifying grain movement pathways,ensuring regional grain security,and guiding sustainable management strategies.However...A comprehensive assessment of grain supply,demand,and ecosystem service flows is essential for identifying grain movement pathways,ensuring regional grain security,and guiding sustainable management strategies.However,current studies primarily focus on short-term grain provision services while neglecting the spatiotemporal variations in grain flows across different scales.This gap limits the identification of dynamic matching relationships and the formulation of optimization strategies for balancing grain flows.This study examined the spatiotemporal evolution of grain supply and demand in the Beijing-Tianjin-Hebei(BTH)region from 1980 to 2020.Using the Enhanced TwoStep Floating Catchment Area method,the grain provision ecosystem service flows were quantified,the changes in supply–demand matching under different grain flow scenarios were analyzed and the optimal distance threshold for grain flows was investigated.The results revealed that grain production follows a spatial distribution pattern characterized by high levels in the southeast and low levels in the northwest.A significant mismatch exists between supply and demand,and it shows a scale effect.Deficit areas are mainly concentrated in the northwest,while surplus areas are mainly located in the central and southern regions.As the spatial scale increases,the ecosystem service supply–demand ratio(SDR)classification becomes more clustered,while it exhibits greater spatial SDR heterogeneity at smaller scales.This study examined two distinct scenarios of grain provision ecosystem service flow dynamics based on 100 and 200 km distance thresholds.The flow increased significantly,from 2.17 to 11.81million tons in the first scenario and from 2.41 to 12.37 million tons in the second scenario over nearly 40 years,forming a spatial movement pattern from the central and southern regions to the surrounding areas.Large flows were mainly concentrated in the interior of urban centers,with significant outflows between cities such as Baoding,Shijiazhuang,Xingtai,and Hengshui.At the county scale,supply–demand matching patterns remained consistent between the grain flows in the two scenarios.Notably,incorporating grain flow dynamics significantly reduced the number of grain-deficit areas compared to scenarios without grain flow.In 2020,grain-deficit counties decreased by28.79 and 37.88%,and cities by 12.50 and 25.0%under the two scenarios,respectively.Furthermore,the distance threshold for achieving optimal supply and demand matching at the county scale was longer than at the city scale in both grain flow scenarios.This study provides valuable insights into the dynamic relationships and heterogeneous patterns of grain matching,and expands the research perspective on grain and ecosystem service flows across various spatiotemporal scales.展开更多
Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture...Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture models often struggle to capture the complexities introduced by coarse boulders and multi-phase interactions,while strong-coupling methods can be computationally prohibitive for practical hazard assessments.In this study,we propose a semi-hybrid,fully resolved coupling numerical framework for modeling boulder-laden debris flows.This framework conceptualizes debris flows as a composite system comprising a continuous viscous fluidphase(including finesediments)and a discrete phase of arbitrarily shaped coarse particles.The continuous phase is treated as a generalized nonlinear Coulomb-viscoplastic fluidusing the smoothed particle hydrodynamics(SPH)method,while coarse particles are modeled via the distributed contact discrete element method(DCDEM).These two phases are coupled through an efficienttwo-way resolved scheme,ensuring accurate simulation of flow-boulder interactions within a unifiedtimeframe.We validate the proposed method against two physical experiments:(1)gravity-driven concrete flows and(2)debris flowinteracting with slit-type barriers.Results confirmthe method's robustness in accurately capturing fluid-solid-structureinteractions and deposition processes.Its capabilities are further showcased through the simulation of a stony debris-flowevent inWenchuan County,China,highlighting its promise for real-world engineering applications and validating the effectiveness of the existing cascade dam system in mitigating debrisflowimpact and energy dissipation.展开更多
Hypersonic flow over a canonical 25°−55°double-cone configuration with a freestream Mach number of 10.38 is numerically investigated.In contrast to a relatively stable and moderate separation region in exper...Hypersonic flow over a canonical 25°−55°double-cone configuration with a freestream Mach number of 10.38 is numerically investigated.In contrast to a relatively stable and moderate separation region in experiments,axisymmetric calculations in the literature demonstrated that the separation bubble continuously grew and notably exceeded the experimental result.To explain the causes of the discrepancy,time-accurate axisymmetric and Three-Dimensional(3D)simulations are conducted for the double-cone flow to investigate the evolution of three-dimensionality and unsteadiness.Both the axisymmetric calculation and the 3D simulation without external disturbances predict a significantly larger separation region than that in experiments and misrepresent the distributions of surface pressure and heat flux.The random forcing approach with two levels of noise amplitude is then applied to 3D simulations.A better agreement with the measured data is observed for the time-averaged heat flux and pressure when the white noise is enforced.As the forcing amplitude is increased,the agreement is slightly improved.However,discrepancies between the 3D results and experimental data still exist in the prediction of the heat flux and pressure distributions,indicating the essential difference between the injected white noise and the wind tunnel freestream disturbances.Realistic noise models are required to reveal the sources of such discrepancies.展开更多
Objective:To analyze the diagnostic value of transvaginal three-dimensional ultrasound(3D-TVS)in evaluating endometrial receptivity(ER)for ovulation disorder infertility(ODI),and to investigate the impact of subendome...Objective:To analyze the diagnostic value of transvaginal three-dimensional ultrasound(3D-TVS)in evaluating endometrial receptivity(ER)for ovulation disorder infertility(ODI),and to investigate the impact of subendometrial endometrial vascular index(VI)and endometrial vascular flow index(VFI)levels on ODI.Methods:A total of 110 patients diagnosed with ODI admitted between January 2023 and June 2024 were selected.All patients underwent ovulation induction therapy,3D-TVS examination,and sex hormone testing.Based on pregnancy outcomes,patients were divided into a successful pregnancy group(73 cases)and an unsuccessful pregnancy group(37 cases).ER parameters,sex hormone levels,and endometrial blood flow patterns were compared between the two groups.Receiver operating characteristic(ROC)curves were plotted to evaluate the predictive value of ER for ODI.Results:The spiral artery peak systolic velocity(PSV),endometrial volume(EMV),endometrial flow index(FI),and VFI in the successful pregnancy group were significantly higher than those in the unsuccessful pregnancy group(p<0.05).No significant differences were observed in other ER parameters between the two groups(p>0.05).There was no significant difference in sex hormone levels between the two groups on the day of human chorionic gonadotropin(hCG)treatment(p>0.05).Among the endometrial blood flow classifications in the pregnant group,the proportion of Type II was lower than that in the nonpregnant group(p<0.05).The Receiver Operating Characteristic(ROC)curve demonstrated that the area under the curve(AUC)for Endometrial Volume(EMV)in predicting pregnancy after Ovarian Dysfunction Infertility(ODI)treatment was 0.854,with a sensitivity of 92.61%and a specificity of 71.75%.The AUC for Vascularization Index(VI)was 0.771,with a sensitivity of 52.18%and a specificity of 88.70%.The AUC for Vascularization Flow Index(VFI)of the endometrium was 0.887,with a sensitivity of 80.01%and a specificity of 69.20%.Conclusion:Three-dimensional transvaginal sonography(3D-TVS)assessment of endometrial receptivity(ER)can effectively detect ODI,and the levels of subendometrial VI and VFI demonstrate superior predictive performance for pregnancy outcomes in this condition,serving as commonly used predictive indicators for the disease.展开更多
Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock s...Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock solutions in rectilinear ducts with square cross-sections.In this article,we are devoted to proving rigorously that a large class of these transonic shock solutions are stable,under multidimensional small perturbations of the upcoming supersonic flows and back pressures at the exits of ducts in suitable function spaces.This manifests that frictions have a stabilization effect on transonic shocks in ducts,in consideration of previous works which shown that transonic shocks in purely steady Euler flows are not stable in such ducts.Except its implications to applications,because frictions lead to a stronger coupling between the elliptic and hyperbolic parts of the three-dimensional steady subsonic Euler system,we develop the framework established in previous works to study more complex and interesting Venttsel problems of nonlocal elliptic equations.展开更多
In the past decades,physical modeling has been widely used in hydrogeology for teaching,studying and exhibition purposes.Most of these models are used to illustrate hydrogeological profiles,but few can depict three-di...In the past decades,physical modeling has been widely used in hydrogeology for teaching,studying and exhibition purposes.Most of these models are used to illustrate hydrogeological profiles,but few can depict three-dimensional groundwater flows,making it impossible to validate groundwater flows simulated by numerical methods with physical modeling.展开更多
Long-duration energy storage has become critical for renewable energy integration.While redox flow batteries,especially vanadium-based systems,are scaling up in capacity,their performance at the stack level remains in...Long-duration energy storage has become critical for renewable energy integration.While redox flow batteries,especially vanadium-based systems,are scaling up in capacity,their performance at the stack level remains insufficiently optimized,demanding more profound mechanistic studies and engineering refinements.To address the difficulties in resolving the flow inhomogeneity at the stack scale,this study establishes a multi-physics field coupling model and analyzes the pressure distributions,flow rate differences,active substance concentration,and electrochemical characteristics.The results show that the uneven cell pressure distribution is a key factor affecting the consistency of the system performance,and the increase in the flow rate improves the reactant homogeneity,with both the average concentration and the uniformity factor increasing with the flow rate.In contrast,high current densities lead to an increased imbalance between electrochemical depletion and reactant replenishment,resulting in a significant decrease in reactant concentration in the under-ribs region.In addition,a higher flow rate can expand the high-current-density region where the stack operates efficiently.This study provides a theoretical basis for optimizing the design of the stack components.展开更多
The investigation of Beltrami flows is important for the research on the mechanism of turbulent structure.In this paper the general solutions of the Beltrami flows are given,which depend explicitly on the solutions of...The investigation of Beltrami flows is important for the research on the mechanism of turbulent structure.In this paper the general solutions of the Beltrami flows are given,which depend explicitly on the solutions of three independent Helmholtz equations with scalar unknowns.Velocity fields of Beltrami flows can then be obtained explicitly after the application of some curl operations on the solutions of Helmholtz equa- tions.On the basis of the exact solutions of Euler equations given above,we obtain one kind of exact solu- tions of non-steady Navier-Stokes equations which are also the Beltrami flows.Some interesting examples of Beltrami flows other than'ABC flows','Kolmogolov flows','Rayleigh-Bernard flows','Q-flows' are given.The detailed analytic results of these examples will be published in the near future.展开更多
Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealize...Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealized numerical experiment using the large eddy simulation method was conducted to analyze the three-dimensional flow field structure and the mechanism for dense water subsidence.Specifically,a negative salt flux is set at the sea surface,in which salt flux enters the sea surface to simulate the icing and salting-out phenomena that occur at high latitudes.Results show that the mean-state 3D flow field of dense water subsidence exhibits a hollow conical distribution.The horizontal flow field is characterized by a cyclonic vortex that driven primarily by the pressure gradient and influenced by the Coriolis effect.Moreover,the inverse vertical pressure gradient generated by this vortex inhibits the sinking of the plume,leading to its off-axis deflection and the development of an anticyclonic precession.In addition,the impact of rotation on the structure of a sinking plume within a stratified environment is discussed.Both horizontal vortex intensity and cone angle of the hollow cone flow field are increased with increasing rotation rate,resulting in a decrease in the plume’s maximum sinking depth.Variances in rotation direction cause the horizontal vortex and sinking plumes of dense water in the northern and southern hemispheres to rotate in opposite directions.展开更多
In order to accurately simulate strong three-dimensional （3-D） free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equati...In order to accurately simulate strong three-dimensional （3-D） free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional （2-D） grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.展开更多
Recently,the Liutex-based force field models are proposed and applied to some two-dimensional flows to explore the possibility of direct vortex control strategies.In this paper,we applied the Liutex-based centripetal ...Recently,the Liutex-based force field models are proposed and applied to some two-dimensional flows to explore the possibility of direct vortex control strategies.In this paper,we applied the Liutex-based centripetal force model to three-dimensional flow around a circular cylinder at Re=3900.The turbulent flows around the three-dimensional cylinder are modeled by the delayed detached-eddy simulations based on the two-equation k-ω shear stress transport(SST)model.The cylinder flow without adding the force field model is simulated in the first place.The statistical flow quantities are compared with experimental data to validate the accuracy of the current numerical models.Then the force field model is applied to the momentum equation in a specific control region near the cylinder to study the change of flow field and hydrodynamic performance.The effect of different control region locations with the same force field strength is discussed to find out the most influenced locations on the drag and lift.Based on the screened control region,the influence source term strength is further studied.Variations on the global statistical flow quantities,local instantaneous three-dimensional flow structures are analyzed and discussed.It is concluded that the Liutex-based force field model can effectively change the hydrodynamic behavior of flow past a cylinder and may provide a new direction for performance improvement and optimization.展开更多
In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by mean...In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.展开更多
A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （200...A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.展开更多
A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-co...A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.展开更多
文摘This paper presents a numerical study of the flow topologies of three-dimensional （3D） flows in a high pressure compressor stator blade row without and with boundary layer aspiration on the hub wall. The stator blade is representative of the first stage operating under transonic inlet conditions and the blade design encourages development of highly complex 3D flows. The blade has a small tip clearance. The computational fluid dynamics （CFD） studies show progressive increase of hub corner stall with the increase in incidence. Aspiration is implemented on the hub wall via a slot in the comer between the hub wall and the suction surface. The CFD studies show aspiration to be sensitive to the suction flow rate; lower rate leads to very complex flow struc- tures and increased level of losses whereas higher rate renders aspiration effective for control of hub comer separation. The flow topologies are studied by trace of skin friction lines on the walls. The nature of flow can be explained by the topological rules of closed separation. Furthermore, a deeper analysis is done for a particular case with advanced criterion to test the non-degeneracy of critical points in the flow field.
基金supported by the National Natural Science Foundation of China(No.52276025)the Science Center for Gas Turbine Project of China(Nos.P2022-A-Ⅱ-001-001,P2022-A-Ⅱ-002-001 and P2022-B-Ⅱ-002-001)。
文摘The increasing performance demands of modern aero engines necessitate the integrated design of compressor transition ducts with upstream components to reduce the axial length of the engine.However,this design approach narrows the spacing between the stator and the strut,making traditional research on transition ducts only with struts unsuitable.The numerical results and experimental oil flow visualization results were utilized to reconstruct the three-dimensional flow structures in the stator passages under various operating conditions.Additionally,numerical methods were employed to analyze the mechanisms of the strut's effect on the upstream stator in an aggressive transition duct.The results show that the strut potential field increases the load on the upstream stator,leading to severe blade surface separation and corner separation/stall,and redistributes the inflow angle of the upstream stators circumferentially,resulting in significant differences in the flow structures within the stator passages on both sides.The separation flows within the stator passages mainly manifest in five types:pressure surface separation vortex,suction surface concentrated shedding vortex,suction surface separation vortex,suction surface-corner stall separation vortex,and suction surface separation vortex pair.Under different operating conditions,the separation flows within the stator passages are always composed of a part of these five types or a transitional state between two of them.
文摘The newly formulated non-Newtonian rivulet flows streaming down an inclined planar surface,with additional periodic perturbations arising from the application of the 2nd Stokes problem to the investigation of rivulet dynamics,are demonstrated in the current research.Hereby,the 2nd Stokes problem assumes that the surface,with a thin shared layer of the fluid on it,oscillates in a harmonic manner along the x-axis of the rivulet flow,which coincides with the main flow direction streaming down the underlying surface.We obtain the exact extension of the rivulet flow family,clarifying the structure of the pressure field,which fully absorbs the arising perturbation.The profile of the velocity field is assumed to be Gaussian-type with a non-zero level of plasticity.Hence,the absolutely non-Newtonian case of the viscoplastic flow solution,which satisfies the motion and continuity equations,is considered(with particular cases of exact solutions for pressure).The perturbed governing equations of motion for rivulet flows then result in the Riccati-type ordinary differential equation(ODE),describing the dynamics of the coordinate x(t).The approximated schematic dynamics are presented in graphical plots.
文摘Particle-and droplet-laden flows are central to many problems in mechanics and transport.They occur in sedimentladen boundary layers,gas-solid and gas-liquid dispersions,and surface-water films driven by external forcing.They also underpin practical applications ranging from environmental transport to high-speed and aerothermal systems.Despite decades of progress,prediction remains difficult.The physics spans a wide range of scales and often couples turbulence,interphase momentum exchange,collisions,and interfacial transport.Reliable computation therefore requires both robust numerical methodology and careful physical interpretation.
基金support of the National Natural Science Foundation of China(No.12172023)。
文摘According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.
基金Supported by China Postdoctoral Science Foundation (No.2004036050) and Open Fund of Nanjing Hydraulic Research Institut(No. Yk90504)
文摘Large eddy simulations were used to model the three-dimensional flows around a non-submerged spur dike. Based on the rigid lid assumption, all these vortex flows around the spur dike, which probably affected the whole flow field, were numerically simulated and analyzed. The largest circulating flow region caused by the spur dike behind it was weak, which would naturally lead to sedimentation for silt-laden twophase flows. Based on the Smagorinsky model, the finite volume method was used to discretize the NavierStokes equations, and the SIMPLEC algorithm was used to solve them. Meanwhile, these flows were investigated experimentally in a 0.5-m wide flume with a dike placed as a barrier. The computational results are in a fairly good agreement with the experimental data.
基金supported by the National Natural Science Foundation of China(42471336,52379021 and 42201278)the Hebei Province Backbone Talent Program,China(Returnee Platform for Overseas Study)(A20240028)+2 种基金the Hebei Province Statistical Science Research Project,China(2024HZ04)the Hebei Province Graduate Education and Teaching Reform Research Project,China(YJG2024046)the Innovation Ability Training Program for Postgraduate Students of Hebei Provincial Department of Education,China(CXZZSS2025048)。
文摘A comprehensive assessment of grain supply,demand,and ecosystem service flows is essential for identifying grain movement pathways,ensuring regional grain security,and guiding sustainable management strategies.However,current studies primarily focus on short-term grain provision services while neglecting the spatiotemporal variations in grain flows across different scales.This gap limits the identification of dynamic matching relationships and the formulation of optimization strategies for balancing grain flows.This study examined the spatiotemporal evolution of grain supply and demand in the Beijing-Tianjin-Hebei(BTH)region from 1980 to 2020.Using the Enhanced TwoStep Floating Catchment Area method,the grain provision ecosystem service flows were quantified,the changes in supply–demand matching under different grain flow scenarios were analyzed and the optimal distance threshold for grain flows was investigated.The results revealed that grain production follows a spatial distribution pattern characterized by high levels in the southeast and low levels in the northwest.A significant mismatch exists between supply and demand,and it shows a scale effect.Deficit areas are mainly concentrated in the northwest,while surplus areas are mainly located in the central and southern regions.As the spatial scale increases,the ecosystem service supply–demand ratio(SDR)classification becomes more clustered,while it exhibits greater spatial SDR heterogeneity at smaller scales.This study examined two distinct scenarios of grain provision ecosystem service flow dynamics based on 100 and 200 km distance thresholds.The flow increased significantly,from 2.17 to 11.81million tons in the first scenario and from 2.41 to 12.37 million tons in the second scenario over nearly 40 years,forming a spatial movement pattern from the central and southern regions to the surrounding areas.Large flows were mainly concentrated in the interior of urban centers,with significant outflows between cities such as Baoding,Shijiazhuang,Xingtai,and Hengshui.At the county scale,supply–demand matching patterns remained consistent between the grain flows in the two scenarios.Notably,incorporating grain flow dynamics significantly reduced the number of grain-deficit areas compared to scenarios without grain flow.In 2020,grain-deficit counties decreased by28.79 and 37.88%,and cities by 12.50 and 25.0%under the two scenarios,respectively.Furthermore,the distance threshold for achieving optimal supply and demand matching at the county scale was longer than at the city scale in both grain flow scenarios.This study provides valuable insights into the dynamic relationships and heterogeneous patterns of grain matching,and expands the research perspective on grain and ecosystem service flows across various spatiotemporal scales.
基金supported by the Japan Society for the Promotion of Science(JSPS)KAKENHI(Grant Nos.JP23KK0182,JP23K26356,and JP24K00971).
文摘Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture models often struggle to capture the complexities introduced by coarse boulders and multi-phase interactions,while strong-coupling methods can be computationally prohibitive for practical hazard assessments.In this study,we propose a semi-hybrid,fully resolved coupling numerical framework for modeling boulder-laden debris flows.This framework conceptualizes debris flows as a composite system comprising a continuous viscous fluidphase(including finesediments)and a discrete phase of arbitrarily shaped coarse particles.The continuous phase is treated as a generalized nonlinear Coulomb-viscoplastic fluidusing the smoothed particle hydrodynamics(SPH)method,while coarse particles are modeled via the distributed contact discrete element method(DCDEM).These two phases are coupled through an efficienttwo-way resolved scheme,ensuring accurate simulation of flow-boulder interactions within a unifiedtimeframe.We validate the proposed method against two physical experiments:(1)gravity-driven concrete flows and(2)debris flowinteracting with slit-type barriers.Results confirmthe method's robustness in accurately capturing fluid-solid-structureinteractions and deposition processes.Its capabilities are further showcased through the simulation of a stony debris-flowevent inWenchuan County,China,highlighting its promise for real-world engineering applications and validating the effectiveness of the existing cascade dam system in mitigating debrisflowimpact and energy dissipation.
基金supported by the Hong Kong Research Grants Council,China(Nos.15206519,15204322 and 25203721)the National Natural Science Foundation of China(No.12102377)。
文摘Hypersonic flow over a canonical 25°−55°double-cone configuration with a freestream Mach number of 10.38 is numerically investigated.In contrast to a relatively stable and moderate separation region in experiments,axisymmetric calculations in the literature demonstrated that the separation bubble continuously grew and notably exceeded the experimental result.To explain the causes of the discrepancy,time-accurate axisymmetric and Three-Dimensional(3D)simulations are conducted for the double-cone flow to investigate the evolution of three-dimensionality and unsteadiness.Both the axisymmetric calculation and the 3D simulation without external disturbances predict a significantly larger separation region than that in experiments and misrepresent the distributions of surface pressure and heat flux.The random forcing approach with two levels of noise amplitude is then applied to 3D simulations.A better agreement with the measured data is observed for the time-averaged heat flux and pressure when the white noise is enforced.As the forcing amplitude is increased,the agreement is slightly improved.However,discrepancies between the 3D results and experimental data still exist in the prediction of the heat flux and pressure distributions,indicating the essential difference between the injected white noise and the wind tunnel freestream disturbances.Realistic noise models are required to reveal the sources of such discrepancies.
文摘Objective:To analyze the diagnostic value of transvaginal three-dimensional ultrasound(3D-TVS)in evaluating endometrial receptivity(ER)for ovulation disorder infertility(ODI),and to investigate the impact of subendometrial endometrial vascular index(VI)and endometrial vascular flow index(VFI)levels on ODI.Methods:A total of 110 patients diagnosed with ODI admitted between January 2023 and June 2024 were selected.All patients underwent ovulation induction therapy,3D-TVS examination,and sex hormone testing.Based on pregnancy outcomes,patients were divided into a successful pregnancy group(73 cases)and an unsuccessful pregnancy group(37 cases).ER parameters,sex hormone levels,and endometrial blood flow patterns were compared between the two groups.Receiver operating characteristic(ROC)curves were plotted to evaluate the predictive value of ER for ODI.Results:The spiral artery peak systolic velocity(PSV),endometrial volume(EMV),endometrial flow index(FI),and VFI in the successful pregnancy group were significantly higher than those in the unsuccessful pregnancy group(p<0.05).No significant differences were observed in other ER parameters between the two groups(p>0.05).There was no significant difference in sex hormone levels between the two groups on the day of human chorionic gonadotropin(hCG)treatment(p>0.05).Among the endometrial blood flow classifications in the pregnant group,the proportion of Type II was lower than that in the nonpregnant group(p<0.05).The Receiver Operating Characteristic(ROC)curve demonstrated that the area under the curve(AUC)for Endometrial Volume(EMV)in predicting pregnancy after Ovarian Dysfunction Infertility(ODI)treatment was 0.854,with a sensitivity of 92.61%and a specificity of 71.75%.The AUC for Vascularization Index(VI)was 0.771,with a sensitivity of 52.18%and a specificity of 88.70%.The AUC for Vascularization Flow Index(VFI)of the endometrium was 0.887,with a sensitivity of 80.01%and a specificity of 69.20%.Conclusion:Three-dimensional transvaginal sonography(3D-TVS)assessment of endometrial receptivity(ER)can effectively detect ODI,and the levels of subendometrial VI and VFI demonstrate superior predictive performance for pregnancy outcomes in this condition,serving as commonly used predictive indicators for the disease.
基金This work was supported in part by National Nature Science Foundation of China(11371141 and 11871218)by Science and Technology Commission of Shanghai Municipality(18dz2271000).
文摘Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock solutions in rectilinear ducts with square cross-sections.In this article,we are devoted to proving rigorously that a large class of these transonic shock solutions are stable,under multidimensional small perturbations of the upcoming supersonic flows and back pressures at the exits of ducts in suitable function spaces.This manifests that frictions have a stabilization effect on transonic shocks in ducts,in consideration of previous works which shown that transonic shocks in purely steady Euler flows are not stable in such ducts.Except its implications to applications,because frictions lead to a stronger coupling between the elliptic and hyperbolic parts of the three-dimensional steady subsonic Euler system,we develop the framework established in previous works to study more complex and interesting Venttsel problems of nonlocal elliptic equations.
基金supported by the State Key Program of National Natural Science of China(Grant No.41130637)
文摘In the past decades,physical modeling has been widely used in hydrogeology for teaching,studying and exhibition purposes.Most of these models are used to illustrate hydrogeological profiles,but few can depict three-dimensional groundwater flows,making it impossible to validate groundwater flows simulated by numerical methods with physical modeling.
基金supported by National Natural Science Foundation of China(No.524B2078,12426307,51906203)Guangdong Major Project of Basic and Applied Basic Research(2023B0303000002)+6 种基金Guangdong Basic and Applied Basic Research Foundation(2023B1515120005)Natural Science Foundation of Shenzhen(JCYJ20241202125327036,JCYJ20240813100103005)Shenzhen Engineering Research Center of Redox Flow Battery for Energy Storage(XMHT20230208003)Research Project on Medium-and Long-Duration Flow Battery Energy Storage Technology(2024KJTW0015)China Association for Science and Technology(OR2308010)High level of special funds(G03034K001)supported by the Center for Computational Science and Engineering at the Southern University of Science and Technology.
文摘Long-duration energy storage has become critical for renewable energy integration.While redox flow batteries,especially vanadium-based systems,are scaling up in capacity,their performance at the stack level remains insufficiently optimized,demanding more profound mechanistic studies and engineering refinements.To address the difficulties in resolving the flow inhomogeneity at the stack scale,this study establishes a multi-physics field coupling model and analyzes the pressure distributions,flow rate differences,active substance concentration,and electrochemical characteristics.The results show that the uneven cell pressure distribution is a key factor affecting the consistency of the system performance,and the increase in the flow rate improves the reactant homogeneity,with both the average concentration and the uniformity factor increasing with the flow rate.In contrast,high current densities lead to an increased imbalance between electrochemical depletion and reactant replenishment,resulting in a significant decrease in reactant concentration in the under-ribs region.In addition,a higher flow rate can expand the high-current-density region where the stack operates efficiently.This study provides a theoretical basis for optimizing the design of the stack components.
文摘The investigation of Beltrami flows is important for the research on the mechanism of turbulent structure.In this paper the general solutions of the Beltrami flows are given,which depend explicitly on the solutions of three independent Helmholtz equations with scalar unknowns.Velocity fields of Beltrami flows can then be obtained explicitly after the application of some curl operations on the solutions of Helmholtz equa- tions.On the basis of the exact solutions of Euler equations given above,we obtain one kind of exact solu- tions of non-steady Navier-Stokes equations which are also the Beltrami flows.Some interesting examples of Beltrami flows other than'ABC flows','Kolmogolov flows','Rayleigh-Bernard flows','Q-flows' are given.The detailed analytic results of these examples will be published in the near future.
基金Supported by the National Natural Science Foundation of China(Nos.42250710152,42192562)the Southern Laboratory of Ocean Science and Engineering(Guangdong Zhuhai)(No.SML 2020 SP 007)。
文摘Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealized numerical experiment using the large eddy simulation method was conducted to analyze the three-dimensional flow field structure and the mechanism for dense water subsidence.Specifically,a negative salt flux is set at the sea surface,in which salt flux enters the sea surface to simulate the icing and salting-out phenomena that occur at high latitudes.Results show that the mean-state 3D flow field of dense water subsidence exhibits a hollow conical distribution.The horizontal flow field is characterized by a cyclonic vortex that driven primarily by the pressure gradient and influenced by the Coriolis effect.Moreover,the inverse vertical pressure gradient generated by this vortex inhibits the sinking of the plume,leading to its off-axis deflection and the development of an anticyclonic precession.In addition,the impact of rotation on the structure of a sinking plume within a stratified environment is discussed.Both horizontal vortex intensity and cone angle of the hollow cone flow field are increased with increasing rotation rate,resulting in a decrease in the plume’s maximum sinking depth.Variances in rotation direction cause the horizontal vortex and sinking plumes of dense water in the northern and southern hemispheres to rotate in opposite directions.
基金financially supported by the Science and Technology Project of the Ministry of Transport(Grant No.2013328352570)
文摘In order to accurately simulate strong three-dimensional （3-D） free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional （2-D） grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.
基金Projects supported by the National Natural Science Foundation of China(Grant Nos.51909160,51879159)the National Key Research and Development Program of China(Grant Nos.2019YFB1704200,2019YFC0312400).
文摘Recently,the Liutex-based force field models are proposed and applied to some two-dimensional flows to explore the possibility of direct vortex control strategies.In this paper,we applied the Liutex-based centripetal force model to three-dimensional flow around a circular cylinder at Re=3900.The turbulent flows around the three-dimensional cylinder are modeled by the delayed detached-eddy simulations based on the two-equation k-ω shear stress transport(SST)model.The cylinder flow without adding the force field model is simulated in the first place.The statistical flow quantities are compared with experimental data to validate the accuracy of the current numerical models.Then the force field model is applied to the momentum equation in a specific control region near the cylinder to study the change of flow field and hydrodynamic performance.The effect of different control region locations with the same force field strength is discussed to find out the most influenced locations on the drag and lift.Based on the screened control region,the influence source term strength is further studied.Variations on the global statistical flow quantities,local instantaneous three-dimensional flow structures are analyzed and discussed.It is concluded that the Liutex-based force field model can effectively change the hydrodynamic behavior of flow past a cylinder and may provide a new direction for performance improvement and optimization.
文摘In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.
基金Supported by the Science Foundations of Laboratory of Computational PhysicalScience Foundation of China Academy of Engineering Physics under Grant Nos. 2009A0102005, 2009B0101012National Natural Science Foundation under Grant Nos. 10775018, 11074300, and 1107521 of China
文摘A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.
文摘A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.
