• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.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
• /
• v.17 no.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
• /
• v.51 no.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.

• Journal of the Korean Solar Energy Society
• /
• v.31 no.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.11a
• /
• 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.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1079-1089
• /
• 2009

The objective of this study is to develop the accurate, robust and high resolution two-dimensional numerical model that solves the computationally difficult hydraulic problems, including the wave front propagation over dry bed and abrupt change in bathymetry. The developed model in this study solves the conservative form of the two-dimensional shallow water equations using an unsplit finite volume scheme and HLLC approximate Riemann solvers to compute the interface fluxes. Bed-slope term is discretized by the divergence theorem in the framework of FVM for application of unsplit scheme. Accurate and stable SGM, in conjunction with the MUSCL which is second-order-accurate both in space and time, is adopted to balance with fluxes and source terms. The exact C-property is shown to be satisfied for balancing the fluxes and the source terms. Since the spurious oscillations in second-order schemes are inherent, an efficient slope limiting technique is used to supply TVD property. The accuracy, conservation property and application of developed model are verified by comparing numerical solution with analytical solution and experimental data through the simulations of one-dimensional dam break flow without bed slope, steady transcritical flow over a hump and two-dimensional dam break flow with a constriction.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.12
• /
• pp.5-15
• /
• 2012

In this study, a series of CFD (Computational Fluid Dynamics) numerical analyses were performed in order to evaluate the thermal performance of six full-scale closed-loop vertical ground heat exchangers constructed in a test bed located in Wonju. The circulation HDPE pipe, borehole and surrounding ground formation were modeled using FLUENT, a finite-volume method (FVM) program, for analyzing the heat transfer process of the system. Two user-defined functions (UDFs) accounting for the difference in the temperatures of the circulating inflow and outflow fluid and the variation of the surrounding ground temperature with depth were adopted in the FLUENT model. The relevant thermal properties of materials measured in laboratory were used in the numerical analyses to compare the thermal efficiency of various types of the heat exchangers installed in the test bed. The simulation results provide a verification for the in-situ thermal response test (TRT) data. The CFD numerical back-analysis with the ground thermal conductivity of 4 W/mK yielded better agreement with the in-situ thermal response tests than with the ground thermal conductivity of 3 W/mK.