• 대한토목학회논문집
• /
• 제32권6C호
• /
• pp.267-274
• /
• 2012

오염물질의 거동은 간극 배열과 연결망에 의해 결정되지만 흙에서 오염물질의 이동을 계산하는 전통적인 접근은 거시 스케일에서 적용된다. 나비에-스토크스 방정식을 풀기 위해 소요되는 컴퓨팅 비용 때문에 간극 스케일에서 이동과 결과를 관찰하기 쉽지 않다. 본 논문에서는 단일상 다성분 유체유동에서 국부적인 유속과 밀도를 효과적으로 평가할 수 있는 격자 볼츠만 방법에 대해 설명한다. 오염물질의 시공간적 거동은 유체 유동의 이송에 의해 명시적으로 결정된다. 두 가지 형태의 이상화된 간극은 유체의 경로를 제공한다. 또한 오염물질 이동, 유속장, 오염물질의 평균 농도는 정상상태의 유동에서 계산된다. 굴곡비와 같은 간극 형상은 오염물질 거동의 영향을 준다. 이는 흙이나 암반의 불연속면에서 동일한 간극률를 가져도 간극의 배열과 형상의 중요함을 강조한다.

• 한국지하수토양환경학회:학술대회논문집
• /
• 한국지하수토양환경학회 2002년도 총회 및 춘계학술발표회
• /
• pp.122-126
• /
• 2002

Contaminant transport in porous media is characterized by solving an advection-dispersion equation(ADE). The ADE can cover equilibrium phenomena of interest, which include sorption, decay, and chemical reactions. Among these phenomena, sorption mechanism is described by several types of sorption isotherm. If we assume the sorption isotherm as linear, the solution of ADE can be easily procured. However, if we consider the sorption isotherm as non-linear isotherm like a Dual Reactive Domain Model (DRDM), the resulting differential equation becomes non-linear. In this case, the solution of ADE cannot be easily acquired by an analytic method. In this paper, we present the numerical analysis of ADE using a DRDM. The results reveal that even if sorption data may be fitted well using linear or non-linear isotherm, the characteristics of contaminant transport of the two cases are different from each other. To be concrete, the retardation of linear isotherm has stronger effect than that of the DRDM. As the non-linearity of sorption isotherm increases, the difference of retardation effects of the two cases becomes larger. For a pulse source, the maximum concentration of the linear model is higher than that of the DRDM, but the plume of the DRDM moves faster than that of the linear model. Behaviors of contaminant transport using the DRDM are consistent with common features of a linear model. For instance, biodegradation effect becomes larger as time goes by The faster the seepage velocity is, the faster the plume of contaminant moves. The plume of the contaminant is distributed evenly over overall domain in the event of high dispersion coefficient.

• 대한토목학회논문집
• /
• 제34권3호
• /
• pp.859-871
• /
• 2014

본 연구에서는 염분과 온도차에 의한 밀도류를 해석하기 위하여 기존의 3차원 파동장 모델(LES-WASS-3D ver. 1.0)을 토대로 염분과 온도에 관한 상태방정식 및 3차원 이류-확산 방정식을 도입하여 새로운 수치모델(LES-WASS-3D ver. 2.0)을 개발하였다. 새롭게 개발한 수치모델의 적용성을 검토하기 위하여 기존의 수리모형실험결과와 비교 분석하였다. 그 결과, 이류-확산하는 해수의 형태 및 연직 염분농도 뿐만 아니라 방출하는 온수의 유속 및 온도를 매우 잘 재현하는 것으로 나타났다. 이로써 본 연구에서 개발한 수치모델의 타당성 및 유효성이 검증되었다.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2008년도 정기총회 및 학술발표대회
• /
• pp.799-802
• /
• 2008

General behaviors based on hydraulic characteristics of natural streams and channels have been recently analyzed and developed via various numerical models. However in the states of natural hydraulics, an experimental research must be performed simultaneously with the mathematical analysis due to effects of hydraulic properties such as meander, sediment, and so on. In this study based on 2-D advection-dispersion equation, flow and tracer experiments were performed in the S-curved meandering laboratory channel with a rectangular cross-section. The channel was equipped with instrument carriages which was equipped with an auto-traversing system to be used with velocity measuring sensors throughout the depth and breadth of the flow field. To measure concentration distribution of the salt solution was adjusted to that of the flume water by adding methanol and a red dye (KMnO4) was added to aid the visualization of the tracer cloud, the tracer was instantaneously injected into the flow as a full-depth vertical line source by the instantaneous injector and the initial concentration of the tracer was 100,000 mg/l. The secondary current as well as the primary flow pattern was analyzed to investigate the flow distribution in the meandering channels. The velocity distribution of the primary flow for all cases skewed toward the inner bank at the first bend, and was almost symmetric at the crossovers, and then shifted toward the inner bank again at the next alternating bend. Thus, one can clearly notice that the maximum velocity occurs taking the shortest course along the channel, irrespective of the flow conditions. The result of the tracer tests shows that pollutant clouds are spreading following the maximum velocity lines in each cases with various mixing patterns like superposition, separation, and stagnation of pollutant clouds. Flow characteristics in each cases performed in this study can be compared with tracer dispersion characteristics with using evaluation of longitudinal and transverse dispersion coefficients(LDC, TDC). As expected, LDC and TDC in meandering parts have been evaluated with increasing distribution and straight parts have effected to evaluate minimum of LDC and TDC due to symmetric flow patterns and attenuations of secondary flow.

• 대한기계학회논문집B
• /
• 제30권8호
• /
• pp.796-805
• /
• 2006

In this study, the problem of using the mixing index as a measure of the mixing performance for a certain flow field has been discussed. The flow model subjected to this study is the two-dimensional unsteady lid-driven cavity flow. The transport equation for the concentration within the cavity was solved by using the finite volume method where the convective terms are discretized with the central difference scheme. It was shown that both the concentration dispersion and the mixing index depend highly on the initial distribution of the concentration, and therefore the mixing index obtained from the concentration dispersion equation loses its universal applicability.

• 한국농공학회지
• /
• 제42권2호
• /
• pp.71-77
• /
• 2000

A multi-region model for solute transport through saturated soils has been developed to describe preferential flow. The model consists of numerous discrete pore groups, which are characterized by a discrete dispersion coefficient, flow velocity, and porosity . The hydraulic properties for each pore group are derived from a soil's hydraluic conductivity and soil water characteristic functions . Flow in pore group is described by the classical advection-disersion equation (ADE). An implict finite difference scheme was applied to the governing equation that results in a block-tridiagonal system of equations that is very efficient and allows the soil to be divided into any number of pore groups. The numerical technique is derived from methods used to solve coupled equations in fluid dynamics problems and can also be applied to the transport of interacting solutes. The results of the model are compared to the experimental data from published papers. This paper contributes on the characteristics of the method when applied to the parallel porosity model to describe preferential flow of solutes in soil.

• Korean Journal of Hydrosciences
• /
• 제5권
• /
• pp.85-97
• /
• 1994

A hybrid finite difference method for the longitudinal dispersion equation, which is based on combining the Holly-Preissmann scheme with fifth-degree Hermite interpolating polynomial and the generalized Crank-Nicholson scheme, is described and comparatively evaluated with other characteristics-based numerical methods. Longitudinal dispersion of an instantaneously-loaded pollutant source is simulated, and computational results are compared with the exact solution. The present method is free from wiggles regardless of the Courant number, and exactly reproduces the location of the peak concentration. Overall accuracy of the computation increases for smaller value of the weighting factor, $\theta$of the model. Larger values of $\theta$ overestimates the peak concentration. Smaller Courant number yields better accuracy, in general, but the sensitivity is very low, especially when the value of $\theta$ is small. From comparisons with the hybrid method using cubic interpolating polynomial and with splitoperator methods, the present method shows the best performance in reproducing the exact solution as the advection becomes more dominant.

• 대한기계학회논문집B
• /
• 제27권3호
• /
• pp.347-356
• /
• 2003

The characteristics of near wakes of circular cylinders with serrated fins are investigated experimentally using a hot-wire anemometer for various freestream velocities. Near wake structures of the fin tubes are observed using a phase average technique. With increasing fin height and decreasing fin pitch. oscillation of streamwise velocity increases. It file oscillation of lateral velocity decreases. The time averaged V-component velocity distribution of the finned tube is contrary to that of the circular cylinder due to the different strength of entrainment flow. This strength is affected by the distance of (equation omitted) = 1.0 contour lines. (equation omitted) = 1.0 contour line approaches to the wake center line when the fin density is increased. When the distance between (equation omitted) = 1.0 contour lines comes close the shear force should be increased and the flow toward the wake center line can be more strengthened because of the shear force. Factors related to the velocity recovery in the near wake of the finned tube are attributed to tile turbulent intensity, the boundary layer thickness. the position and strength of entrainment process.

• 대한기계학회논문집B
• /
• 제32권6호
• /
• pp.429-445
• /
• 2008

Two-phase compressible flow fields of air-water are investigated numerically in the fixed Eulerian grid framework. The phase interface is captured via volume fractions of each phase. A way to model two phase compressible flows as a single phase one is found based on an equivalent equation of states of Tait's type for a multiphase cell. The equivalent single phase field is discretized using the Roe‘s approximate Riemann solver. Two approaches are tried to suppress the pressure oscillation phenomena at the phase interface, a passive advection of volume fraction and a direct pressure relaxation with the compressible form of volume fraction equation. The direct pressure equalizing method suppresses pressure oscillation successfully and generates sharp discontinuities, transmitting and reflecting acoustic waves naturally at the phase interface. In discretizing the compressible form of volume fraction equation, phase interfaces are geometrically reconstructed to minimize the numerical diffusion of volume fraction and relevant variables. The motion of a projectile in a water-filled tube which is fired by the release of highly pressurized air is simulated presuming the flow field as a two dimensional one, and several design factors affecting the projectile movement are investigated.

• 대한토목학회논문집
• /
• 제11권4호
• /
• pp.103-112
• /
• 1991

본(本) 연구(硏究)는 일정균등(一定均等)한 열적(熱的) 특성(特性)을 가지고 있으며 상변화(相變化) 없는 등방(等方) 이질성(異質性)의 3차원(次元) 대수층계(帶水層系)의 열(熱)흐름과 정상상태(定常狀態)의 지하수(地下水)흐름을 모의발생(模擬發生)할 수 있는 유한차분(有限差分) 모형(模型)을 확립(確立)한 것이다. 이 모형(模型)은 대규모(大規模) 지하수(地下水) 흐름체계(體系)에서 폐기물(廢棄物)의 지하저류시(地下貯溜時) 지하수(地下水) 흐름과 발생(發生) 혹은 주입(注入)된 열(熱)의 흐름을 예측(豫測) 분석(分析)하기 위하여 확립(確立)된 것이다. 이러한 대수층계(帶水層系)의 지하수(地下水) 흐름에 작용(作用)하는 조건(條件)으로는 강우주입(降雨注入)으로 인한 수문학적(水文學的) 조건(條件)과 고정(固定)된 수리수두(水理水頭) 경계조건(境界條件) 등(等)이 포함(包含)되고, 열(熱)흐름에는 지열(地熱)의 흐름, 지표면(地表面)으로의 전도(傳導), 주입(注入)에 의한 이류(移流), 고정(固定) 수두경계(水頭境界)로 향(向)한 또는 고정수두경계(固定水頭境界)로 부터의 이류(移流) 등(等)이 포함(包含)된다. 본(本) 모형(模型)에서는 지하수(地下水)흐름과 열(熱)흐름 방정식(方程式)을 번갈아 푸는 교대반복과정(交代反復過程)을 사용(使用)하고, 두 방정식(方程式)의 계산(計算)에는 직접해법(直接解法)을 사용(使用)한다. 이동시간(移動時間)은 모형공간(模型空間)에서 입자추적(粒子追跡)으로 결정(決定)되며, 분할(分轄)된 구역내(區域內)의 지하수(地下水) 유속(流速)은 구역내(區域內)의 유속(流速)을 선형(線形)으로 보간(補間)하여 계산(計算)한다. 본(本) 모형(模型)을 경상북도(慶尙北道) 영일군(迎日郡) 송라면(松羅面) 지경리(地境里) 일대(一帶)의 지하수계(地下水系)에 적용(適用)하여 이 일대(一帶) 지하(地下) 암반층(岩盤層)의 수두분포(水頭分布), 유동로(流動路), 이동시간(移動時間) 및 지하수온분포(地下水溫分布)를 계산(計算)하여 지하수(地下水) 유동체계(流動體系)를 분석(分析)하였다.