Issues in Eulerian–Lagrangian modeling of sediment transport under saltation regime     

Reza Barati  Seyed Ali Akbar Salehi Neyshabouri  Goodarz Ahmadi 《国际泥沙研究》2018,33(4):441-461

The saltation regime is very important for understanding the sediment transport mechanism. However,there is no consensus on a model for the saltation regime. This study answers several questions raised with respect to the Eulerian-Lagrangian modeling of sediment transport. The first question is why the previous saltation models that use different combinations of hydrodynamic forces yielded acceptable results? The second question is which shear lift model(i.e. a shear lift expression and its coefficient) is more appropriate? Another important question is which hydrodynamic forces have greater contributions to the saltation characteristics of a sediment particle? The last question is what are the contributions of the turbulence fluctuations as well as effects of using two-and three-dimensional(2 D and 3 D) models on the simulation results? In order to fairly answer these questions, a systematic study was done by considering different scenarios. The current study is the first attempt to clearly discuss these issues. A comprehensive 3 D saltation model for non-cohesive sediment was developed that includes all the hydrodynamic forces acting on the particle. The random nature of sediment transport was included using turbulent flow and bed-particle collision models. The eddy interaction model was applied to generate a3 D turbulent flow field. Bed-particle collisions were considered using the concept of a contact zone and a corresponding contact point. The validation of the model was done using the available experimental data for a wide range of sediment size(0.03 to 4.8 cm). For the first question, the results indicated that some of the hydrodynamic effects show opposing trends and some have negligible effects. With these opposing effects it is possible to adjust the coefficients of different models to achieve acceptable agreement with the same experimental data while omitting some aspects of the physics of the process. A suitable model for the shear lift force was developed by linking the lift coefficient to the drag coefficient and the contributions of the hydrodynamic forces and turbulence fluctuations as well as the consequences of using of 2 D and 3 D models were studied. The results indicate that the shear lift force and turbulent flow fluctuations are important factors for the saltation of both sand and gravel, and they cannot be ignored.  相似文献   

An Actively Controlled Wind Tunnel And Its Application To The Reproduction Of The Atmospheric Boundary Layer     

Cao Shuyang  Akira Nishi  Kimitaka Hirano  Shigehira Ozono  Hiromori Miyagi  Hiromori Kikugawa  Yuji Matsuda  Yasuo Wakasugi 《Boundary-Layer Meteorology》2001,101(1):61-76

An actively controlled wind tunnel equipped with multiple fansand airfoils has been developed, mainly for the purpose of reproducing the atmospheric boundary layer (ABL) for wind engineering applications. Various fluctuating flows can be achieved in this wind tunnel by altering the input data of the fans and airfoils through computer control. In this study, the ABL is physically simulated in this wind tunnel, and particular attention ispaid to the simulation of the profile of Reynolds stress. The method of generating the fluctuating flow and the experimental results of reproducing the ABL are presented. As the results show, the spatial distribution of Reynolds stress is satisfactorily simulated, and the profiles of other statisticalturbulent parameters, such as mean velocity, turbulent intensity, integral scale and power spectrum are successfully reproduced simultaneously.  相似文献   

Turbulence Spectra And Dissipation Rates Above And Within A Forest Canopy   总被引：4，自引：0，他引：4  

Shuhua Liu  Heping Liu  M. Xu  M. Y. Leclerc  Tingyao Zhu  Changjie Jin  Zhongxiang Hong  Jun Li  Huizhi Liu 《Boundary-Layer Meteorology》2001,98(1):83-102

Three velocity componentsand temperature were measured usingthree-dimensional sonic anemometers/thermometers attwo levels, above and within a forest canopy, in theChangbai Mountains of northeast China. Turbulencespectral structure, local isotropy anddissipation rates above and within the forest canopywere calculated using the eddy correlation method.Results show that the normalized turbulent spectralcurves have -2/3 slopes in the inertial subrange.While the shapes of the spectra are in good agreementwith the Kansas flat terrain results, the atmosphericturbulence is anisotropic above the forest canopy. Dueto breaking down of large eddies by the foliage,branches and trunks, the spectral peak frequencies forvelocity and temperature are higher withinthan above the forest canopy. Compared withmeasurements from previous studies over flat terrain,the velocity and temperature spectra above andinside the forest canopy appear to shift toward higherfrequencies. The turbulence is approximately isotropicin the inertial subrange within the forest canopy, and isanisotropic above the forest canopy. The turbulentkinetic energy and heat energy dissipation rates aboveand inside the forest canopy are much larger thanthose obtained by Kaimal and Hogstrom over grasslandand grazing land. The distinct features in the resultsof the present experiment may be attributed to thedynamic forcing caused by the rough surface of the forestcanopy.  相似文献   

WRF模式中QNSE方案的湍流长度尺度系数的调整试验研究   总被引：2，自引：0，他引：2  

徐慧燕  徐亚钦  王智  朱佩君  李小凡  翟国庆 《大气科学》2017,41(2):357-371

边界层参数化方案中湍流混合对数值模拟起着重要的作用,湍流混合作用的恰当描述对于温度、湿度、风场以及降水的准确模拟至关重要。我国长江中下游流域人口密集,暴雨灾害频发,很有必要寻找一种适合该地区降水模拟的边界层参数化方案。本文运用WRF（Weather Research and Forecasting）中尺度数值模式,以QNSE（Quasi-Normal Scale Elimination）边界层参数化方案为基础,将其中湍流混合长度尺度系数调整为可变参数。本文将Noh et al.（2003）提出的Prandtl公式与Janji?提出的修正湍流长度尺度系数的方法相结合,通过考虑非局地项的强迫、地表稳定度与边界层高度对湍流长度尺度系数的影响,强调大气的动力结构特征与热力结构特征对湍流混合的共同影响,从而提高QNSE边界层参数化方案在不同地理环境下的模拟能力。文章通过进行长江中下游地区的典型暴雨试验,将调整参数后的QNSE方案与原方案进行比较,重点分析调整参数后的方案与原方案对关键基本气象要素场、边界层结构特征以及降水的模拟能力,并将模拟结果与观测结果进行对比,结果表明:调整参数后的方案一定程度上改进了地表温度、边界层结构以及降水的模拟效果。进一步研究表明,调整参数后的方案主要通过改变边界层混合缓解水汽混合比、位温模拟方面的误差。  相似文献   

Bottom friction and its effects on periodic long wave propagation     

A. Orfila  G. Simarro  P.L.F. Liu 《Coastal Engineering》2007

A new set of Boussinesq-type equations describing the free surface evolution and the corresponding depth-integrated horizontal velocity is derived with the bottom boundary layer effects included. Inside the boundary layer the eddy viscosity gradient model is employed to characterize Reynolds stresses and the eddy viscosity is further approximated as a linear function of the distance measured from the seafloor. Boundary-layer velocities are coupled with the irrotational velocity in the core region through boundary conditions. The leading order boundary layer effects on wave propagation appear in the depth-integrated continuity equation to account for the velocity deficit inside the boundary layer. This formulation is different from the conventional approach in which a bottom stress term is inserted in the momentum equation. An iterative scheme is developed to solve the new model equations for the free surface elevation, depth-integrated velocity, the bottom stress, the boundary layer thickness and the magnitude of the turbulent eddy viscosity. A numerical example for the evolution of periodic waves propagating in one-dimensional channel is discussed to illustrate the numerical procedure and physics involved. The differences between the conventional approach and the present formulation are discussed in terms of the bottom frictional stress and the free surface profiles.  相似文献   

沿海市区近地面大气层风、温、湿湍流特征的观测和分析          下载免费PDF全文

赵永平  陈永利  杨连素  范继铨 《海洋科学》1993,17(2):52-56

用三向风速仪和高灵敏温湿度脉动观测装置,测量沿海城市上空大气表面层风速、温度和湿度脉动。并用涡动相关法计算了大气表面层动量、热量和水汽湍流通量,给出了复杂地形上各物理量的谱特征,讨论了各湍流通量的特性。结果指出,风速、温度和湿度谱密度分布基本上都服从—5/3律。由于建筑物的影响,在f≥0.2Hz频段内拢动能量明显偏大。动量和热量输送主要由f<0.2Hz的扰动完成的,高频扰动的贡献可以忽略不计。  相似文献   

Microstructure measurements in natural waters: Methodology and applications   总被引：2，自引：0，他引：2  

Elena Roget  Iossif Lozovatsky  Xavier Sanchez  Manuel Figueroa   《Progress in Oceanography》2006,70(2-4):126

Modern approaches to microstructure data processing, including wavelet denoising, are discussed. The wavelet procedure is applied to small-scale shear signals before estimating the dissipation rate ε and to the temperature/density profiles used to calculate Thorpe scales. Microstructure data obtained on the Mediterranean shelf of Catalonia are used to illustrate various approaches to the Thorpe displacement calculations. It is suggested that the Weibull probability function is an appropriate model for the Thorpe scale distribution. Microstructure measurements from the upper layer of the Boadella reservoir (Catalonia, Spain) support this finding.A new analytical approximation for the 1D Panchev–Kesich spectrum is deduced and the results of ε computation are compared with spectral fitting by the widely used Nasmyth spectrum. Applying the Kraichnan spectral model to compute ε from temperature spectra in the convective-viscous sub-range is examined as an alternative to the Batchelor spectrum. Microstructure measurements taken in Lake Banyoles (Catalonia, Spain) and in the North Atlantic were used for spectral calculations.Statistical analysis of eddy K_b and thermal K_θ diffusivities measured on a shallow shelf of the Black Sea shows the importance of process-orientated domain averaging of the diffusivities in obtaining good correspondence between K_b and K_θ in active turbulent regions. In weakly turbulent, stratified interior layers, the averaged K_b and K_θ differ significantly, which may point to the inapplicability of isotropic formulae used for ε and temperature dissipation χ_θ estimates, as well as to a dependence of the mixing efficiency γ on the Richardson number or in some cases on regions of fossil turbulence.  相似文献   

Calculation of diffusivities of passive fields in turbulent media     

N. A. Silant'ev 《地球物理与天体物理流体动力学》2013,107(2-4):183-197

Abstract

The exact numerical and approximate analytical solutions of the simplest nonlinear integral equation with second order nonlinearity for the averaged Green function are presented. It is assumed that the turbulence is stationary, homogeneous, isotropic and incompressible. Numerous examples of turbulent spectra are considered (peak-like spectrum, spectra of Kolmogorov's type with different forms of “pumping” regions, stepwise spectra etc.). Special emphasis is given to investigating the case of so called “frozen” turbulence when the parameter ξ =u ₀τ/R→∞ where uτ₀,R ₀ are characteristic velocity, lifetime and space scale of turbulent pulsations, respectively. It is shown that these solutions allow us to calculate the turbulent diffusivities accurately for arbitrary spectra with any values of the parameter ξ. The results take into account the possible helicity of turbulence concerned only with scalar passive fields (number density and temperature).  相似文献   

A depth-integrated model for weakly dispersive, turbulent, and rotational fluid flows   总被引：2，自引：0，他引：2  

Dae-Hong Kim  Patrick J. Lynett  Scott A. Socolofsky 《Ocean Modelling》2009,27(3-4):198-214

A set of weakly dispersive Boussinesq-type equations, derived to include viscosity and vorticity terms in a physically consistent manner, is presented in conservative form. The model includes the approximate effects of bottom-induced turbulence, in a depth-integrated sense, as a second-order correction. Associated with this turbulence, vertical and horizontal rotational effects are captured. While the turbulence and horizontal vorticity models are simplified, a model with known physical limitations has been derived that includes the quadratic bottom friction term commonly added in an ad hoc manner to the inviscid equations. An interesting result of this derivation is that one should take care when adding such ad hoc models; it is clear from this exercise that (1) it is not necessary to do so – the terms can be included through a consistent derivation from the viscous primitive equations – and (2) one cannot properly add the quadratic bottom friction term without also adding a number of additional terms in the integrated governing equations. To solve these equations numerically, a highly accurate and stable model is developed. The numerical method uses a fourth-order MUSCL-TVD scheme to solve the leading order (shallow water) terms. For the dispersive terms, a cell averaged finite volume method is implemented. To verify the derived equations and the numerical model, four cases of verifications are given. First, solitary wave propagation is examined as a basic, yet fundamental, test of the models ability to predict dispersive and nonlinear wave propagation with minimal numerical error. Vertical velocity distributions of spatially uniform flows are compared with existing theory to investigate the effects of the newly included horizontal vorticity terms. Other test cases include comparisons with experiments that generate strong vorticity by the change of bottom bathymetry as well as by tidal jets through inlet structures. Very reasonable agreements are observed for the four cases, and the results provide some information as to the importance of dispersion and horizontal vorticity.  相似文献   

Characteristics of turbulent boundary layers over a rough bed under saw-tooth waves and its application to sediment transport   总被引：1，自引：0，他引：1  

Suntoyo  Hitoshi Tanaka  Ahmad Sana   《Coastal Engineering》2008,55(12):1102-1112

A large number of studies have been done dealing with sinusoidal wave boundary layers in the past. However, ocean waves often have a strong asymmetric shape especially in shallow water, and net of sediment movement occurs. It is envisaged that bottom shear stress and sediment transport behaviors influenced by the effect of asymmetry are different from those in sinusoidal waves. Characteristics of the turbulent boundary layer under breaking waves (saw-tooth) are investigated and described through both laboratory and numerical experiments. A new calculation method for bottom shear stress based on velocity and acceleration terms, theoretical phase difference, φ and the acceleration coefficient, a_c expressing the wave skew-ness effect for saw-tooth waves is proposed. The acceleration coefficient was determined empirically from both experimental and baseline k–ω model results. The new calculation has shown better agreement with the experimental data along a wave cycle for all saw-tooth wave cases compared by other existing methods. It was further applied into sediment transport rate calculation induced by skew waves. Sediment transport rate was formulated by using the existing sheet flow sediment transport rate data under skew waves by Watanabe and Sato [Watanabe, A. and Sato, S., 2004. A sheet-flow transport rate formula for asymmetric, forward-leaning waves and currents. Proc. of 29th ICCE, ASCE, pp. 1703–1714.]. Moreover, the characteristics of the net sediment transport were also examined and a good agreement between the proposed method and experimental data has been found.  相似文献