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

• Journal of the Korean Mathematical Society
• /
• v.58 no.3
• /
• pp.553-569
• /
• 2021

In this paper, we present and analyze a fully discrete numerical method for solving the time-fractional diffusion wave equation: ∂^β_tu - div(a∇u) = f, 1 < β < 2. We first construct a difference formula to approximate ∂^β_tu by using an interpolation of derivative type. The truncation error of this formula is of O(△t^2+δ-β)-order if function u(t) ∈ C^2,δ[0, T] where 0 ≤ δ ≤ 1 is the Hölder continuity index. This error order can come up to O(△t^3-β) if u(t) ∈ C³ [0, T]. Then, in combinination with the linear finite volume discretization on spatial domain, we give a fully discrete scheme for the fractional wave equation. We prove that the fully discrete scheme is unconditionally stable and the discrete solution admits the optimal error estimates in the H¹-norm and L₂-norm, respectively. Numerical examples are provided to verify the effectiveness of the proposed numerical method.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.6
• /
• pp.7-15
• /
• 2003

A Digital Wave Tank simulation technique, based on a finite-difference method and a modified marker-and-cell (MAC) algorithm, is applied in order to investigate the characteristics of nonlinear Tsunami propagations and their interactions with a 2D sloping beach, Ohkushiri Island, and to predict maximum wove run-up around the island. The Navier-Stokes (NS) and continuity equation are governed in the computational domain, and the boundary values are updated at each time step, by a finite-difference time-marching scheme in the frame of the rectangular coordinate system. The fully nonlinear, kinematic, free-surface condition is satisfied by the modified marker-density function technique. The near shore Tsunami is assumed to be a solitary wave, and is generated from the numerical wave-maker in the developed Digital Wave Tank. The simulation results are compared with the experiments and other numerical methods, based on the shallow-water wave theory.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.05a
• /
• pp.231-239
• /
• 2003

A Digital Wave Tank simulation technique based on a finite-difference method and a modified marker-and-cell (MAC) algorithm is applied to investigate the characteristics of nonlinear Tsunami propagations and their interactions with a 2D sloping beach and Ohkushiri island, and to predict maximum wave run-up around the island. The Navier-Stokes (NS) and continuity equation are governed in the computational domain and the boundary values updated at each time step by a finite-difference time-marching scheme in the frame of rectangular coordinate system. The fully nonlinear kinematic free-surface condition is satisfied by the modified marker-density function technique. The Nearshore Tsunami is assumed to be a solitary wave and generated from the numerical wavemaker in the developed Digital Wave Tank. The simulation results are compared with the experiments and other numerical methods based on the shallow-water wave theory.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.192-198
• /
• 1997

In the concurrent engineering, the CAD-based design model is necessary for multidisciplinary analysis and for computer-aided manufacturing (CAM). A shape and configuration design velocity field computation of structure has been developed using a computer-aided design (CAD) tool, Pro/ENGINEER. The design Parameterization with CAD tool is to characterize the change in dimensions and movements of geometric control points that govern the shape/orientation of the structural boundary. The boundary velocity is obtained by using a CAD-based finite difference method and the domain velocity field is obtained from finite element analysis (FEA) using the boundary displacement method. In this paper, the continuum configuration DSA for NVH problem, which requires the shape velocity field and the orientation velocity field at the same time, is developed using linear shape functions. For validation of continuum design sensitivity coefficients, design sensitivity coefficients are compared with the coefficients computed using by the finite difference method.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.25 no.9B
• /
• pp.1527-1533
• /
• 2000

Basically, a general characterization procedure based on the extraction of the characteristic impedance and propagation constant for analyzing a single MIS(Metal-Insulator-Semiconductor) transmission line is used. In this paper, an analysis for a new substrate shielding MIS structure consisting of grounded cross-bars at the interface between Si and SiO2 layer using the Finite-Difference Time-Domain (FDTD) method is presented. In order to reduce the substrate effects on the transmission line characteristics, a shielding structure consisting of grounded cross bar lines over time-domain signal has been examined. The extracted distributed frequency-dependent transmission line parameters and corresponding equivalent circuit parameters as well as quality factor have been examined as functions of cross-bar spacing and frequency. It is shown that the quality factor of the transmission line can be improved without significant change in the characteristic impedance and effectve dielectric constant.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.13 no.6
• /
• pp.590-598
• /
• 2002

We have numerically analyzed the electromagnetic wave scattering from randomly rough dielectric surfaces by using the finite volume time domain (FVTD) method. We have then shown that the present method yields a reasonable solution even at low-grazing angle (LGA). It should be noted that the number of sampling points per wavelength should be increased when more accurate numerical results are required, which fact makes the computer simulation impossible at LGA for a stable result. However, when the extrapolation is used for calculating the scattered field, an accurate result can be estimated. If we want to obtain the ratio of backscattering between the horizontal and vertical polarization, we do not need the large number of sampling points. The results are compared with the experimental data.

• Computational Structural Engineering
• /
• v.10 no.3
• /
• pp.217-223
• /
• 1997

This paper deals with the basic concepts of sensitivity analysis in a time domain for the transient response of a lumped mass system. Sensitivity analysis methods in thme domain for determining the effects of parameter changes on the response of a dynamic system by external excitation are presented. The parametric sensitivity of a lumped mass system in time domain can be investigated using different types of sensitivity functions, including first order standard and percentage sensitivity functions. These sensitivity functions are determined as a function of partial derivatives of system variables taken with respect to system parameters. In addition, we compared the results of the analytical method by direct method and those of numerical methods.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.39 no.4
• /
• pp.63-71
• /
• 2002

In this paper, we present microwave characteristics of traveling-wave photodetectors (TWPD) using the finite-difference time-domain method (FDTD). Current and voltage in the time domain are calculated by the FDTD. Also, characteristic impedance and propagation constant in frequency domain are obtained from the time-domain data. As the thickness of i-layer gets thicker and the waveguide width gets narrower, TWPD's show less microwave loss and higher velocity. The 50Ω impedance matching design is achieved for 2.4${\mu}{\textrm}{m}$ waveguide width and 1.2${\mu}{\textrm}{m}$ thickness of i-layer at 100 GHz.

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