• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2008년도 춘계학술대회 논문집
• /
• pp.33-36
• /
• 2008

The transient behavior of a passive air breathing direct methanol fuel cell (DMFC) operated on vapor-feeding mode is studied in this paper. It generally takes 30 minutes after starting for the cell response to come to its steady-state and the response is sometimes unstable. A mathematical dynamic one-dimensional model for simulating transient response of the DMFC is presented. In this model a DMFC is decomposed into its subsystems using lumped model and divided into five layers, namely the anodic diffusion layer, the anodic catalyst layer, the proton exchange membrane (PEM), the cathodic catalyst layer and the cathodic diffusion layer. All layers are considered to have finite thickness, and within every one of them a set of differential-algebraic governing equations are given to represent multi-components mass balance, such as methanol, water, oxygen and carbon dioxide, charge balance, the electrochemical reaction and mass transport phenomena. A one-dimensional, isothermal and mass transport model is developed that captures the coupling between water generation and transport, oxygen consumption and natural convection. The single cell is supplied by pure methanol vapor from a methanol reservoir at the anode, and the oxygen is supplied via natural air-breathing at the cathode. The water is not supplied from external source because the cell uses the water created at the cathode using water back diffusion through nafion membrane. As a result of simulation strong effects of water transport were found out. The model analysis provides several conclusions. The performance drop after peak point is caused by insufficiency of water at the anode. The excess water at the cathode makes performance recovery impossible. The undesired crossover of the reactant methanol through the PEM causes overpotential at the cathode and limits the feeding methanol concentration.

• Journal of Mechanical Science and Technology
• /
• 제17권6호
• /
• pp.906-913
• /
• 2003

Lean premixed combustion has been considered as one of the promising solutions for the reduction of NOx emissions from gas turbines. However, unstable combustion of lean premixed flow becomes a real challenge on the way to design a reliable, highly efficient dry low NOx gas turbine combustor. Contrary to a conventional diffusion type combustion system, characteristics of premixed combustion significantly depend on a premixing degree of combusting flow. Combustion behavior in terms of stability has been studied in a model gas turbine combustor burning natural gas and air. Incompleteness of premixing is identified as significant perturbation source for inducing unstable combustion. Application of a simple convection time lag theory can only predict instability modes but cannot determine whether instability occurs or not. Low frequency perturbations are observed at the onset of instability and believed to initiate the coupling between heat release rate and pressure fluctuations.

• 한국화재소방학회:학술대회논문집
• /
• 한국화재소방학회 1996년도 학술발표회
• /
• pp.5-9
• /
• 1996

In foreign country such as U.S.A and Japan, considerable research has been done regarding the spread of smoke in room of fire involvement by using computer. But, in our country it has not been. So, this paper presents a detailed qualitative description of phenomena which occures during typical fire scenarios through numerical analysis. This research, in the view of field model, is focused on finding out the smoke movement and temperature distribution. And it is planned to analyze governing equation including smoke diffusion equation by numerical analysis with finite volume method and non-staggered grid system. The SIMPLE method for pressure-velocity couple and power-law scheme for convection terms are used. It shows that a plume is formed, hot plume is formed, hot plume gases impinge on the ceiling and they spread across it. then, it eventually reaches the bounding walls of the enclosure. It takes 60s for smoke to fill the enclosure.

• 대한수학회지
• /
• 제34권4호
• /
• pp.1049-1064
• /
• 1997

We study the existence and uniqueness of nonnegative singular solution u(x,t) of the semilinear parabolic equation $$ u_t = \Delta u - a \cdot \nabla(u^q) = u^p, $$ defined in the whole space $R^N$ for t > 0, with initial data $M\delta(x)$, a Dirac mass, with M > 0. The exponents p,q are larger than 1 and the direction vector a is assumed to be constant. We here show that a unique singular solution exists for every M > 0 if and only if 1 < q < (N + 1)/(N - 1) and 1 < p < 1 + $(2q^*)$/(N + 1), where $q^* = max{q, (N + 1)/N}$. This result agrees with the earlier one for N = 1. In the proof of this result, we also show that a unique singular solution of a diffusion-convection equation without absorption, $$ u_t = \Delta u - a \cdot \nabla(u^q), $$ exists if and only if 1 < q < (N + 1)/(N - 1).

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2000년도 추계 학술대회논문집
• /
• pp.129-134
• /
• 2000

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.

• 한국안전학회지
• /
• 제13권4호
• /
• pp.3-8
• /
• 1998

The vortex particle method, which includes viscous effects, consists of diffusion of boundary vorticity and creation of the vortex particles, convection, particle strength exchange, and particle redistribution. Accuracy of the boundary element method is very important since it creates the particles around the body at every time step. A boundary element method based on source panel was investigated as part of computation of unsteady separated flows by rising the vortex particle method. The potential flows were computed around a circular cylinder and a square cylinder. The results around the circular cylinder were compared with the exact solution, and the distribution of vorticity, in particular near the sharp comers of the square cylinder, is scrutinized for different number of panels.

• 설비공학논문집
• /
• 제6권1호
• /
• pp.20-28
• /
• 1994

A numerical study has been performed on enhancement of heat transfer in a forced convection of the modified driven-cavity flow which was previously found by the author to give a regular or chaotic stirring depending on the parameter value. It is found that for the present case wherein heat is transmitted between fluid and the surrounding walls, the chaotic stirring enhances the heat transfer at high Peclet numbers. The optimal condition of the flow modulation for the best heat transfer can be predicted by purely investigating the hydrodynamic facet, i.e. the stirring effect.

• 대한기계학회논문집B
• /
• 제26권11호
• /
• pp.1521-1530
• /
• 2002

This article proposes a pressure based method for predicting flows at all speeds. The compressible SIMPLE algorithm is extended to unstructured grid framework. Convection terms are discretized using second-order scheme with deferred correction approach. Diffusion term discretization is based on structured grid analogy that can be easily adopted to hybrid unstructured grid solver. This method also uses node centered scheme with edge based data structure for memory and computing time efficiency of arbitrary grid types. Both incompressible and compressible benchmark problems are solved using the above methodology. The demonstration of this method is extended to slip flow problem that has low Reynolds number but compressibility effect. It is shown that the proposed method can improve efficiency in memory usage and computing time without losing any accuracy.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.21-28
• /
• 1997

This paper presents a two-dimensional horizontal implicit model to general circulation in estuaries and coastal seas. The model is developed in non-orthogonal curvilinear coordinates system, using the Interpolating Matrix Method (IMM), in combination with a technique of multi-domain. In the propose model, the Saint-Venant equations are solved by a splitting-up technique, in the successive steps; convection, diffusion and wave propagation. The ability of the proposed model to deal with full scale nature is illustrated by the interpretation of a dye-tracing experiment in the Gironde estuary.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1997년도 춘계학술대회 논문집
• /
• pp.327-331
• /
• 1997

As demand for various kinds and small lot products has been increasing, batch annealing has been appreciated for its small restiction for the opteration. The cold spot of the coil is very important in the BAF(Batch annealing furnace) annealing process. Because of the annealing cycle time in the BAF, annealing was decided on the cold spot of the coil. So, we tested the effect,variation of cold spot, for hydrogen contents of atmospheric gas at the annealing furnace. As a result of several investigations. We confirmed the following characteristics ; after the heating and soaking,the cold spot of coil moved to 1/3 of coil thickness in the NHx atmospheric gas, but the mid point of the coil thickness is the cold spot in the Ax or .H/sub2. atmospheric gas. Therefore, the use of hydrogen instead of nitrogen as the protective gas,combined with high convection in batch annealing furnaces, has shown that considerable increases in furnace output and material quality are attainable. Owing to the low density, high diffusion and reducing character of hydrogen, a better transfer resulting in uniform material temperatures and improved coil surfaces can be achieved.