不同湍流模式下钱塘江涌潮水流三维模拟

钱塘江涌潮具有动力作用强和流速变化快等特点。涌潮水流紊动复杂,流速的垂向分布和紊动强度息息相关。通过涌潮水流实测资料的分析可以发现,涌潮作用下流速垂向分布在底部和上层存在差异。为研究涌潮作用下流速垂向分布的特征,应用基于非结构网格下有限体积法模型FVCOM对钱塘江涌潮河段水流运动进行三维数值模拟。考虑到涌潮紊动作用复杂且对流速的垂向分布起着重要影响,采用不同的湍流模式对涌潮传播过程中水流的运动特征开展研究。通过与涌潮河段实测资料的验证,复演涌潮到达前后水流运动特征,给出涌潮水流湍动能的变化过程。研究成果有助于深入认识涌潮水流紊动特征和流速的分布规律,为涌潮作用下物质输运的研究提供基础。     

A three-dimensional multi-level turbulence model for tidal motion

A three-dimensional multi-level turbulence model is developed to simulate tide induced circulation in coastal waters. Based on the bathymetry data, the coastal waters are divided into a number of layers. In every layer, the velocities are integrated along the layer depth. The eddy viscosity and diffusivity are computed from the Prandtl mixing length turbulence model. This multi-level model solves for the water surface elevations and currents in different water depths. Comparison of numerical results with the measured data shows good conformity.     

Combined tidal and wind driven flows and bedload transport over a flat bottom

The combined tidal and wind driven flow and resulting sediment transport in the ocean over a flat bottom at intermediate water depth has been investigated, using a simple one dimensional two-equation turbulence closure model. This model has been verified against field measurements of a tidal flow in the Celtic Sea. The tidal velocity ellipses and the time series of the horizontal velocity components at given elevations above the bottom are well predicted through the water column although there are some deviations between the predicted and measured velocities near the bottom due to the uncertainty of the bottom roughness. For the combined tidal and wind driven flows the velocity profiles, turbulent kinetic energy profiles and surface particle trajectories are predicted for weak and strong winds. Furthermore, the bottom shear stress and the resulting bedload transport have been predicted; the parts of the particle trajectories in the close vicinity of the bottom where the bedload transport exists are displayed. Finally, the direction and magnitude of the surface drift, the depth-averaged mean velocity and the mean bedload transport are given, and the effect of the bottom roughness on the sea surface drift is investigated.     

涌潮水沙动力过程现场观测研究进展

涌潮是潮波传播过程中产生的一种自然现象,是潮波非线性畸变的结果。涌潮有波状和破碎形态之分,波状涌潮是一系列平行向前传播的涌波构成的波列,破碎涌潮则是前锋陡立向前推进的水滚。基于国内外涌潮水沙动力过程现场观测的主要成果,从形成机理出发,归纳涌潮生成的必要条件,剖析潮波运动非线性和摩擦效应对涌潮生成的影响;针对典型的波状和破碎涌潮,总结潮头的自由表面特征参数、流动结构和传播演化特征;回顾涌潮局部湍流和混合过程、泥沙输运和沉积的研究进展,评述涌潮脉冲过程对河口生态环境的影响。涌潮的周期性传播引起自然系统的大规模混合,对潮汐河口区域的生态环境平衡具有重要意义。涌潮现象的研究推动潮汐学的发展,现场观测是涌潮研究的基础。随着仪器设备和分析手段的进步,涌潮多尺度生成与演化机制、涌潮多物理过程耦合作用机理和涌潮脉冲过程的生态效应定量评价是今后需要深入研究的问题。     

河口底边界层湍流观测后处理技术方法分析

河口底边界层过程是河口海岸研究与工程应用中的重要内容。三维点式高频流速仪(ADV)已经成为湍流现场观测的最有效的工具之一,然而受测量状态、复杂的波流环境、底床几何结构等因素的影响,湍流观测的后处理目前还不成熟。在前人工作的基础上,提出了河口底边界层湍流观测后处理的综合技术方法,包括测量状态判断、数据质量检测、坐标系旋转、去除毛刺及滤波,探讨了这些处理方法中的某些步骤及处理顺序对湍流参数估算可能产生的影响,提出了综合后处理技术的准确性评估方法。该研究对于近岸海洋湍流混合、泥沙输运等重要问题的解决可以提供较为扎实的技术支持。     

A Numerical Model for Mixing Areas of Pollutant Discharging into Tidal Flows

- A combined numerical model for computing mixing areas of pollution vertical jet discharging into tidal flows has been developed. This numerical model is composed of a 2-D depth-averaged dynamic far-field numerical model and a 2-D vertical dynamic near-field numerical model. The former uses finite node method to compute velocity, and improved finite node method to compute pollutant concentration distribution; the latter is a k-e turbulence model, and uses SIMPLE (Semi-Implicit Method for Pressure Linked Equations) method to compute velocity. At the same time, the former provides boundary conditions for the latter. This model can simulate both far- field pollutant concentration distribution and near-field vertical recirculation quickly and precisely. This model has been verified by measured data of pollutant F of the Dachang reach of the Changjing River and test data presented by Chen el al. (1992). On the basis of verification, the authors use a designed hydrograph to compute this mixing area for a cer     

Turbulent mixing in a small estuary: Detailed measurements

Between 2003 and 2007, a series of field studies were performed in a typical small coastal plain type estuary (Eprapah Creek) located on the Southeast coast of Australia. The aim of these field studies was to investigate the turbulence and turbulent mixing properties in the estuarine zone. During these studies, high frequency turbulence and physio-chemistry data were collected continuously over a relatively long duration (up to 50 h). This article provides a summary of the key outcomes of these studies, highlighting the implications that these findings have on the modelling of small estuaries. The studies showed that the response of the turbulence and water quality properties were distinct under spring and neap tidal forcing and behaved differently in the middle and upper estuarine zones. The behaviour of turbulence properties to spring tidal forcing differed from that observed in larger estuaries and seemed unique to small estuarine systems. An investigation of several key turbulence parameters used in the modelling of estuarine mixing showed that many assumptions used in larger estuaries must be applied with caution or are simply untrue in small estuaries. For example the assumption that the mixing coefficient parameters are constant over the tidal cycle in a small estuary is simply untrue. These distinctions between the turbulence and mixing properties in small and large estuarine systems highlight the need for the continued study of small estuaries, so this type of system can be properly understood.     

Field studies of velocity,salinity and suspended solids concentration in a shallow tidal channel near tidal flap gates

The design and operation of mathematical models of solute mixing and sediment transport in estuaries rely heavily on the provision of good-quality field data. We present some observations of salinity, suspended sediment concentration and velocity at one of the tidal limits of a semi-enclosed tidal lagoon in Southern England (Pagham Harbour, West Sussex, UK) where the natural processes of tidal incursion and solute mixing have been heavily modified as a result of the construction of sea walls dating back to the 18th Century. These observations, made immediately downstream of two parallel tidal flap gates by conductivity-temperature-depth (CTD) profiler, and also using bed-mounted sensor frames to measure velocity at 2 fixed depths, have yielded a set of results covering 11 tidal cycles over the period 2002–04. It is clear from the results obtained that over a typical tidal cycle, the greatest vertical salinity gradients occur in the 1–2 h immediately after the onset of the flood tide, and that subsequently, energetic mixing acts to rapidly break down this stratification. Under moderate-to-high fresh water flows (>0.5 m³/s), the break-down in vertical salinity gradient is more gradual, while under low fresh water flows (<0.2 m³/s), the vertical salinity gradient is generally less pronounced. Estimates of Richardson number during the early flood-tide period reveal values that vary rapidly between <1 and about 20, with lower values occurring after around 1.5–2 h after low water. Observations of suspended sediment concentration vary widely even for similar tidal and fresh water flow conditions, revealing the possible influence of wind speed, the storage effects of the water in the lagoon downstream of the observation site, and the complexity of the hydrodynamics downstream of tidal flap gates. The data also show that most of the sediment transport is landward, and occurs during flood tides, with estimated total tidal landward flood tide flux of fine sediment of the order of 50–120 kg under low fresh water flow conditions. These observations, which reinforce the results presented in Warner et al. (2004) and elsewhere, can help to provide information about the appropriate techniques for managing sediments and pollutants, including nutrients from sewage effluent waters, in estuaries where hydraulic flap gates are used to control the entry of fresh water over the tidal cycle.     

Numerical modelling of suspended sediment fluxes in estuarine waters

A two-dimensional finite difference numerical model, capable of predicting depth-averaged tidal flow fields in coastal and estuarine waters, has been extended to include tide-induced non-cohesive sediment transport processes. The partial differential equations governing the conservation of mass, momentum and suspended sediment in an incompressible turbulent flow are included in a depth-integrated form in the model. For the representation of the processes of erosion and deposition of sediment from the bed an empirically based source-sink term was refined, based on the results of three mobile bed flume studies. The model has been tested by simulating tidal flows and suspended sediment fluxes in two estuaries, with particular application to the Humber estuary in the U.K. The model was calibrated and found to produce an encouraging degree of agreement between the numerical predictions and corresponding field measurements for this estuary. Furthermore, the predicted gross deposition and erosion features of the estuary were found to be in close agreement with interpretations from Eulerian tidal residual predictions.     

Effects of time dependence in unstratified tidal boundary layers: results from large eddy simulations

A three-dimensional Large Eddy Simulation (LES) model is used to simulate oscillating tidal boundary layers and test previous results obtained from one-dimensional boundary layer models and turbulence measurements in tidal channels. The LES model produces low-order turbulence statistics in agreement with the semi-analytic theory and observations. It shows a logarithmic layer in the mean velocity profile and a linear distribution of Reynolds stress with water depth. However, the eddy viscosity profile predicted by the LES model is not parabolic but better matches a parabolic profile modified by wake effect observed in the outer part of depth-limited steady boundary layers. Low-order turbulence statistics can be scaled by the instantaneous friction velocity at the bottom boundary. Although turbulence intensities in three directions fluctuate over a tidal cycle, their normalized values are in good agreement with those determined from laboratory experiments of steady open-channel flows. The LES model confirms that tidal turbulence is in quasi-equilibrium. However, it also demonstrates the importance of flow acceleration/deceleration term in the depth-integrated momentum balance for the mean flow. Phase differences are found between flows at different heights above the bottom boundary.     

Mixing parameterizations in ocean climate modeling

Results of numerical experiments with an eddy-permitting ocean circulation model on the simulation of the climatic variability of the North Atlantic and the Arctic Ocean are analyzed. We compare the ocean simulation quality with using different subgrid mixing parameterizations. The circulation model is found to be sensitive to a mixing parametrization. The computation of viscosity and diffusivity coefficients by an original splitting algorithm of the evolution equations for turbulence characteristics is found to be as efficient as traditional Monin–Obukhov parameterizations. At the same time, however, the variability of ocean climate characteristics is simulated more adequately. The simulation of salinity fields in the entire study region improves most significantly. Turbulent processes have a large effect on the circulation in the long-term through changes in the density fields. The velocity fields in the Gulf Stream and in the entire North Atlantic Subpolar Cyclonic Gyre are reproduced more realistically. The surface level height in the Arctic Basin is simulated more faithfully, marking the Beaufort Gyre better. The use of the Prandtl number as a function of the Richardson number improves the quality of ocean modeling.     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

Cohesive sediment transport in the 3D-hydrodynamic-baroclinic circulation model,: study case for idealized tidal inlet

This work provides a general hydrodynamic circulation model that can be used to understand density driven flows, which may arise in the case of suspension of fine-grained materials. The research is expected to provide a better understanding of the characteristics of spatial and temporal variability of current, which is associated with the period of ebb and flood tidal cycles.The model development includes extending the existing three-dimensional (3D) ADCIRC model with (1) baroclinic forcing term and (2) transport module of suspended and soluble materials. The transport module covers the erosion, material suspension and deposition processes for cohesive type sediment. In the case of an idealized tidal inlet in stratified water, the inclusion of baroclinic term can demonstrate the prevailing longshore sediment transport. It is shown that the model has application to the transport of the cohesive sediments from the mouth of the Mississippi River along the north shore of the Gulf of Mexico towards and along the Texas coast.     

The effect of standing biomass on flow velocity and turbulence in Spartina alterniflora canopies

Plant-flow interactions on the surface of tidal wetlands result in flow characteristics that are profoundly different from non-vegetated flows. Reductions in mean flow velocity and turbulence, especially the vertical components, limit vertical mixing and may impact a wide range of processes including geochemical exchanges at the sediment water interface, larval recruitment and dispersion, and sediment deposition and retention. The goal of this paper is to quantify horizontal and vertical components of velocity, turbulence intensity and total turbulent kinetic energy in Spartina alterniflora canopies in southeastern North Carolina and to relate flow characteristics to particulate transport on the marsh surface. Another aim of this paper is to assess the extent to which the distribution of standing biomass affects mean flow and turbulence by comparing S. alterniflora data to other canopy types and through a series of canopy manipulations which altered canopy height and stem densities.The results of this study indicate that flow velocity, turbulence intensity, and total turbulent kinetic energy (TKE) are significantly reduced within the vegetated canopy and that this reduction is inversely related to the amount of biomass present in the water column. Within the canopy, approximately 50% of the initial mean velocity and TKE is reduced within 5 m of the canopy edge. Within the canopy, mean velocity and TKE_horiz usually exceeded vertical velocity or TKE_vert and the vertical components of flow were attenuated more strongly than the horizontal. These results suggest that within the vegetation, turbulence contributes more to lateral advection than to vertical mixing. As a result, total suspended solid concentrations were shown to decrease logarithmically with distance from the canopy edge and to decrease at a faster rate in more densely vegetated regions of the canopy (i.e. lower TKE_vert) as compared to areas of sparser vegetation (i.e. higher TKE_vert).     

潮汐河口环流、湍流、混合与层化的物理学

简要回顾了潮汐河口环流、湍流、混合与层化的基本物理概念、内涵、研究方法、研究成果,指出了主要的研究进展,最后,展望了今后的研究方向。本文不考虑悬沙和风浪的影响。经典的河口环流也因潮汐应变的出现而受到挑战,河口环流由重力环流和潮汐应变环流构成。"涡黏度-剪切协方差(ESCO)"概念的提出,又区分出重力ESCO环流与潮汐ESCO环流。横向环流,尤其具有曲率的弯道中的横向环流,也得到进一步的理解。涡度方法的应用,揭示横向环流不仅由各种不同物理机制造成,而且对纵向河口重力环流有重要的影响。分层流中剪切湍流的理论加深了人们对潮汐河口湍流、混合的物理学的认识,势能差异方程更是使得定量理解潮汐河口混合与层化的三维时间、空间变化成为可能。     

On the turbulent prandtl number in a stably stratified atmospheric boundary layer

The results obtained from both atmospheric and laboratory measurements and from LES data show that, in the stably stratified flows of the atmospheric boundary layer, turbulent mixing occurs at gradient Richardson numbers Ri _g that significantly exceed one: the inverse turbulent Prandtl number Pr _t ⁻¹ decreases with an increase in the thermal flow stability. The decreasing trend of the inverse turbulent Ptandtl number is reproduced in a stably stratified atmospheric boundary layer in agreement with measurement data with the aid of an improved three-parameter turbulence model. In this model, a modified model that takes into account the effect of stratification in the expression for the time scale of the scalar field is used for the pressure-scalar correlation.     

Tidal Mixing in the Kuril Straits and Its Impact on Ventilation in the North Pacific Ocean

A numerical study using a 3-D nonhydrostatic model has been applied to baroclinic processes generated by the K ₁ tidal flow in and around the Kuril Straits. The result shows that large-amplitude unsteady lee waves are generated and cause intense diapycnal mixing all along the Kuril Island Chain to levels of a maximum diapycnal diffusivity exceeding 10³ cm²s⁻¹. Significant water transformation by the vigorous mixing in shallow regions produces the distinct density and potential vorticity (PV) fronts along the Island Chain. The pinched-off eddies that arise and move away from the fronts have the ability to transport a large amount of mixed water (∼14 Sv) to the offshore regions, roughly half being directed to the North Pacific. These features are consistent with recent satellite imagery and in-situ observations, suggesting that diapycnal mixing within the vicinity of the Kuril Islands has a greater impact than was previously supposed on the Okhotsk Sea and the North Pacific. To examine this influence of tidal processes at the Kurils on circulations in the neighboring two basins, another numerical experiment was conducted using an ocean general circulation model with inclusion of tidal mixing along the islands, which gives a better representation of the Okhotsk Sea Mode Water than in the case without the tidal mixing. This is mainly attributed to the added effect of a significant upward salt flux into the surface layer due to tidal mixing in the Kuril Straits, which is subsequently transported to the interior region of the Okhotsk Sea. With a saline flux into the surface layer, cooling in winter in the northern part of the Okhotsk Sea can produce heavier water and thus enhance subduction, which is capable of reproducing a realistic Okhotsk Sea Mode Water. The associated low PV flux from the Kuril Straits to the open North Pacific excites the 2nd baroclinic-mode Kelvin and Rossby waves in addition to the 1st mode. Interestingly, the meridional overturning in the North Pacific is strengthened as a result of the dynamical adjustment caused by these waves, leading to a more realistic reproduction of the North Pacific Intermediate Water (NPIW) than in the case without tidal mixing. Accordingly, the joint effect of tidally-induced transport and transformation dominating in the Kuril Straits and subsequent eddy-transport is considered to play an important role in the ventilation of both the Okhotsk Sea and the North Pacific Ocean. This revised version was published online in July 2006 with corrections to the Cover Date.     

Validation of turbulence closure parameterisations for stably stratified flows using the PROVESS turbulence measurements in the North Sea

A number of parameterisations for the simulation of mixing processes in the thermocline are compared and tested against the microstructure data of the PROVESS campaigns, conducted in the northern part of the North Sea during the autumn of 1998. The transport term in the turbulent kinetic energy equation is parameterised via the introduction of a third stability function S_k for turbulent energy diffusion. The formulations are compared with a simpler scheme based upon limiting conditions for turbulence variables. Improved results are obtained with a new form of S_k. The best agreement is, however, found with the simpler limiting scheme. This is explained in terms of a turbulence length scale theory for stably stratified turbulence. In agreement with previous laboratory and ocean data it is found that the ratios of the Thorpe and Kolmogorov scales to the Ozmidov length scale approach critical limiting values in the thermocline. The first of these conditions is satisfied when limiting conditions are implemented into the scheme, providing the necessary minimum value for the dissipation rate, whereas the schemes without limiting conditions fail to produce this critical ratio. The basic reason for this failure is that the Thorpe scale is overestimated, which is shown to be connected to an even larger overprediction of the dissipation rate of temperature variance. To investigate the impact of non-resolved advective processes and salinity stratification on the turbulence predictions, additional numerical experiments were conducted using a simple scheme for data assimilation. The best agreement is found again with the limiting scheme, which is able to make reasonable predictions for the dissipation rate without knowing the detailed shape of the mean stratification profile. It is shown that advective transport due to tidally and wind-driven motions has a non-negligible impact on vertical mixing. This is seen in the data and the models by periodic enhancements of turbulent mixing inside the thermocline.     

Nutrient transport from an artificial upwelling of deep sea water

The transport of nutrient-rich, deep sea water from an artificial upwelling pipe has been simulated. A numerical model has been built within a commercial Computational Fluid Dynamics (CFD) package. The model considers the flow of the deep sea water after it is ejected from the pipe outlet in a negatively buoyant plume (densimetric Froude number = −2.6), within a stably stratified ocean environment subject to strong ocean current cross flow. Two cross-flow profiles were tested with momentum flux ratios equal to 0.92 and 3.7. The standard k-ε turbulence model has been employed and a range of turbulent Schmidt and Prandtl numbers tested. In all cases the results show that the rapid diffusion of heat and salinity at the pipe outlet causes the plume to attain neutral buoyancy very rapidly, preventing strong fountain-like behavior. At the higher momentum flux ratio fountain-like behavior is more pronounced close to the pipe outlet. The strong cross-current makes horizontal advection the dominant transport process downstream. The nutrient plume trajectory remains largely within one relatively thin stratified layer, making any ocean cross-flow profile less important. Very little unsteady behavior was observed. The results show that the nutrient is reduced to less than 2% of its inlet concentration 10 meters downstream of the inlet and this result is largely independent of turbulent Prandtl or Schmidt number. Initial results would suggest that if such an artificial upwelling were to be viable for an ocean farming project, a large number of upwelling pipes would be necessary. Further work will have to determine the minimum nutrient concentration required to sustain a viable phytoplankton population and the required spacing between upwelling pipes.     

Tidal sea mathematical modelling: Current applications and future developments

The purpose of this paper is to describe the contribution of the LNH/EDF, France, which work on the two following topics in sea modelling: (a) calibration and use of mathematical tools in order to perform numerical forecasts and simulation of sea patterns; and (b) important research effort, in order to develop ever more accurate numerical algorithms and to understand ever better the physical processes.Simulating marine flows means in any case filtering of the sea-phenomenon on both time and space scales. The use of such filters related to the preponderant scales in each study leads to parameterize the unsolved scales through specific behaviour laws such as turbulence and dispersion.Specific properties of sea modelling at LNH are based on the splitting-up technique applied to each term of the shallow water equations. The advective transport is solved through a Lagrangian method. The diffusive and propagative part is solved through two different algorithms depending on which method is used, a finite-difference method or a finite-element method.The current applications of tidal computations are numerous and various, the examples given in this paper concern on one hand the simulations of uncovering tidal flats in the estuary of the River Canche and on the other hand of storm surge and tide interaction in the North Sea and English Channel. Some future developments among others are relative to the improvement of the properties of the solver of the advective term, and the treatment of uncovering and overcovering tidal areas. About the mixing properties of fluid, some research work concerns the better understanding and modelling of the big vortex shedding behind a hindrance cross a flow; therefore an experiment is done at this time which will provide a data collection useful to the calibration of a 2D-turbulence model.Other simulations and research works are performed in fields connected to tidal current computations such as thermal impact or bed-load transport studies which lead to the problems of long term simulation and modelling of residual currents which are of a lower order of accuracy.In order to conclude, the necessity is stressed of a double effort, scientific on one hand and algorithmic on the other hand.