• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.15 no.4
• /
• pp.274-286
• /
• 2003

Numerical simulations are performed to investigate the turbulent convective heat transfer of the supercritical carbon dioxide flows in vertical and horizontal square ducts. The gas cooling process at the supercritical state experiences a sudden change in thermodynamic and transport properties. This results in the extraordinary variations of the heat transfer coefficients in the supercritical state, which are much different from those of single or two phase flows. Algebraic second moment closure which can include the effects of large thermophysical property variations of carbon dioxide and of buoyancy is employed to model the Reynolds stresses and turbulent heat fluxes in the governing equations. The previous correlations for the turbulent heat transfer coefficient for the supercritical carbon dioxide flows couldn't reflect the buoyancy effect. The present results are used to establish a new heat transfer coefficient correlation including the effects of large thermophysical property variation and buoyancy on in-duct cooling process of supercritical carbon dioxide.

• Journal of Korean Association for Spatial Structures
• /
• v.15 no.1
• /
• pp.75-84
• /
• 2015

This paper presents the dimensionless wall distance, y+ effect on SST turbulent model for wind turbine blade. The National Renewable Energy Laboratory (NREL) Phase VI wind turbine was used for the study, which the wind tunnel and structural test data has publicly available. The near wall treatment and turbulent characteristics have important role for proper CFD simulation. Most of the CFD development in this area is focused on advanced turbulence model closures including second moment closure models, and so called Low-Reynolds (low-Re) number and two-layer turbulence models. However, in many cases CFD aerodynamic predictions based on these standard models still show a large degree of uncertainty, which can be attributed to the use of the $\epsilon$-equation as the turbulence scale equation and the associated limitations of the near wall treatment. The present paper demonstrates the y+ definition effect on SST (Shear Stress Transport) turbulent model with advanced automatic near wall treatment model and Gamma theta transitional model for transition from lamina to turbulent flow using commercial ANSYS-CFX. In all cases the SST model shows to be superior, as it gives more accurate predictions and is less sensitive to grid variations.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.264-267
• /
• 2010

This research focuses on the development of numerical code to deal with compressible two phase flow around three dimensional objects combined with cavitation model suggested by Weishyy et al. with k-e turbulent model. The cryogenic cavitation is carried out by considering the thermodynamic effect on physical properties of cryogenic fluids in physical point of view and implementing the temperature sensitivity in the energy equation of the government equations in numerical point of view, respectively. The formulation has been extensively validated for both liquid nitrogen and liquid hydrogen by simulating the experiments of Hord on hydrofoils. Then, simulations of cavitating turbopump inducers at their design flow rate are presented. Results over a broad range of Nss numbers extending from single-phase flow conditions through the critical head break down point are discussed. In particular, thermal depression effects arising from cavitation in cryogenic fluids are identified and their impact on the suction performance of the inducer quantified.

• Transactions of the Korean hydrogen and new energy society
• /
• v.26 no.6
• /
• pp.581-589
• /
• 2015

In this study, the viscosity of a liquid carbon dioxide ($LCO_2$) that can potentially be used in a wet feed coal gasifier was evaluated. A pressurized capillary viscometer was employed to obtain the viscosity data of $LCO_2$ using two different methods. During the first method, the measurements were conducted under quasi-steady and high pressure flow conditions where two-phase flow was greatly minimized. The viscosity of $LCO_2$ was determined using turbulent friction relationship. At the second flow condition where unsteady flow is induced, the viscosity of $LCO_2$ was measured using the half-time pressure decay data and was further compared with values calculated by the first method.

• Journal of the Korean Society of Visualization
• /
• v.9 no.4
• /
• pp.63-68
• /
• 2011

The flow-field of a liquid-metal system is very important for the safety analysis and the design of the steam generator of liquid-metal fast breeder reactor. Dynamic neutron radiography (DNR) is suitable for a visualization and measurement of a liquid metal flow and a two-phase flow in a metallic duct. However, the three dimensional DNR techniques is not enough to obtain the velocity information in the wide channel up to now. In this research, a high speed DNR technique was applied to visualize the heavy liquid-metal flow field in the narrow channel with the HANARO-beam facility. The images were taken with a high frame-rate neutron radiography at 250 fps and analyzed with a Particle Image Velocimetry(PIV) method. The images were compared with the results of the commercial CFX code to study the feasibility of DNR technique for the measuring the heavy liquid-metal flow field. The PIV images could discern the turbulent vortex flow in the two-dimensional narrow channel.

• Journal of Energy Engineering
• /
• v.15 no.1 s.45
• /
• pp.1-7
• /
• 2006

The multi-dimensional two-phase flow models were developed for analyze the multi-dimensional behaviors or nuclear systems. To verify the simple turbulence model, The single phase mixing problem in a rectangular slab was calculated and compared with the commercial CFD code results. That result shows a good agreement with the CFD result. And the RPI Air-water experiments were simulated to assess the two-phase turbulence model in the multi-dimensional component. The first calculated distribution or void-fraction is highly dispersed and diffusive. It was revealed that the main reason is undesirable stratification force in a horizontal stratified flow regimes. Therefore the horizontally stratified flow regime is deleted because the stratified flow regime is not expected in multi-dimensional flow. With the modification of the flow regime, the predicted flow patterns and void fraction profiles are in good agreement with the measured data.

• 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.

• Nuclear Engineering and Technology
• /
• v.38 no.8
• /
• pp.763-778
• /
• 2006

Injection of microbubbles within the turbulent boundary layer has been investigated for several years as a method to achieve drag reduction. However, the physical mechanism of this phenomenon is not yet fully understood. Experiments in a channel flow for single phase (water) and two phase (water and microbubbles) flows with various void fraction values are studied for a Reynolds number of 5128 based on the half height of the channel and bulk velocity. The state-of-the art Particle Tracking Velocimetry (PTV) measurement technique is used to measure the instantaneous full-field velocity components. Comparisons between turbulent statistical quantities with various values of local void fraction are presented to elucidate the influence of the microbubbles presence within the boundary layer. A decrease in the Reynolds stress distribution and turbulence production is obtained with the increase of microbubble concentration. The results obtained indicate a decorrelation of the streamwise and normal fluctuating velocities when microbubbles are injected within the boundary layer.

• Ocean Science Journal
• /
• v.43 no.2
• /
• pp.81-90
• /
• 2008

A study of estuarine flows during a neap tide was performed using 13-hour roving acoustic Doppler current profiles (ADCP) and conductivity-temperature-depth (CTD) profiles in the Altamaha River estuary, Georgia, U.S.A. The least-squared harmonic analysis method was used to fit the tidal ($M_2$) component and separate the flow into two components: the tidal and residual ($M_2$-removed) flows. We applied this method to depth-averaged data. Results show that the $M_2$ component demonstrates over 95% of the variability of observation data. As the flow was dominated by the $M_2$ tidal component in a narrow channel, the tidal ellipse distribution was essentially a back-and-forth motion. The amplitude of $M_2$ velocity component increased slightly from the river mouth (0.45 m/sec) to land (0.6 m/sec) and the phase showed fairly constant values in the center of the channel and rapidly decreasing values near the northern and southern shoaling areas. The residual flow and transport calculated from depth-averaged flow shows temporal variability over the tidal time scale. Strong landward flows appeared during slack waters which may be attributed to increased baroclinic forcing when turbulent mixing decreases.

• International Journal of Air-Conditioning and Refrigeration
• /
• v.16 no.3
• /
• pp.89-99
• /
• 2008

Indoor ozone has received attention because of its well-documented adverse effects on health. In addition to the inherently harmful effects of ozone, it can also initiate a series of reactions that generate potentially irritating oxidation products, including free radicals, aldehydes, organic acids and secondary organic aerosols (SOA). Especially, ozone reacts actively with terpene. The overarching goal of this work was to better understand ozone and terpene distributions within rooms. Towards this end, the paper has two parts. The first describes the development of a cylindrical test chamber that can be used to obtain the second order rate constant $(k_b)$ for the bi-molecular chemical reaction of ozone and terpene in the air phase. The second consists of model room experiments coupled with Computational Fluid Dynamics (CFD) analysis of the experimental scenarios to obtain ozone and terpene distributions in various turbulent flow fields. The results of CFD predictions were in reasonable agreement with the experimental measurements.