• 한국해안·해양공학회논문집
• /
• 제29권2호
• /
• pp.92-108
• /
• 2017

본 연구에서는 Navier-Stokes solver에 기초한 3차원혼상류해석법인 TWOPM-3D를 적용하여 고립파와 연안교량과의 상호작용을 수치해석적으로 수행하였다. 이 때, 고립파는 수치파동수로에서 수치조파되었고, 작용파력의 타당성은 본 수치해석결과와 기존의 실험결과와의 비교 분석으로부터 검증되었다. 이로부터 고립파의 파고, 교량의 종류와 위치 및 거더 수 등의 변화에 따른 수위변화와 유속변화를 포함하여 수평파력과 연직파력(연직상방파력과 연직하방파력)의 변동특성을 면밀히 검토하였다. 거더교의 경우 입사파고의 증가에 따라 수평파력과 연직파력의 각 최대치는 증가하는 경향을 나타낸 반면, 거더 수의 변화에 따라서는 서로 상이한 변동특성을 나타내었으며, 연직파력에 수중에 연행되는 공기가 큰 영향을 미치는 것 등을 결과로부터 알 수 있었다.

• 한국수자원학회논문집
• /
• 제56권3호
• /
• pp.173-184
• /
• 2023

기후 변화로 인하여 국내뿐만 아니라 전 세계적으로 녹조류 번무 현상이 잦아지면서 유해 녹조의 독성에 대한 위험성이 최근 대두된 바 있다. 대표적인 유해 조류인 남조류는 수온, 체류시간 및 영양염류 이렇게 3가지 요소의 교집합으로 발생한다고 알려져 있다. 본 연구에서는 내성천 내 영주댐의 완전개방으로 체류시간 감소에 따른 수질 변화 분석을 위해 준3차원 수치모델인 EFDC를 활용하여 수질모의를 수행하였다. 실제 2021년 여름철 영주댐 내 엽록소-a의 농도는 상당하여 과거 조류경보제 기준으로 '조류경보'를 오랜 기간 동안 상회하였다. 반면 댐 수문을 완전 개방한 조건으로 모의를 수행한 결과, 모의기간 동안 엽록소-a 농도가 대부분 '조류경보' 수준 이하로 저감되었다. 내성천의 흐름을 복원하여 체류시간을 감소시키는 것이 영주댐 내 녹조를 즉각적으로 저감시키는 방안임을 확인하였다.

• 대한환경공학회지
• /
• 제32권7호
• /
• pp.690-698
• /
• 2010

유기물에 의한 지하수 자원의 오염이 중요한 환경 문제로 대두되고 있다. 계면활성제 증진 대수층 복원(SEAR)은 유기 오염물에 대하여 가장 유망한 원위치 복원법 중 하나로 인식되고 있다. 계면활성제 또는 계면활성제/폴리머 용액은 복원 프로세스를 촉진시켜 물 세척 시보다 처리 시간을 크게 감소시킨다. 유기물로 오염된 대수층 복원을 위하여 계면활성제 기반의 처리법을 설계할 때 충분한 수치 시뮬레이션을 통한 사전 분석이 매우 중요하다. 본 연구에서는 3차원 다상 다성분 유한 차분 모델인 UTCHEM을 이용하여 유기물 오염 대수층에 대한 SEAR 과정 중 미세에멀젼의 형성과 유동을 분석하였다. 본 연구에서 고려한 복원 프로세스 변수로는 주입 유체의 배열, 유체 주입 및 추출율, 계면활성제 슬러그 및 추격수 내 폴리머 농도, 계면활성제 슬러그 주입 시간 등이다. 각 변수에 대하여 주입 및 생산정에서의 변화와 유기상의 공간적 분포를 비교하였다. 이 결과로부터 정화 시간 및 누적 유기물 회수량을 검토하였다. 본 연구의 결과는 비수상 액체로 오염된 대수층의 정화 시 설계 전략을 수립하는데 중요한 정보를 제공한다.

• 한국가스학회지
• /
• 제2권1호
• /
• pp.28-39
• /
• 1998

본 연구에서는 $20\%$ 천장개구부가 있는 정사각형 밀폐공간내의 순수자연대류와 자연대류 -복사가 고려된 복합열전달을 순차해석과 실험을 통하여 비교 분석하였다. 수치해석은 순수자연대류에 대하여 SIMPLE 알고리즘을 사용하였고, 복사열전달에 대해서는 S-N 구분 종좌표법을 이용하였으며 난류유동의 경계조건은 벽함수를 적용하였다. 실험은 수치해석의 결과와 비교하기 위하여 동일한 조건에 대하여 수행되었다. 그 결과 순수자연대류와 복합열전달의 유동장, 온도장의 형상은 유사한 유선함수를 보이고 있으며, 유동가시화를 통한 실험결과와 잘 일치하고 있음을 보여준다. 수치해석과 실험의 온도분포를 비교한 결과 평균 $8.5\%$의 오차를 보였다.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2011년도 정기총회 및 추계학술대회 논문집
• /
• pp.1199-1213
• /
• 2011

이수가압식 쉴드공법은 사질토에서 적용성이 우수하지만, 이수압이 낮으면 이수 유출 및 지반변형이 발생하기도 한다. 따라서 이수가압식 쉴드에서는 초기 막장압보다 큰 초과 이수압을 가하여 막장 안정을 유지한다. 그러나 이수압이 너무 높으면 전방 지반의 수동 파괴를 유발하므로 수동 파괴 위험성을 줄이고 이수압을 증가시키는 방법으로 막장 전방에 수평 차수층을 설치하는 방안이 있으나 위치와 규모, 효과가 잘 알려져 있지 않다. 따라서 본 연구에서는 포화 사질토에서 막장전방에 차수 그라우팅 적용시 발생하는 효과를 규명하기 위하여 모형실험을 수행하였다. 실험 결과, 차수층이 있는 경우가 없는 경우보다 이수압을 크게 가할 수 있어서 전방 차수층이 막장 안정성을 증대시킬 수 있다는 것을 확인하였다. 막장 안정성 증대에 가장 큰 영향을 주는 적정 차수층은 길이 1.0~1.5D, 설치높이 1.0D로 나타났다. 초기 막장압 대비 최대 이수압의 비로 막장의 자립 안전율(F)을 제안할 수 있으며, 전방 차수층을 적정 위치에 설치할 경우 초기 막장압보다 3.5~4.0배 크게 이수압을 가할 수 있는 것으로 나타났다.

• 한국방사성폐기물학회:학술대회논문집
• /
• 한국방사성폐기물학회 2009년도 학술논문요약집
• /
• pp.187-187
• /
• 2009

Opening of fractures induced by shear dilation or normal deformation can be a significant source of fracture permeability change in fractured rock, which is important for the performance assessment of geological repositories for spent nuclear fuel. As the repository generates heat and later cools the fluid-carrying ability of the rocks becomes a dynamic variable during the lifespan of the repository. Heating causes expansion of the rock close to the repository and, at the same time, contraction close to the surface. During the cooling phase of the repository, the opposite takes place. Heating and cooling together with the, virgin stress can induce shear dilation of fractures and deformation zones and change the flow field around the repository. The objectives of this work are to examine the contribution of thermal stress to the shear slip of fracture in mid- and far-field around a KBS-3 type of repository and to investigate the effect of evolution of stress on the rock mass permeability. In the first part of this study, zones of fracture shear slip were examined by conducting a three-dimensional, thermo-mechanical analysis of a spent fuel repository model in the size of 2 km $\times$ 2 km $\times$ 800 m. Stress evolutions of importance for fracture shear slip are: (1) comparatively high horizontal compressive thermal stress at the repository level, (2) generation of vertical tensile thermal stress right above the repository, (3) horizontal tensile stress near the surface, which can induce tensile failure, and generation of shear stresses at the comers of the repository. In the second part of the study, fracture data from Forsmark, Sweden is used to establish fracture network models (DFN). Stress paths obtained from the thermo-mechanical analysis were used as boundary conditions in DFN-DEM (Discrete Element Method) analysis of six DFN models at the repository level. Increases of permeability up to a factor of four were observed during thermal loading history and shear dilation of fractures was not recovered after cooling of the repository. An understanding of the stress path and potential areas of slip induced shear dilation and related permeability changes during the lifetime of a repository for spent nuclear fuel is of utmost importance for analysing long-term safety. The result of this study will assist in identifying critical areas around a repository where fracture shear slip is likely to develop. The presentation also includes a brief introduction to the ongoing site investigation on two candidate sites for geological repository in Sweden.

• 한국지반공학회논문집
• /
• 제30권3호
• /
• pp.59-72
• /
• 2014

모래지반 표면에 위치한 강체 원형기초를 대상으로 수치해석을 통하여 수직-모멘트 조합하중 조건에서의 지지력을 구하였다. 지반은 Mohr-Coulomb 소성모델을 이용하여 모델링하였으며 관련흐름법칙을 적용하였고, 거친 기초 바닥면 조건에 대하여 검토하였다. 적은 수의 해석으로 조합하중 상관도를 산출할 수 있는 swipe 재하 방법과 통상적인 재하실험에서 적용되는 probe 재하 방법을 적용하여 비교한 결과, 두 방법은 유사한 결과를 나타내었다. 모멘트하중을 고려하기 위하여 전통적으로 사용되는 유효폭 및 유효면적 개념을 사용한 결과와 편심계수($e_{\gamma}$)를 사용한 방법들을 비교하였으며, 기존의 제안식들과 수치모델링으로 구해진 본 연구의 결과를 비교하였다. 수직-모멘트 조합하중 지지력의 내부마찰각에 따른 변화는 미미한 것으로 나타났으며, 유효폭 개념은 편심계수의 형태로 변환하여 원형기초에도 그대로 적용이 가능한 것으로 나타났다. 본 연구의 수치모델링 결과는 기존의 실험에 기반한 결과들에 비해 다소 작은 값을 주는 것으로 나타났으며, 편심 및 모멘트하중이 증가할수록 그 차이는 증가하였다. 수치모델링과 실험 결과가 차이를 나타내는 요인과 향후 연구 방향에 대하여 고찰하였다.

• 한국인터넷방송통신학회논문지
• /
• 제19권1호
• /
• pp.203-209
• /
• 2019

본 연구에서는 WENNER 4전극법 기반에 GMD(신규 대지 고유 저항 측정 장치)와 측정용 Probe(접지동봉)가 PLC(전력선 통신)로 연결된다. 측정용 Probe는 2개(P1,P2)가 1조로 모두 5조 10개의 Probe가 직렬로 각각 1m, 2m, 4m, 8m, 16m 간격으로 대지(토양)에 설치되어 있다. GMD에서 보낸 PLC 신호를 측정용 Probe 1조(P1)의 수신기가 감지하면 Probe에 부착된 PSD(전력 공급 장치)에서 측정용 미세 전압과 전류가 대지로 흐르게 되고 P1과 P2 사이의 토양을 거쳐 Probe 1조(P2)에 유입 된다. 이때 대지 저항으로 인해 전압 강하가 발생되는 원리로 저항값을 측정하게 된다. 이렇게 1~5조까지 T초 간격으로 대지 저항을 측정하고 측정된 데이터는 메인 장비에 탑재된 Arduino Server에 저장 한다. 저장된 측정 데이터는 옴의 법칙(Ohm's Law)에 의한 수식 R=E/I와 고유저항 ${\rho}=6.28aR$ (여기서, R: 측정저항, E: 측정전압, I: 측정전류, a:Probe 간격, ${\rho}$: 고유저항 )를 통해 고유저항을 얻을 수 있다. 실시간으로 얻어진 데이터를 Main PC에 설치된 CDGES 프로그램과 연동되어 데이터 분석이 가능하게 되고 대지(토양)의 접지 환경을 실시간 모니터링 할 수 있게 된다. 또한, 대지(토양)의 온도, 습도 등 계절의 특성을 파악하여 3D 그래프 지원으로 입체적인 Display가 가능하다. 연구의 한계점은 실험적으로 개발 운용한 모델로 상업적인 접근을 위해 Test Bed의 구체적인 적용 방안이 필요할 것이다.

• 한국해안·해양공학회논문집
• /
• 제25권2호
• /
• pp.63-75
• /
• 2013

평택 당진항 내항 2공구 수역은 3개의 투수성 호안으로, 그리고 내항 2공구 동측의 투기장 수역은 2개의 투수성 호안으로 둘러싸여 있다. 2010년 5월에 관측된 내항 2공구 외곽호안 내측 수역과 내항 2공구 투기장 내측 수역의 최대조차는 각각 4.70 m와 2.32 m로서, 동시에 호안 외측에서 관측된 최대조차 8.74 m의 54%와 27%에 달한다. 호안 내 외수위차와 내측 수용적 변화율간의 회귀식을 도출하고, 이 식을 이용하여 투수성 호안의 해수 유통량을 매 계산시간마다 산정하는 모듈을 EFDC 모델에 추가하여 아산만의 3차원 해수유동 수치모형을 구축하였다. 2010년 5월 13~27일의 모의기간에 대하여 주요 5개 분조($M_2$, $S_2$, $K_1$, $O_1$, $N_2$)의 합성조석과 아산, 삽교, 남양, 석문방조제의 담수방류량을 실시간으로 입력하여 해수유동을 모의하였다. 2공구 내측과 2공구 투기장 내측에서 평균고조위, 평균해면과 평균저조위의 실측치에 대한 모델치의 skill score는 96~100%로서 매우 양호한 재현율을 보인다. 투수성 호안의 해수유통을 차단한 모의결과와 비교하면, 최강유속은 주수로를 따라 0.05~0.10 m/s 증가하고, 2공구 외곽호안 외측에서 국지적으로 0.1~0.2 m/s 증가한다. 해저면 전단응력은 유속이 강한 주수로에서 0.1~0.4 $N/m_2$의 범위로 증가하고, 2공구 외곽호안 우각부 주변에서 국지적으로 0.4 $N/m_2$ 이상 증가한다. 본 연구에서 적용한 투수성 호안의 해수유통 모의기법은 대규모 투수성 호안이 유지되는 해역에서 물질의 이류 확산과 해저지형의 침식 퇴적 및 호안 주변의 국부 세굴 등을 모의 예측하는데 유용하게 적용될 수 있다.

• Steel and Composite Structures
• /
• 제46권1호
• /
• pp.53-73
• /
• 2023

Man-made structure materials like concrete usually contain inclusions. These inclusions affect the mechanical properties of concrete. In this investigation, the influence of inclusion length and inclination angle on three-dimensional failure mechanism of concrete under uniaxial compression were performed using experimental test and numerical simulation. Approach of acoustic emission were jointly used to analyze the damage and fracture process. Besides, by combining the stress-strain behavior, quantitative determination of the thresholds of crack stress were done. concrete specimens with dimensions of 120 mm × 150 mm × 100 mm were provided. One and two holes filled by gypsum are incorporated in concrete samples. To build the inclusion, firstly cylinder steel tube was pre-inserting into the concrete and removing them after the initial hardening of the specimen. Secondly, the gypsum was poured into the holes. Tensile strengths of concrete and gypsum were 2.45 MPa and 1.5 MPa, respectively. The angle bertween inclusions and axial loadind ary from 0 to 90 with increases of 30. The length of inclusion vary from 25 mm to 100 mm with increases of 25 mm. Diameter of the hole was 20 mm. Entirely 20 various models were examined under uniaxial test. Simultaneous with experimental tests, numerical simulation (Particle flow code in two dimension) were carried out on the numerical models containing the inclusions. The numerical model were calibrated firstly by experimental outputs and then failure behavior of models containing inclusions have been investigated. The angle bertween inclusions and axial loadind vary from 0 to 90 with increases of 15. The length of inclusion vary from 25 mm to 100 mm with increases of 25 mm. Entirely 32 various models were examined under uniaxial test. Loading rate was 0.05 mm/sec. The results indicated that when inclusion has occupied 100% of sample thickness, two tensile cracks originated from boundaries of sample and spread parallel to the loading direction until being integrated together. When inclusion has occupied 75% of sample thickness, four tensile cracks originated from boundaries of sample and spread parallel to the loading direction until being integrated together. When inclusions have occupied 50% and 25% of sample thickness, four tensile cracks originated from boundaries of sample and spread parallel to the loading direction until being integrated together. Also the inclusion was failed by one tensile crack. The compressive strength of samples decease with the decreases of the inclusions length, and inclusion angle had some effects on that. Failure of concrete is mostly due to the tensile crack. The behavior of crack, was affected by the inclusion length and inclusion number.