• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.10
• /
• pp.2418-2425
• /
• 1993

A two-dimensional transient heat transfer model for reactive gas assisted laser cutting process with a moving Gaussian heat source is developed using a numerical finite difference technique. The kerf width, melting front shape and temperature distribution were calculated by using the boundary-fitted coordinate system to handle the ejection of workpiece material and heat input from reaction and evaporation. An analytical solution for cutting front movement was adopted and numerical simulation was performed to calculate the temperature distribution and melting front thickness. To calculate the moving velocity of cutting front, the normal distribution of the cutting gas velocity was used. The kerf width was revealed to be dependent on the cutting velocity, laser power and cutting gas velocity.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.5
• /
• pp.1338-1346
• /
• 1994

The objective of this study is to investigate the effect of various parameters, such as temperature, mean current density and voltage on the performance of phosphoric acid fuel cell (PAFC) by numerical analysis. Two types of flow passages, which are Z-parallel type and Z-counter type, are evaluated to obtain the best current density and temperature distribution. Parametric studies and sensitivity analysis of the PAFC system's operation in single cell are accomplished. A steady state simulation of the entire system is developed using nonlinear ordinary differential equations. The finite difference method and trial and error procedures are used to obtain a solution.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.12 no.12
• /
• pp.2253-2258
• /
• 2008

To design micro flow sensor varying depending on temperature of driving heater in the detector of Oxide semiconductor, Markov chain MCM(MCMCM), which is a kind of stochastic and microscopic method, was introduced. The formulation for the thermal transfer equation based on the FDM to obtain the MCMCM solution was performed and investigated, in steady state case. MCMCM simulation was successfully applied, so that its application can be expanded to a three-dimensional model with inhomogeneous material and complicated boundary.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.1149-1154
• /
• 2005

The finite difference lattice Boltzmann method(FDLBM) is a quite recent approach for simulating fluid flow, which has been proven as a valid and efficient tool in a variety of complex flow problems. It is considered an attractive alternative to conventional FDM and FVM, because it recovers the Navier-Stokes equations and is computationally more stable, and easily parallelizable to simulate for various laminar flows and a direct simulation of aerodynamics sounds. However, the research of a numerical simulation of turbulent flow by FDLBM, which is important to analyze the structure of turbulent flow in engineering fields, is not carried out. In this research, the FDLBM built in the turbulent model is applied, and a flowfield around 2-dimensional square to validate the applied model with 2D9V is simulated. Besides, 2D computation of the cavity noise generated by flow over a cavity at a Mach number of 0.1 and a Reynolds number based on cavity depth of 5000 is calculated. The computation result is well presented a understanding of the physical phenomenon of tonal noise occurred primarily by well-jet shear layer and vortex shedding and an aeroacoustic feedback loop.

• Computers and Concrete
• /
• v.16 no.5
• /
• pp.669-682
• /
• 2015

Degradation of reinforced concrete (RC) structures due to chloride penetration followed by reinforcement corrosion has been a serious problem in civil engineering for many years. The numerical simulation methods at present are mainly finite element method (FEM) and finite difference method (FDM), which are based on mesh. Mesh generation in engineering takes a long time. In the present article, the numerical solution of chloride transport in concrete is analyzed using radial point interpolation method (RPIM) and element-free Galerkin (EFG). They are all meshless methods. RPIM utilizes radial polynomial basis, whereas EFG uses the moving least-square approximation. A Galerkin weak form on global is used to attain the discrete equation, and four different numerical examples are presented. MQ function and appropriate parameters have been proposed in RPIM. Numerical simulation results are compared with those obtained from the finite element method (FEM) and analytical solutions. Two case of chloride transport in full saturated and unsaturated concrete are analyzed to test the practical applicability and performance of the RPIM and EFG. A good agreement is obtained among RPIM, EFG, and the experimental data. It indicates that RPIM and EFG are reliable meshless methods for prediction of chloride concentration in concrete structures.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2009.09a
• /
• pp.1441-1446
• /
• 2009

Numerical analysis is a powerful method in evaluating the soil-pile-structure interaction under the dynamic loading, and this approach has been applied to the practical area due to the development of computer technology. Finite Difference Method, one of the most popular numerical methods, is sensitive to the shape and the number of mesh. However, the trial and error approach is conducted to obtain the accurate results and the reasonable simulation time because of the lack of researches about mesh size and the number. In this study, FLAC 3D v3.1 program(FDM) is used to simulate the dynamic pile model tests, and the numerical results are compared with the 1G shaking table tests results. With the different size and shape of mesh, the responses of pile behavior and the simulation time are estimated, and the optimum mesh sizes in dynamic analysis of single pile is studied.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.24 no.4
• /
• pp.309-313
• /
• 2011

In this study, we proposed a novel electrode structure for the fringe field switching (FFS) mode LCD and performed a three-dimensional computer simulation to calculate the optical transmittance for the new structure. In the simulation Erickson-leslie equation and Berreman $4{\times}4$ matrix were used for obtaining the director distribution profiles of liquid crystal molecules and the electro-optical characteristics, respectively. Considering the complexity of the motional equation of the liquid crystal molecules, FDM (finite difference method) was used as a numerical method. From the results, We revealed that the light transmission of the newly designed pixel structure is expended to the edge of the pixel electrode. We also confirmed that the light transmittance increased more than 13% compared to that of the conventional electrode structure.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.26 no.1
• /
• pp.48-61
• /
• 1989

In this paper, two-dimensional numerical simulation of GaAs MESFFT with 0.7${\mu}m$ gate length is perfomed. Drift-diffusion model which consider that mobility is a function of local electric field, is used. As a discretization method, instead of FDM (finite difference method) and FEM (finite element method), the Control-Volume Formulation (CVF) is used and as a numerical scheme current hybrid scheme or upwind scheme is replaced by power-law scheme which is very approximate to exponential scheme. In the process of numerical analysis, Peclet number which represents the velocity ratio of drift and diffusion, is introduced. And using this concept a current equation which consider numerical scheme at the interface of control volume, is proposed. The I-V characteristics using the model and numerical method has a good agreement with that of previous paper by others. Therefore, it is confined that it may be useful as a simulator for GaAs MESFET. Besides I-V characteristics, the mechanism of both velocity saturation in drift-diffusion model is described from the view of velocity and electric field distribution at the bottom of the channel. In addition, the relationship between the mechanism and position of dipole and drain current, are described.

• Journal of the Institute of Electronics Engineers of Korea CI
• /
• v.48 no.1
• /
• pp.90-100
• /
• 2011

We present a parallel bi-conjugate gradient (Bi-CG) matrix solver for large scale Bio-FET simulations based on recent graphics processing units (GPUs) which can realize a large-scale parallel processing with very low cost. The proposed method is focused on solving the Poisson equation in a parallel way, which requires massive computational resources in not only semiconductor simulation, but also other various fields including computational fluid dynamics and heat transfer simulations. As a result, our solver is around 30 times faster than those with traditional methods based on single core CPU systems in solving the Possion equation in a 3D FDM (Finite Difference Method) scheme. The proposed method is implemented and tested based on NVIDIA's CUDA (Compute Unified Device Architecture) environment which enables general purpose parallel processing in GPUs. Unlike other similar GPU-based approaches which apply usually 32-bit single-precision floating point arithmetics, we use 64-bit double-precision operations for better convergence. Applications on the CUDA platform are rather easy to implement but very hard to get optimized performances. In this regard, we also discuss the optimization strategy of the proposed method.

• Smart Structures and Systems
• /
• v.21 no.5
• /
• pp.623-635
• /
• 2018

In order to study the mechanical behavior of shield tunnel segments during assembly stage, the in-situ tests and FDM numerical simulation were conducted based on the Foguan Shiziyang Tunnel with large cross-section. Analysis for the load state of the assembling segments in different assembly steps as well as the investigation for the changing of inner forces and longitudinal stress of segments with assembling steps were carried out in this paper. By comparing the tested results with the simulated results, the conclusions and suggestions could be drawn as follows: (1) It is the most significant for the effects on axial force and bending moment caused by the assembly of adjacent segment, followed by the insertion of key segment while the effects in the other assembly steps are relative smaller. With the increasing value of axial force, the negative bending moment turns into positive and remains increasing in most monitored sections, while the bending moment of segment B1and B6 are negative and keeping increasing; (2) The closer the monitored section to the adjacent segments or the key segment, the more significant the internal forces response, and the monitored effects of key segment insertion are more obvious than that of calculation; (3) The axial forces are all in compression during assembling and the monitored values are about 1.5~1.75 times larger than the calculated values, and the monitored values of bending moment are about 2 times the numerical calculation. The bending moment is more sensitive to the segments assembly process compared with axial force, and it will result in the large bending moment of segments during assembling when the construction parameters are not suitable or the assembly error is too large. However, the internal forces in assembly stage are less than those in normal service stage; (4) The distribution of longitudinal stress has strong influence on the changing of the internal forces. The segment side surface and intrados in the middle of two adjacent jacks are the crack-sensitive positions in the early assembly stage, and subsequently segment corners far away from the jacks become the crack-sensitive parts either.