• 한국전자파학회논문지
• /
• 제21권7호
• /
• pp.805-814
• /
• 2010

기존 Crank-Nicolson FDTD 기법(CN FDTD 기법)의 비등방성 분산 특성을 개선하기 위한 CN ID-FDTD 기법을 제안하였다. 제안한 CN ID-FDTD 기법은 공간 미분 연산을 위해 기존 CN FDTD 기법의 centered 유한 차분식 (Finite Difference equation: FD 연산식)이 아닌 isotropic-dispersion 유한 차분식(ID-FD 연산식)$^{[1],[2]}$을 이용한다. 본 논문에서는 손실 매질에 대한 CN ID-FDTD 기법의 분산 관계식을 유도하였고, 이 분산 관계식을 이용해 ID-FD 연산식에서 분산 오차(dispersion error)를 줄이는 가중치(weighting factor)와 보정값(scaling factor)을 제시하였다. 그리고 해석 결과의 정확성 비교를 통해 CN ID-FDTD 기법에서는 기존 CN FDTD 기법의 단점이었던 비등방성 분산 오차가 확연하게 감소하는 것을 확인하였다.

• 대한전자공학회논문지
• /
• 제27권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년도 추계 학술대회논문집
• /
• 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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• 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.

• 한국전산유체공학회지
• /
• 제6권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년도 춘계 학술대회논문집
• /
• 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.

• 한국해안해양공학회지
• /
• 제6권3호
• /
• pp.260-265
• /
• 1994

대규모 흐름이 존재하는 불규칙한 해역에서 새로운 타원형 파동방정식을 유도하고, 유한차분법을 이용한 효율적인 수치모형을 개발하였다. 이때 청원형 방정식은 초기식 문제의 해법과 유사한 방법을 사용하여 해를 구하였다. 이 방법은 불규칙파의 변형을 계산하는 데 특히 효과적이며 수리모형 실험결과(Hiraishi, 1991)와 잘 일치하였다. 마지막으로 수중천퇴가 존재하는 완경사 해역에서 파랑과 흐름의 상호작용에 의한 수치해를 예시하였다.

• International Journal of Aeronautical and Space Sciences
• /
• 제7권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.

• 전자공학회논문지D
• /
• 제34D권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.

• 전기학회논문지
• /
• 제60권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.