TITAN2D를 이용한 제주도에서 발생 가능한 용암돔 붕괴에 의한 화쇄류 수치모의

장철우(Cheolwoo Chang),윤성효(Sung-Hyo Yun) 한국암석학회 2017 암석학회지 Vol.26 No.1

본 연구에서는 제주도에서 화산분화에 의해 발생 가능한 화쇄류의 영향 범위를 파악하기 위하여 백록담 외륜산 외측 사면 8개 지점에서 새로운 용암돔이 발생하는 것을 가정하고 이 과정에서의 붕괴로 발생 가능한 화쇄류를 수치모의 하였다. 수치모의에는 TITAN2D 수치모의 프로그램을 사용하였다. TITAN2D프로그램에서 화쇄류 흐름의 속도를 제어하는 변수로 사용되는 내부마찰각과 층저마찰각 중 내부마찰각을 30°, 층저마찰각을 20°로 설정하였다. 붕괴되는 돔의 높이와 반경, 초기 붕괴 속도, 수치모의 시간 등을 선정하여 총 96개 시나리오에 대하여 수치모의를 실시하였다. 수치모의 결과 돔 붕괴에 의한 화쇄류는 최대 13.4 km까지 도달하여 한라산 정상부 및 제주시와 서귀포시까지 도달하는 것으로 나타났다. 본 연구를 통해 구축된 시뮬레이션 자료들을 이용하여 제주도의 화산 분화 시 발생 가능한 화쇄류의 피해범위 DB를 구축하고 이를 기반으로 재해도를 작성하여 화쇄류로 인한 피해를 최소화하는데 활용할 수 있을 것으로 기대된다. In order to determine the runout range of pyroclastic density currents on Jeju island, lava dome collapse on 8 locations of outer rim of Baekrokdam crater were simulated by TITAN2D numerical simulation program. We set parameters as internal friction angle as 30° and bed friction angle as 20° to control velocity of currents occurred by lava dome collapse. Then we set the height and radius of lava dome, initial speed of collapse and simulation times. And we carried out numerical simulations for a total of 96 scenarios. The result shows that the maximum runout distance was 13.4 km in case of lava dome collapse. This study can be used database for manufacturing of hazard map to minimize damages caused by pyroclastic density currents occurred on Jeju island.

둥지격자기법에 따른 영동지방의 겨울철 강수량에 대한 WRF 모델의 민감도: 사례 연구

이재규,성현민,Lee. Jae Gyoo,Sung. Hyun Min 한국방재학회 2013 한국방재학회논문집 Vol.13 No.1

기후변화로 인해 전 세계적으로 홍수, 폭염, 한파, 가뭄의 피해가 빈번하게 발생하고 있으며, 특히, 강원도 지역은 태백산맥과 동해의 복합적인 영향으로 다른 지역에 비해 재해기상이 빈번히 발생하는 지역으로 2011년 2월에 발생한 대설로 인해 농 어업시설 파괴, 도로 등을 마비시켜 막대한 경제적 피해를 초래한바 있다. 따라서 기상재해를 유발시키는 겨울철 강원 영동지역의 대설사례를 대상으로 둥지격자기법(One-way nesting 기법과 Two-way nesting 기법)을 각각 적용해 WRF 수치모델을 구동하여 그 결과를 분석해 봄으로써, 영동지역의 정확도 높은 대설 예측을 위해 어떠한 둥지격자기법을 사용하여야 하는지 알아보고 자 하였다. 먼저, 영동지역의 대표적인 두 가지 형태의 강수사례, 즉 해안형 강수사례(event A)와 산악형 강수사례(event B)를 각각 선정하여, 지표바람장과 발산장 그리고 누적 강수량의 분포를 분석하였다. 그리고 강수 예측의 정확성을 알아보고자 oneway nesting 기법과 two-way nesting 기법을 사용하여 모의된 3시간 누적 강수량에 대한 Threat Score를 계산하여 분석하였다. 먼저, one-way nesting 기법과 two-way nesting 기법의 사용에 따른 수치모의 결과의 차이를 상호 비교해보면, 해안형 강수사례(event A)의 경우, two-way nesting 기법을 사용한 모의가 one-way nesting 기법을 사용한 모의보다 영동 해안지역에서 서풍 계열의 바람이 부는 영역을 보다 넓게, 그리고 서풍 계열의 바람을 다소 강하게 모의하여 one-way nesting 기법을 사용한 모의보다 two-way nesting 기법을 사용한 모의가 해상쪽(동쪽)으로 갈수록 누적 강수량을 많이 모의하였다. 결과적으로 two-way nesting 기법이 one-way nesting 기법보다 영동지역에 모의되는 강수분포를 다소 동쪽으로 치우치는 경향을 보였다. 산악형 강수사례(event B)의 경우, two-way nesting 기법이 one-way nesting 기법보다 산악지역에서 동풍기류를 상대적으로 약하게 모의하였고 이에 따라, 산악지형에 의한 강제 상승효과가 상대적으로 작았기 때문에 산악지역에서 two-way nesting 기법을 이용하여 모의한 누적 강수량이 one-way nesting 기법을 사용하여 모의한 누적 강수량보다 적게 나타났다. 한편, 수치모의 결과에 대한 상호비교가 아닌 실제 관측과 비교하기 위하여, event A와 event B사례에 대하여 one-way nesting 기법과 two-way nesting 기법을 각각 적용하여 얻어진 수치모의 결과자료와 AWS 관측자료를 비교한 Threat score를 구하였다. 비록 one-way nesting 기법을 사용한 경우보다 two-way nesting 기법을 사용한 경우, 강수역이 보다 해상쪽(동쪽)으로 치우쳐 모의되는 경향을 보였지만, 두 기법에 대한 TS값의 차이(평균 약 0.01 정도)는 매우 작게 나타나, 어느 기법을 사용하든 강수 예측의 정확성은 거의 비슷한 수준이라고 볼 수 있었다. WRF (the Weather Research and Forecasting) model was run with one-way grid nesting or two-way grid nesting, respectively for Yeongdong province`s two typical precipitation cases of the coastal precipitation case (event A) and the mountain precipitation case (event B) which were synoptically under the influence of the northeasterly to examine whether the choice of nesting method can have a significant impact on the precipitation distribution. For the coastal precipitation case (event A), the average accumulated-precipitation from simulation using two-way nesting was approximately 7.78 mm bigger than that using one-way nesting in the Yeongdong coast area. However, that from simulation using two-way nesting was approximately 4.23 mm smaller than that using one-way nesting in the mountain area. Similarly, for the mountain precipitation case (event B), the average accumulated-precipitation from simulation using two-way nesting was approximately 1.56 mm bigger than that using one-way nesting in the Yeongdong coast area. In addition, that from simulation using two-way nesting was approximately 4.91 mm smaller than that using one-way nesting in the mountain area. That is, for both the coastal precipitation case (event A) and the mountain precipitation case (event B), the average accumulated-precipitation from simulation using two-way nesting was bigger than that using one-way nesting in the coast area, but in the mountain area that was the reverse.

강원 연안 통합수치모의시스템 구축 연구

이정렬,김정원,김인호 한국방재학회 2011 한국방재학회 학술발표대회논문집 Vol.10 No.-

강원 동해안은 해안침식 및 개발 후유증으로 2007년 이후 모니터링 체계가 차근히 구축되고 있으며 연안관리체계도에 따라 각 해안의 표사계를 정립하고 임시기준점을 설정하여 연구조사 및 모니터링 단계의 기초를 마련하고 있다. 통합수치모의시스템은 모니터링 시스템과 수치모의 시스템 간의 상호 결합 관계를 나타낸 것으로 수치모의시스템에 지속적인 변화 자료의 제공과 시스템 검증, 그리고 시스템으로부터 산출된 재생산 자료의 장기간 축적 과정이 도시된다. 표사계로부터 모의 시스템의 하천토사유입량, 연안토사유출입량 등이 제공되며 임시기준점, 해안선 및 해빈단면 측량 자료로부터 해저면 자료가 지속적으로 반영된다. 통합수치모의시스템은 재난대응을 위한 실시간 예보, 원인규명을 위한 진단, 개발계획의 영향평가를 목적으로 광역 또는 모니터링 자료의 축적에 따른 Database로부터 중간역, 상세 연안역을 거쳐 좀 더 상세한 연안 현상이 모의되도록 하며 본 논문에서는 일례로 태풍 TALAS의 동해 진입에 따른 상세역 수치모의를 GW10표사계인 강원도 문암-교암 일대에 적용하였다.

3차원 수치모의를 이용한 보 접속부 세굴 특성 분석

정석일,여창건,윤광석,이승오,Jung,Suk Il,Yeo,Chang Geon,Yoon,Gwang Seok,Lee,Seung Oh 한국방재학회 2011 한국방재학회 학술발표대회논문집 Vol.2011 No.1

보는 하천에서 취수나 하상유지를 위한 하천횡단구조물로 일반적으로 본체, 물받이, 바닥보호공 등으로 구성되며 제방 연결부 호안 및 밑다짐 등에 대한 설계기준이 하천설계기준에 제시되어 있다. 그러나 태풍 루사나 매미에 의한 피해사례를 보면 하천 횡단 구조물 본체가 파괴되는 피해 뿐만 아니라, 구조물과 제방과의 연결부가 세굴이 발생되어 붕괴되는 사례가 많이 발생하고 있다. 하천설계기준 해설(2009)에는 이러한 보와 제방의 연결부 부분을 연결호안이라하여 관련 기준을 제시하고 있으나, 설치구간의 길이를 정할 때 하천의 규모나 하도의 특성을 고려하지 못하고 일률적으로 결정하도록 하고 있다. 이에 건설기술연구원의 `보 및 낙차공 설계기술 개발 연구보고서`(윤광석 등, 2006)에서는 고정상 실험을 통해 연결호안 설치구간에 대한 실험식을 제시하였다. 본 연구에서는 3차원 수치모의를 통하여 건설기술연구원에서 수행하였던 실험을 재현하고, 제방을 이동상으로 하여 인자들 간의 상호 관계를 밝히며, 그 특성을 분석하고자 한다. 그리고 윤광석 등(2006)에 의해 제시된 실험식에 적용하여 검증하였다. 수치모의는 유량을 <TEX>$0.7{\sim}2.8m^3/sec$</TEX>까지 변화하며 수행하였으며, 유사의 대표입경은 0.63mm로, 상류수심은 1.0m로 일정하게 유지하였다. 수치모의는 평형세굴 발생 후, 최소 모의시간의 10%정도 지난 시간까지 하였다. 수치모의 결과는 다음과 같다. (1) 유량이 증가함에 따라 유속 및 Froude 수가 증가하여 상 하류부 세굴 발생 범위와 폭은 증가하였다. (2) 상류는 하류에 비해 유량에 따른 세굴발생 범위가 상대적으로 작게 나타났으며, 이는 하류단에 비해 상류단의 유속 및 Froude 수의 차이가 작았기 때문인 것으로 판단된다. (3) 세굴의 폭을 측정함으로써 세굴에 가장 취약한 부분을 짐작할 수 있으며, 설계에 반영되어 호안이나 옹벽의 두께결정에 적용한다면 세굴에 대해 좀 더 안전한 설계가 될 수 있을 것이다. (4) 건설기술연구원(2006)에서 제시한 식(1)과의 비교를 통해 수치모의 결과가 식(1)로부터 계산된 값보다 작음을 알 수 있으며, 그 이유는 식(1)의 범위는 와류영역구간을 나타내서 연결호안 설치구간 길이를 제시하고 있는 반면, 본 연구의 수치모의 결과는 세굴이 발생한 범위를 제시하였기 때문이다. 향후 보 높이와 좀 더 다양한 유량에 대한 경우를 수치모의하고, 이동상 제방에 대한 실험을 통해 명확한 식을 제안 할 것이다.

잠긴수제가 설치된 만곡수로에서의 이차류 거동 수치모의

이정섭,박상덕,최철희,백중철 한국수자원학회 2019 한국수자원학회논문집 Vol.52 No.10

The flow in the meandering channel is characterized by the spiral motion of secondary currents that typically cause the erosion along the outer bank. Hydraulic structures, such as spur dike and groyne, are commonly installed on the channel bottom near the outer bank to mitigate the strength of secondary currents. This study is to investigate the effects of submerged vanes installed in a 90° meandering channel on the development of secondary currents through three-dimensional numerical modeling using the hybrid RANS/LES method for turbulence and the volume of fluid method, based on OpenFOAM open source toolbox, for capturing the free surface at the Froude number of 0.43. We employ the second-order-accurate finite volume methods in the space and time for the numerical modeling and compare numerical results with experimental measurements for evaluating the numerical predictions. Numerical results show that the present simulations well reproduce the experimental measurements, in terms of the time-averaged streamwise velocity and secondary velocity vector fields in the bend with submerged vanes. The computed flow fields reveal that the streamwise velocity near the bed along the outer bank at the end section of bend dramatically decrease by one third of mean velocity after the installation of vanes, which support that submerged vanes mitigate the strength of primary secondary flow and are helpful for the channel stability along the outer bank. The flow between the top of vanes and the free surface accelerates and the maximum velocity of free surface flow near the flow impingement along the outer bank increases about 20% due to the installation of submerged vanes. Numerical solutions show the formations of the horseshoe vortices at the front of vanes and the lee wakes behind the vanes, which are responsible for strong local scour around vanes. Additional study on the shapes and arrangement of vanes is required for mitigate the local scour. 만곡수로에서의 흐름은 나선형 운동 형태의 이차류가 지배적이며, 이로 인해 일반적으로 만곡 외측을 따라 침식 현상이 발생하게 된다. 이러한 이차류를 약화시키기 위해서 보통 만곡수로 외측을 따라서 수제와 같은 수공구조물을 설치한다. 이 연구에서는 OpenFOAM 오프소스 소프트웨어를 토대로 난류 해석을 위한 하이브리드 RANS/LES 기법과 자유수면 해석을 위한 VoF기법을 이용한 3차원 수치모의를 통해서 90° 만곡수로에 설치된 잠긴수제가 후루드수가 0.43인 조건에서 이차류의 발달에 미치는 영향을 분석하였다. 시간과 공간에 대해서 2차 정확도의 유한체적법을 이용하여 수치모의를 수행하였으며, 수치해석 결과는 실험결과와 비교하여 수치모의의 정확도를 평가하였다. 잠긴수제가 설치된 경우의 수치모의 결과를 흐름방향 유속 분포와 횡방향 순환 유속벡터장을 중심으로 수리실험 관측값들과 비교할 때 수치모의 결과는 수리실험에서 관측된 주요 이차류 흐름 거동과 현상 들을 대부분 양호한 정확도로 잘 재현하는 것으로 나타났다. 수치모의 결과를 비교해보면, 잠긴수제 설치로 인해서 만곡이 끝나는 단면 외측 하상부근에서의 유속은 약 평균유속의 1/3 정도 감소하는 반면에 수제 상단부에서의 전단층 발달에 따른 흐름 가속으로 자유수면 부근까지 유속이 증가하고 만곡 수충부에서는 수면 부근 유속이 약 20% 증가하는 것으로 나타났다. 결과적으로 잠긴수제는 만곡부에서 발생하는 이차류의 강도를 약화시켜 만곡부 외측 하상의 안정에 도움이 될 것으로 판단된다. 한편, 각 잠긴수제 전면부에서 말발굽와가 그리고 후면부에서는 후류가 형성되면서 수제 구조물 주변에서 강한 국부세굴이 발생하는 것으로 나타남으로, 국부세굴을 최소화할 수 있는 수제의 형상 및 배열에 대한 추가 연구가 요구된다.

수치모형을 이용한 보의 개방구성에 따른 흐름모의 실험

강태운(Tae Un Kang),장창래(Chang-Lae Jang) 응용생태공학회 2020 Ecology and resilient infrastructure Vol.7 No.3

본 연구는 수치모형을 기반으로 하는 예측방법론을 위한 수치기법개발의 첫 단계로서 제방형 보에 의한 자유형 월류흐름에 대한 수치해석실험을 수행하였다. 이를 위해 2차원 흐름모형인 Nays2DH를 이용하여 기존에 수행된 실험자료를 바탕으로 수치실험을 수행하여 모형을 비교검증하였다. 우선은 위어의 형태에 의한 불연속 흐름의 수치적 재현성에 대해서 고찰했으며 하도의 전처리를 통해 모의흐름을 보정하였다. 더 나아가, 수문갯수에 따른 월류흐름의 예측모의를 수행하여 결과를 비교하였다. 모의결과, 수문의 갯수가 증가할 경우 동일한 통수단면적에서 보 하류부의 최대유속이 증가할 수 있다는 것을 알 수 있었으며 이는 하상변동에도 영향을 줄 것으로 판단된다. 이러한 결과들을 통해, 본 연구는 추후의 보의 운영에 따른 흐름 및 유사이송 등을 고려하는 수리학적 연구의 기초자료를 제공할 수 있었다. This study investigated that the numerical experiment for analysis on free overtopping flow by a weir of levee type, as the first stage of the development of a numerical technique for prediction methodology based on a numerical model. Using 2-dimensional flow models, Nays2DH, we conducted numerical simulations based on existing experimental data to compare and verify the models. We firstly discussed the numerical reproducibility for the discontinued flow by weir shape, and calibrated the computational flow through preprocessing of channel bed. Further, we carried out and compared the simulations for prediction on the overtopping flow by the number of weir gates. As a result of simulations, we found that the maximum flow velocity of downstream of weir increases when the number of weir gates increases under the same cross sectional area of flow. Through such results, this study could present basic data for hydraulic research to consider the water flow and sediment transport depending on weir operation in the future work.

우리나라 영향 태풍에 따른 호우와 강풍에 대한 수치예보모델 모의 성능 비교 및 검증

김은진,박병권,안용준 한국기상학회 2021 한국기상학회 학술대회 논문집 Vol.2021 No.4

북서 태평양 적도 부근에서 주기적으로 발생하는 태풍은 여름철 우리나라를 내습하여 큰 피해를 야기시키고 있다. 특히, 태풍은 우리나라 한 지역에서 짧은 시간 동안 많은 양의 강우와 강한 바람을 발생시킨다. 이로 인한 피해규모는 우리나라의 자연재해 피해규모 중 약 54%를 차지한다. 지금까지의 선행연구는 주로 태풍 이동 경로에 따른 강우 분포와 재해 특성에 관한 연구가 수행되었으나 태풍으로 인한 우리나라에서의 지역별 최대 위험시점과 그 시점의 강수량, 풍속에 관한 연구는 미미하다. 이는 태풍 내습 시 발생할 수 있는 태풍의 피해를 대비하기 위한 중요한 정보가 될 수 있다. 따라서 본 연구에서는 기상청에서 현업 운영중인 통합모델(UM) 기반의 전지구 수치예보시스템(GDPS)과 유럽중기예보센터(ECMWF)에서 생산하는 전지구 수치모델, 그리고 한국형수치예보모델개발사업단에서 개발한 한국형수치예보모델(KIM)을 이용하여 우리나라에 영향을 미치는 태풍에 동반하는 호우와 강풍의 모의 성능과 최대 위험시점을 비교 분석한다. 또한, 연구에 사용한 수치예보모델들에 대하여 태풍 기간 내 누적 강수량, 풍속의 RMSE, MAE, BIAS 검증을 실시한다. 분석을 위해 5km 해상도로 내삽된 모델 격자자료를 사용하였으며, 718개 자동방재기상관측장비(AWS)에서 수집한 관측자료는 객관분석을 통해 격자화하여 사용하였다. 본 실험을 통하여 수치예보모델 모의 경향을 이해하고, 우리나라를 내습하는 태풍으로 인한 호우와 강풍 피해를 예방할 수 있도록 방재유관기관 업무 의사결정에 배경자료로써 활용할 수 있을 것이다.

수치모의를 통한 해운대 이안류에 대한 해상요소의 영향 연구: 파주기, 파고, 파향, 조위

최준우,신충훈,윤성범,Choi, Junwoo,Shin, Choong Hun,Yoon, Sung Bum 한국수자원학회 2013 한국수자원학회논문집 Vol.46 No.2

Boussinesq 모형인 FUNWAVE를 사용한 수치모의 결과를 바탕으로, 파고, 주기, 파향 및 조위와 같은 파랑특성 변수에 따라 해운대 이안류의 발생가능 정도를 추정하였다. 2011년 6월 12일에 발생한 이안류를 수치모의하여 이안류 발생 가능정도를 확인함으로써 정량화 기법을 제시한다. 이 이안류 발생정도의 정량화 기법은 다음의 절차를 따른다. 첫째로, 광범위하게 파고, 주기, 파향, 조위를 변화시키며 불규칙파를 해운대 수심지형에 적용하여 해운대 연안흐름을 수치모의 하였다. 둘째로, 이 수치모의 결과로부터 연안의 직각방향(cross shore direction)의 2주기 시간평균된 유속성분을 추출하고, 연안영역에 대하여 이 유속값의 심해방향 최대값을 그 영역의 '이안류 유속'으로 정의하였다. 셋째로, 전체 수치모의 시간에 대한 '이안류 유속' 시계열이 임의의 위험유속을 초과한 시간의 비를 위험 이안류 발생 가능 정도로 정의하여 정량화하였다. 추정된 이안류 발생정도에 따른 분석으로부터, 파고가 클수록, 주기가 길수록, 조위가 낮을수록 이안류가 더 잘 발생하는 것으로 나타났으며, 파향은 남동쪽보다 남서쪽 파향에서 조금 더 이안류가 잘 발생하는 것으로 나타났다. The likelihood of rip current at Haeundae beach according to wave parameters, such as wave height, period, direction, and tidal elevation, was estimated by using numerical simulations with a Boussinesq model, FUNWAVE. To examine the estimation, the rip current occurred on 12th June, 2011 at Haeundae beach was simulated based on observations. For the estimation, the following procedure was carried out. First, extensive numerical simulations of nearshore circulations are performed under various random sea conditions according to the wave parameters. Second, from the numerical results, cross shore components of two-wave-period averaged velocities over the nearshore area were computed, and their seawardly maximum was defined as rip current velocity of the area. Third, using time series of the rip current velocity, we computed the ratio of the simulation time and the time period in which the rip current velocity exceed a threshold velocity for rip-current accidents, and thus the ratio was quantified as the likelihood of rip current at Haeundae beach for the input wave parameters. From the resultant estimations, it was found that the rip current likelihood increases as wave height and period increase, and tidal elevation decreases.

수치실험을 통한 제방붕괴부 내 수리학적 특성 검토

양지로,김경모,김수영,이승오 한국방재학회 2014 한국방재학회 학술발표대회논문집 Vol.2014 No.-

기존 수행되는 제방붕괴관련 연구는 댐붕괴 형태로 가정하여 해석하기 때문에 하도 내의 흐름을 반영하지 못하였다. 이러한 단점을 보완하기 위해 하도 내의 흐름을 반영할 수 있는 수리모형실험이 진행되었다. 하지만 붕괴부 내의 흐름을 파악하기에는 붕괴부가 발달함에 따라 측점이 이동하기 때문에 계측장비를 이용한 측정이 어려웠기 때문에 수리모형실험에서는 붕괴부의 유량이나 폭에 대한 연구가 대부분이였다. 따라서 본 연구에서는 하도흐름과 붕괴부 형상을 수치모의를 통해 재현한 후 붕괴부의 흐름특성을 분석하고자 하였다. 수치모의는 기존에 수행된 수리모형실험과 동일한 조건으로 수행하였으며 제방의 제원은 제방고 0.3 m 제방경사 1:2(H:L) 마루폭 0.3 m이다. 수리모형실험 결과를 이용하여 시간에 따른 붕괴형상을 FLOW-3D의 GMO(General Moving Object)의 기능을 이용하여 구현하였다. 수치모의와 수리모형실험에 대해서 붕괴유량과 하도수위를 비교하여 수치모의가 수리모형실험을 유사하게 재현함을 확인하였고 수치모의를 통해 도출된 유속분포를 이용하여 제체에 가해지는 외력을 분석함으로써 제방 붕괴의 발달 매커니즘을 분석할 수 있었다. 본 연구에서 도출된 결과를 이용하여 제방붕괴를 발달시키는 주요변수를 파악하면 제방붕괴에 대한 이론을 확립하는데 중요한 자료로 활용될 수 있을 것으로 판단된다. 추후 붕괴부내 흐름뿐만 아니라 홍수터를 전파하는 흐름까지 재현하여 보다 정확한 침수면적을 예측하여 홍수 피해를 감소시키는 데 주요한 역할을 할 수 있을 것이다.

수치모의를 통한 해운대 이안류에 대한 해상요소의 영향 연구: 파주기, 파고, 파향, 조위

최준우,신충훈,윤성범 한국수자원학회 2013 한국수자원학회논문집 Vol.46 No.2

Boussinesq 모형인 FUNWAVE를 사용한 수치모의 결과를 바탕으로, 파고, 주기, 파향 및 조위와 같은 파랑특성 변수에 따라 해운대 이안류의 발생가능 정도를 추정하였다. 2011년 6월 12일에 발생한 이안류를 수치모의하여 이안류 발생 가능정도를 확인함으로써 정량화 기법을 제시한다. 이 이안류 발생정도의 정량화 기법은 다음의 절차를 따른다. 첫째로, 광범위하게 파고, 주기, 파향, 조위를 변화시키며 불규칙파를 해운대 수심지형에 적용하여 해운대 연안흐름을 수치모의 하였다. 둘째로, 이 수치모의 결과로부터 연안의 직각방향(cross shore direction)의 2주기 시간평균된 유속성분을 추출하고, 연안영역에 대하여 이 유속값의 심해방향 최대값을 그 영역의 ‘이안류 유속’으로 정의하였다. 셋째로, 전체 수치모의 시간에 대한 ‘이안류 유속’ 시계열이 임의의 위험유속을 초과한 시간의 비를 위험 이안류 발생 가능 정도로 정의하여 정량화하였다. 추정된 이안류 발생정도에 따른 분석으로부터, 파고가 클수록, 주기가 길수록, 조위가 낮을수록 이안류가 더 잘 발생하는 것으로 나타났으며, 파향은 남동쪽보다 남서쪽 파향에서 조금 더 이안류가 잘 발생하는 것으로 나타났다. The likelihood of rip current at Haeundae beach according to wave parameters, such as wave height, period, direction, and tidal elevation, was estimated by using numerical simulations with a Boussinesq model, FUNWAVE. To examine the estimation, the rip current occurred on 12th June, 2011 at Haeundae beach was simulated based on observations. For the estimation, the following procedure was carried out. First, extensive numerical simulations of nearshore circulations are performed under various random sea conditions according to the wave parameters. Second, from the numerical results, cross shore components of two-wave-period averaged velocities over the nearshore area were computed, and their seawardly maximum was defined as rip current velocity of the area. Third, using time series of the rip current velocity, we computed the ratio of the simulation time and the time period in which the rip current velocity exceed a threshold velocity for rip-current accidents, and thus the ratio was quantified as the likelihood of rip current at Haeundae beach for the input wave parameters. From the resultant estimations, it was found that the rip current likelihood increases as wave height and period increase, and tidal elevation decreases.