• Proceedings of the Korean Nuclear Society Conference
• /
• 1998.05a
• /
• pp.520-525
• /
• 1998

난류혼합율에 대한 예측은 원자로의 노심 열수력 설계에 있어 매우 중요한 일이다. 봉다발 구조에서 난류혼합의 주요 원인으로 지목되고 있는 유동액동(flow pulsation) 현상에 대한 척도평가(scale analysis)틀 통해 봉다발 유동장을 흐르는 저 Prandtl 수 유채에 대판 난류혼합율 평가식을 유도하였다. 난류혼합에 기여하는 인자가 분자운동, 등방성 난류운동(유동맥동 효과률 배제한 난류운동), 그리고 유동맥동의 세 부분으로 구성되어 있다고 가정하고, 각각에 대한 길이 및 속도척도를 평가하여 난류혼합율을 유도하였다. 평가식에는 P/D, Re수 P${\gamma}$ 수 등의 인자가 고려되어 있어 다양한 기하학적, 수력학적 조건과 유체의 물리적 특성이 반영되어 있다. 유도원 난류혼합율 평가식을 실험 상관식과 비교하였으며, 비교 결과 만족스러운 것으로 나타났다.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.13-21
• /
• 2013

Cavitating flow is usually formed on the surface of a high speed underwater object. When a object moves near a free surface at very high speed, the cavity signature becomes one of the major factors to be overcome by sensors of military satellite. The present work was to study the free surface effect on the ventilated cavitation. The governing equations were Navier-Stokes equations based on a homogeneous mixture model. The multiphase flow solver used an implicit preconditioning method in the curvilinear coordinate system. The cavitation model used here was the one first presented by Merkle et al.(2006) and redeveloped by Park & Ha(2009). Computations considered the free surface effects were carried out with a NACA0012 hydrofoil and the corresponding results were compared with the experimental data to have a good agreement. Calculations were then performed considering the ventilated cavitation, including the effect of non-condensable gas under the free surface effects.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.26 no.4
• /
• pp.27-34
• /
• 1989

A potential-based panel method is formulated for the analysis of a partially cavitating 2-dimensional hydrofoil. The method employs dipoles and sources distributed on the foil surface to represent the lifting and cavity problems, respectively. The kinematic boundry condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the inner flow region of the foil. The dynamic boundary condition on the cavity surface is satisfied by requiring that the potential vary linearly, i.e., the velocity be constant. Green's theorem then results in a potential-based boundary value problem rather than a usual velocity-based formulation. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with more improved accuracy than the zero-thickness hydrofoil theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. It was found that five iterations are necessary to obtain converged values, while only two iterations are sufficient for engineering purpose.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.2
• /
• pp.156-162
• /
• 1995

A numerical method for a 2-D oil boom model considering the flexibility of skirt has been developed The neater is assumed rigid and the skirt is tensioned membrane having a point mass at its end The fluid motion is potential. The kinematic condition which demands the continuity of the displacement is imposed at the joint between the floater and the skirt. The dynamic condition for the point mass is imposed at the bottom end of the skirt. The numerical method is based on the Green's function method in the frame of linear potential theory. It finds it's solution simultaneously from the total system of three equations, integral equation, the equation of motion of the floater and the equilibrium equation of the deformation of the skirt. Integral equation is derived by applying the Green's theorem to radiation potential and Green's function. Proper descretization of those three equations leads to the system of a linear algebraic equation. Due to the flexibility of skirt the motion of floater can be diminished in some range of wave frequency and furthermore the mechanism of resonance of the oil boom can be changed. The motion responses of various oil booms have been compared varying the length of the skirt and the point mass. The numerical method has been validated indirectly from the good correspondence between the motion responses of the flexible skirt model and the rigid skirt model in low frequency limit.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.702-715
• /
• 1990

This paper describes behaviors of two-phase open channel flow inside the flash chamber of a horizontal Multi-Stage-Flash evaporator numerically along with the experimental observations. Bubble trajectories and the velocity and temperature distributions of the liquid phase were predicted by using the particle-source-in-cell(PSI-Cell) method with the appropriate bubble motion/growth equations. Size and number of bubble nuclei embedded in the incoming liquid(brine) were taken into account as important parameters in addition to the conventional ones such as the velocity, degree of inlet superheat, inlet opening height, and the liquid level. Bubble motions, which are unsteady, appeared to be mostly determined by the buoyancy and the drag forces. The calculations, though a number of simplifying assumptions were made, reasonably simulated the hydrodynamic behaviors of the two-phase horizontal stream observed in the experiments. The simulated temperature distributions also agreed fairly well with the other's measurements. Non-equilibrium allownaces, evaluated from the simulated temperature distributions, were within the range of those obtained from the existing correlations, and reduced with the increases of the number and size of incoming bubble nuclei due to vigorous flashing.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.21 no.1
• /
• pp.45-50
• /
• 2013

In this study, lumped-vortex element method and thin airfoil theory were used to analyze aerodynamic characteristics of airfoils with relative motion that had camber lines of NACA $44{\times}{\times}$ airfoil in 2-dimensional unsteady incompressible potential flow. Velocity disturbance due to airfoil was calculated by lumped-vortex element model and force distribution on airfoil by unsteady Bernoulli's equation. Variables in relative motion were considered the period p, the amplitude of flapping $A_f$ and pitching $A_p$, and the phase difference between flapping and pitching ${\phi}_p$ and the angle of attack ${\alpha}$. Due to movement of an airfoil, dag was induced in 2-dimensional unsteady incompressible potential flow. The numerical results show that the aerodynamic characteristics of the airfoil with flapping and pitching at the same time are illustrated. Especially the mean lift coefficient became smaller, but drag coefficient became larger.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.766-774
• /
• 2017

Software has been developed for simulating particle dispersion in a circular pipe with sudden-expansion, which models the fuel feeding system of a combustor that uses metal powder like aluminum as fuel. The Lagrangian based discrete tracer point method was employed for a plug flow of particles that satisfies local turbulent velocity fluctuations. A radial velocity component was created to improve the flow turning outwards in the recirculation zone. The particle distribution patterns from both with and without the component were directly compared with the experiments near the reattachment.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.1189-1196
• /
• 1988

Thermal entrance lengths of turbulent tube flow for viscoelastic polymer solutions are investigated experimentally in the recirculating flow system with tubes of inside diameters 8.5mm(L/D=710) and 10.3mm(L/D=1158), respectively. In the present system, the hydrodynamic and thermal boundary layers develop simultaneously from the beginning of the test section. To provide the boundary condition of constant heat flux at the wall, the test tubes are heated directly by electricity. The polymer solution used in the current study is 1000 wppm aqueous solution of polyacrylamide(Separan AP-273). The apparent viscosity of the polymer solutions circulating in the flow system are measured by the capillary tube viscometer at regular time intervals. Thermal entrance lengths vary due to the rate of degradation. The entrance lengths of degraded polymer solutions are about 500~600 times the diameter. However, the entrance lengths of fresh polymer solutions are greater than the lengths of the test tubes used in this study suggesting that thermal entrance lengths for viscoelastic polymer solutions are greater than 1100 tube times the diameters. Friction factor is almost insensitive to the degradation, but the heat transfer $j_{H}$-factor is affected seriously by degradation. Based on the present experimental data of fresh solutions a correlation for the heat transfer $j_{H}$-factor is presented.ted.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.1 no.1
• /
• pp.15-22
• /
• 2000

There are lots of industrial machines of which functions are achieved by operation of multi-phase fluids. Some of them take advantage of the characteristics of counter-current two-phase flow The maximum flow rates of gas and liquid phases which flow in opposite-directions (counter-current flow) are limited by a phenomenon known as a Counter-Current Flow Limitation (CCFL or Flooding) The mass and momentum conservation equations for each Phase were established to build a first-order hyperbolic partial derivative equations system. A new CCFL model is developed based on the characteristic equation of the hyperbolic PDE system. The present model has its applicationto the case in which a non-uniform flow is developed around a square or sharp-edged entrance of liquid phase. The model is able to he used to Predict the operating-limit of components in which mass and heat transfer are taking place between liquid and gas phases.