• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.819-822
• /
• 2002

It is necessary for the numerical simulation of 3-dimensional incompressible isotropic decaying turbulence to construct 3-dimensional initial velocity field which resembles the fully developed turbulence. Although the previous velocity field generation method proposed by Rogallo(1981) satisfies continuity equation and 3-dimensional energy spectrum, it has limitation, as indicated in his paper, that it does not produce the higher velocity moments(e. g. velocity derivative skewness) characteristic of real turbulence. In this study, a new velocity field generation method which is able to control velocity derivative skewness of initial velocity field is proposed. Brief descriptions of the new method and a few parameters which is used to control velocity derivative skewness are given. A large eddy simulation(LES) of isotropic decaying turbulence using dynamic subgrid-scale model is carried out to evaluate the performance of the initial velocity field generated by the new method. It was shown that the resolved turbulent kinetic energy decay curve and the resolved enstrophy decay curve from the initial field of new method were more realistic than those from the initial field of Rogallo's method. It was found that the dynamic model coefficient from the former was initially half the stationary value and experienced relatively short transition period, though that from the latter was initially zero and experienced relatively longer transition period.

• Journal of Mechanical Science and Technology
• /
• 제14권12호
• /
• pp.1386-1395
• /
• 2000

The mear field structure of round turbulent jets with initially asymmetric velocity distributions is investigated experimentally. Experiments are carried out using a constant temperature hot-wire anemometry system to measure streamwise velocity in the jets. The measurements are undertaken across the jet at various streamwise stations in a range starting from the jet exit plane and up to a downstream location of twelve diameters. The experimental results include the distributions of mean and instantaneous velocities, vorticity field, turbulence intensity, and the Reynolds shear stresses. The asymmetry of the jet exit plane was obtained by using circular cross-section pipes with a bend upstream of the exit. There pipes used here include a straight pipe, and 90 and 160 degree-bend pipes. Therefore, at the upstream of the upstream of the pipe exit, secondary flow through the bend mean streamwise velocity distribution could be controlled by changing the curvature of pipes. The jets into the atmosphere have two levels of initial velocity skewness in addition to an axisymmetric jet from a straight pipe. In case of the curved pipe, a six diameter-long straight pipe section follows the bend upstream of the exit. The Reynolds number based on the exit bulk velocity is 13,400. The results indicate that the near field structure is considerably modified by the skewness of an initial mean velocity distribution. As the skewness increases, the decay rate of mean velocity at the centerline also increases.

• 천문학회지
• /
• 제37권4호
• /
• pp.237-241
• /
• 2004

We perform numerical experiments on supernova-driven turbulent flows in order to see whether or not supernovae playa major role in driving turbulence in the interstellar medium. In a $(200pc)^3$ computational box, we set up, as initial conditions, uniformly magnetized gas distributions with different pairs of hydrogen number densities and magnetic field strengths, which cover the observed values in the Galactic midplane. We then explode supernovae at randomly chosen positions at a Galactic explosion rate and follow up the evolution of the supernova-driven turbulent flows by integrating numerically the ideal MHD equations with cooling and heating terms. From the numerical experiments we find that the density-weighted velocity dispersions of the flows are in the range of 5-10 km $s^{-l}$, which are consistent with the observed velocity dispersions of cold and warm neutral media. Additionally, we find that strong compressible flows driven by supernova explosions quickly change into solenoidal flows.

• Journal of Mechanical Science and Technology
• /
• 제15권1호
• /
• pp.117-124
• /
• 2001

This paper presents turbulent characteristics of an impinging F-O-O-F type injector in which fuel ad oxidizer impinge on each other to atomize under the different momentum ratio. Water was used as an inert simulant liquid instead of fuel and oxidizer. The droplet size and velocity in the impinging spray flow field were measured using a PDPA. The gradient of the spray half-width(b$_2$) along the long-axis direction declined throughout the entire spray flow field with increasing the momentum ratio from 1.19 to 6.48. However, the gradient of the half-width(b$_1$) along the short-axis direction decreased with increasing the momentum ratio. The turbulence intensity and turbulent kinetic energy were converged into the center of the center of the initial region with increasing the momentum ratio. As the momentum ratio increased from MR=1.19 to MR=6.48, the turbulent shear stress decreased. The results of this study can be used for the design of an impinging type injector for liquid rackets.

• 대한기계학회논문집
• /
• 제9권4호
• /
• pp.403-413
• /
• 1985

본 논문에서는 수직벽 하류에 형성되는 박이 전단층의 발전과 재부착 그리고 재발전 경계층에 대해 평균 속도, 벽면의 압력 분포, 난류 강도, 레이놀즈 전단 압력 및 아직 수직벽에 대해서는 보고된바 없는 난류 떨림 속도 성분들의 3승곱 통계치를 측정하여 난류 구조의 변화를 분석하고 이를 수치적 계산 모델개발의 자 료로 제공하고자 함이 이 연구의 목적이다.

• Nuclear Engineering and Technology
• /
• 제55권12호
• /
• pp.4561-4569
• /
• 2023

A correct understanding of vibration-based degradation is crucial from the standpoint of maintenance for Steam Generators (SG) as crucial mechanical equipment in nuclear power plants. This study has established a novel approach to developing a model for investigating tube bundle degradation according to crack growth caused by two-phase Flow-Induced Vibration (FIV). An important step in the approach is to calculate the two-phase flow field parameters between the SG tube bundles in various zones using the porous media model to determine the velocity and vapor volume fraction. Afterward, to determine the vibration properties of the tube bundles, the Fluid-Solid Interaction (FSI) analysis is performed in eighteen thermal-hydraulic zones. Tube bundle degradation based on crack growth using the sixteen most probable initial cracks and within each SG thermal-hydraulic zone is performed to calculate useful lifetime. Large Eddy Simulation (LES) model, Paris law, and Wiener process model are considered to model the turbulent crossflow around the tube bundles, simulation of elliptical crack growth due to the vibration characteristics, and estimation of SG tube bundles degradation, respectively. The analysis shows that the tube deforms most noticeably in the zone with the highest velocity. As a result, cracks propagate more quickly in the tube with a higher height. In all simulations based on different initial crack sizes, it was observed that zone 16 experiences the greatest deformation and, subsequently, the fastest degradation, with a velocity and vapor volume fraction of 0.5 m/s and 0.4, respectively.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 1998년도 제17회 KOSCI SYMPOSIUM 논문집
• /
• pp.139-154
• /
• 1998

Under certain circumstance, premixed turbulent flame can be treated as wrinkled thin laminar flame and its motion in a hydrodynamic flow field has been investigated by employing G-equation. Past studies on G-equation successfully described certain aspects of laminar flame propagation such as effects of stretch on flame speed. In those studies, flames were regarded as a passive interface that does not influence the flow field. The experimental evidences, however, indicate that flow field can be significantly modified by the propagation of flames through the volume expansion of burned gas. In the present study, a new method to be used with G -equation is described to include the effect of volume expansion in the flame dynamics. The effect of volume expansion on the flow field is approximated by Biot-Savart law. The newly developed model is validated by comparison with existing analytical solutions of G -equation to predict flames propagating in hydrodynamic flow field without volume expansion. To further investigate the influence of volume expansion, present method was applied to initially wrinkled or planar flame propagating in an imposed velocity field and the average flame speed was evaluated from the ratio of flame surface area and projected area of unburned stream channel. It was observed that the initial wrinkling of flame cannot sustain itself without velocity disturbance and wrinkled structure decays into planar flame as the flame propagates. The rate of decay of the structure increased with volume expansion. The asymptotic change in the average burning speed occurs only with disturbed velocity field. Because volume expansion acts directly on the velocity field, the average burning speed is affected at all time when its effect is included. With relatively small temperature ratio of 3, the average flame speed increased 10%. The combined effect of volume expansion and flame stretch is also considered and the result implied that the effect of stretch is independent of volume release.

• 대한기계학회논문집B
• /
• 제22권12호
• /
• pp.1669-1680
• /
• 1998

A method is developed to include the effect of volume expansion in the description of the flame dynamics using G-equation. Line volume-source is used to represent the effect of the exothermic process of combustion with source strength assigned by the density difference between the burned and the unburned region. The present model provides good agreement with the experimental results. Including volume expansion, the flow field is adjusted to accommodate the increased volume flow rate which crossing the flame front and the result predicts the same behavior of measured velocity field qualitatively. The effect of increasing volume expansion does not change the initial growth rate of flame area but increase the residence time. Consequently this effect increases the maximum area of flame front. The flame propagation in varying flow field due to volume expansion provides a promising way to represent the wrinkled turbulent premixed flames in a numerically efficient manner.

• 한국안전학회지
• /
• 제37권2호
• /
• pp.1-9
• /
• 2022

Most factories deal with toxic or flammable chemicals in their industrial processes. These hazardous substances pose a risk of leakage due to accidents, such as fire and explosion. In the event of chemical release, massive casualties and property damage can result; hence, quantitative risk prediction and assessment are necessary. Several methods are available for evaluating chemical dispersion in the atmosphere, and most analyses are considered neutral in dispersion models and under far-field wind condition. The foregoing assumption renders a model valid only after a considerable time has elapsed from the moment chemicals are released or dispersed from a source. Hence, an initial dispersion model is required to assess risk quantitatively and predict the extent of damage because the most dangerous locations are those near a leak source. In this study, the dispersion model for initial consequence analysis was developed with three-dimensional unsteady advective diffusion equation. In this expression, instantaneous leakage is assumed as a puff, and wind velocity is considered as a coordinate transform in the solution. To minimize the buoyant force, ethane is used as leaked fuel, and two different diffusion coefficients are introduced. The calculated concentration field with a molecular diffusion coefficient shows a moving circular iso-line in the horizontal plane. The maximum concentration decreases as time progresses and distance increases. In the case of using a coefficient for turbulent diffusion, the dispersion along the wind velocity direction is enhanced, and an elliptic iso-contour line is found. The result yielded by a widely used commercial program, ALOHA, was compared with the end point of the lower explosion limit. In the future, we plan to build a more accurate and general initial risk assessment model by considering the turbulence diffusion and buoyancy effect on dispersion.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2002년도 추계 학술대회논문집
• /
• pp.133-138
• /
• 2002

A 3-D axisymmetric computer program is developed to predict the NO behavior in SNCR system for the stoker incinerator with the waste treatment capacity, 200ton/day. To this end a turbulent reacting flow field calculation is made using proper assumption and empiricism. The stoker bed is assumed to be a homogeneous waste-volatilized gaseous state. The initial composition or reactants are assumed based on the data of the ultimate analysis. Turbulent is resolved by k-e model and turbulent reaction is handled by eddy-breakup model harmonized with empirical chemistry data for gaseous combustion, NO and urea reaction. The liquid droplet is traced by Lagrangian method incorporated by aerodynamic drag, Coriolis and crntrifugal forces. Radiation is treated by sensible heat loss model. Calculation results are in good agreement with experimental data at the outlet of post combustion chamber in Daejon 4th industrial complex. The flue gas shows the temperature range of $900\sim1000^{\circ}C$, velocity of 5m/s and NO concentration of 140ppm at the exit while the measured temperature, flue gas velocity and NO concentration are $967^{\circ}C$, $3\sim4m/s$ and $100\sim200ppm$respectively. Using the developed computer program a parametric study has been made with the variation of heat content of waste, castable length and SNCR variables for the determination of proper injector location. In general, the calculated results are consistent and physically acceptable.