the inviscid case. Thus, converged NS runs require about three hours on eight processors.

The establishment of an efficient baseline multiblock flow solver is extremely important because it enables both sea-keeping analysis and optimization to become practical tools in the design process.

For a viscous incompressible fluid moving under the influence of gravity, the differential form of the continuity equation and the Reynolds Averaged Navier-Stokes equations (RANS) in a Cartesian coordinate system can be cast, using tensor notation, in the form,

Here, is the mean velocity components in the *x _{i}* direction, the mean pressure, and the gravity force acting in the

where is the Froude number and the Reynolds number *Re* is defined by where *ν* is the kinematic viscosity, and is a dimensionless form of the Reynolds stress.

Figure 1 shows the reference frame and ship location used in this work. A right-handed coordinate system *Oxyz,* with the origin fixed at the intersection of the bow and the mean tree surface is established. The *z* (*x*_{3}) direction is positive upwards, *y* (*x*_{2}) is positive towards the starboard side and *x* (*x*_{1}) is positive in the aft direction. The free stream velocity vector is parallel to the *x* axis and points in the same direction. The ship hull pierces the uniform flow and is held fixed in place, ie. the ship is not allowed to sink (translate in *z* direction) or trim (rotate in *x−z* plane).

It is well known that the closure of the Reynolds averaged system of equation requires a model for the Reynolds stress. There are several alternatives of increasing complexity. Generally speaking, when the flow remains attached to the body,

a simple turbulence model based on the Boussinesq hypothesis and the mixing length concept yields predictions which are in good agreement with experimental evidence. For this reason a Baldwin and Lomax turbulence model has been initially implemented and tested [5]. On the other hand, more sophisticated models based on the solution of additional differential equations for the component of the Reynolds stress may be required. Notice that when the Reynolds stress vanishes, the form of the equation is identical to that of the Navier Stokes equations. Also, the inviscid form of the Euler equations is recovered in the limit of high Reynolds numbers. Thus, a hierarchy of mathematical model can be easily implemented on a single computer code, allowing study of the controlling mechanisms of the flow. For example, it has been shown in reference [9] that realistic prediction of the wave pattern about an advancing ship can be obtained by using the Euler equations as the mathematical model of the bulk flow, provided that a non-linear evolution of the free surface is accounted for. This is not surprising, since the typical Reynolds number of an advancing vessel is of the order of 10^{8}.

When surface tension as well as tangential stresses are neglected, the boundary condition on the free surface consists of two equations. The first, the dynamic condition, states that the pressure acting on the free surface is equal to the stresses normal to the free surface. It was found that the inclusion of the viscous stresses had little to no effect on the solution since they are of the order . Thus the dynamic condition that pressure is a constant is a good approximation. The kinematic condition states that the free surface is a material surface: once a fluid particle is on the free surface, it forever remains on the surface. The dynamic and kinematic boundary conditions may be expressed as

*p=constant*