• Nuclear Engineering and Technology
• /
• v.16 no.2
• /
• pp.47-57
• /
• 1984

Mathematically-rigorous time-volume averaged conservation equations were simplified to established the differential equations of THERMIT-6S, which is a two-fluid 3-D code. The difference equations of THERMIT-6S were obtained by discretizing the proceeding set of differential equations. The spatial discretization is characterized by a first-order spatial scheme, donor cell method, and staggered mesh layout. For time discretization, a first order semi-implicit scheme treats implictly sonic terms and terms relating to local transport phenomena and explicitly convective terms. The results were linearized by the Newton-Raphson method. In order to construct the reduced pressure equation, the linearized equations were manipulated so that all variables are coupled between mesh cells through only the pressure variable. By simulating numerically the OPERA-15 experiment, it was found that THERMIT-6S is a very powerful code in predicting reactor behavior after sodium boiling including flow coastdown, reversal flow and flow oscillation.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.22 no.6
• /
• pp.750-757
• /
• 2016

This study presents the wave run-up height around single and multiple surface-piercing cylinders according to wave period and steepness. In order to simulate 3D incompressible viscous two-phase turbulent flow, the present study employed a volume of fluid (VOF) method with realizable $k-{\varepsilon}$ turbulence model based on commercial Computational Fluid Dynamics (CFD) software, "STAR-CCM". The wave periods at model scale were 1.269s and 1.692s for a single cylinder and 1.716s for multiple cylinders. In each case, wave steepness of has 1/30 and 1/16 were used, respectively. Consequently, the results for wave run-up height with regard to wave steepness and period were compared with those of relevant previous experimental studies. The numerical simulation results showed a good qualitative agreement with experiments.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.35 no.6
• /
• pp.1237-1246
• /
• 2015

A novel technique was proposed to calculate fluxes accurately by separation of flow area into a part of step face which is dominated by flow resistance of it and an upper part which is relatively less affected by the step face in analyzing shallow-water flows over discontinuous topography. This technique gives fairly good agreement with exact solutions, 3D simulations, and experimental results. It has been possible to directly analyze shallow-water flows over discontinuous topography by the technique developed in this study. It is expected to apply the developed technique to accurate evaluation of overflows over weirs or retaining walls (riverside roads) and areas flooded by the inundation in the city covered in discontinuous topography.

• Journal of Korea Water Resources Association
• /
• v.48 no.12
• /
• pp.1023-1035
• /
• 2015

This study applies 3-D N-S solver based on PBM (Porous Body Model), LED-WASS-3D ver 2.0 to directly analyze non-linear interaction of seawater-freshwater-coastal aquifer in order to simulate the seawater infiltration into coastal aquifer. This numerical simulation is the first trial in Korea, as well as unusual and new numerical analysis abroad. Firstly, to validate the applied numerical model, the validity and effectiveness was verified for the numerical model by comparing and considering it with the result of laboratory experiment for seawater-freshwater interface in coastal aquifer. And then it simulated the seawater infiltration into coastal aquifer considering the changed levels of seawater and groundwater in order to analyze the distribution characteristics of flow field and seawater-freshwater interface of coastal aquifer as the level difference between seawater and groundwater and rate of seawater level (${\Delta}h/h$) increased. In addition, the characteristics of seawater infiltration were analyzed from the vertical salinity in the coastal aquifer by ${\Delta}h/h$, which cannot be obtained from existing non-diffusion numerical models. Finally, it analyzed the effect of ${\Delta}h/h$ on the seawater infiltration distance in coastal aquifer, which was indexed.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.20 no.2
• /
• pp.63-75
• /
• 2017

In offshore environment, an accurate estimation of a wave-structure interaction has been an important issue for safe and cost effective design of fixed and floating offshore structures exposed to a harsh environment. In this study, a wave-structure interaction around a circular column was investigated with regular waves. To simulate 3D two-phase flow, open source computational fluid dynamics libraries, called OpenFOAM, were used. Wave generation and absorption in the wave tank were activated by the relaxation method, which implemented in a source term. To validate the numerical methods, generated Stokes 2nd-order wave profiles were compared with the analytic solution with deep water condition. From the validation test, grid longitudinal and vertical sizes for wave length and amplitude were selected. The simulated wave run-up and wave loads on the circular column were studied and compared with existing experimental data.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.320-334
• /
• 2020

This study presents an application of hydro-mechanical coupled Particle Flow Code 3D (PFC3D) to simulation of fluid injection induced fault slip experiment conducted in Mont Terri Switzerland as a part of a task in an international research project DECOVALEX-2019. We also aimed as identifying the current limitations of the modelling method and issues for further development. A fluid flow algorithm was developed and implemented in a 3D pore-pipe network model in a 3D bonded particle assembly using PFC3D v5, and was applied to Mont Terri Step 2 minor fault activation experiment. The simulated results showed that the injected fluid migrates through the permeable fault zone and induces fault deformation, demonstrating a full hydro-mechanical coupled behavior. The simulated results were, however, partially matching with the field measurement. The simulated pressure build-up at the monitoring location showed linear and progressive increase, whereas the field measurement showed an abrupt increase associated with the fault slip We conclude that such difference between the modelling and the field test is due to the structure of the fault in the model which was represented as a combination of damage zone and core fractures. The modelled fault is likely larger in size than the real fault in Mont Terri site. Therefore, the modelled fault allows several path ways of fluid flow from the injection location to the pressure monitoring location, leading to smooth pressure build-up at the monitoring location while the injection pressure increases, and an early start of pressure decay even before the injection pressure reaches the maximum. We also conclude that the clay filling in the real fault could have acted as a fluid barrier which may have resulted in formation of fluid over-pressurization locally in the fault. Unlike the pressure result, the simulated fault deformations were matching with the field measurements. A better way of modelling a heterogeneous clay-filled fault structure with a narrow zone should be studied further to improve the applicability of the modelling method to fluid injection induced fault activation.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2006.06a
• /
• pp.45-46
• /
• 2006

Experimental and numerical analysis on the heat transfer behaviors and the associated frictional loss in a rectangular channel with two inclined perforated baffles($\;5^{\circ}$) mounted on the bottom plate has been systematically performed. The parametric effects of perforated baffles (3, 6, 9 holes) and flow Reynolds number on heat transfer characteristics of the heated target surface are explored. A combination of two baffles of same overall size was considered and the flow Reynolds number for this study is varied between 28,900 and 61,800. Comparisons of the experimental data with the numerical results by commercial code CFX 5.7 are made. As for the investigation of heat transfer behaviors on local Nusselt number with the two baffles installed at $x/D_h=0.8\;and\;x/D_h=8.0$, it is evident that there exist an optimum perforation density to maximize heat transfer coefficients; i.e., the maximum Nusselt number decreases with increasing number of holes.

• Journal of Fisheries and Marine Sciences Education
• /
• v.28 no.3
• /
• pp.831-840
• /
• 2016

Hydrodynamic condition can be used to predict particle movement within water column and the results used to optimize environmental conditions for effective site selection, setting of environmental quality standard, waste dispersion, and pathogen transfer. To predict the extent of movement of particle from land, 3D hydrodynamic model that includes particle tracking module was applied to Geoje Bay and to calibrate particle tracking model, floating buoy measurement is operated. The model results show that short time is required for particles released into system from river to be transported to the shellfish farming area. It takes about 2 days for the particles to shellfish farming area under mean flow condition. It meant Geoje Bay, especially shellfish farming area is vulnerable to anthropogenic waste from river.

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

최근 기상변화로 인한 게릴라성 집중호우 등의 증가로 재해 발생이 빈번하게 발생하고 있으며, 특히 도시지역에서의 홍수피해로 막대한 인명과 재산피해가 발생하고 있는 실정이다. 도시지역 내로 유입되는 홍수는 매우 복잡한 흐름특성을 나타내는 것으로 알려져 있다. 홍수의 대부분은 도로와 합류부를 통해 흐르므로 도시지역 내에서 홍수로 인한 피해를 저감하기 위해서는 도로 교차로에서의 홍수파의 전파특성을 분석하는 것이 중요하다. 교차로에서 흐름은 유입되는 유량이 상대적으로 적고 도로의 경사가 비교적 완만한 경우 상류 특성을 나타내나 이와 반대일 경우에는 상류 및 사류가 공존하는 복잡한 수면형을 나타낸다. (Mignot, 2005) 현재까지 국내외적으로 홍수흐름 분석에 관한 대부분의 수치적 연구에서는 주로 수직방향으로의 홍수파 변화양상에 대해서는 무시하며, 물의 압력이 정수압이라는 2차원적 가정으로 수행되어져 왔다. 그리고 대부분의 연구는 도시지역을 구성하는 건물 등이 홍수흐름에 어떤 영향을 미치는가 및 직각 교차로에 대한 분석을 수행하였으나 사각 교차로에 대한 홍수흐름에 대한 연구는 아직 미비한 실정이다. 본 연구는 3차원 수직모형을 이용하여 사각 수로 합류부에서의 상류 및 사류에 대한 흐름특성을 수리모형실험 결과와 비교분석하였다.

• Journal of Korea Water Resources Association
• /
• v.41 no.12
• /
• pp.1187-1198
• /
• 2008

Many reservoirs in Korea and their downstream environments are under increased pressure for water utilization and ecosystem management from longer discharge of turbid flood runoff compared to a natural river system. Turbidity($C_T$) is an indirect measurement of water 'cloudiness' and has been widely used as an important indicator of water quality and environmental "health". However, $C_T$ modeling studies have been rare due to lack of experimental data that are necessary for model validation. The objective of this study is to validate a coupled three-dimensional(3D) hydrodynamic and particle dynamics model (ELCOM-CAEDYM) for the simulation of turbid density flows in stratified Daecheong Reservoir using extensive field data. Three different groups of suspended solids (SS) classified by the particle size were used as model state variables, and their site-specific SS-$C_T$ relationships were used for the conversion between field measurements ($C_T$) and state variables (SS). The simulation results were validated by comparing vertical profiles of temperature and turbidity measured at monitoring stations of Haenam(R3) and Dam(R4) in 2004. The model showed good performance in reproducing the reservoir thermal structure and propagation of stream density flow, and the magnitude and distribution of turbidity in the reservoir were consistent with the field data. The 3D model and turbidity modeling framework suggested in this study can be used as a supportive tool for the best management of turbidity flow in other reservoirs that have similar turbidity problems.