• Journal of Ocean Engineering and Technology
• /
• v.23 no.6
• /
• pp.12-16
• /
• 2009

In this study, the wave profiles and kinematics of highly nonlinear waves at various water depths were calculated using a 2D fully nonlinear Numerical Wave Tank (NWT). The NWT was developed based on the Boundary Element Method (BEM) with the potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme by 4th-order Runge-Kutta time integration. The spatial variation of intermediate-depth waves along the direction of wave propagation was caused by the unintended generation of 2nd-order free waves, which were originally investigated both theoretically and experimentally by Goda (1998). These free waves were induced by the mismatch between the linear motion of wave maker and nonlinear displacement of water particles adjacent to the maker. When the 2nd-order wave maker motion was applied, the spatial modulation of the waves caused by the free waves was not observed. The respective magnitudes of the nonlinear wave components for various water depths were compared. It was found that the high-order wave components greatly increase as the water depth decreases. The wave kinematics at various locations were calculated and compared with the linear and the Stokes 2nd-order theories.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.3
• /
• pp.258-268
• /
• 2008

The radiation potential of a ship advancing in waves is studied using the 3D time-domain forward-speed free-surface Green function and the Green integral equation. Numerical solutions are obtained by making use of the 2nd order BEM(Boundary Element Method) which make it possible to take account of the line integral along the waterline in a rigorous manner. The 6 degree of freedom motion memory functions of a hemisphere and the Wigley seakeeping model obtained by direct integration of the time-domain 3D potentials over the wetted surface are presented for various Froude numbers.

• 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 the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.291-305
• /
• 2010

A three-dimensional time-domain calculation method is of crucial importance in prediction of the motions and wave loads of a ship advancing in a severe irregular sea. The exact solution of the free surface wave-ship interaction problem is very complicated because of the essentially nonlinear boundary conditions. In this paper, an approximate body nonlinear approach based on the three-dimensional time-domain forward-speed free-surface Green function has been presented. The Froude-Krylov force and the hydrostatic restoring force are calculated over the instantaneous wetted surface of the ship while the forces due to the radiation and scattering potentials over the mean wetted surface. The time-domain radiation and scattering potentials have been obtained from a time invariant kernel of integral equations for the potentials which are discretized according to the second-order boundary element method (Hong and Hong 2008). The diffraction impulse-response functions of the Wigley seakeeping model advancing in transient head waves at various Froude numbers have been presented. A simulation of coupled heave-pitch motion of a long rectangular barge advancing in regular head waves of large amplitude has been carried out. Comparisons between the linear and the approximate body nonlinear numerical results of motions and wave loads of the barge at a nonzero Froude number have been made.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.6
• /
• pp.16-22
• /
• 2010

In this paper, the diffraction problems for fixed offshore structures are solved using a hybrid scheme. In this hybrid scheme, potential-based solutions and the Navier-Stokes-based finite volume method (FVM) with a volume-of-fluid (VOF) method are combined. We introduce a buffer zone for efficient wave-making and damping. In this buffer zone, the near field solution from FVM-VOF is gradually changed to Stokes' 2nd order wave solutions. Three different models, including the truncated cylinder, sphere, and wigleyIII model, are numerically investigated in regular waves with a wave steepness of 1/30. The efficiency and accuracy of the hybrid scheme are numerically validated from results using different domain sizes and buffer zones. The wave exciting forces from the FVM-VOF simulations are compared with experiments and potential-based solutions from the higher-order boundary element method (HOBEM). This comparison shows good agreement between the hybrid scheme and potential-based solutions.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.34 no.1
• /
• pp.43-50
• /
• 2001

In order to analyze the dynamic motions of fish cage systems made of a frame and a netting under the conditions of waves and current, the hydraulic model experiment at towing tank and the numerical computation using boundary integral element method based on linear potential theory were carried out on a square and a circular type of fish cage, The computed and measured results for the dynamic motions of model fish cage systems showed that the heave and pitch motions were almost unaffected by the inclusion of nets, while the surge motions were very reduced by drag force acting on them. In addition, irregular wave-induced motions of fish cages included non-negligible 2nd order harmonic components at high frequency nearly twice the wave frequency. The reason why these motions were considered was due to resonance or structural components of frames being overflown and out of water during a wave cycle. It was found that circular type was more desirable structure in the open sea than square one only in the respect of dynamic motions due to waves and current. Further verifications were needed considering hydrodynamic forces, fatigue life, and structure analysis based on long term stochastic waves including frequency and time domain for the purpose of analyzing and designing fish cage systems.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.8 no.4
• /
• pp.110-120
• /
• 2000

In order to improve the vibration characteristics of mid sized passenger car automatic transmission at idle experimental and theoretical studies have been carried out. Idle shake in "D" range occurs by various reasons such as characteristics of body bending resonance between subsystems and engine mounts etc. Using full vehicle finite element analyses and modal tests we introduce the way to reduce the idle shake in the early design stage. It shows that the exciting forces are the 2nd order torque and force of engine. A powertrain system modes in "D" range are entirely effected by the additional boundary conditions of drive line. As a result the frequencies of subsystems are arranged to be lined up at the idle frequency range in order to avoid the resonances with subsystems To reduce the idle shake mounts of radiator are tuned to act as a dynamic damper to 1st bending frequency of the body. In addition a hydraulic mount which is optimized by Phase Shift Method is applied to the rear engine mount.e rear engine mount.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.3 s.46
• /
• pp.291-302
• /
• 2000

In recent years the development of high modulus, high strength and low density boron and graphite fibers bonded together has brought renewed interestes in structural elements. When a plate with arbitrarily oriented layers and clamped boundary conditions is subjected to uniform loading, it is difficult to analyze and apply, compared with isotropic and orthotropic cases. Therefore the numerical methods, such as finite difference method or finite element method, should be emloyed to analyse such problems. In this study the finite difference technique is used to formulate the bending analysis of symmetric composite laminated skew plates. When this technique is used to solve the problem, it is desirable to reduce the order of the derivatives in order to minimize the number of the pivotal points involved in each equation. The 4th order partial differential equations of laminated skew plates are converted to an equivalent three of 2nd order partial differential equations with three dependant variables.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.12
• /
• pp.889-898
• /
• 2022

In this paper, parallel computing method on the three-dimensional electromagnetic field is proposed. The present electromagnetic scattering analysis is conducted based on the time-harmonic vector wave equation and the finite element method. The edge-based element and 2^nd -order absorbing boundary condition are used. Parallelization of the elemental numerical integration and the matrix assemblage is accomplished by allocating the partitioned finite element subdomain for each processor. The graph partitioning library, METIS, is employed for the subdomain generation. The large sparse matrix computation is conducted by MUMPS, which is the parallel computing library based on the multi-frontal method. The accuracy of the present program is validated by the comparison against the Mie-series analytical solution and the results by ANSYS HFSS. In addition, the scalability is verified by measuring the speed-up in terms of the number of processors used. The present electromagnetic scattering analysis is performed for a perfect electric conductor sphere, isotropic/anisotropic dielectric sphere, and the missile configuration. The algorithm of the present program will be applied to the finite element and tearing method, aiming for the further extended parallel computing performance.