• Title/Summary/Keyword: Spectral finite difference scheme

Search Result 12, Processing Time 0.026 seconds

A Study on Natural Convection from Two Cylinders in a Cavity

  • Mochimaru Yoshihiro;Bae Myung-Whan
    • Journal of Mechanical Science and Technology
    • /
    • v.20 no.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.

Numerical analysis of natural convection heat transfer from a fin in parallel enclosure

  • Bae, Myung-Whan;Mochimaru, Yoshihiro
    • Journal of Advanced Marine Engineering and Technology
    • /
    • v.40 no.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.

2D Crank-Nicolson FDTD Method Based on Isotropic-Dispersion Finite Difference Equation for Lossy Media (손실 매질에 대한 Isotropic-Dispersion 유한 차분식의 2D Crank-Nicolson FDTD 기법)

  • Kim, Hyun;Koh, Il-Suek;Yook, Jong-Gwan
    • 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.

On the Suitability of Centered and Upwind-Biased Compact Difference Schemes for Large Eddy Smulations (II) - Static Error Analysis - (LES에서 중심 및 상류 컴팩트 차분기법의 적합성에 관하여 (II) - 정적 오차 해석 -)

  • Park, No-Ma;Yoo, Jung-Yul;Choi, Hae-Cheon
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.27 no.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.

Secondary Instability in the Wake of a Circular Cylinder (원주 후류에서의 2차적 불안정성)

  • KNAG S. J.;TANAHASHI M.;MIYAUCHI T.;LEE Y. H.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 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.

  • PDF

Three-Dimensional Transition in the Wake of a Circular Cylinder By Direct Numerical Simulation (DNS에 의한 원주 후류에서의 3차원 천이)

  • Knag, S.J.;Tanahashi, M.;Miyauchi, T.;Mo, J.O.;Lee, Y.H.
    • 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.

  • PDF

CALCULATION METHODS OF SOLAR ATMOSPHERIC MODEL (태양대기모델 계산법)

  • KIM KAP-SUNG
    • Publications of The Korean Astronomical Society
    • /
    • v.15 no.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.

  • PDF

Enhancing the Performance of Coherent Sources SAC OCDMA Networks via Spatial Multiplexing

  • Alhassan, Ahmed M.;Badruddin, Nasreen;Saad, Naufal M.;Aljunid, Syed A.
    • Journal of the Optical Society of Korea
    • /
    • v.17 no.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.

Direct Numerical Simulation of the Flow Past an Oscillating Circular Cylinder (진동하는 원주주위 유동의 직접수치해석)

  • Kang S. J.;Tanahashi M.;Miyauchi T.;Lee Y. H.
    • 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.

  • PDF

Direct Numerical Simulation of the Flow Past an Oscillating Circular Cylinder (진동하는 원주주위 유동의 직접수치해석)

  • KANG Shin-Jeong;TANAHASHI Mamoru;MIYAUCHI Toshio;NAM Cheong-Do;LEE Young-Ho
    • 한국전산유체공학회:학술대회논문집
    • /
    • 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.

  • PDF