• Journal of Korean Society for Atmospheric Environment
• /
• v.8 no.3
• /
• pp.162-168
• /
• 1992

Numerical solutions to the advection equations used for long-range transport air pollution models are calculated using three numerical methods; Antidiffusion correction method(Smolarkiewicz, 1983), Positive definite advecton scheme obtained by nonlinear renormalization of the advective fluxes(Bott, 1989), and Positive definite pseudospectral method(Bartnicki, 1989). Accuracy, numerical diffusion and computational time requirement are compared for two-dimensional transport calculations in a uniform rotational flow field. The solutions from three methods are positive definite. Bartnicki(1989)'s method is most conservative but requires approximately 10 times as much computational time as Smolarkiewicz(1983)'s method of which numerical diffusion is the largest. All three methods are more conservative for a cone shape initial condition than for a rectangular block initial condition with a steep gradient.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.1
• /
• pp.29-41
• /
• 1996

The two-stage numerical model was used to study the relation between three-dimensional local wind model, advection/diffusion model of random walk method and second moment method on Pusan coastal area. The first stage is three dimensional time-dependent local wind model which gives the wind field and vertical dirrusion coefficient. The second stage is advection/diffusion model which uses the results of the first stage as input data. First, wind fields on Pusan coastal area for none synoptic scale wind showed typical land and sea breeze circulation, and convergence zone occured at 1200LST in northern of domain, in succession, moved northward of domain. Emissions from Sinpyeong industrial district were trasnported toward the inland by sea breeze during daytime, and reached the end part of domain about 1800LST. During nighttime, emissions return to sea by land breeze and vertical diffusion also contributes to upward transport. In order to use this model for forecast of air pollution concentration on the Pusan coastal area, it is necessary that computed value must be compared with measured value and wind fields model must also be dealt in detail.

• Proceedings of the Korea Air Pollution Research Association Conference
• /
• 2003.11a
• /
• pp.431-432
• /
• 2003

최근에 와서 연안에 위치한 공단지역과 대도시지역에서 배출된 대기오염물질들의 이류 및 확산에 대한 연구들이 많이 이루어짐에 따라 중규모 기상장 모델링에 관한 연구도 함께 이루어지고 있다(Pielke and Uliasz, 1998). 이런 지역의 오염도는 배출원과 오염물질의 종류 및 특성에 따라 좌우되지만, 그때의 기상 조건과 지형에도 큰 영향을 받는다. 그러므로 대기 오염 농도를 정확히 예측하기 위해서는 대기 오염 물질의 수송, 확산 및 침적에 큰 영향을 미치는 대기 유동장을 표현할 수 있는 수치 모델이 필요하다(오은주, 1995). (중략)

• Proceedings of the Korean Geotechical Society Conference
• /
• 2008.10a
• /
• pp.1512-1517
• /
• 2008

Laboratory column tests were conducted in this study using $Cl^-$ tracers on Jumunjin sand to analyze contaminant transport in mixed contaminated soils. Results obtained from clean soils and soils containing residual diesel verified heterogeneous distribution of residual diesel, and clear acceleration of solute movement. In addition, asymmetric breakthrough curves indicated development of immobile region where solute movement becomes stagnant and creates tailing phenomenon.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.1
• /
• pp.62-72
• /
• 2008

In this study, an analytical model has been developed to predict the build-up of heat field due to a point heat source in the presence of sharp cross-flow depth change. The model has been applied to investigate the effect of the depth change and flow pattern on the heat field. Model results show that, when there is a sharp depth change in cross-flow direction, the heat transport across the boundary of the depth change is enhanced or diminished according to the increasing or decreasing of the horizontal diffusion flux. Including residual components as well as tidal currents give rise to reduce the effect of the horizontal diffusion on the heat transport because of increasing the advection of heat.

• Journal of Korean Society of Environmental Engineers
• /
• v.27 no.10
• /
• pp.1043-1051
• /
• 2005

In this stduy, adsorption and transport characteristics of Fe in the soil were investigated using convection-dispersion local equilibrium sorption model and two-site non-equilibrium sorption model. In batch experiments with different Fe concentration, characteristics of Fe adsorption was investigated using Freundlich and linear isotherm. Column experiments with different flow rate, organic matter content md Fe concentration were also carried out. We measured Fe concentrations in injection-liquid and in effluent, and then applied them to CXTFIT program. As a result of column experiments, some parameters(D, R, ${\beta}$, ${\omega}$) used in two-site non-equilibrium adsorption model were obtained. Characteristics of Fe transport were analyzed using the parameters(D, R, ${\beta}$, ${\omega}$) obtained from the CXTFIT program, Consequently, characteristics of Fe transport in the soil were predicted through two-site non-equilibrium adsorption model.

• Nuclear Engineering and Technology
• /
• v.25 no.1
• /
• pp.154-165
• /
• 1993

A stochastic method using continuous time Markov process is presented to model the one-dimensional convective nuclide transport in geologic media, which have usually heterogeneous feature in physical/geochemical parameters such as velocity, dispersion coefficient, and retardation factor resulting poor description by conventional deterministic advection-dispersion model. The primary desired quantities from a stochastic model are the mean values and variance of the state variables as a function of time. The time-dependent probability distributions of nuclides are presented for each discretized compartment given the volumetric groundwater flux and the intensity of transition. Since this model is discrete in medium space, physical/geochemical parameters which affect nuclide transport can be easily incorporated for the heterogeneous media as well as remarkably layered media having spatially varied parameters. Even though the Markov process model developed in this study was shown to be sensitive to the number of discretized compartments showing numerical dispersion as the number of compartments are increased, this could be easily calibrated by comparing with the analytical deterministic model.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.1
• /
• pp.1-10
• /
• 2020

The application of multiphase flows is increasingly being applied to analyze phenomena such as single phase flows where the fluid boundary changes continuously over time or the problem of mixing a liquid phase and a gas phase. In particular, multiphase flow models that take into account incompressible Newtonian fluids for liquid and gas are often applied to solve the problems of the free water surface such as wave fields. In general, multi-phase flow models require time-based the surface tracking of each fluid's phase boundary, which determines the accuracy of the final calculation of the model. This study evaluates the advection performance of representative VOF-type boundary tracking techniques applied to various CFD numerical codes. The effectiveness of the FCT method to control the numerical flux to minimize the numerical diffusion in the conventional VOF-type boundary tracking method and advection calculation was mainly evaluated. In addition, the possibility of tracking performance of free surface using CIP method (Yabe and Aoki, 1991) was also investigated. Numerical results show that the FCT-VOF method introducing an anti-diffusive flux to precent excessive diffusion is superior to other methods under the confined conditions in this study. The results from this study are expected to be used as an important basic data in selecting free surface tracking techniques applied to various numerical codes.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.105-105
• /
• 2012

최근 우리나라는 집중호우로 인한 토석류의 발생이 현저하게 증가하고 있으며, 2002년 태풍 루사와 2003년 매미 그리고 2006년 7월 집중호우 등으로 인한 피해의 규모는 여러 조사와 문헌에서도 확인할 수 있다. 2011년 7월, 서울 우면산 일대에서는 집중호우로 인한 토석류가 발생하였으며, 16명의 인명손실을 포함한 큰 피해가 있었다. 우면산 토석류는 십여 개 지역에서 동시다발적으로 발생하였으며, 토석류 발생 유역 특성과 토석류의 유동특성을 분석하기 위한 현장 조사가 진행되었다. 조사한 자료에 따르면 우면산 일대의 일일 강우량은 서초 기상측정소를 기준으로 최대 24시간 누적 324mm 그리고 시간당 최대 68.5mm/hr를 기록하였다. 상업용 소프트웨어인 FLO-2D는 유사농도의 함수로서 점성(viscous)응력, 항복(yield)응력, 난류 및 분산(dispersive) 응력항을 포함하는 2차 유변학(quadratic rheology) 모델을 기본으로 사용하여 이류(mudflow)와 토석류(debris flow)를 모의할 수 있다. FLO-2D는 흐름의 운동량 및 에너지 보존을 고려하여 격자와 시간에 관계없이 유동심도, 속도, 압력을 예측할 수 있으며, 격자 기반의 모델로서 GIS 및 기타 응용 프로그램들과 연동이 쉽다는 장점이 있다. 그러나 하상침식에 의해 유발된 토석류의 체적 증가는 고려 할 수 없으므로 토석류의 전파 및 퇴적영역에서의 토석류 모의에만 사용할 수 있는 단점이 있다. 이 연구의 목적은 FLO-2D 소프트웨어를 이용하여 우면산에서의 토석류 현상을 재현하는 것이다. 우면산 일대에서 발생한 토석류 중 서초구 방배동의 래미안 아파트 부근에서 발생한 토석류에 대하여 수치지도(DEM)와 현장조사를 통해 얻은 지형자료, 해당 지역의 강우량 및 지질 특성 자료 등을 토대로 FLO-2D 모델을 적용하여 토석류의 흐름특성을 검토한다. 토석류 유동 및 퇴적에 대한 가용한 현장관측 자료와의 비교 분석을 통하여 토석류 특성 값을 산정하고, 모델의 적용성을 검증한다.

• Economic and Environmental Geology
• /
• v.54 no.4
• /
• pp.465-473
• /
• 2021

This numerical study aims at analyzing the effect of flow characteristics, caused by geometrical features such as intersecting angles, on solute mixing at fracture junctions. It showed that not only Pe, the ratio of advection to diffusion, but also the intersecting angles played an important role in solute mixing at the junction. For the intersection angles less than 90°, the fluid flowed to the outlet in the same direction as the injected flow direction, which increased the contact at the junction with the streamlines coming from the different inlets. On the other hand, for the intersecting angles greater than 90°, the fluid flowed out to the outlet opposite to the flow direction in the inlet, leading to minimizing the contact at the junction. Therefore, in the former case, solute mixing occurred even at high Pe, and in the latter case, solutes transport along the streamlines even at low Pe. For Pe < 1, the complete mixing model was known to occur, but for the intersecting angle greater than 150°, no complete solute mixing occurred. Overall, the transition from the complete mixing model to the streamline-routing model occurred for Pe = 0.1 - 100, but it highly depended on the intersecting angles. Specifically, the transition occurred at Pe = 0.1 - 10 for intersecting angles ≧ 150° and at Pe = 10 - 100 for intersecting angles ≦ 30°. For Pe > 100, the streamline-routing model was dominant regardless of intersecting angles. For Pe > 1,000, the complete streamline-routing model appeared only for the intersecting angles greater than 150°. For the intersecting angles less than 150°, the streamline-routing model dominated over the complete solute mixing, but solute mixing still occurred at the fracture junction.