• Proceedings of the Korea Electromagnetic Engineering Society Conference
• /
• 2005.11a
• /
• pp.183-186
• /
• 2005

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.

• Journal of electromagnetic engineering and science
• /
• v.8 no.4
• /
• pp.139-143
• /
• 2008

In this paper, a new scheme to calculate the weighting factor of the 2-D isotropic-dispersion finite difference time domain(ID-FDTD) is proposed. The weighting factor in [1] was formulated in free space, so that it may not be optimal in dielectric media. Therefore, the weighting factor was reformulated by considering the material properties and using the least mean square method. As a result, a minimum numerical dispersion error for any dielectric media is guaranteed.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.21 no.7
• /
• pp.805-814
• /
• 2010

The Crank-Nicolson isotropic-dispersion finite difference time domain(CN ID-FDTD) scheme is proposed based on isotropic-dispersion finite difference(ID-FD) $equation^{[1],[2]}$. The dispersion relation of CN ID-FDTD is derived for lossy media by solving the eigenvalue problem of iteration matrix in spatial spectral domain, in addition, the weighting factors and scaling factors of the CN ID-FDTD scheme are presented for low dispersion error. The CN ID-FDTD scheme makes the dispersion error drastically reduced and shows accurate numerical results compared to the conventional Crank-Nicolson FDTD method.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.22 no.4
• /
• pp.407-415
• /
• 2011

The stability condition and wideband characteristics of 3D ID-FDTD algorithm which has low dispersion error with isotropic dispersion are presented in this paper. 3D ID-FDTD method was proposed to improve the defect of the Yee FDTD such as the anisotropy and large dispersion error. The published paper calculated the stability condition of 3D ID-FDTD algorithm by using numerical method, however, it is thought that the examples were not sufficient to verify the stability condition. Thus, in this paper, various simulations are included in order to hold reliability under the conditions that the plane wave propagation is assumed with a single frequency and a wideband frequency. Also, the 3D ID-FDTD algorithm is compared to those that have the similar FDTD algorithm with ID-FDTD such as Forgy's method and non-standard FDTD method in a wideband. Finally, the radar cross section(RCS) for the large sphere with high dielectric constant is calculated.

• 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.

• Journal of electromagnetic engineering and science
• /
• v.13 no.3
• /
• pp.158-164
• /
• 2013

This paper expands a previously proposed optimized higher order (2, 4) finite-difference time-domain scheme (H(2, 4) scheme) for use with lossy material. A low dispersion error is obtained by introducing a weighting factor and two scaling factors. The weighting factor creates isotropic dispersion, and the two scaling factors dramatically reduce the numerical dispersion error at an operating frequency. In addition, the results confirm that the proposed scheme performs better than the H(2, 4) scheme for wideband analysis. Lastly, the validity of the proposed scheme is verified by calculating a scattering problem of a lossy circular dielectric cylinder.

• Smart Structures and Systems
• /
• v.21 no.2
• /
• pp.195-205
• /
• 2018

Research on Lamb wave-based damage identification in plate-like structures depends on precise knowledge of dispersive wave velocity. However, boundary reflections with the same frequency of interest and greater amplitude contaminate direct waves and thus compromise measurement of Lamb wave dispersion in different materials. In this study, non-reflecting boundaries were proposed in both numerical and experimental cases to facilitate time-frequency characterization of Lamb wave dispersion. First, the Lamb wave equations in isotropic and laminated materials were analytically solved. Second, the non-reflecting boundaries were used as a series of frames with gradually increased damping coefficients in finite element models to absorb waves at boundaries while avoiding wave reflections due to abrupt property changes of each frame. Third, damping clay was sealed at plate edges to reduce the boundary reflection in experimental test. Finally, the direct waves were subjected to the slant-stack and short-time Fourier transformations to calculate the dispersion curves of phase and group velocities, respectively. Both the numerical and experimental results suggest that the boundary reflections are effectively alleviated, and the dispersion curves generated by the time-frequency analysis are consistent with the analytical solutions, demonstrating that the combination of non-reflecting boundary and time-frequency analysis is a feasible and reliable scheme for characterizing Lamb wave dispersion in plate-like structures.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.25 no.3
• /
• pp.222-227
• /
• 2005

An efficient technique fur the calculation of guided wave dispersion curves in composite pipes is presented. The technique uses a forward-calculating variational calculus approach rather than the guess and iterate process required when using the more traditional partial wave superposition technique. The formulation of each method is outlined and compared. The forward-calculating formulation is used to develop finite element software for dispersion curve calculation. Finally, the technique is used to calculate dispersion curves for several structures, including an isotropic bar, two multi-layer composite bars, and a composite pipe.

• Earthquakes and Structures
• /
• v.4 no.1
• /
• pp.1-10
• /
• 2013

The present investigation deals with the analysis of wave motion in the layer of an anisotropic, initially stressed, fiber reinforced thermoelastic medium. Secular equations for symmetric and skew-symmetric modes of wave propagation in completely separate terms are derived. The amplitudes of displacements and temperature distribution were also obtained. Finally, the numerical solution was carried out for Cobalt and the dispersion curves, amplitudes of displacements and temperature distribution for symmetric and skew-symmetric wave modes are presented to evince the effect of anisotropy. Some particular cases are also deduced.

• Journal of electromagnetic engineering and science
• /
• v.10 no.2
• /
• pp.45-49
• /
• 2010

This paper presents an analytical formulation of the numerical reflection coefficient of the ID-FDTD scheme at the planar dielectric boundary for a TM wave incidence. The reflection coefficient is formulated in an approximate manner, and the accuracy of this method is numerically verified. The effective dielectric constant for a grid on the interface is obtained, and then reduced to that of the Yee scheme for a small cell size.