백두산 화산분출물의 이류와 확산 예측 특성을 살펴보기 위하여 3차원 대기역학 모형인 WRF와 입자 확산 모형인 FLEXPART를 결합하여 다양한 수치실험을 실시하였다. 기상자료의 시간 평균화에 따른 영향을 보기 위하여 네가지 서로 다른 평균화 기간을 가진 수치 실험을 실시하였다. 본 연구에 적용된 기상자료의 평균화 기간은 각각 1개월, 10일, 3일이다. 또한 평균화를 실시하지 않은 실시간자료를 이용한 수치실험도 실시하였다. 기상자료의 평균화 시간이 길어질수록 주풍성분인 동서 방향의 운동이 뚜렷해지고, 짧을수록 주풍의 법선 방향인 남북 운동이 명확하게 나타나며, 화산분출물을 가정한 입자의 확산 역시 동일한 특성이 나타난다. 상하층의 바람은 강도의 차이가 명확하고, 평준화 기간에 따른 영향이 다르게 나타나기 때문에 상층입자와 하층입자의 이동과 지상 침적이 다양하게 나타난다. 또한 한반도의 지상 침적은 방출고도 2 km 이하의 입자가 주로 영향을 미친다. 따라서 기상자료의 평균화 시간 간격이 화산분출물의 지상 침적을 결정하는 하나의 요인으로 작용하기 때문에, 확산 실험 전에 적절한 기상자료의 시간평균을 결정할 필요가 있다.

Stokes drift(SD) and Lagrangian discharge(LD) are important factors for analysis of flushing time, tidal exchange, solute transport and pollutant dispersion. The factors should be calculated using the approached method to flow phenomena. The aim of this paper re-examines the previous procedures for computing the SD and LD, and is to propose the new method approached to stratified flow field in the cross-section of coastal region, e.g. Masan Bay. The intensity of velocity near the bottom boundary layer(BBL) depends on the sea-bed irregularity in the coastal estuaries. So we calculated the depth mean velocity(DMV) considering that of BBL omitted in Kjerfve's calculation method. It revealed that BBL effect resulting in application of the bay acts largely on DMV in half more among 1l stations. The new expression of SD and LD per unit width in the cross-section using the developed DMV and proposed decomposition procedure of current were derived as follow : Q=u0+½H1U1cos(Φh-Φu)+U3cos(Φh-Φud) LD ED SD(Qskim+Qsk2) The third term, Qsk2, on the right-hand of the equation is showed newly and arise from vertical oscillatory shear. According to the results applied in 3 cross-sections including 11 stations of the bay, the volume difference between proposed and previous SD was founded to be almost 2 times more at some stations. But their mean transport volumes over all stations are 18% less than the previous SD. Among two terms of SD, the flux of second term, Qskim, is larger than third term, Qsk2, in the main channel of cross-section, so that Qskim has a strong dependence on the tidal pumping, whereas third term is larger than second in the marginal channel. It means that Qsk2 has trapping or shear effect more than tidal pumping phenomena. Maximum range of the fluctuation in LD is 40% as compared with the previous equations, but mean range of it is showed 11% at all stations, namely, small change. It mean that two components of SD interact as compensating flow. Therefore, the computation of SD and LD depend on decomposition procedure of velocity component in obtaining the volume transport of temporal and spacial flow through channels. The calculation of SD and LD proposed here can separate the shear effect from the previous SD component, so can be applied to non-uniform flow condition of cross-section, namely, baroclinic flow field.

본 연구에서는 3차원 이산 균열망 수치모형을 이용하여 균열망을 구성하는 균열요소의 길이분포가 유체 흐름 특성이 미치는 영향에 대해 수치적으 로 분석하였다. 균열요소의 길이분포의 생성을 위해 절단멱분포법칙을 적용하였으며, 지수 β₁을 1.0에서부터 6.0까지 변화시키면서 유체 흐름 모의를 수행하였다. 모의결과 지수 β₁이 증가함에 따라 균열요소들의 길이분포는 점차적으로 작아지며, 이로 인해 균열망의 투수성에 영향을 미치는 균열요소들 간의 연결성은 취약해지는 것으로 나타났다. 각각의 지수 β₁에 대해 균열요소 각각에서 계산된 유량분포를 분석하였을 때 β₁= 1.0에 서의 평균유량이   6.0에 비해 약 447배 크게 산정되었으며, 균열망의 유출경계에서 계산된 유량의 경우 β₁일 때가 6.0에 비해 약 6,440배 크게 산정되었다.