• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.3
• /
• pp.62-74
• /
• 1993

The vortex shedding from a circular cylinder placed in a steady uniform stream is simulated by the vortex cloud model of the discrete vortex method. The vorticity created at the cylinder surface is discretely represented by a number of nascent vortices at each time step and the motion of these cumulative vortices is monitored to produce the evolution of the vortex distribution pattern. Convection of vortices was traced by the vortex-in-cell technique and the force coefficients were calculated by both Sarpkaya's formulae and Lee's formulae for comparison. Discussions concerning the interrelation between the computational parameters and some principles for choosing the suitable values are included.

• Solar Energy
• /
• v.12 no.3
• /
• pp.70-83
• /
• 1992

The interaction between the surface radiation and the mixed convection transport from an isolated thermal source, with a uniform surface heat flux input and located in a rectangular enclosure, is stuied numerically. The enclosure simulates a practical system such an air cooled electric device, where an air-stream flows through the openings on the two vertical walls. The heat source represents an electric component located in such an enclosure. The size of this cavity is $0.1[m]{\times}0.1[m]$. The inlet velocity is assumed as 0.07[m/s] and the inlet temperature is maintained as $27^{\circ}C$. The inflow is kept at a fixed position. Laminar, two dimensional flow is assumed, and the problem lies in the mixed convection regime, governed by buoyancy force and surface readiation. The significant variables include the location of the out-flow opening, of the heat source and the wall emissivity. The basic nature of the resulting interaction betwwn the externally induced air stream and the buoyancy-driven flow generated by the source is investigated. As a result, the best location of the heat source to make the active heat transfer is 0.075[m] from the left wall on the floor. The trends observed are also discussed in terms of heat removal from practical systems such as electric circuitry.

• Tribology and Lubricants
• /
• v.32 no.1
• /
• pp.1-8
• /
• 2016

This paper focuses on the thermal deformation induced preload change in the tilting pad journal bearing, using a three-dimensional (3D) thermo-hydro-dynamic (THD) approach. Preload is considered as a critical factor in designing the tilting pad journal bearing. The initial preload measured under nil external load and nil thermal gradient is influenced by two factors, namely, the thermal deformation and elastic deformation. Thermal deformation is due to a temperature distribution in the bearing pads, whereas the elastic deformation is due to fluid forces acting on the pads. This study focuses on the changes induced in preload and film clearance due to thermal deformation. The generalized Reynolds equation is used to evaluate the force of the fluid and the 3D energy equation is used to calculate the temperature of the lubricant. The abovementioned equations are combined by establishing a relationship between viscosity and temperature. The heat transfer within the bearing pads, the lubricant, and the spinning journal is calculated using the heat flux boundary condition. The 3D Finite Element Method (FEM) is used in modeling the (1) heat conduction in the spinning journal and bearing pads, (2) thermal gradient induced thermal distortion of the spinning journal and pads, and (3) viscous shearing, and heat conduction and convection in a thin film. This evaluation method has an increased fidelity, and it can prove to be a cost-effective tool that can be used by designers to predict the dynamic behavior of a bearing.

• Ocean and Polar Research
• /
• v.30 no.3
• /
• pp.261-275
• /
• 2008

The profile of a fixed site at station M ($34.77^{\circ}N,\;129.13^{\circ}E$) in the Korea Strait was studied from March 2006 to February 2007. The aim was to understand the relationship between the annual thermal stratification pattern and seasonal variation in phytoplankton community structure. Physicochemical factors including temperature, salinity and nutrient concentrations, which strongly influence the proliferation and diversity of phytoplankton, were measured. The study period was divided into three due to the characteristic of thermohaline structures; mixed I (March-May 2006), stratified (June-November 2006) and mixed II(December 2006-Feburuary 2007). Diatoms dominated during the mixed I (89%) and II (48%) periods, while nanoplankton group occupied over 83% of total population during the stratified period. The dominant species during the mixed I and II was Chaetoceros socialis (47% and 29%, respectively), while during the stratified period Gyrodinium sp.(4%) was the most dominant. Averaged total chl a concentrations during the mixed I and II periods were 0.61 mg $m^{-3}$ and 0.72 mg $m^{-3}$, respectively, which were at least two-fold higher than that during the stratified period (0.30 mg $m^{-3}$). The vertical mixing and convection process of the water column induced nutrient supply from the bottom layer to the euphotic zone. It also led to the dominance of diatoms during the mixed periods, whereas small phytoplankton prevailed over large phytoplankton as stratification blocked the upward movement of nutrients to subsurface during the stratified period. During the mixed I and II periods, microplanktonic chl a dominated concentrations (50% and 48%, respectively), while picoplanktonic chl a occupied over 37% of total chl a during the stratified period.

• Journal of Welding and Joining
• /
• v.17 no.1
• /
• pp.104-115
• /
• 1999

There have been many studies on the two dimensional thermo-elasto-plastic analysis in welding process, mostly from viewpoint of residual stresses. In this study, the temperature distribution, transient thermal stress, and angular distortion during bead-on-plate gas metal arc welding of rectangular plates were analyzed by using the finite element method. A nonlinear heat transfer analysis was first performed by taking account of the temperature-dependent material properties and convection heat losses on the surface. This was followed by a thermo-elasto-plastic stresses and distortion analysis that incorporates the constrained boundary condition of the two dimensional solution domain to get the three dimensional size effect of the plate. The constrained boundary conditions adopted in this study were the constant displacement condition over the whole two dimensional section for axial movement in the welding direction, and the force boundary condition for rotational movementof the domain around the axis of the welding direction. It could be revealed that the theoretical predictions of the angular distortion have an improved agreement with the experimentally obtained data presented in the previous study.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.6 s.144
• /
• pp.619-630
• /
• 2005

A two-frame PIV (Particle Image Velocimetry) technique is used to investigate the wake characteristics behind a marine propeller with 4 blades at high Reynolds number. For each of 9 different blade phases from $ 0^{\circ} $ to $ 80^{\circ} $, one hundred and fifty instantaneous velocity fields are measured. They are ensemble averaged to study the spatial evolution of the propeller wake in the region ranging from the trailing edge to one propeller diameter (D) downstream location. The phase-averaged mean velocity shows that the trailing vorticity is related to radial velocity jump, and the viscous wake is affected by boundary layers developed on the blade surfaces and centrifugal force. Both Galilean decomposition method and vortex identification method using swirling strength calculation are very useful for the study of vortex behaviors En the propeller wake legion. The slipstream contraction occurs in the near-wake region up to about X/D : 0.53 downstream. Thereafter, unstable oscillation occurs because of the reduction of interaction between the tip vortex and the wake sheet behind the maximum contraction point.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• Proceedings of the KIEE Conference
• /
• 2007.07a
• /
• pp.1185-1186
• /
• 2007

The UTM-01 developed in 1998 was the first maglev vehicle in Korea for the urban transit maglev (UTM) system. Through the improvement of UTM-01 and development of UTM02, the commercialization of the UTM system is being prepared now. In order to prepare for the commercialization of maglev, it is necessary that an optimal design of the levitation magnet should be provided for the safe operation of the vehicle. The levitation force is formed through the function of magnetic flux density on the top of magnet poles and gap between magnet pole and guide rail. To generate a magnetic field that is high enough to levitate the vehicle, ferromagnetic materials, such as pure iron for magnet pole and SS400 for guide rail, were used. The heat generated by $I^2R$ loss of magnet conductor makes the thermal convection on the surface of magnet including coil and poles. As these two characteristics are nonlinear phenomena, this paper deals with the nonlinear analysis on the magnetic and thermal properties of the U-type levitation magnet by using 3-D finite element method (FEM). Base on the analysis results, a small scale U-type magnet was designed, manufactured, and tested and it was verified that the magnet manufactured was satisfactory to all the design specifications.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.2
• /
• pp.123-136
• /
• 1996

A three-dimensional Galerkin-FEM model which can handle the temporal and spatial variation of density is presented. The hydrostatic approximation is used and density effects are included by means of conservation equation of heat and the equation of state. The finite difference grids are used in the horizontal plane and a set of linear-shape functions is used for the vertical expansion. The similarity transform is introduced to solve resultant matrix equations. The proposed model was first applied to the density-driven circulation in an idealized basin in the presence of the heat exchange between the air and the sea. The advection terms in the momentum equation were ignored, while the convection terms were retained in the heat equation. Coefficients of the vertical eddy viscosity and diffusivity were fixed to be constant. Calculation in a non-rotating idealized basin shows that the difference in heat capacity with depth gives rise to the horizontal gradient of temperature. Consequently, there is a steady new in the upper layer in the direction of increasing depth with compensatory counter flow .in the lower layer. With Coriolis force, geostrophic flow was predominant due to the balance between the pressure gradient and the Coriolis force. As a test in region of irregular topography, the model is applied to the Yellow Sea. Although the resultant flow was very complex, the character of the flow Showed to be geostrophic on the whole.

• Applied Chemistry for Engineering
• /
• v.2 no.2
• /
• pp.97-107
• /
• 1991

Highly-siliceous dealuminated zeolite, silicalite(end member of ZSM-5) was synthesized from a batch composition of 2.55 $Na_2O-5.0$ TPABr-$100SiO_2-2800H_2O $ at $180^{\circ}C$ and at times ranging from one to seven days of reaction time. Autoclaves containing the synthesis mixture were centrifuged within the specially-equipped convection oven to provide an elevated gravitational force field like 30 and 50 G. Tests were also conducted at normal gravity. For synthesis performed under elevated gravities, average and maximum crystal sizes were substantially greater than those synthesized under normal gravity and product yields were also found to be affected by elevated gravity ; that is, product yields were substantially enhanced under elevated gravity from 4 % to 55 % with respect to normal gravity. The average crystal sizes of silicalite synthesized at normal gravity were 50 to $70{\mu}m$ over an entire range of reaction time, one to seven days while the average crystal sizes synthsized under elevated gravities, 30 and 50 G, were 160 to $190{\mu}m$ respectively. For the elevated gravity, in particular, two separate nucleations and growths were observed. For examples, at 50G, large crystals of $200{\mu}m$ were produced through the second growing stage after 5 days of reaction following the rapid first growing stage where fairly large crystals of $135{\mu}m$ were produced only in 2 days of reaction. The maximum crystal sizes obtained through the above two growing stages were 190 and $300{\mu}m$, respectively. A discussion of how elevated gravity affects nucleation, growth, yield and crystal size of silicalite is presented.