• 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 Institute of Telematics and Electronics
• /
• v.27 no.12
• /
• pp.1925-1932
• /
• 1990

In this paper, we consider the effect of diffraction due to the finite-sized aperture and propose the new correction method of errors in measurements of attenuation coefficient owing to the diffraction effect. In the existing correction method, we obtained the attenuation coefficient after correct the spectrum at each depth. However, this paper obtain the attenuation coefficient using lg-spectral difference approach and then correct errors. As a proposed method is not correction for the spectrum at each depth but the difference spectrum, we reduce the calculation. Also the correction is performed through the total frequency range, the accurate attenuation coefficient in whole bandwidth is produced.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.26-34
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Navier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to 25% of the cylinder diameter and in the case of the lock-in region it is 60%. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.181-188
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Wavier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to $25\%$ of the cylinder diameter and in the case of the lock-in region it is $60\%$. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.3
• /
• pp.260-265
• /
• 1994

A new set of elliptic wave equations describing the deformations of irregular waves on a large-scale current field in water of irregular depth is given, and using finite difference scheme an efficient numerical method is also presented. The elliptic equations are solved in a similar way to Initial value problem. This method is extensively used for the calculation of wave spectral transformation. and computation results agree very well with experimental data (Hiraishi, 1991). Finally numerical examples are presented concerning the interactions between waves and currents over a mildly sloping beach and also over a mound.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.155-174
• /
• 2006

Parabolized stability equations for compressible flows in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Compressible and incompressible flat plate flow stability under two-dimensional and three¬dimensional disturbances has been investigated to test the present code. Results of the present computation are found to be in good agreement with the multiple scale analysis and DNS data. Stability calculation results by the present PSE code for compressible boundary layer at Mach numbers ranging from 0.02 to 1.5 are also presented and are again seen to be as accurate as the spectral method.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.34D no.7
• /
• pp.79-87
• /
• 1997

Untill now, the most well-known cofiguration for waveguide grating router(WGR) is composed of radiative star couplers and arrayed waveguide grating(AWG), which usually suffer form the rdiation loss of around 3dB or more. Therefore, te improved design of WGRs is needed to reduce the loss. In ths paper, we propose a novel WGR composed of multimode interference couplers which have good unifiormity, fabrication tolerance, and very low excess loss, and suggest the efficient algorithm to find the proper path length differences of AWG for given channel spacing and channel assignment to each output prot. The simulated spectral responses of the proposed WGR using the finite difference beam propagation method (BPM) show that the excess loss is less than 0.3dB and the crosstalk less than -25dB in case of 4x4 WGR, and the excess loss less than 0.4dB and the crosstalk less than -25dB in case of 8x8 WGR for all the channel wavelengths.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.60 no.4
• /
• pp.817-822
• /
• 2011

We propose a novel ring structure based on the stadium-shaped ring resonator (SSRR) with dual photonic crystal microcavity (DPCM) for biosensor and analyzed the sensing characteristics. The Q-factor of the photonic crystal microcavity (PCM) can be significantly enhanced when the PCM or DPCM has the same resonance condition as the SSRR. The simulation results show that the Q-factor of the SRR with DPCM was increased by three times in comparison with single PCM structure. We also defined a mutual interference between two PCMs. Assuming a detectable spectral resolution of 10 picometers, a refractive index resolution of $3.03\times10-5$ can be measured on the SSRR-DPCM.