水动力学研究与进展:英文版最新文献

In this paper, the 3-D turbulent flow around two side-by-side circular cylinders of different diameters, at sub-critical Reynolds number (Re = 3900), is numerically simulated by the large eddy simulation (LES). The spacing ratios (T/D) between the two cylinders are considered in four cases (T/D = 1.2, 1.5, 1.8 and 2.7) to study the vortex shedding and turbulent properties in the flow field. The main results are focused on the drag and lift coefficients, the vortex shedding frequency, the coherent structure, and the scale properties. It is shown that when T/D is equal to 1.2, the vortex shedding of the main cylinder is strongly suppressed by the small cylinder, the drag and lift coefficients of the main cylinder are smaller than those in other three cases. While T/D is equal to 1.5, the vortex shedding of the main cylinder can be improved, the drag and lift coefficients of the main cylinder are larger than those in other three cases. The empirical mode decomposition (EMD) method is applied to decompose the velocity signals traced by the LES. It is shown that there is a linear relationship between the mean period and the mode in the semi-log coordinates. The vortex shedding period of the main cylinder is consistent with the period of the restructured coherent structures quantitatively.

The slit-type energy dissipater (STED) is widely used in hydraulic projects of high water head, large discharge, and narrow river valley, thanks to its simple structure and high efficiency. However, the water wing caused by the shock waves in the contraction section of the STED may bring about harmful effects. A coefficient is introduced for the application of Ippen's theory in the STED. The expression of the coefficient is experimentally obtained. Simplified formulas to calculate the shock wave angle and the water wing scope are theoretically derived, with relative errors within 5%.

The capability of the smoothed-particle hydrodynamics (SPH) method to model supercritical flow in circular pipe bends is considered. The standard SPH method, which makes use of dynamic boundary particles (DBP), is supplemented with the original algorithm for the treatment of open boundaries. The method is assessed through a comparison with measured free-surface profiles in a pipe bend, and already proposed regression curves for estimation of the flow-type in a pipe bend. The sensitivity of the model to different parameters is also evaluated. It is shown that an adequate choice of the artificial viscosity coefficient and the initial particle spacing can lead to correct presentation of the flow-type in a bend. Due to easiness of its implementation, the SPH method can be efficiently used in the design of circular conduits with supercritical flow in a bend, such as tunnel spillways, and bottom outlets of dams, or storm sewers.

The best hydraulic channel section makes the maximum flow capacity for the same flow cross-area, and the minimum cross-area and wetted perimeter for the same discharge. The construction cost can be reduced nearly to the minimum at the same time. The horizontal bottom parabolic section (HBP section) is a composite section. It is important for design to find the best combination form of the horizontal bottom and the parabolic sides. This paper studies the best hydraulic section and its hydraulic characteristics. The explicit formulae are proposed to determine the dimensions and the best combination form of the horizontal bottom and the parabolic sides. These explicit formulae and the parameters make it easy to design the channel. It is shown that the ratios of the surface width to the depth and the bottom width to the depth are constant for the best hydraulic section. The comparisons with the classic parabolic, rectangular, trapezoid, triangular, semi-cubic and horizontal-bottomed semi-cubic sections show that the HBP section has the largest flow capacity and the shortest wetted perimeter for the same flow area, and has the smallest flow area for the same discharge. It is indicated that the parabolic side parts of the best hydraulic HBP section are different from those of the classic section. The results of the best hydraulic section of the classic parabolic channel cannot be applied directly to the HBC section.

The cavitation shedding flow around a 3-D Clark-Y hydrofoil is simulated by using an improved filter-based model (FBM) and a mass transfer cavitation model with the consideration of the maximum density ratio effect between the liquid and the vapor. The unsteady cloud cavity shedding features around the Clark-Y hydrofoil are accurately captured based on an improved FBM model and a suitable maximum density ratio. Numerical results show that the predicted cavitation patterns and evolutions compare well with the experimental visualizations, and the prediction errors of the time-averaged lift coefficient, drag coefficient and Strouhal number for the cavitation number, the angle of attack at a Reynolds number are only 3.29%, 2.36% and 9.58%, respectively. It is observed that the cavitation shedding flow patterns are closely associated with the vortex structures identified by the criterion method. The predicted cloud cavitation shedding flow shows clearly three typical stages: (1) Initiation of the attached sheet cavity, the growth toward the trailing edge. (2) The formation and development of the re-entrant jet flow. (3) Large scale cloud cavity sheds downstream. Numerical results also indicate that the non-uniform adverse pressure gradient is the main driving force of the re-entrant jet, which results in the U-shaped cavity and the 3-D bubbly structure during the cloud cavity shedding.

A parallel direct-forcing fictitious domain method is applied in fully-resolved numerical simulations of particle-laden turbulent flows in a square duct. The effects of finite-size heavy particles on the mean secondary flow, the mean streamwise velocity, the root-mean-square velocity fluctuation, and the particle concentration distribution are investigated at the friction Reynolds number of 150, the particle volume fraction of 2.36%, the particle diameter of 0.1 duct width, and the Shields number ranging from 1.0 to 0.2. Our results show that the particle sedimentation breaks the up-down symmetry of the mean secondary vortices, and results in a stronger secondary-flow circulation which transports the fluids downward in the bulk center region and upward along the side walls at a low Shields number. This circulation has a significant impact on the distribution of the mean streamwise velocity, whose maximum value occurs in the lower half duct, unlike in the plane channel case. The flow resistance is increased and the turbulence intensity is reduced, as the Shields number is decreased. The particles accumulate preferentially at the face center of the bottom wall, due to the effect of the mean secondary flow. It is observed that the collision model has an important effect on the results, but does not change the results qualitatively.

In this work, we perform a numerical study of a water flow over a stepped spillway. This flow is described by the Reynolds averaged Navier-Stokes equation (RANS) associated with the turbulence model. These equations are solved using a commercial software based on the finite volume scheme and an unstructured mesh. The air-water flow was modeled using volume of fluid (VOF) and multiphasic methods. The characteristics of the flow were investigated including the total pressure, the velocity profile, etc.. We analyze the effects on the flow structure of the steps and countermarch inclination, the air injection through the countermarch into the water flow and the dynamics water discharges. Results show that the inclination of the countermarch relative to the vertical and the air injection into the water flow increase the total pressure in the neighbourhood of the steps.

A wall motion of a bubble freely oscillating in a liquid is studied from the point of view of energy conversions at different instants. It is shown that the time of the bubble oscillation can be divided into two distinct intervals. In the first long interval (here called PK and KP intervals) the prevailing energy conversion is between the potential energy of the bubble and the kinetic energy of the liquid. In the second short interval (here called KI and IK intervals) the kinetic energy of the liquid is transformed into the internal energy of the gas/vapor in the bubble interior and into some other forms of energy. By observing the bubble wall motion in the PK and KP intervals, it is shown that only the value of the maximum bubble radius in the corresponding oscillation can be determined. However, only the knowledge of the maximum bubble radii is insufficient for formulation of a correct theoretical model. Unfortunately this fact is often not noticed in the literature.

In this study, three-dimensional liquid-metal magnetohydrodynamic flows in a conduit with multiple channels under a uniform magnetic field are numerically investigated. The geometry of the conduit is of a four-parallel-channels system including one inflow channel and three outflow channels. The liquid-metal flows into the inflow channel, then turns through in the transition segment, finally flows into three different outflow channels. This kind of channel system can induce counter flow and co-flow, which is rarely investigated before though the conceptual designs of duct flow in the blanket have suggested this type of flow. A structured grid system is chosen after a series of mesh independence tests in the present study. The axial velocity in the side layer near the first partitioning wall, located between the inflow channel and the first outflow channel, is the highest with the lowest electric potential formed therein. The pressure almost linearly decreases in the main flow direction, except in the transition segment. Moreover, the pressure gradient in the first outflow channel is the largest among the three outflow channels. The interdependency of the current, fluid velocity, pressure, electric potential is examined in order to describe the electromagnetic characteristics of the liquid-metal flows.

A wing-body junction flow of a navigating underwater vehicle is considered to be a crucial source of the flow radiating acoustic noise, which attracts much research interest. In this paper, wing-plate junction flows are experimentally investigated in a low-speed wind tunnel by smoke-wire flow visualizations and time-resolved PIV measurements. To reveal the physical behavior of such flows, smoke-wire flow visualizations are conducted for a laminar wing-plate junction. A novel control strategy is proposed, to accurately locate the suction openings where the streamline is about to roll up to form a vortex in the turbulent junction flows. The control effect is discussed in perspectives of both the time-averaged and instantaneous flow fields.