• Ocean Systems Engineering
• /
• 제7권4호
• /
• pp.399-412
• /
• 2017

The Gauss-Legendre integral method is applied to numerically evaluate the Green function and its derivatives in finite water depth. In this method, the singular point of the function in the traditional integral equation can be avoided. Moreover, based on the improved Gauss-Laguerre integral method proposed in the previous research, a new methodology is developed through the Gauss-Legendre integral. Using this new methodology, the Green function with the field and source points near the water surface can be obtained, which is less mentioned in the previous research. The accuracy and efficiency of this new method is investigated. The numerical results using a Gauss-Legendre integral method show good agreements with other numerical results of direct calculations and series form in the far field. Furthermore, the cases with the field and source points near the water surface are also considered. Considering the computational efficiency, the method using the Gauss-Legendre integral proposed in this paper could obtain the accurate numerical results of the Green function and its derivatives in finite water depth and can be adopted in the near field.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제11권1호
• /
• pp.37-44
• /
• 2006

본 연구는 오염물 거동에 대한 수치해석을 위해 보편적으로 사용되고 있는 수치 방법들의 장단점을 총괄적으로 나타내고, 효율적인 수치모델링 기법 개발을 위해 ELLAM과 LEZOOMPC를 비교분석하였다. 지하수 분야에서 가장 많이 사용되는 수치 방법은 Eulerian-Lagrangian 방식과 Eulerian 방식인데, Eulerian-Lagrangian 방식은 수치영역 내에서 일반적으로 질량을 보존하지 못하고, 경계조건을 체계적으로 처리하지 못하는 한계를 갖고 있다. 반면에 Eulerian 빙식은 시간 및 공간 절삭 오차로 인해서 시간 간격 및 격자 크기를 극히 줄여야 하는 제약을 갖고 있다. 최근 10 년간 지하수 분야에서 크게 대두되고 있는 수치기법인 ELLAM(Eulerian Lagrangian Localized Adjoint Method)은 Eulerian-Lagrangian 방식과 Eulerian 방식에서 나타나는 수치 제약점이나 한계점을 동시에 해결하는 수치기법으로 알려져 왔다. 그러나 본 연구에서는 ELLAM의 장단점을 파악하고 보완점을 제안한다. ELLAM의 단점을 파악하기 위해, mesh Peclet number가 다른 예제들을 설정하고, 그 예제들에 대한 ELLAM, LEZOOMPC(Lagrangian-Eulerian ZOOMing Peak and valley Capturing)와 visual MODFLOW의 수치결과들을 해석해와 비교하였다. Mesh Peclet number가 무한대일 때 ELLAM의 수치결과는 수치진동으로 인해 해석해와 일치하지 않았으나, LEZOOMPC의 수치 결과는 해석해와 일치했다. 위의 결과는 ELLAM의 수치오차가 LEZOOMPC의 특성을 이용하여 개선 및 보완될 수 있는 가능성을 시사해 준다. 따라서 ELLAM에 LEZOOMPC의 후향 입지추적, 전향 입지추적, 선택적 국부 격자 세립화 과정과 최고/최저 농도점 이동 추적 과정을 결합하면 ELLAM의 수치적 장점을 유지하면서 mesh Peclet number에 제약을 받지 않는 효율적인 수치모델링 기법을 개발할 수 있을 것으로 판단된다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제22권1호
• /
• pp.75-89
• /
• 2018

This paper presents the numerical valuation of the arithmetic Asian option by using the operator-splitting method (OSM). Since there is no closed-form solution for the arithmetic Asian option, finding a good numerical algorithm to value the arithmetic Asian option is important. In this paper, we focus on a two-dimensional PDE. The OSM is famous for dealing with plural-dimensional PDE using finite difference discretization. We provide a detailed numerical algorithm and compare results with MCS method to show the performance of the method.

• 한국유체기계학회 논문집
• /
• 제1권1호
• /
• pp.72-80
• /
• 1998

Computational analysis of incompressible flows by numerically solving Navier-Stokes equations using multi-block finite volume method is conducted on a parallel computing system. Numerical algorithms adopted in this study $include^{(1)}$ QUICK upwinding scheme for convective $terms,^{(2)}$ central differencing for other terms $and^{(3)}$ the second-order Euler differencing for time-marching procedure. Structured grids are used on the body-fitted coordinate with multi-block concept which uses overlaid grids on the block-interfacing boundaries. Computational code is parallelized on the MPI environment. Numerical accuracy of the computational method is verified by solving a benchmark test case of the flow inside two-dimensional rectangular cavity. Computation in the axial compressor cascade is conducted by using 4 PE's md, as results, no numerical instabilities are observed and it is expected that the present computational method can be applied to the turbomachinery flow problems without major difficulties.

• 한국해안·해양공학회논문집
• /
• 제20권5호
• /
• pp.461-471
• /
• 2008

WKB 방법에 의한 일차원 완경사 파랑식의 해석해를 유도하였고 이는 Porter(2003)의 해와 유사한 형태를 갖는다. 적용적 측면에서 본 해석해는 관련 수치해에 견줄 만큼 해석상의 일반성을 갖는다. 유도과정에서 해면함수로 표현된 굴절 방정식의 해도 얻었다. Bremmer 방법을 이용한 본 해석해에 대한 수치계산 결과를 제시하였고 이들은 기존 결과와 일치한다.

• 한국군사과학기술학회지
• /
• 제11권6호
• /
• pp.164-171
• /
• 2008

A modified generalized particle algorithm, MGPA, was suggested to improve the computational efficiency of standard SPH method in numerical analysis of high speed impact behavior. This method uses a numerical failure mechanism than material failure models to describe the target penetration. MGPA algorithm was more effective to describe the impact phenomena and new boundaries produced during the calculation process were well recognized and treated in the target penetration problem of a bullet. When bullet perforation problems were analyzed by this method, MGPA algorithm calculation gives the stable numerical solution and stress oscillation or particle penetration phenomena were not shown. The error range in ballistic velocity limit is less than $2{\sim}13%$ for various target thickness.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2010년도 춘계 학술발표회
• /
• pp.852-855
• /
• 2010

Numerical simulations for large scale triaxial tests with large diameter rockfill materials are conducted using distinct element method. For generation of compacted assembly with specific grain size distribution and initial material porosity, the clump logic method and expansion of generated particles are adapted. With micro parameters which are chosen by calibration process, discrete particle modelling of triaxial test in case of other confining stress and cyclic loading condition were conducted. Also numerical simulations of fluid injection into particulate materials were conducted to observe cavity initiation and propagation using distinct element method. The fluid scheme solves the continuity and Navior-Stokes equations numerically, then derives pressure and velocity vectors for fixed grid by considering the existence of particles within the fluid cell.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2000년도 춘계학술대회 논문집
• /
• pp.140-145
• /
• 2000

The spectral method which is composed of an eigenfunction expansion of free modes in the wave number domain is used to produce two dimensional unsteady inviscid wave simulation such as progressive waves in a numerical pneumatic wave tank. A spatial and time dependent free surface elevation and the potential are calculated by integrating ODE derived from fully nonlinear kinematic and dynamic free surface boundary condition at each time step. The nonlinear characteristics in the waves by this method were notable as increasing wave steepness. This method is very useful and powerful in terms of saving computational time caused by rapid convergence exponentially with increasing number of nodes, even preserving accurate numerical results. Moreover, it will given us many possibilities to apply to naval and ocean engineering fields.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1997년도 추계학술대회 논문집
• /
• pp.799-803
• /
• 1997

Smell and sound are closely connected with our life. But there is no way to represent the exact numerical value of smell. On the contary, there is a method of represent the exact numerical value of sound. In this thesis, using and adapting this method about sound, a new method of expressing the numercal value of smell is going to be derived. These two methods are similar, but there are so many substances which emit smells and the methods of measuring the density of the substances are various according to kinds of the substances. So the new method about smell will be derived by a new idea.

• Journal of applied mathematics & informatics
• /
• 제28권3_4호
• /
• pp.611-624
• /
• 2010

A Numerical scheme based on collocation method using quartic B-spline functions is designed for the numerical solution of one-dimensional modified equal width wave (MEW) wave equation. Using Von-Neumann approach the scheme is shown to be unconditionally stable. Performance of the method is validated through test problems including single wave, interaction of two waves and use of Maxwellian initial condition. Using error norms $L_2$ and $L_{\infty}$ and conservative properties of mass, momentum and energy, accuracy and efficiency of the suggested method is established through comparison with the existing numerical techniques.