• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5A
• /
• pp.457-465
• /
• 2009

This study presents a moving least squares (MLS) finite difference method for solving elastic crack problems with stress singularity at the crack tip. Near-tip functions are intrinsically employed in the MLS approximation to model near-tip field inducing singularity in stress field. employment of the functions does not lose the merit of the MLS Taylor polynomial approximation which approximates the derivatives of a function without actual differentiating process. In the formulation of crack problem, computational efficiency is considerably improved by taking the strong formulation instead of weak formulation involving time consuming numerical quadrature Difference equations are constructed on the nodes distributed in computational domain. Numerical experiments for crack problems show that the intrinsically enriched MLS finite difference method can sharply capture the singular behavior of near-tip stress and accurately evaluate stress intensity factors.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.24 no.7
• /
• pp.748-751
• /
• 2013

We perform a wideband radiated pulse coupling analysis of simple building structures using the finite-deference time-domain(FDTD) method. Toward this purpose, the building structures composed of concrete and window materials are assumed and we numerically model the electrical properties of each material. In this work, we apply a dispersive FDTD algorithm for the electromagnetic analysis of building structures and investigate their shielding effectiveness in the frequency range of 50 MHz to 1 GHz.

• Fire Science and Engineering
• /
• v.6 no.1
• /
• pp.17-22
• /
• 1992

We have derived the general transfer equation for governing the continuity, energy transfer, mass and momentum transfer, and turbulent energy dissipation rate within the fire compartment which has the 800t fire source at the center of the floor. The governing transfer equations have been descretized using the finite volume approach and numerically experimented under the SIMPLE algorithm. In order for the SIMPLE algorithm approach to be physically reliable, the test results are compared with those of Morita's SOR Method using Conjugate Residual Method and found to be close to physical values though the computational convergence time still remains to be upgraded. The treatment of source terms in the system of finite difference equations has been critical in order to converge the governing equations within the appropriate time steps. The criteria of convergence allowance for the whole domain have been checked and the sudden change of the non-linear effects from the source term have been avoided. The criteria has been allowed to be for 5$\times$10$^{-5}$ .

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.17 no.2 s.105
• /
• pp.158-164
• /
• 2006

In this paper, we proposed a new algorithm of the inverse scattering for the reconstruction of unknown dielectric scatterers using the finite-difference time-domain method and the design sensitivity analysis. We introduced the design sensitivity analysis based on the gradient information for the fast convergence of the reconstruction. By introducing the adjoint variable method for the efficient calculation, we derived the adjoint variable equation. As an optimal algorithm, we used the steepest descent method and reconstructed the dielectric targets using the iterative estimation. To verify our algorithm, we will show the numerical examples for the two-dimensional $TM^2$ cases.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.17 no.2 s.105
• /
• pp.165-170
• /
• 2006

A two dimensional(2-D) finite-difference time-domain(FDTD) method based on a novel finite difference scheme is developed to eliminate the numerical dispersion errors. In this paper, numerical dispersion and stability analysis of the new scheme are given, which show that the proposed method is nearly dispersionless, and stable for a larger time step than the standard FDTD method.

• Proceedings of the IEEK Conference
• /
• 2002.07a
• /
• pp.137-140
• /
• 2002

In this paper, a new algorithm for two-dimensional (2-D) finite-difference time-domain (FDTD) method is presented. By this new method, the maximum time step size can be increased over the Courant-Friedrich-Levy (CFL) condition restraint. This new algorithm is adapted from an Alternating-Direction Implicit FDTD (ADI-FDTD) method. However, unlike the ADI-FDTD algorithm. the alternation is performed with respect to the blocks of fields rather than with respect to each respective coordinate direction. Moreover. this method can be efficiently simulated with parallel computation. and it is more efficient than the conventional FDTD method in terms of CPU time. Numerical formulations are shown and simulation results are presented to demonstrate the effectiveness and efficiency of our proposed method.

• Journal of applied mathematics & informatics
• /
• v.35 no.3_4
• /
• pp.277-302
• /
• 2017

In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of reaction-diffusion type of third order Ordinary Differential Equations (ODEs). The SPBVP is reduced into a weakly coupled system of one first order and one second order ODEs, one without the parameter and the other with the parameter ${\varepsilon}$ multiplying the highest derivative subject to suitable initial and boundary conditions, respectively. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. The weakly coupled system is decoupled by replacing one of the unknowns by its zero-order asymptotic expansion. Finally the present numerical method is applied to the decoupled system. In order to get a numerical solution for the derivative of the solution, the domain is divided into three regions namely two inner regions and one outer region. The Shooting method is applied to two inner regions whereas for the outer region, standard finite difference (FD) scheme is applied. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing. The main advantage of this method is that due to decoupling the system, the computation time is very much reduced.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11a
• /
• pp.620-628
• /
• 2005

Spatial control of sound is essential to deliver better sound to the listener's position in space. As it can be experienced in many listening environments, the quality of sound can not be manifested over every position in a hall. This motivates us to control sound in a region we select. The primary focus of the developed method has to do with the brightness and contrast of acoustic image in space. In particular, the acoustic brightness control seeks a way to increase loudness of sound over a chosen area, and the contrast control aims to enhance loudness difference between two neighboring regions. This enables us to make two different kinds of zone - the zone of quiet and the zone of loud sound - at the same time. The other perspective of this study is on the direction of sound. It is shown that we can control the direction of perceived sound source by focusing acoustic energy in wavenumber domain. To begin with, the proposed approaches are formulated for pure-tone case. Then the control methods are extended to a more general case, where the excitation signal has broadband spectrum. In order to control the broadband signal in time domain, an inverse filter design problem is defined and solved in frequency domain. Numerical and experimental results obtained in various conditions certainly validate that the acoustic brightness, acoustic contrast, direction of wave front can be manipulated for some finite region in space and time.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.10 no.6
• /
• pp.847-854
• /
• 1999

In this paper, we have analysed frequency-dispersion characteristics of higher-order modes for uniform planar transmission lines, using the 2-dimensional finite-difference time-domain method. The method presented in this paper uses both informations of amplitude and phase of the electromagnetic spectrum to determine resonant frequencies, while methods previously reported use the magnitude only. This algorithm is very useful when a resonant mode has a relatively small magnitude, where the identification of the resonant mode is quite difficult. Numerical results show that a strip line supports few higher-order modes within the frequency range of 20 GHz, but there occur many higher-order modes in the structure of grounded coplanar waveguide, where resonant frequencies of the first higher-order mode is very close to those of the fundamental mode and there occur lots of very adjacent higher-order modes. As in this example, for the analysis of planar transmission lines which support many resonant modes very close each other, the method presented in this paper can be applied very efficiently.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.9 no.5
• /
• pp.678-686
• /
• 1998

This paper describes an anisotropic perfectly matched layer (APML) for mesh truncation of the Finite Difference Time Domain(FDTD) method. The APML method can classified by a split type and an unsplit type, in case of the split type be made up 12 equations or 8 equations largely, and in case of the unsplit type be made of 6 equations. Therefore the latter is more simple as compare with the former. For the APML method presented in this paper is the latter, is directly derived from the time domain equations of Maxwell and extend to the three dimensional problem for the method suggested by Chen. Especially, in the edge and corner parts except the planes, the APML method effectively treated as compound with the Chen's method and Gedney's method newly. The results of the numerical method in this paper show the radiation patterns and the time responses of electromagnetic fields of the wire antennas according to wavelengths and the APML results are compared with Mur's first order absorbing boundary condition and Kraus's analytical results. Eventually, Mur's first order condition have many errors at the edge and corner. On the other hand, in comparison with Kraus's analytical results, it is quite good agreement, and the validity of present method is confirmed.