• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.160-161
• /
• 2003

Thermocapillary convection is a surface tension driven flow due to a temperature gradient along an interface. It occurs during a crystal-growth process and therefore understanding the convection is important to material processing in microgravity. Although modelling of the float-zone crystal-growth process has been of interest for a few decades, most studies of liquid bridges assumed non-deformable flat surfaces. In reality, the surface profile, g(t,z), is unknown and should be obtained as a solution to the coupled transport equations along with the surface force balance. Here we report on a numerical study of axisymmetric thermocapillary convection in liquid bridges with deformable surfaces. The interface is determined as part of the complete solution. The influence of the capillary number (Ca), Reynolds number (Re), Prandtl number (Pr) and aspect ratio(Ar) on the dynamics is explored.

• 한국결정성장학회:학술대회논문집
• /
• 한국결정성장학회 1996년도 The 9th KACG Technical Annual Meeting and the 3rd Korea-Japan EMGS (Electronic Materials Growth Symposium)
• /
• 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

• International Journal of Aeronautical and Space Sciences
• /
• 제15권1호
• /
• pp.91-96
• /
• 2014

The aim of this work is to determine supersonic nozzle profiles that are used in propulsion, for launchers or embarked with satellites. This design has the role of firstly, providing important propulsion, i.e. with uniform and parallel flow at exit; and secondly, to find short length profiles, without modification of the flow in the nozzle. The first elaborate program is used to determine the profile of the divergent, by using the characteristics method for an axisymmetric flow. The second program is conceived by using the finite volume method, to determine and test the profile found connected to a convergent.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 2006년 제4회 한국유체공학학술대회 논문집
• /
• pp.199-202
• /
• 2006

The objective of this study is to propose methods of controlling the wake behind a sphere for drag reduction using the suboptimal control theory and surrogate management framework, respectively. The Reynolds numbers considered is 300 at which the base flow is unsteady planar symmetric. Given the cost function defined as the square of the difference between the target pressure (potential-flow pressure) and real flow pressure on the sphere surface, the suboptimal control makes the flow steady axisymmetric and produces drag reduction. Based on the actuation profile from the suboptimal control, the optimal wavy actuation profile is obtained using the surrogate management framework and produces drag reduction.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집A
• /
• pp.470-474
• /
• 2000

Cold expansion of fastener holes is a mechanical process widely used in the aerospace industry. This treatment leads to an improvement of fatigue behavior due to the developed compressive residual stresses on the hole surface. The residual stress profile depends on the parameters of cold expansion, which are, expanding rate, inserting direction of mandrel, material properties etc. and the method to confirm this profile is only measurement by X-ray diffractometer. Despite its importance to aerospace industries, little attention has been devoted to the accurate modelling of the process. In this paper, Two-dimensional axisymmetric finite element simulations have been conducted for the cold expansion in an aluminium plate in order to predict the magnitude and distribution of the residual stress and plastic deformation. Maximum compressive residual stress could be increase about 7 percent using the 2-step cold expansion method.

• 한국전산유체공학회지
• /
• 제18권4호
• /
• pp.69-73
• /
• 2013

Mass manufacturing of optical fiber includes the process of very thin glass fiber drawing by heating and softening the high purity silica preform and applying the draw tension on the softened tip of preform neck-down profile in a draw furnace. In this computational study, this process is numerically modeled with simplified geometry of the draw furnace which is comprised of essential parts such as concentric graphite heater, muffle tube, and insulation surrounding the heater. The iterative computational scheme is employed between one-dimensional model of neck-down profile prediction and two-dimensional axisymmetric thermo-fluid CFD computation of radiative heating and working gas convection. The computational results show the experimentally observed neck-down profile in heated section of preform, while yielding the reasonable values of draw tension and heater wattage. Also, this study analyzes and discusses the effects of heating conditions such as heater length and temperature on several important aspects of glass fiber drawing process.

• 한국정밀공학회지
• /
• 제12권8호
• /
• pp.63-72
• /
• 1995

Closed-die forging of spur gears with hollow cylindrical billet has been analysed by using the upper-bound method. A kinematically admissible velocity field has been developed, wherein, an involute curve has been introduced to represent the forging die profile. In the analysis, the deformation region has been divided into nine zones. A constant frictional stress has been assumed on the contacting surfaces. Utilizing the formulated velocity field, numerical calculations have been carried out to investigate the effects of various parameters, such as module, number of teeth and friction factor, on the forging of spur gears. Hardness and accuracy of forged gears are measured. The following results have been obtained: (1) It is verified that an axisymmetric deformation zone exists between root circle and center of gear through forged gears. (2) The average relative forging pressure is predominantly dependent on the number of teeth and increases near the final filling stage as the addendum modification coefficient increases. (3) Close agreement was found between the predicted values of forging load and those obtained from experimental results.

• 한국결정성장학회:학술대회논문집
• /
• 한국결정성장학회 1997년도 Proceedings of the 12th KACG Technical Meeting and the 4th Korea-Japan EMGS (Electronic Materials Growth Symposium)
• /
• pp.197-201
• /
• 1997

The momentum, heat and mass trasfer phenomena in the silicon melt of the Czochralki crystal growth system are calculated using a three dimensional numerical simulation thechnique. Even though axisymmetrical boundary conditions are imposed to all calculations in a 3cm diameter crucible, several types of non-axisymmetric profiles of velocities, temperature and oxygen concentration appeared in the melt. Because of the non-axisymmetric profiles of velocities, temperature and oxygen concentration appeared in the melt. Because of the non-axisymmetric profiles and rotations of fluid induced by the crucible rotation, temperatures and oxygen concentrations in the silicon melt fluctuate. The rotating velocity of the profile is calculated from the phase shift of the data of temperature or oxygen at two different points which have same radius from center but 90 degree angular difference. From this calculation, it is found that the rotating veolocity of the oxygen and temperature is different from the crucible rotation rates. Therefore the frequencies of the oscillating temperature and oxygen concentrations are not same to the frequencies of the crucible rotations. Futhermore, the components of the frequencies of the temperature and oxygen concentration at the same point are not same. The fluctuation behaviors of the temperature or oxygen themselves are also different when the points are different. The calculation show that the temperature and the oxygen concentration near the interface also fluctuate. The results suggest that the striation pattern found in the grown silicon single crystals may ben generated by the oxygen concentration and the temperature oscillations of the melt occurred near the interface.

• 대한기계학회논문집B
• /
• 제21권11호
• /
• pp.1422-1436
• /
• 1997

Numerical investigations were conducted to study the convective phenomena of an initially stably stratified salt water solution with lateral heating in a uniformly rotating annulus. The method of investigation is the finite difference analysis of the basic conservation equation for an axisymmetric, unsteady, double-diffusive convection and calculation is made for R $a_{\eta}$=2*10$^{5}$ and Ta=10$^{7}$ ~ 2.5*10$^{8}$ . Formation of layered flow structure, merging process of layers, the corresponding temperature and concentration distributions, Nu variation with time are examined. Numerical results show that in each layer, the temperature profile looks 'S'-shaped and the concentration profile is uniform due to the convective mixing. At the interface between adjacent layers, the temperature changes smoothly but the concentration changes rapidly. As the effect of the rotation increases, the generation of rolls at hot wall, the formation and merging of layers are delayed. The average Nu shows the trend of conduction heat transferees the effect of the rotation increases.n increases.

• Korea-Australia Rheology Journal
• /
• 제21권1호
• /
• pp.39-46
• /
• 2009

Atherosclerosis is a disease that narrows, thickens, hardens, and restructures a blood vessel due to substantial plaque deposit. The geometric models of the considered stenotic blood flow are three different types of constriction of cross-sectional area of blood vessel; 25%, 50%, and 75% of constriction. The computational model with the fluid-structure interaction is introduced to investigate the wall shear stresses, blood flow field and recirculation zone in the stenotic vessels. The velocity profile in a compliant stenotic artery with various constrictions is subjected to prescribed physiologic waveform. The computational simulations were performed, in which the physiological flow through a compliant axisymmetric stenotic blood vessel was solved using commercial software ADINA 8.4 developed by finite element method. We demonstrated comparisons of the wall shear stress with or without the fluid-structure interaction and their velocity profiles under the physiological flow condition in the compliant stenotic artery. The present results enhance our understanding of the hemodynamic characteristics in a compliant stenotic artery.