• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.195-199
• /
• 2005

Wave breaking phenomenon near the fore body of a ship is numerically simulated. The ship advance with uniform velocity in calm water. For the simulation, incompressible Navier-Stokes equations and continuity equation are adopted as governing equations. The simulation is carried out in staggered variable mesh system with finite difference method. Marker and Cell(MAC) method and Marker-Density method are employed to track the free surface. Body boundary conditions are satisfied with the adoption of porosity method and no-slip condition on the hull surface. The ship model has a wedge type fore-body, and the computational domain is an appropriate region around the fore-body. The computation results are compared with some experimental results. Also the difference of the free surface tracking methods are discussed.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.4
• /
• pp.18-22
• /
• 2011

The Ranking source panel method was used to predict the flow phenomenon of a ship with a goose-neck type bulbous bow penetrating the free surface. The non-linearity of the free surface boundary condition was fully satisfied using an iterative calculation method, and the raised panel method was adopted to obtain a more stable solution at each iteration step. The panel cutting method was applied to generate a hull calculation grid at each iteration step, including the first step. At that time, the nose of the goose-neck type bulbous bow was divided by the free surface and the free surface panel was modified at each iteration step using the variable free surface panel method. Numerical calculations were performed to investigate the validity and efficiency of the applied numerical algorithm using the 3600 TEU container carrier. The computed wave resistance coefficients were compared with the experimentally achieved residual resistance coefficients.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.2 s.42
• /
• pp.71-80
• /
• 2005

The solution to a liquid sloshing problem is challenge to the field of engineering. This is not only because the dynamic boundary condition at the free surface is nonlinear, but also because the position of the free surface varies with time in a manner not known a priori. Therefore, this nonlinear phenomenon, which is characterized by the oscillation of the unrestrained free surface of the fluid, is a difficult mathematical problem to solve numerically and analytically. In this study, three-dimensional boundary element method(BEM), which is based on the so-called an arbitrary Lagrangian-Eulerian(ALE) approach for the fluid flow problems with a free surface, was formulated to solve the behavior of the nonlinear free surface motion. An ALE-BEM has the advantage to track the free surface along any prescribed paths by using only one displacement variable, even for a three-dimensional problem. Also, some numerical examples were presented to demonstrate the validity and the applicability of the developed procedure.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.80-89
• /
• 1992

This paper deals with the open boundary problems, and two numerical schemes are used for the implementation of open boundary condition. One is to add the artificial damping term to dynamic free-surface boundary condition. Determination of suitable damping coefficient and the damping cone is the most important in this scheme. The other scheme is a modified Orlanski's method. This will be useful for the problems with unidirectional waves. A few typical free-surface wave problems are modeled for the numerical test. Method of solution is fundamental source-distribution method and the fully nonlinear boundary conditions are applied. The computed results are compared with those of others for the proof of practicality of these schemes.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.4
• /
• pp.12-20
• /
• 1997

A potential-based panel method has been developed to investigate the resistance characteristics of a high speed catamaran advancing on the free surface. Normal dipoles and sources are distributed on the body surface while sources are distributed on the free surface. Linearised free surface conditions are used in the present analysis. To avoid the instabilities due to the velocity difference between inner and outer flow of a high speed catamaran, Kutta condition has been applied at the stern. Model test has been carried out not only to validate the numerical results but to confirm the capabilities of a CWC(Circulating Water Channel). It is believed that we can obtain the qualitatively reasonable results in the CWC. Computed results are compared with those of experiments and Insel's experimental values. Since the Kutta condition is applied at the stern, stable solutions are obtained at the high speed range. The present method, using linearised free surface conditions at the high speed range, seems to be a useful tool in the hull form design of a high speed catamaran.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.21-32
• /
• 1992

Numerical solutions are presented for solving the free surface flow created by a three-dimensional body moving beneath the free surface with constant velocity at an angle of attack. The solution is obtained using a panel method based on the perturbation potential, which employs Havelock sources and normal dipoles distributed on the body surface and Havelock normal dipoles in the wake downstream of the trailing edge. A pressure Kutta condition with an iterative solution procedure is implemented to satisfy equal pressure condition on the upper and lower surfaces at the trailing edge. Numerical calculation examples in the present paper include an ellipsoid at zero angle of attack, a rectangular planform wing at a small angle of attack in the limit of zero Froude number and then free surface flows and hydrodynamic forces acting on the submerged spheroid and parabolic strut are calculated. Discussions are made about the validity of the present method.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.2
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.3
• /
• pp.204-211
• /
• 1993

Boundary element method is applied to simulate nonlinear water waves using Green's identity formula in a numerical wave flume. A system of linear equations is formulated from the governing equation and free surface boundary conditions in order to calculate velocity potential and water surface elevation at each nodal point. The velocity square terms are included in the dynamic free surface boundary condition. The free surface is treated as a moving boundary. the vertical variation of velocity potential being considered in calculating the time derivative of the velocity potential at the free surface. The present method is applied to simulate solitary wave and Stokes 2nd order wave, and shows excellent agreements with their theoretical values.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.35 no.2
• /
• pp.41-48
• /
• 2023

In this study, simplified linear numerical method that can simulate wave generation and transformation by a moving bottom is introduced. Numerical analysis is conducted in wave number domain after continuity equation, linear dynamic and kinematic free surface boundary conditions and linear kinematic bottom boundary condition are Fourier transformed, and the results are expressed in space domain by an inverse Fourier transform. In the wavenumber domain, the dynamic free water surface boundary condition and the kinematic free water surface boundary condition are numerically calculated, and the velocity potential in the mean water level (z = 0) satisfies the continuity equation and the kinematic bottom boundary condition. Wave generation and transformation are investigated when the triangular and rectangular shape of bottoms move periodically. The results of the simplified numerical method are compared with the results of previous analytical solutions and agree well with them. Stability of numerical results according to the calculation time interval (Δt) and the calculation wave number interval (Δk) was also investigated. It was found that the numerical results were appropriate when Δt ≤ T(period)/1000 and Δk ≤ π/100.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.56-72
• /
• 1990

The hydrodynamic aspects of bow flare impact are discussed and the assumptions implicit in typical flare impact calculations are examined. Calculations based upon a pressure release free surface boundary condition are presented. While this simplified boundary condition eliminates the possibility of important factors such as splash-up and above-surface wetting, it has the significant advantage of being computationally simple. Both three-dimensional and two-dimensional results are compared with experiments. Errors resulting from the approximate free surface condition are examined. Also included in the two-dimensional and three-dimensional analysis are longitudinal distributions of the vertical bending moment and vertical shear force.