• 수산해양기술연구
• /
• 제53권4호
• /
• pp.309-316
• /
• 2017

In this paper, numerical modeling is conducted to analyze the tension of an anchor line by varying the size and drag coefficient of a buoy when the trapnet is influenced by the wave and the current simultaneously. A mass-spring model was used to analyze the behavior of trapnet underwater under the influence of waves and current. In the simulation of numerical model, wave height of 3, 4, 5 and 6 m, a period of 4.4 s, and the flow speed of 0.7 m/s were used for the wave and current condition. The drag coefficients of buoy were 0.8, 0.4 and 0.2, respectively. The size of buoy was 100, 50 and 25% based on the cylindrical buoy ($0.0311m^3$) used for swimming crab trap. The drag coefficient of the trapnet, the main model for numerical analysis, was obtained by a circular water channel experiment using a 6-component load cell. As a result of the simulation, the tension of the anchor line decreased proportional to buoy's drag coefficient and size; the higher the wave height, the greater the decrease rate of the tension. When the buoy drag coefficient and size decreased to one fourth, the tension of the anchor line decreased to a half and the tension of the anchor line was lower than the holding power of the anchor even at 6 m of wave height. Therefore, reducing the buoy drag coefficient and size appropriately reduces the trapnet load from the wave, which also reduces the possibility of trapnet loss.

• 한국수자원학회논문집
• /
• 제43권9호
• /
• pp.801-811
• /
• 2010

급경사에서의 고립파의 처오름을 예측하기 위해 3차원 수치모형에 내부조파기법을 도입하여 수치모형실험을 수행하였다. 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 푸는 LES(large-eddy-simulation) 기반의 수치모형을 사용한다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two-step projection 기법을 사용하여 SANS 방정식을 풀었으며, Poisson 방정식을Bi-CGSTAB 기법을 이용하여 풀고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 고립파를 3차원 공간의 일정 수심상에서 내부조파하여 해석해와 비교한 후 분산오차에 대해 분석하였다. 그리고 고립파를 내부조파하여 급경사에서의 고립파의 처오름 및 처내림 현상을 예측하고 수리모형 실험결과와 비교 및 분석하였다.

• 한국해양공학회지
• /
• 제11권2호
• /
• pp.122-130
• /
• 1997

Coastal diaster induced by waves and countermeasures were investigated in the viewpoint of reduction of overtopping rate with enviroment in fishing port. The reduction method of wave overtopping rate using ecosystem control structures was proposed and studied on the efficiency by hydraulic and numerical experiments. The estimation models on wave overtopping rate was proposed after comparing previous models with dimensional analysis and experimental results. Control function o fwave overtopping by use of ecosystem controlstructures was simulated and discussed with combining wave shoaling-dissipation-breaking deformation model around ecosystem control structures and newly proposed calculation model for wave overtopping rate. Feasiblilty of ecosystem control structures could be confirmed for reduction of wave overtopping and fisheries-based multipurpose development of coastal zone.

• 비파괴검사학회지
• /
• 제34권6호
• /
• pp.457-466
• /
• 2014

This paper presents an analytic and numerical simulation of the generation and propagation of pico-second ultrasound with nano-scale wavelength, enabling the production of bulk waves in thin films. An analytic model of laser-matter interaction and elasto-dynamic wave propagation is introduced to calculate the elastic strain pulse in microstructures. The model includes the laser-pulse absorption on the material surface, heat transfer from a photon to the elastic energy of a phonon, and acoustic wave propagation to formulate the governing equations of ultra-short ultrasound. The excitation and propagation of acoustic pulses produced by ultra-short laser pulses are numerically simulated for an aluminum substrate using the finite-difference method and compared with the analytical solution. Furthermore, Fourier analysis was performed to investigate the frequency spectrum of the simulated elastic wave pulse. It is concluded that a pico-second bulk wave with a very high frequency of up to hundreds of gigahertz is successfully generated in metals using a 100-fs laser pulse and that it can be propagated in the direction of thickness for thickness less than 100 nm.

• 한국해안·해양공학회논문집
• /
• 제31권1호
• /
• pp.17-27
• /
• 2019

Bound mode의 외중력파와 이로 인한 보강간섭이 이상파 생성 기작이라는 가설을 확인하기 위해 OpenFoam 기반 Tool box인 ihFoam과 Bi-spectrum에 기초한 수치모의를 수행하였다. 수치모의는 건설 과정에서 이상파가 관측된 삼척 LNG 생산기지 전용부두를 대상으로 수행되었으며, Bound mode의 외중력파는 Bi-spectrum에 기반하여 출현빈도가 높은 주기가 7초인 국지 풍성파와 11.4초인 너울의 difference interaction으로 생성하였다. 또한 비교를 위해 비선형 Cnoidal wave, linear wave를 대상으로 한 수치모의도 병행하여 수행하였다. 모의결과 N 계열 bound mode의 외중력파에 의해 삼척 LNG 생산기지 전용부두를 따라 진행되는 연파가 생성되며, 이상파는 전술한 연파와 남측 도류제로부터의 반사파가 더해져 출현하는 것으로 모의되었다.

• Ocean and Polar Research
• /
• 제29권1호
• /
• pp.9-31
• /
• 2007

Most seismic sea waves in the East Sea originate from earthquakes occurring near the Japanese west coast. While the waves propagate in the East Sea, they are deformed by refraction, diffraction and scattering. Though the Boussinesq equation is most applicable for such wave phenomena, it was not used in numerical modelling of seismic sea waves in the East Sea. To examine characteristics of seismic sea waves in the East Sea, numerical models based on the Boussinesq equation are established and used to simulate recent tsunamis. By considering Ursell parameter and Kajiura parameter, it is proved that Boussinesq equation is a proper equation for seismic sea waves in the East Sea. Two models based on the Boussinesq equation and linear wave equation are executed with the same initial conditions and grid size ($1min{\times}1min$), and the results are compared in various respects. The Boussinesq equation model produced better results than the linear model in respect to wave propagation and concentration of wave energy. It is also certified that the Boussinesq equation model can be used for operational purpose if it is optimized. Another Boussinesq equation model whose grid size is $40sec{\times}30sec$ is set up to simulate the 1983 and 1993 tsunamis. As the result of simulation, new propagation charts of 2 seismic sea waves focused on the Korean east coast are proposed. Even though the 1983 and 1993 tsunamis started at different areas, the propagation paths near the Korean east coast are similar and they can be distinguished into 4 paths. Among these, total energy and propagating time of the waves passing over North Korea Plateau(NKP) and South Korea Plateau(SKP) determine wave height at the Korean east coast. In case of the 1993 tsunami, the wave passing over NKP has more energy than the wave over SKP. In case of the 1983 tsunami, the huge energy of the wave passing over SKP brought about great maximum wave heights at Mukho and Imwon. The Boussinesq equation model established in this study is more useful for simulation of seismic sea waves near the Korean east coast than it is the Japanese coast. To improve understanding of seismic sea waves in shallow water, a coastal area model based on the Boussinesq equation is also required.

• 한국항해학회지
• /
• 제17권1호
• /
• pp.61-78
• /
• 1993

Large offshore structure are to be considered for oil storage facilities , marine terminals, power plants, offshore airports, industrial complexes and recreational facilities. Some of them have already been constructed. Some of the envisioned structures will be of the artificial-island type, in which the bulk of structures may act as significant barriers to normal waves and the prediction of the wave intensity will be of importance for design of structure. The present study deals wave scattering problem combining reflection and diffraction of waves due to the shape of the impermeable rigid upright structure, subject to the excitation of a plane simple harmonic wave coming from infinity. In this study, a finite difference technique for the numerical solution is applied to the boundary integral equation obtained for wave potential. The numerical solution is verified with the analytic solution. The model is applied to various structures, such as the detached breakwater (3L${\times}$0.1L), bird-type breakwater(318L${\times}$0.17L), cylinder-type and crescent -type structure (2.89L${\times}$0.6L, 0.8L${\times}$0.26L).The result are presented in wave height amplification factors and wave height diagram. Also, the amplification factors across the structure or 1 or 2 wavelengths away from the structure are compared with each given case. From the numerical simulation for the various boundary types of structure, we could figure out the transformation pattern of waves and predict the waves and predict the wave intensity in the vicinity of large artificial structures.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 제26회 KOSCO SYMPOSIUM 논문집
• /
• pp.57-62
• /
• 2003

A general procedure of obtaining reliable one-step kinetics model for hydrocarbon mixture from the fully detailed chemistry is described iin this study. One-step theoretical formulation of the induction parameter model IPM uses a theoretical reconstruction of the induction time database obtained from a detailed kinetics library. Non-dimensional induction time calculations is compared with that of detailed kinetics. The IPM was latter implemented to fluid dynamics code and applied for the numerical simulation of detonation wave propagation. The numerical results including the numerical smoked-foil record show the all the details of the detonation wave propagation characteristics at the cost around 1/100 of the detailed kinetics calculation.

• 대한기계학회논문집
• /
• 제18권3호
• /
• pp.612-623
• /
• 1994

Shock-wave focusing from a two-dimensional parabolic reflector was simulated using an explicit finite volume upwind TVD scheme. Computations were performed for three different incident shock speeds of $M_s$ = 1.1, 1.2 and 1.3, corresponding to the relatively weak, intermediate, and strong shock waves, respectively. Numerical solutions nicely resolved all the waves evolving through the focusing process. As the incident shock strength increase, a transition was observed in the shock-fronts geometry that was caused by the change in the reflection type of converging shock fronts on the axis of symmetry, from regular-type to Mach-type reflection. The computed maximum on-axis pressure amplification and the trajectories of three-wave intersections showed good agreement with experimental results. The strong nonlinear effect near the focal region which determines the shock-fronts geometries at and behind the focus and at the same time confines the pressure amplification at the focus was clearly revealed from the present numerical simulation.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2000년도 춘계 학술대회논문집
• /
• pp.51-59
• /
• 2000

The Interactions of shock wave with turbulent boundary layers in high-speed flows cause complex flowfields which result in increased adverse pressure gradients, skin friction and temperature. Accurate and reliable prediction of such phenomena is needed in designing high-speed propulsion systems. Such analyses of the complex flowfields require sophisticated numerical scheme that can resolve interactions between shock wave and boundary layers accurately. Therefore the purpose of the present. article is to introduce an accurate and efficient mixed explicit-implicit generalized Galerkin finite element method. To demonstrate the validity of the theory and numerical procedure, several benchmark cases are investigated.