• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2008년도 공동학술대회
• /
• pp.57-62
• /
• 2008

탄성파 수치 모형 계산에 있어서 널리 사용되는 엇갈린 격자 방법이 아니라 회전된 엇갈린 격자 방법을 사용하여 탄성파 수치 모사를 수행하였다. 표준 엇갈린 격자 방법에서는 자유 경계조건을 적용하여야 하는 단점이 있지만 회전된 격자 방법에서는 자유 경계조건을 적용할 필요 없이 물성을 부여함으로써 해결될 수 있다. 파동전파에 있어서 유한 경계 조건에서 발생하는 인공 반사파를 제거하기 위하여 최근에 뛰어난 성능을 보이는 PML(Perfectly Matched Layer)법이 많이 적용되고 있다. 이 연구에서는 PML보다 성능이 뛰어나기도 하지만 차분에 있어 파동식을 분리할 필요가 없는 CPML(Convolutional Perfectly Matched Layer)법을 회전된 엇갈린 격자법(RSG : Rotatged Staggered Grid)에 적용하고자 하였다. Cerjan법의 감쇠법과 비교한 결과 매우 효과적으로 인공 반사파를 제거할 수 있음을 알 수 있었다.

• 한국전산유체공학회지
• /
• 제7권1호
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• pp.1821-1826
• /
• 2003

In this study, we discuss CIP method, which can treat compressible/incompressible problem and multi-phase problem. We can apply this method to the general equations by using CCUP. In this paper, non-staggered grid arrangement and predictor-corrector method are used to enhance later development and the solution accuracy and convergence performance. To validate the numerical algorithm proposed in this paper, the two-dimensional unsteady flow in the driven cavity is simulated. The numerical results of this subject using the present algorithm are compared with other numerical results. As a result, it is prived that the present scheme gives stable and improved solutions, and the results even coarse grid are in good agreement with other result using a fine grid arrangement.

• 지구물리와물리탐사
• /
• 제12권2호
• /
• pp.183-191
• /
• 2009

탄성파 수치 모형 계산에 있어서 널리 사용되는 엇갈린 격자 방법이 아니라 회전된 엇갈린 격자 방법을 사용하여 탄성파 수치 모사를 수행하였다. 표준 엇갈린 격자 방법에서는 특별한 자유 경계조건을 적용하여야 하는 단점이 있지만 회전된 격자 방법에서는 물성으로 진공 또는 공기층을 부여함으로써 자유 경계조건을 실현가능하다는 것을 확인할 수 있었다. 파동전파에 있어서 유한 경계 조건에서 발생하는 인공 반사파를 제거하기 위해 PML (Perfectly Matched Layer)의 파동식 분해라는 단점을 극복할 수 있고 좋은 성능을 보이는 CPML (Convolutional Perfectly Matched Layer)법을 회전된 엇갈린 격자법(RSG: Rotatged Staggered Grid)에 적용하였다. 회전된 격자 유한 차분법에서 CPML의 고주파수 흡수 특성과 에너지 흡수율 조사, Cerjan법의 감쇠를 비교한 결과 흡수경계조건으로 좋은 성능을 확인하였다. 유체와 고체의 모형에 대한 경계에 대하여서도 매우 효과적으로 경계면에서 발생하는 반사파를 제거할 수 있음을 알 수 있었다.

• Journal of applied mathematics & informatics
• /
• 제22권1_2호
• /
• pp.67-82
• /
• 2006

In this paper, we present a finite difference method for the implementation of the rotation of a circular cylinder in the incompressible flow field by solving the two-dimensional unsteady Navier-Stokes equations. The approach is to use staggered grid method so that the accuracy and order of convergence of the associated algorithms can be maintained. The proposed method is easy to be implemented and is effective. A set of simulations for the flow dynamics is provided to show the computational results.

• 한국전산유체공학회지
• /
• 제12권1호
• /
• pp.1-8
• /
• 2007

A code is developed using the hybrid Cartesian/immersed boundary method and it is applied to simulate flows around a three-dimensional deforming body. A new criterion is suggested to distribute the immersed boundary nodes based on edges crossing a body boundary. Velocities are reconstructed at the immersed boundary nodes using the interpolation along a local normal line to the boundary. Reconstruction of the pressure at the immersed boundary node is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other experimental and numerical results for the velocity profiles around a circular cylinder under the forced in-line oscillation and the pressure coefficient distribution on a sphere. The code is applied to simulate the flow fields around a plate whose tail is periodically flapping under a translation. The effects of the velocity and acceleration due to the deformation on the periodic shedding of pairs of tip vortices are investigated.

• 한국전산유체공학회지
• /
• 제12권2호
• /
• pp.37-45
• /
• 2007

The hybrid Cartesian/immersed boundary method is applied to simulate fluid-structure interaction of a two-dimensional orbiting flexible foil. The elastic deformation of the flexible foil is modelled based on the dynamic equation of a thin-plate. At each time step, the locations and velocities of the Lagrangian control points on the flexible foil are used to reconstruct the boundary conditions for the flow solver based on the hybrid staggered/non-staggered grid. To test the developed code, the flow fields around a flapping elliptical wing are calculated. The time history of the vertical force component and the evolution of the vorticity fields are compared with recent other computations and good agreement is achieved. For the orbiting flexible foil, the vorticity fields are compared with those of the case without the deformation. The combined effects of the angle of attack and the orbit on the deformation are investigated. The grid independency study is carried out for the computed time history of the deformation at the tip.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.96-101
• /
• 1998

Various discretiation methods of Laplacian operator in the Pressure-Poisson equation are investigated for the solution of incompressible Navier-Stokes equations using the non-staggered grid. Laplacian operators previously proposed by other researchers are applied to a Driven-Cavity problem. The computational results are compared with those of Ghia. The results show the characteristics of the discrete Laplacian operators.

• 한국지구과학회지
• /
• 제39권4호
• /
• pp.317-326
• /
• 2018

Effect of the nonuniform grid on the two-dimensional transport equation was investigated in terms of theoretical analysis and finite difference method (FDM). The nonuniform grid having a typical structure of the numerical weather forecast model was incorporated in the vertical direction, while the uniform grid was used in the zonal direction. The staggered and non-staggered grid were placed in the vertical and zonal direction, respectively. Time stepping was performed with the third-order Runge Kutta scheme. An error analysis of the spatial discretization on the nonuniform grid was carried out, which indicated that the combined effect of the nonuniform grid and advection velocity produced either numerical diffusion or numerical adverse-diffusion. An analytic function is used for the quantitative evaluation of the errors associated with the discretized transport equation. Numerical experiments with the non-uniformity of vertical grid were found to support the analysis.

• 지구물리와물리탐사
• /
• 제20권3호
• /
• pp.121-128
• /
• 2017

국내에서 시간영역 전자탐사(time-domain/transient electromagnetic, TEM) 자료의 해석은 1차원 알고리듬에 주로 의존하고 있는 실정이기 때문에 정밀한 해석을 위하여 3차원 모델링 및 역산 해석 프로그램의 개발이 필요한 상황이다. 이 연구에서는 엇갈린 격자를 이용한 시간영역 유한차분(staggered-grid finite-difference time-domain, FDTD)법에 기초하여 3차원 TEM 반응 모델링 알고리듬을 개발하였다. 시간영역 전자탐사의 모델링을 위해 맥스웰 방정식을 현시적 중앙점 FDTD법을 이용하여 이산화하였으며 수치 안정성을 높이기 위해 가상 변위전류항을 도입하였다. 일반적으로 많이 활용되는 소형 코일 송신원을 수치적으로 구현하여 균질 반무한 공간에서의 해석해와 비교 검증하고 3차원 이상체에 대한 반응을 분석하였다. 이 연구에서 개발된 모델링 프로그램은 향후 TEM 전자탐사 자료의 정밀 해석에 기초가 될 것으로 기대한다.