• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.6
• /
• pp.399-409
• /
• 2017

This study numerically investigates the characteristics of the velocity field including the average flow velocity, longshore current and turbulent kinetic energy acting as the main external forces of the salient formed behind the permeable submerged breakwaters. Shoreline response is also predicted by the longshore-induced flux. In this paper, a three-dimensional numerical wave tank based on the OLAFOAM, CFD open source code, is utilized to simulate the velocity field around permeable submerged breakwaters under the formation condition of salient. The characteristics of the velocity field around permeable submerged breakwaters with respect to the gap width between breakwaters and the installing position away from the shoreline under a range of regular waves for different wave height are evaluated. The numerical results revealed that as the gap width between breakwaters increases, the longshore currents become stronger. Furthermore, as the gap width becomes narrower, the point where flow converges moves from the center of the breakwater to the head part. As a result, it is possible to understand the formation of the salient formed behind the submerged breakwaters. In addition, it was found that the longshore currents caused by the gap width between breakwaters and the installation position away from the shoreline are closely related to the turbulent kinetic energy.

• Journal of Ocean Engineering and Technology
• /
• v.30 no.2
• /
• pp.91-99
• /
• 2016

The global demand for oil and natural gas has increased, and resource development is moving to the deep sea. Floating and flexible offshore structures such as semi-submersible, spar, and FPSO structures have been widely used. The major equipment of floating structures is always exposed to waves, currents, and other marine environmental factors, which cause structural damage. Moreover, flexible risers are susceptible to an exciting force due to the motion of the floating body. The inline and transverse responses from the three-dimensional behavior of a floating structure occur because of various forces. Typical risers are made of steel pipe and applied in the oil and gas development field, but flexible materials such as polyethylene are suitable for OTEC risers. Consequently, the optimal design of a flexible offshore plant requires a dynamic behavior analysis of slender bodies made of the different materials commonly used for offshore flexible risers. In this study, a three-dimensional motion measurement device was used to analyze the displacements of riser models induced by external force factors, and forced oscillation of a riser was linked to forced oscillation under a steady flow and regular wave condition.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.3 no.2
• /
• pp.18-24
• /
• 2000

A two-dimensional numerical method for inviscid two-fluid flows with a significant entrainment into both directions is established, and the oil leakage from a non-pressurized underwater pipe is studied. The interface between two fluids is modeled at a vortex sheet. The flow field and the subsequent interface evolution are solved by using the vortex-in-cell method. For longer flow simulation with a realistic two fluids interaction, an efficient merging scheme is introduced. In the Boussinesq limit, the speed of the external fluid intrusion into the pipe is very close to the existing mathematical models, and the lock exchange is observed in spite of a significant roll-up of the interface and entrainments. It is believed that the developed method can be utilized effectively for further detailed studies on various two-fluid flows which are encountered in many different marine oil spill problems.

• Wind and Structures
• /
• v.15 no.4
• /
• pp.331-343
• /
• 2012

Thunderstorm downbursts are responsible for numerous structural failures around the world. The wind characteristics in thunderstorm downbursts containing vortex rings differ with those in 'traditional' boundary layer winds (BLW). This paper initially performs an unsteady-state simulation of the flow structure in a downburst (modelled as a impinging jet with its diameter being $D_{jet}$) using a computational fluid dynamics (CFD) method, and then analyses the pressure distribution on a solar updraft tower (SUT) in the downburst. The pressure field shows agreement with other previous studies. An additional pair of low-pressure region and high-pressure region is observed due to a second vortex ring, besides a foregoing pair caused by a primary vortex ring. The evolutions of pressure coefficients at five orientations of two representative heights of the SUT in the downburst with time are investigated. Results show that pressure distribution changes over a wide range when the vortices are close to the SUT. Furthermore, the fluctuations of external static pressure distribution for the SUT case 1 (i.e., radial distance from a location to jet center x=$D_{jet}$) with height are more intense due to the down striking of the vortex flow compared to those for the SUT case 2 (x=$2D_{jet}$). The static wind loads at heights z/H higher than 0.3 will be negligible when the vortex ring is far away from the SUT. The inverted wind load cases will occur when vortex is passing through the SUT except on the side faces. This can induce complex dynamic response of the SUT.

• The Journal of the Acoustical Society of Korea
• /
• v.37 no.5
• /
• pp.263-270
• /
• 2018

In most industry fields concerning external flow noise problems, the hybrid computational aeroacoustic techniques based on the FW-H (Ffowcs Williams and Hawkings) equation are widely used for its numerical efficiency. However, when the surface integral form of FW-H equation is used without volume quadrupole sources, it is known to generate significant non-physical noise in a certain case. Especially, in the case of a flow in which the tip vortex cavitation is formed in the distant downstream direction such as flow driven by an underwater propeller, the accuracy in noise prediction becomes poor unless it is not properly modelled. Therefore, in this study, the nonphysical acoustic waves caused by the surface integral form of FW-H equation is reduced by adding the quadrupole correction term. First, to verify the accuracy of the in-house code of FW-H equation, the noise by an axial fan used in the outdoor unit of air conditioner was calculated and compared with the results of ANSYS Fluent. In order to verify the effects of the quadrupole correction term, the noise prediction for isentropic vortex convection is performed and it is confirmed that the error is reduced by the quadrupole correction term. Finally, the noise prediction is performed for the flow field generated by the Clark-Y hydrofoil in underwater. It is confirmed that the error caused by the cavitation passing through the integral surface can be reduced by the quadrupole correction term.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.43 no.1
• /
• pp.139-147
• /
• 2023

As heavy precipitation rates have increased due to climate change, the risk of landslides has also become greater. Studies in the field of disaster risk assessment predominantly focus on evaluating intrinsic importance represented by the use or role of facilities. This work, however, focused on evaluating risks according to the external conditions of facilities, which were presented via debris flow simulation. A random walk model (RWM) was partially improved and used for the debris flow simulation. The existing RWM algorithm contained the problem of the simulation results being overly concentrated on the maximum slope line. To improve the model, the center cell height was adjusted and the inertia application method was modified. Facility information was collected from a digital topographic map layer. The risk level of each object was evaluated by combining the simulation result and the digital topographic map layer. A risk assessment technique suitable for the polygon and polyline layers was applied, respectively. Finally, by combining the evaluated risk with the attribute table of the layer, a system was prepared that could create a list of objects expected to be damaged, derive various statistics, and express the risk of each facility on a map. In short, we used an easy-to-understand simulation algorithm and proposed a technique to express detailed risk information on a map. This work will aid in the user-friendly development of a debris flow risk assessment system.

• The Journal of the Korea Contents Association
• /
• v.14 no.1
• /
• pp.435-443
• /
• 2014

The aim of this study was to look for the cause of shortage of K-League spectators from a consumer's perspective and to offer solutions. For this, we offered enhancement method of the K-League spectators through dual path model based on existing relative researches. The first path involved consumer reaction to external factors of the K-League game. This path explained relational view(=we-ness) of consumers on K-League game(Hypothesis 1, 2). The second path explained personal view on K-League game. This path means consumer reaction to K-League game itself(Hypothesis 3, 4). The empirical study was based on a field survey and structural equation model. The results showed that we-ness positively affected team loyalty(Hypothesis 1) and the team loyalty positively affected spectating intentions(Hypothesis 2). It also revealed that flow of K-League game positively affected consumer satisfaction(Hypothesis 3). Finally consumer satisfaction was found to affect spectating intentions(Hypothesis 4). The significance of present study is to extend scope of research which is related with cheerleading of sports game. And this study emphasize the importance of we-ness and flow in sports game.

• Progress in Superconductivity and Cryogenics
• /
• v.15 no.2
• /
• pp.34-38
• /
• 2013

AC Losses for face to face stacks of four identical coated conductors (CCs) were numerically calculated using the H-formulation combined with the E-J power law and the Kim model. The motive sample was the face to face stack of four 2 mm-wide CC tapes with 2 ${\mu}m$ thick superconducting layer of which the critical current density, $J_c$, was $2.16{\times}10^6A/cm^2$ on IBAD-MgO template, which was suggested for the mitigation of ac loss as a round shaped wire by Korea Electrotechnology Research Institute. For the calculation the cross section of the stack was simply modeled as vertically aligned 4 rectangles of superconducting (SC) layers with $E=E_o(J(x,y,t)/J_c(B))^n$ in x-y plane where $E_o$ was $10^{-6}$ V/cm, $J_c$(B) was the field dependence of current density and n was 21. The field dependence of the critical current of the sample measured in four-probe method was employed for $J_c$(B) in the equation. The model was implemented in the finite element method program by commercial software. The ac loss properties for the stacks were compared with those of single 4 cm-wide SC layers with the same critical current density or the same critical current. The constraint for the simulation was imposed in two different ways that the total current of the stack obtained by integrating J(x,y,t) over the cross sections was the same as that of the applied transport current: one is that one fourth of the external current was enforced to flow through each SC. In this case, the ac loss values for the stacks were lower than those of single wide SC layer. This mitigation of the loss is attributed to the reduction of the normal component of the magnetic field near the SC layers due to the strong expulsion of the magnetic field by the enforced transport current. On the contrary, for the other case of no such enforcement, the ac loss values were greater than those of single 4cm-wide SC layer and. In this case, the phase difference of the current flowing through the inner and the outer SC layers of the stack was observed as the transport current was increased, which was a cause of the abrupt increase of ac loss for higher transport current.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.15 no.1
• /
• pp.593-601
• /
• 2014

In this study, we investigated the methods of improving the capturing ability of acid fume by assessing the performance of slot-type external hood installed on both sides of an open surface tank for acid washing process. A field survey and the results of computational fluid dynamics revealed that capturing performance of existing hoods is very poor. To solve such problem, 'push-pull hood' that pushes from one side of an open surface tank and pulls on the other side was suggested. The initial prediction was that if a push-pull hood is used, the acid fume of an acid-washing tank surface could be moved towards the hood through the push flow. However, this study has confirmed that if the push flow velocity becomes too high, it could spread to other areas due to flooding from the hood. Therefore, if the push air supply is maintained at around 25 $m^3/min$(push 10 m/s), proper control flow is formed on the surface of a tank and acid fume that stayed at the upper part of the tank is smoothly captured toward the hood, significantly enhancing the capturing performance.