On the numerical accuracy of the prediction of resistance coefficients in ship stern flow calculations

This article presents a study on the accuracy of the numerical determination of the friction and pressure resistance coefficients of ship hulls. The investigation was carried out for the KVLCC2 tanker at model- and full-scale Reynolds numbers. Gravity waves were neglected, i.e., we adopted the so-called double-model flow. Single-block grids with H–O topology were adopted for all the calculations. Three eddy viscosity models were employed: the one-equation eddy viscosity and the two-equation models proposed by Menter and the TNT version of the two-equation k-ω model. Verification exercises were performed in sets of nearly geometrically similar grids with different densities in the streamwise, normal, and girthwise directions. The friction and pressure resistance coefficients were calculated for different levels of the iterative error and for computational domains of different size. The results show that on the level of grid refinement used, it is possible to calculate the viscous resistance coefficients in H–O grids that do not match the ship contour with a numerical uncertainty of less than 1%. The differences between the predictions of different turbulence models were larger than the numerical uncertainty; however, these differences tended to decrease with increases in the Reynolds number. The pressure resistance was remarkably sensitive to domain size and far-field boundary conditions. Either a large domain or the application of a viscous–inviscid interaction procedure is needed for reliable results. This work was presented in part at the International Conference on Computational Methods in Marine Engineering—MARINE 2007, Barcelona, June 3–4, 2007.     

Sidescan sonar survey of a dumping site in the Mecklenburg Bight (south-western Baltic Sea)

Two dumping test sites of dredged sediment (glacial till, mixed sediment with sand) in the south-western Baltic Sea were repeatedly investigated with sidescan sonar. The first survey was conducted before dumping, the second survey 1 week after dumping, and eight more surveys were run during the following three and a half years. Sidescan mosaics were calculated from raw data. Comparing the mosaics, it becomes obvious that the initial strong microrelief of the dumping sites vanishes with time. The heaps of dumped material were eroded. Coarse material remains at the surface, fine material fills in the gaps between the heaps. Fine sediment structures (filaments and aureoles) created by the dumping process, and elongated traces of dumped material outside the dumping places disappeared with time.     

The use of polycyclic aromatic hydrocarbons as a particulate tracer in the water column of Gaoping (Kaoping) Submarine Canyon

Time-series samples of settling particles were collected in the water column of Gaoping (formerly spelled Kaoping) Submarine Canyon (KPSC) with two sediment traps on taut-line moorings deployed at two different depths (60 and 280 m) between May 26 and June 27, 2004. Average total polycyclic aromatic hydrocarbon (PAH) concentrations of upper and lower trap array samples were 310 ± 61 ng g^− 1 dw (range: 200–440) and 240 ± 36 ng g^− 1 dw (range: 180–290), respectively. Principal component analysis results suggest that PAH sources in the trap-collected particles included diesel vehicle/coal burning, diagenetic sources, and petroleum release. PAH downward fluxes based on settling particles were estimated to be 12–44 μg m^− 2 d^− 1. These values are higher than those reported in the literature for most coastal areas. During the sampling period, both traps were significantly tilted by tidal current and fluctuated vertically. The upper traps experienced greater vertical movements, thus their particle characteristics (e.g., POC, particle mass, and fine particle fraction) varied more than those of the lower traps. Hourly depth variations of the tilted sediment trap array were echoed by the corresponding total PAH concentrations. Moreover, the PAH composition of the collected particles was related to the flow direction and speed. These observations suggest that PAHs can be used as an effective chemical tracer for the transport of terrestrial and marine particulates in a complex aquatic environment like Gaoping (Kaoping) Submarine Canyon.     

Water entry of a bow-flare ship section with roll angle

The two-dimensional water entry of a bow-flare ship section with constant roll angle, or heel angle, was studied by using a boundary element method. The fully nonlinear free surface conditions and exact body boundary conditions were satisfied. Nonviscous flow separation from the knuckles of the section or from the curved bottom could be simulated. The numerical calculations were compared with existing experimental results. First, the effects of roll angle were investigated and then the characteristics associated with large roll angles were examined in particular. The evolution of the free surfaces and the pressure distributions on the section surface are illustrated and the influence of nonviscous flow separation from the leeward section surface is discussed.     

Characteristic flow patterns generated by macrozoobenthic structures

A laboratory flume channel, equipped with an acoustic Doppler flow sensor and a bottom scanning laser, was used for detailed, non-intrusive flow measurements (at 2 cm s^− 1 and 10 cm s^− 1) around solitary biogenic structures, combined with high-resolution mapping of the structure shape and position. The structures were replicates of typical macrozoobenthic species commonly found in the Mecklenburg Bight and with a presumed influence on both, the near-bed current regime and sediment transport dynamics: a worm tube, a snail shell, a mussel, a sand mound, a pit, and a cross-stream track furrow. The flow was considerably altered locally by the different protruding structures (worm tube, snail, mussel and mound). They reduced the horizontal approach velocity by 72% to 79% in the wake zone at about 1–2 cm height, and the flow was deflected around the structures with vertical and lateral velocities of up to 10% and 20% of the free-stream velocity respectively in a region adjacent to the structures. The resulting flow separation (at flow Reynolds number of about 4000 and 20,000 respectively) divided an outer deflection region from an inner region with characteristic vortices and the wake region. All protruding structures showed this general pattern, but also produced individual characteristics. Conversely, the depressions (track and pit) only had a weak influence on the local boundary layer flow, combined with a considerable flow reduction within their cavities (between 29% and 53% of the free-stream velocity). A longitudinal vortex formed, below which a stagnant space was found. The average height affected by the structure-related mass flow rate deficit for the two velocities was 1.6 cm and 1.3 cm respectively (80% of height and 64%) for the protruding structures and 0.6 cm and 0.9 cm (90% and 127% of depth) for the depressions. Marine benthic soft-bottom macrozoobenthos species are expected to benefit from the flow modifications they induce, particularly in terms of food particle capture due to altered particle pathways and residence times, but also for the exchange of gases, solutes and spawn. The present results confirm previous studies on flow interaction effects of various biogenic structures, and they add a deeper level of detail for a better understanding of the fine-scale effects.     

On the observability of the fortnightly cycle of the Tagus estuary turbid plume using MODIS ocean colour images

Using three years (2003 to 2005) of MODIS-Aqua normalized water-leaving radiance at 551 nm this paper shows a fortnightly cycle in the Tagus estuary turbid plume. The Tagus estuary is one of the largest estuaries of the west coast of Europe and is located in the most populated area of Portugal, including the capital Lisbon. The turbid plume has been detected by the backscattering characteristics of the surface waters in the vicinity of the estuary mouth. In fortnightly scales, the turbid plume has smaller dimensions during and after neap tides and higher dimensions during and after spring tides. This is most probably associated with the fortnightly spring–neap tidal cycle and the consequent increase in turbidity inside the estuary during spring tides. During the summer weak spring tides (tidal amplitude approximately 2.5 m) no turbid plume is observed for an entire fortnightly cycle. Outside the summer months, precipitation, river discharge and winds, were found to increase the turbid area, but the fortnightly cycle appears to be superimposed on the large time-scale variability, and present throughout the year.     

The structure of dissipation in the western Irish Sea front

We report on an intensive campaign in the summer of 2006 to observe turbulent energy dissipation in the vicinity of a tidal mixing front which separates well mixed and seasonally stratified regimes in the western Irish Sea. The rate of turbulent dissipation ε was observed on a section across the front by a combination of vertical profiles with the FLY dissipation profiler and horizontal profiles by shear sensors mounted on an AUV (Autosub). Mean flow conditions and stratification were obtained from a bed mounted ADCP and a vertical chain of thermistors on a mooring. During an Autosub mission of 60 h, the vehicle, moving at a speed of ~ 1.2 m s^− 1, completed 10 useable frontal crossings between end points which were allowed to move with the mean flow. The results were combined with parallel measurements of the vertical profile of ε which were made using FLY for periods of up to 13 h at positions along the Autosub track. The two data sets, which show a satisfactory degree of consistency, were combined to elucidate the space–time variation of dissipation in the frontal zone. Using harmonic analysis, the spatial structure of dissipation was separated from the strong time dependent signal at the M₄ tidal frequency to yield a picture of the cross-frontal distribution of energy dissipation. A complementary picture of the frontal velocity field was obtained from a moored ADCP and estimates of the mean velocity derived from the thermal wind using the observed density distribution. which indicated the presence of a strong (0.2 m s^− 1) jet-like flow in the high gradient region of the front. Under neap tidal conditions, mean dissipation varied across the section by 3 orders of magnitude exceeding 10^− 2 W m^− 3 near the seabed in the mixed regime and decreasing to 10^− 5 W m^− 3. in the strongly stratified interior regime. The spatial pattern of dissipation is consistent in general form with the predictions of models of tidal mixing and does not reflect any strong influence by the frontal jet.     

Quasi-two-dimensional flow on the polar β-plane: Laboratory experiments

Geophysical turbulence is strongly affected by the variation of the Coriolis parameter with latitude. This variation results in the so-called β-effect, which forces energy from small-scales to be transferred preferentially into zonal motions. This effect results in the formation of narrow jet-like zonal flows that dominate the dynamics and act as transport barriers. Here, laboratory experiments are used to reproduce this effect in decaying turbulent flows. An electromagnetic cell is used to generate an initial field of vorticity in a rotating tank. Under conditions of quasi-geostrophic flow, the β-effect is produced by depth variation of the flow instead of variation of the Coriolis parameter. The effects of changing the container geometry and the overall fluid depth on the production of jets are investigated. The results suggest that this laboratory configuration can be used to model jet formation in the oceans and that increasing fluid depth is a practical way to decrease viscous effects.     

Non-local simulation of the stable boundary layer with a third order moments closure model

A new higher order closure model for the stable boundary layer is presented and compared with Large Eddy Simulation data. The model includes numerical solutions for the mean values, second and third order moments equations. A satisfactory agreement is found between the calculated vertical profiles of the turbulent quantities with those provided by the LES. Furthermore the new model results are compared with profiles obtained with a lower order closure model in order to verify the effective importance of including third order dynamical equations in the model.     

Analysis of velocity fluctuations and their intermittency properties in the surf zone using empirical mode decomposition

We consider here surf zone turbulence measurements, recorded in the Eastern English Channel using a sonic anemometer. In order to characterize the intermittent properties of their fluctuations at many time scales, we analyze the experimental time series using the Empirical Mode Decomposition (EMD) method. The series is decomposed into a sum of modes, each one narrow-banded, and we show that some modes are associated with the energy containing wave breaking scales, and other modes are associated with small-scale intermittent fluctuations. We use the EMD approach in association with a newly developed method based on Hilbert spectral analysis, representing the probability density function in an amplitude–frequency space. We then characterize the fluctuations in a stochastic framework using a cumulant generating function for all scales, and compare the results obtained from direct and classical structure function analysis, to EMD–Hilbert spectral analysis results, showing that the former method saturates at large scales, whereas the latter method is more precise in its scale approach. These results show the strength of the new EMD–Hilbert spectral analysis method for data presenting a strong forcing such as found in shallow water, wave dominated situations.