• Title/Summary/Keyword: a finite difference approximation

검색결과 126건 처리시간 2.623초

삼중대각행렬 시스템 풀이의 빠른 GPU 구현 (Fast GPU Implementation for the Solution of Tridiagonal Matrix Systems)

  • 김영희;이성기
    • 한국정보과학회논문지:시스템및이론
    • /
    • 제32권11_12호
    • /
    • pp.692-704
    • /
    • 2005
  • 컴퓨터 하드웨어의 급속한 발전으로 그래픽 프로세서 유닛(Graphics Processor Units : GPUs)은 굉장한 메모리 대역폭과 산술 능역을 보유하게 되어 범용 계산에 많이 활용되고 있으며, 특히 계산 집약적인 물리 기반 시뮬레이션(physics based simulation)의 GPU 구현이 활발하게 연구되고 있다. 물리 기반 시뮬레이션의 기본이 되는 미분방정식 풀이 과정에서 삼중대각행렬(tridiagonal matrix) 시스템은 유한차분(finite-difference) 근사에 의해서 자주 나타나는 선형시스템으로 물리 기반 시뮬레이션 관점에서 삼중대각행렬 시스템의 빠른 풀이는 중요한 연구 분야이다. 본 논문에서는 GPU에서 삼중대각행렬 시스템 풀이를 빠르게 구현할 수 있는 방법을 제안한다. 벡터 프로세서(vector processor) 계산에서 삼중대각행렬 시스템 풀이 방법으로 널리 사용되는 cyclic reduction 또는 odd-even reduction 알고리즘을 GPU에서 구현하였다. 본 논문에서 제안한 방법을 삼중대각행렬 시스템 풀이 방법으로 잘 알려져 있는 Thomas 방법과 GPU를 이용한 선형시스템 풀이에서 좋은 성과를 보이고 있는 conjugate gradient 방법과 비교할 때 상당한 성능 향상을 얻을 수 있었다. 또한, 열전도(heat conduction) 방정식, 이류 확산(advection-diffusion) 방정식, 얕은 물(shallow water) 방정식에 의한 물리 기반 시뮬레이션의 GPU 구현에 본 논문에서 제안한 방법을 사용하여 1024x1024 격자의 계산 영역에서 초당 35프레임 이상의 놀라운 성능을 보여주었다.

하이브리드 발견적 탐색기법을 이용한 천부 굴절법 자료의 파형역산 (Waveform inversion of shallow seismic refraction data using hybrid heuristic search method)

  • Takekoshi, Mika;Yamanaka, Hiroaki
    • 지구물리와물리탐사
    • /
    • 제12권1호
    • /
    • pp.99-104
    • /
    • 2009
  • 본 연구에서는 천부 토양층의 2차원 불균질 S파 단면을 결정하기 위해 천부 굴절법 탐사로부터 얻은 SH파 자료에 대한 파형역산 기법을 제안한다. 2차원 매질에서 SH파의 전파를 모사하기 위해 2.5차원 파동 방정식을 사용한다. 합성탄성파를 계산하기 위해 공간축으로 4차, 시간축으로 2차 근사한 staggered grid 유한차분법을 사용하여 파동 방정식을 푼다. 계산된 파형과 측정 파형의 잔여오차로 정의되는 목적함수를 하이브리드 발견적 탐색기법에 의해 최소화한다. 2차원 지하구조 모형은 각기 다른 심도 경계면을 갖는 블록과 블록 내부의 S파 속도에 의해 매개화한다. 수치실험은 암영층파 불규칙한 경계를 갖는 모델에 대해 백색잡음을 추가한 합성 SH파 자료를 이용하여 수행하였다. 지표 굴절법 자료로부터 암영층을 포함한 구조를 적절한 계산시간 내에 영상화할 수 있었다

차분 간섭도로부터 지표변위의 시계열 관측을 위한 개선된 Small Baseline Subset (SBAS) 알고리즘 (Improvement of Small Baseline Subset (SBAS) Algorithm for Measuring Time-series Surface Deformations from Differential SAR Interferograms)

  • 정형섭;이창욱;박정원;김기동;원중선
    • 대한원격탐사학회지
    • /
    • 제24권2호
    • /
    • pp.165-177
    • /
    • 2008
  • 최근 spatial decorrelation을 극복하기 위하여 비교적 짧은 기선(baseline)을 지니는 여러 시기의 차분간섭도(differential interferogram)로부터 시계열 지표변위(time-series surface deformation)를 관측할 수 있는 small baseline subset(SBAS) 알고리즘이 개발되었다. 이 SBAS 알고리즘은 singluar value decomposition(SVD)을 이용하여 시간별로 완벽하게 연결되지 않는 차분간섭도로부터 시계열 지표변위를 관측하였을 뿐 아니라 공간적으로 저주파 필터와 시간적으로 고주파 필터를 이용하여 대기효과를 보정하였다. 그러나 이 알고리즘은 초기 관측시 시계열 지표변위를 선형으로 가정하였기 때문에 각 차분간섭도의 phase unwrapping 오차를 정확하게 보정하기 어려웠을 뿐 아니라 시계열의 지표변위에 존재하는 노이즈 성분을 완화시키지 못했다. 이와 같은 단점을 보완하기 위하여 이 연구에서는 기존의 SBAS 알고리즘을 개선하였다. 이 개선된 SBAS 알고리즘은 각 차분간섭도의 phase unwrapping 오차를 최소화하기 위하여 반복적으로 시계열 지표면 변위를 개선하였고, 시계열 지표변위의 노이즈를 제거하기 위하여 유한차분근사법(finite difference approximation)을 이용하였다. 서로 다른 지역의 26개의 ERS-12자료와 21개의 RADARSAT-1 fine beam (F5) 자료를 이용하여 개선된 SBAS 알고리즘을 실험하고 분석하였다. ERS-1/2자료에서는 LOS(line-of-sight) 지표변위가 약 13년 동안 최대 -40cm가 관측되었고, RADARSAT-1 fine beam 자료에서는 약 2년 동안 최대 -3cm의 LOS 지표변위가 관측되었다.

107의 Rayleigh 수에서 원형 실린더가 존재하는 사각형 실린더 내부의 자연대류에 관한 수치적 연구 (Numerical Study of Natural Convection in a Square Enclosure with an Inner Circular Cylinder for Rayleigh Number of 107)

  • 유동훈;윤현식;하만영
    • 대한기계학회논문집B
    • /
    • 제34권8호
    • /
    • pp.739-747
    • /
    • 2010
  • 고온의 내부 원형 실린더가 존재하는 저온의 사각 밀폐계 사이의 온도차이에 의해 발생하는 자연대류에 관해 수치해석을 수행하였다. 본 연구는 내부의 원형 실린더를 표현하기 위해 유한 체적법에 기초한 가상 경계법을 사용하여 Ra = $10^7$ 에서 2 차원 비정상상태의 자연대류 현상에 대한 해를 얻었다. 더욱이 밀폐계의 수직 중심선에 따른 내부 원형 실린더의 위치 변화에 의한 열전달과 유동에 관한 영향을 연구하였다. 내부원형 실린더의 위치 변화에 따라 자연대류 현상은 비정상상태에서 정상상태로 변화 되었다. 연구 결과 두가지의 임계 위치를 얻을 수 있었다. 하부 경계는 0.05 인 지점이고 상부경계는 0.18 인 지점이다. 하부 경계와 상부 경계의 사이의 위치에서는 열 및 유동장이 정상상태임을 알 수 있었다.

지표 부근에서의 노멀전기검층 수치 모델링 (Numerical Simulation of Normal Logging Measurements in the Proximity of Earth Surface)

  • 남명진;황세호
    • 자원환경지질
    • /
    • 제43권3호
    • /
    • pp.259-267
    • /
    • 2010
  • 국내에서 노멀전기검층은 지반조사, 지하수 환경조사, 지열조사, 지질조사, 광물자원 평가 등의 다양한 분야에서 널리 이용되고 있다. 노멀전기검층은 지표 전기비저항탐사법과는 달리 완전공간에 대한 자료를 취득하기 때문에 지표부근에 수위가 위치하는 시추공에서 자료를 취득하는 경우, 수위 및 지표에 대한 영향을 고려해야 한다. 이 연구는 노멀전기검층 존데, 전류리턴전극, 기준점전위전극, 시추공내 지하수위 등을 포함하는 실제 물리검층 환경과 동일한 조건에서 노멀전기검층을 시물레이션하여 지표 및 지하수 수위가 노멀전기검층에 미치는 영향을 분석하였다. 수치 모델링은 2차원 목표지향 자기적응 고차 hp 유한요소법(2D goal-oriented high-order self-adaptive hp finite element method)을 이용하였다(여기서 h는 요소의 크기, p는 노드에서의 근사 차수). 이 알고리듬을 이용하여 측정한 겉보기비저항의 오차가 1%보다 작도록 계산할 수 있는 최적 hp 격자를 구성함으로써 매우 정밀한 결과를 얻었다. 수치실험결과, 지표부근에서 취득한 노멀전기검층 자료는 기준점전위전극이 시추공에 가까울수록 자료의 왜곡이 증가하며 장노멀전기검층에서의 왜곡이 단노멀전기검층에서 보다 심함을 알 수 있었다.

밀도 변화를 포함하는 3차원 연직함수 전개모형 (A Three-Dimensional Galerkin-FEM Model with Density Variation)

  • 이호진;정경태;소재귀;강관수;정종율
    • 한국해안해양공학회지
    • /
    • 제8권2호
    • /
    • pp.123-136
    • /
    • 1996
  • 본 연구에서는 시 공간적으로 변화하는 밀도구조를 고려한 3차원 연직함수 전개모형을 제안한다. 정수압을 가정하며 열보존식과 상태방정식의 도입으로 밀도변화가 고려된다 수평방향으로의 변화는 유한 차분 격자상에서 계산되며 유속 및 수온의 연직구조는 선형보간함수를 사용하여 계산된다. 구성된 행렬방정식은 유사변환 개념을 도입하여 시간적분된다. 개발된 모델의 테스트를 위해 간단한 이상해역과 황해역에서 대기와 해양간의 열교환에 따른 수온구조 변동을 실험하였다. 이류효과는 열수송방정식에서만 고려하였으며 연직 와동점성계수와 와동확산계수는 시 공간적으로 일정한 값을 사용하였다. 이상해역에서의 수치실험결과 모델영역의 수심의 차이에 따른 열저장의 차이로 인해 수온의 수평적 구배가 발생하였다. 결과적으로 전향력을 고려하지 않을 경우에는 상층에서는 수심이 증가하는 방향으로 흐름이 발생하고 하층에는 반대방향의 흐름이 유도된 반면 전향력을 고려할 경우에는 수온차에 의한 압력구배력과 전향력이 균형을 이루면서 지형류가 뚜렷하게 나타났다. 황해역에서는 복잡한 흐름이 나타났지만 전체적으로는 지형류의 특성이 우세하게 나타났다.

  • PDF