• 대한지하수환경학회지
• /
• 제4권1호
• /
• pp.54-59
• /
• 1997

하천 주변의 쓰레기 매립장 등과 같은 오염원은 하천수와 지하수의 오염을 야기시킨다. 하천 수위의 영향을 받는 하천구역에서 오염원의 이동과 대수층 특성과의 반응을 고찰하였다. 이를 위해 하천 수위의 변화에 따른 변동속도 성분을 지하수의 지배방정식에 대한 해석해를 이용하여 유도하였다. 변동 속도 성분을 고려하여 이송확산방정식에 대한 수치모형을 구성하였으며 구성된 수치 모형의 타당성을 위해 해석해와 비교한 결과 거의 일치하였다. 수치실험을 위한 자료는 난지도의 대수층 특성값을 이용하였다. 대수층 특성값에 대한 범위를 설정하고 이를 이용하여 하천 수위와 오염원의 이동에 대한 반응을 고찰한 결과, 저류계수를 제외한 수리전도도와 유효공극율은 오염원의 이동에 대한 반응이 미소하게 나타났다. 그리고 하천 수위의 변화에 민감한 반응을 나타낸 경우는 저류계수가 $10^{-2}$ 차수임을 알았다.

• 한국통신학회논문지
• /
• 제19권8호
• /
• pp.1595-1611
• /
• 1994

본 논문에서는 전자파의 전파현상의 불연속모델로서 시간영역 유한 차분법의 수치적 성질이 연구된다. 시간 공간의 차원에서 막스웰 방정식을 개구리뜀 근사식으로 나타내므로 수치적인 특성과 의존 영역의 항으로 전자파의 전파현상을 모사한다. 시간영역 유한차분법의 수치적모사과정이 기하학적으로 설명된다. 개구리뜀 근사법의 채용으로 인한 수치적인 분산현상이 예시된다. 개구리뜀 근사법을 기초로 한 시간영역 유한차분법은 원래 계산 결과만을 산출하는 모델이 아니고 묘사적인 모델이므로 전자파 전파현상에 대한 몰리적인 현상을 묘사할 뿐만 아니라 이러한 묘사직언 결과로부터 푸리에 변환을 통하여 주파수 영역에서의 결과를 추출할 수 잇는 매우 유연한 수치해석 방법이다. 그래서 본 수치해석 방법을 이용하여 WR-28과 WR-90 도파관의 E-평면 휠터와 인턱티브 아이리스의 특성성분적 결과를 포함시킨다.

• 대한토목학회논문집
• /
• 제31권6B호
• /
• pp.507-514
• /
• 2011

본 논문에서는 하천에 유입되는 열 오염물질의 혼합거동을 분석하기 위해 2차원 수심 평균된 이송-분산방정식에 유한요소법을 적용하여 수치 모형을 개발하였다. 유한요소법의 여러 수치기법 중 SUPG법을 적용하였으며, 복잡한 하천경계를 보다 정확히 재현할 수 있도록 삼각 및 사각 요소망의 혼용이 가능하도록 하였다. 열 오염물질의 거동을 표현하기 위해서 열 교환을 묘사하는 반응항을 평형온도와 수온과의 차이에 비례하는 식으로 나타내고, 열교환 계수 및 평형온도에 따라 수온의 변화가 적용되도록 방정식을 구성하였다. 모형의 검증을 위하여 직사각형 수로에 선원으로 연속주입하여 얻은 수치해와 1차원 정상상태의 해석해를 비교하였다. 비교결과 수치해와 해석해의 결과가 서로 일치하는 것으로 밝혀졌다. 모형의 현장적용을 위해 상수원 보호구역인 팔당댐 하류부터 잠실수중보까지 22.5 km 구간을 대상영역으로 하였다. 구리하수처리장 방류수에 의한 수온 변화를 모의한 결과 수질측정망 측정자료와 비교적 비슷한 경향이 나타났다. 본 연구에서 개발한 수치모형이 열 오염원 유입으로 인한 수온 변화를 잘 표현하는 것을 알 수 있었다.

• 대한토목학회논문집
• /
• 제7권1호
• /
• pp.11-22
• /
• 1987

흐름수역(水域)에서 연직상향으로 방류되는 평면부력(平面浮力)?의 거동이 연속방정식(連續方程式), 운동량방정식(運動量方程式) 및 추적물수송식(追跡物輸送式)의 기본방정식을 수치적(數値的)으로 풀음으로서 해석(解析)된다. 난류확산(亂流擴散)에는 Prandtl의 혼합거리이론(混合距離理論)을 도입한 난류수송모형(亂流輸送模型)이 이용된다. 수치해과정(數値解過程)은 기본방정식을 유함수(流凾數)(stream function)식(式)과 골도수송(滑度輸送)(vorticity transport)식을(式) 이용하여 변환(變換)한 후, ?방류속도(放流速度), ?방류구폭(放流口幅) 등(等)으로 표현되는 변수(變數)와 흐름을 지배(支配)하는 무차원매개변수(無次元媒介變數)를 도입하여 무차원화(無次元化)하고 successive under-relaxation을 이용하여 Gauss-Seidal 반복법(反復法)으로 해를(解) 구(求)하는 것이다. 수치실험(數値實驗)은 방류(放流)Froude수(數)가 4~32, 방류속도(放流速度)와 가로흐름속도와의 비로(比) 정의되는 속도비가 8~15 의 범위의 흐름영역(領域)에서 수행되었다. 부력(浮力)?으로 인한 주변(周邊)흐름수역(水域)의 속도변화(速度變化), 온도상승(溫度上昇)범위, 흐름상태 및 골도(滑度)가 조사되었으며, ?의 경로에 대한 속도비와 방류밀도Froude 수의 영향이 또한 조사되었다. ?중심선의 속도, 온도변화, 국부밀도(局部密度)Froude 수(數)의 변화가 계산되며 퍼짐율(spreading rate)과 확산비(擴散比)(dispersion ratio)가 방류밀도(放流密度)Froude 수, 국부밀도(局部密度)Froude 수(數) 및 속도비(速度比)의 항(項)으로 해석되었다. 또한 속도와 온도분포를 상사(相似)(similarity)로 나타낼 수 있음이 밝혀졌으며, Gaussian 분포(分布)를 이용한 적분형해석(積分型解析)(integral type analysis)이 가능한 것으로 사료된다.

• 대한지하수환경학회지
• /
• 제3권2호
• /
• pp.51-59
• /
• 1996

지하 LPG 저장시설과 같은 지하수에 의존하는 시설물 인접지역에 폐기물매립장이 설치될 경우, 침출수 누출로 인한 지하수오염 예측 및 오염방지 대책수립이 요구되고, 이를 위한 지하수유동 및 용질 이동(Mass Transport) 해석을 실시하게 된다. 본 연구에서는 현재 U 지역 석유화학공단내에서 발생된 상기 와 같은 내용의 문제점을 해결하기 위하여 해석학적 해석 및 수치해석에 의한 침출수이동 예측을 실시하였으며 침출수제어 공법을 검토하였다. 해석학적 해석은 입력요소를 이송(Advenction)요인과 분산 (Dispersion)요인으로 구분하여, Peclet 수로 부터 각 요인의 침출수이동에 미치는 영향을 분석 검토하였으며, 보존성용질 이송확산 방정식으로부터 입력요소 변화율과 침출수이동속도 변화율의 관계를 함수식으로 도출하여, 입력요소 변화에 따른 침출수 이동속도 변화를 간단히 예측할 수 있었다. 수치해석은 지하수유동 및 용질이동 해석용 FEM 프로그램인 AQUA2D를 이용하여 침출수이동 예측을 위한 수치모의를 실시하였다. 침출수제어공법 검토를 위해서 3D 불연속체모형을 모의하고 모의된 모형을 대상으로 지하수유동해석을 실시하여 암반내 열극(Fracture)의 상호 연결성 분석을 통한 수벽(Water Curtain) 시스템의 타당성을 분석하였다. 해석학적 해석 및 수치모의 결과, 폐기물매립장으로 부터 지하 LPG 저장기지로 침출수가 30년 이내로 유입되는 것으로 예측되었으며, 이를 방지하기 위한 침출수 제어공법으로는 수직수벽 시스템이 효율적인 것으로 예측되었다.

• 대한환경공학회지
• /
• 제30권11호
• /
• pp.1146-1153
• /
• 2008

본 연구에서는 지하수 내 투수성 반응벽체(permeable reactive barrier, PRB)의 TCE 처리에 관한 모델링을 수행하여 trichloroethylene (TCE)의 농도, 컬럼의 단위 부피당 철 매질의 질량, 철환원균(iron-reducing bacteria, IRB)의 농도에 대하여 각각의 유기적인 관계를 고찰하였다. 1차원 이송 확산 반응 방정식을 MATLAB을 이용하여 이송, 확산, 그리고 분해 반응 등을 컬럼의 길이, 실험 수행 시간에 따라 모델하였으며, 유한차분법(finite differential method, FDM)으로 수치해를 구하였다. 영가철 및 2가 산화철은 TCE에 의한 반응항과 철환원균에 의한 반응 항으로 나누어서 식을 정리했다. TCE 주입농도는 10 mg/L로 설정하여 영가철 및 2가 산화철에 의한 각각의 관계를 모델링했다. 또한, 철환원균 농도와 산화철 환원 모델을 통해 철환원균의 농도에 따른 산화철 환원 효율을 해석했고, 이것이 전체 TCE 분해에 어떤 영향을 주는지 모델로 나타냈다. 영가철 컬럼에서는 TCE 제거 효율이 60시간에서 235시간 동안 99% 이상을 나타냈고, 1,365시간 이후에 1% 이하로 떨어졌다. 2가 산화철 컬럼의 경우 TCE와 반응을 시작한 210시간 이후에 평형을 이루었고, 85.3%의 일정한 제거 효율을 나타냈다. 모델의 결과에 따르면, 철환원균에 의한 2가 산화철의 경우 영가철보다 TCE 제거 효율이 떨어지지만 더 높은 제거수명을 가질 수 있는 것으로 나타났다.

• 대한원격탐사학회지
• /
• 제30권6호
• /
• pp.731-742
• /
• 2014

본 연구에서는 건물 꼭지점의 위 경도 좌표를 제공하는 GIS로부터 수치 모델의 건물 정보를 구축할 수 있는 알고리즘을 개선하였다. 이 알고리즘은 인접한 건물 꼭지점 위 경도 좌표를 지나는 선분을 순차적으로 연결하여 건물 외곽선을 구성하고, 외곽선 내부의 지점을 건물로 인식하기 때문에, Lee et al. (2009)에서 개발한 알고리즘의 한계를 개선하였고, 복잡한 형태의 건물을 실제에 가깝게 재현할 수 있었다. 본 연구에서는 GIS로부터 수치 건물을 구축할 때, 알고리즘 한계에 의해 발생한 건물 형태의 변화가 건물 주변의 흐름과 오염물질 확산에 미치는 영향을 조사하였다. 이를 위해, 알고리즘에 의한 건물 변형이 나타날 수 있는 세 가지 형태의 건물을 대상으로 전산유체역학 모델을 이용한 수치 실험을 수행하였다. 알고리즘 한계에 의해 발생한 건물 변형은 풍하영역의 흐름 패턴에는 상대적으로 작은 영향을 미쳤으나, 건물 사이의 공간에 나타나는 소용돌이와 같은 건물 규모 대기 현상의 수치 모의에는 매우 중요한 영향을 미쳤다. 건물 변형에 따른 건물 사이 공간의 축소는 건물 주위에 나타날 수 있는 소용돌이를 전혀 모의하지 못하거나 소용돌이 규모를 과소 모의 하는 등의 결과를 초래하였다. 건물 변형에 따른 평균 바람장 변화는 건물 주변 지역에서 배출된 스칼라 오염물질의 확산 패턴에도 큰 영향을 미침을 확인할 수 있었다.

• 한국해안해양공학회지
• /
• 제14권4호
• /
• pp.308-318
• /
• 2002

본 논문에서는 Korteweg de Vries(이하 KdV)방정식의 골격내에서 천해의 풍파의 수치모형요구조건에 대한 토의를 수행하였다. KdV식을 실험자료를 포함하는 발달된 해상상태, 비정상적 풍파와 나블상황에 적용시의 비선형성과 분산성의 취약점을 논하였다. 최극해파통계의 충분한 평가를 위해서는 파고기록이 적어도 10.000개 정도의 해파를 다루어야 하는데 이는 숫적으로 다루기 힘들다. 따라서 유의파의 2-2.3배에 상응하는 최극해파의 통계적 특성을 효과적으로 평가할 수 있는 가능성을 제시하는 불규칙해파의 수치적 모형을 제시하였다.

• Geomechanics and Engineering
• /
• 제2권2호
• /
• pp.141-156
• /
• 2010

In the framework of meshfree methods, a new methodology is developed based on radial point interpolation method (RPIM). This methodology is applied to a one-dimensional contaminant transport modelling in the saturated porous media. The one-dimensional form of advection-dispersion equation involving reactive contaminant is considered in the analysis. The Galerkin weak form of the governing equation is formulated using 1D meshfree shape functions constructed using thin plate spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Numerical examples representing various phenomena, which occur during migration of contaminants, are presented to illustrate the applicability of the proposed method and the results are compared with those obtained from the analytical and finite element solutions. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints. In order to test the practical applicability and performance of the RPIM, three case studies of contaminant transport through the landfill liners are presented. A good agreement is obtained between the results of the RPIM and the field investigation data.

• 한국수자원학회논문집
• /
• 제32권5호
• /
• pp.593-601
• /
• 1999

하천의 유속구조에 관한 지식은 하천 양안의 침식, 사행현상, 유사 및 오염물질의 확산과정을 규명하는데 필수적인 요소이다. 본 연구에서는 적절한 수리량을 사용하여 자연하천에서 종방향유속의 횡분포를 예측할 수 있는 이론식으로 수치모형을 구성하였다. 그리고 이를 낙동강에서 실측된 유속자료에 적용시켰다. 그 결과 수치모형은 실측된 유속분포의 일반적인 형태를 대체적으로 잘 근사하였다. 특히 모형은 유속의 횡분포에 있어서 왜곡도 및 편평도를 적절하게 재현하였다.