• 한국수자원학회논문집
• /
• 제31권5호
• /
• pp.599-612
• /
• 1998

수송방정식의 양면적은 특성으로 인하여 이송항이 지배적인 흐름에 있어서 수송방정식의 수채해석은 매우 난해하다. 특히 유한요소법을 사용하여 수치해석할 때, 상류방향으로 더 많은 가중치를 두기 위하여 변화된 시행함수를 사용한다. 이러한 방법을 Petrov-Galerkin 방법이라고 한다. 본 논문에서는 N+1 과 N+2 Petrov-Galerkin 방법을 격자 Peclet 수가 큰 수송문제에 적용하였다. 주파수맞춤 기법을 사용하여 N+2 Petrov- Galerkin 방법을 격자 Peclet 수가 큰 소송문제에 적용하였다. 주파수맞춤 기법을 사용하여 N+2 Petrov-Galerkin 방법의 적정 풍상정도를 찾아내었고, 그 결과를 토의하였다. 이 기법의 시행함수는 이송항과 확산항의 상대적 크기에 따라 그 모양이 변화된다. 수치실험을 통하여 제시된 새로운 수치해석기법의 우수성을 설명하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2019년도 학술발표회
• /
• pp.98-98
• /
• 2019

중력류 또는 밀도류는 주변 유체에 비해 상대적으로 밀도가 큰 유체가 밀도차에 의한 추진력으로 흐르는 것이다. 중력류의 수치모델링에는 두 가지 어려움이 있다. 즉, 적합한 지배방정식을 구성하여 적용하는 것 그리고 난류의 영향을 합리적으로 반영하는 것이다. 기존 중력류 해석을 위한 지배방정식들은 유체의 연속방정식과 운동량 방정식 그리고 밀도 또는 농도의 이송방정식을 조합하여 구성된다. 이들 지배방정식을 이용한 연구들은 대부분 두 유체 사이의 밀도차가 충분히 작아서 밀도 변동(variations)의 영향은 오로지 부력항에서만 유지된다는 Boussinesq 근사에 근거를 둔다. 그리고 이송방정식에서 밀도 또는 농도의 확산계수을 점성계수의 함수로 표현하기 위해서 Schmidt 수를 이용한다. 수치모델링에서 Schimdt 수는 상수값을 적용하지만, 이 값은 밀도의 연직방향 경사에 근거한 부력빈도(buoyancy frequency)와 난류량의 따라 큰 차이를 보이는 것으로 알려져있다. 한편, 표준 통계학적 난류모델과 벽함수를 적용한 수치모델링은 초기 중력에 의해서 무너지는(slumping) 단계를 넘어 관성력으로 추진되는 단계와 점성 효과가 지배적인 단계에서는 정확도에 현저히 낮아지기 때문에 대부분 큰와모의(large-eddy simulation, LES) 또는 DNS(direct numerical simulation)수준의 고해상도(high-resolution) 해석기법을 적용하여 공학적인 문제에 적용하는 데는 한계가 있다. 이 연구에서는 Boussinesq 근사와 Schmidt 수를 사용하지 않으며, LES 보다 적용이 용이한 DES (detached-eddy simulation)기법을 조합한 다상흐름 수치모델을 적용하여 중력류를 해석을 시도하였다. 수치해석결과를 실험값과 함께 기존 수치모델링 기법으로 구한 수치해와 비교분석하여 이 연구에서 개발 및 적용된 수치모델링 기법의 적용성을 평가한다.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2011년도 정기 학술발표대회
• /
• pp.44-44
• /
• 2011

저수지와 같은 갇혀진 수체는 상류에서 유입되는 오염물질 뿐만 아니라 성층현상에 의해서도 오염될 수 있다. 안정된 성층은 혐기성 조건을 제공하여 바닥에서의 과도한 조류성장을 유발하고 연직순환흐름을 방해한다. 갇혀진 수체에서의 연직순환은 환경적인 문제를 감소시킬 수 있도록 중요한 역할을 한다. 갇혀진 수체에서의 연직순환은 이러한 오염을 줄이는데 중요한 역할을 하는데, 연직순환을 일으키는 인자로는 빛의 입사, 바람, 물의 온도 및 열의 확산 등이 있으며, 그중에서도 가장 중요한 것은 바람의 영향이다. 그러므로 성층화된 흐름에서 바람에 의해 발생하는 연직순환에 대한 수치모형을 개발하고 적용하는 것이 필요하다. 본 연구에서는 수온성층흐름에서의 순환흐름을 해석하기 위하여 유속성분을 계산하고 자유수면 변위와 온도, 염도등과 같은 스칼라양의 해석을 위해 다음과 같은 3단계의 방법을 이용하였다. 첫 번째 단계(정수압 계산단계)에서는 운동량 방정식의 경사항을 음해적으로 해석하고, 두 번째 단계(자유수면 보정단계)에서는 자유수면의 변화를 계산하고 수평방향 유속성분을 계산한다. 예측-수정자 방법(predictor-corrector step method)을 이용하여 자유수면변위와 유속을 구하였다. 마지막으로 세 번째 단계(이송-확산단계)에서는 이송-확산 방정식을 이용하여 스칼라양을 계산하였다. 본 연구에서 제시한 모형의 정확도를 검증하기 위하여 정사각형수조에서 진동하는 자유수면의 해석해와 비교하였고, 성층화된 흐름에서 발생하는 연직순환에 대하여 수치모의를 실시하였다. 그 결과, 본 연구에서 개발된 수치모형이 흐름 내부의 현상을 잘 묘사함을 알 수 있었다.

• 한국수자원학회논문집
• /
• 제30권2호
• /
• pp.143-154
• /
• 1997

이송항에는 5차 보간다항식을 사용하는 Holly-Pressmann 기법을, 확산항에는 Hobson 등이 제안한 양해법을 사용하는 연산자 분리기법을 사용하여 1차원 이송-확산방정식의 수치모형을 제안하였다. 제안된 모형을 검정하기 위하여 일정한 유속과 종확산계수를 갖는 순간적으로 부하된 오염원의 경우와 상류단에 연속적인 오염원을 갖는 경우에 대하여 본 모형의 해를 해석해와 기존의 모형으로부터 구한 해를 비교검토하였다. Courant 수와 Peclet 수를 가진 경우에 대한 수치해석을 통하여, 본 모형이 Courant 수가 1보다 큰 경우에 대해서도 안정된 해를 제공함을 알 수 있었으며, 해석해가 존재하는 경우에 본 모형을 적용하여 얻은 수치해와 비교한 바 전반적으로 잘 일치하였다. 본 모형의 확산항에 사용된 양해법에서는 일반적인 양해법의 단점인 계산시간간격의 제약이 상당히 완화되어 상대적으로 큰 계산시간간격에 대해서도 양호한 결과를 보였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• pp.1219-1223
• /
• 2007

수온은 수계에서 가장 중요한 물리적 특성 중 하나로서, 생물 군집, 특히 어류와 무척추 동물에 관련된 많은 수질 인자에 영향을 미친다. 하천의 생태학적 모습을 개선하기 위한 하천 복원 사업 수행에 있어서, 서식처 및 산란처 조건으로서의 수온 조사 및 모델링의 필요성이 점차 증가하고 있는 상황이다. 본 연구는 낙동강의 중상류에 위치한 반변천 합류부에 평면 이차원 비정상 수치모형인 KU-RLMS 모형을 적용하여 수온의 변화 특성을 규명할 목적으로 수행하였다. KU-RLMS 모형은 하천 및 저수지의 국부적인 수리, 수질, 유사이동 해석을 위해 개발된 평면 이차원 비정상 수치모형이다. 직사각형 격자를 사용하는 유한차분법의 단점을 보완하기 위해, 수심적분된 2차원 연속방정식, 운동량방정식, 이송확산방정식을 불규칙한 경계를 현실적으로 모사할 수 있는 직교곡선 좌표계로 변환한 방정식을 사용한다. 이 모형은 흐름, 농도, 지형변화를 조합하여 계산할 수 있는 모형으로서 점착성 및 비점착성 유사의 이동, 보존성 및 비보존성 오염물질의 이동, 수온 변화를 모의할 수 있다. 수치모형 적용을 위한 현황분석으로 안동 및 임하 조정지댐의 방류량, 안동 수위관측소의 수위, 법흥교 및 포진교 지점의 수온 자료를 분석하였다. 이송확산모형의 보정을 위해, 안동대교 지점의 수온 횡분포 측정자료를 사용하여 확산계수에 대한 매개변수 추정 및 검증을 수행하였다. 또한, 안동조정지댐과 임하조정지댐의 방류량 및 방류수온을 고려하여 수치모의조건을 결정하였으며, 각 조건에 대한 수온 변화 특성을 분석하였다.

• 한국수자원학회논문집
• /
• 제45권2호
• /
• pp.179-189
• /
• 2012

오염물질의 이동 현상을 모의하기 위하여, 감쇠항이 있는 3차원 이송-확산 방정식의 수치모형이 개발되었다. 개발된 모형은 유한차분 모형으로서 시간단계의 가중치 ${\alpha}$를 포함하는 음해법(implicit finite difference method)과, 반복법인 Gauss-Seidel SOR(successive over relaxation)이 사용되었다. 모형은 보다 단순화된 가정 하에서 존재하는 두 가지의 해석적인 해와 비교되었다. 그 결과 Peclet number가 5~20 이하에서는 수치 분산의 영향이 크지 않았고 작은 오차범위 내에서 해석적인 해와 동일하였다. 또한 가중치 ${\alpha}$의 변화에 대한 모형의 거동은 Crank-Nicolson 모형(${\alpha}$=0.5)이 fully-implicit 모형(${\alpha}$=1)보다 해석적인 해에 접근함을 보여주었다. 모형의 검증과 실효성 제고를 위하여, mass balance를 검토하였다. 즉, 이송, 확산 및 감쇠항 각각에 대한 질량 이동을 산출하였으며, 그 결과 질량 이동의 계산 오차는 약 3% 이내였다. 본 모형은 감쇠 과정이 수반되는 3차원 이송-확산의 농도분포와 질량이동을 산출할 수 있으며 다양한 경계조건을 설정함으로서 현장조건을 반영할 수 있다. 그러나본 모형은 고정격자를 기반으로하는 유한차분 모형이므로 Peclet number가 비교적 작게 나타날 수 있는 토양 및 지하수계의 오염물질 이동 등의 문제에서 유용하게 적용될 수 있을 것으로 사료된다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2006년도 학술발표회 논문집
• /
• pp.1741-1745
• /
• 2006

본 연구에서는 하천에 유입되는 오염물질 중 대부분을 차지하는 비보존성 오염물질의 확산거동을 예측하기 위해 2차원 수심평균된 이송분산방정식에 유한요소모형을 적용하였다. 수치모형 구성을 위해 Galerkin법을 이용한 가중잔차법을 사용하였으며, 복잡한 하천경계를 보다 정확히 재현할 수 있도록 삼각 및 사각요소망의 혼용이 가능하도록 하였다. 모의대상 수질인자는 BOD, DO, 질소화합물, 인화합물, 수온, pH 및 대장균군수이며, 이 가운데 BOD와 DO는 상호 쌍을 이루는 방정식을 풀어야 하는 특수한 형태이므로 별도로 취급하였다. 순간주입 및 연속주입에 의한 비보존성 오염물질의 확산거동을 모의하였으며, 시간에 따른 민감도를 분석하기 위해 보존성 오염물질의 확산거동과 비교하였다. 해석해를 이용해 순간주입된 오염물질의 오염운을 구하는 식을 유도하고 해석해와 본 연구에서 개발한 수치모형에 의한 수치해를 비교하였다.

• 물과 미래
• /
• 제28권4호
• /
• pp.125-136
• /
• 1995

특정곡선을 고려한 ELM을 이송-확산방정식에 적용하여 그 결과를 Eulerian 기법(Stone-Brian, QUICKEST)과 비교하였다. 이송항의 계산을 위해서는 Lagrangian 보간법과 Cubic spline 보간법을 이용하였고 확산항의 계산에 있어서는 Crank-Nicholson 방법을 이용하였다. 수치모형의 적용결과는 다음과 같다. (1) Gaussian hill에의 적용:Lagrangian 보간법을 사용하여 계산한 경우가 가장 정확한 결과를 보였다. Cubic spline 보간법을 사용한 경우와 QUICKEST 방법의 경우에는 Peclet수가 50인 경우에 감쇠현상을 보였다. Stone-Brian방법은 Peclet수 10,50에서 위상오차가 발생하였다. (2) Advanced front에의 적용: 모든 방법이 Peclet수 1,4에서 정확한 결과를 얻었다. Peclet수가 50인 경우에 Lagrangian 보간법을 사용하여 계산한 경우와 Stone-Brian 방법은 증폭오차가 발생하였고 Cubic spline 보간법을 사용한 경우와 QUICKEST 방법의 경우는 수치진동 현상을 보였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2015년도 학술발표회
• /
• pp.169-169
• /
• 2015

유사 입자의 크기는 유사의 특성 및 그에 따른 거동 변화에 중요한 영향을 미친다. Stokes 침강 속도 모형에서 유사의 침강 속도에 가장 많은 영향을 주는 인자는 유사의 크기인 것 또한 확인된다. 유사 입자의 크기가 약 $60{\mu}m$보다 작은 유사들은 알갱이 사이의 점착력을 무시할 수 없다. 이로 인해 유사들은 응집 현상을 겪으며 입자 본래의 크기보다 크기가 큰 플럭을 형성하는 점착성 유사로 분류된다. 응집 현상이란, 흐름 내 점착성을 띠는 일차입자(Primary Particle)가 응집과 파괴를 반복하며 플럭을 형성하는 현상을 뜻한다. 입자 간의 충돌을 통해 응집이 진행되며 난류 전단으로 인해 형성된 플럭의 파괴가 발생한다. 많은 연구에서 점착성 유사의 충돌을 야기하는 가장 지배적인 원리는 난류라 알려져 있다. 이러한 응집 현상으로 인하여 플럭의 크기와 밀도는 지속적으로 변화를 겪으며 비점착성 유사와 다른 특징들을 보인다. 흐름에 존재하는 유사의 이동은 이송-확산 방정식을 통해 표현된다. 이송-확산 방정식은 시간 변화에 따른 농도의 변화를 입자의 침강과 난류 및 유사 자체의 특징에 의한 확산으로 해석한다. 침강속도로 대변되는 이송과 달리, 확산은 난류흐름 내에서 유사가 확산되는 정도를 정량화하기 위한 인자가 요구된다. 난류에 의한 유사의 확산은 유사 자체 특성에 따른 물질 확산에 비하여 매우 큰 값을 가지며, 이를 확산 계수로 개념화 한다. 확산계수는 와점성계수와 Schmidt 수(${\sigma}_c$)의 비로 정의된다. ${\sigma}_c$는 난류의 점성과 난류로 인한 부유과정에 의해 유사가 확산되는 정도를 나타낸다. 이에 따라 ${\sigma}_c$의 변화가 유사의 부유 및 침강거동에 많은 영향을 미칠 것이라 판단되나, 국내외에서 수행된 연구 동향에서는 ${\sigma}_c$를 0.5부터 1.0 사이의 상수를 적용하여 수행되었다. 이에 본 연구에서는 ${\sigma}_c$의 크기에 따라 달라지는 유사의 부유 및 침강 변화에 의한 총 부유량을 살펴보고자 한다. 유사의 점착성을 고려할 수 있는 1DV 수치 모형을 이용하여 비점착성 유사와 점착성 유사를 대상으로 수치연구를 수행하며, 유사의 크기 및 ${\sigma}_c$의 변화에 따른 총 부유량 경향을 살펴본다. 그 결과, 점착성 유사는 ${\sigma}_c$의 증가에 따라서 유사의 총 부유량이 증가하는 현상이 나타난 반면 비점착성 유사는 ${\sigma}_c$의 증가에 따라 유사의 총 부유량이 감소하는 경향이 나타났다. 그러나 크기가 아주 작은 비점착성 유사를 대상으로 수치 연구를 수행한 결과, ??에 따른 총 부유량의 경향은 유사의 점착성에서 기인하는 것이 아닌 입자의 크기로부터 야기되는 특성이라는 결론이 도출되었다.

• Ecology and Resilient Infrastructure
• /
• 제1권1호
• /
• pp.1-7
• /
• 2014

본 연구에서는 파랑과 흐름이 공존하고 있는 해안지역에 이용이 가능한 물에 용해된 물질의 정확한 이동을 예측하기 위한 수심적분형 수치모의 기법을 제시한다. 대상 영역에 일반적으로 발생하는 파랑의 전파와 변형 과정 및 쇄파와 흐름의 발달 과정에 대한 모의가 가능한 boussinesq equation 흐름모형과 동일한 과정을 거쳐 유도된 수심적분형 물질수송모형을 지배방정식으로 이용한다. 지배방정식은 approximate riemann solver를 이용하는 유한체적법을 이용하여 해석한다. 제시된 수치모형을 이용하여 해일발생에 의한 흐름양상을 계측한 실험을 재현하였으며, 해당 수역에 가상의 물질의 이송과 확산에 대한 수치모의를 수행하고 그 결과를 분석하였다.