• 한국수자원학회논문집
• /
• 제52권5호
• /
• pp.361-372
• /
• 2019

본 연구의 목적은 운동량방정식에서 이송가속도항을 제외한 지배방정식을 이용하여 정형 사각 격자 기반의 2차원 지표면 침수해석 모형을 개발하는 것이다. 공간적 이산화는 유한체적법을 이용하였으며, 시간적 이산화는 음해법을 적용하였다. 모형의 실행시간을 단축하기 위해서 CPU를 이용한 병렬계산 기법을 적용하였다. 개발된 모형의 검증을 위해서 해석해와 비교하고, 가상 도메인에서 수치실험을 통해 모형의 거동을 평가하였다. 또한 국내의 장호원 지역과 모로코의 Sebou 강 지역에 대해서 각기 다른 공간해상도로 침수해석을 수행하고, 그 결과를 CAESER-LISFLOOD (CLF) 모형을 이용한 해석 결과와 비교하였다. 모형의 검증 결과 해석해와 잘 일치된 모의 결과를 나타내었고, 가상 도메인에서의 흐름 해석도 타당한 것으로 평가되었다. 장호원 지역과 Sebou 강 지역에 대한 본 연구와 CLF 모형의 침수모의 결과는 침수심과 침수범위에서 서로 유사하게 나타났으며, 장호원 지역의 경우 홍수위험지도의 침수범위와도 유사한 값을 보였다. 본 연구와 CLF 모형의 모의결과에서 상이한 부분에 대해서는 각각의 모의결과를 비교 평가하였다. 연구결과 본 연구에서 제시된 모형은 홍수터에서의 침수 양상을 잘 모의할 수 있는 것으로 평가되었다. 그러나 본 연구에서 제시된 모형을 이용하여 침수해석을 할 경우에는 도메인 구성 방법과 지배방정식 및 해석 방법에 의한 모형의 특징과 한계점을 충분히 고려해야 할 것이다.

• 대한원격탐사학회지
• /
• 제36권6_3호
• /
• pp.1653-1667
• /
• 2020

본 연구에서는 국지예보시스템(LDAPS)과 전산유체역학(CFD) 모델을 접합하여, 부산 중구 광복동에 소재한 건물 밀집 지역의 상세 흐름과 PM2.5 농도 분포를 조사하였다. 도로 배출이 건물 밀집 지역의 PM2.5 농도에 미치는 영향을 분석하기 위해, PM2.5의 연간 시·군·구별, 배출 원소 별, 연료 별 도로이동오염원·비산먼지 배출량 자료와 월별·일별·시간 별 배출 계수를 이용하여 부산의 단위 면적당 시간별 PM2.5 배출량을 산정하였다. 본 연구에서는 건물 옥상과 도로변에서 수행된 특별 측정 자료를 이용하여 수치 모의 결과를 검증하고, 도로배출 유·무에 따른 PM2.5 농도 분포 특성을 분석하였다. 대상 기간(2020년 06월 22일) 동안 대상 지역에서는 바람이 약하게 나타났다. 새벽 시간에는 북동풍과 북서풍이 불고 주간에는 주로 남동풍이 불었다. 도로 배출을 고려하지 않은 경우에 LDAPS-CFD 접합 모델은 측정 지점(PKNU-AQ Sensor)의 PM2.5 농도를 과소모의 하였으나, 도로 배출을 고려하여 수치 모의한 PM2.5 농도는 도로 배출의 영향으로 PM2.5 농도가 증가하여 측정 결과와 유사하게 나타났다. 2020년 6월 22일 07시와 19시의 유입 풍향은 각각 북동풍과 남동풍이지만, 주변 지형과 건물에 의해 흐름이 변화되어, 두 시각 모두 측정 지점 주변에서는 주로 남풍 계열의 흐름이 나타났다. 07시와 19시의 유사한 흐름에 의해, 두 시각의 PM2.5 농도 분포도 매우 유사하게 나타났다. 건물 옥상 측정 지점에서 수치 모의된 PM2.5 농도는 도로 배출 영향을 크게 받지 않았으나, 도로변 에서는 도로 배출 영향을 상대적으로 크게 받았다. 도로 배출을 고려한 경우, 풍속이 약한 북쪽 도로와 긴 도로 협곡에 위치한 서쪽 도로에서 PM2.5 농도가 높고, 상대적으로 건물의 밀집도가 낮은 동쪽 도로에서는 PM2.5 농도가 낮게 나타났다. LDAPS-CFD 접합모델은 모든 도로에서 배출량이 동일하게 적용되기 때문에, 좁은 골목과 건물 밀도가 낮은 지역의 지형 특성이 반영되어 도로 별 PM2.5 농도 특성이 다양하게 나타났다.

• 한국전산구조공학회논문집
• /
• 제37권1호
• /
• pp.67-76
• /
• 2024

본 논문에서는 역학적 변수들을 측정하는 방안으로 디지털 이미지 프로세싱과 강형식 기반의 MLS 차분법을 융합한 DIP-MLS 시험법을 소개하고 추적점의 위치와 이미지 해상도에 대한 영향을 분석하였다. 이 방법은 디지털 이미지 프로세싱을 통해 시료에 부착된 표적의 변위 값을 측정하고 이를 절점만 사용하는 MLS 차분법 모델의 절점 변위로 분배하여 대상 물체의 응력, 변형률과 같은 역학적 변수를 계산한다. 디지털 이미지 프로세싱을 통해서 표적의 무게중심 점의 변위를 측정하기 위한 효과적인 방안을 제시하였다. 이미지 기반의 표적 변위를 이용한 MLS 차분법의 역학적 변수의 계산은 정확한 시험체의 변위 이력을 취득하고 정형성이 부족한 추적 점들의 변위를 이용해 mesh나 grid의 제약 없이 임의의 위치에서 역학적 변수를 쉽게 계산할 수 있다. 개발된 시험법은 고무 보의 3점 휨 실험을 대상으로 센서의 계측 결과와 DIP-MLS 시험법의 결과를 비교하고, 추가적으로 MLS 차분법만으로 시뮬레이션한 수치해석 결과와도 비교하여 검증하였다. 이를 통해 개발된 기법이 대변형 이전까지의 단계에서 실제 시험을 정확히 모사하고 수치해석 결과와도 잘 일치하는 것을 확인하였다. 또한, 모서리 점을 추가한 46개의 추적점을 DIP-MLS 시험법에 적용하고 표적의 내부 점만을 이용한 경우와 비교하여 경계 점의 영향을 분석하였고 이 시험법을 위한 최적의 이미지 해상도를 제시하였다. 이를 통해 직접 실험이나 기존의 요소망 기반 시뮬레이션의 부족한 점을 효율적으로 보완하는 한편, 실험-시뮬레이션 과정의 디지털화가 상당한 수준까지 가능하다는 것을 보여주었다.

• 대한원격탐사학회지
• /
• 제36권5_1호
• /
• pp.763-774
• /
• 2020

수치예보모델의 예측 바람장은 대기운동벡터 알고리즘의 표적 추적 과정에서 추적 정확도 향상이나 계산 시간 단축을 위해 초기 추정치로 사용된다. 대기운동벡터는 수치예보모델의 자료동화 시 활용가치가 높다고 알려졌으나, 초기 추정치로 사용된 수치예보모델 바람장이 대기운동벡터의 검증 과정에 참 값으로 사용된다는 모순이 있다. 이를 해결하기 위해서는 수치예보모델로부터 독립적인 초기 추정치가 필요하다. 본 연구에서는 Lucas and Kanade 옵티컬 플로우 방법을 적용하여 바람장을 도출한 후 이를 초기 추정치로 사용함으로써 표적 추적과정에서의 모델 의존성을 제거하고 계산 속도를 향상시키고자 하였다. 대기운동벡터 산출에는 2015년 8월 18일 ~ 9월 5일 00, 06, 12, 18시 동안의 정지궤도 위성 Himawari-8/AHI의 14번 채널 Level 1B 자료를 사용하였다. 옵티컬 플로우 방법이 대기운동벡터 산출에 미치는 영향을 평가하기 위하여 다음과 같은 세가지 방법으로 교차 검증을 수행 하였다. (1) 초기 추정치 없이, (2) KMA/UM 예보바람장을 초기 추정치로 사용하여, 그리고 (3) 옵티컬 플로우 방법으로 계산된 바람장을 초기 추정치로 사용하여 대기운동벡터를 산출하고 ECMWF ERA-Interim 재분석장과 비교 검증한 결과, 옵티컬 플로우 기반 바람장을 초기 추정치로 사용한 경우에 가장 높은 정밀도를 보였다(RMSVD: 5.296-5.804 ms^-1). 계산 속도는 초기 추정치를 사용하지 않은 경우에 가장 느렸고, 나머지 테스트는 유사한 속도를 보였다. 그러므로 대기운동벡터 알고리즘의 표적 추적 과정에 옵티컬 플로우 방법을 적용하면, 모델 의존성 없는 고품질 바람벡터의 산출이 가능할 것으로 사료된다.

• 한국방재학회 논문집
• /
• 제3권3호
• /
• pp.151-163
• /
• 2003

본 연구에서는 한국형 지하수 프로그램 개발(3-DFM, 3-Dimensional Finite Difference Method)을 위하여 대수층에 있어 지형 지질상태가 지하수유동시스템내에서 동적거동을 하는 것으로 취급하여 유동과정의 알고리즘을 확립토록 하였다. 본 연구에서 개발된 3-DFM모델은 입력변수 자료에 대한 설정이 모두 한글로 구성되어 있으며, 각 입력자료와 매개변수들의 이해와 적용치에 대한 도움말을 설정하여 두었다. 따라서, 입력변수에 대해서는 아이콘을 입력변수에 두면 각각에 대한 상세한 정보를 알 수 있도록 설계하였다. 또한, 각 지층의 지질경계 상태나 초기수위자료를 지정할 때는 work sheet상에서 간단히 지정할 수 있도록 설계되어 있다. 그리고 각 대수층의 특성과 더불어 정류 및 부정류 해석시에 각 매개변수들에 대한 입력은 기존의 모델과 같이 복잡하지 않도록 활성칸이 설정되도록 설계되어 있다. 최종 입력자료를 이용한 분석결과에서는 우측에 입력자료에 대하여 설명과 더불어 좌측에 분석 결과치를 나타나게 하였으며 이에 대한 결과는 TXT파일로도 출력할 수 있도록 설정하였다. 본 연구에서 개발된 모델은 유한차분법을 이용한 수치모델이며, 실제 함양량을 적용하고 매개변수들을 결정하여 관측 지하수두치와 모의발생으로 얻은 계산 지하수두치를 비교 분석하여 개발모델의 적용성을 검토하였다. 본 연구에서는 제주도 세화리 및 송당리일대의 양수에 따른 지하수 유동시스템 해석을 위하여 3-DFM모델을 적용 분석한 결과, 정류상태에서 따른 관측치와 계산된 지하수두와의 상대오차백분율(E.P.)이 $0.03{\sim}0.07$의 범위로서 관측치와 거의 일치하였다. 그리고 분석유역의 양수 전의 모의발생분석 결과를 이용하여 지하등수두분포와 유속벡터를 산정한 결과 지하수 유동분포는 높은오름과 문석이오름 등에서 월랑봉, 용눈이오름 및 손자봉 등 각 방향으로 고르게 유출되고 있는 것으로 분석되었다. 이러한 분석결과는 MODFLOW모델과 비교할 때 일치된 결과를 나타내었다.

• 한국의학물리학회지:의학물리
• /
• 제19권1호
• /
• pp.80-88
• /
• 2008

최근의 방사선 치료선량 계획시스템은 대체로 커널빔을 컨볼루션하여 조직선량을 구하고 있다. 본 연구에서는 광자선 빔에 따른 심부선량과 임의의 깊이에서 프로파일 선량을 구하기 위하여 반복적 수치해석을 통해 투과 필터에 의한 감쇠선량으로부터 에너지 스펙트럼을 구성하였다. 실험은 15 MV X선(Oncor, Siemens사)과 이온선량계 0.125 cc (PTW T31010)을 이용하여 납필터를 투과한 선량을 측정하여 이루어졌다. 15 MV X선의 에너지스펙트럼은 0.25 MeV 간격으로 납필터 0.51 cm에서 8.04 cm의 감쇠선량으로 실측치와 비교하여 구하였다. 실험 연산에서 15 MV X선의 최대유량은 3.75 MeV에서 나타났으며, 평균에너지는 4.639 MeV를 보였으며, 투과선량은 평균 0.6%의 오차인 반면에 최대오차는 납두께 5 cm에서 2.5%를 보였다. 조직선량은 에너지에 크게 의존하므로, 평탄형 필터의 중심과 Tangent 0.075와 0.125인 가장자리의 에너지를 구하였으며, 각각 4.211 MeV와 3.906 MeV로 나타났다. 심부선량과 프로파일 선량은 상업화로 공급되고 있는 선량계획시스템에 중심 선속과 가장자리의 각 에너지스펙트럼을 적용하여 구하여 실측선량률과 비교하였다. 생성된 심부선량 곡선은 조사면 $6{\times}6cm^2$에서 $30{\times}30cm^2$까지 실측치와 비교한 결과 1% 이내의 거의 일치하는 값을 얻었으며, 프로파일 곡선은 $10{\times}10cm^2$에서 1% 이내의 오차를 보였으나, $30{\times}30cm^2$와 같이 큰 조사면의 얕은 깊이에서는 2%의 오차를 보였다. 따라서 투과선량을 연산으로 구한 에너지 스펙트럼이 조직선량을 평가하는 데 상당히 적은 오차범위 내에서 정량적이고 정성적으로 얻을 수 있음을 알 수 있다.