• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.14 no.1
• /
• pp.38-44
• /
• 2010

In this paper, the higher order FDTD(Finite-Difference Time-Domain) method is used to obtain the frequency response characteristics of the cylindrical metamaterial slab. FDTD method is one of strongest electromagnetic numerical method which is widely used to analyze the metamaterial structure because of its simplicity and the dispersive FDTD equation which has the dispersive effective dielectric constant and permeability are derived to analyze the metamaterials. This derived dispersive FDTD equation has no errors in analyzing the dielectric materials but there are some time and frequency errors in case of analyzing the metamaterials. We used the higher order FDTD method to obtain the accurate frequency response of the metamaterials. Comparisons between the dispersive FDTD method and the higher order FDTD method are performed in this paper also. From the results, we concluded that more accurate frequency response for various metamaterials applications can be obtained using the proposed method in this paper.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.411-420
• /
• 2009

This study presents an extended finite difference method based on moving least squares(MLS) method for solving potential problems with interfacial boundary. The approximation constructed from the MLS Taylor polynomial is modified by inserting of wedge functions for the interface modeling. Governing equations are node-wisely discretized without involving element or grid; immersion of interfacial condition into the approximation circumvents numerical difficulties owing to geometrical modeling of interface. Interface modeling introduces no additional unknowns in the system of equations but makes the system overdetermined. So, the numbers of unknowns and equations are equalized by the symmetrization of the stiffness matrix. Increase in computational effort is the trade-off for ease of interface modeling. Numerical results clearly show that the developed numerical scheme sharply describes the wedge behavior as well as jumps and efficiently and accurately solves potential problems with interface.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.19 no.4
• /
• pp.427-435
• /
• 2008

In this paper, a stable marching-on in time(MOT) method with a time domain combined field integral equation(CFIE) is presented to obtain the transient scattering response from arbitrarily shaped three-dimensional conducting bodies. This formulation is based on a linear combination of the time domain electric field integral equation(EFIE) with the magnetic field integral equation(MFIE). The time derivatives in the EFIE and MFIE are approximated using a central finite difference scheme and other terms are averaged over time. This time domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. Numerical results with the proposed MOT scheme are presented and compared with those obtained from the conventional method and the inverse discrete Fourier transform(IDFT) of the frequency domain CFIE solution.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.9 no.6
• /
• pp.761-767
• /
• 1998

The Perfectly Matched Layer(PML) provide good performance in absorption over a wide frequency range and is an appropriate ABC for waveguides with high dispersion. In this paper, a novel algorithm is proposed to improve the computational efficiency of the PML. In the input and output ports, the fields are decomposed into a series of modes, and then an appropriate ABC is applied to each mode. CPU time and memory storage requirements are greatly reduced, since the computational region is analyzed in one dimension. A WG-90 rectangular waveguide with a thick asymmetric iris is analyzed by Finite-Difference Time-Domain(FDTD) simulations with the conventional PML and the proposed one-dimensional (1-D) PML. Numerical results show that the computational efficiency is significantly improved by the proposed method.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.4
• /
• pp.38-47
• /
• 1995

The three-dimensional multi-layer computer simulation technology for tidal current developed in the previous study is updated to a new version. many improvements are achieved by following changes : (1) No-reflection condition is adopted instead of no-gradient condition as an open boundary condition. (2) Time marching algorithm is changed so that velocity and pressure(surface movement) might be salved in turn at different time step (3) Convection term in equation of motion is estimated by upwind differencing scheme instead of central differencing. The stability is improved considerably and the steady state is achieved within 2 tidal periods which is about 3 times shorter than that of the old version. Moreover, fluctuations in time disappeared by introducing the new time marching technique. An application to the real near shore area(near Inchon harbor) is performed by the new version. Simulated results are compared with those by the simulation total developed in the University of Tokyo. Validity and effectiveness of the two simulation technologies are chocked through the comparative research works.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.16 no.2
• /
• pp.211-221
• /
• 2018

Nodal transport methods are proposed for solving the simplified even-parity neutron transport (SEP) equation. These new methods are attributed to the success of existing nodal diffusion methods such as the Polynomial Expansion Nodal and the Analytic Function Expansion Nodal Methods, which are known to be very effective for solving the neutron diffusion equation. Numerical results show that the simplified even-parity transport equation is a valid approximation to the transport equation and that the two nodal methods developed in this study also work for the SEP transport equation, without conflict. Since accuracy of methods is easily increased by adding node unknowns, the proposed methods will be effective for coarse mesh calculation and this will also lead to computation efficiency.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.1
• /
• pp.45-53
• /
• 1991

The nonlinear forward diffraction of a modulated wave train by a thin wedge has been studied analytically. Since the physical variables involved in the problem have vastly different scales, the multiple scale expansion method has been used to obtain an approximate solution. To simplify the problem. the wedge is assumed to be thin and the parabolic approximation is utilized. The wave evolution can be described by a kind of the cubic Schrodinger equation. which consists of the linear time evolution. the lateral dispersion and the nonlinearity. Numerical results indicate that the nonlinearity. which it defined by the ratio of the ratio of the incident wave to the wedge angle. governs the amplitude and the stability of diffracted waves. The instability of dirffracted waves becomes more pronounced as the nonlinearity increases and the modulation ratio decreases. It is also found that the stem waves. initially developed along the wedge. can not sustain for a long time.

• Journal of the Korean Geotechnical Society
• /
• v.24 no.8
• /
• pp.111-123
• /
• 2008

Geotechnical engineering problems are characterized by many sources of uncertainty. Some of these sources are connected to the uncertainties of soil properties involved in the analysis. In this paper, a numerical procedure for a probabilistic analysis that considers the spatial variability of soil properties is presented to study the response of spatially random soil. The approach integrates a commercial finite difference method and random field theory into the framework of a probabilistic analysis. Two-dimensional non-Gaussian random fields are generated based on a Karhunen-$Lo{\grave{e}}ve$ expansion in a fashion consistent with a specified marginal distribution function and an autocorrelation function. A Monte Carlo simulation is then used to determine the statistical response based on the random fields. A series of analyses were performed to study the effects of uncertainty due to the spatial heterogeneity on the settlement and bearing capacity of a rough strip footing. The simulations provide insight into the application of uncertainty treatment to the geotechnical problem and show the importance of the spatial variability of soil properties with regard to the outcome of a probabilistic assessment.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.39 no.2
• /
• pp.14-26
• /
• 2002

A new model for degenerate semiconductor quantum well devices was developed. In this model, the multi-subband Boltzmann transport equation was formulated by applying the Pauli exclusion principle and coupled to the Schrodinger and Poisson equations. For the solution of the resulted nonlinear system, the finite difference method and the Newton-Raphson method was used and carrier energy distribution function was obtained for each subband. The model was applied to a Si MOSFET inversion layer. The results of the simulation showed the changes of the distribution function from Boltzmann like to Fermi-Dirac like depending on the electron density in the quantum well, which presents the appropriateness of this modeling, the effectiveness of the solution method, and the importance of the Pauli -exclusion principle according to the reduced size of semiconductor devices.

• Journal of the Korean Electrochemical Society
• /
• v.15 no.3
• /
• pp.165-171
• /
• 2012

Here is implemented MATLAB program to analyze the characteristic curves of cyclic voltammetry which involves the multi-electron electrode reaction considered as key processes in electrochemical systems. For the electrochemical mass-transfer system, Fick's concentration equation subject to semi-infinite diffusion model for the boundary condition was discretized and solved by the explicit finite difference method. The resulting concentration values were converted into currents at each node by using Butler-Volmer equation. Based on the good agreement between the present numerical solution and the existing experimental results, effects of kinetic constants and CV scan rates on the reaction mechanism in multi-electron transfer processes were investigated effectively.