• Journal of Korean Port Research
• /
• v.6 no.2
• /
• pp.33-42
• /
• 1992

In this paper we investigate the motion response of a moored ship in the fluid region sheltered by inclined breakwaters. The matched asymptotic expansion technique is employed to analyze the wave fields scattered by the inclined breakwaters. Fluid domain is subdivided into the ocean, entrance and sheltered regions. Unknown coefficients contained in each region can be determined by matching at the intermediate zone between two neighboring regions. The wave field generated by the ship motion can be analyzed in terms of Green's function method. To obtain the velocity jump across the ship associated with the symmetric motion modes, the sheltered region is further divided into near field of the ship and the rest field. The image method is introduced to consider the effect of the pier near the ship. The integral equation for the velocity jump is derived by the flux matching between the inner region and the outer region of a moored ship. Throughout the numerical calculation it is found that the inclined angle width of entrance of breakwaters as well as the location of moored vessel play an important role in the motion response of a moored ship.

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.6
• /
• pp.31-37
• /
• 1999

This paper describes the far-field estimation using the near-field measurement data. Measurement in far-field region gives us the acoustical characteristics of the source but in general measurement is made in near-field such as acoustic water tank or anechoic chamber, so far-field acoustical characteristics of the source should be predicted from near-field data. In this case, the number of measurement points in the near field which relates to the accuracy of the predicted field and the amount of data processing, should be optimized. Existing papers say that measurement points is proportional to kL and depends on geometry and directivity of the source. But they do not give us any definite criterion for the required number of measurement points. Boundary Collocation Method which is one of the far-field prediction methods, is analyzed based on Helmholtz integral equation and Green function and it has been found that the number of measurement points is optimized as 0.54kL which is about one half of the existing results.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.4
• /
• pp.25-34
• /
• 1986

A two-dimensional linear method has been developed for the motion and the second-order steady force arising from the hydrodynamic coupling between waves and currents in the presence of a body of arbitrary shape. Interaction between the incident wave and current in the absence of the body lies in the realm beyond our interest. A Fredholm integral equation of the second kind is employed in association with the Haskind's potential for a steadily moving source of pulsating strength located in or below the free surface. The numerical calculations at the preliminary stage showed a significant fluctuation of the hydrodynamic forces on the surface-piercing body. The problem is approximately solved by using the asymptotic Green function for $U^2{\rightarrow}0$. The original Green function, however, is applied for the fully submerged body. Numerical calculations are made for a submerged and for a half-immersed circular cylinder and extensively for the mid-ship section of a Lewis-form. Some of the results are compared with other analytical results without any available experimental data. The current has strong influence on roll motion near resonance. When the current opposes the waves, the roll response are generally negligible in the low frequency region. The current has strong influence on roll motion near resonance. When the current opposes the wave, the roll response decreases. When the current and wave come from the same direction, the roll response increases significantly, as the current speed increases. The mean drift forces and moment on the submerged body are more affected by current than those on the semi-immersed circular cylinder or on the ship-like section in the encounter frequency domain.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.12-18
• /
• 1991

Three-dimensional nonlinear sloshing effects due to tank motions are simulated by solving boundary value problem using the panel method based on boundary integral technique. While Shinkai used boundary elements on which source strengths vary linearly between nodes, the source of constant strength is distributed on each triangular panel in the present study. The source strength at each time step is determined by solving the Fredholm integral equation of the second kind obtained from Green's theorem. To avoid cumulative numerical errors as time elapses, Adam-Bashforth-Moulton method is employed. Numerical examples for the case of partially filled spherical tank on board oscillating in harmonic sway mode or pitch mode are solved. The elevation of the free surface is compared with the result by Shinkai and confirmed in good agreement during early time. The input and the output energy are comparatively evaluated to check the overall accuracy of the present numerical scheme. Although some leakage of energy are found as time marches, it is plausible when we take into account nonlinearities of the problem and the number of panels of the model.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.4
• /
• pp.15-26
• /
• 1990

This paper describes a potential-baoed panel method formulated for the analysis cf a supercavitating two-dimensional symmetri strut. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring that the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type, With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lifting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.3
• /
• pp.35-47
• /
• 1996

A surface panel method treating a boundary-value problem of the Dirichlet type is presented to design a three dimensional body with free surface corresponding to a prescribed pressure distribution. An integral equation is derived from Green's theorem, giving a relation between total potential of known strength and the unknown local flux. Upon discretization, a system of linear simultaneous equations is formed including free surface boundary condition and is solved for an assumed geometry. The pseudo local flux, present due to the incorrect positioning of the assumed geometry, plays a role f the geometry corrector, with which the new geometry is computed for the next iteration. Sample designs for submerged spheroids and Wigley hull and carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. For the calculation of the wave resistance, normal dipoles and Rankine sources are distributed on the body surface and Rankine sources on the free surface. The free surface boundary condition is linearized with respect to the oncoming flow. Four-points upwind finite difference scheme is used to compute the free surface boundary condition. A hyperboloidal panel is adopted to represent the hull surface, which can compensate the defects of the low-order panel method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but by experiment.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.159-173
• /
• 1991

This paper describes a potential-based panel method formulated for the analysis of a super-cavitating two-dimensional hydrofoil. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring thats the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lilting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged. Characteristics of iteration and discretization of the present numerical method are much faster and more stable than the existing nonlinear theories. The theory shows good correlations with the existing theories and experimental results for the super-cavitating flow. In the region of small angles of attack, the present prediction shows and excellent comparison with the Geurst's linear theory. For the long cavity, the method recovers the trends of the Wu's nonlinear theory. In the intermediate regions of the short super-cavitation, the method compares very well with the experimental results of Parkin and also those of Silberman.