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

In this paper a new finite-volume non-hydrostatic and shock-capturing three-dimensional model for the simulation of wave-structure interaction and hydrodynamic phenomena (wave refraction, diffraction, shoaling and breaking) is proposed. The model is based on an integral formulation of the Navier-Stokes equations which are solved on a time dependent coordinate system: a coordinate transformation maps the varying coordinates in the physical domain to a uniform transformed space. The equations of motion are discretized by means of a finite-volume shock-capturing numerical procedure based on high order WENO reconstructions. The solution procedure for the equations of motion uses a third order accurate Runge-Kutta (SSPRK) fractional-step method and applies a pressure corrector formulation in order to obtain a divergence-free velocity field at each stage. The proposed model is validated against several benchmark test cases.

The spatial relationship between the energy dissipation slabs and the vortex tubes is investigated based on the direct numerical simulation (DNS) of the channel flow. The spatial distance between these two structures is found to be slightly greater than the vortex radius. Comparison of the core areas of the vortex tubes and the dissipation slabs gives a mean ratio of 0.16 for the mean swirling strength and that of 2.89 for the mean dissipation rate. These results verify that in the channel flow the slabs of intense dissipation and the vortex tubes do not coincide in space. Rather they appear in pairs offset with a mean separation of approximately 10η.

The flow pattern in the confluent meander bend channel under the conditions of different discharge ratios and junction angles is numerically simulated by means of the large eddy simulation (LES), and the characteristics of the flow separation zone are analyzed. Numerical results are well validated by experimental data with a good agreement. Analysis of the vertical confinement shows that the turbulence within the separation zone can be characterized as quasi-2-D. Details of the separation zone characteristics are revealed as shown by mean velocity isolines. According to the analysis of numerical results, the length and the width of the separation zone generally increase with the increase of the discharge ratio and the junction angle. However, the width of the separation zone keeps substantially constant when the junction angle increases from to 60° to 90°. The dimensionless shape of the separation zone is nearly the same for three discharge ratios and three junction angles. The formulas of the relative width and the relative length of the separation zone are obtained by means of the polynomial fit method.

Compared with conventional channels, experiments of microchannel often exhibit some controversial findings and sometimes even opposite trends, most notably the effects of the Reynolds number and the scaled channel height on the Poiseuille number. The experimental method has still been constrained by two key facts, firstly the current ability to machine microstructures and secondly the limitation of measurement of parameters related to the Poiseuille number. As a consequence, numerical method was adopted in this study in order to analyze a flow in two-dimensional rectangular microchannels using water as working fluid. Results are obtained by the solution of the steady laminar incompressible Navier-Stokes equations using control volume finite element method (CVFEM) without pressure correction. The computation was made for channel height ranging from 50 μm to 4.58 μm and Reynolds number varying from 0.4 to 1 600. The effect of Reynolds number and channel heights on flow characteristics was investigated. The results showed that the Poiseuille numbers agree fairly well with the experimental measurements proving that there is no scale effect at small channel height. This scaling effect has been confirmed by two additional simulations being carried out at channel heights of 2.5 μm and 0.5 μm, respectively and the range of Reynolds number was extended from 0.01 up to 1 600. This study confirm that the conventional analysis approach can be employed with confidence for predicting flow behavior in microchannels when coupled with carefully matched entrance and boundary conditions in the dimensional range considered here.

In this paper, we investigate the verification and validation (V&V) procedures for the URANS simulations of the turbulent cavitating flow around a Clark-Y hydrofoil. The main focus is on the feasibility of various Richardson extrapolation-based uncertainty estimators in the cavitating flow simulation. The unsteady cavitating flow is simulated by a density corrected model (DCM) coupled with the Zwart cavitation model. The estimated uncertainty is used to evaluate the applicability of various uncertainty estimation methods for the cavitating flow simulation. It is shown that the preferred uncertainty estimators include the modified Factor of Safety (FS1), the Factor of Safety (FS) and the Grid Convergence Index (GCI). The distribution of the area without achieving the validation at the V_v level shows a strong relationship with the cavitation. Further analysis indicates that the predicted velocity distributions, the transient cavitation patterns and the effects of the vortex stretching are highly influenced by the mesh resolution.

A hybrid approach coupled with a surface panel method for the propeller and a Reynolds averaged Navier-Stokes (RANS) model for the hull with the propeller body forces are presented for predicting the self-propulsion performance and the effective wake field of underwater vehicles. To achieve a high accuracy and simplicity, a radial basis function (RBF) based approach is proposed for mapping the force field from the blade surface panels to the RANS model. The effective wake field is evaluated in two ways, i.e., by extrapolation from the flat planes upstream of the propeller disk, and by direct computation in a curved surface upstream of and parallel to the blade leading edges. The hull-propeller system of a real propeller geometry is further simulated with the sliding mesh model to numerically verify the hybrid approach. Numerical simulations are conducted for the fully appended SUBOFF submarine model. The high accuracy of the RBF-based interpolation scheme is confirmed, and the effective wake fraction predicted by the hybrid approach is found consistent with that obtained by the sliding mesh model. The effective wake fractions predicted by the two methods are, respectively, 4.6% and 3% larger than the nominal one.

The cavitation in a mechanical heart valve (MHV) is a serious concern. In most of the investigations of the MHV cavitation in vitro, the tap water, the distilled water, or the glycerin are used as the test liquids, instead of the real blood. Therefore, the effects of the liquid properties on the cavitation can not be well revealed. In this paper, the cavitation erosion in the porcine bloods is experimentally investigated as well as in the tap water and the distilled water by means of a vibratory apparatus. The results show that the blood produces a weaker intensity of the cavitation erosion than the tap water or the distilled water. The cavitation erosion decreases with the decrease of the liquid temperature or with the increase of the concentration of the blood, especially with the increase of the liquid viscosity. It is the viscosity that could be a major dominant factor affecting this erosion. The temperature or the concentration of the blood changes the viscosity, and in turns changes the intensity of the cavitation erosion.

This paper investigates the effect of initial volume fraction on the runout characteristics of collapse of granular columns on slopes in fluid. 2-D sub-grain scale numerical simulations are performed to understand the flow dynamics of granular collapse in fluid. The discrete element method (DEM) technique is coupled with the lattice Boltzmann method (LBM), for fluid-grain interactions, to understand the evolution of submerged granular flows. The fluid phase is simulated using multiple-relaxation-time LBM (LBM-MRT) for numerical stability. In order to simulate interconnected pore space in 2-D, a reduction in the radius of the grains (hydrodynamic radius) is assumed during LBM computations. The collapse of granular column in fluid is compared with the dry cases to understand the effect of fluid on the runout behaviour. A parametric analysis is performed to assess the influence of the granular characteristics (initial packing) on the evolution of flow and run-out distances for slope angles of 0°, 2.5°, 5° and 7.5°. The granular flow dynamics is investigated by analysing the effect of hydroplaning, water entrainment and viscous drag on the granular mass. The mechanism of energy dissipation, shape of the flow front, water entrainment and evolution of packing density is used to explain the difference in the flow characteristics of loose and dense granular column collapse in fluid.

The turbulent flow over a channel bed roughened by three layers of closely packed spheres with a Reynolds number of Re = 15000 is investigated using the large eddy simulation (LES) and the double-averaging (DA) method. The DA velocity is compared with the results of the corresponding laboratory experiments to validate the LES results. The existence of the types of vortex structures is demonstrated by the Q - criterion above the permeable bed. The turbulent kinetic energy (TKE) fluxes and budget are quantified and discussed. The results show that the TKE fluxes are directed downward and downstream near the virtual bed level. In the TKE budget, the form-induced diffusion rate is significant in the vicinity of the crest bed level, and the TKE production rate and the dissipation rate attain their peaks at the crest bed level and decrease sharply below it.