• 지구물리와물리탐사
• /
• 제1권1호
• /
• pp.25-30
• /
• 1998

공동충전 효과를 검증하기 위하여 실시한 시차 공대공 탄성파 탐사자료로부터 지하공동 부존 지역에서 충전 전과 후에 매질의 탄성파 전파속도의 변화를 확인하였다. 시차 공대공 탄성파 탐사자료에 나타난 반응과 시추조사 결과에 의하면 본 지역의 공동은 규모가 극히 소규모이거나 또는 폐석 등으로 충전된 것으로 보인다. 공동충진 효과는 토모그래피로부터 도출된 속도단면상의 탄성파 속도의 증가량을 분석함으로써 평가하였다. 시추공용 에어건을 진원으로 24-채널 하이드로폰을 수진기로 하여 자료를 취득하였다. 취득한 자료에는 무시할 수 없을 정도의 source statics를 확인할 수 있었다. 본 논문에서 제시한 보정방법은 2단계로; 1) 불규칙한 발파시점에 의한 영향 보정과 2) 잔여 정보정으로 이는 진원의 부정확한 위치에 대한 정보정이다. 본 논문에서는 고주파수 성분의 수치잡음이 억제되고 관심대상 부분에서 비교적 고분해능 영상을 도출할 수 있는 다단계 역산 방안을 제시하였다. 일반적으로 최소자승 주시토모그래피로는 평활화된 속도 영상을 얻을 수 있다. 따라서 이러한 역산으로는 비교적 소규모의 구간에서 발생한 적은 속도변화를 영상화하기에는 어려운 면이 있다. 본 논문에서는 속도모델의 파라메터를 변화시킨 2단계 제어 역산법으로 도출한 시차 토모그램으로부터 채굴 영향대에서 발생한 매질의 속도변화를 시각화 할 수 있었다. 2단계 역산법은 1-단계에서는 적정한 크기의 균일 격자로 구성된 모델을 사용하여 토모그램을 작성하고 이 토모그램에 2차원 중위수 필터를 적용하여 대략적인 속도구조 모델을 작성한다. 2-단계 역산시는 1-단계에서 작성한 속도모델을 수정하여 초기 모델로 한다. 모델 수정은 관심대상 부분만을 작은 크기의 균일격자로 재구성하는 것이다. 기준조사 토모그램을 2차 조사자료 역산의 초기 속도모델로 사용하였다. 속도변화는 공동대 부근에서만 예상되므로 그 이외 부분의 속도는 기준 토모그램과 동일하게 고정시키고 역산을 수행하였다.

• 대한토목학회논문집
• /
• 제33권6호
• /
• pp.2211-2221
• /
• 2013

4대강 살리기 사업을 통해 가뭄 및 홍수문제를 예방하고 수질개선 및 종합적인 물 관리를 위해 4대강에 다수의 보가 설치되었다. 그러나 보의 설치로 인해 상류로부터 유입되는 유속이 감소되어 유입된 유사의 이송능력 또한 저하됨으로써 보 상류구간에서의 하상변화가 발생한다. 특히 함안보는 낙동강에 설치된 8개의 보 중 낙동강 최하류단에 설치되어 완만한 경사와 유속 감소로 인해 하상변동 및 지속적인 퇴사문제와 더불어 하상의 불안정화가 더욱 클 것으로 예상된다. 이에 본 연구에서는 2차원 모형인 CCHE2D를 이용하여 함안보 설치에 따른 보 상류에서의 흐름 및 하상변동 분석을 실시하고 이를 통해 하상의 안정화를 위한 방안들의 정량적인 평가를 실시하였다. 함안보 설치 후 흐름 및 하상변동 모의 결과, 모든 모의유량 조건에서 초기하상에서의 유속이 일정 지속기간이 지난 후 하상에서의 유속보다 빠르게 나타났으며, 전체 모의구간에서 하상변동이 크게 발생한 지점의 위치가 모든 적용유량에 대해 동일하게 나타났다. 이에 따라 하상안정화를 위해 하류단 수위를 함안보 관리수위 5.0 m에서 4.5 m로 저하시킬 경우, 유속분포는 관리수위일 때보다 전반적으로 빠르게 나타났으며 하폭이 가장 좁은 지점에서의 침식현상은 하류단 수위저하에 관계없이 지속적으로 나타났다. 이에 본 연구에서는 하폭이 가장 좁은 지형의 하폭을 확대시키는 방법을 제안하였으며 수치모의 분석 결과, 하폭 확대 후 지형에서 지속적인 침식이 예상되는 구간에서의 하상변동은 거의 발생하지 않아 하상의 안정화 효과가 있는 것으로 나타났다.

• 한국지하수토양환경학회지:지하수토양환경
• /
• 제9권1호
• /
• pp.79-86
• /
• 2004

암석 내 균열을 따른 수리전도도는 균열의 기하학적 요소, 즉 방향, 간극, 거칠기 그리고 상호 연결도에 주로 좌우된다. 따라서, 균열 내 투수계수를 정확하게 계산하기 위해서는 이와 같은 기하 요소들을 최대한 계산모델에 반영할 필요가 있다. 이 연구에서는 균열 기하양상을 최대한 정확히 반영한 균열모델에서 기존 수치해석과는 다른 새로운 방법인 균질화 해석법(homogenization analysis method)을 이용하여 균열을 따른 투수계수를 구하기 위해 수치해석을 수행하였다. 먼저, 공초점 레이저 스캔 현미경(Confocal Laser Scanning Microscope)을 이용하여 암석시료의 균열 조도와 균열에 가한 수직압축력의 변화에 따른 간극 변화량을 직접 측정하고, 이와 같이 획득한 자료는 균열모델 재현을 위한 입력자료로 사용되었다. 재현된 균열모델을 토대로 한 균질화 해석법은 미시규모(microscale) 매질특성과 거시규모(macroscale) 매질특성을 동시에 고려하여 투수계수를 계산할 수 있는 것이다. 즉, 균질화 해석법은 주기적 미세구조(microstructure)를 갖는 미소 불균질 물질의 거동특성을 구명하기 위해 개발된 새로운 형태의 섭동(perturbation) 이론이다. 이는 균질한 미시규모에서 미시 투수특성을 계산한 후, 거시규모에서의 균질화 투수계수를 계산하게 된다. 그러므로, 이 방법은 균열 기하양상의 국부적 영향을 고려한 투수특성을 정확히 해석할 수 있다. 균질화법을 이용한 투수계수 산정결과를 기존 연구에서 제안한 경험식과 비교하여 그 타당성을 검증하기 위해 전술한 2차원 균열모델을 이용한 투수계수 계산을 수행하였다. 균열모델은 거칠기(roughness)를 반영하고 동일한 간극을 할당한 평행판 모델을 가정하였다. 계산결과에 의하면, 균질화 해석법에 의해 계산한 C-투수계수는 실내투수시험에 의해 구한 투수계수와 같은 범위의 값을 가지거나 $10^1$ 정도의 차이를 보여, 그 계산결과는 타당하다고 볼 수 있다. 그러나, 균질화 해석법은 국부적으로 불균질한 균열 기하양상과 물질특성이 미시규모와 거시규모에서 모두 고려되므로, 이들 특성을 정확히 알고 있을 경우 기존에 제안된 경험식들에 의한 계산결과 보다 균질화 해석법의 결과가 훨씬 정확함을 주목하여야 한다.

• 한국해양학회지:바다
• /
• 제18권3호
• /
• pp.111-121
• /
• 2013

지역규모의 정확한 일기예보와 해양생태계 변화 이해에 있어서 수온 예측은 매우 중요하다. 황해는 조류가 매우 빠르고 탁도가 높다. 이러한 해역에서는 수치 모델의 수직 혼합 기법 및 해수의 탁도에 따른 수형(water type)이 수온 구조 결정에 많은 영향을 미친다. 수직 혼합 기법 변화와 탁도의 변화에 따른 황해 수온 모사의 민감도를 알아보기 위해 3차원 해양 순환 모델인 Regional Ocean Modeling System (ROMS)을 사용하여 수치 실험을 수행하였다. 수직 혼합 기법은 해양 순환 모델에서 많이 사용되는 Mellor-Yamada level 2.5 closure(M-Y)와 K-Profile parameterization (KPP)을 사용하고, 탁도는 Jerlov의 분류에 따른 수형 1, 3, 5를 사용하여 수치 실험을 수행하고 그 결과를 국립수산과학원에서 제공하는 정선 해양 관측 자료와 비교, 분석하였다. M-Y 기법은 수직적 혼합을 상대적으로 강하게 모의하였으며 그 결과로 저층수온이 높게 형성되었다. 높은 저층 수온은 탁도를 높게 설정하면 완화되지만 표층 수온이 높아지는 단점이 있다. KPP 기법은 M-Y 기법보다는 수직 혼합을 약하게 모의하고 이 약한 수직 혼합 때문에 황해 연안을 따라 형성되는 조석전선을 잘 재현하지 못하였으나, 저층 수온은 관측 수온에 더 가깝게 재현하였다. 결과적으로 황해 3차원 해양순환 모델실험에서 M-Y 기법은 수직 혼합이 잘 되어 표층과 저층의 수온 차이가 작게 나타나고, KPP 기법은 이와 반대로 모의하였다. 탁도의 영향을 표현하는 Jerlov 수형은 높을수록 일사량이 낮은 수심까지만 투과되어 성층을 잘 표현하였고, 낮을수록 깊은 수심까지 일사량이 투과되어 표층과 저층의 수온차를 작게 모의하였다.

• 대한환경공학회지
• /
• 제36권8호
• /
• pp.570-578
• /
• 2014

수심이 낮은 하구에서 바람이 국부적인 해수면 상승/하강에 미치는 영향을 연구하기 위해 Mobile Bay에 3차원 해수 순환모델을 적용하였다. Mobile Bay의 남단 경계면, 즉 northern Gulf of Mexico에서 시작된 조위는 Mobile River system 북부까지 직접적인 영향을 준다. 그러나 Mobile Bay 남단에서 발생한 조위변화는 Mobile Bay 북부로 이동하면서 Mobile River system으로부터 들어오는 담수와 국부적 바람의 영향으로 왜곡된다. Mobile Bay 남단에 위치한 기상관측소에서의 바람정보를 Mobile Bay 전체에 적용하였을 경우 Mobile Bay 북부에서 실제보다 강한 바람의 영향으로 과대한 수위 상승과 하강 현상이 발생하였다. 그러나 Mobile Bay 남단과 중단에 위치한 두 개의 관측소에서 측정된 바람 정보를 활용하였을 경우 Mobile Bay 북부에서의 조위 변화를 보다 정확하게 재현하는 것으로 나타났다. 특히 바람의 강도가 센 경우 Mobile Bay 남단과 북단의 풍속이 현저하게 차이나는 것으로 나타났으며(~ 88%), 이는 Mobile Bay 북단에서 나무와 건물 등의 영향으로 바람의 강도가 급격하게 줄어들었기 때문으로 판단된다. 따라서 Mobile Bay와 같이 수심이 낮고 국부적으로 풍속이 다른 하구 또는 만에서의 수위변화 재현 또는 예측을 위해서는 국부적인 바람 정보가 매우 중요한 것으로 연구되었다.

• 한국수산과학회지
• /
• 제34권1호
• /
• pp.43-50
• /
• 2001

규칙파, 불규칙파 및 파랑과 흐름이 공존하는 상태에서 프레임에 우리형 그물을 부착하여 구성되어 있는 가두리 시설의 운동특성을 분석하기 위하여 정사각형 및 원형 구조의 가두리 시설을 대상으로 예인 수조에서 수리 모형실험을 실시하였으며, 규칙파중 운동 특성은 선형 포텐셜 이론에 의한 수치 해석의 그것과 비교하였다. 그 결과는 다음과 같이 요약할 수 있다. 1) 규칙파중 모형 가두리 시설의 상하 및 종 동요는 전후 동요와는 달리 그물의 영향을 거의 받지 않았으므로 이 시설물의운동 특성중 가장 중요한 상하 동요를 해석하는 경우에는 그물을 제외하고 프레임만을 고려해도 된다는 것을 확인하였다. 2) 불규칙파중 및 파랑과 흐름이 공존하는 상태에서 모형 가두리 시설의 운동 특성은 입사파 주기의 2배 되는 고주파수에서 시설물의 고유 주기 등에 의해 동적 운동의 peak frequency가 나타남으로써 비선형 즉, 2nd order harmonic 성분이 존재하였다 3) 파동 및 흐름장에서 동요 특성이 정사각형에 비해 비교적 작게 나타난 원형이 파도, 조류 등 환경 조건이 거친 해역에 보다 적합한 것으로 판단되었다.

• 한국전산구조공학회논문집
• /
• 제26권2호
• /
• pp.99-112
• /
• 2013

이 논문은 부분공간 시스템 확인기법을 이용하여 전단빌딩의 강성행렬과 부재의 강성을 추정하는 기법을 소개한다. 시스템 행렬은 입력-출력 데이터로 구성된 행켈행렬을 LQ 분해와 특이치 분해를 통해 추정한다. 추정된 시스템 행렬은 닮음 변환을 통해 실제 좌표축으로 변환하고, 변환된 시스템 행렬로부터 강성행렬을 계산한다. 추정된 강성행렬의 정확성과 안정성은 행켈행렬의 크기에 따라 변한다. 전단빌딩의 기저 유한요소 모델을 이용하여 행켈행렬의 크기에 따른 강성행렬의 추정 오차 곡선을 구한다. 오차 곡선을 이용하여 목표 정확도 수준에 부합하는 행켈행렬의 크기들을 결정한다. 이렇게 선택된 행렬의 크기들 중에서 부분공간 시스템 확인의 계산비용을 고려하여 보다 적절한 행렬의 크기를 결정할 수 있다. 결정된 크기의 행켈행렬을 이용하여 강성행렬을 추정하고 추정된 강성행렬로부터 부재의 강성을 추정한다. 제안된 방법을 손상 전후의 5층 전단빌딩 수치 예제에 적용하여 타당성을 검증한다.

• 한국안광학회지
• /
• 제7권2호
• /
• pp.1-11
• /
• 2002

눈물 층을 사이에 두고 각막 위에 부착되어 있는 콘택트 렌즈(하드렌즈)에는 모세관 현상에 따른 장력이 렌즈 가장 자리에 균일하게 방사형으로 작용한다. 이러한 장력은 렌즈가 각막에 부착되는 힘의 근원이 된다. 렌즈착용 시 자체 무게에 의한 렌즈의 회전과 렌즈가 평형상태에 도달 한 후 각막에 부착되는 힘을 계산할 수 있는 방정식과 이의 수치 계산 컴퓨터 프로그램을 수립하였으며 이 컴퓨터 모델을 사용하여 렌즈 구경, 렌즈 B.C, 눈물의 점성, 눈물 층의 두께 등의 변수가 렌즈의 평형 위치 및 부착력에 미치는 영향을 모사(simulation)하였다. 일정 각막 B.C에서 렌즈 B.C 증가에 따라 눈물 충의 간격은 급격히 증가하는데 비해 렌즈의 직경 증가는 완만한 간격 증가를 수반한다. 렌즈의 B.C가 증가할수록, 즉 렌즈의 곡률이 완만할수록, 렌즈가 각막에 fitting되는 부착력은 급격하게 감소하며 직경 감소는 완만한 부착력 감소를 초래한다. 렌즈가 각막에 부착되는 순간 렌즈의 무게로 인하여 아래 방향으로의 회전 모멘트가 발생하게 되며 이와 동시에 렌즈 위 부분 눈물 층의 간격이 아래 부분보다 좁아지게 되어 위 방향으로 알짜 힘이 작용하게 되며 이에 따른 역방향 모멘트는 중력에 의한 모멘트를 상쇄하게 되어 렌즈는 일정 위치에서 평형 상태에 도달하게 된다. 렌즈 착용에 수반되는 초기 회전각 또는 회전 변위는 렌즈의 B.C가 증가할수록 급격히 증가하는데 비해 렌즈 구경 증가는 변위 증가율이 완만하다.

• 대한토목학회논문집
• /
• 제33권1호
• /
• pp.297-304
• /
• 2013

기후변화에 대응하여 수송부문의 탄소배출 저감노력이 세계적인 이슈가 되고 있다. 특히, 지역단위의 수송부문 탄소배출 저감전략의 수립 및 효율 집행을 위해서는 정밀한 탄소배출 인벤토리 구축은 필수적이나, 현재 도로교통의 특성과 수단의 다양성, 통행패턴의 복잡성 등 자료수집 한계로 신뢰성 있는 지역단위 인벤토리는 개략적 산정에 그치고 있는 실정이다. 본 연구는 이와 같은 문제점을 극복하기 위해 국가교통DB의 자료를 활용하여 지역단위 도로교통 탄소배출 인벤토리 구축 방법론을 제시하였다. 이를 위해 국가교통DB포함 자료 중 $CO_2$ 관련 속성을 파악하고, 이에 적합한 $CO_2$ 산정식을 도출하여 인벤토리 모형을 구축하였다. 구축된 모형의 적용성을 검토하기 위해 부산시를 사례로 2008년도 기준 도로부문의 탄소배출량을 산정하여 기존 연구결과와 비교분석하였다. 또한, 연구의 결과를 기초로 더욱 정밀한 인벤토리 구축을 위한 향후 연구방향을 제시하였다.

• 한국전산구조공학회논문집
• /
• 제33권3호
• /
• pp.171-178
• /
• 2020

본 연구에서는 다물체 페리다이나믹 해석 코드의 MPI-OpenMP 혼합 병렬화를 수행하였다. 페리다이나믹 해석 모델은 복잡한 동적파괴 거동 및 불연속 특성을 모사하는데 적합하지만, 비국부 영역을 통한 절점 간 상호작용을 계산하기 때문에 유한요소 모델에 비해 계산 시간이 많이 소요된다. 또한 다중적층구조물의 다물체 페리다이나믹 해석에서 추가된 비국부 접촉 모델과 가상 층간 결합 모델을 통한 여러 물체 간 상호작용으로 계산 부담이 증가한다. 더불어 고속 충돌 파괴와 같은 복잡한 동적 파괴 거동 해석을 위해 세밀한 절점 간격과 작은 시간 간격이 요구되기 때문에 코드 최적화와 병렬화를 통한 고성능 해석 코드 개발이 필수적이다. 해석 코드는 Intel Fortran MPI compiler와 OpenMP를 사용하여 개발되었으며, 한국과학기술정보원(KISTI)의 슈퍼컴퓨팅센터 누리온(Nurion)으로 실행되었다. 다물체 해석 코드를 최적화하기 위한 핵심 요소들을 분석하고, 모델 의존성 발생 서브루틴 분석 및 프로세스 통신 데이터 분별을 통해 MPI-OpenMP 혼합 병렬 처리 구조를 적용하였다. 다물체 충돌 파괴 현상 시뮬레이션을 통해 개발된 병렬 처리 코드의 성능을 확인하였다.