• 대한기계학회논문집B
• /
• 제21권9호
• /
• pp.1149-1164
• /
• 1997

In this study, the applicability of the RNG k-.epsilon. model to the analysis of the complex flows is studied. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretized by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-.epsilon. model of three complex flows, i.e., the flow over a backward-facing step and a blunt flat plate, the flow around a 2D model car are compared to these from the standard k-.epsilon. model and experimental data. That of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly and the spray characteristics within a chamber of direct injection model engine are compared to these from the standard k-.epsilon. model and experimental data. The results of reattachment length, separated eddy size, average surface pressure distribution using the RNG k-.epsilon. model show more reasonable trends comparing with the experimental data than those using the modified k-.epsilon. model. Although the predicted rms velocity using the modified k-.epsilon. model is lower considerably than the experimental data in incylinder flow with poppet valve, predicted axial and radial velocity distributions at the valve exit and in-cylinder region show good agreements with the experimental data. The spray tip penetration predicted using the RNG k-.epsilon. model is more close to the experimental data than that using the modified k-.epsilon. model. The application of the RNG k-.epsilon. model seems to have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly and the spray characteristics over the modified k-.epsilon. model.

• Wind and Structures
• /
• 제17권6호
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.

• Nuclear Engineering and Technology
• /
• 제51권6호
• /
• pp.1487-1503
• /
• 2019

The methods and performance of a pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT employing the drift-flux model are presented. This code aims at providing an accurate yet fast core thermal-hydraulics solution capability to high-fidelity multiphysics core analysis systems targeting massively parallel computing platforms. The four equation drift-flux model is adopted for two-phase calculations, and numerical solutions are obtained by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models involving turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analyses. ESCOT is parallelized by a domain decomposition scheme that involves both radial and axial decomposition to enable highly parallelized execution. The ESCOT solutions are validated through the applications to various experiments which include CNEN $4{\times}4$, Weiss et al. two assemblies, PNNL $2{\times}6$, RPI $2{\times}2$ air-water, and PSBT covering single/two-phase and unheated/heated conditions. The parameters of interest for validation include various flow characteristics such as turbulent mixing, spacer grid pressure drop, cross-flow, reverse flow, buoyancy effect, void drift, and bubble generation. For all the validation tests, ESCOT shows good agreements with measured data in the extent comparable to those of other subchannel-scale codes: COBRA-TF, MATRA and/or CUPID. The execution performance is examined with a mini-sized whole core consisting of 89 fuel assemblies and for an OPR1000 core. It turns out that it is about 1.5 times faster than a subchannel code based on the two-fluid three field model and the axial domain decomposition scheme works as well as the radial one yielding a steady-state solution for the OPR1000 core within 30 s with 104 processors.

• 한국태양에너지학회 논문집
• /
• 제31권3호
• /
• pp.23-28
• /
• 2011

In this study, the ejector design was modeled using Fluent 6.3 of FVM(Finite Volume Method) CFD(Computational Fluid Dynamics) techniques to resolve the flow dynamics in the ejector. A vacuum system with the ejector has been widely used because of its simple construction and easy maintenance. Ejector is the main part of the desalination system, of which designs determine the efficiency of system. The effects of the ejector was investigated geometry and the operating conditions in the hydraulic characteristics. The ejector consists mainly of a nozzle, suction chamber, mixing tube (throat), diffuser and draft tube. Liquid is supplied to the ejector nozzle, the fast liquid jet produced by the nozzle entrains and the non condensable gas was sucked into the mixing tube. The multiphase CFD modeling was carried out to determine the hydrodynamic characteristics of seawater-air ejector. Condition of the simulation was varied in entrance mass flow rate (1kg/s, 1.5kg/s, 2kg/s, 2.5kg/s, 3kg/s), and position of driving nozzle was located from the central axis of the suction at -10mm, 0mm, 10mm, 20mm, 30mm.. Asaresult, suction flow velocity has the highest value in central axis of the suction.

• 한국태양에너지학회:학술대회논문집
• /
• 한국태양에너지학회 2011년도 추계학술발표대회 논문집
• /
• pp.68-73
• /
• 2011

This study was performed numerical analysis in order to analyze liquid film flow of heat exchanger tube arrangement and configuration of evaporative multi effect distillation system using medium-temperature. Simulation was accomplished the two-dimensional calculations using commercial analyses program FLUENT based on the FVM(finite volume method). Fresh water generator of this study used Shell & Tubes heat exchanger with Cu_Ni tube, configuration of tube used bare tube and corrugated tube, and arrangement of tube used in-line array and staggered array. Performance of heat exchanger through the formation of liquid film was compared and analyzed. Liquid film flow occurred that falling on heat exchanger tube wall. Result of simulation showed that liquid film thickness of in-line arrangement was found 0.57mm with bare tube and 0.67mm with corrugated tube, respectively. And liquid film thickness of staggered arrangement was found 0.39mm with bare tubes and 0.62mm with corrugated tubes, respectively. Liquid film thickness of corrugated tube showed thicker than bare tube, but heat transfer rates of corrugated tube showed higher than bare tube. The reason was considered that surface area of corrugated tube was wider than bare tube. And liquid film thickness of staggered arrangement showed thinner than in-line arrangement, so thermal performance of staggered arrangement showed higher than in-line arrangement.

• 한국수자원학회논문집
• /
• 제42권12호
• /
• pp.1079-1089
• /
• 2009

본 연구에서는 마른하도 및 복잡한 지형에서의 파의 전파와 같은 수공학 분야에서 해결하기 어려운 문제를 해석하기 위한 고정확도 2차원 수치모형을 개발하기 위해, unsplit 유한체적기법과 HLLC Riemann 해법을 이용한 흐름율 계산으로 쌍곡선형 적분 보존형의 2차원 천수방정식을 해석하였다. Unsplit 기법의 적용을 위해 하상경사항은 발산정리를 이용하여 이산화한 형태를 적용하였으며, 흐름율과 생성항의 균형을 이루기 위해 수면경사법을 시간과 공간에 대해 2차정확도를 가지는 MUSCL 기법과 연계하였다. 그리고 적용한 생성항 처리기법과 흐름율과의 보존특성이 만족함을 보였다. 2차정확도의 사용으로 불연속 지점에서 발생할 수 있는 수치진동을 제거하기 위해서 경사제한자를 사용한 TVD 기법을 적용하였다. 개발모형을 정확해가 존재하는 생성항이 없는 1차원 댐 붕괴 흐름에 적용하여 흐름율 계산의 정확성을 검증하였고, 하상융기를 가진 하도의 정상류 및 천이류 모의를 통해 개발모형의 보존특성을 검증하였으며, 하상경사 및 단면의 확대/축소구간이 존재하는 2차원 댐 붕괴 흐름에 적용하여 개발모형의 적용성을 검증하였다.

• 한국지반공학회논문집
• /
• 제28권12호
• /
• pp.5-15
• /
• 2012

본 연구에서는 일련의 현장 열응답 시험결과를 동일한 지중열교환기와 지반 조건에 대한 CFD(Computational Fluid Dynamics) 수치해석 결과와 비교하고 역해석을 통해 지반의 열전도도를 평가하였다. 총 6개의 보어홀을 원주에 소재하고 있는 시험시공 현장에 설치하였으며 순환 파이프의 형상과 그라우트 재료에 대한 수직 밀폐형 지중열교환기의 성능을 비교하기 위해 일반적인 U형 순환 파이프와 새롭게 개발된 3공형 순환 파이프를 보어홀 내 시공하였다. 수치해석은 CFD 해석 프로그램인 FLUENT를 적용하여 3차원 열전달 거동 해석을 수행하였으며 각각의 보어홀에 대해 시간에 따른 순환수의 유입, 유출 온도 차이와 지반의 깊이별 온도변화를 User Define Function (UDF)을 이용하여 실제 조건을 모사하였다. 주어진 보어홀 조건과 실내시험을 통해 시험시공 현장의 열 물성을 입력치로 적용하여 수치 해석을 수행하였으며, 현장 열응답 시험에서 측정된 시간에 따른 유입, 유출 순환수의 온도 변화를 모사하였다. 수치해석 결과, 지반의 열전도도를 3W/mK로 적용하였을 때 보다 4W/mK일 때 현장 열응답 시험과 유사한 결과를 얻었다.