• Nuclear Engineering and Technology
• /
• v.16 no.1
• /
• pp.1-10
• /
• 1984

General time-volume averaged conservation equations and jump conditions for two-phase flows are derived here. The time-averaged equations for a single phase region in two-phase flow are obtained from local instant balance equations by a technique often used for single phase turbulent flow equations. The results obtained by integrating the time averaged equations over a flow volume are spatially averaged twice; first, they are averaged over a single phase region of the k-th phase and then averaged over the total volume of the k-th phase, in a flow volume. The mass, momentum, and energy conservation equations are obtained from the general time-volume averaged equations. The advantages of the present model are explained by comparing it with Ishii's model (1) and Banerjee's model (2). Finally, the assumptions and approximate terms of the equations of the THERMIT-6S are clarified.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1327-1339
• /
• 1999

Characteristics of two-phase flow and heat transfer were numerically investigated in a submerged gas Injection system. Effects of both the gas flow rate and bubble size were investigated. In addition, heat transfer characteristic and effects of heat transfer were investigated when temperature of the injected gas was different from that of the liquid. The Eulerian approach was used for the formulation of both the continuous and the dispersed phases. The turbulence in the liquid phase was modeled by the use of the standard $k-{\varepsilon}$ turbulence model. The interphase friction and heat transfer coefficient were calculated by means of correlations available in the literature. The turbulent dispersion of the phases was modeled by introducing a "dispersion Prandtl number". The plume region and the axial velocities are increased with increases in the gas flow rate and with decreases in the bubble diameter. The turbulent flow field grows stronger with the increases in the gas flow rate and with the decreases in the bubble diameter. In case that the heat transfer between the liquid and the gas is considered, the axial and the radial velocities are decreased in comparison with the case that there is no temperature difference between the liquid and the gas when the temperature of the injected gas is higher than the mean liquid temperature. The results in the present research are of interest in the design and the operation of a wide variety of material and chemical processes.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05b
• /
• pp.241-246
• /
• 1996

An analytical model to calculate rate of vapor generation due to heterogeneous wall nucleation in flashing flow is developed. In the present model, an important parameter of the vapor generation term, i.e. nucleation site density is calculated by integrating its probability distribution function with respect to active cavity radius. The limits of integration are minimum and maximum active cavity radii, and these are formulated using an active cavity model for nucleate boiling. This formulation, therefore. can statistically account for the effect of surface specific thermo-physical and geometric conditions on the vapor generation rate and flashing inception. For verifying the adequacy of the present model, steady state two-fluid and the bubble transport equations are solved with applicable constitutive equations. The applicable region of the bubble transport equation is also extended to churn-turbulent flow regime to predict interfacial area concentration at high void fraction. Predicted results in terms of axial pressure and void fraction profiles along the channels are compared with experimental data of Super Moby Dick and BNL Reasonable agreements have been achieved and this shows the applicability of the present model to flashing flow analysis.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.11
• /
• pp.3639-3648
• /
• 1996

The effect of 15% pinched diffuser in a centrifugal air compressor with a cascade airfoil diffuser on the aerodynamic performance is investigated using a numerical approach. The commercial CFD code for three-dimensional, turbulent, compressible flow fields is executed for various mass flow rates at a design speed which can be obtained as long as the calculation succeeds. The pinched diffuser is found to help improve the instability of flow within vaneless diffuser space, especially the reverse flow near shroud, and to change both stall/surge line and choking line to increase the surge margin. It is also found to generate more favorable increase of static pressure in diffuser region, and to increase the resulting pressure ratio and efficiency.

• Journal of Advanced Marine Engineering and Technology
• /
• v.28 no.8
• /
• pp.1225-1238
• /
• 2004

The present study has been performed for obtaining the melting heat transfer enhancement characteristics of water mixture slurries of plural microcapsules having different diameters encapsulated with solid-liquid phase change material(PCM) flowing in a pipe heated under a constant wall heat flux condition. In the turbulent flow region, the friction factor of the present PCM slurry was to be lower than that of only water flow due to the drag reducing effect of the PCM slurry. The heat transfer coefficient of the PCM slurry flow in the pipe was increased by both effects of latent heat involved in phase change process and microconvection around plural microcapsules with different diameters. The experimental results revealed that the average heat transfer coefficient of the PCM slurry flow was about 2～2.8 times greater than that of a single phase of water.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.358-361
• /
• 2010

A Reynolds-Averaged Navier-Stokes (RANS) code with transition prediction model is developed and the computational results on an oscillating airfoil are compared with the experimental data for OA209 airfoil. An approximated eN method that can predict transition onset points and the length of transition region is directly applied to the RANS code. The hysteresis loop in dynamic stall is compared for the computational results using transition prediction and fully turbulent models with the experimental data. Results with transition prediction show more correlation with the experimental data than the fully turbulent computation.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.2
• /
• pp.91-98
• /
• 1993

A two-dimensional turbulence model for the surface discharge of heated water into cross-flow field has been developed. The depth-averaged continuity, momentum and temperature equations, are solved by an efficient finite-difference procedure known as SIMPLE. Turbulent stresses and heat fluxes are determined from a depth-averaged version of the $textsc{k}$-$\varepsilon$ equation. Results of test run clearly demonstrate its effectiveness in handling strong turbulent phenomena in very shallow near-field region.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.39-45
• /
• 2010

The turbulent flows in a tunnel mock-up($10L{\times}0.5W{\times}0.25H$ m3 : scale reduction 1/20) with rectangular cross section were investigated. The instantaneous velocity fields of Re = 49,029, 89,571 were measured by the 2-D PIV system which is consisted of double pulsed Nd:Yag laser and the tracer particles in the straight-duct mock-up where the flows were fully developed. The mean velocity profiles were taken from the ensemble averages of 1,000 instantaneous velocity fields. Simultaneously, numerical simulations(RANS) were performed to compare with experimental data using STREAM code. Non-linear eddy viscosity model (NLEVM : Abe-Jang-Leschziner Eddy Viscosity Model) was employed to resolve the turbulent flows in the duct. The calculated mean velocity profiles were well compared with PIV results. In the log-law profiles, the experimental data were in good agreement with numerical simulations all the way to the wake region except the viscous sub-layer (near wall region).

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.2
• /
• pp.470-483
• /
• 1990

Turbulent jet in a crossflow, issuing from a row of holes on a convex surface of 90 .deg. bend duct, is predicted by a 3-dimensional numerical method. The Cartesian coordinate system in adopted in upstream and downstream tangents and the cylindrical polar coordinate system in curved region. The Reynolds stresses and heat fluxes are obtained from a standard k-e model in the core region and van Driest model in the vicinity of the wall. The governing equations are discretized by a finite volume method and solutions are obtained by a locally elliptic calculation procedure. Pressure and convective terms are treated by SIMPLE algorithm and hybrid scheme respectively. A vortex initially induced by the injected jet has been built up due to the interaction with the secondary flow caused by pressure gradient and centrifugal force. The vortex structure has a strong influence on the wall cooling effectiveness. Another vortex like horseshoe is formed in the vicinity of the injection hole and its strength is getting weak as it moves downward.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.8
• /
• pp.1927-1939
• /
• 1995

An experimental investigation is conducted to find out the large scale coherent structures in the intermediate wake past a rotating cylinder with a single tripping wire attached. Relation between the vortex shedding frequency and the spin rate of rotating cylinder and the effects of the tripping wire on the flow characteristics were studied by using spectral analysis and conditional phase average technique, respectively. It is found that the vortex shedding frequency is bound to a certain range and varies regularly as spin rate increases. The coherent structures are compared with those of the plain rotating cylinder in the case of spin rate of 1.0. Distance between the upper and lower center of vortices increase and the vortex shedding time is delayed, the velocity fluctuation energy decreases near the center line of vortices and it spreads out to the outer region. The Reynolds shear stress increases highly in the upper region and the turbulent wake width expands with strong entrainment process.