• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.34D no.8
• /
• pp.15-25
• /
• 1997

The integration method of carier concentrations to redcue the discretization error of th box integratio method used in the discretization of the three-dimensional poisson's equation is presented. The carrier concentration is approximated in the closed form as an exponential function of the linearly varying potential in the element. The presented method is implemented in the three-dimensional poisson's equation solver running under the windows 95. The accuracy and the convergence chaacteristics of the three-dimensional poisson's equation solver are compared with those of DAVINCI for the PN junction diode and the n-MOSFET under the thermal equilibrium and the DC reverse bias. The potential distributions of the simulatied devices from the three-dimensional poisson's equation solver, compared with those of DAVINCI, has a relative error within 2.8%. The average number of iterations needed to obtain the solution of the PN junction diode and the n-MOSFET using the presented method are 11.47 and 11.16 while the those of DAVINCI are 21.73 and 23.0 respectively.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.10
• /
• pp.947-954
• /
• 2009

Starting from the Eu's GH(Generalized Hydrodynamic) theory, the multi-species GH numerical solver is developed in this research and its computatyional behaviors are examined for the hypersonic rarefied flow over an axisymmetric body. To improve the accuracy of the developed multi-species GH solver, various slip-wall boundary conditions are tested and the computed results are compared. Additionally, in order to validate the entropy characteristics of the GH equation, the entropy production and entropy generation rates of the GH equation are investigated in the 1-dimensional normal shock structure test at a high Knudsen number.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1992.10a
• /
• pp.349-353
• /
• 1992

The purpose of this paper is to develop methods for the dynamic analysis of multibody system that consist of interconnected rigid and deformable component. The equations of motion are derived by using the Lagrange's equation and finite element theory for the elastic mechanism systems. The type of equation of motion is the differential algebraic equation included kinematic nonlinear algebraic equation. The generalized coordinate partitioning method is used for solving this equation. To show the validity of this analysis solver, couple of models were canalized and those results were compared with the commercial package(ADAMS).

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.36D no.11
• /
• pp.45-55
• /
• 1999

The small signal analyzer for the stationary drift-diffusion equation is developed. The slotboom variables of the potential, electron and hole concentrations for the response of applied small signal are defined and the stationary drift-diffusion equation is linearlized on DC operation point by $S^3A$ method. Frontal solver, which is used to solve the global matrix, progresses the accuracy of the solution in high frequency and minimizes the requirement of the memory. The simulations are executed on the structure of 3 dimensional N'P junction diode and 2 dimensional n-MOSFET to verify the proposed algorithm. The average relative errors of the conductance and the capacitance compared with MEDICI are about 26% and 0.67 for N'P junction diode and 7.75% and 2.24% for n-MOSFET. The simulation by the proposed algorithm can analyze the stationary drift-diffusion equation for applied small signal in high frequency region about 100GHz.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.31-35
• /
• 2006

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

• Journal of electromagnetic engineering and science
• /
• v.19 no.1
• /
• pp.6-12
• /
• 2019

An IE-FFT algorithm is implemented and applied to the electromagnetic (EM) solution of perfect electric conducting (PEC) scattering problems. The solution of the method of moments (MoM), based on the magnetic field integral equation (MFIE), is obtained for PEC objects with closed surfaces. The IE-FFT algorithm uses a uniform Cartesian grid to apply a global fast Fourier transform (FFT), which leads to significantly reduce memory requirement and speed up CPU with an iterative solver. The IE-FFT algorithm utilizes two discretizations, one for the unknown induced surface current on the planar triangular patches of 3D arbitrary geometries and the other on a uniform Cartesian grid for interpolating the free-space Green's function. The uniform interpolation of the Green's functions allows for a global FFT for far-field interaction terms, and the near-field interaction terms should be adequately corrected. A 3D block-Toeplitz structure for the Lagrangian interpolation of the Green's function is proposed. The MFIE formulation with the IE-FFT algorithm, without the help of a preconditioner, is converged in certain iterations with a generalized minimal residual (GMRES) method. The complexity of the IE-FFT is found to be approximately $O(N^{1.5})$and $O(N^{1.5}logN)$ for memory requirements and CPU time, respectively.

• Journal of applied mathematics & informatics
• /
• v.7 no.3
• /
• pp.905-915
• /
• 2000

In this article we discuss some aspects of operational Tau Method on delay differential equations and then we apply this method on the differential delay equation defined by $\omega(u)\;=\frac{1}{u}\;for\;1\lequ\leq2$ and $(u\omega(u))'\;=\omega(u-1)\;foru\geq2$, which was introduced by Buchstab. As Khajah et al.[1] applied the Recursive Tau Method on this problem, they had to apply that Method under the Mathematica software to get reasonable accuracy. We present very good results obtained just by applying the Operational Tau Method using a Fortran code. The results show that we can obtain as much accuracy as is allowed by the Fortran compiler and the machine-limitations. The easy applications and reported results concerning the Operational Tau are again confirming the numerical capabilities of this Method to handle problems in different applications.

• Journal of Ocean Engineering and Technology
• /
• v.7 no.2
• /
• pp.150-161
• /
• 1993

원초변수를 이용한 Navier-Stokes 방정식의 수치계산기법을 개발하고, 이를 응용하여 backward facing step의 층류 유동을 계산하였다. 직교좌표계에서의 비압축성 Navier-Stokes방정식을 풀기위해 시간과 공간항을 2차 정도의 유한 차분을 사용하여 이산화하였고 비교차격자계를 사용하여 양해법으로 수치 계산하였다. 운동량방정식과 연속방정식으로 부터 유도된 압력방정식(pressure-poisson equation)을 이용하여 무발산 조건을 만족시켰ㄲ다. Backward facing step의 층류 유동을 100.$\leq$R$_e$$\leq$1000 범위에 대해서 수치 계산하였으며 실험결과와 잘 일치하는 결과를 구할 수 있었다. 특히 step뒤에서 생기는 박리구간의 길이는 다른 계산결과들보다 실험치에 가까운 값을 얻을 수 있었으며, Re가 600보다 클때는 위쪽 벽에 또 다른 박리 유동이 발생되는 현상이 예측되었다.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.84-91
• /
• 2008

Generalized hydrodynamic (GH) theory for multi-species gas and the computational models are developed for the numerical simulation of hypersonic rarefied gas flow on the basis of Eu's GH theory. The rotational non-equilibrium effect of diatomic molecules is taken into account by introducing excess normal stress associated with the bulk viscosity. The numerical model for the diatomic GH theory is developed and tested. Moreover, with the experience of developing the dia-tomic GH computational model, the GH theory is extended to a multi-species gas including 5 species; O$_2$, N$_2$, NO, O, N. The multi-species GH model includes diffusion relation due to the molecular collision and thermal phenomena. Two kinds of GH models are developed for an axisymmetric flow solver. By compar-ing the computed results of diatomic and multi-species GH theories with those of the Navier-Stokes equations and the DSMC results, the accuracy and physical consistency of the GH computational models are examined.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.84-91
• /
• 2008

Generalized hydrodynamic (GH) theory for multi-species gas and the computational models are developed for the numerical simulation of hypersonic rarefied gas flow on the basis of Eu's GH theory. The rotational non-equilibrium effect of diatomic molecules is taken into account by introducing excess normal stress associated with the bulk viscosity. The numerical model for the diatomic GH theory is developed and tested. Moreover, with the experience of developing the dia-tomic GH computational model, the GH theory is extended to a multi-species gas including 5 species; $O_2,\;N_2$, NO, O, N. The multi-species GH model includes diffusion relation due to the molecular collision and thermal phenomena. Two kinds of GH models are developed for an axisymmetric flow solver. By compar-ing the computed results of diatomic and multi-species GH theories with those of the Navier-Stokes equations and the DSMC results, the accuracy and physical consistency of the GH computational models are examined.