Ultimate strength of ship hulls under torsion

For a ship hull with large deck openings such as container vessels and some large bulk carriers, the analysis of warping stresses and hatch opening deformations is an essential part of ship structural analyses. It is thus of importance to better understand the ultimate torsional strength characteristics of ships with large hatch openings. The primary aim of the present study is to investigate the ultimate strength characteristics of ship hulls with large hatch openings under torsion. Axial (warping) as well as shear stresses are normally developed for thin-walled beams with open cross sections subjected to torsion. A procedure for calculating these stresses is briefly described. As an illustrative example, the distribution and magnitude of warping and shear stresses for a typical container vessel hull cross section under unit torsion is calculated by the procedure. By theoretical and numerical analyses, it is shown that the influence of torsion induced warping stresses on the ultimate hull girder bending strength is small for ductile hull materials while torsion induced shear stresses will of course reduce the ship hull ultimate bending moment.     

Factors influencing the oxidation,reduction, methylation and demethylation of mercury species in coastal waters

The objective of this study was to examine the redox reactions and other transformations of mercury (Hg) species in surface waters, and the factors determining the rates of these reactions. For the redox studies completed at the Chesapeake Biological Laboratory (CBL), two isotopes (¹⁹⁹Hg^II and ²⁰²Hg⁰) were added into different types of filtered water (fresh to seawater) to examine the oxidation and reduction reactions. Further studies of both the redox reactions and methylation/demethylation reactions of Hg were conducted with unfiltered water on board research vessels during cruises in May and July 2005 on the Chesapeake Bay and shelf. While CH₃¹⁹⁹Hg^II was added to allow the examination of demethylation, ²⁰¹Hg^II was used to examine both reduction and methylation, and ²⁰²Hg⁰ was used to examine oxidation. Overall, the results showed that both Hg oxidation and reduction were simultaneously occurring and were photochemically mediated in the waters investigated. In contrast to the previously assumed “unreactive” nature of Hg⁰, the studies found that the magnitude of the rate constant for Hg⁰ oxidation was greater than that for reduction, indicating its importance in estuarine and coastal waters. In addition, both experiments at CBL and on board ship showed that Hg^II reduction was similar in magnitude, suggesting that biotic processes were relatively unimportant. While no measurable methylation occurred during the incubation period during the on board studies, concentration of CH₃¹⁹⁹Hg^II decreased over the time during the experiments. It appeared that the demethylation processes were not dominantly photochemically driven, but could be microbially mediated. Further studies are needed in order to help better understand Hg redox and transformations in natural water systems.     

Physiological and behavioral responses of Bathynerita naticoidea (Gastropoda: Neritidae) and Methanoaricia dendrobranchiata (Polychaeta: Orbiniidae) to hypersaline conditions at a brine pool cold seep

Bathynerita naticoidea (Gastropoda: Neritidae) and Methanoaricia dendrobranchiata (Polychaeta: Orbiniidae) are two of the most abundant invertebrates associated with cold‐seep mussel beds in the Gulf of Mexico. At the methane seep known as Brine Pool NR‐1 (27 °43.415 N, 91 °16.756 W; 650 m depth), which is surrounded by a broad band of mussels (Bathymodiolus childressi), these species have distinctly different patterns of abundance, with the gastropod being found mostly at the outer edge of the mussel bed (average density in November 2003: 817 individuals·m⁻² in outer zone, 20·m⁻² in inner zone) and the polychaete being found almost exclusively near the inner edge (average density in November 2003: 3155 individuals·m⁻² in inner zone, 0·m⁻² in outer zone), adjacent to the brine pool itself. The salinity of the brine pool exceeds 120, so we hypothesized that M. dendrobranchiata should be more tolerant of high salinities than B. naticoidea. The opposite proved to be true. The gastropods were capable of withstanding salinities at least as high as 85, whereas the polychaetes died at salinities higher than 75. Both species were osmoconformers over the range of salinities (35–75) tested. Behavioral responses of B. naticoidea to salinities of 50, 60, and 70 were investigated in inverted vertical haloclines. Gastropods generally did not enter water of salinity greater than 60, but tolerated short periods at 60. Behavioral avoidance of brine should limit the vertical distribution of B. naticoidea in the inner zone to the top 2.5–5 cm of the mussel bed. Behavior is also a likely (though unproven) mechanism for controlling horizontal distribution of this species across the mussel bed. Methanoaricia dendrobranchiata can tolerate short excursions into the brine, but probably avoids hypersaline conditions by aggregating on the tops of the mussels.     

The Vema Transverse Ridge (Central Atlantic)

Transverse ridges are elongate reliefs running parallel and adjacent to transform/fracture zones offsetting mid-ocean ridges. A major transverse ridge runs adjacent to the Vema transform (Central Atlantic), that offsets the Mid-Atlantic Ridge by 320 km. Multibeam morphobathymetric coverage of the entire Vema Transverse ridge shows it is an elongated (300 km), narrow (<30 km at the base) relief that constitutes a topographic anomaly rising up to 4 km above the predicted thermal contraction level. Morphology and lithology suggest that the Vema Transverse ridge is an uplifted sliver of oceanic lithosphere. Topographic and lithological asymmetry indicate that the transverse ridge was formed by flexure of a lithospheric sliver, uncoupled on its northern side by the transform fault. The transverse ridge can be subdivided in segments bound by topographic discontinuities that are probably fault-controlled, suggesting some differential uplift and/or tilting of the different segments. Two of the segments are capped by shallow water carbonate platforms, that formed about 3–4 m.y. ago, at which time the crust of the transverse ridge was close to sea level. Sampling by submersible and dredging indicates that a relatively undisturbed section of oceanic lithosphere is exposed on the northern slope of the transverse ridge. Preliminary studies of mantle-derived ultramafic rocks from this section suggest temporal variations in mantle composition. An inactive fracture zone scarp (Lema fracture zone) was mapped south of the Vema Transverse ridge. Based on morphology, a fossil RTI was identified about 80 km west of the presently active RTI, suggesting that a ridge jump might have occurred about 2.2 m.a. Most probable causes for the formation of the Vema Transverse ridge are vertical motions of lithospheric slivers due to small changes in the direction of spreading of the plates bordering the Vema Fracture Zone.     

A dynamic response analysis of tension leg platforms including hydrodynamic interaction in regular waves

A numerical procedure is described for predicting the motion and structural responses of tension leg platforms (TLPs) in waves. The developed numerical approach, in a TLP is assumed to be flexible instead of rigid, is based on a combination of the three dimensional source distribution method and the finite-element method. The hydrodynamic interactions among TLP members, such as columns and pontoons, are included in the motion and structural response analysis. Numerical results are compared with the experimental and numerical ones. The results of comparison confirmed the validity of the proposed approach.     

A B-spline based higher order panel method for analysis of steady flow around marine propellers

A higher order panel method based on B-spline representation for both the geometry and the solution is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then shown to be evaluated using conventional numerical quadrature. By selecting a proper order for numerical quadrature, the accuracy of the present method can be increased to the machine limit. The far- and near-field influences are shown to be evaluated based on the same far-field approximation, but the near-field solution requires subdividing the panels into smaller subpanels continuously, which can be effectively implemented due to the B-spline representation of the geometry. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution on the blade surface, including very close to the tip and trailing edge regions, with far fewer panels than existing low-order panel methods.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

3D gravity modelling for Anyongbok Seamount in the East Sea

A seamount chain with an approximately WNW trend is observed in the northeastern Ulleung Basin. It has been argued that these seamounts, including two islands called Ulleung and Dok islands, were formed by a hotspot process or by ridge related volcanism. Many geological and geophysical studies have been done for all the seamounts and islands in the chain except Anyongbok Seamount, which is close to the proposed spreading ridge. We first report morphological characteristics, sediment distribution patterns, and the crustal thickness of Anyongbok Seamount using multibeam bathymetry data, seismic reflection profiles, and 3D gravity modeling. The morphology of Anyongbok Seamount shows a cone shaped feature and is characterized by the development of many flank cones and flank rift zones. The estimated surface volume is about 60 km³, and implies that the seamount is smaller than the other seamounts in the chain. No sediments have been observed on the seamount except the lower slope, which is covered by more than 1,000 m of strata. The crustal structure obtained from a 3D gravity modeling (GFR = 3.11, SD 3.82 = mGal) suggests that the seamount was formed around the boundary of the Ulleung Plateau and the Ulleung Basin, and the estimated crustal thickness is about 20 km, which is a little thicker than other nearby seamounts distributed along the northeastern boundary of the Ulleung Basin. This significant crustal thickness also implies that Anyongbok Seamount might not be related to ridge volcanism.     

Three‐dimensional simulation of tsunami run‐up around conical island

At the circular Babi Island in the Flores tsunami (1992) and pear shaped island in the Okushiri event (1993), unexpectedly large tsunami run‐up heights in the lee of conic islands were observed. The flume and basin physical model studies were conducted in the Coastal Hydraulic Laboratory, Engineering Research and Development Center, U.S. Army Corps of Engineers to provide a better understanding of the physical phenomena and verify numerical models used in predicting tsunami wave run‐up on beaches, islands, and vertical walls. Reasonably accurate comparison of run‐up height of solitary waves on a circular island has been obtained between laboratory experimental results and two‐dimensional computation model results. In this study we apply three‐dimensional RANS model to simulate wave run‐up on conical island. In the run‐up computation we obtain that 3D calculations are in very good comparison with laboratory and 2D numerical results. A close examination of the three‐dimensional velocity distribution around conical island to compare with depth‐integrated model is performed. It is shown that the velocity distribution along the vertical coordinate is not uniform: and velocity field is weaker in the bottom layer and higher on the sea surface. The maximum difference (about 40%) appears at the time when solitary wave reached the circular island.