The motion of the fins and control surfaces of underwater vehicles in a fluid is an interesting and challenging research subject.Typically the effect of fin oscillations on the fluid flow around such a body is highly ...The motion of the fins and control surfaces of underwater vehicles in a fluid is an interesting and challenging research subject.Typically the effect of fin oscillations on the fluid flow around such a body is highly unsteady, generating vortices and requiring detailed analysis of fluid-structure interactions.An understanding of the complexities of such flows is of interest to engineers developing vehicles capable of high dynamic performance in their propulsion and maneuvering.In the present study, a CFD based RANS simulation of a 3-D fin body moving in a viscous fluid was developed.It investigated hydrodynamic performance by evaluating the hydrodynamic coefficients (lift, drag and moment) at two different oscillating frequencies.A parametric analysis of the factors that affect the hydrodynamic performance of the fin body was done, along with a comparison of results from experiments.The results of the simulation were found in close agreement with experimental results and this validated the simulation as an effective tool for evaluation of the unsteady hydrodynamic coefficients of 3-D fins.This work can be further be used for analysis of the stability and maneuverability of fin actuated underwater vehicles.展开更多
This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensiona...This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensional(3-D)steady Euler equations can be projected onto the unit sphere and the state of fluid can be characterized by the polar and azimuthal angles.Given a segment smooth curve as a conical-sonic line in the polar-azimuthal angle plane,we construct a classical conical-supersonic solution near the curve under some reasonable assumptions.To overcome the difficulty caused by the parabolic degeneracy,we apply the characteristic decomposition technique to transform the Euler equations into a new degenerate hyperbolic system in a partial hodograph plane.The singular terms are isolated from the highly nonlinear complicated system and then can be handled successfully.We establish a smooth local solution to the new system in a suitable weighted metric space and then express the solution in terms of the original variables.展开更多
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.展开更多
In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since...In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180v bend, 0.6 m to 0.45 m wide, was simulated using the SSI1M 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widtbwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel's centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.展开更多
According to turbulent theory and characteristics of beam blank continuous casting, 3-D model to represent the flow of beam blank mould is established. The predicted results indicate that the exit obliquity of up 15...According to turbulent theory and characteristics of beam blank continuous casting, 3-D model to represent the flow of beam blank mould is established. The predicted results indicate that the exit obliquity of up 15°(+15°) should be adopted, which will benefit the floatation of non-metallic inclusion and purification of the molten steel. When the nozzle angle is 120°, the flow pattern is reasonable. Proper nozzle depth can be 200mm. Turbulent kinetic of meniscus can be reduced by adopting the square nozzle and suitable area of side outlet when casting speed increases. The results are consistent with those of water model experiment, so the model is exact and reasonable. The model can provide important information for design of SEN and defining of immersion depth.展开更多
A three-dimensional viscous code has been developed to solve Reynolds-averaged Navier-Stokes equations. The governing equations in finite volume form are solved by two-step Runge-Kutta scheme with implicit residual sm...A three-dimensional viscous code has been developed to solve Reynolds-averaged Navier-Stokes equations. The governing equations in finite volume form are solved by two-step Runge-Kutta scheme with implicit residual smoothing. The eddy viscous is obtained using the Baldwin-Lomax model. A prediction of the 3-D turbulent flow and the performance in the “all-over controlled vortex distribution” centrifugal impeller with a vaneless diffuser has been made for the compressor at design and off-design condition. The predicted effi-ciency is a little higher than the experiment data. These results suggest that the present calculation code is able to determine the flow development in the impeller and also the turbulence model in the centrifugal im-peller should be improved.展开更多
A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing t...A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.展开更多
In the current study, a numerical investigation of three-dimensional combined convection-radiation heat transfer over an inclined forward facing step (FFS) in a horizontal rectangular duct is presented. The fluid is t...In the current study, a numerical investigation of three-dimensional combined convection-radiation heat transfer over an inclined forward facing step (FFS) in a horizontal rectangular duct is presented. The fluid is treated as a gray, absorbing, emitting and scattering medium. To simulate the incline surface of FFS, the blocked-off method is employed in this study. The set of governing equations for gas flow are solved numerically using the CFD technique to obtain the temperature and velocity fields. Since the gas is considered as a radiating medium, all of the convection, conduction and radiation heat transfer mechanisms are presented in the energy equation. For computation of radiative term in energy equation, the radiative transfer equation (RTE) is solved numerically by the discrete ordinates method (DOM) to find the divergence of radiative heat flux distribution inside the radiating medium. The effects of optical thickness, radiation-conduction parameter and albedo coefficient on heat transfer behavior of the system are carried out.展开更多
Various models exist to explain the formation of the Tibetan Plateau,including“tectonic escape”,“pure shear thickening”,“convective removal of the lithospheric mantle”,and“lower crustal flow”model.The first tw...Various models exist to explain the formation of the Tibetan Plateau,including“tectonic escape”,“pure shear thickening”,“convective removal of the lithospheric mantle”,and“lower crustal flow”model.The first two models are primarily constructed on pure mechanical models but are unable to reasonably explain the tension and shear phenomena inside the plateau.The latter two are rheological dynamic models based on deep geophysical observations.However,the spatial range of the lower crustal flow and its role in the plateau formation/uplift remain controversial.Five multi-terrane viscoplastic thermomechanical models were constructed to simulate the uplift and lithospheric structure change of the Tibetan Plateau during the post-collision stage(since 35 Ma)under the convergence of the Indian Plate.Results show that the plateau's formation begins with crustal thickening,blocked by strong terranes at the northern plateau,and expanded laterally to the east.The lithosphere thickens gradually and experiences delamination at its base,elevating temperature within the crust and forming partial melting layers in the central plateau.As convergence persists on the southern side,the northern plateau's lithosphere bends downward and undergoes delamination,further heating the crust and promoting the northward and eastward flow of partial melting layers,leading to secondary uplift around the plateau.展开更多
Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package PO...Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package POLYFLOW. Based on the velocity fields calculated, the particle trajectories in both machines are visualized using particle tracking technique. The numerical results indicate that the flow patterns in three-screw extruders are similar to those in twin-screw extruders. The triangularly-arranged three-screw extruder has the largest pumping capacity and also the highest extrusion stability in terms of flowrate fluctuation with screw rotation. The instantaneous mixing and cumulative residence time distribution (RTD) characteristics are also analyzed and compared with traditional intermeshing co-rotating twin-screw extruders. It is shown that the start section of the cumulative RTD curve for the triangularly-arranged machine has a small shoulder, which is attributed to the faster flow in the central region of this type of extruder.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
The generalized lattice Boltzmann equation(GLBE),with the addition of the standard Smagorinsky subgrid-stress(SGS) model,has been proved that it is more suitable for simulating high Reynolds number turbulent flows whe...The generalized lattice Boltzmann equation(GLBE),with the addition of the standard Smagorinsky subgrid-stress(SGS) model,has been proved that it is more suitable for simulating high Reynolds number turbulent flows when compared with the lattice BGK Boltzmann equation(LBGK).However,the computing efficiency of lattice Boltzmann method(LBM) is too low to make it for practical applications,unless using a massive parallel computing clusters facility.In this study,the massive parallel computing power from an inexpensive graphic processor unit(GPU) and a typical personal computer has been developed for improving the computing efficiency,more than 100 times.This developed three-dimensional(3-D) GLBE-SGS model,with the D3Q19 scheme for simplifying collision and streaming courses,has been successfully used to study 3-D rectangular cavity flows with Reynolds number up to 10000.展开更多
A generalized hend flow model, treating a 90° single bend and 60° continuous hends, was designed to quantitatively describe 3-D turhulenee mechanism of circulating notfully-developed flow in open channels wi...A generalized hend flow model, treating a 90° single bend and 60° continuous hends, was designed to quantitatively describe 3-D turhulenee mechanism of circulating notfully-developed flow in open channels with hends. The 3-D fluctuating veloeities of turbulent flow were measured and analyzed with a 3 D acoustic-Doppler velocimeter. Formula for 3 D turbulent intensity was derived using the dimension analysis approaeh. Expressions of vertical turbulent intensity distributions were obtained with the multivariant-rcgression theo ry, whieh agree with experiment data. Distrihutions of turbulent intensity and turbulent stress were characterized, and their relationships were concluded. In the bend-turbulent flow core region, longitudinal and lateral turbulent-intensity distri hutions are coincident with linear distribution, hut in nearwall region are coincident with the Gamma distribution. Verotical turbulent intensity distributions are coincident with the Rayleigh distribution. Herein, it is concluded that the bend turbulence is anisotropic.展开更多
An efficient numerical method for calculating the three-dimensional transonic flows in turbomachinery is proposed. Instead of the Euler equation, streamsurface-governing equations are deduced in the generalized von Mi...An efficient numerical method for calculating the three-dimensional transonic flows in turbomachinery is proposed. Instead of the Euler equation, streamsurface-governing equations are deduced in the generalized von Mises coordinate system to reflect the flow feature in turbomachinery. Its main advantage is that it is easier to specify more reasonable initial values, i.e. initial streamsurface position, thus accelerating the convergence rate of the iteration process. Moreover, to use the generalized von Mises coordinates makes the present method capable of incorporating the calculation of the flow field, design and modification of the blade contour into a unified algorithm. A rotated finite difference scheme for the streamsurface-governing equations is constructed, and a new measure is presented to deal with the double-value problem of the velocity and density caused by the application of the stream functions as coordinates in the transonic flow. Three test cases were considered with the present approach to demonstrate the solution method for 3-D inviscid flow analysis. Numerical results confirm that the present method has rapid convergence and high accuracy.展开更多
Most natural rivers are curved channels, where the turbulent flows have a complex helical pattern, as has been extensively studied both numerically and experimentally. The helical flow structure in curved channels ha...Most natural rivers are curved channels, where the turbulent flows have a complex helical pattern, as has been extensively studied both numerically and experimentally. The helical flow structure in curved channels has an important bearing on sediment transport, riverbed evolution, and pollutant transport study. In this article, different turbulence closure schemes i.e., the mixing-length model and the k-ε model with different pressure solution techniques i. e., hydrostatic assumptions and dynamic pressure treatments are applied to study the helical secondary flows in an experiment curved channel. The agreements of vertically-averaged velocities between the simulated results obtained by using different turbulence models with different pressure solution techniques and the measured data are satisfactory. Their discrepancies with respect to surface elevations, superelevations and secondary flow patterns are discussed.展开更多
This study is aimed at using blade 3-D optimization to control corner flows in the high through-flow fan/booster of a high bypass ratio commercial turbofan engine. Two kinds of blade 3-D optimization, end-bending and ...This study is aimed at using blade 3-D optimization to control corner flows in the high through-flow fan/booster of a high bypass ratio commercial turbofan engine. Two kinds of blade 3-D optimization, end-bending and bow, are focused on. On account of the respective operation mode and environment, the approach to 3-D aerodynamic modeling of rotor blades is different from stator vanes. Based on the understanding of the mechanism of the corner flow and the consideration of intensity problem for rotors, this paper uses a variety of blade 3-D optimization approaches, such as loading distribution optimization, perturbation of departure angles and stacking-axis manipulation, which are suitable for rotors and stators respectively. The obtained 3-D blades and vanes can improve the corner flow features by end-bending and bow effects. The results of this study show that flows in corners of the fan/booster, such as the fan hub region, the tip and hub of the vanes of the booster, are very complex and dominated by 3-D effects. The secondary flows there are found to have a strong detrimental effect on the compressor performance. The effects of both end-bending and bow can improve the flow separation in corners, but the specific ways they work and application scope are somewhat different. Redesigning the blades via blade 3-D optimization to control the corner flow has effectively reduced the loss generation and improved the stall margin by a large amount.展开更多
In this article,the anaglyph video maker is employed for generating realistic 3-D flows and the software FlowAnimator is developed using that technology.Based on Microsoft Windows Presentation Foundation(WPF),the re...In this article,the anaglyph video maker is employed for generating realistic 3-D flows and the software FlowAnimator is developed using that technology.Based on Microsoft Windows Presentation Foundation(WPF),the real 3-D scene is set up and marker particles are distributed in it randomly in order to create a more natural flow scenario.The trajectory of the particle motion is calculated with Lagrangian description in 3-D space.During the simulation,the viewport can be changed in order to focus on different parts of the model by panning,zooming,rotating and inclination variation etc.Marker particles may appear in different shapes:spheres,tracking-balls,cylinders and ribbons in order to fit different flows.It is the first time that the video anaglyph technology is employed in the 3-D hydrodynamic simulation,which removes the obstacles for 3-D scenes to be rendered on a flat-panel display.展开更多
The 3-D characteristics of the water-air flow patterns in a corridor-shaped air-cushion surge chamber during hydraulic transients need to be considered in the shape optimization. To verify the reliability of the water...The 3-D characteristics of the water-air flow patterns in a corridor-shaped air-cushion surge chamber during hydraulic transients need to be considered in the shape optimization. To verify the reliability of the water-air two-phase model, namely, the volume of fluid model, the process of charging water into a closed air chamber is successfully simulated. Using the model, the 3-D flow characteristics under the load rejection and acceptance conditions within the air-cushion surge chamber of a specific hydro- power station are studied. The free surface waves, the flow patterns, and the pressure changes during the surge wave process are analyzed in detail. The longitudinal flow of water in the long corridor-shaped surge chamber is similar to the open channel flow with respect to the wave propagation, reflection and superposition characteristics. The lumped parameters of the 3-D numerical simulation agree with the results of a 1-D calculation of hydraulic transients in the whole water conveying system, which validates the 3-D method. The 3-D flow structures obtained can be applied to the shape optimization of the chamber.展开更多
A 3-D model based on the Reynolds equations with closed k-ε turbulence model is presented in this paper,which can be used to predict surface water flow in open channels.In- stead of the“rigid lid”approximation,the ...A 3-D model based on the Reynolds equations with closed k-ε turbulence model is presented in this paper,which can be used to predict surface water flow in open channels.In- stead of the“rigid lid”approximation,the solution of the free surface equation is implemented in the velocity-pressure iterative procedure on the basis of the conventional SIMPLE method.This model was used to compute the flow in rectangular channels with trenches dredged across the bot- tom.The velocity,eddy viscosity coefficient,turbulent shear stress,turbulent kinetic energy and elevation of the free surface over the trenches dredged in the main channel,can be obtained. The computed results are in good agreement with existing experimentaing data.展开更多
基金Supported by the National Natural Science Foundation of China under Grant No.50879014
文摘The motion of the fins and control surfaces of underwater vehicles in a fluid is an interesting and challenging research subject.Typically the effect of fin oscillations on the fluid flow around such a body is highly unsteady, generating vortices and requiring detailed analysis of fluid-structure interactions.An understanding of the complexities of such flows is of interest to engineers developing vehicles capable of high dynamic performance in their propulsion and maneuvering.In the present study, a CFD based RANS simulation of a 3-D fin body moving in a viscous fluid was developed.It investigated hydrodynamic performance by evaluating the hydrodynamic coefficients (lift, drag and moment) at two different oscillating frequencies.A parametric analysis of the factors that affect the hydrodynamic performance of the fin body was done, along with a comparison of results from experiments.The results of the simulation were found in close agreement with experimental results and this validated the simulation as an effective tool for evaluation of the unsteady hydrodynamic coefficients of 3-D fins.This work can be further be used for analysis of the stability and maneuverability of fin actuated underwater vehicles.
基金the two referees for very helpful comments and suggestions to improve the quality of the paper.This work was partially supported by the Natural Science Foundation of Zhejiang province of China(LY21A010017)the National Natural Science Foundation of China(12071106,12171130).
文摘This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensional(3-D)steady Euler equations can be projected onto the unit sphere and the state of fluid can be characterized by the polar and azimuthal angles.Given a segment smooth curve as a conical-sonic line in the polar-azimuthal angle plane,we construct a classical conical-supersonic solution near the curve under some reasonable assumptions.To overcome the difficulty caused by the parabolic degeneracy,we apply the characteristic decomposition technique to transform the Euler equations into a new degenerate hyperbolic system in a partial hodograph plane.The singular terms are isolated from the highly nonlinear complicated system and then can be handled successfully.We establish a smooth local solution to the new system in a suitable weighted metric space and then express the solution in terms of the original variables.
基金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.
文摘In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180v bend, 0.6 m to 0.45 m wide, was simulated using the SSI1M 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widtbwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel's centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.
文摘According to turbulent theory and characteristics of beam blank continuous casting, 3-D model to represent the flow of beam blank mould is established. The predicted results indicate that the exit obliquity of up 15°(+15°) should be adopted, which will benefit the floatation of non-metallic inclusion and purification of the molten steel. When the nozzle angle is 120°, the flow pattern is reasonable. Proper nozzle depth can be 200mm. Turbulent kinetic of meniscus can be reduced by adopting the square nozzle and suitable area of side outlet when casting speed increases. The results are consistent with those of water model experiment, so the model is exact and reasonable. The model can provide important information for design of SEN and defining of immersion depth.
文摘A three-dimensional viscous code has been developed to solve Reynolds-averaged Navier-Stokes equations. The governing equations in finite volume form are solved by two-step Runge-Kutta scheme with implicit residual smoothing. The eddy viscous is obtained using the Baldwin-Lomax model. A prediction of the 3-D turbulent flow and the performance in the “all-over controlled vortex distribution” centrifugal impeller with a vaneless diffuser has been made for the compressor at design and off-design condition. The predicted effi-ciency is a little higher than the experiment data. These results suggest that the present calculation code is able to determine the flow development in the impeller and also the turbulence model in the centrifugal im-peller should be improved.
基金financially supported by the Science and Technology Project of the Ministry of Transport (Grant No. 2011329224170)
文摘A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.
文摘In the current study, a numerical investigation of three-dimensional combined convection-radiation heat transfer over an inclined forward facing step (FFS) in a horizontal rectangular duct is presented. The fluid is treated as a gray, absorbing, emitting and scattering medium. To simulate the incline surface of FFS, the blocked-off method is employed in this study. The set of governing equations for gas flow are solved numerically using the CFD technique to obtain the temperature and velocity fields. Since the gas is considered as a radiating medium, all of the convection, conduction and radiation heat transfer mechanisms are presented in the energy equation. For computation of radiative term in energy equation, the radiative transfer equation (RTE) is solved numerically by the discrete ordinates method (DOM) to find the divergence of radiative heat flux distribution inside the radiating medium. The effects of optical thickness, radiation-conduction parameter and albedo coefficient on heat transfer behavior of the system are carried out.
基金sponsored by the National Key R&D Program of China(No.2021YFA0715100)the Shenzhen Fundamental Research Program,China(No.JCYJ20220818102601004)+1 种基金the National Natural Science Foundation of China(No.41774145)the Pre-research Project on Civil Aerospace Technologies(No.D020101)of CNSA。
文摘Various models exist to explain the formation of the Tibetan Plateau,including“tectonic escape”,“pure shear thickening”,“convective removal of the lithospheric mantle”,and“lower crustal flow”model.The first two models are primarily constructed on pure mechanical models but are unable to reasonably explain the tension and shear phenomena inside the plateau.The latter two are rheological dynamic models based on deep geophysical observations.However,the spatial range of the lower crustal flow and its role in the plateau formation/uplift remain controversial.Five multi-terrane viscoplastic thermomechanical models were constructed to simulate the uplift and lithospheric structure change of the Tibetan Plateau during the post-collision stage(since 35 Ma)under the convergence of the Indian Plate.Results show that the plateau's formation begins with crustal thickening,blocked by strong terranes at the northern plateau,and expanded laterally to the east.The lithosphere thickens gradually and experiences delamination at its base,elevating temperature within the crust and forming partial melting layers in the central plateau.As convergence persists on the southern side,the northern plateau's lithosphere bends downward and undergoes delamination,further heating the crust and promoting the northward and eastward flow of partial melting layers,leading to secondary uplift around the plateau.
文摘Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package POLYFLOW. Based on the velocity fields calculated, the particle trajectories in both machines are visualized using particle tracking technique. The numerical results indicate that the flow patterns in three-screw extruders are similar to those in twin-screw extruders. The triangularly-arranged three-screw extruder has the largest pumping capacity and also the highest extrusion stability in terms of flowrate fluctuation with screw rotation. The instantaneous mixing and cumulative residence time distribution (RTD) characteristics are also analyzed and compared with traditional intermeshing co-rotating twin-screw extruders. It is shown that the start section of the cumulative RTD curve for the triangularly-arranged machine has a small shoulder, which is attributed to the faster flow in the central region of this type of extruder.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
基金supported by the Virginia Institute of Marine Science,College of William and Mary for the Study Environmentthe National Natural Science Foundation of China(Grant No.50679008)
文摘The generalized lattice Boltzmann equation(GLBE),with the addition of the standard Smagorinsky subgrid-stress(SGS) model,has been proved that it is more suitable for simulating high Reynolds number turbulent flows when compared with the lattice BGK Boltzmann equation(LBGK).However,the computing efficiency of lattice Boltzmann method(LBM) is too low to make it for practical applications,unless using a massive parallel computing clusters facility.In this study,the massive parallel computing power from an inexpensive graphic processor unit(GPU) and a typical personal computer has been developed for improving the computing efficiency,more than 100 times.This developed three-dimensional(3-D) GLBE-SGS model,with the D3Q19 scheme for simplifying collision and streaming courses,has been successfully used to study 3-D rectangular cavity flows with Reynolds number up to 10000.
文摘A generalized hend flow model, treating a 90° single bend and 60° continuous hends, was designed to quantitatively describe 3-D turhulenee mechanism of circulating notfully-developed flow in open channels with hends. The 3-D fluctuating veloeities of turbulent flow were measured and analyzed with a 3 D acoustic-Doppler velocimeter. Formula for 3 D turbulent intensity was derived using the dimension analysis approaeh. Expressions of vertical turbulent intensity distributions were obtained with the multivariant-rcgression theo ry, whieh agree with experiment data. Distrihutions of turbulent intensity and turbulent stress were characterized, and their relationships were concluded. In the bend-turbulent flow core region, longitudinal and lateral turbulent-intensity distri hutions are coincident with linear distribution, hut in nearwall region are coincident with the Gamma distribution. Verotical turbulent intensity distributions are coincident with the Rayleigh distribution. Herein, it is concluded that the bend turbulence is anisotropic.
文摘An efficient numerical method for calculating the three-dimensional transonic flows in turbomachinery is proposed. Instead of the Euler equation, streamsurface-governing equations are deduced in the generalized von Mises coordinate system to reflect the flow feature in turbomachinery. Its main advantage is that it is easier to specify more reasonable initial values, i.e. initial streamsurface position, thus accelerating the convergence rate of the iteration process. Moreover, to use the generalized von Mises coordinates makes the present method capable of incorporating the calculation of the flow field, design and modification of the blade contour into a unified algorithm. A rotated finite difference scheme for the streamsurface-governing equations is constructed, and a new measure is presented to deal with the double-value problem of the velocity and density caused by the application of the stream functions as coordinates in the transonic flow. Three test cases were considered with the present approach to demonstrate the solution method for 3-D inviscid flow analysis. Numerical results confirm that the present method has rapid convergence and high accuracy.
基金supported by the National Natural Science Foundation of China (Grant No.50479034)the Natural Science Foundation of Tianjin (Grant No.09YFSZSF02100)+1 种基金the financial support of the USDA Agriculture Research Service under Specific Research Agreement (Grant No. 58-6408-2-0062)(monitored by the USDA-ARS National Sedimentation Laboratory)the US State Department Agency for International Development under Agreement (Grant No.EE-G-00-02-00015-00) and the University of Mississippi
文摘Most natural rivers are curved channels, where the turbulent flows have a complex helical pattern, as has been extensively studied both numerically and experimentally. The helical flow structure in curved channels has an important bearing on sediment transport, riverbed evolution, and pollutant transport study. In this article, different turbulence closure schemes i.e., the mixing-length model and the k-ε model with different pressure solution techniques i. e., hydrostatic assumptions and dynamic pressure treatments are applied to study the helical secondary flows in an experiment curved channel. The agreements of vertically-averaged velocities between the simulated results obtained by using different turbulence models with different pressure solution techniques and the measured data are satisfactory. Their discrepancies with respect to surface elevations, superelevations and secondary flow patterns are discussed.
基金supported by National Natural Science Foundation of China (51006005,50736007)"Fan-Zhou" Youth Foundation(20100401)
文摘This study is aimed at using blade 3-D optimization to control corner flows in the high through-flow fan/booster of a high bypass ratio commercial turbofan engine. Two kinds of blade 3-D optimization, end-bending and bow, are focused on. On account of the respective operation mode and environment, the approach to 3-D aerodynamic modeling of rotor blades is different from stator vanes. Based on the understanding of the mechanism of the corner flow and the consideration of intensity problem for rotors, this paper uses a variety of blade 3-D optimization approaches, such as loading distribution optimization, perturbation of departure angles and stacking-axis manipulation, which are suitable for rotors and stators respectively. The obtained 3-D blades and vanes can improve the corner flow features by end-bending and bow effects. The results of this study show that flows in corners of the fan/booster, such as the fan hub region, the tip and hub of the vanes of the booster, are very complex and dominated by 3-D effects. The secondary flows there are found to have a strong detrimental effect on the compressor performance. The effects of both end-bending and bow can improve the flow separation in corners, but the specific ways they work and application scope are somewhat different. Redesigning the blades via blade 3-D optimization to control the corner flow has effectively reduced the loss generation and improved the stall margin by a large amount.
文摘In this article,the anaglyph video maker is employed for generating realistic 3-D flows and the software FlowAnimator is developed using that technology.Based on Microsoft Windows Presentation Foundation(WPF),the real 3-D scene is set up and marker particles are distributed in it randomly in order to create a more natural flow scenario.The trajectory of the particle motion is calculated with Lagrangian description in 3-D space.During the simulation,the viewport can be changed in order to focus on different parts of the model by panning,zooming,rotating and inclination variation etc.Marker particles may appear in different shapes:spheres,tracking-balls,cylinders and ribbons in order to fit different flows.It is the first time that the video anaglyph technology is employed in the 3-D hydrodynamic simulation,which removes the obstacles for 3-D scenes to be rendered on a flat-panel display.
基金support by the National Natural Science Foundation of China(Grant Nos.51039005,50909076)
文摘The 3-D characteristics of the water-air flow patterns in a corridor-shaped air-cushion surge chamber during hydraulic transients need to be considered in the shape optimization. To verify the reliability of the water-air two-phase model, namely, the volume of fluid model, the process of charging water into a closed air chamber is successfully simulated. Using the model, the 3-D flow characteristics under the load rejection and acceptance conditions within the air-cushion surge chamber of a specific hydro- power station are studied. The free surface waves, the flow patterns, and the pressure changes during the surge wave process are analyzed in detail. The longitudinal flow of water in the long corridor-shaped surge chamber is similar to the open channel flow with respect to the wave propagation, reflection and superposition characteristics. The lumped parameters of the 3-D numerical simulation agree with the results of a 1-D calculation of hydraulic transients in the whole water conveying system, which validates the 3-D method. The 3-D flow structures obtained can be applied to the shape optimization of the chamber.
文摘A 3-D model based on the Reynolds equations with closed k-ε turbulence model is presented in this paper,which can be used to predict surface water flow in open channels.In- stead of the“rigid lid”approximation,the solution of the free surface equation is implemented in the velocity-pressure iterative procedure on the basis of the conventional SIMPLE method.This model was used to compute the flow in rectangular channels with trenches dredged across the bot- tom.The velocity,eddy viscosity coefficient,turbulent shear stress,turbulent kinetic energy and elevation of the free surface over the trenches dredged in the main channel,can be obtained. The computed results are in good agreement with existing experimentaing data.