Large-amplitude motion responses of a Ro–Ro ship to regular oblique waves in intact and damaged conditions

This article presents a nonlinear time-domain simulation method for the prediction of large-amplitude motions of a Ro–Ro ship in regular oblique waves in an intact and a damaged condition. Numerical computations and model tests have been carried out to investigate the dynamic motion responses of Ro–Ro ship Dextra to various wave amplitudes at three different wave headings. The results of numerical and experimental investigations for stern quartering waves are reviewed. Comparisons between predictions and measurements show good agreement except in the roll-resonant region. Nonlinear effects are significant in horizontal modes of motion, and resonant roll motion, and there is strong coupling between all modes of motion in the roll-resonant region for large-amplitude responses. On the other hand, the time-domain simulation technique suffers from numerical drift in horizontal modes of motion as wave amplitude increases. This is due to nonlinear equations of motion and the lack of a restoring force and moment in horizontal motion. Received: April 30, 2002 / Accepted: August 9, 2002 Acknowledgments. II Programme of the European Community Commission under contract No. BRPR-CT97-0513. Address correspondence to: H.S. Chan (hoi-sang.chan@ncl.ac.uk)     

Distributed mass/discrete floe model for pack ice rheology computation

A new model, called the distributed mass/discrete floe model, is proposed for performing practical computations of mesoscale pack ice rheology. This model possesses the advantages of both the continuum and the discrete element models: it can express the discrete nature of pack ice, for which it is difficult to use a continuum model, and requires a much shorter computation time than a discrete element model. The pack ice is divided into ice bunches in which the floes, assumed to be distributed uniformly, are modeled as inelastic disks or rectangles floating on the water. The ice interaction forces are formulated from the relationship between the impulse on the bunch and the variation of momentum in the bunch. The ocean flow is calculated simultaneously with the floe movement using a multilayer model. In a circulating water channel, drift tests of physical model floes were performed in order to investigate the characteristics of their motion and interaction with ocean structure models. Near the structure, the floe motion depends on the floe shape. Disk floes show a lateral motion in front of the structure. They flow out around both sides of the structure and the number of floes in front of the structure decreases with the lapse of time. On the other hand, rectangular floes scarcely flow laterally. The number of floes in front of the structure remains constant over time. These experiments indicate that when the motion of pack ice around a structure is simulated, it is important to consider the floe shape. The disk floe motion and the rectangular floe motion can be regarded as extreme cases of pack ice motion. Actual pack ice motion may be between these two extremes. Computations were carried out using the distributed mass/discrete floe (DMDF) model. Simulation results were compared with the circulating water channel experiment results and sea ice motion in the southern part of the Sea of Okhotsk. The DMDF model predicted the circulating water channel drift test results quite closely. The DMDF model results also compared quite well with the sea ice motion.     

PMM tests in a circulating water channel and nonlinear analysis

As a result of frequent marine disasters leading to the loss of human life and pollution of vast areas of the ocean, ship manoeuvrability has become a very important characteristic of ship design. Among several recent experimental techniques to determine ship manoeuvrability, the most popular is captive model testing using a planar motion mechanism (PMM). This article describes some tests, analyses, and results of PMM tests in a circulating water channel (CWC) using a model of a training ship. The hydrodynamic forces and moments acting on a model of the training ship Shioji Maru in pure yawing motion were measured, and hydrodynamic derivatives were obtained using two different methods of analysis: singular value decomposition (a least-squares fit method) and Fourier analysis. Derivatives obtained from the tests were used to simulate the turning trajectory of the actual ship, and these were compared with the results of sea trials. The results indicate that both methods of analysis yield fairly similar derivatives. The simulation results were also found to be a close match with the trial results. Received: February 7, 2002 / Accepted: May 14, 2002 Address correspondence to: K. Shoji (shoji@ipc.tosho-u.ac.jp)     

基于灰色系统理论的滑行艇运动姿态预报(英文)

In order to minimize the harm caused by the instability of a planing craft, a motion prediction model is essential. This paper analyzed the feasibility of using an MGM(1,N) model in grey system theory to predict planing craft motion and carried out the numerical simulation experiment. According to the characteristics of planing craft motion, a recurrence formula was proposed of the parameter matrix of an MGM(1,N) model. Using this formula, data can be updated in real-time without increasing computational complexity significantly. The results of numerical simulation show that using an MGM(1,N) model to predict planing motion is feasible and useful for prediction. So the method proposed in this study can reflect the planing craft motion mechanism successfully, and has rational and effective functions of forecasting and analyzing trends.     

Broaching prediction in the light of an enhanced mathematical model,with higher-order terms taken into account

 We have attempted to develop a more consistent mathematical model for capsizing associated with surf-riding in following and quartering waves by taking most of the second-order terms of the waves into account. The wave effects on the hull maneuvring coefficients were estimated, together with the hydrodynamic lift due to wave fluid velocity, and the change in added mass due to relative wave elevations. The wave effects on the hydrodynamic derivatives with respect to rudder angles were estimated by using the Mathematical Modelling Group (MMG) model. Then captive ship model experiments were conducted, and these showed reasonably good agreements between the experiments and the calculations for the wave effects on the hull and the rudder maneuvring forces. It was also found that the wave effects on restoring moments are much smaller than the Froude–Krylov prediction, and the minimum restoring arm appears on a wave downslope but not on a wave crest amidship. Thus, an experimental formula of the lift force due to the heel angle of the ship is provided for numerical modelling. Numerical simulations were then carried out with these second-order terms of waves, and the results were compared with the results of free-running model experiments. An improved prediction accuracy for ship motions in following and quartering seas was demonstrated. Although the boundaries of the ship motion modes were also obtained with both the original model and the present one, the second-order terms for waves are not so crucial for predicting the capsizing boundaries themselves. Received: June 20, 2002 / Accepted: October 10, 2002 Acknowledgments. This research was supported by a Grant-in-Aid for Scientific Research of the Ministry of Education, Culture, Sports, Science and Technology of Japan (No. 13555270). The authors thank Prof. N. Rakhmanin of the Krylov Ship Research Institute for providing the Russian literature, as well as Mr. H. Murata of NHK (Japan Broadcasting Corporation) for translating it into Japanese. Address correspondence to: N. Umeda (e-mail: umeda@naoe.eng.osaka-u.ac.jp)     

Winter–spring dynamics in sea-ice carbon cycling in the coastal Arctic Ocean

An understanding of microbial interactions in first-year sea ice on Arctic shelves is essential for identifying potential responses of the Arctic Ocean carbon cycle to changing sea-ice conditions. This study assessed dissolved and particulate organic carbon (DOC, POC), exopolymeric substances (EPS), chlorophyll a, bacteria and protists, in a seasonal (24 February to 20 June 2004) investigation of first-year sea ice and associated surface waters on the Mackenzie Shelf. The dynamics of and relationships between different sea-ice carbon pools were investigated for the periods prior to, during and following the sea-ice-algal bloom, under high and low snow cover. A predominantly heterotrophic sea-ice community was observed prior to the ice-algal bloom under high snow cover only. However, the heterotrophic community persisted throughout the study with bacteria accounting for, on average, 44% of the non-diatom particulate carbon biomass overall the study period. There was an extensive accumulation of sea-ice organic carbon following the onset of the ice-algal bloom, with diatoms driving seasonal and spatial trends in particulate sea-ice biomass. DOC and EPS were also significant sea-ice carbon contributors such that sea-ice DOC concentrations were higher than, or equivalent to, sea-ice-algal carbon concentrations prior to and following the algal bloom, respectively. Sea-ice-algal carbon, DOC and EPS-carbon concentrations were significantly interrelated under high and low snow cover during the algal bloom (r values ≥ 0.74, p < 0.01). These relationships suggest that algae are primarily responsible for the large pools of DOC and EPS-carbon and that similar stressors and/or processes could be involved in regulating their release. This study demonstrates that DOC can play a major role in organic carbon cycling on Arctic shelves.     

Long-term prediction method for the green water load and volume for an assessment of the load line

 To develop a practical prediction method for the green water load and volume on the bow deck in irregular waves, model tests were conducted using a tanker and a cargo ship on a domestic Japanese voyage. The relation between green water load and relative water height at the stem was considered. Based on the finding that the maximum value of the green water load is proportional to the square of the maximum value of the water elevation over the bow top, the probability density functions of the green water load and volume in short-term predictions were proposed. It was verified that the proposed functions show good agreement with the measured distributions, especially in the tails, and were better than conventional functions. Using these functions, long-term predictions of the green water load were carried out. It was confirmed that the present method is more rational than the conventional one for estimating the long-term probability of the green water load. An assessment of the bow height of a domestic Japanese ship from the viewpoint of deck wetness was carried out using these prediction methods. This research was used as the technical background for the revision of domestic rules on load lines, which was enforced in October 2001. Received: July 19, 2002 / Accepted: October 30, 2002 Acknowledgment. Some of the present study was carried out as part of a cooperation project (RR45) with the Shipbuilding Research Association of Japan, supported by the Nippon Foundation. Address correspondence to: Y. Ogawa (e-mail: ogawa@nmri.go.jp)     

Study on real-time estimation of the ship motion cross spectra

 Time-varying coefficient vector autoregressive (T-VVAR) modeling with instantaneous responses is applied to spectrum analysis based on the nonstationary motion data of ships. Because of the ship's maneuvers, changes such as course and speed, the ship motions in waves are regarded as a nonstationary random process, although the seaway can be considered as a stationary stochastic process. The T-VVAR model is transformed into a state space model, and the time-varying coefficients can be evaluated by using the Kalman filter algorithm. Using the estimated time-varying coefficients, the instantaneous cross spectra of the ship motions can be calculated at every moment. In order to examine the reliability of the proposed procedure, on-board tests were carried out. Under stationary conditions, at a constant speed and course, the proposed method shows good agreement with stationary vector autoregressive (SVAR) modeling analysis. Moreover, it is confirmed that the proposed method can estimate the instantaneous cross spectra of the ship motions even under nonstationary conditions, showing that this is a powerful tool for on-line analysis of the nonstationary motion data of ships. Received: August 2, 2002 / Accepted: November 28, 2002 Acknowledgments. The authors thank the captain and crew of the training ship Shioji Maru, Tokyo University of Mercantile Marine. Address correspondence to: T. Iseki (iseki@ipc.tosho-u.ac.jp) Updated from the Japanese original, which won the 2002 SNAJ prize (J Soc Nav Archit Jpn 2001;190:161–168)     

Elastic response characteristics of a very large floating structure in waves moored inside a reef

 This article describes the results of hydraulic model tests of the elastic response of a very large floating structure (VLFS) moored inside a reef in an isolated island. The distributions of strains and vertical displacements due to the elastic response of the VLFS were measured. The response characteristics were strongly affected by deformed nonlinear waves inside the reef. A two-step analytical method to compute the elastic response of a VLFS is proposed, and its validity is verified using the results of the hydraulic model tests. Received: May 2, 2002 / Accepted: March 17, 2003 RID="*" ID="*" Address correspondence to: S. Shiraishi (shiraishi@pari.go.jp) Acknowledgment. This study was supported by the Program for Promotions of Fundamental Transport Technology Research from the Corporation for Advanced Transport and Technology (CATT).     

Mixed analysis of shell and solid elements using the overlaying mesh method

 This article presents a mixed method of analyzing shell elements and solid elements using the overlaying mesh method. In the structural design of a ship's hull, the shell elements are used for the global model. However, the solid elements are necessary to analyze the stress concentration zones or the vicinity of a crack. In such cases, the models are analyzed using zooming analysis, in which the results of a global model analysis are transferred to a local model analysis by imposing boundary conditions. This method is more advantageous than zooming analysis in terms of the accuracy of the solution and the modeling flexibility. Some examples of a plate model with a cracked surface or with a projection are shown in order to demonstrate the effectiveness of the method. Received: August 6, 2002 / Accepted: November 25, 2002 Address correspondence to: S. Nakasumi (sumi@nasl.t.u-tokyo.ac.jp) Updated from the Japanese original, which won the 2002 SNAJ prize (J Soc Nav Archit Jpn 2001;189:219–224; and 190:655–662)     

Carbonate chemistry dynamics and carbon dioxide fluxes across the atmosphere–ice–water interfaces in the Arctic Ocean: Pacific sector of the Arctic

Climatic changes in the Northern Hemisphere have led to remarkable environmental changes in the Arctic Ocean, which is surrounded by permafrost. These changes include significant shrinking of sea-ice cover in summer, increased time between sea-ice break-up and freeze-up, and Arctic surface water freshening and warming associated with melting sea-ice, thawing permafrost, and increased runoff. These changes are commonly attributed to the greenhouse effect resulting from increased atmospheric carbon dioxide (CO₂) concentration and other non-CO₂ radiatively active gases (methane, nitrous oxide). The greenhouse effect should be most pronounced in the Arctic where the largest air CO₂ concentrations and winter–summer variations in the world for a clean background environment were detected. However, the air–land–shelf interaction in the Arctic has a substantial impact on the composition of the overlying atmosphere; as the permafrost thaws, a significant amount of old terrestrial carbon becomes available for biogeochemical cycling and oxidation to CO₂. The Arctic Ocean's role in determining regional CO₂ balance has been ignored, because of its small size (only  4% of the world ocean area) and because its continuous sea-ice cover is considered to impede gaseous exchange with the atmosphere so efficiently that no global climate models include CO₂ exchange over sea-ice. In this paper we show that: (1) the Arctic shelf seas (the Laptev and East-Siberian seas) may become a strong source of atmospheric CO₂ because of oxidation of bio-available eroded terrestrial carbon and river transport; (2) the Chukchi Sea shelf exhibits the strong uptake of atmospheric CO₂; (3) the sea-ice melt ponds and open brine channels form an important spring/summer air CO₂ sink that also must be included in any Arctic regional CO₂ budget. Both the direction and amount of CO₂ transfer between air and sea during open water season may be different from transfer during freezing and thawing, or during winter when CO₂ accumulates beneath Arctic sea-ice; (4) direct measurements beneath the sea ice gave two initial results. First, a drastic pCO₂ decrease from 410 μatm to 288 μatm, which was recorded in February–March beneath the fast ice near Barrow using the SAMI-CO₂ sensor, may reflect increased photosynthetic activity beneath sea-ice just after polar sunrise. Second, new measurements made in summer 2005 beneath the sea ice in the Central Basin show relatively high values of pCO₂ ranging between 425 μatm and 475 μatm, values, which are larger than the mean atmospheric value in the Arctic in summertime. The sources of those high values are supposed to be: high rates of bacterial respiration, import of the Upper Halocline Water (UHW) from the Chukchi Sea (CS) where values of pCO₂ range between 400 and 600 μatm, a contribution from the Lena river plume, or any combination of these sources.     

Comparison of CFD and EFD for the Series 60 C B= 0.6 in steady drift motion

 This paper presents comparisons of computational and experimental fluid dynamics results for boundary layers, wakes, and wave fields for the Series 60 C _B= 0.6 ship model in steady drift motion. The numerical method solves the unsteady Reynolds-averaged Navier–Stokes and continuity equations with the Baldwin–Lomax turbulence model, exact nonlinear kinematic and approximate dynamic free-surface boundary conditions, and a body/free-surface conforming grid. The experimental and computational conditions, i.e., Froude numbers of 0.16 and 0.316 for the experiments, and Froude numbers of 0 and 0.316 for the computations, allow comparisons of low and high Froude number results, respectively, which allows an evaluation of Froude number effects and validation of the computational fluid dynamics at both low and high Froude numbers. This article gives an overview of this numerical approach, and the computational conditions and uncertainty analysis are described. Results are presented for the wave and flow fields, with emphasis on the important flow features of drift- and wave-induced effects in comparison with the experiments. Finally, conclusions from the present study are given, together with recommendations for future work. Received: August 31, 2001 / Accepted: March 25, 2002     

潜艇操纵性水动力数值计算中湍流模式的比较与运用

湍流模式的选取对潜艇操纵性水动力数值计算精度有着重要影响,采用六种湍流模式计算了SUBOFF主艇体及主艇体加指挥室围壳两种模型在一定漂角范围内的潜艇操纵性水动力,并与试验值进行了比较,结果表明SSTk-ω模型更为适合潜艇操纵性水动力计算。在此基础上对SUBOFF全附体模型在一定漂角和攻角范围内的艇体水动力进行预报计算,并对该计算方法应用于潜艇操纵性水动力预报计算的计算精度与适用范围进行了探讨。     

Use of SeaWiFS data for light availability and parameter estimation of a phytoplankton production model of the Bay of Biscay

Processing SeaWiFS (Sea-viewing Wide Field-of-view Sensor) data provides useful information for the observation and modelling of the phytoplankton production of the Bay of Biscay. Empirical algorithms allow the retrieval of chlorophyll a and non-living Suspended Particulate Matter (SPM) concentrations. These data are used to constrain a coupled 3D physical–biogeochemical model of the Bay of Biscay continental shelf. Two issues are investigated, depending on the variable used, to constrain the winter to spring phytoplankton production for the year 2001. First, SPM data is used as forcing data to correct the corresponding state variable of our model. This allows the realistic simulation of the light limited bloom at the end of February 2001, as observed with SeaWiFS chlorophyll a images and from the NUTRIGAS field cruise. Second, chlorophyll a data is used for parameter estimation of the biogeochemical model. The ability of assimilating these data is tested to improve the simulation of strong blooms observed in late May 2001 in the Loire and Gironde plumes. A global optimization method (Evolutive Strategies) is adapted to the complete 3-D coupled model, in order to find the best set of parameters. The hydrological conditions during the bloom can be validated with data from the PEL01 field cruise. After selection of the most sensitive parameters, the method is tested with twin experiments. Then, the use of real SeaWiFS data reduces the model/data misfit by a factor of two, improving the simulation of bloom intensities and extensions. The sets of parameters retrieved in each plume are discussed.     

URANS simulations of static and dynamic maneuvering for surface combatant: part 1. Verification and validation for forces, moment, and hydrodynamic derivatives

Part 1 of this two-part paper presents the verification and validation results of forces and moment coefficients, hydrodynamic derivatives, and reconstructions of forces and moment coefficients from resultant hydrodynamic derivatives for a surface combatant Model 5415 bare hull under static and dynamic planar motion mechanism simulations. Unsteady Reynolds averaged Navier–Stokes (URANS) computations are carried out by a general purpose URANS/detached eddy simulation research code CFDShip-Iowa Ver. 4. The objective of this research is to investigate the capability of the code in regards to the computational fluid dynamics based maneuvering prediction method. In the current study, the ship is subjected to static drift, steady turn, pure sway, pure yaw, and combined yaw and drift motions at Froude number 0.28. The results are analyzed in view of: (1) the verification for iterative, grid, and time-step convergence along with assessment of overall numerical uncertainty; and (2) validations for forces and moment coefficients, hydrodynamic derivatives, and reconstruction of forces and moment coefficients from resultant hydrodynamic derivatives together with the available experimental data. Part 2 provides the validation for flow features with the experimental data as well as investigations for flow physics, e.g., flow separation, three dimensional vortical structure, and reconstructed local flows.     

Self-propulsion computations using a speed controller and a discretized propeller with dynamic overset grids

A method that can be used to perform self-propulsion computations of surface ships is presented. The propeller is gridded as an overset object with a rotational velocity that is imposed by a speed controller, which finds the self-propulsion point when the ship reaches the target Froude number in a single transient computation. Dynamic overset grids are used to allow different dynamic groups to move independently, including the hull and appendages, the propeller, and the background (where the far-field boundary conditions are imposed). Predicted integral quantities include propeller rotational speed, propeller forces, and ship’s attitude, along with the complete flow field. The fluid flow is solved by employing a single-phase level set approach to model the free surface, along with a blended k−ω/k−ɛ based DES model for turbulence. Three ship hulls are evaluated: the single-propeller KVLCC1 tanker appended with a rudder, the twin propeller fully appended surface combatant model DTMB 5613, and the KCS container ship without a rudder, and the results are compared with experimental data obtained at the model scale. In the case of KCS, a more complete comparison with propulsion data is performed. It is shown that direct computation of self-propelled ships is feasible, and though very resource intensive, it provides a tool for obtaining vast flow detail.     

Benthic foraminifera in the surface sediments of the Beaufort Shelf and slope, Beaufort Sea, Canada: Applications and implications for past sea-ice conditions

This paper presents new data on distribution patterns of modern benthic foraminifera and other microfossils from the Canadian Arctic, specifically the Beaufort Shelf and slope. The material was collected in June to August of 2004 and is the first of its kind in this area to be collected since 1970. We examined the smaller sizes (45–63µm) as well as > 63µm and discovered that many species had been severely underrepresented in previous studies. Deep sea forms, that had been overlooked previously, were common on the shelf; two species (Elphidiella arctica and Ammotium cassis) appeared in preliminary results to be indicators of methane seepage; and it was possible to make determinations of sea-ice coverage using a combination of foraminifera and tintinnids (planktic ciliates). Our data indicated the presence of many of the same species as previous studies from this area, but improved techniques of sample processing greatly increased the number of specimens and species found (particularly the small deep sea arctic species Buliminella hensoni and Bolivina arctica) which provide much more reliable data for paleoceanographic determinations. One of the primary objectives for this work was to provide baseline data to help determine paleo-ice cover; these data cover a broad range of conditions on the Beaufort Shelf that make it possible to achieve this objective as well as improving what it is known about the assemblages on this shelf as compared to other arctic shelf areas, such as the Siberian Shelf).     

用naoe-FOAM-SJTU求解器数值模拟船舶在波浪上的操纵运动问题（英文）

Ship maneuvering in waves includes the performance of ship resistance, seakeeping, propulsion, and maneuverability. It is a complex hydrodynamic problem with the interaction of many factors. With the purpose of directly predicting the behavior of ship maneuvering in waves, a CFD solver named naoe-FOAM-SJTU is developed by the Computational Marine Hydrodynamics Lab(CMHL) in Shanghai Jiao Tong University. The solver is based on open source platform OpenFOAM and has introduced dynamic overset grid technology to handle complex ship hull-propeller-rudder motion system. Maneuvering control module based on feedback control mechanism is also developed to accurately simulate corresponding motion behavior of free running ship maneuver. Inlet boundary wavemaker and relaxation zone technique is used to generate desired waves. Based on the developed modules, unsteady Reynolds-averaged Navier-Stokes(RANS) computations are carried out for several validation cases of free running ship maneuver in waves including zigzag, turning circle, and course keeping maneuvers. The simulation results are compared with available benchmark data. Ship motions, trajectories, and other maneuvering parameters are consistent with available experimental data, which indicate that the present solver can be suitable and reliable in predicting the performance of ship maneuvering in waves. Flow visualizations, such as free surface elevation, wake flow, vortical structures, are presented to explain the hydrodynamic performance of ship maneuvering in waves. Large flow separation can be observed around propellers and rudders. It is concluded that RANS approach is not accurate enough for predicting ship maneuvering in waves with large flow separations and detached eddy simulation(DES) or large eddy simulation(LES) computations are required to improve the prediction accuracy.     

珠江口区域海上溢油动态预测信息系统的开发与应用

为了对《珠江口区域溢油应急计划》提供决策支持,在国内外相关研究成果基础上,针对珠江口海域的特点,研究开发了先进实用的"珠江口区域海上溢油动态预报信息系统",综合了三维潮流模型、三维溢油扩散模型、溢油风化模型、应急反应模型、以及电子海图、地理信息系统(GIS)、数据库等关键技术。该系统可以预测模拟并可视化显示海上溢油的漂移扩散和性质变化过程,同时显示环境敏感区和应急人员设备分布等相关信息。实际溢油应用案例表明,该系统的预报模拟结果与现场实际情况完全相符,能有效地提高海上溢油污染事故的应急决策效率。     

URANS simulations of static and dynamic maneuvering for surface combatant: part 2. Analysis and validation for local flow characteristics

Part 2 of this two-part paper presents the analysis and validation results of local flow characteristics for a surface combatant Model 5415 bare hull under static and dynamic planar motion mechanism simulations. Unsteady Reynolds averaged Navier–Stokes (URANS) computations are carried out by a general-purpose URANS/detached eddy simulation research code CFDShip-Iowa Ver. 4. The objective of this research is to investigate the capability of the code in relation to the computational fluid dynamics-based maneuvering prediction method. In the current study, the ship is subjected to static drift, steady turn, pure sway and pure yaw motions at Froude number 0.28. The free surface, three dimensional vortical structure and, the validation of two dimensional local flow quantities together with the available experimental data are of the interest in the current study. Part 1 provides the verification and validation results of forces and moment coefficients, hydrodynamic derivatives, and reconstructions of forces and moment coefficients from resultant hydrodynamic derivatives.