首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
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.  相似文献   

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

3.
Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k  ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k  ϵ model Cfk and C can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient CD is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone.  相似文献   

4.
Experimental studies were carried out to investigate the response features of an inclined flexible bare cylinder as well as a straked cylinder in a towing tank, with the main purpose of further improving the understanding of the effect of yaw angle on vortex-induced vibration (VIV) suppression. Four yaw angles (a = 0°, 15°, 30°, 45°), which is defined as the angle between the cylinder axis and the plane orthogonal to the oncoming fluid flow, were tested. The cylinder model was towed along the tank to generate a uniform fluid flow. The towing velocity was in the range of 0.05–1.0 m/s with an interval of 0.05 m/s. The corresponding Reynolds number ranged from 800 to 16000. The strakes selected for the experiments had a pitch of 17.5D and a height of 0.25D, which is generally considered as the most effective configuration for VIV suppression of a flexible cylinder in water. The experimental results indicate that VIV suppression effectiveness of the inclined flexible straked cylinder is closely related to the yaw angle. The displacement amplitudes are significantly suppressed in both cross-flow (CF) and in-line (IL) directions at a = 0°. However, with increasing yaw angle, the suppression efficiencies of the CF and IL displacement amplitudes gradually decrease. In addition, the CF dominant frequencies of the straked cylinder obviously deviate from those of the bare cylinder at a = 0° and 15°. This deviation is substantially alleviated with increasing yaw angle. The IL dominant frequencies show less dependency on the yaw angle. Similar trends are also observed on the dominant modes of vibration and the mean drag coefficients.  相似文献   

5.
6.
《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.  相似文献   

7.
An autonomous upwardly-moving microstructure profiler was used to collect measurements of the rate of dissipation of turbulent kinetic energy (ε) in the tropical Indian Ocean during a single diurnal cycle, from about 50 m depth to the sea surface. This dataset is one of only a few to resolve upper ocean ε over a diurnal cycle from below the active mixing layer up to the air–sea interface. Wind speed was weak with an average value of ~5 m s−1 and the wave field was swell-dominated. Within the wind and wave affected surface layer (WWSL), ε values were on the order of 10−7–10−6 W kg−1 at a depth of 0.75 m and when averaged, were almost a factor of two above classical law of the wall theory, possibly indicative of an additional source of energy from the wave field. Below this depth, ε values were closer to wall layer scaling, suggesting that the work of the Reynolds stress on the wind-induced vertical shear was the major source of turbulence within this layer. No evidence of persistent elevated near-surface ε characteristic of wave-breaking conditions was found. Profiles collected during night-time displayed relatively constant ε values at depths between the WWSL and the base of the mixing layer, characteristic of mixing by convective overturning. Within the remnant layer, depth-averaged values of ε started decaying exponentially with an e-folding time of 47 min, about 30 min after the reversal of the total surface net heat flux from oceanic loss to gain.  相似文献   

8.
A series of flow induced vibration (FIV) experiments for an equilateral triangle prism elastically mounted in a water channel are performed with different system stiffness at constant damping and mass. An amplitude variation coefficient is proposed to describe FIV stationarity in the present study. The FIV of the prism can be divided into three primary regions based on the amplitude and frequency responses, which are the vortex induced vibration (VIV) branch, the transition branch from VIV to galloping, and the galloping branch. The transition branch occurs at the reduced velocity in the range of 7.8 < Ur = U/(fn,air·D) < 10.4, accompanied with a relatively rapid increase in amplitude and a precipitous drop in frequency and vibration stationarity. In addition, the reduced velocity where the transition region is initiated is independent of the system stiffness. The maximum amplitude reaches 3.17 D in the galloping branch. The ratio of the response frequency to the natural frequency of the prism in air remains locked to approximately 0.65 throughout the fully developed galloping branch. Large amplitude responses in an infinite range of flow velocities, excellent vibration stationarity and steady vibration frequencies, which are characteristics of the galloping of the prism, have a positive impact on improving energy conversion.  相似文献   

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

10.
Salt intrusion in estuaries is important for ecological reasons as well as water extraction purposes. The distance salt intrudes upstream depends on a number of factors, including river discharge, tidal and wind mixing and gravitational circulation. In this paper, an analytical solution is presented for the salt intrusion in a well mixed, funnel-shaped estuary whose cross sectional area decreases exponentially (with decay coefficient β) with distance, x, inland, and in which longitudinal mixing is constant along the length of the estuary. The solution predicts that a graph of the logarithm of salinity against exp (βx) should be a straight line, with slope proportional to the mixing coefficient Kx. The solution is tested against observations from 15 surveys over a four-year period in the Incomati estuary. Good straight line fits, as predicted, are observed on all surveys, with a mean R2 = 0.97. The average value of Kx for all surveys is 38 m2 s−1. The solution is used to make predictions about the minimum river flow required to prevent salt intruding to an extent where it causes a detrimental effect on water extraction. The minimum recommended river flow required to prevent this is 35 m3 s−1. In recent years, flow has fallen below this level for several months each year.  相似文献   

11.
Ocean surface gravity waves play a major role in many engineering and environmental problems, both in the open ocean and in coastal zones. Therefore, it is essential to improve our knowledge on spatial and temporal variability of wave climate. This study aims at investigating this variability in the North-East Atlantic Ocean (25°W–0°W and 30°N–60° N), using a 57-year hindcast (1953–2009) obtained with a spectral wave model forced with reanalysis wind fields. The hindcast analysis reveals firstly strong seasonal fluctuations of wave climate, with winters characterized by large and long-period waves of mean direction spreading from south-west to north-west, and summers characterized by smaller and shorter-period waves originating from norther directions. From northern (55°N) to southern (35°N) latitudes, the significant wave height (Hs) decreases by roughly 40%, the mean wave direction (Mwd) rotates clockwise by about 25% while the peak period (Tp) only grows by 5%. These three parameters also exhibit a strong inter-annual variability, particularly when winter-means (from 1st of December to 1st of April) are considered. Linear trend analysis over the studied period shows spatially variable long-term trends, with a significant increase of Hs (up to 0.02 m yr?1) and a counterclockwise shift of Mwd (up to ?0.1° yr?1) at northern latitude, contrasting with a fairly constant trend for Hs and a clockwise shift of Mwd (up to +0.15° yr?1) at southern latitudes. Long-term trends of Tp are less significant, with still a slight increase in the north-eastern part of the study area (up to +0.01 s yr?1). Eventually, a comparison between the inter-annual variability of the winter-means of the three selected wave parameters and the North Atlantic Oscillation (NAO) reveals: (1) a strong positive correlation between Hs and the NAO index at northern latitudes (correlation coefficient up to R = 0.91) and a significant negative correlation at southern latitudes (up to R = ?0.6); (2) no significant correlation for Mwd north of 40°N and a clear positive correlation southward of 40°N (up to R = 0.8) and (3) a northward increasing correlation for Tp (up to R = 0.8). Long-term trends for Hs, Mwd and Tp are finally explained by a significant increase in the NAO index over the studied period.  相似文献   

12.
We describe experiments with multi-directional focused waves interacted with a vertical circular cylinder in a 3D wave basin. The focus of this study is on the run-up of multi-directional focused waves, wave forces, and wave pressures on the cylinder. Part I, the study on wave run-up, has already been presented by Li et al. (2012). In this paper, the analysis of the wave force on the vertical cylinder is presented.In this experiment, a cylinder with 0.25 m in diameter was adopted and different wave parameters, such as focused wave amplitude, peak frequency, frequency bandwidth and directional spreading index, are considered. The model scale kpa (kp is the wave number corresponding to peak frequency, a is the radium of the cylinder) varies from 0.32 to 0.65. The maximum forces of multi-directional focused wave on cylinder were measured and investigated. The results showed that the wave parameters have a significant influence on the wave force, and that the spatial profile of the surface of multi-directional focused wave can also affect its force on the cylinder, which is different from two-dimensional wave. In addition, the ‘secondary loading cycle’ phenomenon was also observed and discussed. In our experiments, the ‘secondary loading cycles’ occur when kA > 0.36 for all cases. While in some referred small scale experiments, the secondary load cycles are observed even for kA = 0.2, when the waves are longer enough. To larger model scale, the pronounced secondary load cycle occurs with larger wave steepness waves.  相似文献   

13.
14.
The results from a~1 km resolution HYbrid Coordinate Ocean Model (HYCOM), forced by 1/2° Navy Operational Global Atmospheric Prediction System (NOGAPS) atmospheric data, were used in order to study the dynamic response of the Persian Gulf to wintertime shamal forcing. Shamal winds are strong northwesterly winds that occur in the Persian Gulf area behind southeast moving cold fronts. The period from 20 November to 5 December 2004 included a well defined shamal event that lasted 4–5 days. In addition to strong winds (16 m s?1) the winter shamal also brought cold dry air (Ta=20 °C, qa=10 g kg?1) which led to a net heat loss in excess of 1000 W m?2 by increasing the latent heat flux. This resulted in SST cooling of up to 10 °C most notably in the northern and shallower shelf regions. A sensitivity experiment with a constant specific humidity of qa=15 g kg?1 confirmed that about 38% of net heat loss was due to the air–sea humidity differences. The time integral of SST cooling closely followed the air–sea heat loss, indicating an approximate one-dimensional vertical heat balance. It was found that the shamal induced convective vertical mixing provided a direct mechanism for the erosion of stratification and deepening of the mixed layer by 30 m. The strong wind not only strengthened the circulation in the entire Persian Gulf but also established a northwestward flowing Iranian Coastal Current (ICC, 25–30 cm s?1) from the Strait of Hormuz to about 52°E, where it veered offshore. The strongest negative sea level of 25–40 cm was generated in the northernmost portion of the Gulf while the wind setup against the coast of the United Arab Emirates established a positive sea level of 15–30 cm. The transport through the Strait of Hormuz at 56.2°E indicated an enhanced outflow of 0.25 Sv (Sv≡106 m3 s?1) during 24 November followed by an equivalent inflow on the next day.  相似文献   

15.
《Marine Chemistry》2007,103(1-2):30-45
The chemistry of dissolved Fe(III) was studied in the Scheldt estuary (The Netherlands). Two discrete size fractions of the dissolved bulk (< 0.2 μm and < 1 kDa) were considered at three salinities (S = 26, 10 and 0.3).Within the upper estuary, where fresh river water meets seawater, the dissolved Fe concentration decreases steeply with increasing salinity, for the fraction < 0.2 μm from 536 nM at S = 0.3 to 104 nM at S = 10 and for the fraction < 1 kDa from 102 nM to 36 nM Fe. Further downstream, in the middle and lower estuary, this decrease in the Fe concentration continues, but is far less pronounced. For all samples, the traditionally recognised dissolved strong organic Fe-binding ligand concentrations are lower than the dissolved Fe concentrations.Characteristics of dissolved Fe-binding ligands were determined by observing kinetic interactions with adsorptive cathodic stripping voltammetry. From these kinetic experiments we concluded that apart from the well-known strong Fe-binding organic ligands (L, logK = 19–22) also weak Fe-binding ligands (P) existed with an α value (binding potential = K · [P]) varying between 1011.1 and 1011.9. The presence of this relatively weak ligand explained the high concentrations of labile Fe present in both size fractions in the estuary. This weak ligand can retard or prevent a direct precipitation after an extra input of Fe.The dissociation rate constants of the weak ligand varied between 0.5 × 10 4 and 4.3 × 10 4 s 1. The rate constants of the strong organic ligand varied between kd = 1.5 × 10 3–17 × 10 2 s 1 and kf = 2.2 × 108–2.7 × 109 M 1 s 1. The dissociation rate constant of freshly amorphous Fe-hydroxide was found to be between 4.3 × 10 4 and 3.7 × 10 3 s 1, more labile or equal to the values found by Rose and Waite [Rose, A.L., Waite, T.D., 2003a. Kinetics of hydrolysis and precipitation of ferric iron in seawater. Environ. Sci. Technol., 37, 3897–3903.] for freshly precipitated Fe in seawater.Kinetic rate constants of Fe with the ligand TAC (2-(2-Thiazolylazo)-p-cresol) were also determined. The formation rate constant of Fe(TAC)2 varied between 0.1 × 108 and 3.6 × 108 M 1 s 1, the dissociation rate constant between 0.2 × 10 5 and 17 × 10 5 s 1 for both S = 26 and S = 10. The conditional stability constant of Fe(TAC)2 (βFe(TAC)2′) varied between 22 and 23.4 for S = 10 and S = 26 more or less equal to that known from the literature (logβFe(TAC)2 = 22.4; [Croot, P.L., Johansson, M., 2000. Determination of iron speciation by cathodic stripping voltammetry in seawater using the competing ligand 2-(2-Thiazolylazo)-p-cresol (TAC). Electroanalysis, 12, 565–576.]). However, at S = 0.3 the logβFe(TAC)2′ was 25.3, three orders of magnitude higher. Apparently the application of TAC to samples of low salinity can only be done when the correct βFe(TAC)2′ is known.  相似文献   

16.
Recently obtained World Ocean Circulation Experiment (WOCE) sections combined with a specially prepared pre-WOCE South Atlantic data set are used to study the dianeutral (across neutral surface) mixing and transport achieving Antarctic Intermediate Water (AAIW) being transformed to be part of the North Atlantic Deep Water (NADW) return cell. Five neutral surfaces are mapped, encompassing the AAIW from 700 to 1100 db at the subtropical latitudes.Coherent and significant dianeutral upwelling is found in the western boundary near the Brazil coast north of the separation point (about 25°S) between the anticyclonic subtropical and cyclonic south equatorial gyres. The magnitude of dianeutral upwelling transport is 10-3 Sv (1 Sv=106 m3 s-1) for 1°×1° square area. It is found that the AAIW sources from the southwestern South Atlantic and southwestern Indian Ocean do not rise significantly into the Benguela Current. Instead, they contribute to the NADW return formation by dianeutral upwelling into the South Equatorial Current. In other words, the AAIW sources cannot obtain enough heat/buoyancy to rise until they return to the western boundary region but north of the separation point. The basin-wide integration of dianeutral transport shows net upward transports, ranging from 0.25 to 0.6 Sv, across the lower and upper boundary of AAIW north of 40°S. This suggests that the equatorward AAIW is a slow rising water on a basin average. Given one order of uncertainty in evaluating the along-neutral-surface and dianeutral diffusivities from the assumed values, K=103 m2 s-1 and D=10-5 m2 s-1, the integrated dianeutral transport has an error band of about 10–20%. The relatively weak integrated dianeutral upwelling transport compared with AAIW in other oceans implies much stronger lateral advection of AAIW in the South Atlantic.Mapped Turner Angle in diagnosing the double-diffusion processes shows that the salty Central Water can flux salt down to the upper half of AAIW layer through salt-fingering. Therefore, the northward transition of AAIW can gain salt either through along-neutral-surface advection and diffusion or through salt fingering from the Central Water and heat through either along-neutral-surface advection and diffusion or dianeutral upwelling. Cabbeling and thermobaricity are found significant in the Antarctic frontal zone and contribute to dianeutral downwelling with velocity as high as −1.5×10-7 m s-1. A schematic AAIW circulation in the South Atlantic suggests that dianeutral mixing plays an essential role in transforming AAIW into NADW return formation.  相似文献   

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

18.
Yearlong 75 kHz acoustic Doppler current profiler (ADCP) data were obtained well above Reykjanes Ridge (northern extension of the Mid-Atlantic Ridge (MAR)). The area is characterized by relatively large semidiurnal tidal (‘D2’) currents that have (at lunar M2) more than half a decade larger variance than inertial (f) currents. However, despite the relatively weak near-inertial kinetic energy, its vertical current shear shows larger magnitudes than at M2 in an otherwise flat fD2 band limited between frequencies 0.74 and 1.35f, which equals the inertio-gravity wave bounds [σmin, σmax](N=f). N represents the buoyancy frequency. The shear in this band dominates all shear computed at 20 m effective vertical scale. As the kinetic energy spectrum peaks at M2, but not (significantly) at S2 and N2, a difference in tidal (and inertial) scales and hence sources is observed. M2-tides contribute mostly to large-scale coherent motions. The dominant incoherent fD2 shear is highly variable in time (∼2-day periodicity). Furthermore, inertial and tidal shear are more or less completely separated in space and time, each occurring in different layers in the vertical. The thin shear layers reflect the rapidly varying short vertical scale Ns profile, to within the ∼20 m limitation of ADCP data, rather than the large-scale smooth NL. In each of large-Ns layers Ri≈1, probably. The yearlong smoothed shear magnitude follows NL, but only as stable Ri≈5. The shear polarization is more circular than rectilinear, albeit varying with time, and highly symmetric around f. During transitions, e.g., between stratified and homogeneous layers and between waves from varying sources, near-circular motions can generate near-rectilinear shear in the direction of wave propagation (in the direction of the minor axis of the current ellipse). This contrasts with the possibility of near-rectilinear barotropic oscillatory motions generating near-circular shear under viscosity in shallow seas.  相似文献   

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

20.
Wave–current flow is a phenomenon that is present in many practical engineering situations. Over the past several decades, this type of flow has been increasingly investigated under controlled laboratory conditions. This paper presents a numerical study of wave–current flow in the ocean basin of the LabOceano (COPPE/UFRJ). A homogeneous multiphase model based on the RANS equations and the kɛ turbulence model implemented in ANSYS-CFX code were used. A cross section of the ocean basin was represented. A regular wave with a height of 0.08 m and a period of 1.80 s (i.e., a wave steepness of H/L = 0.016), propagating on favourable currents, was simulated. The behaviour of the free surface elevation over time and the streamlines along the basin for wave and wave–current flows were presented. The numerical results were compared to the non-viscous theory given by the Rayleigh equation applied to the problem of wave–current interaction. Good agreement was found between the wave length estimated by the numerical results and the analytical solutions, with a deviation of less than 2%.  相似文献   

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

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