• The Journal of the Acoustical Society of Korea
• /
• v.18 no.4
• /
• pp.71-77
• /
• 1999

Exact forward modeling of acoustic propagation is crucial in MFP such as inverse problems and various other acoustic applications. As acoustic propagation in shallow water environments become important, range dependent modeling has to be considered of which PE method is considered as one of the most accurate and relatively fast. In this paper higher order numerical rode employing the PE method is developed. To approximate the depth directional operator, Galerkin's method is used with partial collocation to lessen necessary calculations. Linearization of tile depth directional operator is achieved via expansion into a multiplication form of (equation omitted) approximation. To approximate the range directional equation, Crank-Nicolson's method is used. Final1y, numerical self stater is employed. Numerical tests are performed for various occan environment scenarios. The results of these tests are compared to exact solutions, OASES and RAM results.

• Water for future
• /
• v.16 no.4
• /
• pp.221-236
• /
• 1983

In semi-enclosed shallow sea areas typified by the Yellow sea and the East China Sea, currents and sea surface variations are predominantly tidal. During the recent years two-dimensional numerical hydrodynamic model of the Yellow Sea and the East China Sea has been developed, based on the vertically-integrated equations of motion and continuity, capable of reproducing amplitudes and phases of the principal components of tides to satisfiable accuracy. As a subsequent development a three-dimensional hydrodynamical nymerical model covering the Yellow Sea and the East China Sea has been formulated to investigate the vertical distribution of horizontal tidal current and the response of the continented to investigate the vertical distribution of horizontal tidal current and the response of the continental shelf sea to steady uniform wind stress field imposed over the surface. Features of the M2 tidal current and the wind-induced three-dimensional current structure determined from the computation have been examined and discussed.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.177-184
• /
• 1997

In this paper, numerical methods are developed for solving the elastica of simple beams with variable-arc-length subjected to a point loading. The beam model is based on Bernoulli-Euler beam theory. The Runge-Kutta and Regula-Falsi methods, respectively, are used to solve the governing differential equations and to compute the beam's rotation at the left end of the beams. Extensive numerical results of the elastica responses, including deflected shapes, rotations of cross-section and bending moments, are presented in non-dimensional forms. The possible maximum values of the end rotation, deflection and bending moment are determined by analyzing the numerical data obtained in this study.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.2
• /
• pp.75-84
• /
• 2002

A numerical analysis of unsteady motion in solid rocket motors with a nozzle has been conducted. The numerical formulation including modified $\kappa$-$\varepsilon$ turbulence model treats the complete conservation equation for the gas phase and the one-dimensional equations in the radial direction for the condensed phase. A fully coupled implicit scheme based on a dual time-stepping integration algorithm has been adopted to solve the governing equations. After obtaining a steady state solution, pulse and periodic oscillations of pressure are imposed at the head-end to simulate acoustic oscillations of a travelling-wave motion in the combustion chamber. Various steady and unsteady state features in the combustion chamber of a rocket motor has been analyzed as results of numerical calculations.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.7 no.2
• /
• pp.82-88
• /
• 2004

To evaluate the tidal flow around the transmission type harbour structures, the numerical models considering hydraulic resistance are used. The traditional governing equations of depth-integratef hydrodynamics are modified by using of the drag force term. As a numerical model to evaluate the rate of tidal exchange, the Random-walk method is used. The models are applied in a design of the pile-supported pier structures in Busan harbour site, Korea, where the flow speed and the tidal exchange are more promotive than the case of non-transmission structure. The developed model will be applied usefully to design the transmission type structure in the ocean space.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.237-237
• /
• 2022

우리나라 하천 내에 설치된 횡단구조물은 보 35,000개소 이상으로 전국 하천의 0.6 km마다 수공구조물이 위치해 있으며 이로 인해 종적(longitudinal) 연결성에 영향을 미치는 것으로 나타났다. 농어촌연구원의 조사에 따르면 하천 횡단구조물이 설치된 구간에 어도 설치율은 14.9%로 나타났으며 이는 수생태 관점에서 연속성이 매우 열악한 상태임을 확인 할 수 있다. 이에 우리나라 각 부처 국토교통부, 환경부 등에서는 수생태 연속성을 확보하고자하는 노력이 지속되고 있으며 대표적으로 어류의 이동통로인 어도를 설치하거나 기존 어도의 효율을 향상시키기 위하여 개보수 작업을 지속적으로 실시하고 있다. 어도의 기능을 평가하기 위해서는 어도 내의 수리특성을 정확히 파악하는 것이 중요한데 연속적인 구조물로 구성된 어도 내 흐름은 매우 복잡하다. 특히 어도 내 구조물간 상호작용에 의하여 비정상 흐름이 발생하며 구조물 뒤에서는 사수역(dead zone)이 형성된다. 사수역에서 나타나는 와류의 거동은 구조물의 기하학적 특성에 따라 변화한다. 본 연구에서는 2차원 수심적분 수치모형을 활용하였으며 벽면 근처 점성저층의 유속분포을 재현하기 위하여와 점성항에 감쇠함수(damping function)를 고려하였다. 수치모형의 검증을 위해 실내 실험수로의 직선 개수로에서 PIV(particle image velocimetry)를 활용하여 연속적으로 배치된 구조물에 의한 유속자료를 활용하였다. 이 결과는 향후 새로운 어도설계 혹은 기 설치된 어도의 수리학적 기능을 평가하는데 활용이 가능할 것으로 보인다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.3
• /
• pp.265-276
• /
• 1995

In this Paper, the Performance of the convolution method has been investigated as an effort to develop a simple system of predicting wind-driven surface current on a real time basis. In this approach wind stress is assumed to be spatially uniform and the effect of atmospheric pressure is neglected. The discrete convolution weights are determined in advance at each point using a linear three-dimensional Galerkin model with linear shape functions(Galerkin-FEM model). Four directions of wind stress(e.g. NE, SW, NW, SE) with unit magnitude are imposed in the model calculation for the construction of data base for convolution weights. Given the time history of wind stress, it is then possible to predict with-driven currents promptly using the convolution product of finite length. An unsteady wind stress of arbitrary form can be approximated by a series of wind pulses with magnitude of 6 hour averaged value. A total of 12 pulses are involved in the convolution product To examine the accuracy of the convolution method a series of numerical experiments has been carried out in the idealized basin representing the scale of the Yellow Sea and the East China Sea. The wind stress imposed varies sinusoidally in time. It was found that the predicted surface currents and elevation fields were in good agreement with the results computed by the direct integration of the Galerkin model. A model with grid 1/8$^{\circ}$ in latitude, l/6$^{\circ}$ in longitude was established which covers the entire region of the Yellow Sea and the East China Sea. The numerical prediction in terms of the convolution product has been carried out with particular attention on the formation of upwind flow in the middle of the Yellow Sea by northerly wind.

• Proceedings of the Acoustical Society of Korea Conference
• /
• 1992.06a
• /
• pp.91-95
• /
• 1992

개영역(open boundary) 내의 변환기에 대한 음장 특성은 경계 조건이 좌표계와 정합되는 특수한 경우 헤름홀쯔 방정식풀이로 그렇지 않은 복잡한 형태의 음장에 대해서는 경계적분법이나 근사해법으로 얻어진다. 그러나 최근 이러한 문제들은 발달된 컴퓨터의 계산능력으로 유한요소법이나 경계요소법을 이용한 수치해법으로 해결되고 있다. 본 연구에서는 경계조건 변화에 따른 원통형 압전변환기의 방사임피던스 및 지향특성을 유한요소법 및 하이브리드형 무한요소법을 이용하여 구하고 기존의 해석적 결과와 비교하여 수중용 변환기 제작의 설계에 대한 기초자료를 제공하고자 한다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.4
• /
• pp.279-287
• /
• 1993

Three dimensional linear transient wave propagation over a constant depth is presented using Green's function method. Present solution is proved to be a more general solution from which the existing solutions are obtained by using the appropriate assumptions. The effect of dispersion on wave height attenuation is shown and discussed on the numerical results of the solution.

• Proceedings of the Korea Air Pollution Research Association Conference
• /
• 1999.10a
• /
• pp.194-196
• /
• 1999

대기나 수용액 속에 부유 입자는 서로 충돌하여 합쳐져서 그 크기가 커지게 된다. 이러한 과정을 응집(Coagulation)이라고 하며, 이는 대기중 부유입자의 농도 및 크기분포의 변화, 구름 속에서의 빗방울형성 등에 매우 중요한 기작 중의 하나이다. 응집방정식은 일반적으로 비선형 편미적분 방정식으로 표현되어 일반 해를 구하는 것은 불가능하다. 이러한 이유로 응집방정식을 풀 때에는 수치 해석적인 방법이 주로 이용되고 있다.(Tolof, 1977; Gelbard and Seinfeld, 1978; Reed ea al., 1980; Mick et al., 1991).(중략)