• 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 the Computational Structural Engineering Institute of Korea
• /
• v.30 no.6
• /
• pp.495-500
• /
• 2017

In the boundary layer of supersonic or hypersonic vehicles, there is the conversion from kinetic energy to thermal energy, called aerodynamic heating. Aerodynamic heating has to be considered to design supersonic vehicles, because it induces severe heat flux to surface. Transient heat transfer analysis with CFD is used to predict thermal response of vehicles, however transient heat transfer analysis needs excessive computing powers. Loosely coupled method is widely used for evaluating thermal response, however it needs to be revised for overestimated heat flux. In this research, quasi-transient method, which is combined loosely coupled method and conjugate heat transfer analysis, is proposed for evaluating thermal response with efficiency and reliability. Defining reference time of splitting flight scenario for transient simulation is important on accuracy of quasi-transient method, however there is no algorithm to determine. Therefore the research suggests the algorithm with various flow conditions to define reference time. Supersonic flow field of blunt body with constant acceleration is calculated to evaluate quasi-transient method. Temperature difference between transient and quasi-transient method is about 11.4%, and calculation time reduces 28 times for using quasi-transient method.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.62 no.6
• /
• pp.131-142
• /
• 2020

The objective of this study was to evaluate future inundation risk of farmland according to the application of coupled model intercomparison project phase 5 (CMIP5) and coupled model intercomparison project phase 6 (CMIP6). In this study, future weather data based on CMIP5 and CMIP6 general circulation model (GCM) were collected, and inundation was simulated using the river modeling system for small agricultural watershed (RMS) and GATE2018 in the Tanjung district of the Moohan stream watershed. Although the average probable rainfall of CMIP5 and CMIP6 did not show significant differences as a result of calculating the probability rainfall, the difference between the minimum and maximum values was significantly larger in CMIP6. The results of the flood discharge calculation and the inundation risk assessment showed similar to trends to those of probability rainfall calculations. The risk of inundation in the future period was found to increase in all sub-watersheds, and the risk of inundation has been analyzed to increase significantly, especially if CMIP6 data are used. Therefore, it is necessary to consider climate change effects by utilizing CMIP6-based future weather data when designing and reinforcing water structures in agricultural areas in the future. The results of this study are expected to be used as basic data for utilizing CMIP6-based future weather data.

• Tunnel and Underground Space
• /
• v.26 no.2
• /
• pp.59-67
• /
• 2016

In this paper, we perform a numerical analysis to evaluate the potential of fault reactivation during gas production from hydrate bearing sediments and the moment magnitude of induced seismicity. For the numerical analysis, sequential coupling of TOUGH+Hydrate and FLAC3D was used and the change in effective stress and consequent geomechanical deformation including fault reactivation was simulated by assuming that Mohr-Coulomb shear resistance criterion is valid. From the test production simulation of 30 days, we showed that pore pressure reduction as well as effective stress change hardly induces the fault reactivation in the vicinity of a production well. We also investigated the influence of stress state conditions to a fault reactivation, and showed that normal fault stress regime, where vertical stress is relatively greater than horizontal, may have the largest potential for the reactivation. We tested one simulation that earthquake can be induced during gas production and calculated the moment magnitude of the seismicity. Our calculation presented that all the magnitudes from the calculation were negative values, which indicates that induced earthquakes can be grouped into micro-seismic and as small as hardly perceived by human beings. However, it should be noted that the current simulation was carried out using the highly simplified geometric model and assumptions such that the further simulations for a scheduled test production and commercial scale production considering complex geometric conditions may produce different results.

• Journal of Electrical Engineering and Technology
• /
• v.13 no.6
• /
• pp.2402-2411
• /
• 2018

Fluoronitriles-$CO_2$ gas mixtures are promising alternatives to $SF_6$ in environmentally-friendly gas-insulated transmission lines (GILs). Insulating gas heat transfer characteristics are of major significance for the current-carrying capacity design and operational state monitoring of GILs. In this paper, a three-dimensional calculation model was established for a GIL using the thermal-fluid coupled finite element method. The calculated results showed close agreement with experimentally measured data. The temperature distribution of a GIL filled with the Fluoronitriles-$CO_2$ mixture was obtained and compared with those of GILs filled with $CO_2$ and $SF_6$. Furthermore, the effects of the mixture ratio of the component gases and the gas pressure on the temperature rise and current-carrying capacity of the GIL were analyzed. Results indicated that the heat transfer performance of the Fluoronitriles-$CO_2$ gas mixture was better than that of $CO_2$ but worse than that of $SF_6$. When compared with $SF_6$, use of the Fluoronitriles-$CO_2$ gas mixture caused a reduction in the GIL's current-carrying capacity. In addition, increasing the Fluoronitriles gas component ratio or increasing the pressure of the insulating gas mixture could improve the heat dissipation and current-carrying capacity of the GIL. These research results can be used to design environmentally-friendly GILs containing Fluoronitriles-$CO_2$ gas mixtures.

• Nuclear Engineering and Technology
• /
• v.52 no.4
• /
• pp.721-733
• /
• 2020

In nuclear engineering, the occurrence of critical heat flux (CHF) is complicated for rod bundle, and it is much more difficult to predict the CHF when it is in natural circulation under motion condition. In this paper, the dryout-type CHF is investigated for the rod bundle in a natural circulation loop under rolling motion condition based on the coupled analysis of subchannel method, a one-dimensional system analysis method and a CHF mechanism model, namely the three-fluid model for annular flow. In order to consider the rolling effect of the natural circulation loop, the subchannel model is connected to the one-dimensional system code at the inlet and outlet of the rod bundle. The subchannel analysis provides the local thermal hydraulic parameters as input for the CHF mechanism model to calculate the occurrence of CHF. The rolling motion is modeled by additional motion forces in the momentum equation. First, the calculation methods of the natural circulation and CHF are validated by a published natural circulation experiment data and a CHF empirical correlation, respectively. Then, the CHF of the rod bundle in a natural circulation loop under both the stationary and rolling motion condition is predicted and analyzed. According to the calculation results, CHF under stationary condition is smaller than that under rolling motion condition. Besides, the CHF decreases with the increase of the rolling period and angular acceleration amplitude within the range of inlet subcooling and mass flux adopted in the current research. This paper can provide useful information for the prediction of CHF in natural circulation under motion condition, which is important for the nuclear reactor design improvement and safety analysis.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.22 no.6
• /
• pp.79-88
• /
• 2008

In this study electromagnetic fields near electric power lines were derived by dipole antenna theory and electromagnetic fields near 3 phase power lines with vertical configurations were formulated and could be computed easily using these formula. It seems that those formula could be applicable to the consideration of electromagnetic fields during the design of transmission and distribution lines. Those formulated equations on elements of electromagnetic fields were applied to the model of a transmission-line system and were calculated by Matlab programs. The calculation results are follows. For variation of horizontal distance profiles of $E_y$ and $B_z$ are same each other, and also those of $B_y$ and $E_z$ are same each other. This means that coupled elements of E and B are perpendicular each other and have the propagation direction of the right-hand system such as $x{\rightarrow}E_y{\rightarrow}B_z$. Resultant electric field E is dominated by the element $E_y$ and resultant magnetic field B is dominated by the element $B_z$.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.5
• /
• pp.575-582
• /
• 2010

In recent years, several rapid-mold-heating techniques that can be used for the injection molding of thin-walled parts or micro/nano structures have been developed. High-frequency induction heating, which involves heating by electromagnetic induction, is an efficient method for the rapid heating of mold surfaces. The present study proposes an integrated numerical model of the high-frequency induction heating process and the resulting injection molding process. To take into account the effects of thermal boundary conditions in induction heating, we carry out a fully integrated numerical analysis that combines electromagnetic field calculation, heat transfer analysis, and injection molding simulation. The proposed integrated simulation is extended to the injection molding of a thin-wall part, and the simulation results are compared with the experimental findings. The validity of the proposed simulation is discussed according to the ways of the boundary condition imposition.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.20 no.1
• /
• pp.29-35
• /
• 2010

This paper studies causes of the L-1 blade damage of a low pressure turbine, which was found during the scheduled maintenance, in 500 MW fossil power plants. Many failures of turbine blades are caused by the coupling of aerodynamic forcing with bladed-disk vibration characteristics. In this study the coupled vibration characteristics of the L-1 turbine bladed-disk in a fossil power plant is shown for the purpose of identifying the root cause of the damage and confirming equipment integrity. First, analytic and experimental modal analysis for the bladed-disk at zero rpm as well as a single blade were performed and analyzed in order to verify the finite element model, and then steady stresses, natural frequencies and corresponding mode shapes, dynamic stresses were calculated for the bladed-disk under operation. Centrifugal force and steady steam force were considered in calculation of steady and dynamic stress. The proximity of modes to sources of excitation was assessed by means of an interference diagram to examine resonances. In addition, fatigue analysis was done for the dangerous modes of operation by a local strain approach. It is expected that these dynamic characteristics will be used effectively to identify the root causes of blade failures and to perform prompt maintenance.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.291-305
• /
• 2010

A three-dimensional time-domain calculation method is of crucial importance in prediction of the motions and wave loads of a ship advancing in a severe irregular sea. The exact solution of the free surface wave-ship interaction problem is very complicated because of the essentially nonlinear boundary conditions. In this paper, an approximate body nonlinear approach based on the three-dimensional time-domain forward-speed free-surface Green function has been presented. The Froude-Krylov force and the hydrostatic restoring force are calculated over the instantaneous wetted surface of the ship while the forces due to the radiation and scattering potentials over the mean wetted surface. The time-domain radiation and scattering potentials have been obtained from a time invariant kernel of integral equations for the potentials which are discretized according to the second-order boundary element method (Hong and Hong 2008). The diffraction impulse-response functions of the Wigley seakeeping model advancing in transient head waves at various Froude numbers have been presented. A simulation of coupled heave-pitch motion of a long rectangular barge advancing in regular head waves of large amplitude has been carried out. Comparisons between the linear and the approximate body nonlinear numerical results of motions and wave loads of the barge at a nonzero Froude number have been made.