• 제목/요약/키워드: 유한차분식 모형

검색결과 58건 처리시간 0.025초

지진해일 전파모의를 위한 유한차분모형 개발 (Development of Finite Difference Model for Tsunami Propagation)

  • 안성호;하태민;조용식
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2010년도 학술발표회
    • /
    • pp.52-56
    • /
    • 2010
  • 본 연구에서는 지진해일의 전파 과정을 모의함에 있어 선형 천수방정식의 수치분산을 이용하는 기법이 아닌 선형 Boussinesq 방정식을 직접 차분하는 유한차분기법을 제안하였다. 지배방적식과 차분식의 일치성을 해석하기 위해 이산화 오차를 확인하고, 수치해의 안정적 수렴여부를 판단하기 위해 Von neumann 안정성 해석을 수행하였다. 또한 기법의 정확성을 검증하기 위하여 Gauss 분포의 초기 자유수면변위를 갖는 문제에 적용하여 선형 Boussinesq 방정식의 해석해와 비교하였다. 그 결과 기존의 선형 천수방정식을 차분화한 수치모형에 비하여 정확한 결과를 제공하였고 분산보정기법을 이용한 수치모형과 동일한 정확도를 보였으나 본 수치모형을 이용했을 때 비교적 넓은 범위의 조건에서 정확도 높은 결과를 제공하였다.

  • PDF

지진해일 전파 모의를 위한 유한차분모형 (A Finite Difference Model for Tsunami Propagation)

  • 안성호;하태민;조용식
    • 한국방재학회:학술대회논문집
    • /
    • 한국방재학회 2010년도 정기 학술발표대회
    • /
    • pp.107.2-107.2
    • /
    • 2010
  • 본 연구에서는 지진해일의 전파 과정을 모의함에 있어 선형 천수방정식의 수치분산을 이용하는 기법이 아닌 선형 Boussinesq 방정식을 직접 차분하는 유한차분기법을 제안하였다. 기법의 정확성을 검증하기 위하여 Gauss 분포의 초기 자유수면변위를 갖는 문제에 착용하여 선형 Boussinesq 방정식의 해석해와 비교하였다. 그 결과 기존의 선형 천수방정식을 차분화한 수치모형에 비하여 정확한 결과를 제공하였고 분산보정기법을 이용한 수치모형과 동일한 정확도를 보였으나 본 수치모형을 이용했을 때 계산 효율이 개선되었다.

  • PDF

Boussinesq 방정식을 이용한 지진해일 전파모형 (Tsunami Propagation Model Using Boussinesq Equation)

  • 송민종;하태민;조용식
    • 한국방재학회:학술대회논문집
    • /
    • 한국방재학회 2011년도 정기 학술발표대회
    • /
    • pp.57-57
    • /
    • 2011
  • 지진해일은 진행속도가 빠르고 파장이 길며 파형의 변화 없이 먼 거리를 진행 할 수 있어 주변지역은 물론 멀리 떨어진 지역에도 심한 범람피해를 야기시킨다. 지진해일의 일반적인 특징으로 장파와 단파가 합성되어 있고 먼 거리를 전파할 경우 분산효과의 역할이 중요하게 된다. 특히 우리나라의 동해안에 영향을 주는 지진해일은 단주기파 성분이 강하고 파장에 비해 먼 거리를 전파하기에 분산을 고려하는 선형 Boussinesq 방정식을 지배방정식으로 사용하는 것이 바람직하다. 하지만 지금까지의 지진해일 전파모의를 위한 모형은 선형 Boussinesq 방정식의 복잡한 계산과 계산시간이 길다는 단점 때문에 선형 천수방정식을 지배방정식으로 사용하고 분산효과는 수치분산을 이용하여 고려해왔다. 지진해일 해석 시 일반적으로 사용되어 오던 기존의 leap-frog 유한차분 모형(Imamura et al., 1988; 조용식, 1996)은 지배방정식으로 선형 천수방정식을 사용하고 파의 분산효과는 수치분산을 이용하여 고려하므로 정해진 시간 간격에 대해 수심에 따라 격자 간격을 적절히 선택해야 하는데 수심이 복잡하게 변하는 경우 격자간격 조정이 불가능하여 분산효과를 정도 높게 고려할 수 없다. 이 문제점을 해결하기 위하여 윤성범 등(2004)은 파동방정식의 인위적인 분산항을 이용하여 Boussinesq 방정식의 분산효과를 고려할 수 있는 능동적인 분산보정기법을 제안하였고 Cho et al.(2007)는 일정한 수심에서 수치적인 분산오차가 Boussinesq 방정식의 물리적인 분산항을 대체하도록 수심, 격자 간격 및 계산 시간 간격 사이의 관계식을 유도하고 Boussinesq 방정식의 분산항과 일치하는 수치분산을 이용하여 실용적인 분산보정기법을 개발하였다. 이에 Ahn(2010)은 현재 컴퓨터의 계산 능력이 향상되어 선형 Boussinesq 방정식을 직접 차분하여 계산하는데 무리가 없다고 판단하여 선형 Boussinesq 방정식을 직접 차분한 모형을 개발하였다. 본 연구에서는 기존의 원해 지진해일 전파모의에 이용되어왔던 선형 천수방정식에 수치분산을 고려한 모형 대신 선형 Boussinesq 방정식의 유한차분 모형을 제안하였으며 기존의 선형 Boussinesq 방정식 모형의 격자와 수심간의 제약을 없애기 위해 차분 기법을 달리 한 2차 정확도의 유한차분 모형을 제안하였다. 검증을 위하여 선형 Boussinesq 방정식의 해석해(Carrier, 1991)와 비교하였다.

  • PDF

곡선좌표계에서의 2차원 흐름의 수치해석 (Numerical Analysis of Two-dimenstional Flow in Curvilinear Coordinate System)

  • 김형준;조용식;김수진
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2006년도 학술발표회 논문집
    • /
    • pp.402-407
    • /
    • 2006
  • 본 연구는 곡선좌표계에서 유한차분기법(finite difference method)을 이용하여 2차원 흐름이 모의가능한 수치모형을 개발하는 것이다. 기존의 연구는 대부분 직교좌표계(cartesian coordinate system)에서의 격자망을 대상으로 개발되고 적용되었기 때문에 불규칙한 흐름의 경계 및 형상을 올바로 표현하기 어려웠다. 유한요소법이나 유한체적법같은 수치모의기법들이 개발되어 비구조격자체계를 구성하고 자연현상에 가까운 경계 표현할 수 있도록 개발되었다. 하지만 위의 기법들은 질량과 운동량과 같은 물리량을 보존하기 위해서 매우 조밀한 격자체계를 가져야만 한다. 이에 본 연구에서는 기존의 문제점들을 해결하기 위하여 곡선좌표계(curvilinear coordinate system)를 이용하여 지배방정식을 표현하고 2차원 흐름을 모의할 수 있는 모형을 구축한다. 수치모형은 leap-frog기법과 1차 정확도의 풍상차분기법(upwind scheme)을 사용하여 구성하였다. 본 연구에서 개발된 모형을 사각수조 및 만곡수로흐름에 적용하여 모의결과를 해석해 및 실험관측값과 비교하였다. 이로부터 본 수치모형이 해석해 및 실측치와 잘 일치하고 있음을 알 수 있었다.

  • PDF

옵션에 대한 수치해법상의 초기값 불연속성 문제에 관한 연구

  • 김동석;변석준
    • 한국경영과학회:학술대회논문집
    • /
    • 한국경영과학회 1998년도 추계학술대회 논문집
    • /
    • pp.97-100
    • /
    • 1998
  • 옵션의 가격을 계산하기 위한 수치해법은 크게 격자모형, 유한차분법, 그리고 몬테카를로 시뮬레이션의 세 가지로 분류된다. 유한차분법은 옵션가격함수가 만족하는 편미분 방정식의 모든 편도함수를 유한 차분식으로 근사하여 옵션을 평가하는 방법이다. 본 연구에서는 유한차분법을 이용하여 옵션을 평가 할 때 발생하는 가격계산 오차의 가장 큰 원인이 옵션 만기 손익구조(payoff)의 비선형성에 있음을 보인다. 특히, 옵션 시장에서 가장 거래가 많이 이루어지는 손익분기옵션(at the money option) 그리고 손익분기점에 가까운 옵션(around the money option)에서 가장 큰 오차가 발생함을 보인다. 또한 본 연구에서는 이러한 오차를 효율적으로 줄이기 위하여 행사가격 근처의 일부 구간에서만 구간점 사이의 간격을 변화시키는 수정된 유한차분법을 제시하고 오차의 크기와 계산의 효율성 측면에서 기존의 유한차분법과 비교·분석한다.

  • PDF

자유수면흐름에 적용가능한 동수압 모형의 개발 (Development of Hydrodynamic Pressure Model for Free Surface Flow)

  • 이진우;조용식
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2011년도 학술발표회
    • /
    • pp.64-68
    • /
    • 2011
  • 본 연구에서는 일반 좌표계를 사용함으로써 자유수면이 매시간 변화하는 계산격자를 구성해야하는 문제점을 해결하였고, 또한 유한차분법의 단점인 지형변화의 적용성을 보완할 수 있었다. 자유수면과 동수압을 고려하기 위하여 계산단계를 3단계(정수압 계산단계, 동수압 보정단계, 자유수면 보정단계)로 나누어 해석하였고, 개발된 수치모형을 이용하여 실험값이 존재하는 수중방파제를 지나는 비선형파의 전파문제에 대해 수치모의를 실시하였다. 전반적으로 수치모의에 의한 결과는 실험값과 일치하는 경향을 나타내었다.

  • PDF

유속투사법을 이용한 동수압 모형의 개발 (Development of Hydrodynamic Pressure Models with Velocity Projection Method)

  • 이진우;김주영;이정규;조용식
    • 한국방재학회:학술대회논문집
    • /
    • 한국방재학회 2010년도 정기 학술발표대회
    • /
    • pp.52.2-52.2
    • /
    • 2010
  • 본 연구에서는 자유수면의 계산을 위해 동수압을 해석하는 수치모형을 제시하였다. 동수압과 자유수면을 고려하기 위해, 비점성 Navier-Stokes 방정식을 3단계로 나누어 해석하였다. 제1단계에서는 동수압과 자유수면을 전번 시간단계에서 계산된 값으로 대입하여 차분하였으며 이 차분식은 NGMRES(Newton-Generalized Minimal Residual) 방법을 이용해 음해적으로 해석되었다. 이때 계산된 유속장은 연속방정식을 고려하지 않았으므로 각 계산격자에서 질량보존법칙을 만족하지 않을 수도 있다. 제2단계에서 유속과 동수압 보정항으로 이루어진 식을 연속방정식에 대입하여 얻어지는 타원형 방정식인 동수압-포와송 방정식을 해석하므로 여기서 얻어지는 유속은 질량보존법칙을 만족하게 된다. 마지막 3 단계에서는 자유수면과 최종 유속을 계산하였다. 새로 개발된 수치모형을 검증하기위해 정사각형 탱크에서 수면의 자유 진동 문제에 적용한 결과 수치해는 해석해와 잘 일치하였다.

  • PDF

지진해일 전파 수치모형의 고도화를 위한 격자접속기법 (Grid Nesting Scheme for High Accuracy of Tsunami Propagation Numerical Model)

  • 임채호;배재석;전영준;윤성범
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2007년도 학술발표회 논문집
    • /
    • pp.2130-2134
    • /
    • 2007
  • 선형 Boussinesq 형태의 파동방정식을 지배방정식으로 사용하는 능동적인 분산보정 지진해일 전파 유한 차분모형의 정확도를 향상시키기 위한 새로운 동시격자접속기법을 개발하였다. 이 격자접속기법은 공간에 대해 3차 보간식을 사용하므로 짧은 파에 대해서 보간에 따른 오차를 최소화할 수 있고, 시간에 대해 2차 보간식을 사용하기 때문에 기존 기법에 비해 더 정확한 값을 얻을 수 있다. 개발된 격자접속기법의 정확성을 평가하기 위해 수중 원형천퇴상을 전파하는 Gaussian 형상의 가상지진에 대해 격자접속기법을 적용하지 않고 계산한 결과와 격자접속기법을 적용하여 계산한 결과를 FUNWAVE에 포함되어 있는 선형화된 Boussinesq방정식에 의해 계산된 수치해와 비교하였다. 그 결과 개발된 동시격자접속기법의 우수함이 검증되었다.

  • PDF

잠제에 의한 파랑변형 수치모의 (Numerical Simulation of Wave Deformation due to a Submerged Breakwater)

  • 하태민;조용식
    • 한국방재학회:학술대회논문집
    • /
    • 한국방재학회 2010년도 정기 학술발표대회
    • /
    • pp.63.1-63.1
    • /
    • 2010
  • 수중구조물에 의한 파랑의 변형을 예측하기 위해 3차원 수치모형을 도입하여 수치모형 실험을 수행하였다. 본 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(Spatially Averaged Navier-Stokes) 방정식의 해를 구하는 LES(large-eddy-simulation) 기반의 수치모형이다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two Step projection 기법을 사용하여 SANS 방정식을 계산하였으며, Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 선형파를 일정 수심상에서 조파시켜 해석해와 비교한 후 수중구조물이 설치된 지형에 적용하여 파랑의 변형을 수치모의하여 수리모형 실험 결과와 비교 및 분석하였다.

  • PDF

포물형방정식을 이용한 파의 굴절.회절모형 (Parabolic Model for Wave Refraction-Diffraction)

  • 정신택;채장원;안수한;정원무
    • 한국해안해양공학회지
    • /
    • 제2권1호
    • /
    • pp.28-33
    • /
    • 1990
  • 해저면 경사가 비교적 완만한 해역에서의 수심과 흐름에 의한 파랑변형 현상을 정확하고 효율적으로 계산할 수 있도록 포물형 유한차분 수치모형을 수립하였다. 모형의 기본식은 쌍무형 완경사 파동방정식으로부터 Pade'근사를 이용하여 유도한 포물형방정식이며 수치계산은 Kirby(1986) 모형의 차분식을 수정한 다음 계산치의 정확도를 수치실험을 통하여 비교ㆍ분석하였다.

  • PDF