• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.89-99
• /
• 1998

한국형 차세대원자로에서는 비상노심 안전주입수가 저온관을 통하지 않고 원자로용기에 직접 주입된다. 원자로용기의 가압열충격과 열수력적 관점에서 최적의 노즐위치를 결정하기 위해서 전산유체역학을 활용하였다. 상용 전산유체코드인 CFX를 이용하여 원자로 하향유로를 모사하는 해석대상 격자를 다중불록으로 형성한 다음 유동장을 비압축성 Navier-Stokes 운동량 방정식, 에너지 방정식과 표준 k-ε 난류모형 등으로 모형화하여 3차원 비정상상태 계산을 수행하였다. CFX에서는 경계 밀착좌표계, 비엇물림격자와 SIMPLE 알고리즘을 사용한다. 본 연구결과 원자로용기의 가압열충격 관점에서 가장 보수적인 사고인 증기관 파단사고시에도 열적혼합이 잘 일어나 가압열충격이 발생할 가능성이 없는 것으로 판단되며 안전주입수 노즐이 저온관 바로 위에 위치할 때 원자로 하향유로 내의 온도 분포가 가장 균일하여 열적 혼합 관점에서는 최적의 위치로 판단된다.

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

하천의 지형을 조사하고 계측하는 것은 하천을 연구하는 전문가들에게 필수적인 일이다. 하지만 하천의 지형을 계측하는 것은 쉽지 않으며, 조사를 하여도 유사의 이송으로 인하여 하천의 지형은 시간이 지남에 따라 변하게 된다. 그러므로 실험이나 모델링을 통하여 하천의 지형을 예측하고 모의하는 것은 중요한 연구이다. 모델링을 이용하여 유사이송에 의한 하상변동을 잘 예측하기 위해서는 하천의 복잡한 흐름을 정확히 모의하는 것이 중요하며 유사를 발생시키는 힘인 하상전단응력을 정확히 산정하는 것 또한 중요하다. 하상의 전단응력을 산정하는 방법으로는 대표적으로 로그법칙에 의한 방법, 레이놀즈응력 분포를 이용한 방법, 난류운동에너지를 이용한 방법 등이 있다. 앞서 말한 방법으로 산정된 전단응력 값은 차이를 보이며, 이는 하상변동을 정확히 모의하는 것에 문제를 발생시킬 수 있다. 따라서 본 연구에서는 곡선좌표계를 이용하여 3차원 유동 및 하상변동을 모의할 수 있는 수치모형을 이용하여 전단응력 산정 방법에 따른 하상변동량을 분석하는 것이다. 하천의 복잡한 흐름을 정확히 모의하기 위하여 본 연구에서는 RANS (Reynolds Averaged Navier-Stokes) 방정식을 3차원으로 해석하여 흐름 계산을 하였고 유사량 산정공식과 Exner 방정식을 이용하여 유사이송에 의한 하상변동을 계산하였다. 흐름 계산의 검증을 위하여 선행 연구의 실험을 대상으로 모의하였다. 그리고 곡선으로 된 실험 수로를 대상으로 전단응력 산

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.2
• /
• pp.13-29
• /
• 1996

This paper presents calculated results of the free-surface flow around a Series 60($C_B=0.6$) model. Three-dimensional Navier-Sotkes equations are solved and Baldwin-Lomax algebraic turbulence model is adopted to simulate the high Reynolds-number flow. To reduce computational efforts, velocity components near the wall are extrapolated with a the solved by using the Implicit Approximate Factorization method[2]. The successive-over-relaxation method is used for solving pressure-Poisson equation when obtaining the pressure field projecting the divergence-free velocity field. To simulate the free-surface flows more precisely, the numerical scheme solving the equation for the kinematic boundary condition is very important. In this paper, there numerical schemes are employed and the results are compared with the available experimental data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3B
• /
• pp.261-270
• /
• 2008

Even though there is a great deal of progress in a numerical method of high caliber like SPH, it is very rarely deployed in a water resources community. Despite the great stride in computing environment, depth averaged approach like a nonlinear shallow equation is still efficient tool for flood routing in large watershed, but it can give some misleading information like the inundation height of flood. In this rationale, we numerically simulate the flow into the dry channel, dry channel with an obstacle triggered by the collapse of a two dimensional water column using SPH (Smoothed Particle Hydrodynamics) in order to boost the application of numerical method of high caliber like SPH in a water resources community. As a most severe test of the robustness of SPH, we also carry out the simulation of the flow through a clearance into the wet channel driven by the rapid removal of a water gate. As a hydrodynamic model, we used the Navier-Stokes equation, a numerical integration of which was carried out using SPH. To verify the validity of newly proposed numerical model, we compare the numerically simulated flow with the others in the literature mainly from VOF and MAC, and hydraulic experiments by Martin and Moyce (1952), Koshizuka et al. (1995) and Janosi et al. (2004). It was shown that agreements between the numerical results in this study and hydraulic experiments are remarkable.

• Journal of Wetlands Research
• /
• v.14 no.1
• /
• pp.75-87
• /
• 2012

Hyporheic zone is a region beneath and alongside a stream, river, or lake bed, where there is mixing of shallow groundwater and surfacewater. Hyporheic exchange controls a variety of physical, biogeochemical and thermal processes, and provides unique ecotones in a aquatic ecosystem. Field and experimental observations, and modeling studies indicate that hyporheic exchange is mainly in response to pressure gradients driven by the geomorphological features of stream beds. In the reach scale of a stream, pool-riffle structures dominate the exchange patterns. Flow over a pool-riffle sequence develops recirculation zones and stagnation points, and this flow structures make irregular pressure gradient which is driving force of the hyporheic exchange. In this study, 3 D hydro-dynamic model solves the Reynolds-averaged Navier-Stokes equations for the surface water and Darcy's Law and the continuity equation for ground water. The two sets of equations are coupled via the pressure distribution along the interface. Simulation results show that recirculation zones and stagnation points in the pool-riffle structures dominantly control the upwelling and downwelling patterns. With decrease of recirculation zones, length of donwelling zone formed in front of riffles is reduced and position of maximum downwelling point moves downward. The numerical simulation could successfully predict the behavior of hyporheic exchange and contribute the field study, river management and restoration.

• Journal of Korea Water Resources Association
• /
• v.41 no.2
• /
• pp.197-212
• /
• 2008

Recently, stream flow analysis has been accomplished by one or two dimensional equations and was applied by simple momentum equations and fixed energy conservations which contain many condition limits. In this study, FLOW-3D using CFD (Computational Fluid Dynamics) was applied to stream flow analysis which can solve three dimensional RANS (Reynolds Averaged Navier-Stokes Equation) control equation to find out physical behaviors and the effect of hydraulic structures. Numerical simulation accomplished those results was compared by using turbulence models such as ${\kappa}-{\varepsilon}$, RNG (Renormalized Group) ${\kappa}-{\varepsilon}$ and LES (Large Eddy Simulation). Numerical analysis results have been illustrated by the turbulence energy effects, velocity of flow, water level pressure and eddy flows around the piers and transverse overflow type structures. These results will be able to used by basis data that catch hold of effects on long-term bed elevation changes, sediment accumulations, scours and water aggravations by removal of obsolete transverse over flow type structures in urban stream.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.14 no.4
• /
• pp.273-279
• /
• 2011

Turbulent flow analysis around a wind turbine blade was performed to evaluate the power performance of offshore wind turbine. Fluent package was utilized to solve the Reynolds-averaged Navier-Stokes equations in non-inertial rotating coordinates. The realizable k-$\varepsilon$ model was used for turbulence closure and the grid system combining structured and unstructured grids was generated. In the first, lift and drag forces of 2-D foil section were calculated and compared with existing experimental data for the validation. Then torque and thrust of the wind turbine blade having NACA 4-series sections were calculated with fixed pitch angle and rpm. Tip speed ratio was varied by changing wind speed. In the next, three kinds of end plate were attached at the tip of blade in order to increase the power of the wind turbine. Among them the end plate attached at the suction side of the blade was found to be most effective. Furthermore, performance analysis with tilt angle and rake was also performed.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.3
• /
• pp.72-82
• /
• 1995

In the present paper, flow characteristics and free surface waves generated by a submerged hydrofoil advancing with an uniform speed are calculated. Using a numerical method based on a MAC(Marker And Cell) method, the Navier-Stokes and the continuity equations are solved to simulate flow fields around the hydrofoil. Computations are carried out in a rectangular grid system in which grids are concentrated near the foil and the free surface to improve numerical accuracies. Viscous flow phenomenas including pressure distributions are computed. Moreover, the influences of submerged depths upon the generated wave profiles and the wave breaking phenomena are also investigated. Experiments are performed at the towing tank of Inha University to measure free surface wave elevations due to the advancing hydrofoil. The computational results are compared with the present and the other available experimental data to show the accuracy of the numerical method developed.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.484-491
• /
• 2007

Numerical experiments using a numerical wave tank have been performed to verier the nonlinear interaction between monochromatic waves and a submerged horizontal plate. As a model for numerical wave tank, we used a higher-order Boundary Element Method(BEM) based on fully nonlinear potential flow theory and CADMAS-SURF for solving Navier Stokes equations and exact free surface conditions. Both nonlinear models are able to predict the higher harmonic generation in the shallow water region over a submerged horizontal plate. CADMAS-SURF, which involves the viscous effect, can evaluate the higher harmonic generation by flow separation and vortices at the each ends of plate. The comparison of reflection and transmission coefficients with experimental results(Patarapanich and Cheong, 1989) at different lengths and submergence depths of a horizontal plate are presented with a good agreement. It is found that the transfer of energy from the incident fundamental waves to higher harmonics becomes larger as the submergence depth ratio decreases and the length ratio increases.