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

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.11 no.6
• /
• pp.1054-1061
• /
• 2007

In this paper, the comparison between UPML and Liao's absorbing boundary condition in the FDTD(Finite-Difference Time-Domain) method was performed for the analysis of the 2D cylindrical coordinate system. Generally, it is known as the absorbing characteristics of the UPML is bro than Liao's absorbing boundary condition in the 2D rectangular coordinate. The simulation results in this paper showed that Liao's original absorbing boundary condition is better than other two absorbing boundary conditions, Liao's modified condition and UPML. We concluded that more numerical, theoretical studies, simulations and verifications for various absorbing boundary conditions will be needed to get more accurate results for the design of useful 2D cylindrical microwave circuits.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.3
• /
• pp.23-27
• /
• 2008

The new conjunctive surface-subsurface flow model at a large scale was developed by using a 1-D Diffusion Wave (DW) model for surface flow interacting with the 3-D Volume Averaged Soil-moisture Transport (VAST) model for subsurface flow for the comprehensive terrestrial water and energy predictions in Land Surface Models (LSMs). A selection of numerical implementation schemes is employed for each flow component. The 3-D VAST model is implemented using a time splitting scheme applying an explicit method for lateral flow after a fully implicit method for vertical flow. The 1-D DW model is then solved by MacCormack finite difference scheme. This new conjunctive flow model is substituted for the existing 1-D hydrologic scheme in Common Land Model (CLM), one of the state-of-the-art LSMs. The new conjunctive flow model coupled to CLM is tested for a study domain around the Ohio Valley. The simulation results show that the interaction between surface flow and subsurface flow associated with the flow routing scheme matches the runoff prediction with the observations more closely in the new coupled CLM simulations. This improved terrestrial hydrologic module will be coupled to the Climate extension of the next-generation Weather Research and Forecasting (CWRF) model for advanced regional, continental, and global hydroclimatological studies and the prevention of disasters caused by climate changes.