• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.290-297
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.

• Journal of Korea Water Resources Association
• /
• v.43 no.10
• /
• pp.911-920
• /
• 2010

The closed water bodies, such as reservoirs and lakes, could be contaminated by an inflow of pollutants in the upstream as well as a stratification caused by seasonal natural phenomena. The vertical circulation particularly plays an important role in reduction of environmental pollutants. The factors of the vertical circulation are the temperature, wind, thermal diffusivity and sunlight. The wind is probably the most significant factor among them. Thus, it is necessary to describe the validation and application of a three-dimensional numerical model of wind-induced circulation in a thermally stratified flow. In this paper, a three-dimensional numerical model for the thermally stratified flows is presented. The model is conducted in three steps to calculate the velocity components from the momentum equations in x- and y- axis directions, the elevations from the free surface equation and the temperature from the scalar transport equation. Numerical predictions are compared with available analytical solutions for the sloshing free surface movement in a rectangular basin. The numerical results generally show a reasonable agreement with analytical solutions. And the model is applied to the circulation for the wind induced flow in a thermally stratification. Consequently, the developed model is validated by two verifications and phenomena of the internal flow.

• Journal of Korea Water Resources Association
• /
• v.37 no.12
• /
• pp.979-991
• /
• 2004

In this study, RMA2 and RMA4, the 2-D depth-averaged models, were employed to simulate the two-dimensional mixing characteristics of the pollutants in the natural streams. The velocity and depth were first calculated using RMA2, 2-D hydrodynamic model, and then the resulting flow field was inputted to RMA4, 2-D water quality model, to compute the concentration field. RMA models were verified using the velocity and concentration data measured in S-curved meandering channel. The results showed that the RMA2 model simulated well the phenomenon that the maximum velocity line is located at the Inner bank of meandering channel, and the RMA4 model was well adapted to reproduce the general mixing behavior and the separation of tracer clouds. Comparing model simulations with measured data in the field experiments, RMA2 model simulated well general flow field and tendency that the maximum velocity line skewed toward the outer bank which were found in field experiments. The simulations of RMA4 model showed that the center of the tracer cloud tends to follow the path in which the maximum velocity occurs. In this study, the dispersion coefficients are fine-tuned based on the measured coefficients calculated using field concentration data, and the results show reasonable agreement with predictive equations.

• Journal of Korea Water Resources Association
• /
• v.42 no.10
• /
• pp.891-896
• /
• 2009

The puff models have been developed to simulate the advection-diffusion processes of dredging suspended solids, either alone or in combination with Eulerian models. Computational efficiency and accuracy are of prime importance in designing these hybrid approaches to simulate a pollutant discharge, and we characterize two relatively simple Lagrangian techniques in this regard: forward Gaussian puff tracking (FGPT), and backward Gaussian puff tracking (BGPT). FGPT and BGPT offer dramatic savings in computational expense, but their applicability is limited by accuracy concerns in the presence of spatially variable flow or diffusivity fields or complex no-flux or open boundary conditions. For long simulations, particle and/or puff methods can transition to an Eulerian model if appropriate, since the relative computational expense of Lagrangian methods increases with time for continuous sources. Although we focus on simple Lagrangian models that are not suitable to all environmental applications, many of the implementation and computational efficiency concerns outlined herein would also be relevant to using higher order particle and puff methods to extend the near field.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.26 no.1
• /
• pp.48-61
• /
• 1989

In this paper, two-dimensional numerical simulation of GaAs MESFFT with 0.7${\mu}m$ gate length is perfomed. Drift-diffusion model which consider that mobility is a function of local electric field, is used. As a discretization method, instead of FDM (finite difference method) and FEM (finite element method), the Control-Volume Formulation (CVF) is used and as a numerical scheme current hybrid scheme or upwind scheme is replaced by power-law scheme which is very approximate to exponential scheme. In the process of numerical analysis, Peclet number which represents the velocity ratio of drift and diffusion, is introduced. And using this concept a current equation which consider numerical scheme at the interface of control volume, is proposed. The I-V characteristics using the model and numerical method has a good agreement with that of previous paper by others. Therefore, it is confined that it may be useful as a simulator for GaAs MESFET. Besides I-V characteristics, the mechanism of both velocity saturation in drift-diffusion model is described from the view of velocity and electric field distribution at the bottom of the channel. In addition, the relationship between the mechanism and position of dipole and drain current, are described.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.4B
• /
• pp.335-344
• /
• 2006

In this study, sequential mixing model (SMM) was proposed based on the Taylor's theory which can be summarized as the fact that longitudinal advection and transverse diffusion occur independently and then the balance between the longitudinal shear and transverse mixing maintains. The numerical simulation of the model were performed for cases of different mixing time and transverse velocity distribution, and the results were compared with the solutions of 1-D longitudinal dispersion model (1-D LDM) and 2-D advection-dispersion model (2-D ADM). As a result it was confirmed that SMM embodies the Taylor's theory well. By the comparison between SMM and 2-D ADM, the relationship between the mixing time and the transverse diffusion coefficient was evaluated, and thus SMM can integrate 2-D ADM model as well as 1-D LDM model and be an explanatory model which can represents the shear flow dispersion in a visible way. In this study, the predicting equation of the longitudinal dispersion coefficient was developed by fitting the simulation results of SMM to the solution of 1-D LDM. The verification of the proposed equation was performed by the application to the 38 sets of field data. The proposed equation can predict the longitudinal dispersion coefficient within reliable accuracy, especially for the river with small width-to-depth ratio.

• Journal of Korea Water Resources Association
• /
• v.43 no.10
• /
• pp.921-932
• /
• 2010

In this study, the adjustment processes of the disturbed channels by sand or gravel mining were investigated by a two dimensional numerical model in the generalized coordinate system. As a numerical scheme, the CIP (cubic interpolated pseudoparticle method) method was used to calculate the advection term in the flow field and central difference method was used to the diffusion term in it. The pit of the channel was partially filled with sediment at the toe of the pit upstream. As time increased, the headcut erosion upstream in the pit was decreased due to the sediment inflow. The almost inflow sediment upstream was trapped into the pit and the sediment deposit wedge migrated downstream in the pit with the steep submerged angle of repose. The numerical model was reproduced well the evolution processes of the channel. The mining pit migrated with speed as the channel was steep, and the numerical results were in overall agreement with the experimental results.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.1
• /
• pp.57-69
• /
• 1995

The development of a finite difference model for cross-shore sediment transport prediction in the surf tone due to the storm surge event is presented in this paper. Using the inhomogeneous diffusion equation with moving boundaries. the present numerical model is found to be robust and efficient and does not possess a number of restrictions imposed in Kriebel and Dean's(1985) numerical model. Our numerical model is validated through comparison with the analytical solution. the data of a large-scale experiment and the field data of Hurricane Eloise. The Present model if able to predict the averaged volumetric erosion rate of a beach due to the time-varying real storm surge hydrographs and satisfies the conservation of sediment between eroded volume in the onshore region and deposited volume in the offshore region. In addition. the present model is able to reasonably predict the recession of a beach with wide berm and dune. and can describe the change of a breaking point by the offshore deposition. From the sensitivity analysis or the present numerical model with various input parameters, it is concluded that the present numerical model is able to analyze the beach change in a reliable manner including the effects of different sizes of sediments.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.395-395
• /
• 2015

하천 내 식생은 수리학, 지형학 및 생태학적으로 매우 중요하다. 식생은 하천 수생물들의 서식처를 제공할 뿐만 아니라 필터와 같은 역할을 함으로써 부유사에 의한 하천오염물의 퇴적을 유발하여 하천의 수질을 개선시킨다. 더욱이, 하천 내 흐름 및 난류구조를 변경시킴으로써 식생주변의 유사 퇴적량 및 분포에 크게 영향을 미치고, 결국 하천의 지형을 변화시킨다. 개수로의 식생에 대한 영향은 주로 실험 및 수치모델을 이용하여 연구되었고 전단면이 식재된 조건에서 식생의 항력계수, 식생역내의 부유사 및 확산에 관한 연구가 진행되어왔다. 이러한 연구를 통해 식생역 내의 전단력이 감소하여 부유사퇴적이 증가하고 식생역과 비식생역 사이의 운동량 교환에 의해 부유사 퇴적이 증가함을 보였다. 그러나 개수로에서 존재하는 유한한 크기의 식생에 의한 흐름 및 유사분포에 관한 연구는 아직 미흡하다. 이에 본 연구에서는 침수하지 않은 원형 식생 주변에서 발생하는 흐름특성을 수치모의 하였다. 침수하지 않은 원형식생 하류에서 발생하는 흐름을 계산하기 위해 2차원 수치모형을 적용하였다. 식생에 의한 저항을 고려하기 위해 운동량 방정식에 식생항을 추가하였고 $k-{\varepsilon}$ 난류모형을 적용하였다. 수치모의 조건은 Zong and Nepf (2012)의 수리실험을 참고하여 수로의 길이는 12 m, 폭은 1.2 m로 설정하였다. 0.13 m 수심을 갖는 개수로에 0.22 m 지름을 갖는 원형식생을 상류경계로부터 1.0 m 떨어진 곳에 설정하였다. 식생의 밀도($6{\sim}77m^{-1}$)를 변화시키면서 원형식생 하류의 흐름거동을 분석하였다. 식생밀도가 높은 경우에는 원형식생 양 측면에서 유발된 전단층들의 상호작용에 의해 하류에서 와류가 발생하였다. 와류가 발생하는 위치에서 난류강도가 가장 크게 나타났다. 그러나 식생밀도가 일정 값보다 낮아지면 와류가 발생하지 않는 것으로 나타났다.

• Journal of the Korean Geotechnical Society
• /
• v.37 no.12
• /
• pp.25-31
• /
• 2021

Boundary conditions for thermal-hydraulic problems of soils play an essential role in the numerical accuracy. This study presents a boundary condition considering the thermo-hydraulic interaction between the ground and the atmosphere. Ground surface energy balance consists of solar radiation, ground radiation, wind convection, latent heat from water evaporation, and heat conduction to the ground. Equations for each heat flux are presented, and numerical analyses are performed in conjunction with the FEM program for the thermal-hydraulic phenomenon of unsaturated soils. Numerical results using the weather data at the Ulsan Meteorological Observatory are similar to the measured surface temperature. Latent heat caused by water evaporation during the daytime lowers the surface temperature of the bare soil, and a thermal equilibrium is reached at nighttime when the effect of the ground condition is significantly reduced. The temperature change of the surface ground is diminished at the deeper ground due to its thermal diffusion. Numerical analysis where the surface ground temperature is the primary concern requires considering the thermo-hydraulic interaction between the ground and the atmosphere.