首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Previous studies of combined wave and current bottom boundary layer flow have concentrated on the final converged state of the flow following the addition of waves to a current. While this final state is of primary interest to modellers and engineers, it pre-supposes that such a state is actually attained in reality, and this may not always be the case. In addition, it overlooks the interesting and complicated transient effects which occur as a wave-current flow evolves from one state to another. The present study concentrates attention on the transient effects predicted by a “one-equation” turbulence closure model. Results of case studies are presented in which waves are superimposed co-linearly on a current (“forward problem”), and are then removed from the converged wave-current flow (“backward problem”). Two formulations of the “forward” and “backward” problems are discussed. In the first the steady component of the pressure gradient driving the mean flow is held constant throughout, and in the second the steady component of the mass flux is held constant. In each case the detailed evolution of the profiles of mean velocity, turbulent energy, mixing length, eddy viscosity and shear stress are discussed. More generally, the question of the convergence timescale of a combined wave-current flow is considered, and a convergence criterion is proposed.  相似文献   

2.
Pengzhi Lin  C. W. Li   《Ocean Engineering》2003,30(7):855-876
A three-dimensional numerical model is developed in this study to investigate the problem of wave–current–body interaction. The model solves the spatially averaged Navier–Stokes equations. Turbulence effects are modeled by a subgrid-scale (SGS) model using the concept of large eddy simulation (LES). The model is employed to study the wave–current interaction with a square cylinder that is mounted on the bottom and vertically pierces the free surface. The force analysis demonstrates that the presence of waves can reduce both the strength and frequency of vortex shedding induced by a uniform current due to the nonlinear wave–current interaction. The free surface elevation, strain rates of the mean flow, and eddy viscosity are found to closely correlate with the mechanism of vortex shedding. It is also shown that when the vortex shedding is neglected in the calculation such as by the potential flow approach, one may significantly underestimate the magnitude of in-line force. The energy spectral analysis reveals that there exist initiating, growing, and decaying regions for shedding vortices around the cylinder. In the vortex initiating region, both coherent and turbulent structures are nearly two-dimensional that become three-dimensional in the vortex growing region. The kinetic energy of both coherent and turbulent motions is dissipated in the vortex decaying region, within which the mean flow gradually returns back to two-dimensional.  相似文献   

3.
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.  相似文献   

4.
An analytical theory which describes the motion in a turbulent wave boundary layer near a rough sea bottom by using a two-layer time invariant eddy viscosity model is presented. The eddy viscosity in the inner layer increases quadratically with the height above the sea bottom. In the outer layer the eddy viscosity is taken as a constant. The mean velocity and shear stress profiles, the bottom shear stress and the bottom friction coefficient are presented, and comparisons are made with experimental results.  相似文献   

5.
An approximate steady solution of the wave-modified Ekman current is presented for gradually varying eddy viscosity by using the WKB method with the variation of parameters technique. The parameters involved in the solution can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water. The solution reduces to the exact solution when the eddy viscosity is taken as a constant. As illustrative examples, for a fully developed wind-generated sea with different wind speeds and a few proposed gradually varying eddy viscosities, the current profiles calculated from the approximate solutions are compared with those of the exact solutions or numerical ones by using the Donelan and Pierson wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. It is shown that the approximate solution presented has an elegant form and yet would be valid for any given gradually varying eddy viscosity. The applicability of the solution method to the real ocean is discussed following the comparisons with published observational data and with the results from a large eddy simulation of the Ekman layer.  相似文献   

6.
Scaling of the equations of motion of the Antarctic Circumpolar Current indicates that the Rossby number and the Ekman number are 10−4 to 10−5 but the vertical Ekman number may reach unity in the bottom boundary layer. The equations of motion are integrated vertically from the surface to the bottom and averaged over a latitude circle. The resulting equation in the meridional direction is predominantly geostrophic, whereas the main terms of the equation in the zonal direction are the wind stress and the bottom stress. When the vertical eddy viscosity near the bottom is of the order of 102cm2/sec, the total zonal transport through the Drake Passage computed from the balance of the wind stress and the bottom stress equals 260×106m3/sec, the amount determined byReid andNowlin (1970) from observations. The northward transport reduces the eastward transport corresponding to the wind stress of the westerlies in the A. C. C. through the Coriolis' term in the vertically integrated equation of motion of the zonal direction. South of the Drake Passage, such reduction reaches about ten percent of the wind-driven transport mainly due to the peripheral water discharge. North of the Drake Passage, the northward transport may be generated by the effect of the South American coast which prevents free eastward movement of the A. C. C., causing a wake to the east. This transport may contribute to a part of the northward transport of the bottom water postulated byMunk (1966). The effect of the horizontal eddy viscosity in the zonal transport equation is negligible except near the Antarctic coast, if the eddy viscosity is less than 109cm2/sec.  相似文献   

7.
A three-dimensional tidal current model is developed and applied to the East China Sea (ECS), the Yellow Sea and the Bohai Sea. The model well reproduces the major four tides, namely M2, S2, K1 and O1 tides, and their currents. The horizontal distributions of the major four tidal currents are the same as those calculated by the horizontal two-dimensional models. With its high resolutions in the horizontal (12.5 km) and the vertical (20 layers), the model is used to investigate the vertical distributions of tidal current. Four vertical eddy viscosity models are used in the numerical experiments. As the tidal current becomes strong, its vertical shear becomes large and its vertical profile becomes sensitive to the vertical eddy viscosity. As a conclusion, the HU (a) model (Davieset al., 1997), which relates the vertical eddy viscosity to the water depth and depth mean velocity, gives the closest results to the observed data. The reproduction of the amphidromic point of M2 tide in Liaodong Bay is discussed and it is concluded that it depends on the bottom friction stress. The model reproduces a unique vertical profile of tidal current in the Yellow Sea, which is also found in the observed data. The reason for the reproduction of such a unique profile in the model is investigated.  相似文献   

8.
9.
黄大吉 《海洋与湖沼》1993,24(4):385-392
用一垂向二维的诊断模式,研究了浅水区锋面的流场结构。参照James模式并结合杭州湾具体情形来确定垂向涡动粘性系数,得到了与实测结果相一致的浅水区(杭州湾)锋面的流场结构。研究表明,在探讨浅水区锋面流场时必须选取合理的模式;而深水锋面流场结构对所采用的模式依赖性不大;锋面的流场结构不仅与锋面的强度有关,而且与锋面的结构有关。  相似文献   

10.
11.
The seakeeping characteristics of various boom geometries in irregular waves and currents are investigated. The response of a floating boom section on the open sea is a function of a number of parameters, such as boom geometry, distribution of mass, buoyancy/weight ratio, and wave and current characteristics. To understand the relationship between these design parameters more clearly, a series of regular and irregular wave tests were conducted with six different 1:4 scale models for three current velocities and six different wave conditions. To simplify the problem, only rigid boom sections consisting of a buoyancy cylinder and vertical skirt were used. In parallel with this experimental program, a numerical model for the responses of two-dimensional floating boom sections in small-amplitude waves is also developed. The numerical results are compared with our large-scale experimental results. The boom effectiveness on the open sea is evaluated based on the concept of “effective draft” and “effective freeboard” assuming that drainage and oversplashing failures are the prime mechanisms of containment failure. Using the present results, a guideline for the optimum design/selection of future booms is developed.  相似文献   

12.
超浅海风暴潮模型提出后[2],对渤海风潮,作为超浅海问题,进行了数值研究[1]。其结果的分析和观测资料的比较都表明了该模型有一定的应用价值;故,对超浅海风暴潮模型作进一步的探讨是有一定意义的。尤其因为我国是一个多浅水域和多风暴潮的国家,这种研究就具有更重要的意义。  相似文献   

13.
Analysis of measurements from two long-term moored arrays in and near the Gulf Stream suggests a simple parameterization of eddy spatial covariance statistics: a parameterization that can be referred to as “quasi-homogeneous and isotropic”. Taking the normalized covariance function (i.e. the correlation function) for streamfunction to be homogeneous and isotropic and assuming motions to be horizontally nondivergent and hydrostatic permit the velocity and temperature covariances to be derived from the streamfunction covariance. Statistical tests indicate that deviations from these assumptions are indistinguishable from Gaussian random noise. The spatial correlation function used in Gaussian with a decay scale of about 140 km, which is only weakly depth dependent. A simple form is also suggested for the vertical lag dependence. This parameterization permits calculation of derived quantities such as the eddy vorticity flux divergence which is discussed in the context of the mean potential vorticity balances for the depth integrated circulation and for the subthermocline layer. The divergence of the relative vorticity flux is found capable of driving two counter-rotating gyres of strength 30–40 Sv on either side of the Stream, as are observed. The “thickness flux” dominates the lower layer eddy potential vorticity flux and is of the correct sign to make the recirculation more barotropic. The lower layer eddy forcing is weak and the gyres exist in a region of nearly uniform mean potential vorticity.  相似文献   

14.
Morphodynamics of a bar-trough surf zone   总被引:3,自引:0,他引:3  
A field study was made of the distinguishing morphodynamic processes operating in a surf zone which perennially exhibits accentuated bar-trough topography (the “longshore-bar-trough” and “rhytmic-bar-and-beach” states as described by Wright and Short, 1984). Characteristic features of the morphology include a shallow bar with a steep shoreward face, a deep trough, and a steep beach face. This morphology, which is favored by moderate breaker heights and small tidal ranges, strongly controls the coupled suite of hydrodynamic processes. In contrast to fully dissipative surf zones, the bar-trough surf zone is not at all saturated and oscillations at incident wave frequency remain dominant from the break point to the subaerial beach. The degree of incident wave groupiness does not change appreciably across the surf zone. Infragravity standing waves which, in dissipative surf zones, dominate the inshore energy, remain energetically secondary and occur at higher frequencies in the bar trough surf zone. Analyses of the field data combined with numerical simulations of leaky mode and edge wave nodal—antinodal positions over observed surf-zone profiles, indicate that the frequencies which prevail are favored by the resonant condition of antinodes over the bar and nodes in the trough. Standing waves which would have nodes over the bar are suppressed. Sediment resuspension in the surf zone appears to be largely attributable to the incident waves which are the main source of bed shear stress. In addition, the extra near-bottom eddy viscosity provided by the reformed, non-breaking waves traversing the trough significantly affects the vertical velocity profile of the longshore current. Whereas the bar is highly mobile in terms of onshore—offshore migration rates, the beach face and inner regions of the trough are remarkably stable over time.  相似文献   

15.
The aim of the present paper is to explain some of the differences between previously published analytical and numerical models of combined wave and current bottom boundary layer flow. To this end, the Grant and Madsen (1979) model for wave-current, rough turbulent flow is modified to include both first and second harmonic time variations in the eddy viscosity (K). The functional form of the coefficients controlling the amount of time variation is established by analysing the numerical model results of Davies (1990). The addition of time variation in K reduces the strong non-linearity exhibited by the mean stress in the original Grant and Madsen model for current dominated cases, and reproduces the veering of the current predicted by numerical turbulence closure models.  相似文献   

16.
We have investigated the three-dimensional Lagrangian motion of water particles related with tidal exchange between two basins with a constant depth connected through a narrow strait and the effects of density stratification on the exchange processes by tracking a number of the labeled particles. Tide-induced transient eddies (TITEs), which are similar to those in two-dimensional basin, are generated behind the headlands. Upwelling appears around the center of the eddy and sinking around the boundary. When the basins are filled with homogeneous water, a pair of vortices are produced in the vertical cross section of the strait due to bottom stress, with upwellings along the side walls of the strait and sinking in the center of the strait. These circulations form the horizontally convergent field in the cross-strait direction in the upper layers while the horizontal divergence takes place in the bottom layer. These vertical water-motions produce the three-dimensional distribution of velocity shear and phase lag of the tidal current around the strait, and the Lagrangian drifts of water particles become large. As a result, water exchange through the strait is greatly enhanced: The water exchange rate reaches 94.1% which is much larger than that obtained in the vertically integrated two-dimensional model. When the basins are stratified, the stable stratification suppresses the vertical motion so that a pair of vertical vortices are confined in the lower layers. This leads to a decrease in the exchange rate, down to 88.6%. Our numerical results show that the three-dimensional structure of tidal currents should be taken into account in tidal exchange through a narrow strait.  相似文献   

17.
A theory is outlined for time-dependent currents induced near the sea surface in deep water, away from coastal boundaries, by a variable wind stress and deep-water wave field. It is based on the theory of Weber (1983) which uses a second-order perturbation expansion of the Navier-Stokes equations in Lagrangian coordinates and includes the Coriolis effect. It uses an eddy viscosity formulation for both wave dissipation and momentum transfer within the current field: the eddy viscosity ν may be allowed to vary with depth. The wind stress may be time-varying and the wave field may vary in both space and time.For the case of a constant ν, the results agree with those of Ursell (1950), Hasselmann (1970) and Pollard (1970) in the limit ν→0, and the steady-state results agree with those of Weber. For a particular case of depth-varying ν, results (obtained from numerical simulations) are in better general agreement with observations of wind-induced surface drift than when a constant ν is used.An outline is given of the application of the theory to the case of a random sea state. There are good prospects for using output data from numerical wave prediction models to drive the equations of this near-surface current model.  相似文献   

18.
纪艳菊  刘淑波  齐震 《海洋科学》2014,38(12):120-127
本文通过假定底边界层湍黏性的三次多项式参数化形式,基于简化的Navier–Stokes方程,并利用超几何方程的性质,推导出了湍流粗糙底边界层的速度解析解。同时,得到了底边界层内其他的动力参数,如底剪应力、Ekman传输、Ekman抽吸及近底部速度分布场,从理论上讨论了均匀混合底边界层特征量分布特征。通过数值结果分析,进一步得出底边界层的总速度、亏损速度及其剪应力受平均流的角频率和地球自转影响比较大;而底边界层的动力结构对于底边界层顶部粗糙度不敏感。该涡黏性模式从理论上丰富了底边界层涡黏性的形式,为底边界层的动力系统研究提供了借鉴和理论参考。  相似文献   

19.
A new model is presented for the propagation of monochromatic surface waves over a region of arbitrary, one-dimensional bottom topography. The smoothly varying bed profile is divided into a series of shelves separated by abrupt steps. The wave fields on either side of each step are related by a “transfer matrix”, and the propagation of waves along the shelf between adjacent steps is described by a “rotation matrix”. Starting from a point where the surface wavefield is known, the step by step application of the appropriate combination of these matrices allows computation of the wavefield over the region of interest. If the individual steps are small then the transfer matrix reduces to a simpler plane-wave form, with considerable savings in computational effort. Comparisons are made with an exact potential solution for single and double steps in order to investigate the accuracy and validity of the matrix method. Finally, the model is applied to the case of wave reflection by fixed sinusoidal bottom undulations, and good agreement is found between its predictions and existing laboratory data.  相似文献   

20.
The normal force coefficient on a flat planing surface having arbitrary heave and pitch motion in two-dimensional flow is compared with the lift coefficient of a thin wing in an infinite fluid. Despite the totally different derivations, they are found to be identical (at large Froude numbers and low trim angles and allowing for the wing's interaction with twice as much fluid) at low reduced frequencies. For higher frequency motions, the wing's angle of attack induced lift and its pitch and heave damping are less than those of a planing surface, but the acceleration terms remain identical. The differences at the higher reduced frequencies are due to the fact that, in invisad irrotational flow, the planning plate cannot leave a vortex wake, whereas a wing does.It seems to follow that the “virtual mass” planing hull analysis can be applied to “quasi-static” problems involving wings and bodies in an infinite fluid without the slenderness restriction originally imposed by Jones (1946). Certainly, it is remarkable that the so called “quasi-steady” forces on a two-dimensional wing can be obtained in a few lines of elementary analysis. On the other hand, the method fails entirely when used to compute the pitching moment on a two-dimensional plate, even though it has been found to give good results for the three-dimensional case (Payne, 1981c).This work is offered as a very incomplete study of an intriguing relationship between two very different bodies of analysis. Much more work will need to be done before the relationship between the two approaches will be fully understood.  相似文献   

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

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