共查询到20条相似文献，搜索用时 15 毫秒
1.
In this paper, the evolution of focused waves using different paddle displacements (piston type) under laboratory conditions is presented. It is well known that in intermediate water depths, linear paddle displacements will generate spurious, free, sub and super harmonics. Thus, a second order correction to suppress these spurious free sub and super harmonics was used to generate the focused waves. The focused waves were generated in the laboratory using a linear superimposition principle, in which the wave paddle displacement is derived based on the sum of a number of sinusoidal components at discrete frequencies, whose phases are accordingly set to focus at a particular location. For this method of generation, the second order wave maker theory proposed by Schäffer [24] can be easily adopted and was used in the present study. Two different centre frequencies (fc = 0.68 Hz and 1.08 Hz) with three different bandwidth ratios (Δf/fc = 0.5, 0.75 and 1.0) were tested in a constant water depth, to consider both narrow and broadband spectra. These test cases correspond to wave focusing packets propagating in intermediate and deep water regions. Further, for each wave packet, two different amplitudes were considered in order to analyze non-breaking and breaking cases. Thus, by systematically generating the wave packets using the linear and second order paddle displacements, the analysis was carried out for the spectral and temporal evolution of selected long waves. The temporal evolution of the selected harmonics was analyzed using the Inverse Fast Fourier Transform (IFFT), to show the propagation of the spurious, free, long waves. Further, the variations in energy for the lower, higher and primary frequency ranges are reported for different wave paddle displacements. The analysis revealed that for the broadband spectrum the differences are more pronounced when using linear paddle displacements. We have also noticed a shift in focusing/breaking location and time (i.e. premature) due to the increase in crest height using linear displacements. The experiment data used in this paper has been provided as a supplementary, which can be used to validate the numerical models. 相似文献
2.
The formation time of alongshore morphological variability in surf zone sand bars has long been known to differ from one beach to the other and from one post-storm period to another. Here we investigate whether the type of sea state, i.e. distant swell waves or locally generated short period wind sea, affects the formation time of the emerging alongshore topographic variability.A numerical modeling approach is used to examine the emergence of alongshore variability under different shore-normal wave forcing. A research version of Delft3D, operating on the time-scale of wave groups, is applied to a schematised bathymetry with a single bar. The model is then used to investigate several wave scenarios, examining the impact of peak period, frequency spread and directional spread on the formation time of alongshore variability.Results show that an increase in wave period has a large effect, changing the formation time up to O (250%) in case the wave period is changed from a representative value for the Dutch coast (Tp ~ 5–6 s) to an Australian South East coast value (Tp ~ 10–12 s). In contrast, modifications in the directional and frequency spread of the wave field result only in a minor change in the formation time.Examination of hydrodynamics and potential sediment transport shows that the variations in formation time are primarily related to changes in the magnitude of the time-averaged flow conditions. Variations in the magnitude of very low frequency (f < 0.004 Hz) or infragravity (0.004 < f < 0.04 Hz) surf zone flow velocities do not affect the mean sediment transport capacity. Consequently the formation speed of patterns is primarily governed by positive feedback between mean flow and morphology, and low frequency flow fluctuations are of minor importance.These findings indicate that the development of alongshore topographic variability may be faster at swell dominated open coasts, primarily due to the occurrence of longer period swell. Also, at a given site, the arrival of a long wave period swell after a storm can accelerate the emergence of variability. 相似文献
3.
We have employed laboratory and numerical experiments in order to investigate propagation of waves in both long and short-crested wave fields in deep water. For long-crested waves with steepness, ϵ = kcac = 0.1 (a fairly extreme case), reliable prediction can be performed with the modified nonlinear Schrödinger equation up to about 40 characteristic wavelengths. For short-crested waves the accuracy of prediction is strongly reduced with increasing directional spread. 相似文献
4.
A lift based cycloidal wave energy converter (WEC) was investigated using potential flow numerical simulations in combination with viscous loss estimates based on published hydrofoil data. This type of wave energy converter consists of a shaft with one or more hydrofoils attached eccentrically at a radius. The main shaft is aligned parallel to the wave crests and submerged at a fixed depth. The operation of the WEC as a wave-to-shaft energy converter interacting with straight crested waves was estimated for an actual ocean wave climate. The climate chosen was the climate recorded by a buoy off the north-east shore of Oahu/Hawaii, which was a typical moderate wave climate featuring an average annual wave power PW = 17 kWh/m of wave crest. The impact of the design variables radius, chord, span and maximum generator power on the average annual shaft energy yield, capacity factor and power production time fraction were explored. In the selected wave climate, a radius R = 5 m, chord C = 5 m and span of S = 60 m along with a maximum generator power of PG = 1.25 MW were found to be optimal in terms of annual shaft energy yield. At the design point, the CycWEC achieved a wave-to-shaft power efficiency of 70%. In the annual average, 40% of the incoming wave energy was converted to shaft energy, and a capacity factor of 42% was achieved. These numbers exceeded the typical performance of competing renewables like wind power, and demonstrated that the WEC was able to convert wave energy to shaft energy efficiently for a range of wave periods and wave heights as encountered in a typical wave climate. 相似文献
5.
The main purpose of this study is to establish a better understanding of the relationship between drag reduction and surface roughness. Experiments were conducted to measure the force and flow characteristics of a circular cylinder with different types of artificial surface roughness over the range 6 × 103 < Re < 8 × 104 (Re is based on the cylinder diameter D). The roughness cylinder was formed by covering the exterior surface of the cylinder with uniformly distributed (1) sandpaper, (2) netting, and (3) dimples. The roughness coefficient ranged from k/D = 0.0028 to 0.025 (k is the roughness height). A detailed quantitative measurement of the flow field around the cylinder using Particle Imaging Velocimetry (PIV) was carried out. The hydrodynamic force coefficients (drag and lift) of the rough cylinders are compared against those of a smooth cylinder measured under the same flow conditions. It is found that certain configuration of surface roughness significantly reduces the mean drag coefficient of the cylinder, particularly at large Reynolds numbers. In addition, the root-mean-square (r.m.s.) lift coefficient of the rough cylinders is considerably lower than that of a smooth cylinder. 相似文献
6.
We here study the scouring processes that evolve around a submarine pipeline placed on a weakly cohesive seabed. We first analyze some laboratory tests carried out by Vijaya Kumar et al. [21], Xu et al. [25] and Zhou et al. [28] that focused on the scouring around a horizontal cylinder lying on a cohesive bed, subject to waves and currents. The specific purpose is that of finding a new formula for the prediction of the equilibrium scour depth under submarine pipelines. After a theoretical analysis of the main parameters, the sought formula has been found to be a function of: (i) the hydrodynamic forces acting on the cylinder (through the Keulegan–Carpenter parameter KC), (ii) the clay content of the soil Cc, and (iii) the burial depth e0/D. In the presence of small amounts of clay (Cc< 5 %), the scour depth depends directly on KC (as confirmed by many literature works for pipelines lying on sandy soils, e.g. [18]) and inversely on Cc (as already seen for bridge abutments on cohesive soils, e.g. [1]), the best-fit law being characterized by a coefficient of determination R2 = 0.62. If some burial depth is accounted for, this being a novelty of the present work, a more general formulation can be used, valid in the presence of weakly-cohesive soils and with burial depths of the pipe smaller than 0.5 (R2 = 0.79). For large clay-content ranges (2% < Cc < 75 %), the scour depth depends directly on both KC and Cc, this giving R2 = 0.79 (no burial depth) and 0.91 (some burial depth). However, this finding is at odds with the main literature, because, for large amounts of clay, it is fundamental to consider the liquidity index LI, which accounts for some important clay properties, like the plasticity. We argue that the absence of LI is balanced by the direct dependence of the scour depth on Cc. Notwithstanding the small number of available data, a formula for the prediction of the scour depth under pipelines lying on cohesive soils is fundamental for several engineering applications. The present contribution represents the first attempt to build such a formula, when the pipeline is subject to the wave-current forcing and the seabed is characterized by a relatively small clay content. 相似文献
8.
The wave transmission, reflection, and energy dissipation of the double rows of vertical piles suspending horizontal steel C shaped bars are experimentally and theoretically studied under normal regular waves. Different wave and structural parameters are investigated e.g. the wave length, the C shaped bars draft and spacing, the supporting piles diameter and spacing, and the space between the double rows. Also, the theoretical model based on an eigenfunction expansion method is developed to study the hydrodynamic breakwater performance. In order to examine the validity of the theoretical model, the theoretical results are compared with the experimental and theoretical results obtained by different authors. Comparison between experiments and predictions showed that theoretical model provides a good estimate to the different hydrodynamic coefficients when the friction factors of the upper and the lower parts are fU = 1.5 and fL = 0.75. The present breakwater physical model gives efficiency near other similar systems of different shapes. 相似文献
9.
A renewable energy harvester using the piezoelectric effect is developed for the ocean tidal and wind flow. The harvester is made of connected driving blades to an octo-generator, which has a rotator with n blades and a stator attached by eight mass-spring-piston-cylinder-piezoelectricity devices. The resonance and force magnification are utilized to increase the power output of the harvester. A corresponding mathematical model is developed to calculate the root mean square of the generated electric power. The simulation results indicate that the generated power is largely enhanced when the near-resonant condition is established. The power increases with increases in the magnetic flux density, the large-to-small diameter ratio of the cylinder, the size of magnetic bar face, and decreases in the gap between two magnetic faces and the size of the piezoelectric bar face. A generated power of 5 kW is realized by the harvester working under an ocean tidal speed, V = 1.75 m/s, and its geometric and material properties of driving length L = 7.5 m, spring constant kv = 65000 N/m, gap between the two magnets s = 0.0015 m, large to small diameter ratio of the cylinder z = 6, and magnetic flux density Br = 1.45 T. 相似文献
10.
《Coastal Engineering》2006,53(9):781-792
This paper presents results of a series of detailed measurements of geometric and migrating characteristics of ripples superimposed upon sandwaves under the action of combined waves and currents. Velocity measurements within the fluid, surface wave characteristics and 3D mapping of the bottom were recorded with an Acoustic Doppler Velocimeter (ADV), an acoustic water level sensor and a 32 composite element array of sub-aquatic acoustic sensors, respectively. Bottom records were statistically analyzed to obtain height, length and migration rates of ripples. Experiments examined ripple heights and wavelengths for the mobility factors (as defined in Eq. (4)) and the Reynolds wave number within the ranges 10 < ψ < 88 and 16 × 103 < Rew < 5 × 105, respectively. Measured values were compared with laboratory and field data together with semi-empirical and analytical formulae from the literature. Good correlation was obtained when plotting measured ripple length and length in dimensionless form as a function of the Reynolds wave number Rew. Under a given hydraulic condition, it was observed that ripples with different geometric characteristics may coexist at different locations over the sandwave. Ripple steepness is presented as a function of the Shields parameter although characterized with rather large scatter. Finally, average ripple migration speed is presented as a function of the Shields parameter and the mobility number. 相似文献
11.
《Coastal Engineering》2005,52(7):633-645
New experimental laboratory data are presented on swash overtopping and sediment overwash on a truncated beach, approximating the conditions at the crest of a beach berm or inter-tidal ridge-runnel. The experiments provide a measure of the uprush sediment transport rate in the swash zone that is unaffected by the difficulties inherent in deploying instrumentation or sediment trapping techniques at laboratory scale. Overtopping flow volumes are compared with an analytical solution for swash flows as well as a simple numerical model, both of which are restricted to individual swash events. The analytical solution underestimates the overtopping volume by an order of magnitude while the model provides good overall agreement with the data and the reason for this difference is discussed. Modelled flow velocities are input to simple sediment transport formulae appropriate to the swash zone in order to predict the overwash sediment transport rates. Calculations performed with traditional expressions for the wave friction factor tend to underestimate the measured transport. Additional sediment transport calculations using standard total load equations are used to derive an optimum constant wave friction factor of fw = 0.024. This is in good agreement with a broad range of published field and laboratory data. However, the influence of long waves and irregular wave run-up on the overtopping and overwash remains to be assessed. The good agreement between modelled and measured sediment transport rates suggests that the model provides accurate predictions of the uprush sediment transport rates in the swash zone, which has application in predicting the growth and height of beach berms. 相似文献
12.
Laminar flow past a circular cylinder with 3 small control rods is investigated by numerical simulation. This study is concerned with the suppression efficacy of vortex induced vibration by small control rods located around a main cylinder. The effects of the attack angle and rod-to-cylinder gap ratio on the hydrodynamics and vibration responses of the main cylinder are investigated. The attack angle of α = 45° is performed as the critical angle for VIV suppression of 3 control rods. The 3 control rods have no effect on VIV suppression when the attack angle is less than the critical angle. The 3 control rods have an excellent VIV suppression efficacy when the attack angle is larger than the critical angle. The transverse vibration frequency of the cylinder with 3 control rods is less than that for an isolated cylinder for all the configurations. The numerical results for the configurations of α = 45° & 60°, G/D = 0.6–1.2 show excellent suppression efficient among the cases investigated in this study. The best suppression efficient is found at α = 45°, G/D = 0.9 for 3 control rods. 2 rods in behind of the main cylinder perform more efficient than that of 1 rod in front for VIV suppression as the gap ratio of G/D less than 1.0. 相似文献
13.
《Coastal Engineering》2006,53(11):929-945
A finite difference model based on a recently derived highly-accurate Boussinesq-type formulation is presented. Up to the third-order space derivatives in terms of the velocity variables are retained, and the horizontal velocity variables are re-formulated in terms of a velocity potential. This decreases the total number of unknowns in two horizontal dimensions from seven to five, simplifying the implementation, and leading to increased computational efficiency. Analysis of the embedded properties demonstrates that the resulting model has applications with errors of 2 to 3% for (wavenumber times depth) kh ≤ 10 in terms of dispersion and kh ≤ 4 in terms of internal kinematics. The stability and accuracy of the discrete linearised systems are also analysed for both potential and velocity formulations and the advantages and disadvantages of each are discussed. The velocity potential model is then used to study physically demanding problems involving highly nonlinear wave run-up on a bottom-mounted (surface-piercing) plate. New cases involving oblique incidence are considered. In all cases, comparisons with recent physical experiments demonstrate good quantitative accuracy, even in the most demanding cases, where the local wave steepness can exceed (waveheight divided by wavelength) H / L = 0.20. The velocity potential model is additionally shown to have numerical advantages when dealing with wave–structure interactions, requiring less smoothing around exterior structural corners. 相似文献
14.
This paper provides an overview of a new large scale laboratory data set on the kinematics of breaking tsunami wavefronts. The aim of the experiments was to provide an open access data set for model testing, calibration and verification, with particular emphasis on fluid kinematics in the wave breaking and run-up (swash) zones. The experiments were performed over a composite slope in the tsunami wave basin at the O. H. Hinsdale Wave Research Laboratory at Oregon State University. Data for ten different wave conditions were collected, including non-breaking and breaking waves, and both shore breaks and fully developed long bores.Surface elevation and fluid kinematics were measured with a closely spaced array of surface piercing wave gauges, non-contact ultrasonic wave gauges and four 3-D side-looking Acoustic Doppler Velocimeters. The array was traversed from the nearshore (depth = 0.2 m) to the middle and upper run-up zone, providing kinematic data at 30 cross-shore locations. Video was also recorded from 4 cameras covering the propagation, breaking and run-up zones. Surface elevation, flow velocities and the wave maker displacement were also recorded to provide offshore boundary conditions.The experiments include conditions with wave heights up to 0.55 m, notional wave periods up to 20 s and run-up lengths of up to 15.2 m on a 1/30 slope. In terms of the slope in the shoaling and breaker zones, the data correspond to Iribarren numbers in the range of 0.26–5.6. Raw, calibrated and processed data are stored with open access within the OSU Tsunami Wave Basin Experiment Notebook, which provides full access to all the wave maker control signals, data, instrument coordinates, and processing and plotting software. This paper serves as an introduction to the data set, demonstrates data quality and provides an initial analysis of some key parameters that govern the impact of tsunami events, including run-up versus offshore wave conditions and nearshore bore height, the maximum inundation depths at the original shoreline position, and the time to maximum inundation depth and flow reversal. Examples of temporal and convective accelerations and turbulent flow components are also presented to illustrate further details of the kinematics. 相似文献
15.
《Coastal Engineering》2005,52(9):745-770
New experiments were carried out in the Large Oscillating Water Tunnel of WL|Delft Hydraulics (scale 1:1) using asymmetric 2nd-order Stokes waves. The main aim was to gain a better understanding of size-selective sediment transport processes under oscillatory plane-bed/sheet-flow conditions. The new data show that for uniform sand sizes between 0.2 < D < 1.0 mm, measured net transport rates are hardly affected by the grain size and are proportional to the third-order velocity moment. However for finer grains (D = 0.13 mm) net sand transport rates change from the ‘onshore’ direction into the ‘offshore’ direction in the high velocity range. A new measuring technique for sediment concentrations, based on the measurement of electro-resistance (see [McLean, S.R., Ribberink, J.S., Dohmen-Janssen, C.M. and Hassan, W.N.M., 2001. Sediment transport measurements within the sheet flow layer under waves and currents. J. Waterw., Port, Coast., Ocean Eng., ISSN 0733-950X]), was developed further for the improved measurement of sediment dynamics inside the sheet-flow layer. This technique enabled the measurements of particle velocities during the complete wave cycle. It is observed that for long period waves (T = 12.0 s), time-dependent concentrations inside the sheet-flow layer are nearly in phase with the time-dependent flow velocities. As the wave period decreases, the sediment entrainment from the bed as well as the deposition process back to the bed lags behind the wave motion more and more. The new data show that size-gradation has almost no effect on the net total transport rates, provided the grain sizes of the sand mixture are in the range of 0.2 < D < 1.0 mm. However, if very fine grains (D = 0.13 mm) are present in the mixture, net total transport rates of graded sand are generally reduced in comparison with uniform sand with the same D50. The transport rates of individual size fractions of a mixture are strongly influenced by the presence of other fractions in a mixture. Fine particles in sand mixtures are relatively less transported than in that uniform sand case, while the opposite occurs for coarse fractions in a mixture. The relative contribution of the coarse grains to the net total transport is therefore larger than would be expected based on their volume proportion in the original sand mixture. This partial transport behaviour is opposite to what is generally observed in uni-directional (e.g. river) flows. This is caused by vertical sorting of grain sizes in the upper bed layer and in the sheet flow and suspension layers. Kinematic sorting is believed to be responsible for the development of a coarse surface layer on top of a relatively fine sub-layer, providing in this way a relatively large flow exposure for the coarser sizes. Furthermore fine grains are suspended more easily than coarse grains to higher elevations in the flow where they are subject to increasing phase-lag effects (settling lags). The latter also leads to reduced net transport rates of these finer sizes. 相似文献
16.
Pinna nobilis is the largest endemic bivalve of the Mediterranean Sea, declared protected since 1992. Although hydrodynamic stress induced by waves is known to influence density, size and orientation of P. nobilis, the effect of other hydrological features is unknown. This paper considers a P. nobilis population living within a Posidonia oceanica meadow in the Gulf of Oristano (Sardinia, Italy). We hypothesize that spatial differences in density and orientation of P. nobilis may be related to significant wave height (HS), wave direction (DW), bottom current direction (DBC) and bottom current speed (SBC). A population of P. nobilis was investigated at different sites and its distribution was correlated to hydrodynamics by means of a numerical modeling approach. The spatial distribution was patchy, with a density of 0.06–6.7 ind. 100 m− 2. A non-uniform distribution of shell orientations (OS) was demonstrated in 4 sites out of 6. DBC and SBC were the main factors affecting OS, while waves had little influence. A SBC of 0.07 m s− 1 appears to be the threshold for inducing specimen directionality with shells aligned to the current and the ventral side exposed to the flow. This suggests that feeding strategy is a key factor in determining OS, in addition to drag minimization. We also highlighted the role of adjacent lagoons in supporting high densities as a result of high food availability. These findings demonstrate the usefulness of modeling techniques in explaining the spatial distribution pattern of P. nobilis and in contributing to our knowledge of its ecological traits. 相似文献
17.
《Deep Sea Research Part I: Oceanographic Research Papers》2006,53(3):449-473
Mode-1 internal tides were observed the western North Atlantic using an ocean acoustic tomography array deployed in 1991–1992 centered on 25°N, 66°W. The pentagonal array, 700-km across, acted as an antenna for mode-1 internal-tides. Coherent internal-tide waves with O(1 m) displacements were observed traveling in several directions. Although the internal tides of the region were relatively quiescent, they were essentially phase locked over the 200–300 day data record lengths. Both semidiurnal and diurnal internal waves were detected, with wavenumbers consistent with those calculated from hydrographic data. The M2 internal-tide energy flux was estimated to be about 70 W m−1, suggesting that mode-1 waves radiate 0.2 GW of energy, with large uncertainty, from the Caribbean island chain at this frequency. A global tidal model (TPXO 5) suggested that 1–2 GW is lost from the M2 barotropic tide over this region, but the precise value was uncertain because the complicated topography makes the calculation problematic. In any case, significant conversion of barotropic to baroclinic tidal energy does not occur in the western North Atlantic basin. It is apparent, however, that mode-1 internal tides have very weak decay and retain their coherence over great distances, so that ocean basins may be filled up with such waves. Observed diurnal amplitudes were an order of magnitude larger than expected. The amplitude and phase variations of the K1 and O1 constituents observed over the tomography array were consistent with the theoretical solutions for standing internal waves near their turning latitude. The energy densities of the resonant diurnal internal waves were roughly twice those of the barotropic tide at those frequencies. 相似文献
18.
《Ocean Modelling》2009,26(3-4):154-171
Ocean surface mixing and drift are influenced by the mixed layer depth, buoyancy fluxes and currents below the mixed layer. Drift and mixing are also functions of the surface Stokes drift Uss, volume Stokes transport TS, a wave breaking height scale Hswg, and the flux of energy from waves to ocean turbulence Φoc. Here we describe a global database of these parameters, estimated from a well-validated numerical wave model, that uses traditional forms of the wave generation and dissipation parameterizations, and covers the years 2003–2007. Compared to previous studies, the present work has the advantage of being consistent with the known physical processes that regulate the wave field and the air–sea fluxes, and also consistent with a very large number of in situ and satellite observations of wave parameters. Consequently, some of our estimates differ significantly from previous estimates. In particular, we find that the mean global integral of Φoc is 68 TW, and the yearly mean value of TS is typically 10–30% of the Ekman transport, except in well-defined regions where it can reach 60%. We also have refined our previous estimates of Uss by using a better treatment of the high frequency part of the wave spectrum. In the open ocean, Uss ≃ 0.013U10, where U10 is the wind speed at 10 m height. 相似文献
19.
When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 104 and 7.57 × 105. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem. 相似文献
20.
The ∼8.15 ka Storegga submarine slide was a large (∼3000 km3), tsunamigenic slide off the coast of Norway. The resulting tsunami had run-up heights of around 10–20 m on the Norwegian coast, over 12 m in Shetland, 3–6 m on the Scottish mainland coast and reached as far as Greenland. Accurate numerical simulations of Storegga require high spatial resolution near the coasts, particularly near tsunami run-up observations, and also in the slide region. However, as the computational domain must span the whole of the Norwegian-Greenland sea, employing uniformly high spatial resolution is computationally prohibitive. To overcome this problem, we present a multiscale numerical model of the Storegga slide-generated tsunami where spatial resolution varies from 500 m to 50 km across the entire Norwegian-Greenland sea domain to optimally resolve the slide region, important coastlines and bathymetric changes. We compare results from our multiscale model to previous results using constant-resolution models and show that accounting for changes in bathymetry since 8.15 ka, neglected in previous numerical studies of the Storegga slide-tsunami, improves the agreement between the model and inferred run-up heights in specific locations, especially in the Shetlands, where maximum run-up height increased from 8 m (modern bathymetry) to 13 m (palaeobathymetry). By tracking the Storegga tsunami as far south as the southern North sea, we also found that wave heights were high enough to inundate Doggerland, an island in the southern North Sea prior to sea level rise over the last 8 ka. 相似文献