• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.18 no.2
• /
• pp.125-129
• /
• 2005

Two kinds of physical treatments were examined for the analysis both of intrinsic surface and interior nature of CuInS $e_2$[CIS] and CuGaS $e_2$[CGS] films grown in separated systems. For the first method, a selenium protection layer which was immediately deposited after the growth of the CIS was investigated. The Se cap layer protects CISe surface from oxidation and contamination during the transport under ambient atmosphere. The Se cap was removed by thermal annealing at temperature above 15$0^{\circ}C$. After the decapping treatment at 2$25^{\circ}C$ for 60 min, ultraviolet photoemission and inverse photoemission measurements of the CIS film showed that its valence band maximum(VBM) and conduction band minimum (CBM) are located at 0.58 eV below and 0.52 eV above the Fermi level $E_{F}$, respectively. For the second treatment, an Ar ion beam etching was exploited. The etching with ion kinetic energy $E_{k}$ above 500 eV resulted in broadening of photoemission spectra of core signals and occasional development of metallic feature around $E_{F}$. These degradations were successfully suppressed by decreasing $E_{k}$ below 400 eV. CGS films etched with the beam of $E_{k}$ = 400 eV showed a band gap of 1.7 eV where $E_{F}$ was almost centered.st centered.

• Journal of Hydrogen and New Energy
• /
• v.19 no.5
• /
• pp.423-429
• /
• 2008

This work was carried out to improve the performance of anodic electrocatalysts in direct ethanol fuel cell(DEFC). PtRu and $Pt_5Ru_4M$(M= Ni, Sn, Mo and W) electrocatalysts were prepared by using a $NaBH_4$ reduction method. Alloy crystal structure and particle size of electrocatalysts were characterized by X-ray diffraction(XRD) and transmission electron microscopy(TEM). The XRD analysis of the electrocatalysts revealed that the face-centered cubic(fcc) peaks shifted to slightly higher diffraction angles when third metals were added. Average size of the uniform particles was observed to be approximately $3{\sim}3.5\;nm$ from the TEM image. The electrochemical measurements were carried out in the solution 1M $H_2SO_4$ and 1M $C_2H_5OH$ at room temperature. Cyclic-voltammogram results showed that $Pt_5Ru_4W$ electrocatalyst exhibited much higher current density for ethanol oxidation of $2.73\;mA/cm^2$ than PtRu electrocatalyst of $0.73\;mA/cm^2$.

• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2005.06a
• /
• pp.231-238
• /
• 2005

On the basis of flow resistance theory the conceptual model and related mathematical descriptions is proposed for resistance modeling of groundwater flow in CPM(continuum Porous medium), DFN(discrete fracture network) and fractured-porous medium. The proposed model is developed on the basis of finite volume method assuming steady-state, constant density groundwater flow. The basic approach of the method is to evaluate inter-block flow resistance values for a staggered grid arrangement, i.e. fluxes are stored at cell walls and scalars at cell centers. The balance of forces, i.e. the Darcy law, is utilized for each control volume centered around the point where the velocity component is stored. The transmissivity (or permeability) at the interface is assumed to be the harmonic average of neighboring blocks. Flow resistance theory was utilized to relate the fluxes between the grid blocks with residual pressures. The flow within porous medium is described by three dimensional equations and that within an individual fracture is described by a two dimensional equivalent of the flow equations for a porous medium. Newly proposed models would contribute to develop flow simulation techniques with various matrix characteristics.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.7
• /
• pp.723-733
• /
• 2011

A numerical method of the CSF(Continuum Surface Force) model is presented for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.7
• /
• pp.721-730
• /
• 2012

The present paper presents a numerical study on multiphase flows induced by wall adhesion. The continuum surface force (CSF) model with the wall adhesion boundary condition model is used for calculating the surface tension force; this model is implemented in an in-house solution code (PowerCFD). The present method (code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with a volume capturing method (CICSAM) in a volume of fluid (VOF) scheme for phase interface capturing. The effects of wall adhesion are then numerically simulated by using the present method for a shallow pool of water located at the bottom of a cylindrical tank with no external forces such as gravity. Two different cases are computed, one in which the water wets the wall and one in which the water does not wet the wall. It is found that the present method efficiently simulates the surface tension-dominant multiphase flows induced by wall adhesion.

• Journal of Ship and Ocean Technology
• /
• v.10 no.4
• /
• pp.1-11
• /
• 2006

The finite volume based multi-block RANS code, WAVIS developed at MOERI is applied to the numerical self-propulsion test. WAVIS uses the cell-centered finite volume method for discretization of the governing equations. The realizable $k-{\epsilon}$ turbulence model with a wall function is employed for the turbulence closure. The free surface is captured with the two-phase level set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. The propeller forces are obtained using an unsteady lifting surface method based on potential flow theory. The numerical procedure followed the self-propulsion model experiment based on the 1978 ITTC performance prediction method. The self-propulsion point is obtained iteratively through balancing the propeller thrust, the ship hull resistance and towing force that is correction for Reynolds number difference between the model and full scale. The unsteady lifting surface code is also iterated until the propeller induced velocity is converged in order to obtain the propeller force. The self-propulsion characteristics such as thrust deduction, wake fraction, propeller efficiency, and hull efficiency are compared with the experimental data of the practical container ship. The present paper shows that hybrid RANS and potential flow based numerical method is promising to predict the self-propulsion parameters of practical ships as a useful tool for the hull form and propeller design.

• Journal of the Korean Electrochemical Society
• /
• v.11 no.2
• /
• pp.109-114
• /
• 2008

Though ethanol can theoretically generate 12 electrons during oxidation to carbon dioxide, the complete oxidation of ethanol is hard to achieve due to the strong bond between the two carbons in its molecular structure. Therefore, development of high activity catalyst for ethanol oxidation is necessary for the commercialization of direct ethanol fuel cell. In this study, some binary and ternary electrocatalysts of PtSn/C and PtSnAu/C have been synthesized and characterized. The catalysts were fabricated with modified polyol method with the amounts of 20 wt%, where the Pt : Sn ratios in the PtSn/C were 1 : 0, 4 : 1, 3 : 1, 2 : 1, 1.5 : 1, 1 : 1, 1 : 1.5 and Pt:Sn:Au ratios in the PtSnAu/C were 5 : 5 : 0, 5 : 4 : 1, 5 : 3 : 2, 5 : 2 : 3. From the XRD and TEM analysis results, the catalysts were found to have face centered cubic structure with particle size of around $1.9{\sim}2.4\;nm$. The activity in the ethanol oxidation was examined with cyclic voltammetry and the results indicated that PtSn(1.5 : 1)/C and PtSnAu(5 : 2 : 3)/C had the highest activity in each catalyst system. Further tests with single cell were performed with those catalysts. It was found that PtSn/C(1.5 : 1) exhibited the best performance while the long term stability of PtSnAu/C(5 : 2 : 3) is better than PtSn/C(1.5 : 1).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.4 s.247
• /
• pp.337-342
• /
• 2006

Natural convection flows in a cubical air-filled cavity that has one pair of opposing faces isothermal at different temperatures, $T_h\;and\;T_c$, respectively, the remaining four faces having a linear variation from $T_c\;to\;T_h$ are numerically simulated by a solution code(PowerCFD) using unstructured cell-centered method. Special attention is paid to three-dimensional flow and thermal characteristics according to the variation of inclination angle $\theta$ of the isothermal faces from horizontal: namely $\theta=0^{\circ},\;15^{\circ},\;30^{\circ},\;45^{\circ},\;50^{\circ},\;60^{\circ},\;75^{\circ}\;and\;90^{\circ}$. Comparisons of the average Nusselt number at the cold face are made with experimental benchmark results found in the literature. It is demonstrated that the average Nusselt number at the cold face has a maximum value around the inclination angle of $50^{\circ}$. It is also found that the code is capable of producing accurately the nature of the laminar convection in a cubical air-filled cavity with differentially heated walls.

• Nuclear Engineering and Technology
• /
• v.45 no.6
• /
• pp.707-720
• /
• 2013

We analyze piecewise homogenization with flux-weighted cross sections and preservation of averaged currents at the boundary of the homogenized domain. Introduction of a set of flux discontinuity ratios (FDR) that preserve reference interface currents leads to preservation of averaged region reaction rates and fluxes. We consider the class of numerical discretizations with one degree of freedom per volume and per surface and prove that when the homogenization and computing meshes are equal there is a unique solution for the FDRs which exactly preserve interface currents. For diffusion submeshing we introduce a Jacobian-Free Newton-Krylov method and for all cases considered obtain an 'exact' numerical solution (eight digits for the interface currents). The homogenization is completed by extending the familiar full assembly homogenization via flux discontinuity factors to the sides of regions laying on the boundary of the piecewise homogenized domain. Finally, for the familiar nodal discretization we numerically find that the FDRs obtained with no submesh (nearly at no cost) can be effectively used for whole-core diffusion calculations with submesh. This is not the case, however, for cell-centered finite differences.