首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 285 毫秒
1.
Large Eddy Simulation for Plunge Breaker and Sediment Suspension   总被引:1,自引:1,他引:1  
BAI  Yuchuan 《中国海洋工程》2002,16(2):151-164
Breaking waves are a powerful agent for generating turbulence that plays an important role in many fluid dynamical processes, particularly in the mixing of materials. Breaking waves can dislodge sediment and throw it into suspension, which will then be carried by wave-induced steady current and tidal flow. In order to investigate sediment suspension by breaking waves, a numerical model based on large-eddy-simulation (LES) is developed. This numerical model can be used to simulate wave breaking and sediment suspension. The model consists of a free-surface model using the surface marker method combined with a two-dimensional model that solves the flow equations. The turbulence and the turbulent diffusion are described by a large-eddy-simulation (LES) method where the large turbulence features are simulated by solving the flow equations, and a subgrid model represents the small-scale turbulence that is not resolved by the flow model. A dynamic eddy viscosity subgrid scale stress model has been used for the  相似文献   

2.
A boundary layer flow under spilling breakers in a laboratory surf zone with a smooth bottom is investigated using a high resolution particle image velocimetry (PIV) technique. By cross-correlating the images, oscillatory velocity profiles within a viscous boundary layer of O(1) mm in thickness are resolved over ten points. Using PIV measurements taken for an earlier study and the present study, flow properties in the wave bottom boundary layer (WBBL) over the laboratory surf zone are obtained, including the mean velocities, turbulence intensity, Reynolds stresses, and intermittency of coherent events. The data are then used to estimate the boundary layer thickness, phase variation, and bottom shear stress. It is found that while the time averaged mass transport inside the WBBL is onshore in the outer surf zone, it changes to offshore in the inner surf zone. The zero Eulerian mass transport occurs at h/hb ≈ 0.92 in the outer surf zone. The maximum overshoot of the streamwise velocity and boundary layer thickness are not constant across the surf zone. The bottom shear stress is mainly contributed by the viscous stress through mean velocity gradient while the Reynolds stress is small and negligible. The turbulence level is higher in the inner surf zone than that in the outer surf zone, although only a slight increase of turbulent intensity is observed inside the WBBL from the outer surf zone to the inner surf zone. The variation of phase inside and outside the WBBL was examined through the spatial velocity distribution. It is found the phase lead is not constant and its value is significantly smaller than previous thought. By analyzing instantaneous velocity and vorticity fields, a remarkable number of intermittent turbulent eddies are observed to penetrate into the WBBL in the inner surf zone. The size of the observed large eddies is about 0.11 to 0.16 times the local water depth. Its energy spectra follow the − 5/3 slope in the inertial subrange and decay exponentially in the dissipation subrange.  相似文献   

3.
A laboratory study on the turbulence and wave energy dissipations of spilling breakers in a surf zone is presented. Instantaneous velocity fields of propagating breaking waves on a 1/20 slope were measured using Particle Image Velocimetry (PIV). Due to the large region of the evolving wave breaking generated turbulent flow, seven PIV fields of view (FOVs) were mosaicked to form a continuous flow field in the surf zone. Mean and turbulence quantities were extracted by ensemble averaging 25 repeated instantaneous measurements at each FOV. New results for distribution and evolution of turbulent kinetic energy, mean flow energy, and total energy across the surf zone were obtained from analyzing the data. The turbulence dissipation rate was estimated based on several different approaches. It was found that the vertical distribution of the turbulence dissipation rate decays exponentially from the crest level to the bottom. The resulting energy budget and energy flux were also calculated. The calculated total energy dissipation rate was compared to that based on a bore approximation. It was found that the ratio of turbulence dissipation rate to total energy dissipation rate was about 0.01 in the outer surf zone and increased to about 0.1 after the breaking waves transformed into developed turbulent bores in the inner surf zone.  相似文献   

4.
The flocculation of cohesive sediment in the presence of waves is investigated using high-resolution field observations and a newly-developed flocculation model based on artificial neural networks. Vertical profiles of suspended sediment concentration and turbulent intensity are estimated using measurements of current profile and acoustic backscatter. The vertical distribution of floc size is estimated using an artificial neural network (ANN) that is trained and validated using floc size measurements at one vertical level. Data analysis suggests a linear correlation between suspended sediment concentration and turbulence intensity. Observations and numerical simulations show that floc size is inversely related to sediment concentration, turbulence intensity and water temperature. The numerical results indicate that floc growth is supported by low concentration and low turbulence. In the vertical direction, mean size of flocs decreases toward the bottom, suggesting floc breakage due to increasing turbulence intensity toward the bed. A significant decrease in turbulent shear could occur within the bottom few-cm, related to increased damping of turbulence by sediment induced density stratification. The results of the numerical simulations presented here are consistent with the concept of a cohesive sediment particle undergoing aggregation-fragmentation processes, and suggest that the ANN can be a precise tool to study flocculation processes.  相似文献   

5.
Active turbulence in lakes is confined to the surface mixed layer, to boundary layers on the lake sides and bottom, and to turbulent patches in the interior. The density stratification present in most lakes fundamentally alters the pathways connecting external mechanical energy inputs, for example by wind, with its ultimate fate as dissipation to heat; the density stratification supports internal waves and intrusions that distribute the input energy throughout the lake. Intrusions may be viewed as internal waves with zero horizontal wavenumber and are formed each time localised mixing occurs in a stratified fluid. Intrusions are also formed in the epilimnion by differential heating or cooling and by differential deepening. The fraction of lake volume below the diurnal mixed layer that is subject to active turbulence is very small, probably of the order of 1% or less in small to medium‐sized lakes. By contrast, in the surface mixed layer, turbulence is less intermittent and maintains phytoplankton in suspension and controls their exposure to the underwater solar flux. Nutrient transport to individual cells depends not only on the cell Reynolds number but also on the Peclet number, which, if large, implies enhanced mass transfer above purely diffusive values.  相似文献   

6.
《Coastal Engineering》1999,36(1):59-85
Simple theoretical models to determine the equilibrium profile shape under breaking and non-breaking waves are presented. For the case of breaking waves, it is assumed that the seaward transport in the undertow is locally balanced by a net vertical sedimentation, so that no bottom changes occur at equilibrium. The parameterization of the water and sediment flux in the surf zone yields a power curve for the equilibrium profile with a power of 2/3, which is in agreement with previous field investigations on surf zone profile shapes. Three different models were developed to derive the profile shape under non-breaking waves, namely (1) a variational formulation where the wave energy dissipation in the bottom boundary layer is minimized over the part of the profile affected by non-breaking waves, (2) an integration of a small-scale sediment transport formula over a wave period where the slope conditions that yield zero net transport determine equilibrium, and (3) a conceptual formulation of mechanisms for onshore and offshore sediment transport where a balance between the mechanisms defines equilibrium conditions. All three models produced equilibrium profile shapes of power-type with the power typically in the range 0.15–0.30. Comparison with field data supported the results obtained indicating different powers for the equilibrium profile shape under breaking and non-breaking waves.  相似文献   

7.
通过大尺度水槽波浪引起泥沙悬移的动床模型实验,研究了沙坝海岸破波带内水底悬沙浓度形成机理,通过比较时间平均水底悬沙浓度与时间平均水底波浪水质点动能或时间平均水底湍动能之间的相关性,论证了利用时间平均湍动能比利用时间平均波浪水质点动能计算时间平均水底悬沙浓度更为适用,并提出了以上时间平均水底悬沙浓度与水底湍动能之间的关系也可以用来近似表达时间变化的水底悬沙浓度与时间变化的水底湍动能之间的关系。研究针对规则波、波群和不规则波3种波浪形态进行,并分别对破波带内的爬坡区、内破波区和沙坝区3个区域实验结果进行讨论。  相似文献   

8.
Based on a wave bottom boundary layer model and a sediment advection-diffusion model, seven turbulence schemes are compared regarding their performances in prediction of near-bed sediment suspension beneath waves above a plane bed. These turbulence algorithms include six empirical eddy viscosity schemes and one standard two-equation k-ε model. In particular, different combinations of typical empirical formulas for the eddy viscosity profile and for the wave friction factor are examined. Numerical results are compared with four laboratory data sets, consisting of one wave boundary layer hydrodynamics experiment and three sediment suspension experiments under linear waves and the Stokes second-order waves. It is shown that predictions of near-bed sediment suspension are very sensitive to the choices of the empirical formulas in turbulence schemes. Simple empirical turbulence schemes are possible to perform equally well as the two-equation k-ε model. Among the empirical schemes, the turbulence scheme, combining the exponential formula for eddy viscosity and Swart formula for wave friction factor, is the most accurate. It maintains the simplicity and yields identically good predictions as the k-ε model does in terms of the wave-averaged sediment concentration.  相似文献   

9.
渤海海峡沉积物输运的参数化计算   总被引:1,自引:1,他引:0  
本文以2018年冬季渤海海峡两个站位的定点连续观测数据为基础,使用一维参数化方案,计算了观测站位底边界层内的水平悬浮物输运通量以及推移质输运量。在参数化方案中,简化的一维对流扩散方程被用于计算底边界层内的垂向悬浮物浓度。为了验证参数化方案的可靠性,本文基于观测数据对比了两种底剪切应力计算模型、四种临界起动剪切应力计算方法和两种一维对流扩散方程解法。对比结果表明:(1)不同模型计算的底剪切应力结果相近;(2)临界起动剪切应力受到颗粒间黏性作用的影响;(3)一维对流扩散方程的求解过程需要考虑沉积物浓度的分层效应和不同粒级颗粒临界起动剪切应力的差异。基于上述对比结果确定的最优参数化方案,进一步计算了观测站位的沉积物输运量:(1)在有再悬浮的时段,距底5 m内的水平悬浮物通量占全水深悬浮物通量的比例(T01站约为21%,T02站约为17%)显著高于相同层位水通量的占比;(2)依据参数化方案估算的冬季平均的悬浮物通量比忽略底边界层悬浮物浓度垂向变化的传统方法结果高约16%;(3)推移质输运量比悬移质输运量约低两个数量级。  相似文献   

10.
11.
基于2018年8月福建三沙湾湾内外共两个定点站位的船基和座底三脚架观测数据,研究了三沙湾底边界动力过程及悬沙输运特征。结果表明,三沙湾湾内湾外两个站位均表现出涨落潮历时相近但涨落潮流速明显不对称的现象,即湾内涨潮流速大于落潮流速,湾外则相反。湾内水体受淡水输入影响较大,表现出落潮期间显著的温盐层化,而涨潮期间水体混合良好;湾外水体受淡水影响不明显,表现为水体温度主导的层化。通过对底边界层动力过程的分析表明,湾内(距底0.75 m)、湾外(距底0.50 m)站位底边界层的平均摩阻流速分别是0.016 m/s、0.013 m/s,且两个站位拖曳系数基本相等(2.03×10-3),表明在相同流速下湾内站位的底部切应力更大,近底沉积物再悬浮和搬运相对湾外站位更为显著。因此观测期间悬沙浓度最大值出现在湾内站位,为109 mg/L,且悬沙在垂向上的分布可达上层水体;湾外站位悬沙浓度更低,并且底部悬浮泥沙仅能影响至距底5 m的水体。悬沙通量机制分解结果表明,三沙湾夏季的潮周期单宽悬沙从湾外向湾内方向净输运,湾内站位向湾内方向净输运74.88 g/(m·s),平流输沙占主导作用,贡献率41.7%;湾外站位向湾内方向净输运10.57 g/(m·s),主要受平流输沙和垂向净环流的控制,贡献率94.9%  相似文献   

12.
The characteristics of turbulence created by a plunging breaker on a 1 on 35 plane slope have been studied experimentally in a two-dimensional wave tank. The experiments involved detailed measurements of fluid velocities below trough level and water surface elevations in the surf zone using a fibre-optic laser-Doppler anemometer and a capacitance wave gage. The dynamical role of turbulence is examined making use of the transport equation for turbulent kinetic energy (the k-equation). The results show that turbulence under a plunging breaker is dominated by large-scale motions and has certain unique features that are associated with its wave condition. It was found that the nature of turbulence transport in the inner surf zone depends on a particular wave condition and it is not similar for different types of breakers. Turbulent kinetic energy is transported landward under a plunging breaker and dissipated within one wave cycle. This is different from spilling breakers where turbulent kinetic energy is transported seaward and the dissipation rate is much slower. The analysis of the k-equation shows that advective and diffusive transport of turbulence play a major role in the distribution of turbulence under a plunging breaker, while production and dissipation are not in local equilibrium but are of the same order of magnitude. Based on certain approximate analytical approaches and experimental measurements it is shown that turbulence production and viscous dissipation below trough level amount to only a small portion of the wave energy loss caused by wave breaking. It is suggested that the onshore sediment transport produced by swell waves may be tied in a direct way to the unique characteristics of turbulent flows in these waves.  相似文献   

13.
北戴河海滩泥沙捕获实验及其初步结果分析   总被引:2,自引:1,他引:2  
介绍了利用泥沙捕获器观测破波带泥沙垂直分布结构的现场实验方法和基本程序,以及利用实验结果计算泥沙通量的方法。研究表明,近岸带泥沙运移通量及其在垂向上的分布受破波带相对位置和海滩地形变化的影响。在破波点附近,波浪的搅动和流场作用强,泥沙运移通量增大,泥沙在波浪的作用下可以大量进入垂直水体以悬移和跃移的方式运移。在本实验中,破波点附近的泥沙在距海底100cm的垂直水体中运移,通量垂向向上逐渐减小。远离破波点,泥沙运移通量和进入垂直水体的高度明显下降。在地形变化复杂的有坝海滩,沙坝顶部的泥沙运移通量最大,泥沙进入垂直水体运移的机率增加,而在沙坝问的沟槽内,波浪和海流作用减弱,泥沙通量和垂向进入水体运移的比例下降。  相似文献   

14.
A three-dimensional wave radiation stress is introduced into the hydrodynamic sediment coupled model COHERENS-SED,which has been developed through introducing wave-enhanced bottom shear stress,wave dependent surface drag coefficient,wave-induced surface mixing,SWAN,damping function of sediment on turbulence,sediment model and depth-dependent wave radiation stress to COHERENS.The COHERENS-SED is adopted to study the effects induced by wave-induced three-dimensional longshore current on suspended sediment spreading of the Huanghe River (Yellow River) mouth.Several different cases divided by setting different wave parameters of inputting boundary waves are carried out.The modeling results agree with measurement data.In terms of simulation results,it is easy to know that three-dimensional wave radiation stress plays an obvious role when inputting boundary wave height is stronger than 3 m.Moreover,wave direction also affects the sediment spreading rules of the mouth strongly too.  相似文献   

15.
Regional Ocean Modeling System (ROMS v 3.0), a three-dimensional numerical ocean model, was previously enhanced for shallow water applications by including wave-induced radiation stress forcing provided through coupling to wave propagation models (SWAN, REF/DIF). This enhancement made it suitable for surf zone applications as demonstrated using examples of obliquely incident waves on a planar beach and rip current formation in longshore bar trough morphology (Haas and Warner, 2009). In this contribution, we present an update to the coupled model which implements a wave roller model and also a modified method of the radiation stress term based on Mellor (2008, 2011a,b,in press) that includes a vertical distribution which better simulates non-conservative (i.e., wave breaking) processes and appears to be more appropriate for sigma coordinates in very shallow waters where wave breaking conditions dominate. The improvements of the modified model are shown through simulations of several cases that include: (a) obliquely incident spectral waves on a planar beach; (b) obliquely incident spectral waves on a natural barred beach (DUCK'94 experiment); (c) alongshore variable offshore wave forcing on a planar beach; (d) alongshore varying bathymetry with constant offshore wave forcing; and (e) nearshore barred morphology with rip-channels. Quantitative and qualitative comparisons to previous analytical, numerical, laboratory studies and field measurements show that the modified model replicates surf zone recirculation patterns (onshore drift at the surface and undertow at the bottom) more accurately than previous formulations based on radiation stress (Haas and Warner, 2009). The results of the model and test cases are further explored for identifying the forces operating in rip current development and the potential implication for sediment transport and rip channel development. Also, model analysis showed that rip current strength is higher when waves approach at angles of 5° to 10° in comparison to normally incident waves.  相似文献   

16.
《Coastal Engineering》2004,51(10):1021-1049
A numerical process-based model to forecast beach profile morphodynamics has been developed. In the present paper, an analysis of various modelling approaches and key parametrizations involved in the estimation of the wave-driven current and the suspended sediment concentration is carried out.Several resolution techniques for the 1DV horizontal (i.e., in the x-direction perpendicular to coastline) momentum equation governing the Mean Horizontal Velocity (MHV) are analysed. In the first kind of techniques, the mean horizontal velocity is computed from the momentum equation, whereas the Mean Water Level (MWL) is computed using a parametrization of the depth-averaged momentum equation. Two boundary or integral conditions are thus needed. In the second kind, both mean horizontal velocity and mean water level gradient in the x-direction are the unknowns of the momentum equation, thus, three boundary or integral conditions are needed. Various additional conditions are discussed. We show that using a technique of the first kind is equivalent to imposing the difference between the surface and the bottom shear stresses in the 1D vertical equation. Both techniques lead to results that are in good agreement with the Delta Flume experimental data, provided the Stokes drift flow discharge is imposed as an additional condition. The influence of the breaking roller model and of the turbulent viscosity parametrization are also analysed.Suspended sediment transport by the mean current and wave-induced bedload transport are taken into account in the sediment flux. Three turbulent diffusivity parametrizations are compared for suspended sediment concentration estimations. A linear profile for the turbulent diffusivity taking into account the wave bottom shear stress and the surface wave breaking turbulence production is shown to give the best results. Using experimental data, we put forward the poor estimation of the bottom sediment concentration given by the three implemented parametrizations. We thus propose a new parametrization relying on a Shields parameter based on the breaking roller induced surface shear stress. Using this new parametrization, the bottom profile used in the tests keeps its two bars which disappear otherwise. However, the morphodynamical model still overestimates the bars offshore motion, a bias already observed in other models.  相似文献   

17.
底边界层水沙观测系统和应用   总被引:2,自引:0,他引:2  
利用现代先进的光学和声学测量仪器(ADP、PC-ADP、ADV、OBS等),组制四角架观测系统用于底边界层的水沙观测.此观测系统既可获取垂向流速场,更可获取常规水文观测难以获取的高分辨率的近底流速、近底单点紊动过程、近底温盐沙和波浪过程.利用此观测系统于2007年8月在长江口北槽最大浑浊带内进行了实地观测.分析四角架观测数据与同步船测底部数据表明,四角架观测系统在近底(特别是底部明显分层)观测时表现出明显优势.常规水文观测测得底部含沙量一般只相当于距底70~120 cm高度含沙量;船测底部最大含沙量为5.8 kg/m 3,而架测距底30 cm处最大含沙量有25 kg/m 3;近底含沙量较高时往往出现强分层,强分层时实测距底30 cm处含沙量可达距底120 cm处含沙量的5倍;船测底部流速比架测距底30 cm处流速平均偏大45%;但距底30 cm处最大流速仍可达1.2 m/s.船测数据也验证了四角架观测数据的正确性.  相似文献   

18.
According to the mechanism of sediment suspension under waves,namely,the main reason of sediment suspension changes from the turbulent mixing in the bottom boundary layer to the periodic motion of the water particle near the free water surface,a three-hyer model d sediment concentration distribution due to waves is presented along the whole water depth based on the concept of the finite mixing |ength.The determination of the parameters in the model is discussed and an empirical formula is suggested.Comparisons between the calculated results and the measurements indicate that the resuits of the model agree well with the data from both the large and small scale flume experiments.  相似文献   

19.
Hydrodynamics and sediment transport in the nearshore zone were modeled numerically taking into account turbulent unsteady flow. The flow field was computed using the Reynolds Averaged Navier–Stokes equations with a kε turbulence closure model, while the free surface was tracked using the Volume-Of-Fluid technique. This hydrodynamical model was supplemented with a cross-shore sediment transport formula to calculate profile changes and sediment transport in the surf and swash zones. Based on the numerical solutions, flow characteristics and the effects of breaking waves on sediment transport were studied. The main characteristic of breaking waves, i.e. the instantaneous sediment transport rate, was investigated numerically, as was the spatial distribution of time-averaged sediment transport rates for different grain sizes. The analysis included an evaluation of different values of the wave friction factor and an empirical constant characterizing the uprush and backwash. It was found that the uprush induces a larger instantaneous transport rate than the backwash, indicating that the uprush is more important for sediment transport than the backwash. The results of the present model are in reasonable agreement with other numerical and physical models of nearshore hydrodynamics. The model was found to predict well cross-shore sediment transport and thus it provides a tool for predicting beach morphology change.  相似文献   

20.
Novel laboratory experiments and numerical modelling have been performed to study the advection scales of suspended sediment in the swash zone. An experiment was designed specifically to measure only the sediment picked up seaward of the swash zone and during bore collapse. The advection scales and settling of this sediment were measured during the uprush along a rigid sediment-free beach face by a sediment trap located at varying cross-shore positions. Measurements were made using a number of repeated solitary broken waves or bores. Approximately 25% of the pre-suspended sediment picked up by the bores reaches the mid-swash zone (50% of the horizontal run-up distance), indicating the importance of the sediment advection in the lower swash zone. The pre-suspended sediment is sourced from a region seaward of the shoreline (still water line) which has a width of about 20% of the run-up distance. An Eulerian–Lagrangian numerical model is used to model the advection scales of the suspended sediment. The model resolves the hydrodynamics by solving the non-linear shallow water equations in an Eulerian framework and then solves the advection–diffusion equation for turbulence and suspended sediment in a Lagrangian framework. The model provides good estimates of the measured mass and distribution of sediment advected up the beach face. The results suggest that the correct modelling of turbulence generation prior to and during bore collapse and the advection of the turbulent kinetic energy into the lower swash is important in resolving the contribution of pre-suspended sediment to the net sediment transport in the swash zone.  相似文献   

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

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