• 제목/요약/키워드: 유한 차분 방법

검색결과 395건 처리시간 0.023초

DRASTIC 모델 및 지하수 수치모사 연계 적용에 의한 부여읍 일대의 지하수 오염 취약성 평가 (A Joint Application of DRASTIC and Numerical Groundwater Flow Model for The Assessment of Groundwater Vulnerability of Buyeo-Eup Area)

  • 이현주;박은규;김강주;박기훈
    • 한국지하수토양환경학회지:지하수토양환경
    • /
    • 제13권1호
    • /
    • pp.77-91
    • /
    • 2008
  • 본 연구에서는 지하수 오염 취약성 평가 방법 중 가장 활용도가 높은 DRASTIC 기법 및 지하수 수치모사 기법을 병행하여 충청남도 부여군 부여읍 일대의 지하수 오염 취약성 변화를 평가하였다. 특히, DRASTIC 기법에서 사용되는 7가지 수리인자들 중 취약성에 대한 민감도가 상대적으로 높고 비교적 장기간 자료의 획득이 용이하지 않은 '지하수위(depth to water)' 인자에 대하여 수치모사 기법을 적용하여 시-자료에 대한 동적 보간을 실시하였으며, 이를 통하여 시간에 따른 지하수위 변화가 지하수 오염 취약성에 미치는 영향을 검토하였다. 연구의 대상 지역은 부여군 부여읍 일대로 충적 대수층을 대상으로 기 실시된 지하수 수위자료를 보정대상으로 하였으며 연구지역의 지하수위 분포양상을 모사하기 위하여 유한차분모델인 MODFLOW를 이용하였다. 지하수 수치모사에는 정류 및 부정류 모사가 모두 실시되었다. 그 결과 기존 지구통계기법을 통한 지하수위 보간에서 가장 큰 가중치를 갖는 지형적 요인뿐만 아니라 유역 암상이나 지질특성 등과 같은 제반 수리지질학적 요인까지 복합적으로 반영하는 지하수위 자료를 얻을 수 있었다. 기존 연구에서 수위자료를 보완할 때 주로 이용되는 기법인 크리깅은 수치모사와는 달리 강우나 계절적 영향, 암상이나 지질등과 같은 요인들은 반영하지 못한다. 또한 수치모사를 통하여 얻어진 수위자료는 부정류 모사가 가능함으로 인하여 풍수기와 갈수기 그리고 평수기의 지하수위 변화 등도 파악할 수 있다는 장점이 있다. 모사를 통하여 얻어진 지하수위 자료와 공동크리깅 방법을 통해 얻어진 지하수위 자료를 입력자료로 하여 연구지역에 지하수 오염 취약성을 비교 평가한 결과 평수기의 경우 유사한 경향의 지하수 오염 취약성을 보였다. 또한 공동크리깅을 통해 분석이 어려운 분기별 지하수 오염 취약성 평가 결과 강우량이 많은 여름철을 포함한 분기와 강우량이 가장 적은 겨울철을 포함한 분기가 뚜렷한 취약성 차이를 보이는 것으로 분석되었다. 본 연구에서 제시하는 바와 같이 지하수 오염 취약성 평가에 있어 가장 큰 민감도를 갖는 수위자료의 수가 공간적 혹은 시간적으로 부족한 경우 수치모사 적용을 통해 자료의 한계성을 보완할 수 있으며 시간적인 보간이 이루어질 수 있어 유용할 수 있을 것이라 판단된다.

대기 누출 방사성물질 선원 위치 추적을 위한 3차원 궤적모델 개발 (Development of Three-Dimensional Trajectory Model for Detecting Source Region of the Radioactive Materials Released into the Atmosphere)

  • 서경석;박기현;민병일;김소라;양병모
    • Journal of Radiation Protection and Research
    • /
    • 제41권1호
    • /
    • pp.31-39
    • /
    • 2016
  • 연구배경: 우리나라를 포함한 중국, 대만, 북한, 일본 등에서 원전, 재처리시설과 같은 원자력시설의 증가에 따라 주변국 핵활동 분석의 종합적 대책이 필요하다. 우리나라와 포괄적핵실험금지조약기구(Comprehensive Nuclear-Test-Ban Treaty Organization, CTBTO)는 동북아시아 지역에서 핵종 탐지소를 운영 중으로, 핵종탐지 장비에서 특이 값 측정시 모니터링 자료의 분석과 더불어 배출원 탐색모델을 이용하여 핵종의 기원이 어디인지 추정하고 평가하는 것은 주변국 핵활동에 대한 감시 및 안전성 확보 측면에서 중요하다. 재료 및 방법: 주변국의 은밀한 핵활동 시 방사성핵종의 기원을 추정하기 위하여 3차원 전진/후진형 궤적모델을 개발하였다. 개발된 궤적모델은 궤적 미분방정식을 유한차분법을 이용한 방법으로 주어진 바람자료를 이용하여 방사성핵종의 방출지점으로부터 입자의 궤적을 순차적으로 찾아가는 전진형 모델과 시간 역산으로 방출기원을 추정하는 후진형 모델로 구성되었다. 결과 및 논의: 개발된 궤적모델의 검증을 위하여 체르노빌 사고 당시 측정된 농도자료를 이용하였다. 검증결과 관측지점의 농도가 높게 측정된 지점과 방출기원에서 가까운 지역으로부터 시간 역산의 방출지점을 추정한 결과의 정확도가 높았다. 3차원 궤적모델은 방출시간, 방출높이, 방출간격 등의 변수에 의해 계산결과가 달라지는 불확도를 내포하고 있는데, 이러한 궤적모델의 불확도를 최소화하기 위해 한국원자력연구원에서 개발한 대기확산모델(long-range accident dose assessment system, LADAS)를 이용하여 fields of regards (FOR) 기법에 의해 오염물 방출영역을 추정한바 신뢰성 있는 결과를 얻었다. 결론: 본 연구를 통하여 개발된 배출원 탐색모델은 주변국의 은밀한 핵활동 시 핵종 탐지장비와 연계하여 방사성핵종의 방출지역과 기원을 파악하여 우리나라의 핵종탐지 능력을 향상하고 핵활동 및 방사선 안전 분야에서 주도적 역할을 할 수 있을 것으로 생각된다.

SUMT법(法)에 의(依)한 2골절(滑節) I형(形) 강재(鋼材) 아치의 최적설계(最適設計) (Optimum Design of Two Hinged Steel Arches with I Sectional Type)

  • 정영채
    • 대한토목학회논문집
    • /
    • 제12권3호
    • /
    • pp.65-79
    • /
    • 1992
  • 본(本) 연구(硏究)는 2골절(滑節) 강재(鋼材) 아치의 최적설계(最適設計)에 관(關)한 것으로 아치구조(構造)의 정확(正確)한 해석(解析)과 구조(構造)를 안전(安全)하며 경제적(經濟的)으로 설계(設計)하는 데 목적(目的)을 둔다. 구조해석(構造解析) 방법(方法)은 해석과정(解析過程)에서 구조물(構造物)의 처짐울 고려(考慮)하는 유한차분법(有限差分法)을 도입(導入)하므로 해석오차(解析誤差)를 소거(消去)하여 구조물(構造物)의 단면력(斷面力)을 결정(決定)할 수 있는 방법(方法)을 사용(使用)한다. 최적화문제(最適化問題)는 설계변수(設計變數)를 단면(斷面)의 칫수들(B, D, $t_f$, $t_w$)로 하는 목적함수(目的函數)와 제약건식(制約件式)으로 형성(形成)한다. 목적함수(目的函數)는 아치구조(構造)의 총(總) 중량(重量)으로하고 제약조건(制約條件)은 한국(韓國) 도로교(道路橋) 표준시방서(標準示方書)에 규정(規定)된 허용응력(許容應力), 플랜지와 복부(腹部)의 최소칫수에 관한 규준(規準)을 사용(使用)하고 I형(形) 단면(斷面)의 경제적(經濟的) 높이 조건(條件)과 복부(腹部)의 상한계(上限界) 칫수와 플랜지 폭(幅)의 하한계(下限界) 칫수를 포함(包含)하여 유도(誘導)된다. 본(本) 연구(硏究)에서 개발(開發)된 비선형계획문제(非線型計劃問題)를 풀기 위해 수정(修正) Newton Raphson 탐사법(探査法)을 사용(便用)하는 SUMT 기법(技法)을 도입(導入)하여 수치예(數値例) 통(通)하여 시험(試驗) 본다. 본(本) 연구(硏究)에서 개발(開發)된 아치구조(構造)의 최적화(最適化) 프로그램은 여러 아치구조(構造) 수치예(數値例)를 통하여 시행(試行)하고 고찰(考察)한다. 이러한 수치결과(數値結果)를 통(通)하여 본 알고려즘의 최적화(最適化) 가능성(可能性), 적용(適用) 가능성(可能性) 및 수검성(收檢性)과 타(他) 문현(文獻)(30)을 사용(使用)한 수치결과(數値結果)와도 비교분석(比較分析)한다. 본(本) 연구(硏究)의 최적단면적(最適斷面績)과 2차(次)모멘트의 상관관계식(相關關係式)은 많은 수치적(數値的) 최적설계(最適設計) 결과(結果)로부터 도출(導出)한다.

  • PDF

지표 현무암을 통해 전파하는 탄성파의 거동 - 석탄 탄성파탐사에 적용 (Seismic wave propagation through surface basalts - implications for coal seismic surveys)

  • Sun, Weijia;Zhou, Binzhong;Hatherly, Peter;Fu, Li-Yun
    • 지구물리와물리탐사
    • /
    • 제13권1호
    • /
    • pp.1-8
    • /
    • 2010
  • 탄성파 반사법 탐사는 석탄 구조를 묘사하기 위해 그리고 석탄 채굴법 중 장벽식 채단법(longwall mining)의 위험 경감을 위해 가장 널리 이용되며 가장 효과적인 방법 중의 하나이다. 그러나 이 탄성파 탐사법의 분해능은 화산암의 존재에 의해 영향을 받는다. 호주의 보웬분지(Bowen Basin)와 시드니 분지(Sydney Basin) 일부에 이러한 화산암이 분포하고 있으며, 이러한 지역에서 탄성파 탐사를 이용하여 지하지질구조를 정확히 밝혀내는 데는 어려움이 따를 수 있다. 결과적으로 그러한 지역은 석탄광 후보지로 관성을 끌지 못하게 된다 따라서, 이러한 현무암 지역에서 수행되는 탄성파 탐사의 성공률을 높이기 위한 기법들이 필요하다. 이 연구에서는 탄성파동 방정식에 기초한 모델링 기법을 이용하여 현무암의 물성차, 두계, 횡방향 연장성, 송신원에 대한 상대적인 위치, 다양한 형태의 불균질성이 탄성파 전파에 어떠한 영향을 주는가를 살펴본다. 탄성파 모델링 결과로부터 다음과 같은 결론을 얻었다. 1) 임피던스의 차가 큰 현무암과 다중 흐름(multiple flow)이 있는 경우 다중반사파가 강하게 나타나며 일차반사파가 약하게 나타난다. 2) 현무암층이 앓을 경우 현무암의 존재가 파의 전파에 미치는 효과가 두꺼운 현무암층에 의한 효과보다 적은 것으로 나타났으며, 3) 현무암이 부분적으로 덮여 있을 경우 표변전체가 현무암으로 덮여 있는 경우에 비해 그 효과가 미약하게 나타났다. 4) 저주파의 탄성파(특히 원거리 오프셋에서)가 고주파의 탄성파에 비해 현무암층을 더 잘 통과함을 알 수 있었다. 또한 5) 석탄층이 현무암층으로부터 얼리 떨어져 깊은 곳에 존재할수록 현무암층이 파의 전파에 미치는 영향이 적음을 알 수 있었다. 이러한 연구결과들은 현무암층 하부에 대하여 탄성파 탐사를 수행할 때 발생할 수 있는 문제에 대한 통찰력을 제공할 뿐 아니라 탄성파 잡음을 적절히 처리하고 저주파수 대역의 지오폰으로 원거리 오프셋에서 탄성파 탐사자료를 획득함으로써 탄성파 탐사 결과를 개선시킬 수 있다는 것을 보여준다.

소나무재선충병 피해를 받은 곰솔 원목의 열처리 소요시간 예측 (Prediction of Heat-treatment Time of Black Pine Log Damaged by Pine Wilt Disease)

  • 한연중;서연옥;정성철;엄창득
    • Journal of the Korean Wood Science and Technology
    • /
    • 제44권3호
    • /
    • pp.370-380
    • /
    • 2016
  • 제주도 지역에서 소나무재선충병 피해를 받은 곰솔 원목의 이용확대를 위하여 열처리를 수행하였다. 열처리는 처리 원목의 중심부가 소나무재선충의 사멸온도인 $56^{\circ}C$를 30분간 유지하여야 한다. 곰솔 원목의 초기함수율과 말구지름은 각각 46% ~ 141%, 180 mm ~ 500 mm의 범위이고, 기본비중과 전건비중은 각각 0.47, 0.52이었다. $105^{\circ}C$ 조건에서 함수율과 말구지름에 따라 열처리에 소요되는 시간은 7.7 h ~ 44.2 h의 범위로 측정되었다. 다양한 함수율 및 지름을 갖는 곰솔 원목의 열처리 소요시간을 예측하기 위하여 열처리 진행 중 처리목 내부의 온도분포를 유한차분법을 적용한 2차원 열전달 해석을 통하여 제시하였다. 열전달 해석을 위한 목재의 열적 특성은 함수율에 따른 열전도계수와 비열을 적용하였으며, 자연대류와 강제대류를 합한 형태의 혼합대류에 의한 혼합대류계수를 적용하였다. 실험값과 예측 값의 오차는 3 ~ 45%의 범위로 분석되었다. 곰솔 원목에서 초기함수율이 50%이고, 말구지름이 200 mm, 300 mm, 400 mm인 경우, 예측된 열처리 소요시간은 각각 10.9 h, 18.3 h, 27.0 h이었다. 초기함수율이 75%일 때, 지름에 따라 각각 13.6 h, 22.5 h, 32.8 h이고, 초기함수율이 100%일 때, 지름에 따라 각각 16.2 h, 26.5 h, 38.2 h이었다. 이러한 열처리 소요시간의 예측방법에 소나무와 잣나무 등 다른 소나무재선충병 피해목의 물리적 특성을 적용하면, 함수율과 말구지름에 따른 열처리 소요시간을 제시할 수 있을 것으로 판단된다.