• Journal of Mechanical Science and Technology
• /
• 제20권10호
• /
• pp.1773-1778
• /
• 2006

Steady-state natural convection heat transfer characteristics from cylinders in a multiply-connected bounded region are clarified. A spectral finite difference scheme (spectral decomposition of the system of partial differential equations, semi-implicit time integration) is applied in numerical analysis, with a boundary-fitted conformal coordinate system through a Jacobian elliptic function with a successive transformation to formulate a system of governing equations in terms of a stream function, vorticity and temperature. Multiplicity of the domain is expressed explicitly.

• Journal of Advanced Marine Engineering and Technology
• /
• 제40권5호
• /
• pp.412-417
• /
• 2016

A fin of finite width with infinitely small thickness is assumed to be placed horizontally between two horizontal parallel plates of infinite extension in the exactly central position. The lower plate and the half of the upper plate are kept at a constant lower temperature, and the remaining upper plate is kept at a constant higher temperature. The fin is also kept at a constant temperature (variable). Steady-state two-dimensional laminar natural convection is analyzed as a problem of boundary value under a boundary-fitted conformal mapping system, using a spectral finite difference scheme, with a condition of doubly-connectedness. The steady-state solution is obtained as a limit of the transient solution.

• 한국전자파학회논문지
• /
• 제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 기법의 단점이었던 비등방성 분산 오차가 확연하게 감소하는 것을 확인하였다.

• 대한기계학회논문집B
• /
• 제27권7호
• /
• pp.984-994
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a spectral, static error analysis. To investigate the effect of numerical dissipation on LES solutions, power spectra of discretization errors are evaluated for isotropic turbulence models in both continuous and discrete wavevector spaces. Contrary to the common belief, the aliasing errors from upwind-biased schemes are larger than those from comparable non-dissipative schemes. However, this result is the direct consequence of the definition of the power spectral density of the aliasing error, which poses the limitation of the static error analysis for upwind schemes.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 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.

• 천문학논총
• /
• 제15권spc2호
• /
• pp.65-71
• /
• 2000

We have investigated the numerical methods to calculate model atmosphere for the analysis of spectral lines emitted from the sun and stars. Basic equations used in our calculations are radiative transfer, statistical equilibrium and charge-particle conservations. Transfer equation has been solved to get emitting spectral line profile as an initial value problem using Adams-Bashforth-Moulton method with accuracy as high as 12th order. And we have calculated above non linear differential equations simultaneously as a boundary value problem by finite difference method of 3 points approximation through Feautrier elimination scheme. It is found that all computing programs coded by above numerical methods work successfully for our model atmosphere.

• Journal of the Optical Society of Korea
• /
• 제17권6호
• /
• pp.471-480
• /
• 2013

The beating of two or more lasers that have the same or a finite difference in the central frequencies, is the main source of noise in spectral amplitude coding optical code division multiple access (SAC OCDMA) systems. In this paper we adopt a spatial multiplexing (SM) scheme for SAC OCDMA systems to mitigate this beat noise. The results show that for different code weights and different data rates SM SAC can support a larger number of users than the conventional SAC for all different laser source configurations. However, SM SAC requires a more complex system than the conventional SAC, and almost twice as much optical component.

• 한국전산유체공학회지
• /
• 제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.