首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
海洋是多尺度强迫-耗散系统,机械能主要在大尺度输入,在小尺度耗散。在大、中尺度运动的能量向小尺度湍流传递过程中,内波扮演着重要角色。内波的生成和破碎可打破海洋动力平衡,而在陆架区,内波(主要是内孤立波)的浅化演变与耗散则是驱动湍流混合的关键过程。通过长期的理论、观测与数值模拟研究,目前已认识到内波浅化过程中主要发生如下演变:波形调制、极性转变、裂变、破碎与耗散。相较于直接发生破碎,浅化演变过程中的裂变及其引发的剪切不稳定和对流不稳定是内孤立波在陆架区的主要耗散机制,显著调制陆架区的跃层混合。从能量串级的角度讲,内孤立波浅化裂变为动力不稳定的高频内波是潮能串级的重要通道。本文简要回顾南海北部陆架区内波的研究历史,并着重总结内波在陆架区演变与耗散机制的研究进展。  相似文献   

2.
Airy waves have a sinusoidal profile in deep water that can be modeled by a time series at any point x and time t, given by η(x,t) = (Ho/2) cos[2πx/Lo − 2πt/Tw], where Ho is the deepwater height, Lo is the deepwater wavelength, and Tw is the wave period. However, as these waves approach the shore they change in form and dimension so that this equation becomes invalid. A method is presented to reconstruct the wave profile showing the correct wavelength, wave height, wave shape, and displacement of the water surface with respect to the still water level for any water depth.  相似文献   

3.
A nonlinear model for nonbreaking shoaling random waves   总被引:1,自引:0,他引:1  
AnonlinearmodelfornonbreakingshoalingrandomwavesLiuXin'an,HuangPeiji,ChenXueying,HuZejian(ReceivedOctober15,1996,acceptedAugu...  相似文献   

4.
The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.  相似文献   

5.
《Coastal Engineering》2006,53(4):311-318
The extended mild-slope equations of Suh et al. [Suh, K.D., Lee, C., Park, W.S., 1997. Time-dependent equations for wave propagation on rapidly varying topography. Coastal Eng., 32, 91–117] and Lee et al. [Lee, C., Kim, G., Suh, K.D., 2003. Extended mild-slope equation for random waves. Coastal Eng., 48, 277–287] are compared analytically and numerically to determine their applicability to random wave transformation. The geometric optics approach is used to compare the two models analytically. In the model of Suh et al., the wave number of the component wave with a local angular frequency ω is approximated with an accuracy of O(ω  ω¯) at a constant water depth, where ω¯ is the carrier frequency of random waves. In the model of Suh et al., however, the diffraction effects and higher-order bottom effects are considered only for monochromatic waves, and the shoaling coefficient of random waves is not accurately approximated. This inaccuracy arises because the model of Suh et al. was derived for regular waves. In the model of Lee et al., all the parameters of random waves such as wave number, shoaling coefficient, diffraction effects, and higher-order bottom effects are approximated with an accuracy of O(ω  ω¯). This approximation is because the model of Lee et al. was developed using the Taylor series expansion technique for random waves. The result of dispersion relation analysis suggests the use of the peak and weighted-average frequencies as a carrier frequency for Suh et al. and Lee et al. models, respectively. All the analytical results are verified by numerical experiments of shoaling of random waves over a slightly inclined bed and diffraction of random waves through a breakwater gap on a flat bottom.  相似文献   

6.
The three-dimensional numerical model with σ-coordinate transformation in the vertical direction is applied to the simulation of surface water waves and wave-induced laminar boundary layers. Unlike most of the previous investigations that solved the simplified one-dimensional boundary layer equation of motion and neglected the interaction between boundary layer and outside flow, the present model solves the full Navier–Stokes equations (NSE) in the entire domain from bottom to free surface. A non-uniform mesh system is used in the vertical direction to resolve the thin boundary layer. Linear wave, Stokes wave, cnoidal wave and solitary wave are considered. The numerical results are compared to analytical solutions and available experimental data. The numerical results agree favorably to all of the experimental data. It is found that the analytical solutions are accurate for both linear wave and Stokes wave but inadequate for cnoidal wave or solitary wave. The possible reason is that the existing analytical solutions for cnoidal and solitary waves adopt the first-order approximation for free stream velocity and thus overestimate the near bottom velocity. Besides velocity, the present model also provides accurate results for wave-induced bed shear stress.  相似文献   

7.
《Coastal Engineering》1998,35(3):185-209
Two depth inversion algorithms (DIA) applicable to coastal waters are developed, calibrated, and validated based on results of computations of periodic waves shoaling over mild slopes, in a two-dimensional numerical wave tank based on fully nonlinear potential flow (FNPF) theory. In actual field situations, these algorithms would be used to predict the cross-shore depth variation h based on sets of values of wave celerity c and length L, and either wave height H or left–right asymmetry s2/s1, simultaneously measured at a number of locations in the direction of wave propagation, e.g., using video or radar remote sensing techniques. In these DIAs, an empirical relationship, calibrated for a series of computations in the numerical wave tank, is used to express c as a function of relative depth koh and deep water steepness koHo. To carry out depth inversion, wave period is first predicted as the mean of observed L/c values, and Ho is then predicted, either based on observed H or s2/s1 values. The celerity relationship is finally inverted to predict depth h. The algorithms are validated by applying them to results of computations for cases with more complex bottom topography and different incident waves than in the original calibration computations. In all cases, root-mean-square (rms)-errors for the depth predictions are found to be less than a few percent, whereas depth predictions based on the linear dispersion relationship—which is still the basis for many state-of-the-art DIAs—have rms-errors 5 to 10 times larger.  相似文献   

8.
General characteristics of topographically trapped subinertia waves are discussed from the viewpoint of an eigenvalue problem and ray theory. Special attention is paid to the slope parameterS(x) (=(dh/dx)/h, wherex denotes the coordinate perpendicular to the shoreline, increasing seaward, andh(x) is the depth) which is a measure of the strength of the restoring force of the waves. Three cases for theS distribution are considered, in whichS is assumed to be positive at the coast and to tend to zero far from the coast. The first is whereS(x) decreases monotonically towards the open ocean. It is found in this case that waves are trapped near the coast. The second is whereS(x) does not decrease monotonically, but has a maximum. It is concluded that this case may contain two types of waves, i.e., those trapped near the coast and those trapped near the maximum, and the dispersion curves corresponding to different types may nearly intersect, namely, result in “kissing”. The third is whereS(x) has a negative region (corresponding to the presence of a trench). It is found in this case that an infinite sequence of waves is trapped in the negativeS region which propagate with the coast to their left (right) in the northern (southern) hemisphere besides the waves trapped near the coast. The topography in the second case corresponds to a typical continental shelf and a typical continental slope. It is shown by model calculation that trapped waves are present over the continental slope as well as over the continental shelf. Almost the same results are obtained for superinertia waves ifS is replaced by 1/h which is a measure of the restoring force of superinertia waves.  相似文献   

9.
Adjustment of Wind Waves to Sudden Changes of Wind Speed   总被引:1,自引:0,他引:1  
An experiment was conducted in a small wind-wave facility at the Ocean Engineering Laboratory, California, to address the following question: when the wind speed changes rapidly, how quickly and in what manner do the short wind waves respond? To answer this question we have produced a very rapid change in wind speed between U low (4.6 m s?1) and U high (7.1 m s?1). Water surface elevation and air turbulence were monitored up to a fetch of 5.5 m. The cycle of increasing and decreasing wind speed was repeated 20 times to assure statistical accuracy in the measurement by taking an ensemble mean. In this way, we were able to study in detail the processes by which the young laboratory wind waves adjust to wind speed perturbations. We found that the wind-wave response occurs over two time scales determined by local equilibrium adjustment and fetch adjustment, Δt 1/T = O(10) and Δt 2/T = O(100), respectively, in the current tank. The steady state is characterized by a constant non-dimensional wave height (H/gT 2 or equivalently, the wave steepness for linear gravity waves) depending on wind speed. This equilibrium state was found in our non-steady experiments to apply at all fetches, even during the long transition to steady state, but only after a short initial relaxation Δt 1/T of O(10) following a sudden change in wind speed. The complete transition to the new steady state takes much longer, Δt 2/T of O(100) at the largest fetch, during which time energy propagates over the entire fetch along the rays (dx/dt = c g) and grows under the influence of wind pumping. At the same time, frequency downshifts. Although the current study is limited in scale variations, we believe that the suggestion that the two adjustment time scales are related to local equilibrium adjustment and fetch adjustment is also applicable to the ocean.  相似文献   

10.
南海北部陆架区内孤立波向岸传播过程研究   总被引:1,自引:0,他引:1  
南海北部是全球海洋中内孤立波最强和最为活跃的海域。然而,内孤立波在传入陆架区后,其形态发生显著变化,其传播演变过程表现出高度的复杂性。本研究综合卫星图像和数值模式手段研究了内孤立波在向岸传播过程中的空间变化特征。可见光卫星图像研究结果显示,南海北部陆架区存在三种形态的内孤立波,分别为第一模态下凹型内孤立波、第一模态上凸型内孤立波和第二模态内孤立波。受水深和层结变化的控制,它们的分布区域显著不同。基于MITgcm的数值模拟研究表明,上凸型内孤立波由第一模态下凹内孤立波经过极性转换过程发展而来,而第二模态内孤立波由第一模态下凹内孤立波与急剧变浅地形相互作用而产生。  相似文献   

11.
An existing 2D time-domain method for separating irregular incident and reflected waves by wavelet transform [Ma et al., 2010. A new method for separation of 2D incident and reflected waves by the Morlet wavelet transform. Coastal Eng., 57(6):597–603] is extended to account for obliquely incident irregular waves propagating over sloping bottoms. The linear shoaling and refraction coefficients are adopted to determine the amplitude and phase changes of waves. The optimal central frequency of the Morlet wavelet is determined by the minimum Shannon wavelet entropy. Numerical tests show that the present method can accurately separate waves over horizontal depths. For waves at sloping bottoms, however, the separation errors increase as bottom slope increases and are significant for waves with incident angle larger than π/3.  相似文献   

12.
A high resolution modeling study is undertaken, with a 2.5-dimensional nonhydrostatic model, of the generation of internal waves induced by tidal motion over the ridges in Luzon Strait. The model is forced by the barotropic tidal components K1, M2, and O1. These tidal components, along with the initial density field, were extracted from data and models. As the barotropic tide moves over the Luzon Strait sills, there is a conversion of barotropic tidal energy into baroclinic tidal energy. Depressions are generated that propagate towards the Asian Seas International Acoustics Experiment (ASIAEX) test site on the Chinese continental shelf. Nonlinear effects steepen the depressions, frequency and amplitude dispersion set in, and disintegration into large amplitude solitary waves occurs. The effects of varying the initial density field, tidal component magnitudes, as well as adding a steady background current to represent the occasional excursions of the Kuroshio Current into the strait, are considered.Depressions are generated at each of the two sills in Luzon Strait which radiate away, steepening and evolving into internal solitary wave trains. Baroclinic fluxes of available potential energy, kinetic energy and linear are calculated for various parameter combinations. The solitary wave trains produced in the simulations generally consist of large amplitude wave trains alternating with small amplitude wave trains. During strong tidal flow, Kelvin–Helmholtz type instabilities can develop over the taller double-humped sill. The solitary waves propagating towards the ASIAEX test site have been observed to reach amplitudes of 120–250 m, depending on the tidal strength. ASIAEX observations indicate amplitudes up to 150 m and the Windy Island Experiment (WISE) measurements contain magnitudes over 200 m. The model results yield solitary wave amplitudes of 70–300 m and half widths of 0.60–3.25 km, depending on parameter values. These are in the range of observations. Measurements by Klymak et al. (2006), in the South China Sea, exhibit amplitudes of 170 m, half widths of 3 km and phase speeds of 2.9 m s?1. Model predictions indicate that the solitary waves making up the wave packet each experience different background currents with strong near surface shear.The energy in the leading soliton of the large amplitude wave trains ranges between 1.8 and 9.0 GJ m?1. The smaller value, produced using barotropic tidal currents based on the Oregon State University data base, is the same as the energy estimated to be in a solitary wave observed by Klymak et al. (2006). Estimates of the conversion of barotropic tidal energy into radiating internal wave energy yield conversion rates ranging between 3.6% and 8.3%.  相似文献   

13.
《Oceanologica Acta》2002,25(2):51-60
A new composite model, which consists of a generation model of the internal tides and a regularized long wave propagation model, is presented to study the generation and evolution of internal solitary waves in the sill strait. Internal bores in the sill strait are first simulated by the generation model, and then the internal tidal field outside of the sill region is given as input for the propagation model. Numerical experiments are carried out to study the imposing tide, depth profile, channel width and shoaling effect, etc., on the generation and evolution of internal solitary waves. It is shown that only when the amplitude of internal tide at the forcing boundary of the propagation model is large enough that a train of internal solitary waves would be induced. The amplitude of the imposing tide in the generation model, shoaling effect, asymmetry of the depth profile and channel width have some effects on the amplitude of the induced internal solitary wave. The imposing tidal flow superimposed on a constant mean background flow has a great damping effect on the induced internal waves, especially on those propagate against the background flow direction. The generation and propagation of internal solitary waves in three possible straits among the Luzon Strait are simulated, and the reasons for the asymmetry of their propagation are also explained.  相似文献   

14.
The transformation of cylindrical solitary waves is studied theoretically and numerically in the context of the Gardner equation. The amplitudes of solitons are determined as functions of the ray coordinate for different combinations of external parameters.  相似文献   

15.
16.
Solitary waves have been commonly used as an initial condition in the experimental and numerical modelling of tsunamis for decades. However, the main component of a tsunami waves acts at completely different spatial and temporal scales than solitary waves. Thus, use of solitary waves as approximation of a tsunami wave may not yield realistic model results, especially in the coastal region where the shoaling effect restrains the development of the tsunami wave. Alternatively, N-shaped waves may be used to give a more realistic approximation of the tsunami wave profile. Based on the superposition of the sech2(*) waves, the observed tsunami wave profile could be approximated with the N-shaped wave method, and this paper presents numerical simulation results based on the tsunami-like wave generated based on the observed tsunami wave profile measured in the Tohoku tsunami. This tsunami-like wave was numerically generated with an internal wave source method based on the two-phase incompressible flow model with a Volume of Fluid (VOF) method to capture the free surface, and a finite volume scheme was used to solve all the governing equations. The model is first validated for the case of a solitary wave propagating within a straight channel, by comparing its analytical solutions to model results. Further, model comparisons between the solitary and tsunami-like wave are then made for (a) the simulation of wave run-up on shore and (b) wave transport over breakwater. Comparisons show that use of these largely different waveform shapes as inputs produces significant differences in overall wave evolution, hydrodynamic load characteristics as well as velocity and vortex fields. Further, it was found that the solitary wave uses underestimated the total energy and hence underestimated the run-up distance.  相似文献   

17.
A fully nonlinear Boussinessq-type model with several free coefficients is considered as a departure point. The model is monolayer and low order so as to simplify numerical solvability. The coefficients of the model are here considered functions of the local water depth. In doing so, we allow to improve the dispersive and shoaling properties for narrow banded wave trains in very deep waters. In particular, for monochromatic waves the dispersion and shoaling errors are bounded by ~ 2.8% up to kh = 100, being k the wave number and h the water depth. The proposed model is fully nonlinear in weakly dispersive conditions, so that nonlinear wave decomposition in shallower waters is well reproduced. The model equations are numerically solved using a fourth order scheme and tested against analytical solutions and experimental data.  相似文献   

18.
Surface gravity waves are commonly observed to slow down and to stop at a beach without any noticeable reflection taking place. We assume that as a consequence the waves are continuously giving up their linear and angular momenta, which they carry with them, along with energy, as they propagate into gradually decreasing mean depths of water. It takes a force to cause a time rate of decrease in the linear momentum and a torque to produce a time rate of decrease in the angular momentum. Both a force and a torque operate on the shoaling waves, due to the presence of the sloping bottom, to cause the diminution of their linear and angular momenta. By Newton’s third law, action equals reaction, an equal but opposite force and torque are exerted on the bottom. No other mechanisms for transferring linear and angular momenta are included in the model. Since the force on the waves acts over a horizontal distance during shoaling, work is done on the waves and energy flux is not conserved. Bottom friction, wave interaction with a mean flow, scattering from small-scale bottom irregularities and set-up are neglected. Mass flux is conserved, which leads to a shoreward monotonic decrease in amplitude consistent with available swell data. The formula for the time-independent force on the bottom agrees qualitatively with observations in seven different ways: four for swell attenuation and three for sediment transport on beaches. Ardhuin (2006) argues against a mean force on the bottom that is not hydrostatic, mainly by using conservation of energy flux. He also applies the action balance equation to shoaling waves. Action is a difficult concept to grasp for motion in a continuum; it cannot be easily visualized, and it is not really necessary for solving the shoaling wave problem. We prefer angular momentum because it is clearly related to the observed orbital motion of the fluid particles in progressive surface waves. The physical significance of wave action for surface waves has been described recently by showing that in deep water action is equivalent to the magnitude of the wave’s orbital angular momentum (Kenyon and Sheres, 1996). Finally, Ardhuin requires that there be a significant exchange of linear momentum between shoaling waves and an unspecified mean flow, although the magnitude and direction of the exchange are not predicted. No mention is made of what happens to the orbital angular momentum during shoaling. Mass flux conservation is not stated.  相似文献   

19.
20.
An analysis of the spatial and temporal variability of the field of mechanical energy transfer (MET) from the atmosphere into the ocean is based on a separate numerical simulation of evolution for the terms of source function for a wind-wave model conducted in the Indian Ocean area for the period 1998 to 2009. The MET field is described by two integral values calculated per unit area: the total rate of energy flux from the wind to waves, I E (x, t), and the rate of energy-loss flux for the wind waves, D E (x, t). To solve this problem, the wind field W(x, t) is used, downloaded from the NCEP/NOAA archive [1], and the fields I E (x, t) and D E (x, t) were calculated using the numerical model WAM [2] with the modified source function proposed in [3]. Maps for the fields I E (x, t) and D E (x, t) were obtained by calculations with different scales of the space-time averaging, extreme and average values of the MET were found, seasonal and interannual variability was estimated, and the 12-year trend for several mean quantities was obtained.  相似文献   

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

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