• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1996.06a
• /
• pp.179-200
• /
• 1996

The intrinsic instabilities of fluid flow occurred in the melt of the Czochralski crystal growth system Czochralski method, asymmetric flow patterns and temperature profiles in the melt have been studied by many researchers. The idea that the non-symmetric structure of the growing equipment is responsible for the asymmetric profiles is usually accepted at the first time. However further researches revealed that some intrinsic instabilities not related to the non-symmetric equipment structure in the melt could also appear. Ristorcelli had pointed out that there are many possible causes of instabilities in the melt. The instabilities appears because of the coupling effects of fluid flow and temperature profiles in the melt. Among the instabilities, the B nard type instabilities with no or low crucible rotation rates are analyzed by the visualizing experiments using X-ray radiography and the 3-D numerical simulation in this study. The velocity profiles in the Silicon melt at different crucible rotation rates were measured using X-ray radiography method using tungsten tracers in the melt. The results showed that there exits two types of fluid flow mode. One is axisymmetric flow, the other is asymmetric flow. In the axisymmetric flow, the trajectory of the tracers show torus pattern. However, more exact measurement of the axisymmetrc case shows that this flow field has small non-axisymmetric components of the velocity. When fluid flow is asymmetric, the tracers show random motion from the fixed view point. On the other hand, when the observer rotates to the same velocity of the crucible, the trajectory of the tracer show a rotating motion, the center of the motion is not same the center of the melt. The temperature of a point in the melt were measured using thermocouples with different rotating rates. Measured temperatures oscillated. Such kind of oscillations are also measured by the other researchers. The behavior of temperature oscillations were quite different between at low rotations and at high rotations. Above experimental results means that the fluid flow and temperature profiles in the melt is not symmetric, and then the mode of the asymmetric is changed when rotation rates are changed. To compare with these experimental results, the fluid flow and temperature profiles at no rotation and 8 rpm of crucible rotation rates on the same size of crucible is calculated using a 3-dimensional numerical simulation. A finite different method is adopted for this simulation. 50×30×30 grids are used. The numerical simulation also showed that the velocity and flow profiles are changed when rotation rates change. Futhermore, the flow patterns and temperature profiles of both cases are not axisymmetric even though axisymmetric boundary conditions are used. Several cells appear at no rotation. The cells are formed by the unstable vertical temperature profiles (upper region is colder than lower part) beneath the free surface of the melt. When the temperature profile is combined with density difference (Rayleigh-B nard instability) or surface tension difference (Marangoni-B nard instability) on temperature, cell structures are naturally formed. Both sources of instabilities are coupled to the cell structures in the melt of the Czochralski process. With high rotation rates, the shape of the fluid field is changed to another type of asymmetric profile. Because of the velocity profile, isothermal lines on the plane vertical to the centerline change to elliptic. When the velocity profiles are plotted at the rotating view point, two vortices appear at the both sides of centerline. These vortices seem to be the main reason of the tracer behavior shown in the asymmetric velocity experiment. This profile is quite similar to the profiles created by the baroclinic instability on the rotating annulus. The temperature profiles obtained from the numerical calculations and Fourier transforms of it are quite similar to the results of the experiment. bove esults intend that at least two types of intrinsic instabilities can occur in the melt of Czochralski growing systems. Because the instabilities cause temperature fluctuations in the melt and near the crystal-melt interface, some defects may be generated by them. When the crucible size becomes large, the intensity of the instabilities should increase. Therefore, to produce large single crystals with good quality, the behavior of the intrinsic instabilities in the melt as well as the effects of the instabilities on the defects in the ingot should be studied. As one of the cause of the defects in the large diameter Silicon single crystal grown by the

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.575-578
• /
• 2006

In designing smoke-control system of rescue station in train tunnel, a purpose is to prevent a disaster by proposing the jet fan operation together with smoke-control curtain in tunnel fire. This study has investigated the relationship of the Heat Release Rate(HRR) and a adequate ventilation velocity to control the fire propagation in tunnel fire, and has improved the effect of the smoke-control curtain on preventing the flow of pollutants. In this study, Computational Fluid Dynamics(CFD) simulations with ST AR-CD(ver 3.24) were carried out on predicting the fire spreading and the flow of pollutants, considering jet fan operations and effect of smoke-control curtain. Our simulation domain is the full scale model of the 'DAEGWALLYEONG' 1st tunnel. The results represent that ventilation operation can control the fire spreading and pollutants effectively to prevent a disaster.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.141-148
• /
• 2007

We present an immersed boundary (IB) method for 3D simulation of flappingflags in a uniform flow. The proposed formulation is manipulated on the basis of an efficient Navier-Stokes solver adopting the fractional step method and a staggered Cartesian grid system. A direct numerical method is developed to calculate the flag motion, with the elastic force treated implicitly. The fluid motion defined on an Eulerian grid and the flag motion defined on a Lagrangian grid are independently solved and the mass of flag is handled in a natural way. An additional momentum forcing is formulated from the flag motion equation in a way similar with the direct-forcing IB formulation and acts as the interaction force between the flag and ambient fluid. A series of numerical tests are performed and the present results are compared qualitatively and quantitatively with previous studies. The instantaneous flag motion is analyzed under different conditions and surrounding vortical structures are identified. The effects of physical parameters on the flapping frequency are studied.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.557-560
• /
• 2002

In this study, simulation of weak shock waves are peformed by a two-dimensional thermal fluid or compressible fluid model of the lattice Boltzmann method. The shock wave represents an abrupt change in fluids properties, in which finite variations in pressure, internal energies, and density occur over the shock thickness. The characteristics of the proposed model with a simple distribution function is verified by calculation of the sound speeds, and the shock tube problem. The reflection of a weak shock wave by wedge propagating in a channel is performed. The results agree well with those by finite difference method or by experiment. In the simulation of unsteady shock wave diffraction around a sharp corner, we show a flow field of vortical structure near the comer.

• Structural Engineering and Mechanics
• /
• v.31 no.3
• /
• pp.349-365
• /
• 2009

Membrane structures are quite sensitive to wind and therefore the fluid-solid interaction can not be neglected in dynamic analysis. A boundary element method (BEM) for 3D simulation of wind-structure interaction in tensile membrane structures is presented in this paper. The flow is treated as incompressible and potential. The flow field is solved with boundary element method codes and structural simulation is performed by finite element method software ANSYS. The nonlinear equations system is solved iteratively, with segregated treatment of the fluid and structure equations. Furthermore this method has been demonstrated to be effective by typical examples. Besides, the influence of several parameters on the wind-structure interaction, such as rise-span ratio, prestress and the wind velocity are investigated according to this method. The results provide experience in wind resistant researches and engineering.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.867-876
• /
• 2003

The Lagrangian dispersion of fluid particles in inhomogeneous turbulence is investigated by a direct numerical simulation of turbulent channel flow. Four points Hermite interpolation in the homogeneous direction and Chebyshev polynomials in the inhomogeneous direction is adopted to simulate the fluid particle dispersion. An inhomogeneity of Lagrangian statistics in turbulent boundary layer is investigated by releasing many particles at several different wall-normal locations and tracking those particles. The fluid particle dispersions and Lagrangian structure functions of velocity are scaled by the Kolmogorov similarity. The auto-correlations of velocity and acceleration are shown at the different releasing locations. Effect of initial particle location on the dispersion is analyzed by the probability density function at the several downstreams and time instants.

• Journal of Energy Engineering
• /
• v.24 no.2
• /
• pp.115-119
• /
• 2015

Pipe vibration caused great threat to the safety in production. Strong pipeline vibration will line accessories, especially the joints and pipe fittings etc. pipe joints loosening and rupture, causing serious accidents. By the action of the compressor constant fluid flow within the pipe, this process produces pulsating fluid flow may cause vibration of the pipe, thereby reducing the efficiency of the pipeline, structural vibration induced fatigue, thereby resulting in even piping structural damage. This paper studies on the vibration problems caused by fluid, by analyzing the causes of pipeline vibration and factors affecting pipeline vibrations, FEM (Finite Element Method) analysis of modal and enforced vibration.

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

Due to rapid progress in the performance of high-end computers, numerical prediction of fluid flow and flow-induced sound is expected to become a vital tool for aero- and hydro- dynamic design of various flow-related products. This presentation focuses on the applications of large-scale numerical simulations to complex engineering problems with a particular emphasis placed on the low-speed flows. Flow field computations are based on a large eddy simulation that directly computes all active eddies in the flow and models only those eddies responsible for energy dissipations. The sound generated from low-speed turbulent flows are computed either by direct numerical simulation or by decoupled methods, according to whether or not the feedback effects of the generated sound onto the source flow field can be neglected. Several numerical examples are presented in order to elucidate the present status of such computational methods and discussion on the future prospects will also be given.

• International Journal of Fluid Machinery and Systems
• /
• v.4 no.1
• /
• pp.209-216
• /
• 2011

The flow upstream of a centrifugal pump impeller has been investigated by both experimental test and numerical simulation. For experimental study, the flow field at four sections in the pump suction is measured by using PIV method. For calculation, the three dimensional turbulent flow for the full flow passage of the pump is simulated based on RANS equations combined with RNG k-$\varepsilon$ turbulence model. From those results, it is noted that at both design lo ad and quarter load condition, the pre-swirl flow whose direction is the same as the impeller rotation exists at all four sections in suction pipe of the pump, and at each section, the pre-swirl velocity becomes obviously larger at higher rotational speed. It is also indicated that at quarter load condition, the low pressure region at suction surface of the vane is large because of the unfavorable flow upstream of the pump impeller.

• Korea-Australia Rheology Journal
• /
• v.14 no.4
• /
• pp.161-174
• /
• 2002

A new technique for numerical calculation of viscoelastic flow based on the combination of Neural Net-works (NN) and Brownian Dynamics simulation or Stochastic Simulation Technique (SST) is presented in this paper. This method uses a "universal approximator" based on neural network methodology in combination with the kinetic theory of polymeric liquid in which the stress is computed from the molecular configuration rather than from closed form constitutive equations. Thus the new method obviates not only the need for a rheological constitutive equation to describe the fluid (as in the original Calculation Of Non-Newtonian Flows: Finite Elements St Stochastic Simulation Techniques (CONNFFESSIT) idea) but also any kind of finite element-type discretisation of the domain and its boundary for numerical solution of the governing PDE's. As an illustration of the method, the time development of the planar Couette flow is studied for two molecular kinetic models with finite extensibility, namely the Finitely Extensible Nonlinear Elastic (FENE) and FENE-Peterlin (FENE-P) models.P) models.