• 제목/요약/키워드: shallow water equations

검색결과 173건 처리시간 0.022초

도시범람모의를 위한 다공성천수방정식의 고차 정확도 기법 개발 (Development of high-order method of porous shallow water equations for urban inundation modeling)

  • 정재영;황진환
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2022년도 학술발표회
    • /
    • pp.100-100
    • /
    • 2022
  • 일반적으로 유체와 구조물간 상호작용의 수리동역학적 모의에서는 벽경계조건을 통하여 유동에 대한 구조물의 영향이 반영된다. 하지만 도심지에서 발생한 홍수를 예측하려는 경우 이러한 방법으로는 밀집한 구조물들 사이에 형성된 좁은 길들로 인하여 세밀한 격자망을 요하여 큰 계산량을 유발하고 빠른 예측 속도를 기대할 수 없게 한다. 최근 이러한 문제를 극복하기 위해 성긴 격자망에서도 구조물의 유체에 대한 영향을 반영할 수 있도록 하는 방법들이 큰 관심을 받고 있다. 그 중에서도 다공성 천수방정식은 벽경계조건 대신 다공도(posority)의 개념을 이용한 모형으로 도시범람모의에 있어 계산량과 정확도를 가장 적절하게 타협한 모형으로 보고되고 있다. 이러한 흐름에 맞추어 본 연구는 다공도 천수방정식을 해석하는 수치 기법을 개발하였고, 여기에 최근 쌍곡선계 방정식의 수치적 연구들에서 소개된 주요 특징들이 반영되도록 설계하였다. 우선, WENO 기법과 Runge-Kutaa 기법을 통하여 공간과 시간에 대한 고차 정확도를 만족시켰다. 이 때, 재구성 변수와 알고리즘를 새롭게 제시하여 정상흐름조건에 대한 플럭스항과 생성·소멸항간 절단오차에 의한 비물리적인 흐름생성을 억제하였다. 또한, 수치모의 중 음수심의 발생으로 인하여 수치모형이 불안정해지는 현상을 막기 위해, 양-보존성 제한자를 구축하였다. 마지막으로 도심지에서 즐비한 인위적인 구조물에 의해 나타나는 지형적인 불연속의 효과를 적절하게 반영할수 있도록 정상파 재구축의 단계를 구축하여 수치 기법에 반영하였다. 이렇게 구성된 수치기법은 리만문제의 해석해에 기반하여 기존의 주요 연구들의 결과와 비교되었고, 그 결과 본 연구의 방법이 정확성, 수렴성, 안전성의 측면에서 가장 우수함을 수치적으로 증명하였다.

  • PDF

삼중대각행렬 시스템 풀이의 빠른 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프레임 이상의 놀라운 성능을 보여주었다.

가로흐름 수역으로 방출되는 2차원 표면온배수 수치모형 비교연구 (A Comparative Study of Two-Dimensional Numerical Models for Surface Discharge of Heated Water into Crossflow Field)

  • 이남주;최흥식;이길성
    • 한국해안해양공학회지
    • /
    • 제6권1호
    • /
    • pp.40-50
    • /
    • 1994
  • 가로흐름이 존재하는 천해역으로 방출되는 표면온배수에 의한 온도장의 정확한 예측을 위한 근해역 2차원 수치모형을 개발하였다. 개발된 모형은 4-방정식 난류모델로서 열적 시간상수에 대한 정보를 얻을 수 없는 2-방정식 난류모델의 단점을 극복하기 위하여 변동온도 자승 평균항 및 그것의 감쇠율에 대한 전달방정식을 2-방정식 모델의 전달방정식에 추가한 모델이다. 또한, 부력생성항 및 난류 열플럭스항을 도입하여 연직방향 확산현상을 고려하고자 하였다. 개발된 수치모형을 간단한 단면을 갖는 개수로 정상류의 경우에 대하여 적용하였으며, 그 결과를 기존의 실험결과 및 2-방정식 난류모델을 사용한 수치계산결과와 비교하였다. 4-방정식 모형에 의한 계산결과가 2-방정식 모형보다 실험결과와 잘 일치하였으며, 제트 포획 및 안정화 영역에서 온배수의 물리적 특성을 잘 재현하였다.

  • PDF

유한요소법에 의한 호소의 수질예측모형;삽교 담수호에 적용하여 (Water Quality Prediction Model in a Lake by Finite Element Method;Application to Sapkyo Lake)

  • 유병로;안상진
    • 한국환경농학회지
    • /
    • 제8권1호
    • /
    • pp.37-46
    • /
    • 1989
  • 본 연구는 호소의 2차원 수질해석을 위해 흐름의 기본 방정식과 확산이송 방정식에 유한요소법을 도입하여 모의한 것으로, 얻은 결론은 다음과 같다. 1.단순화된 1차원 수로에서의 모형검정 결과 수치해와 정확해가 거의 일치하므로 모형의 알고리즘이 잘 구성 되었음이 입증되며, 적용결과 실측치와의 비교에서도 만족할만한 결과가 도출되어 모형의 적용 가능성이 시사되었다. 2. 삽교호의 BOD, COD, T-N, T-P, 및 SS의 농도분포가 2차원적으로 해석되었다. 3. 비점오염원의 비중이 높아 오염부하량 산정이 어려운 T-N, T-P를 사용하지 않고 COD를 모의 하므로서 부영양상태를 예측할 수 있도록 하였다. 본 연구결과는 하구언 및 방조제 등으로 인한 호수와 같이 수심에 비해 폭이 넓은 자연수역의 수질 해석에 적용될 수 있으며 이를 활용하므로서 용수취수지점의 선택, 하천 방류수의 처리기준 선정등, 수자원의 효율적인 이용 관리와 수질 대책 수립에 이용될 수 있을 것으로 판단된다.

  • PDF

유한요소법을 이용한 천수방정식의 수치실험 (Numerical Experiments of Shallow Water Eqs. by FEM)

  • 최성욱;이길성
    • 대한토목학회논문집
    • /
    • 제10권4호
    • /
    • pp.141-150
    • /
    • 1990
  • 하구나 연안에서 해수의 순환형상을 모사(模寫)하게 위해 천수방정식(淺水方程式)을 여러 가지 경계조건 하에서 수치해석하였다. 공간영역은 Galerkin방법으로 이산화(離散化)하였으며 시간영역에 대해서는 유한차분법(Crank-Nicolson방법)을 사용하였다. 네 가지 검정실험이 해석적인 해가 있는 일차원 수로에서 행하여졌으며, 해석해를 구할 수 없는 이차원 모형에도 적용되었다. 해석해가 있는 경우 수치모사 결과가 이와 잘 일치하였으며, 이차원 모형에서의 결과도 매우 합당함을 알 수 있었다. 또 일차원 문제에서 4점 bilinear요소와 삼각형 요소를 사용한 결과를 각각 비교하였으며 시간적분도 2단계 Lax-Wendroff방법을 사용하여 결과를 비교하였다. 음해법을 사용할 경우 비교적 정확한 결과를 얻을 수 있으나 요소의 갯수가 많아지면 구성되는 대수방정식(代數方程式)이 커지기 때문에 각 시간마다의 계산량이 엄청나게 늘어나게 되며 양해법을 사용할 때는 원하는 만큼의 정확한 결과를 얻기 위하여 시간간격이나 공간격자 간격을 선정하는데 각별히 유의하여야 할 것이다.

  • PDF

파랑 변형 해석을 위한 복합 요소 모형 (Hybrid Element Model for Wave Transformation Analysis)

  • 정태화;박우선;서경덕
    • 한국해안해양공학회지
    • /
    • 제15권3호
    • /
    • pp.159-166
    • /
    • 2003
  • 본 연구에서는 천해 파랑 계산에 널리 사용되어지고 있는 확장형 완경사 방정식과 계산 효율은 같게 유지하면서 Laplace방정식을 직접 풀 수 있는 유한 요소 모형에 대해서 연구하였다. 기존의 확장형 완경사 방정식을 사용하는 경우와 같은 계산효율을 유지하기 위하여 파동장을 수심방향으로 1층인 유한요소로 나누고, 요소내의 포텐셜을 수면에 위치한 절점에 대한 포텐셜만으로 표시하도록 한 후, Galerkin 기법을 적용하여 수치모형을 구성하였다. 요소 내 수평방향에 대해서는 통상의 보간함수를 채택하였으며, 수심방향에 대해서는 진행파의 수심방향 거동인 함수를 사용하여 보간함수를 구성하였다. 모형의 개발은 우선 연직 2차원 문제를 대상으로 하였다. 개발된 모형의 검증을 위하여 연직 2차원에서의 파랑 반사 및 전달문제에 적용한 결과, 개발된 유한 요소 모형은 계산상의 효율면에서나 해의 정확도 면에서 기존의 확장형 완경사 방정식에 기초한 모형과 같은 수준을 보임을 확인하였다.

지진해일 전파 수치해석을 위한 실용적인 분산보정기법 적용 (Application of Practical Dispersion-Correction Scheme for Simulation of Tsunami Propagation)

  • 최문규;하태민;조용식
    • 한국방재학회 논문집
    • /
    • 제8권4호
    • /
    • pp.145-152
    • /
    • 2008
  • 지진해일의 초기 파장은 수십에서 수백 킬로미터에 이른다. 따라서 보다 정확한 지진해일 전파 수치모의를 위해 파장의 크기에 따라 중요성이 변하는 분산효과를 고려하여야 한다. 최근에 분산효과를 고려할 수 있는 실용적인 분산보정기법이 개발되었다(Cho 등, 2007). 분산보정기법은 일정 수심에서 해석해와 비교함으로서 검증되었다. 하지만 아직 수심이 변화하는 실제지형에 적용이 되지 않았다. 본 연구에서는 개발된 모형의 실제 지형에서의 적용성을 검토하였다. 모형의 적용성 검토를 위해 역사 지진해일인 1983 동해 중부 지진해일과 1993 북해도 남서 외해 지진해일을 재현하였다. 대상 지역은 당시 검조기록이 존재하는 속초, 묵호, 포항항이고 지진해일 첫 파 도착시간과 검조소 위치에서의 최대 파고를 비교함으로서 적용성을 검토하였다.

규칙파와 수평판의 비선형 상호작용에 의한 고차 조화항 발생 (Higher Harmonic Generation by Nonlinear Interaction between Monochromatic Waves and a Horizontal Plate)

  • 고혁준;조일형
    • 한국해안해양공학회지
    • /
    • 제19권5호
    • /
    • pp.484-491
    • /
    • 2007
  • 규칙파와 수면아래 잠긴 수평판 사이의 비선형 상호작용을 규명하기 위하여 수치파동수조를 이용하여 수치 실험을 수행하였다. 수치모델로 비선형 포텐셜 이론을 근거로 한 고차 경계요소법과 Navier Stokes 방정식과 엄밀한 자유표면 경계조건식을 푸는 CADMAS-SURF을 사용하였다. 두 모델 모두 수평판 위의 천수역에서 발생하는 고차 조화항들을 예측할 수 있으며 점성효과를 포함하는 CADMAS-SURF는 수평판 양단에서 발생하는 와류와 박리에 의한 고차 조화항의 발생을 재현할 수 있다. 수평판의 잠긴 깊이와 길이 변화에 따른 반사율과 투과율을 Patarapanich and Cheong(1989)의 모형실험결과와 비교하였고 서로 잘 일치함을 확인하였다. 수평판의 길이가 길어지고 잠긴 깊이가 얕아질수록 입사파의 에너지가 고차 조화항으로 더 많이 전이됨을 확인하였다.

BEM-CADMAS-SURF 결합해석법에 기초한 수치조파수조의 응용 (Application based on the strictly combined method of BEM and CADMAS-SURF)

  • 김상호;야마시로마사루;요시다아키노리;신승호;홍기용
    • 한국항해항만학회지
    • /
    • 제33권1호
    • /
    • pp.65-70
    • /
    • 2009
  • 대수심의 유체운동을 포텐셜 운동으로 가정하여 자유수면의 거동을 신속히 해석하는 BEM 해석법과 구조물 근방에서 유체의 자유 수면 변화를 계산하기 위하여 NS방정식의 해석으로 CADMAS-SURF 기법을 결합하여 하이브리드 수치기법을 개발하였다. 하이브리드 해석법에서는 반사파를 고려해야 할 넓은 영역에서, 대수심의 영역은 BEM이, 천수역은 CADMAS-SURF가 계산하게 된다. 특히, 하이브리드 모델은 장시간에 걸친 불규칙파의 운동에 대해서는 단독의 CADMAS-SURF을 이용한 계산에 비해 거의 동일한 정확도로 월등히 신속하게 계산할 수 있다. 본 연구에서는 완경사 해저면을 가진 넓은 해역에서, 호안구조물에 내습하는 파랑의 처오름과 월파와 같은 강비선형 파랑장 계산에 결합해석모델을 적용하였다. 계산결과는 각각 토요시마(풍도(豊島))의 규칙과 처오름 실험과 고다(합전(合田))가 제안한 불규칙파의 월파량 산정도와 비교하였다.

여울-소 구조에서 지표수-지하수 혼합대의 흐름 특성 분석에 관한 수치모의 연구 (Numerical Modeling of Flow Characteristics within the Hyporheic Zones in a Pool-riffle Sequences)

  • 이두한;김영주;이삼희
    • 한국습지학회지
    • /
    • 제14권1호
    • /
    • pp.75-87
    • /
    • 2012
  • 지표수-지하수 혼합대는 하천 및 호소 등에서 지표수와 지하수가 교환되는 공간이다. 지표수-지하수 혼합은 하상의 토층으로 확장되어 다양한 물리적, 생지화화적, 열역학적 교환을 발생시키며 수생태계 내 고유한 생태적 전이대를 형성하는데 주요한 역할을 한다. 과거 실험 및 수치모의 연구에 의하면 혼합대에서 발생하는 물질교환은 하천의 지형적인 특징으로 발생하는 압력분포에 의해 지배된다. 특히 하천의 구간 규모에서 여울-소 구조는 혼합대의 특성을 지배하는 주요 인자로 알려져 있다. 여울-소 연속 구조는 지표수에서 재순환영역과 정체점을 형성하며 이 독특한 흐름 구조에 의해 혼합대의 흐름특성이 영향을 받는다. 본 연구에서는 3차원 동수역학 모형을 이용하여 Reynolds-averaged Navier-Stokes 방정식과 Darcy 방정식을 동시에 해석하여 연속된 여울-소 구조에서 발생하는 지표수의 흐름구조가 혼합대의 흐름에 미치는 영향을 분석하였다. 모의 결과, 여울-소 구조에서 지표수의 재순환영역 및 정체점은 상승류와 하강류 형성과 직접적 연관을 가지며, 재순환영역의 크기가 감소하면 여울 전면부 하강류 형성 구간의 길이가 감소하고 최대 하강류 발생 지점이 하부로 이동하는 특성을 파악하였다. 이와 같은 본 연구의 결과는 혼합대의 현장관측, 하천 관리 및 복원 등의 연구에 활용하여 친환경 하천 조성에 기여할 수 있을 것으로 판단된다.