• Journal of the Korean Society of Industry Convergence
• /
• v.4 no.3
• /
• pp.329-337
• /
• 2001

This simulator is developed for the analysis of a MOSFET based on Thermally Coupled Energy Transport Model(TCETM). The simulator has the ability to calculate not only stationary characteristics but also non - stationary characteristics of a MOSFET. It solves basic semiconductor devices equations including Possion equation, current continuity equations for electrons and holes, energy balance equation for electrons and heat flow equation, using finite difference method. The conventional semiconductor device simulation technique, based on the Drift-Diffusion Model (DDM), neglects the thermal and other energy-related properties of a miniaturized device. I, therefore, developed a simulator based on the Thermally Coupled Energy Transport Model (TCETM) which treats not only steady-state but also transient phenomena of such a small-size MOSFET. In particular, the present paper investigates the breakdown characteristics in transient conditions. As a result, we found that the breakdown voltage has been largely underestimated by the DDM in transient conditions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.103-121
• /
• 2015

In this paper, we briefly review and describe a projection algorithm for numerically computing the two-dimensional time-dependent incompressible Navier-Stokes equation. The projection method, which was originally introduced by Alexandre Chorin [A.J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comput., 22 (1968), pp. 745-762], is an effective numerical method for solving time-dependent incompressible fluid flow problems. The key advantage of the projection method is that we do not compute the momentum and the continuity equations at the same time, which is computationally difficult and costly. In the projection method, we compute an intermediate velocity vector field that is then projected onto divergence-free fields to recover the divergence-free velocity. Numerical solutions for flows inside a driven cavity are presented. We also provide the source code for the programs so that interested readers can modify the programs and adapt them for their own purposes.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.17-29
• /
• 2004

Internal waves are computed using a ghost fluid method on an unstructured grid. Discontinuities in density and dynamic pressure are captured in one cell without smearing or oscillations along a multimaterial interface. A time-accurate incompressible Navier-Stokes/Euler solver is developed based on a three-point backward difference formula for the physical time marching. Artificial compressibility is introduced with respect to pseudotime and an implicit method is used for the pseudotime iteration. To track evolution of an interface, a level set function is coupled with the governing equations. Roe's flux difference splitting method is used to calculate numerical fluxes of the coupled equations. To get higher order accuracy, dependent variables are reconstructed based on gradients which are calculated using Gauss theorem. For each edge crossing an interface, dynamic pressure is assigned for a ghost node to enforce the continuity of total pressure along the interface. Solitary internal waves are computed and the results are compared with other computational and experimental results.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.163-169
• /
• 1999

The numerical procedure has been developed for simulating incompressible viscous flow around a turbine stage with rotor-stator interaction. This study solves 2-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system. The Marker-and-Cell concept is applied to efficiently solve continuity equation. To impose an accurate boundary condition, O-H multiblocked grid system is generated. O-type grid and H-type grid is generated near and outer rotor-stator The cubic-spline interpolation is applied to handle a relative motion of a rotor to the stator. Turbulent flows have been modeled by the Baldwin- Lomax turbulent model. To validate present procedure, the time averaged pressure coefficients around the rotor and stator are compared with experiment and a good agreement obtained.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2005.04a
• /
• pp.137-140
• /
• 2005

Resin transfer molding (RTM) is one of the most popular processes for producing fiber reinforced polymer composites. In the manufacture of complex thick composite structures, analysis on flow front advancement on the resin impregnating the multi-layered fiber preform is helpful for the optimization of the process. In this study, three-dimensional mold filling simulation of RTM is carried out by using CVFEM (Control Volume Finite Element Method). On the assumption of isothermal flow of Newtonian fluid, Darcy’s law and continuity equation are used as governing equations. Different permeability tensors employed in each layer are obtained by experiments. Numerically predicted flow front is compared with experimental one in order to validate the numerical results. Flow simulations are conducted in the two mold geometries, rectangular plate and hollow cylinder. Permeability tensor of each layer preform in Cartesian coordinate system is transformed to cylinder coordinates system so that the flow within the multi-layered preforms of the hollow cylinder can be calculated exactly. Our emphasis is on the three dimensional flow analysis for circular three-dimensional braided preform, which shows outstanding mechanical properties such as high impact strength and toughness compared with other conventional two-dimensional laminar-structured preforms.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.188-192
• /
• 2004

A Numerical analysis method is developed in order to compute the flaw and wave field for 2-dimensional floating body in the free roll motion with 3 degrees of freedom. Navier-Stokes and continuity equations are gaverning equations in tire present study. Finite Difference method is introduced to discretize the governing equation. The free surface is traced by the interface tracking method and the grid system is fitted to boundaries including free surface and body surface, which is moving in the flow field. The numerical scheme is based on Maker and Cell method. For the sake of validation of the numerical method, the computed roll decay factors according to tire midship section shapes are compared with measured results. The numerical results are discussed in order to understand the effect of midship section shape on roll motion.

• Nuclear Engineering and Technology
• /
• v.49 no.1
• /
• pp.92-102
• /
• 2017

The analysis of steam generators as an interface between primary and secondary circuits in light water nuclear power plants is crucial in terms of safety and design issues. VVER-1000 nuclear power plants use horizontal steam generators which demand a detailed thermal hydraulics investigation in order to predict their behavior during normal and transient operational conditions. Two phase flow field simulation on adjacent tube bundles is important in obtaining logical numerical results. However, the complexity of the tube bundles, due to geometry and arrangement, makes it complicated. Employment of porous media is suggested to simplify numerical modeling. This study presents the use of porous media to simulate the tube bundles within a general-purpose computational fluid dynamics code. Solved governing equations are generalized phase continuity, momentum, and energy equations. Boundary conditions, as one of the main challenges in this numerical analysis, are optimized. The model has been verified and tuned by simple two-dimensional geometry. It is shown that the obtained vapor volume fraction near the cold and hot collectors predict the experimental results more accurately than in previous studies.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.4
• /
• pp.327-334
• /
• 1994

The manifold designates a pipe system discharging or absorbing water through multiple holes installed along the finite length of the pipe. The proper design of the manifold requires the pre-analysis of the hydraulic characteristics such as system head and flow rate. head loss and hole discharges, etc. On the contrary to the general pipe systems. the head loss along the pipe is hardly quantified in an explicit way since it is complicately varied by the size and arrangement of the holes. In the present study, both energy and continuity equations are employed to analyze the hydraulic characteristics, constituting nonlinear simultaneous equations which are solved by Newton method. In addition, a hydraulic experiment utilizing the manifold system equipped with an automatically controlled valve is performed to reproduce model tide. The result shows that the manifold system can be effectively used in a tidal basin where water flow should be maintained uniformly over the basin width.

• Journal of Ship and Ocean Technology
• /
• v.5 no.2
• /
• pp.50-61
• /
• 2001

In this paper, nonlinear interactions between water waves and a horizontally submerged circular cylinder are numerically simulated. In this case, the nonlinear interactions between them generated a wave breaking phenomenon. The wave breaking phenomenon plays an important role in the wave farce. Negative drifting forces are raised at shallow submerged cylinders under waves because of the wave breaking phenomenon. For the numerical simulation, a finite difference method based on the unsteady incompressible Navier-Stokes equations and the continuity equation is adopted in the rectangular grid system. The free surface is simulated with a computational simulation method of two-layer flow by using marker density. The results are compared with some existing computational and experimental results.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.16 no.3
• /
• pp.133-140
• /
• 2007

A 2-D model of fluid flow, mass transport and electrochemistry is analysed to examine the effect of current density at the current collector depending on active layer thickness of catlyst in polymer elecrolyte fuel cells. The finite element method is used to solve the continuity, potential and Maxwell-Stefan equations in the flow channel and gas diffusion electrode regions. For the material behavior of electrode reactions in the active catalyst layers, the agglomerate model is implemented to solve the diffusion-reaction problem. The calculated model results are described and compared with the different thickness of active catalyst layers. The significance of the results is discussed in the viewpoint of the current collecting capabilities as well as mass transportation phenomena, which is inferred that the mass transport of reactants dictates the efficiency of the electrode in the present analysis.