• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2294-2303
• /
• 1993

This paper investigated the characteristics of the turbulent incompressible flow past the orifice ring in an axi-symmetric pipe. The flow field was the turbulent pulsatile flow for Reynolds number of $2{\times}10^{5}$ which was defined based on the maximum velocity and the pipe diameter at the inlet, with oscillating frequence $(f_{os})=1/4{\pi}$ which was considered as quasi-steady state frequence. In the present investigation, finite analytic method was used to solve the governing equations in Navier Stokes and turbulent transport formulations. Particularly at high Reynolds number and low oscillation frequency, the effects of orifice ring on the flow were numerically investigated. The separation zone behind the orifice ring during the acceleration phase was found to be decreased. However, during the deceleration phase, the separation behind the orifice ring for pulsatile flow continuously grow to a size even larger than that in steady flow. The pressure drop in steady flow was found to be constant and always positive while for pulsatile flow the pressure drop change with time. And large turbulent kinetic energy, dissipation rate were found to be located in the region where the flow passes through the orifics ring. The maximum turbulent kinetic energy, generally occurs along the shear layer where the velocity gradient is large.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.62-68
• /
• 2008

This paper reports briefly on the computational results of a turbulent Rayleigh-Benard convection with the elliptic-blending second-moment closure (EBM). The primary emphasis of the study is placed on an investigation of accuracy and numerical stability of the elliptic-blending second-moment closure for the turbulent Rayleigh-Benard convection. The turbulent heat fluxes in this study are treated by the algebraic flux model with the temperature variance and molecular dissipation rate of turbulent heat flux. The model is applied to the prediction of the turbulent Rayleigh-Benard convection for Rayleigh numbers ranging from Ra=$2{\times}10^6$ to Ra=$10^9$ and the computed results are compared with the previous experimental correlations, T-RANS and LES results. The predicted cell-averaged Nusselt number follows the correlation by Peng et al.(2006) (Nu=$0.162Ra^{0.286}$) in the 'soft' convective turbulence region ($2{\times}10^6{\leq}Ra{\leq}4{\times}10^7$) and it follows the experimental correlation by Niemela et al. (2000) (N=$0.124Ra^{0.309}$) in the 'hard' convective turbulence region ($10^8{\leq}Ra{\leq}10^9$) within 5% accuracy. This results show that the elliptic-blending second-moment closure with an algebraic flux model predicts very accurately the Rayleigh-Benard convection.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.1063-1071
• /
• 1999

Fine grid calculations are reported for the developing turbulent flow in a curved duct of square cross-section with a radius of curvature to hydraulic diameter ratio ${\delta}=Rc/D_H=3.357 $ and a bend angle of 720 deg. A sequence of modeling refinements is introduced; the replacement of wall function by a fine mesh across the sublayer and a low Reynolds number algebraic second moment closure up to the near wall sublayer in which the non-linear return to isotropy model and the cubic-quasi-isotropy model for the pressure strain are adopted; and the introduction of a multiple source model for the exact dissipation rate equation. Each refinement is shown to lead to an appreciable improvement in the agreement between measurement and computation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.1022-1030
• /
• 1999

A numerical study was carried out to analyze the effect of flow distribution of stirred part and plug flow part on combustion efficiency at the coal gasification process in an entrained bed coal reactor. The model of computation was based on gas phase eulerian balance equations of mass and momentum. The solid phase was described by lagrangian equations of motion. The $k-{\varepsilon}$ model was used to calculate the turbulence flow and eddy dissipation model was used to describe the gas phase reaction rate. The radiation was solved using a Monte-Carlo method. One-step parallel two reaction model was employed for the devolatilization process of a high volatile bituminous Kideco coal. The computations agreed well with the experiments, but the flame front was closer to the burner than the measured one. The flow distribution of a stirred part and a plug flow part in a reactor was a function of the magnitude of recirculation zone resulted from the swirl. The combustion efficiency was enhanced with decreasing stirred part and the maximum value was found around S=1.2, having the minimum stirred part. The combustion efficiency resulted from not only the flow distribution but also the particle residence time through the hot reaction zone of the stirred part, in particular for the weak swirl without IRZ(internal recirculation zone) and the long lifted flame.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.2
• /
• pp.243-250
• /
• 2003

The extinction behavior and the unsteady response of augmented reduced mechanism(ARM) have been investigated by adopting an OPPDIF code and a numerical solver for the flamelet equations. By comparing the performance of the ARM based on Miller and Bowman's mechanism(MB-ARM) with that of the ARM based on GRI-Mech 3.0(GRI-3.0-ARM), it is identified that the MB-ARM is more suitable for the unsteady calculation because it is relatively less stiff than GRI-3.0-ARM during an ignition process. The steady results using the MB-ARM, which is modified to predict reasonably the extinction point of experiment, are in excellent agreement with those from full mechanism. Under the sinusoidal transient disturbances of scalar dissipation rate, the unsteady responses of the flame temperature and species concentrations using a modified MB-ARM show in very close agreement with those from full mechanism. It is presumed that above modified MB-ARM is very suitable for the unsteady simulation of turbulent flames because it gives not only a low computational cost but also a good prediction performance for flame structure, extinction point and unsteady response.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.21-27
• /
• 2012

Reynolds-averaged Navier-Stokes equations solver based on a pressure-based cell-centered finite volume method was developed using OpenFOAM libraries, which was an open source and providing computational continuum mechanics libraries. For the reasonable development of the turbulent boundary layer on the bow of the ship, specified library was developed. Grid sensitivities, such as skewness and aspect ratio of a cell, were tested for the solution convergence. Pressure, turbulent kinetic energy, turbulent dissipation rate contours on the ship surface computed by the developed CFD code were compared with those computed by the commercial CFD code, Fluent.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1368-1379
• /
• 2000

A numerical study on swirling pulverized coal combustion in an axisymmetric enclosure is carried out by applying the 2-phase weighted sum of gray gases model (WSGGM) approach with the discrete ordinate method (DOM) to model the radiative heat transfer equation. In the radiative transfer equation, the same polynomial equation and coefficients for weighting factors as those for gas are adopted for the coal/char particles as a function of partial pressure and particle temperature. The Eulerian balance equations for mass, momentum, energy, and species mass fractions are adopted with the standard and RNG k-${\varepsilon}$ turbulence model, whereas the Lagrangian approach is used for the particulate phase. The eddy-dissipation model is employed for the reaction rate for gaseous mixture, and the single-step and two-step first-order reaction model for the devolatilization process for coal. Special attention is given to establish the thermal boundary conditions on radiative transfer equation By comparing the numerical results with experimental ones, the radiation model used here is confirmed and found to provide an alternative for simulating the radiative transfer.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.7 no.2
• /
• pp.193-203
• /
• 1983

The mixed convective heat transfer problems are characterized by the relatively significant contribution of buoyancy force to the transport processes of momentum and heat. Past analytical studies on this kind of problems have been carried out by employing either the conventional R-.epsilon. turbulence model which includes constant turbulent Prandtl number .sigma.$_{+}$ 1 or an extended R-.epsilon. turbulence model which takes account of the buoyancy effect in appropriate length scale equations. But in the latter case, the temperature variance .the+a.$^{2}$ over bar is approximated by a model under local equilibrium condition and the time scale ratio between velocity and temperature is assumed to be constant. These approximation is known to break down when the buoyancy effect is dominant. The present study is aimed at development of new computational turbulence closure level which can be applied to this rather complex turbulent process. The temperature variance is obtained directly by solving its dynamic transport equation and the time scale ratio which is variable in space is computed by a solution of a dynamic equation for the rate of scalar dissipation .epsilon.$_{\thetod}$ It was found that the computational results are in good agreement with available experimental data of wide range of unstable conditions.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.14 no.3
• /
• pp.38-42
• /
• 2000

Plasma polymerized thin films were manufactured inter-electrode coupled plasma polymerization apparatus. The deposition rate reached its maximum between 40[W] and 100[W]. In the ESCA analysis, peaks revealing -CH2, -CH, -C- were present at 285.4 and 285.5[eV] respectively. The C=O peak at 532.8[eV] and the C-O peak at 533.8[eV], which were grouped with an unignorable amount of oxygen were conformed. In ESR analysis, the curve revealing strong amplification was in saturation, which was affected by weak power. This is considered as a -CH-Ch=Ch- structure containing the Allyl group. The relative permittivity of the plasma polymerized thin films was about 3.5 at a frequency of 100[Hz]∼200[kHz]. The dissipation factor showed allow value of 0.008.

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.52-60
• /
• 2002

A two-dimensional direct numerical simulation is performed to investigate the flame structure of $CH_4/N_2-Air$ counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed reaction mechanism are adopted in this calculation. To quantify the strain on flame induced by a vortex, a scalar dissipation rate (SDR) is introduced. Results show that the fuel and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex is extinguished at much larger SDR than steady flame. It is also found that air-side vortex extinguishes a flame more rapidly than fuel-side vortex. The unsteady effect induced by flame-vortex interaction does not lead to a transient OH overshoot of the maximum steady concentration observed in experiment, while $HO_2$ radical increases more than the maximum steady concentration with increasing SDR. In addition, it is seen that NO and $NO_2$ are not sensitive to the unsteady variation of SDR.