首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
A three-dimensional k-ε-Ap two-fluid turbulence model is proposed to study liquid-particle two-phase flow and bed deformation. By solving coupled liquid-phase and solid-phase governing equations in a finite-volume method, the model can calculate the movement of both water and sediment. The model was validated by water-sediment transport in a 180° channel bend with a movable bed. The validation concerns two-phase time-averaged velocities, bed deformation, water depth, depth-averaged streamwise velocity, cross-stream bed profiles, and two-phase secondary flow velocity vectors. The agreement between numerical results and experimental results was generally good. The comparisons of the numerical results of different models show that the three-dimensional k-ε-Ap two-fluid turbulence model has a relatively higher accuracy than one-fluid model.  相似文献   

2.
Wave breaking and wave runup/rundown have a major influence on nearshore hydrodynamics, morphodynamics and beach evolution. In the case of wave breaking, there is significant mixing of air and water at the wave crest, along with relatively high kinetic energy, so prediction of the free surface is complicated. Most hydrodynamic studies of surf and swash zone are derived from single-phase flow, in which the role of air is ignored. Two-phase flow modeling, consisting of both phases of water and air, may be a good alternative numerical modeling approach for simulating nearshore hydrodynamics and, consequently, sediment transport. A two-phase flow tool can compute more realistically the shape of the free surface, while the effects of air are accounted for. This paper used models based on two-dimensional, two-phase Reynolds-averaged Navier–Stokes equations, the volume-of-fluid surface capturing technique and different turbulence closure models, i.e., kε, kω and re-normalized group (RNG). Our numerical results were compared with the available experimental data. Comparison of the employed method with a model not utilizing a two-phase flow modeling demonstrates that including the air phase leads to improvement in simulation of wave characteristics, especially in the vicinity of the breaking point. The numerical results revealed that the RNG turbulence model yielded better predictions of nearshore zone hydrodynamics, although the kε model also gave satisfactory predictions. The model provides new insights for the wave, turbulence and means flow structure in the surf and swash zones.  相似文献   

3.
A two-dimensional numerical model was presented for the simulation of wave breaking, runup and turbulence in the surf and swash zones. The main components of the model are the Reynolds-Averaged Navier–Stokes equations describing the average motion of a turbulent flow, a kε turbulence closure model describing the transformation and dissipation processes of turbulence and a volume of fluid technique for tracking the free surface motion. Nearshore wave evolution on a sloping bed, the velocity field and other wave characteristics were investigated. First, the results of the model were compared with experimental results for different surf zone hydrodynamic conditions. Spilling and plunging breakers were simulated and the numerical model investigated for different wave parameters. The turbulence field was also considered and the spatial and time-dependent variations of turbulence parameters were discussed. In the next stage of the study, numerical results were compared with two sets of experimental data in the swash zone. Generally, there is good agreement except for turbulence predictions near the breaking point where the model does not represent well the physical processes. On the other hand, turbulence predictions were found to be excellent for the swash zone. The model provides a precise and efficient tool for the simulation of the flow field and wave transformations in the nearshore, especially in the swash zone. The numerical model can simulate the surface elevation of the vertical shoreline excursion on sloping beaches, while swash–swash interactions within the swash zone are accounted for.  相似文献   

4.
The mixing characteristics of dredged sediments of variable size discharged into cross-flow are studied by an Eulerian-Lagrangian method. A three-dimensional (3D) numerical model has been developed by using the modified k-ε parameterization for the turbulence in fluid phase/water and a Lagrangian method for the solid phase/sediments. In the model the wake turbulence induced by sediments has been included as additional source and sink terms in the k-ε model; and the trajectories of the sediments are tracked by the Lagrangian method in which the sediment drift velocities in cross-flow are computed by a multiphase particle-in-cell (MP-PIC) method and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumped sediment cloud is governed by the total buoyancy on the cloud, the drag force on each particle and velocity of cross-flow. The cross-flow destroys more or less the double vortices occurred in stagnant ambience and dominates the longitudinal movement of sediment cloud. The computed results suggest satisfactory agreement by comparison with the experimental results of laboratory.  相似文献   

5.
In this paper, we address the transport of multi-disperse suspended sediment mixtures in open channels, via the use of the two-fluid model. To that end, we extend previously developed frameworks for the dilute and non-dilute transport of suspended sediment. Within the scope of the Reynolds-averaged Navier-Stokes (RANS) equations, these modeling frameworks comprise mass and momentum equations for both phases (water and sediment). Here, we calculate the distribution of total volumetric concentration of sediment using two approaches: (1) by considering the mixture as represented by a single size; we call this approach Partial two-fluid model for uniform sediments (PTFMU); and (2) by combining the volumetric concentration of the sediment corresponding to several particle size classes; we call this approach Partial two-fluid model for non-uniform sediments (PTFMNU). In the second approach, we propose a methodology for the computation of the overall velocity of the disperse phase as a function of the velocities of each size class. k-ε type closures to account for the turbulence in the carrier phase (water) are applied. We also consider the coupling between the two phases through the drag force. Velocities of the carrier and disperse phases, and concentrations for each sediment class size are numerically solved by integrating the differential equations over control volumes. In order to validate our models, we compare numerical results to experimental data of Einstein and Chien [H.A. Einstein, N. Chien, Effects of heavy sediment concentration near the bed on velocity and sediment distribution, MRD sediment series report, University of California, Berkley, 1955] and Taggart et al. [W.C. Taggart, C.A. Yermoli, S. Montes, A. Ippen, Effects of sediment size and gradation on concentration profiles for turbulent flow, Massachusetts Institute of Technology, 1972]. Results of mean velocity of the carrier phase are in close agreement with the experimental data. For the prediction of sediment concentrations, we observe that there is a difference in the results using the two approaches mentioned above. We additionally obtain values of the Schmidt number needed to improve the agreement between predictions of the distribution of suspended sediment and the experimental data, and discuss the effect of sediment size and increasing sediment concentration on the values of the Schmidt number.  相似文献   

6.
This numerical investigation was carried out to advance mechanistic understanding of sediment transport under sheet flow conditions. An Euler–Euler coupled two-phase flow model was developed to simulate fluid–sediment oscillatory sheet flow. Since the concentration of sediment particles is high in such flows, the kinematics of the fluid and sediment phases are strongly coupled. This model includes interaction forces, intergranular stresses and turbulent stress closure. Each phase was modeled via the Reynolds-Averaged Navier–Stokes equations, with interphase momentum conservation accounting for the interaction between the phases. The generation and transformation of turbulence was modeled using the two-equation k–εkε turbulence model. Concentration and sediment flux profiles were compared with experimental data for sheet flow conditions considering both symmetric and asymmetric oscillatory flows. Sediment and fluid velocity variations, concentration profiles, sediment flux and turbulence parameters of wave-generated sheet flow were studied numerically with a focus on sediment transport characteristics. In all applications, the model predictions compared well with the experimental data. Unlike previous investigations in which the flow is driven by a horizontal pressure gradient, the present model solves the Navier–Stokes equations under propagating waves. The model’s ability to predict sediment transport under oscillatory sheet flow conditions underscores its potential for understanding the evolution of beach morphology.  相似文献   

7.
A key problem in computational fluid dynamics (CFD) modelling of gravel‐bed rivers is the representation of multi‐scale roughness, which spans the range from grain size, through bedforms, to channel topography. These different elements of roughness do not clearly map onto a model mesh and use of simple grain‐scale roughness parameters may create numerical problems. This paper presents CFD simulations for three cases: a plane bed of fine gravel, a plane bed of fine gravel including large, widely‐spaced pebble clusters, and a plane gravel bed with smaller, more frequent, protruding elements. The plane bed of fine gravel is modelled using the conventional wall function approach. The plane bed of fine gravel including large, widely‐spaced pebble clusters is modelled using the wall function coupled with an explicit high‐resolution topographic representation of the pebble clusters. In these cases, the three‐dimensional Reynolds‐averaged continuity and Navier–Stokes equations are solved using the standard k ? ε turbulence model, and model performance is assessed by comparing predicted results with experimental data. For gravel‐bed rivers in the field, it is generally impractical to map the bed topography in sufficient detail to enable the use of an explicit high‐resolution topography. Accordingly, an alternative model based on double‐averaging is developed. Here, the flow calculations are performed by solving the three‐dimensional double‐averaged continuity and Navier‐Stokes equations with the spatially‐averaged 〈k ? ε〉 turbulence model. For the plane bed of fine gravel including large, widely‐spaced pebble clusters, the model performance is assessed by comparing the spatially‐averaged velocity with the experimental data. The case of a plane gravel bed with smaller, more frequent, protruding elements is represented by a series of idealized hypothetical cases. Here, the spatially‐averaged velocity and eddy viscosity are used to investigate the applicability of the model, compared with using the explicit high‐resolution topography. The results show the ability of the model to capture the spatially‐averaged flow field and, thus, illustrate its potential for representing flow processes in natural gravel‐bed rivers. Finally, practical data requirements for implementing such a model for a field example are given. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

8.
Suspended sediment is carried by turbulent water flows in rivers.Traditional sediment-laden flow analysis treats the suspension as a mixed liquid,and recent two-phase flow model enables separate velocity measurement of the two coupling phases.A simplified theoretical analysis was presented to discuss the differences between the two models in reporting turbulence intensity in experimental research.The turbulence intensity of the mixture is lower than the weighted average of those of the two phases in mixed-flow experiments.The mixture’s turbulence intensity becomes higher than the average of the two phases in two-phase experiments due to the presence of velocity lag.The same set of data may lead to either an underestimation or an overestimation of actual turbulence levels when different models are used.  相似文献   

9.
《国际泥沙研究》2020,35(4):395-407
A two-dimensional vertical (2DV), Eulerian two-phase model or complete two-fluid model of the free surface flow was developed to simulate water-sediment flow in a local scour hole. In the model, the complete forms of the vertical, two-dimensional, two-fluid Navier-Stokes equations were discretized using a finite volume scheme. This discretization was done based on a standard staggered grid system using a curvilinear network system in compliance with the bed boundaries and water level. At the beginning of the computational cycle, the equations governing the fluid phase were solved based on the two-step projection method with a pressure-correction technique. In the first step, the intermediate fluid velocities were obtained by solving different phases of the momentum equations of the fluid phase using the time-splitting technique. In the second step, pressure was obtained and fluid velocities were updated. In this step a simple discretization method was applied for decreasing the computational complexity. After obtaining all the fluid phase variables at a new time step, the sediment phase momentum equations were solved using the time-splitting technique and sediment velocities were obtained. Then, at the end of the computational cycle, the sediment phase mass equation was solved and the concentrations of both phases were updated. At last, the capacity of the model for simulating of the longitudinal fluid velocity and sediment concentration in a local scour hole was evaluated. Numerical results were found to be in good agreement with experimental data.  相似文献   

10.
Fluvial bedforms generate a turbulent wake that can impact suspended-sediment settling in the passing flow. This impact has implications for local suspended-sediment transport, bedform stability, and channel evolution; however, it is typically not well-considered in geomorphologic models. Our study uses a three-dimensional OpenFOAM hydrodynamic and particle-tracking model to investigate how turbulence generated from bedforms and the channel bed influences medium sand-sized particle settling, in terms of the distribution of suspended particles within the flow field and particle-settling velocities. The model resolved the effect of an engineered bedform, which altered the flow field in a manner similar to a natural dune. The modelling scenarios alternated bed morphology and the simulation of turbulence, using detached eddy simulation (DES), to differentiate the influence of bedform-generated turbulence relative to that of turbulence generated from the channel bed. The bedform generated a turbulent wake that was composed of eddies with significant anisotropic properties. The eddies and, to a lesser degree, turbulence arising from velocity shear at the bed substantially reduced settling velocities relative to the settling velocities predicted in the absence of turbulence. The eddies tended to advect sediment particles in their primary direction, diffuse particles throughout the flow column, and reduced settling likely due to production of a positively skewed vertical-velocity fluctuation distribution. Study results suggest that the bedform wake has a significant impact on particle-settling behaviour (up to a 50% reduction in settling velocity) at a scale capable of modulating local suspended transport rates and bedform dynamics. © 2020 John Wiley & Sons, Ltd.  相似文献   

11.
12.
In this study an incompressible smoothed particle hydrodynamics (ISPH) approach coupled with the sediment erosion model is developed to investigate the sediment bed scour and grain movement under the dam break flows. Two-phase formulations are used in the ISPH numerical algorithms to examine the free surface and bed evolution profiles, in which the entrained sediments are treated as a different fluid component as compared with the water. The sediment bed erosion model is based on the concept of pick-up flow velocity and the sediment is initiated when the local flow velocity exceeds a critical value. The proposed model is used to reproduce the sediment erosion and follow-on entrainment process under an instantaneous dam break flow and the results are compared with those from the weakly compressible moving particle semi-implicit (WCMPS) method as well as the experimental data. It has been demonstrated that the two-phase ISPH model performed well with the experimental data. The study shows that the ISPH modelling approach can accurately predict the dynamic sediment scouring process without the need to use empirical sediment transport formulas.  相似文献   

13.
A 3D Lagrangian model of the saltation of solid spherical particles on the bed of an open channel flow, accounting for turbulence-induced mechanisms, is proposed and employed as the key tool of the study. The differences between conventional 2D models and a proposed 3D saltation model are discussed and the advantages of the 3D model are highlighted. Particularly, the 3D model includes a special procedure allowing generation of 3D flow velocity fields. This procedure is based on the assumption that the spectra of streamwise, vertical and transverse velocity components are known at any distance from the bed. The 3D model was used to identify and quantify effects of turbulence on particle entrainment and saltation. The analysis of particle trajectories focused on their diffusive nature, clarifying: (i) the effect of particle mobility parameter; (ii) the effect of bed topography; and (iii) the effect of turbulence. Specifically, the results of numerical simulations describing the abovementioned effects on the change in time of the variance are presented. In addition, the change in time of the skewness and kurtosis, which are likely to reflect the turbulence influence on the spread of particles, are also shown. Two different diffusion regimes (local and intermediate) for each of the investigated flow conditions are confidently identified.  相似文献   

14.
An in‐house fully three‐dimensional general‐purpose finite element model is applied to solve the hydrodynamic structure in a periodic Kinoshita‐generated meandering channel. The numerical model solves the incompressible Reynolds‐averaged Navier–Stokes equations for mass and momentum, while solving the k ? ε equations for turbulence. The free surface is described by the rigid‐lid approximation (using measured water surface data) for flat (smooth‐bed) and self‐formed (rough‐bed) conditions. The model results are compared against experimental measurements in the ‘Kinoshita channel’, where three‐dimensional flow velocities and turbulence parameters were measured. This validation was carried out for the upstream‐valley meander bend orientation under smooth (flat bed) conditions. After validation, several simulations were carried out to predict the hydrodynamics in conditions where either it was not possible to perform measurements (e.g. applicability of the laboratory acoustic instruments) and to extrapolate the model to other planform configurations. For the flat smooth‐bed case, a symmetric (no skewness) planform configuration was modeled and compared to the upstream‐skewed case. For the self‐formed rough‐bed case, prediction of the hydrodynamics during the progression of bedforms was performed. It appears that the presence of bedforms on a bend has the following effects: (i) the natural secondary flow of the bend is disrupted by the presence of the bedforms, thus depending on the location of the dune, secondary flows might differ completely from the traditional orientation; (ii) an increment on both the bed and bank shear stresses is induced, having as much as 50% more fluvial erosion, and thus a potential increment on the migration rate of the bend. Implications on sediment transport and bend morphodynamics are also discussed in the paper. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

15.
A three-dimensional(3D) numerical model of unstable turbidity currents is developed based on the mechanism of sediment transport and turbulence theory.In this model,numerical simulation of turbidity currents without subsequent supply of muddy water was conducted using the same parameters as were used in the flume experiments.The evolution process of turbidity currents of completely losing supplies observed in the experiment was simulated by the model;validation of the numerical model and the algorithm was conducted.If momentarily interrupted process is regarded as a special case of the gradually interrupted,based on the preceding numerical simulation validity,it is feasible to simulate the motion law of turbidity currents under losing gradually supplies.By this method in this article,the characteristic of sediment-laden flow of losing gradually supplies was obtained,as well as its relationship between front velocity and sediment concentration.  相似文献   

16.
首先以Lorenz混沌方程产生的非线性时间序列为例,讨论了在不同时间序列长度下各种延迟时间算法对噪声的适用性.研究发现,采用C_C算法计算延迟时间的鲁棒性强.在此基础上,给出了垂直上升管中气水两相流电导波动信号混沌表征结果,发现在较低水相表观速度时,随着气相表观速度增加,泡状流及混状流动力学特性变得愈加复杂,而段塞流动力学特性受液相表观速度影响较大;在较高水相表观速度时,随着气相表观速度增加,当流型从泡状流向段塞流转变时,气液两相流动力学特性变得相对简单.但是,由于受液相湍流作用影响,段塞流的动力学特性表现出了涨落现象,呈现不稳定性,当流型从段塞流向混状流转变时,气液两相流动力学特性则变得愈加复杂.研究结果表明:基于电导波动信号的混沌分析可以较好地表征气液两相流流型变化,是理解流型转变机理及其动力学演变特性的有用工具.  相似文献   

17.
A three-dimensional finite-element model is used to investigate the tidal flow around Rattray Island, Great Barrier Reef, Australia. Field measurements and visual observations show both stable eddies developing at rising and falling tide in the wake of the island. The water turbidity suggests intense upwelling able to carry bed sediments upwards. Based on previous numerical studies, it remains unclear at this point whether the most intense upwelling occurs near the centre of the eddies or off the island's tips, closer to the island. All these studies resorted to a very simple turbulence closure, with a zero-equation model whereby the coefficient of vertical viscosity is computed via an algebraic expression. In this work, we aim at studying the influence of the turbulence closure on model results, with emphasis on the prediction of vertical motions. The Mellor and Yamada level 2.5 closure scheme is used and an increase in the intensity of vertical transport is observed. This increase is partly explained by the fact that the Mellor and Yamada model takes into account the hysteresis effect in the time variation of turbulence variables. The influence of the advection of turbulence variables is estimated to be negligible. By a better representation of transient coastal phenomena, the Mellor and Yamada level 2.5 turbulence closure improves the model to a significant degree.  相似文献   

18.
The permeability of river beds is an important control on hyporheic flow and the movement of fine sediment and solutes into and out of the bed. However, relatively little is known about the effect of bed permeability on overlying near‐bed flow dynamics, and thus on fluid advection at the sediment–water interface. This study provides the first quantification of this effect for water‐worked gravel beds. Laboratory experiments in a recirculating flume revealed that flows over permeable beds exhibit fundamental differences compared with flows over impermeable beds of the same topography. The turbulence over permeable beds is less intense, more organised and more efficient at momentum transfer because eddies are more coherent. Furthermore, turbulent kinetic energy is lower, meaning that less energy is extracted from the mean flow by this turbulence. Consequently, the double‐averaged velocity is higher and the bulk flow resistance is lower over permeable beds, and there is a difference in how momentum is conveyed from the overlying flow to the bed surface. The main implications of these results are three‐fold. First, local pressure gradients, and therefore rates of material transport, across the sediment–water interface are likely to differ between impermeable and permeable beds. Second, near‐bed and hyporheic flows are unlikely to be adequately predicted by numerical models that represent the bed as an impermeable boundary. Third, more sophisticated flow resistance models are required for coarse‐grained rivers that consider not only the bed surface but also the underlying permeable structure. Overall, our results suggest that the effects of bed permeability have critical implications for hyporheic exchange, fluvial sediment dynamics and benthic habitat availability. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.  相似文献   

19.
20.
Sediment transport models require appropriate representation of near-bed processes. We aim here to explore the parameterizations of bed shear stress, bed load transport rate and near-bed sediment erosion rate under the sheet flow regime. To that end, we employ a one-dimensional two-phase sheet flow model which is able to resolve the intrawave boundary layer and sediment dynamics at a length scale on the order of the sediment grain. We have conducted 79 numerical simulations to cover a range of collinear wave and current conditions and sediment diameters in the range 210–460 μmμm. The numerical results confirm that the intrawave bed shear stress leads the free stream velocity, and we assess an explicit expression relating the phase lead to the maximum velocity, wave period and bed roughness. The numerical sheet flow model is also used to provide estimates for the bed load transport rate and to inspect the near-bed sediment erosion. A common bed load transport rate formulation and two typical reference concentration approaches are assessed. A dependence of the bed load transport rate on the sediment grain diameter is observed and parameterized. Finally, the intrawave near-bed vertical sediment flux is further investigated and related to the time derivative of the bed shear stress.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号