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

According to Lorenz, chaotic dynamic systems have sensitive dependence on initial conditions (SDIC), i.e., the butterfly-effect: a tiny difference on initial conditions might lead to huge difference of computer-generated simulations after a long time. Thus, computer-generated chaotic results given by traditional algorithms in double precision are a kind of mixture of “true” (convergent) solution and numerical noises at the same level. Today, this defect can be overcome by means of the “clean numerical simulation” (CNS) with negligible numerical noises in a long enough interval of time. The CNS is based on the Taylor series method at high enough order and data in the multiple precision with large enough number of digits, plus a convergence check using an additional simulation with even smaller numerical noises. In theory, convergent (reliable) chaotic solutions can be obtained in an arbitrary long (but finite) interval of time by means of the CNS. The CNS can reduce numerical noises to such a level even much smaller than micro-level uncertainty of physical quantities that propagation of these physical micro-level uncertainties can be precisely investigated. In this paper, we briefly introduce the basic ideas of the CNS, and its applications in long-term reliable simulations of Lorenz equation, three-body problem and Rayleigh-Bénard turbulent flows. Using the CNS, it is found that a chaotic three-body system with symmetry might disrupt without any external disturbance, say, its symmetry-breaking and system-disruption are “self-excited”, i.e., out-of-nothing. In addition, by means of the CNS, we can provide a rigorous theoretical evidence that the micro-level thermal fluctuation is the origin of macroscopic randomness of turbulent flows. Naturally, much more precise than traditional algorithms in double precision, the CNS can provide us a new way to more accurately investigate chaotic dynamic systems.

In this paper, we investigate flows with moving contact lines on curved solid walls on a dual-resolution grid using a diffuse-interface immersed-boundary (DIIB) method. The dual-resolution grid, on which the flows are solved on a coarse mesh while the interface is resolved on a fine mesh, was shown to significantly improve the computational efficiency when simulating multiphase flows. On the other hand, the DIIB method is able to resolve dynamic wetting on curved substrates on a Cartesian grid, but it usually requires a mesh of high resolution in the vicinity of a moving contact line to resolve the local flow. In the present study, we couple the DIIB method with the dual-resolution grid, to improve the interface resolution for flows with moving contact lines on curved solid walls at an affordable cost. The dynamic behavior of moving contact lines is validated by studying drop spreading, and the numerical results suggest that the effective slip length λ_n can be approximated by 1.9Cn, where Cn is a dimensionless measure of the thickness of the diffuse interface. We also apply the method to drop impact onto a convex substrate, and the results on the dual-resolution grid are in good agreement with those on a single-resolution grid. It shows that the axisymmetric simulations using the DIIB method on the dual-resolution grid saves nearly 60% of the computational time compared with that on a single-resolution grid.

Wave breaking plays an important role in wave-structure interaction. A novel control volume finite element method with adaptive unstructured meshes is employed here to study 3-D breaking waves. The numerical framework consists of a “volume of fluid” type method for the interface capturing and adaptive unstructured meshes to improve computational efficiency. The numerical model is validated against experimental measurements of breaking wave over a sloping beach and is then used to study the breaking wave impact on a vertical circular cylinder on a slope. Detailed complex interfacial structures during wave impact, such as plunging jet formation and splash-up are captured in the simulation, demonstrating the capability of the present method.

The distribution of the phosphorus (P) adsorption in a bed sediment at channel confluences is an important issue for understanding the transport of contaminants in channel networks. In this study, the flow structure and its effect on the P distribution in the bed sediment were investigated in a 90° confluence flume. It is shown that the P adsorption amount in the sediment varies significantly in different hydrodynamic zones. The P adsorption amount is high in the flow separation zone where the horizontal velocity of the flow is very small, and it is low in the maximum velocity zone where the flow velocity reaches a maximum. The low P adsorption amount is observed in the downstream portion of the distorted shear layer, while the P enrichment is found in the upstream portion, as is related to the significant downwelling flow in this zone. Thus, the flow structure, especially, the flow velocity, has a significant effect on the distribution of the P adsorption in the bed sediment at the channel confluences.

The air-trapped vertical vortices (ATVVs) are easy to form above the throttled orifices in the widely used long-corridor-shaped surge tanks (LCSSTs), when the tank water level decreases rapidly during hydraulic transients. These ATVVs may jeopardize the operation safety of the hydropower stations and should be avoided. This study elucidates the formation mechanism of the ATVVs and proposes some simple measures to eliminate them. The 3-D CFD model for predicting the ATVVs is validated first by physical model tests in a model tailrace LCSST, and then the formation mechanism is analyzed based on the numerical results. It is shown that the main influence factor for the ATVVs is the critical submergence, which can be reduced by minimizing the velocity circulation around the throttled orifices. Two practical ATVV elimination measures through suppressing the velocity circulation are compared and verified, and the optimized one is recommended.

Hydraulic machinery mainly includes turbine and pump, which is closely related to national economy and people's livelihood involving aerospace industry, marine engineering, hydropower engineering, petroleum industry, chemical industry, mining industry, biomedical engineering, environmental engineering, agricultural water-soil engineering, etc.. The internal flow of hydraulic machinery is extremely complex, and its characteristics can be summarized as high Reynolds number, multi-scales, inhomogeneous and vortex-dominant unsteady turbulence which interact with the rotating dynamic boundary (rotor blade). Based on the analysis of the internal flow characteristics of hydraulic machinery, the author and his research team successively proposed a rotation correction model, a curvature corrected filter-based model, a scalable detached eddy simulation method, and a non-linear hybrid RANS/LES turbulence model to capture unsteady flow structures and then predict hydraulic performance and dynamic characteristics more accurately. According to the analysis on the internal flow, the corresponding flow control measures were put forward. It was verified by experiments that these methods could significantly improve the hydraulic performance, anti-cavitation performance and dynamic characteristics (pressure pulsation and vibration) of hydraulic machinery in a certain range of operating conditions. In addition, the mechanism how flow control measures influence internal flow was analyzed in depth, aiming at finding a feasible and effective way to improve hydraulic performance, anti-cavitation performance and dynamic characteristics of hydraulic machinery.

With the quasi-static analysis method, the terminal floating state of a damaged ship is usually evaluated for the risk assessment. But this is not enough since the ship has the possibility to lose its stability during the transient flooding process. Therefore, an enhanced smoothed particle hydrodynamics (SPH) model is applied in this paper to investigate the response of a simplified cabin model under the condition of the transient water flooding. The enhanced SPH model is presented firstly including the governing equations, the diffusive terms, the boundary implementations and then an algorithm regarding the coupling motions of six degrees of freedom (6-DOF) between the structure and the fluid is described. In the numerical results, a non-damaged cabin floating under the rest condition is simulated. It is shown that a stable floating state can be reached and maintained by using the present SPH scheme. After that, three-dimensional (3-D) test cases of the damaged cabin with a hole at different locations are simulated. A series of model tests are also carried out for the validation. Fairly good agreements are achieved between the numerical results and the experimental data. Relevant conclusions are drawn with respect to the mechanism of the responses of the damaged cabin model under water flooding conditions.

An incompressible smoothed particle hydrodynamics (ISPH) model has been developed to investigate the flow-like landslide phenomena The landslide mass is idealized as rigid and perfectly-plastic material with a constant density. Unlike the widely-used explicit smoothed particle hydrodynamics (SPH) models for landslides, the Chorin's projection method is used herein to implicitly solve the normal stress via a pressure Poisson equation, leading to a realistic distribution of instantaneous stress fields free of spurious fluctuations. The capability of the model is demonstrated through three case studies, including the idealized granular flow, landslide interaction with a rigid barrier, and a historical cut-slope landslide in Hong Kong.

Runoff coefficient is an important parameter for the decision support of urban stormwater management. However, factors like comprehensive land-use type, variable spatial elevation, dynamic rainfall and groundwater elevation, make the direct estimation of runoff coefficient difficult. This paper presented a novel method to estimate the urban runoff coefficient using the inverse method, where observed time-series catchment outfall flow volume was employed as input for the water balance model and runoff coefficients of different catchments were treated as unknown parameters. A developed constrained minimization objective function was combined to solve the model and minimized error between observed and modeled outfall flow is satisfactory for the presenting of a set of runoff coefficients. Estimated runoff coefficients for the urban catchments in Shanghai downtown area demonstrated that practice of low impact design could play an important role in reducing the urban runoff.

The boundary element method (BEM) is a main method for analyzing the interactions between the waves and the marine structures. As with the BEM, a set of linear equations are generated with a full matrix, the required calculations and storage increase rapidly with the increase of the structure scale. Thus, an accelerated method with a low storage is desirable for the wave interaction with a very large structure. A systematic review is given in this paper for the BEM for solving the problem of the wave interaction with a large scale structure. Various integral equations are derived based on different Green functions, the advantages and disadvantages of different discretization schemes of the integral equations by the constant panels, the higher order elements, and the spline functions are discussed. For the higher order element discretization method, the special concerns are given to the numerical calculations of the single-layer potential, the double layer potential and the solid angle coefficients. For a large scale computation problem such as the wave interaction with a very large structure or a large number of bodies, the BEMs with the FMM and pFFT accelerations are discussed, respectively, including the principles of the FMM and the pFFT, and their implementations in various integral equations with different Green functions. Finally, some potential applications of the acceleration methods for problems with large scale computations in the ocean and coastal engineering are introduced.