A low-order potential based boundary element method is applied for the prediction of the performance of flapped rudders as well as all-movable rudders in steady inflow. In order to obtain a reasonable solution at large angles of attack, the location of the trailing wake sheet is determined by aligning freely with the local flow. The effect of the wake sheet roll-up is also included with use of a high order panel method. The flow in the gap of a flapped rudder is modeled as Couette flow and its effect is introduced into the kinematic boundary conditions for flux at both the inlet and the outlet of the gap. In order to validate the present method, the method is applied for a series of rudders and the computational results on forces and moments are compared with experimental data. The effect of the gap size on the forces and moments is also presented.

Reynolds-averaged Navier-Stokes equations are numerically solved using a secondorder finite difference method for the analysis of turbulent flows around single and tandem hydrofoils advancing under the free surface. The location of the free surface, not known a priori, is computed from the kinematic free surface condition and the computational grid is conformed at each iteration to the free surface deformation. The eddy viscosity model of Baldwin-Lomax is employed for the turbulence closure. The method is validated through the comparision of the numerical results with the experimental data for a single hydrofoil of a Joukowski foil section. A computational study is also carried out to investigate the effect of the submergence depth and the Froude number on the lift and the drag of the hydrofoil. For tandem hydrofoils, computations are performed for several separation distances between the forward and aft foils to see the interference effect. The result shows clearly how the lift and drag change with the separation distance.

The slamming force occurring in the free fall impact of cylindrical bodies on the water surface is analyzed in both compressible and incompressible stages. In the compressible phase the hydrodynamic analysis is carried on by the acoustic approximation, obtaining a closed form expression for the maximum impact force. The incompressible analysis is approached through and unsteady boundary element method to compute the free surface evolution and the slamming force on the body. A similar behavior seems to characterize the maximum slamming force versus a dimensionless mass parameter.

Numerical simulation based on the superposition of ring waves generated by the linear periodic source distributions for the plunger type wave maker was accomplished. The characteristics of directional spreading function were investigated. Mirror images are also introduced to consider reflections of side-wall together with the reflection coefficient to account for the imperfect reflection from the real side wall in the long experimental towing tank. Unexpected spurious waves, resulting from the combined effect of finite breadth of segmented wave maker, wavelength and main wave maker, wavelength and main wave propagating direction, were observed in the line source method and also in the analysis of the directivity. The influence of spurious waves to the directional spreading function was also investigated.