Modelling and Simulation in Engineering

Volume 2011, Article ID 356741, 10 pages

http://dx.doi.org/10.1155/2011/356741

## A Fully Nonlinear, Dynamically Consistent Numerical Model for Ship Maneuvering in a Seaway

^{1}Hydromechanics Department, David Taylor Model Basin, NSWCCD, 9500 MacArthur Boulevard, West Bethesda, MD 20817-5700, USA^{2}Planetary Geodynamics Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD, USA

Received 28 March 2011; Accepted 19 July 2011

Academic Editor: Agostino Bruzzone

Copyright © 2011 Ray-Qing Lin and Weijia Kuang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This is the continuation of our research on development of a fully nonlinear, dynamically consistent, numerical ship motion model (DiSSEL). In this paper we report our results on modeling ship maneuvering in arbitrary seaway that is one of the most challenging and important problems in seakeeping. In our modeling, we developed an adaptive algorithm to maintain dynamical balances numerically as the encounter frequencies (the wave frequencies as measured on the ship) varying with the ship maneuvering state. The key of this new algorithm is to evaluate the encounter frequency variation differently in the physical domain and in the frequency domain, thus effectively eliminating possible numerical dynamical imbalances. We have tested this algorithm with several well-documented maneuvering experiments, and our results agree very well with experimental data. In particular, the numerical time series of roll and pitch motions and the numerical ship tracks (i.e., surge, sway, and yaw) are nearly identical to those of experiments.

#### 1. Introduction

Ship motion with fixed heading is an overly simplification of general ship motion dynamics, for example, ship maneuvering in seaway is far more complex than the fixed heading. In our earlier work, Lin and Kuang [1] developed a fully nonlinear, dynamically consistent ship motion model (DiSSEL) that can accurately simulate ship motion with fixed heading. A continuing effort is therefore to expand DiSSEL for simulating ship motions when maneuvering in seaway.

Accurate prediction of ship motion when maneuvering in a seaway is one of the most challenging and important problems in ship motion dynamics. Much effort has been devoted in the past on modeling ship maneuvering in seaway. The earliest studies started from the simplest cases of simulation ship tracks (i.e., surge, sway, and yaw motions) in calm water. For example, Chaplin [2] and Folger et al. [3] predicted ship tracks in calm water based on empirical laws derived from experimental data. Söding [4] made a substantial step forward by including dynamics in ship motion modeling. He introduced a set of hydrodynamic maneuvering coefficients in simulation, though those coefficients are determined mainly from experimental data. Hooft and Pieffers [5] furthered this effort by including the effects of waves on ship motion (e.g., wave induced motions). In 1988, Lee [6] developed a potential flow method to calculate the rudder forces and the moments, thus making it possible to fully simulate ship maneuvering motion state. Nearly a decade late, Bailey et al. [7] proposed a hybrid approach in which experimental data and theoretical formulations (e.g., linearized or parameterized ship response function) are integrated into a single model, thus providing an enhanced description of fluid loading. More recently, several complicated computational models are developed for simulating complex cases, for example, the viscous flow models by Chau [8] and El Moctar [9].

These studies have greatly improved our understandings of ship tracks and ship maneuvering dynamics. The methodologies developed from these efforts are still widely used [10–15]. However, as pointed by Söding [10], there are also many limitations in these methods. For example, the potential flow methods ignore the dissipative effects from (tiny) fluid viscosity, small-scale turbulence, and flow separation. The viscous flow models are capable of evaluating these dissipative effects. But they are computationally difficult and expensive: numerical instabilities and computational capability limitations make direct simulation in the oceanic parameter domain impossible. In addition, inclusion of any a priori condition (e.g., those empirical extrapolations from experimental data) will affect generality of the model.

To utilize the mathematical simplicity of potential flow approaches and to include appropriate dissipative effects, Lin and Kuang [16] developed a physically more realistic and computationally more efficient model in which the “blocking theory” is used to numerically determine the rudder forces, the lift forces, the moments, and the Kutta conditions at the trailing edges. Based on this approach, Lin et al. [17] were able to model accurately ship track and heel motion when maneuvering in calm water. Their results are well benchmarked by experiments data.

The natural extension of the effort by Lin and Kuang [16] and by Lin et al. [17] is to develop a dynamically consistent model for predicting ship maneuvering in seaway. The model by Lin and Kuang [1] can be a very good candidate for this purpose for several reasons: there is no a priori condition nor parameterization in this model, and the model is fully nonlinear and dynamically consistent. In particular, prediction of ship motion in fixed heading from this model agrees well with experiments. For predicting ship maneuvering in seaway, several mechanisms need to be included (or reintroduced) into this model, including change of ship heading, motion acceleration, and propeller forces.

But, one might anticipate a mathematical difficulty of modeling encounter waves in ship maneuvering. The wave frequency and the wave number vector of the encounter waves vary with the ship speed and ship heading. The variation due to the latter is in particular significant in ship maneuvering. As we will find later in this paper, changes in the encounter frequencies do have profound impacts on ship-motion interactions and thus on the ship motion state. In this paper, we will discuss the details of this model expansion, the mathematical/numerical difficulties, and the algorithms resolving these difficulties. This paper is organized as follows: in Section 2 is an overview of DiSSEL. Section 3 describes the details of the relevant algorithm. Model results validation is given in Section 4, followed by the discussions.

#### 2. Fundamentals of DiSSEL Nonlinear Ship Motion Model

DiSSEL is a hybrid-flow ship motion model in time domain: the surface waves are described by velocity potential, but the fluid viscous effect is included in the form of wave breaking and wave blocking (separation effects) [1, 16, 18]. Inclusion of the viscous effect makes DiSSEL satisfying the Kutta condition at a trailing edge, and velocity stagnation points (), thus allowing accurate evaluation of the lift force [17]. This hybrid-flow approach builds a solid foundation of our effort to expand the model for ship maneuvering in a seaway.

In DiSSEL, ship-wave interaction and ship solid body motion dynamics are modeled by two separate components. Interactions between the two dynamical processes are determined through the pressure force on the ship hull surface and the restoring force from the water displaced by underwater ship volume. For the details of the partial differential equations, the boundary conditions and the surface integrations, we refer the reader to Lin et al. [19], Lin and Kuang [1, 20].

To expand DiSSEL for a ship maneuvering in a seaway, it is necessary to include a propulsion mechanism. As the first step, this mechanism can be approximated from experimental data (empirical). This approximation does not affect the dynamical consistency of our model, since the propulsion is determined by the ship engine, independent of the ship-wave interactions. In the following, we focus on the differences in the equations and the physical variables that are different from those of the original DiSSEL model.

Another notice should be made on the definitions of the reference frames used in DiSSEL. In Lin and Kuang [1], two additional reference frames are introduced: the ship reference frame (SRF) and the model reference frame (MRF). The latter is defined to be the reference frame moving horizontally with the ship mass center, but with the origin on the mean free surface. When a ship maneuvers in seaway, its track direction can change. But MRF does not change its orientation. Therefore, the rotational effect from this direction change needs to be accounted for correctly. For details of the reference frames, we refer the reader to Lin and Kuang [1].

##### 2.1. Ship-Wave Interaction Model

In DiSSEL, the surface waves are described by the following set of partial differential equations: where is the pressure, is the water density, is the turbulent viscosity (arising from the wave-breaking processes), is the position vector, is the water depth, is the velocity potential, is the free surface elevation, is the gradient, is the horizontal component part of , and is the ship maneuvering velocity (in SRF, it is the given forward speed). It should be noted that includes both magnitude and direction changes of the maneuvering velocity. In addition, is the horizontal component of translational velocity for surge and sway (absent in the fixed heading scenario).

The impenetrable boundary conditions are at the flat bottom , and on the ship surface . In the above equation, is the normal unit vector of the ship surface , and is the mass center of the ship hull.

##### 2.2. Solid Body Motion Model

The ship solid body motion is decoupled into the translational motion of ship mass center and the rotation about . The translational motion, , is described by the following equation in MRF: where is the total ship mass and is the total force on the ship: where is the (driving) pressure force on the ship hull,

is the buoyancy force (restoring force). This force acts only in the vertical direction. But the rudder force in (4b) is the new addition. We will discuss this in more detail in the next section. is under water volume.

The solid body rotation about the ship mass center is governed by the Liouville equation defined in SRF: where is the moment of inertia tensor about the ship mass center (and its components are denoted by ), and are integrated dissipative effects (e.g., drag) of the fluid on the ship hull, and Γ is the total moment (torque) on the ship: It includes the (driving) pressure torque: (where is total wetted ship surface), the restoring torque:

(where is the position vector in under water volum), and the torque from the rudder (a new addition). It should be pointed out that the restoring torque is always horizontal, thus affecting only roll and pitch motions.

##### 2.3. Rudder Force and Moment

The rudder force in (4b) and the moment in (5b) are the new physical variables included in this expansion. Since the rudder position is often defined relative to the ship hull, it is therefore very convenient to evaluate them in SRF. Therefore, this leads to the following formulations: where is the potential flow pressure; but ( is the wetted rudder surface, and is the blocking effected area) [16], is the unit normal vector of , is the position vector of the center of the circle tangent to the point on the ship trajectory, is the radius vector of ship trajectory to the point , and is the rotating arm vector of the rudder. The force in (7) describes the effect of maneuvering directional changes in SRF and is, therefore, dependant on the changes in the rudder angle (defined relative to the ship hull): where is the drifting angle vector (the angle between the ship’s horizontal velocity vector and the ship’s longitudinal axis). If the rudder can be approximated as a plate rotating about the vertical direction of the ship hull with a rudder deflection angle , the horizontal components of can be substantially simplified: where is the effective inflow speed; and are the and components of the normal vector of the rudder surface, respectively; is pressure on rudder surface due to the incident waves; is the lift curve slope: where and are the length and the width of the rudder, respectively [21]. This theory has been tested by many experiments and is broadly used in ship dynamics. It should be pointed out that the separation effects are included by using (not ) in (6)–(10).

By (10), the vertical moment is since is in general much smaller than . In (12), is the longitudinal distance between the ship mass center and the rudder mass center .

The effective inflow into the rudder is derived from the ship forward speed and the accelerated flow in the propeller slipstream. The latter depends on the propeller’s mechanical design. In the cases discussed in this paper, we choose an empirical estimate from experimental data for naval vessels with Froude numbers ranging from 0.2 to 0.4 [22].

##### 2.4. Ship Track Kinematics

Kinematics of ship track variation can be best described in an inertial reference frame, for example, the earth-fixed reference frame (ERF) (). The ship moving direction angle, , is the sum of the yaw angle (the rotational angle in direction, and obtained from (5a)–(5d)) and the drift angle : where and are ship translation motion in the and directions in ship coordinate, respectively. The sketch map can be find in Figure 5 page 198 [23, Figure 5, page 198]. The track of the ship in ERF is defined as In particular,

If the rudder deflection angle is given, then the rudder force (10) and moment (12) can be determined. If the ship velocity is calculated, for example, via (4a)–(4d) and (5a)–(5d), then the track can be determined via time integration of (16). On the other hand, given the navigation history, that is, the ship track , the total ship horizontal velocity, , and horizontal acceleration can be calculated by simple time differentiation, and the direction angle can be also determined via The the rudder angle can be obtained from horizontal component of (4a)–(4d), vertical component of (5a)–(5d), (10), and (16) using an iterative method. Finally we want to point out that, in DiSSEL ship motion model, hydrodynamic interactions between the hull and the rudder are properly implemented, in addition to the implementation of the interactions of the six-degree freedom of the ship motions. These nonlinear interactions are obtained by solving a set of the nonlinear equations (1)–(7).

##### 2.5. Centrifugal Effect

Ship track direction changes generate a finite centrifugal force. In a very simple mathematical description, the force can be described as a vertical rotation about the arc center of the track curve: where is the rotation speed about and is the position vector of the point in ERF. Therefore, Thus, In other words, the changes in the direction of the ship motion can be described by the effective pressure in (20). Given the ship track , the radius can then be calculated with the standard curvature formulation: Its influence on the solid body rotation of the ship can be computed via integration of the effective pressure over the ship surface.

#### 3. Stable and Conservative Algorithm for Ship Maneuvering

When a ship maneuvers in seaway, its direction and speed change over time. This leads to time-varying encounter waves, thus potentially to numerical difficulties in conserving the momentum of the entire system. This will be a source of numerical instabilities. Therefore, we will discuss first the nature of the numerical difficulties (based on current ship motion models) and then introduce our new stable algorithm that conserves the momentum of the system.

##### 3.1. Dynamical Imbalance from Rapidly Changing Encounter Frequencies

As pointed out by Lin and Kuang [24], the numerical balance between the driving forces and the restoring forces is the key to model correctly ship motions. Lin and Kuang [1] further showed that DiSSEL can well maintain such dynamic balances for fixed heading of various ship hulls in a wide range of environments.

To expand DiSSEL for modeling ship maneuvering, a natural consideration is whether the algorithm used in the model is able to maintain accurately the dynamic balance between the driving forces and the restoring forces on the ship hull when ship heading changes in time. The consideration arises from the very fact that such changes also alter the environment observed on the ship. For example, in SRF, a surface wave can be generally described as where is the amplitude of the mode , is the wave number vector, is the encounter frequency, and is the phase. The integer pair () defines the total number of the wave modes used in numerical simulation. The encounter frequency is defined as where is the phase velocity of the wave mode, is the angle between and , and is the incident wave frequency in ERF. For any given incident wave, changes with the ship heading direction.

One could clearly observe the complications and its numerical difficulties from (23): changes in and (and therefore ) are dynamically determined via the driving forces and the restoring forces on the ship hull which depend on the encounter frequency . But changes in depend also on and . Therefore, errors in the force balances could potentially lead to numerical instabilities and false numerical solutions in ship motion simulation.

To help understand this, we consider a very mathematically simple case: a ship maneuvering in a single environment wave () (and is called Case 1 in the rest of the paper). The ship heading direction changes at a steady rate (with the rudder angle ). The value of the ship speed is also fixed and slow, with the Froude number . We first model this simple system with the existing version of DiSSEL [1]. The numerical driving forces (moments) and the restoring forces (moments) for heave, pitch, and roll are shown in Figure 1. From the figure one can observe that, as soon as the ship direction starts to change, the balance between the forces (moments) is lost numerically.

What is the cause of this imbalance? One possibility is the time variation of the encounter frequency . To verify this, we consider Case 2: an artificial system in which is fixed, but the ship heading direction change remains the same as in the previous case. Numerical results for this case are shown in Figure 2. Clearly, the force balance is well maintained for this artificial system. Of course, this is not a mathematical proof, but it is sufficient to demonstrate the key problem for retaining numerical force balances in ship motion modeling.

##### 3.2. Adaptive Method for the Time-Varying Encounter Frequencies

The above numerical experiments can guide our effort in developing a new algorithm to maintain dynamic balance for time-varying encounter frequencies. In particular, we intend to exploit the fact that the force balances can be well achieved for a steady encounter frequency.

A simple and effective approach is to apply the same numerical algorithm in a “fictitious” reference frame in which the encounter frequency is time invariant. To explain this approach, we will start from a single-wave scenario.

In this case, we introduce a “fictitious” time and its corresponding time step , such that From (23), But (25) could lead to when that is,

To avoid this numerical singularity, we modify (25) by adding a truncation limit: Substituting (28) into (22) and setting (single-model wave), we have where

and , .

To test the applicability of this approach, we reexamine Case 1. The force balance is now well retained with the new algorithm as shown in Figure 3. By comparing this with Figure 1, we can observe clearly the improvement by the new algorithm.

We test further the new algorithm with more complicated cases that have experimental data. Due to space limitations, we show here only two of them: the EHF maneuvering in regular waves with the largest rudder angle (Case 3) and with the fastest forward speed (Case 4). The results for the two cases are shown in Figures 4 and 5. Because in the experiments, the ship motions, wave elevations, and trajectory motions were collected by separate systems, all the data sets could not be aligned in time accurately enough to know the initial conditions of the system. Therefore, the benchmarking comparison focuses only on the overall characteristics of the motion and not a point-by-point prediction comparison. As shown in these figures, our numerical results agree very well with the experimental data. The good agreement is because the dynamical balances are now correctly modeled in the simulation, as shown in Figures 4(c) and 5(c).

This new algorithm can be easily extended to arbitrary two-dimensional irregular waves. The mathematical details are shown as follows: expression (29a) is now modified to irregular waves: where and are defined similarly as in (28) for each wave mode and .

To maintain computational efficiency, we introduce an effective incident wave frequency and wave amplitude: With these definitions, we can then define a single effective “fictitious” time step The notions of the effective incident wave frequency and the “fictitious” time are the same as those in the one-dimensional regular wave formulations. But the actual physical meanings are different. In particular, is given and invariant throughout the simulation. But changes in time because of time-varying . Therefore, (30) can be rewritten as where Again, and (and therefore by definition) must satisfy the CFL condition in the simulation.

It should be pointed that this method is more efficient than iterative methods not only because it calculates the mean once in each time step, but also because it does not impose any convergence requirement and thus ensures a stable solution.

Since there is no experimental data on ship maneuvering in irregular seas with which to compare, we will only examine the dynamical balances in the numerical solutions, with no attempt made for benchmarking purposes. For this, we reconsider Case 4, but replace the regular wave with the irregular wave specified by the spectral shape in Figure 6. The numerical results are summarized in Figure 7. From the figures we can observe that the dynamical balances are well established in the numerical solutions for roll, pitch, and heave motions, that is, the driving forces (moments) are comparable in magnitude and opposite in sign to the restoring forces (moments).

#### 4. Conclusion

When a ship maneuvers in a seaway, arbitrary incident waves measured on the ship are time varying (i.e., time-varying encounter frequencies). This time variation depends on the ship motion state, which, in turn, changes in time under the influence of the incident waves. Therefore, it is critical to maintain the dynamical balances between the driving forces (moments) and the restoring forces (moments) in a numerical simulation for the accurate modeling of ship motion. Such dynamical balance depends on the appropriate evaluation of the incident waves and the corresponding forces on the ship hull and thus requires well-designed numerical approaches to model the nonlinear dependencies between the environment and the ship motion.

We have demonstrated numerically that direct time integration of the ship-wave interactions cannot retain properly the balance when the ship heading changes in time. This heading change results in time-varying encounter wave frequencies. We anticipated that this time-varying encounter frequency could lead to phase misfits between the driving forces and the restoring forces, thus destroying the dynamical balances in the simulation and therefore the numerical stability. This is supported by our numerical simulations of Cases 1 and 2 (shown in Figures 1 and 2).

To eliminate the time variation of the encounter frequency (mainly due to ship heading direction change), we developed a new adaptive algorithm with a mathematical (but “fictitious”) time defined in (32). This is equivalent to rescaling the numerical time (and thus improving numerical temporal resolution) dynamically with the ship motion states. In this new system, the encounter frequencies are steady as measured in . We have tested this new algorithm for many cases. All numerical results show that the dynamical balances between the driving forces (moments) and the restoring forces (moments) are well established.

This new algorithm has been implemented in the DiSSEL nonlinear ship motion model. With this improved model, we are able to reproduce satisfactorily the results from all known laboratory experiments. Results for two of such experiments are shown in Figure 4 (with the largest rudder angle) and Figure 5 (with the largest forward speed) in regular waves. We also tested the improved model for ship maneuvering in irregular waves specified in Figure 6, and the dynamical balances are retained as well, as shown in Figure 7.

#### Acknowledgments

The authors thank Terry Applebee, Y. T. Shen, and J. Gorski for their comments and suggestions to this study. They also thank the Seakeeping Division at NSWCCD for providing the experimental data to benchmark and validate their model. This work is supported by Naval Surface Warfare Center Independent Laboratory In-House Research (ILIR) Program.

#### References

- R.-Q. Lin and W. Kuang, “Solid-body motion in fully nonlinear ship motion model,”
*Proceedings of the Royal Society A*. In press. - H. R. Chaplin, “Wind-tunnel tests of low-aspect-ratio control surfaces with trailing edge to wind,” Tech. Rep. 944, David Taylor Model Basin Report, 1954. View at Google Scholar
- L. W. Folger, D. E. Whicker, and L. F. Fehlner, “Free-stream characteristics of a family of low-aspect-ratio, all-movable control surfaces for application of ship design,” Tech. Rep. 933, David Taylor Model Basin Report, 1958. View at Google Scholar
- H. Söding, “Prediction of ship steering capabilities,”
*Schiffstechnik*, vol. 29, no. 1, pp. 3–29, 1982. View at Google Scholar · View at Scopus - J. Hooft and J. B. Pieffers, “Maneuverability of frigates in waves,”
*Marine Technology*, vol. 25, no. 4, pp. 262–271, 1988. View at Google Scholar · View at Scopus - J.-T. Lee, “A potential based panel method for the analysis of marine propellers in steady flow,” Tech. Rep. 87-13, Department of Ocean Engineering, MIT, 1988. View at Google Scholar
- P. A. Bailey, W. G. Price, and P. Temarel, “A Unified mathematical model descriping the manuoeuvring of a ship travelling in a seaway,”
*The Royal Institution of Naval Architects -RINA*, pp. 131–149, 1997. View at Google Scholar - S.-W. Chau, “Computation of rudder force and moment in uniform flow,”
*Ship Technology Research*, vol. 45, no. 1, pp. 3–13, 1998. View at Google Scholar · View at Scopus - O. M. El Moctar, “Numerical determination of rudder forces,” in
*Proceedings of the European Mechanics Society (Euromech '98)*, vol. 374, Poitiers, France, 1998. - H. Söding, “Limits of potential theory in rudder flow predictions,” in
*Proceedings of the 22nd Symposium on Naval Hydrodynamics*, Weinblum Lecture, pp. 622–637, Washington, DC, USA, 1999. - R. Kimoto, T. Katayama, and Y. Ikeda, “Effects of running attitude on hydrodynamic forces for oblique towed planing craft,” in
*Proceedings of the 2nd Asia-Pacific Workshop on Hydrodynamics (APHydro '04)*, pp. 115–122, Busan, Korea, 2004. - T. Katayama, R. Kimoto, and Y. Ikeda, “Effects of running attitudes on maneuvering hydrodynamic forces for planing hull,” in
*Proceedings of the International Conference on Fast Sea Transportation (FAST '05)*, St. Petersburg, Russia, 2005. - T. Katayama, T. Taniguchi, H. Fujii, and Y. Ikeda, “Development of maneuvering simulation method for high speed craft using hydrodynamic forces obtained from model tests,” in
*Proceedings of the International Conference on Fast Sea Transportation (FAST '09)*, Athens, Greece, 2009. - M. Altosole, M. Figari, and M. Viviani, “6 DOF Simulation of maneuvering and propulsive performance of water jet propelled Mega Yacht,” in
*Proceedings of the International Conference on Fast Sea Transportation (FAST '09)*, Athens, Greece, 2009. - A. C. Hochbaum, F. Stern, K. Agdrup et al., “Final report and recommendations to 25th ITTC,” in
*Proceedings of the 25th Information and Telecommunication Technology Center*, vol. 1, pp. 143–208, 2008. - R.-Q. Lin and W. Kuang, “Modeling nonlinear roll damping with a self-consistent, strongly nonlinear ship motion model,”
*Journal of Marine Science and Technology*, vol. 13, no. 2, pp. 127–137, 2008. View at Publisher · View at Google Scholar · View at Scopus - R.-Q. Lin, M. Hughes, and T. Smith, “Prediction of ship steering capabilities with a fully nonlinear ship motion model. Part 1: maneuvering in calm water,”
*Journal of Marine Science and Technology*, vol. 15, no. 2, pp. 131–142, 2010. View at Publisher · View at Google Scholar · View at Scopus - R.-Q. Lin and W. Kuang, “Modeling the effects of ship appendages on the six-degree of freedom ship motions,” in
*Proceedings of the 2nd International Conference, on Marine Research and Transportation*, pp. 25–32, 2007. - R. Q. Lin, W. Kuang, and A. M. Reed, “Numerical modeling of nonlinear interactions between ships and surface gravity waves I: ship waves in calm water,”
*Journal of Ship Research*, vol. 49, no. 1, pp. 1–11, 2005. View at Google Scholar · View at Scopus - Lin R.-Q and W. Kuang, “Numerical modeling of nonlinear interactions between ships and surface gravity waves II: ship boundary condition,”
*Journal of Ship Research*, vol. 50, no. 2, pp. 181–186, 2006. View at Google Scholar - S. F. Hoerner,
*Fluid Dynamic Lift: Theoretical, Experimental and Statistical Information*, Hoerner Fluid Dynamics, 1965. - T. C. Shen, “Private communication,” 2009.
- E. Lewis,
*Principles of Naval Architecture Second Revision*, The Society of Naval Architects and Marine Engineers, New Jersey, NJ, USA, 1998. - R.-Q. Lin and W. Kuang, “A nonlinear method for predicting motions of fast ships,” in
*Proceeding of 10th International Conference on Fast Sea Transportation*, Athens, Greece, October 2009.