共查询到20条相似文献，搜索用时 15 毫秒
1.
Numerical simulation of two-phase flow for sediment transport in the inner-surf and swash zones 总被引：1，自引：0，他引：1
R. Bakhtyar D.A. Barry A. Yeganeh-Bakhtiary L. Li J.-Y. Parlange G.C. Sander 《Advances in water resources》2010
A two-dimensional two-phase flow framework for fluid–sediment flow simulation in the surf and swash zones was described. Propagation, breaking, uprush and backwash of waves on sloping beaches were studied numerically with an emphasis on fluid hydrodynamics and sediment transport characteristics. The model includes interactive fluid–solid forces and intergranular stresses in the moving sediment layer. In the Euler–Euler approach adopted, two phases were defined using the Navier–Stokes equations with interphase coupling for momentum conservation. The k–ε closure model and volume of fluid approach were used to describe the turbulence and tracking of the free surface, respectively. Numerical simulations explored incident wave conditions, specifically spilling and plunging breakers, on both dissipative and intermediate beaches. It was found that the spatial variation of sediment concentration in the swash zone is asymmetric, while the temporal behavior is characterized by maximum sediment concentrations at the start and end of the swash cycle. The numerical results also indicated that the maximum turbulent kinetic energy and sediment flux occurs near the wave-breaking point. These predictions are in general agreement with previous observations, while the model describes the fluid and sediment phase characteristics in much more detail than existing measurements. With direct quantifications of velocity, turbulent kinetic energy, sediment concentration and flux, the model provides a useful approach to improve mechanistic understanding of hydrodynamic and sediment transport in the nearshore zone. 相似文献
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.
Turbulent flow in a meandering channel is computed with two Computational Fluid Dynamics (CFD) codes solving the Navier–Stokes equations by employing different turbulence closure approaches. The first CFD code solves the steady Reynolds-Averaged Navier–Stokes equations (RANS) using an isotropic turbulence closure. The second code is based on the concept of Large Eddy Simulation (LES). LES resolves the large-scale turbulence structures in the flow and is known to outperform RANS models in flows in which large-scale structures dominate the statistics. The results obtained from the two codes are compared with experimental data from a physical model study. Both, LES and RANS simulation, predict the primary helical flow pattern in the meander as well as the occurrence of an outer-bank secondary cell. Computed primary as well as secondary flow velocities are in reasonably good agreement with experimental data. Evidence is given that the outer-bank secondary cell in a meander bend is the residual of the main secondary cell of the previous bend. However, the RANS code, regardless of the turbulence model employed, overpredicts the size and strength of the outer-bank secondary cell. Furthermore, only LES is able to uphold the outer-bank second secondary cell beyond the bend apex until the exit of the bend as turbulence anisotropy contributes to its persistence. The presence of multiple secondary cells has important consequences for the distribution of shear stresses along the wetted perimeter of the channel, and thereby the sediment transport in meandering channels. Consequently, even though LES is expected to compute the bed-shear stresses along the wetted perimeter of the channel with a higher degree of accuracy than the RANS model, comparisons between LES and RANS computed wall shear stresses agree well. These findings are useful for practitioners who need to rely on RANS model predictions of the flow in meandering channels at field scale. 相似文献
5.
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. 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. 相似文献
6.
This study numerically investigates effects of cutting riparian vegetation on flow characteristics by using a two-dimensional numerical model. The numerical model is based on depth-averaging the time- and volume-averaged Navier–Stokes equation with turbulent effects determined by the standard k–ε turbulence model. Drag forces exerted by the flow on vegetation are considered by adding source terms into momentum equations. In a rectangular channel and compound channel with vegetation along one side, numerical predictions show are in good agreement with those of previous studies. Five cutting scenarios, including the original, cutting along the main channel side, cutting along the bank side, alternative cutting, and reducing vegetative density, are analyzed in this study. The influences of the cutting scenarios on hydrodynamic behaviors are evaluated via numerical simulations. Simulation results suggest that cutting along the main channel side is the most effective scenario for reducing water depth and flow velocities. 相似文献
7.
George S. Constantinescu Witold F. KrajewskiCelalettin E. Ozdemir Talia Tokyay 《Advances in water resources》2007
Wind is responsible for systematic errors that affect rain gauge measurements. The authors investigate the use of computational fluid dynamics (CFD) to calculate airflow around rain gauges by applying a high-resolution large eddy simulation (LES) model to determine the flow fields around a measuring system of two rain gauges. The simulated air flow field is characterized by the presence of massive separation which induces the formation and shedding of highly unsteady eddies in the detached shear layers and wakes. Parts of these detached structures occur over the orifice of the rain gauges and may substantially affect the dynamics of the raindrops in this critical region. Non-dissipative LES methods used with fine enough meshes can successfully predict these eddies and their associated fluctuations. The authors compare statistics from LES with steady-state Reynolds averaged Navier–Stokes (RANS) simulations using the k–ε and shear stress transport k–ω turbulence models. They find that both RANS and LES models predict similar mean velocity distributions around the rain gauges. However, they determine the distribution of the resolved turbulent kinetic energy (TKE) to be strongly dependent on the RANS model used. Neither RANS model predictions of TKE are close to those of LES. The authors conclude that the failure of RANS to predict TKE is an important limitation, as TKE is needed to scale the local velocity fluctuations in stochastic models used to calculate the motion of raindrops in the flow field. 相似文献
8.
Wave-induced, steep vortex ripples are ubiquitous features in shallow coastal seas and it is therefore important to fully understand and model the sediment transport processes that occur over them. To this end, two two-dimensional vertical (2DV) models have been critically tested against detailed velocity and sediment concentration measurements above mobile ripples in regular asymmetric oscillatory flow. The two models are a k–ω turbulence-closure model and a discrete-vortex, particle-tracking (DVPT) model, while the data are obtained in the Aberdeen oscillatory flow tunnel (AOFT). The models and the data demonstrate that the time-dependent velocity and suspended sediment concentration above the ripple are dominated by the generation of lee-side vortices and their subsequent ejection at flow reversal. The DVPT model predicts the positions and strengths of the vortices reasonably well, but tends to overpredict the velocity close to the ripple surface. The k–ω model, on the other hand, underpredicts the height to which the vortices are lifted, but is better able to predict the velocity close to the bed. In terms of the cycle- and ripple-averaged horizontal velocity, both models are able to reproduce the observed offshore flow close to and below the ripple crest and the DVPT model is able to produce the onshore flow higher up. In the vicinity of the vortices, the DVPT model better represents the concentration (because of its better prediction of vorticity). The k–ω model, on the other hand, better represents the concentration close to the ripple surface and higher up in the flow (because of the better representation of the near-bed flow and background turbulence). The measured and predicted cycle- and ripple-averaged suspended sediment concentrations are in reasonable agreement and demonstrate the expected region of exponential decay. The models are able to reproduce the observed offshore cycle- and ripple-averaged suspended sediment flux from the ripple troughs upwards, and as a result, produce net offshore suspended sediment transport rates that are in reasonable agreement. The net measured offshore suspended transport rate, based on the integration of fluxes, was found to be consistent with the total net offshore transport measured in the tunnel as a whole once the onshore transport resulting from ripple migration was taken into account, as would be expected. This demonstrates the importance of models being able to predict ripple-migration rates. However, at present neither of the models is able to do so. 相似文献
9.
Jorge D. Abad Christian E. Frias Gustavo C. Buscaglia Marcelo H. Garcia 《地球表面变化过程与地形》2013,38(13):1612-1622
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. 相似文献
10.
This paper reviews a model, developed by Shiono and Knight [Shiono K, Knight DW. Two-dimensional analytical solution for a compound channel. In: Proceedings of the 3rd international symposium on refined flow modelling and turbulence measurements, Tokyo, Japan, July 1988. p. 503–10; Shiono K, Knight DW. Turbulent open channel flows with variable depth across the channel. J Fluid Mech 1991;222:617–46 [231:693]], which yields analytical solutions to the depth-integrated Navier–Stokes equations, and includes the effects of bed friction, lateral turbulence and secondary flows. Some issues about the original model developed by Shiono and Knight (1988, 1991) are highlighted and discussed. Based on the experimental data concerning the secondary flow, two assumptions are proposed to describe the contribution of the streamwise vorticity to the flow. Two new analytical solutions are compared with the conventional solution for three simple channel shapes and one trapezoidal compound channel to highlight their differences and the importance of the secondary flow and planform vorticity term. Comparison of the analytical results with the experimental data shows that the general SKM predicts the lateral distributions of depth-averaged velocity well. 相似文献
11.
Sheet flows occur widely in natural free-surface flows including rivers in flood, tidal estuaries and coastal waters in storm conditions when bed shear stress becomes sufficiently high. Particle volumetric concentration in sheet flows normally follows a linear distribution with the Rouse [Rouse H. Modern conceptions of the mechanics of fluid turbulence. Trans ASCE, 1937;102:463–543] distribution applicable in the dilute water column above the sheet-flow layer. However, a well-verified formula for determining particle velocity distribution in sheet flows is still lacking. This paper presents formulas to describe the particle velocity profile in steady or oscillatory sheet flows. They compare well with measured data. In particular, the novel formula for determining the particle velocity at the top of bedload–sediment-dominated sublayer in sheet flows is also well verified with measured data. 相似文献
12.
13.
Issues associated with microwave link rainfall estimation such as the effects of spatial and temporal variation in rain, the nonlinearity of R–k relations, temporal sampling, power resolution, and wet antenna attenuation are investigated using more than 1.5 years of data from a high-resolution X-band weather radar. Microwave link signals are generated for different link frequencies and lengths from these radar data, so that retrieved path-averaged rainfall intensities can be compared to true path-averaged values. Results of these simulations can be linked to the space–time structure of rain. A frequency-dependent relation between the rainfall intensity at an antenna and the attenuation caused by its wetting is derived using microwave link and rain gauge data. It is shown that if the correct temporal sampling strategy is chosen, the effects of the degradation of power resolution and of wet antenna attenuation (if a correction is applied) are minor (i.e., MBE and bias-corrected RMSE are >−20% and <20% of the mean rainfall intensity, respectively) for link frequencies and lengths above ∼20 GHz and ∼2 km, respectively. 相似文献
14.
Ship-generated waves and return currents are capable of re-suspending significant quantities of bottom and bank sediments.However,most of the previous studies done on the subject do not show how and where sediment is re-suspended by the wakes and the directions of net transport.In this paper,a 3D numerical model based on hydro-sedimentary coupling is presented to search the relationship between the sediment movement,and the pattern of ship-generated waves around and far away from the vessel and the return currents around the ships.The hydrodynamic model is based on 3D Navier-Stokes equations including the standard k-ε model for turbulence processes,and the sediment transport model is based on a 3D equation for the re-suspended sediment transport.The computation results show that the areas of sediment concentration and transport(whether by resuspension or by the bedload) depend mainly on the position,the speed of the ship in the waterways,the kinematics of ship-generated waves and on the return flows.Thus,a map of sediment distribution and the modes of sediment transport generated by the passage of the ship are presented. 相似文献
15.
This paper presents a 3D model in sigma coordinates. Although the principles it is based on have been established for some time, some original aspects for this type of 3D mode splitting model are presented here. The model was designed to simulate flows in various coastal areas from the regional scale down to the inshore scale of small bays or estuaries where circulation is generally driven by a mix of processes. The processes to be modeled enable simplifications of the Navier–Stokes equations on the classic Boussinesq and hydrostatic hypotheses. These equations are transformed within a sigma framework to make free surface processing easier. The main point of our demonstration focuses on the original aspect of the coupling between barotropic and baroclinic modes especially designed for ADI. It explains how full consistency of the transport calculated within the 2D and 3D equation sets was obtained. Lastly, we describe the physical processes simulated on a realistic configuration at a regional scale in the Bay of Biscay. 相似文献
16.
17.
18.
This study focuses on the development of a next generation multiobjective evolutionary algorithm (MOEA) that can learn and exploit complex interdependencies and/or correlations between decision variables in monitoring design applications to provide more robust performance for large problems (defined in terms of both the number of objectives and decision variables). The proposed MOEA is termed the epsilon-dominance hierarchical Bayesian optimization algorithm (ε-hBOA), which is representative of a new class of probabilistic model building evolutionary algorithms. The ε-hBOA has been tested relative to a top-performing traditional MOEA, the epsilon-dominance nondominated sorted genetic algorithm II (ε-NSGAII) for solving a four-objective LTM design problem. A comprehensive performance assessment of the ε-NSGAII and various configurations of the ε-hBOA have been performed for both a 25 well LTM design test case (representing a relatively small problem with over 33 million possible designs), and a 58 point LTM design test case (with over 2.88×1017 possible designs). The results from this comparison indicate that the model building capability of the ε-hBOA greatly enhances its performance relative to the ε-NSGAII, especially for large monitoring design problems. This work also indicates that decision variable interdependencies appear to have a significant impact on the overall mathematical difficulty of the monitoring network design problem. 相似文献
19.
A 2D depth‐averaged hydrodynamic, sediment transport and bed morphology model named STREMR HySeD is presented. The depth‐averaged sediment transport equations are derived from the 3D dilute, multiphase, flow equations and are incorporated into the hydrodynamic model STREMR. The hydrodynamic model includes a two‐equation turbulence model and a correction for the mean flow due to secondary flows. The suspended sediment load can be subdivided into different size classes using the continuum (two‐fluid) approach; however, only one bed sediment size is used herein. The validation of the model is presented by comparing the suspended sediment transport module against experimental measurements and analytical solutions for the case of equilibrium sediment‐laden in a transition from a rigid bed to a porous bed where re‐suspension of sediment is prevented. On the other hand, the bed‐load sediment transport and bed evolution numerical results are compared against bed equilibrium experimental results for the case of a meander bend. A sensitivity analysis based on the correction for secondary flow on the mean flow including the effect of secondary flow on bed shear stresses direction as well as the downward acceleration effect due to gravity on transverse bed slopes is performed and discussed. In general, acceptable agreement is found when comparing the numerical results obtained with STREMR HySeD against experimental measurements and analytical solutions. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
20.
The boundary layer formed under the footprint of an internal solitary wave is studied by numerical simulation for waves of
depression in a two-layer model of the density stratification. The inviscid outer flow, in the perspective of boundary-layer
theory, is based on an exact solution for the long wave-phase speed, yielding a family of fully nonlinear solitary wave solutions
of the extended Korteweg–de Vries equation. The wave-induced boundary layer corresponding to this outer flow is then studied
by means of simulation employing the Reynolds-averaged Navier–Stokes (RANS) formulation coupled with a turbulence closure
model validated for wall-bounded flows. Boundary-layer characteristics are computed for an extensive range of environmental
conditions and wave amplitudes. Boundary-layer transition, identified by monitoring the eddy viscosity, is correlated in terms
of a boundary-layer Reynolds number. The frictional drag is evaluated for laminar, transitional, and turbulent cases, and
correlations are presented for the friction coefficient plus relevant measures of the boundary-layer thickness. 相似文献