• 제목/요약/키워드: Higher-Order Finite Difference Method

검색결과 50건 처리시간 0.026초

DNS에 의한 원주후류에 대한 유동해석 (The Flow Analysis of Past Flow a Circular Cylinder By Direct Numerical Simulation)

  • 강신정;;;김유택;이영호
    • 한국마린엔지니어링학회:학술대회논문집
    • /
    • 한국마린엔지니어링학회 2001년도 춘계학술대회 논문집
    • /
    • pp.52-57
    • /
    • 2001
  • Laminar two-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for the low Reynolds number (Re=164∼280). 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. The convection term is applied by the 7th order up wind scheme and the pressure and viscosity terms are applied by the 4th order central difference. The grid system makes use of the regular grid system and it is generated by an equation. The calculated results of drag coefficients, lift coefficients, pressure distributions, and vorticity contours and other information are compared with experimental and numerical ones. These results obtained by the present DNS show good agreement with the previous studies.

  • PDF

Non-local 이론을 적용한 단열전단밴드의 국부화에 대한 연구 (A Study of Localization for Adiabatic Shear Band Using Non-local Theory)

  • 이용성;이병섭;황두순;윤수진;홍성인
    • 한국소성가공학회:학술대회논문집
    • /
    • 한국소성가공학회 2001년도 춘계학술대회 논문집
    • /
    • pp.205-208
    • /
    • 2001
  • Localized shear band is investigated through the analysis of one-dimensional model for simple shearing deformation of thermally rate dependent material. Generally mesh size or interval of nodes play an important role in determining the overall flow behavior of the material. In order to observe these size effects we adapted non-local theory by including higher order strain gradients of the equivalent strain into the constitutive equation for the flow stress. for the ease of convergence and numerical stability the inplicit finite difference scheme is employed.

  • PDF

강체모드분리와 급수전개를 통한 준해석적 민감도 계산 방법의 개선에 관한 연구(I) - 정적 문제 - (A Refined Semi-Analytic Sensitivity Study Based on the Mode Decomposition and Neumann Series Expansion (I) - Static Problem -)

  • 조맹효;김현기
    • 대한기계학회논문집A
    • /
    • 제27권4호
    • /
    • pp.585-592
    • /
    • 2003
  • Among various sensitivity evaluation techniques, semi-analytical method(SAM) is quite popular since this method is more advantageous than analytical method(AM) and global finite difference method(FDM). However, SAM reveals severe inaccuracy problem when relatively large rigid body motions are identified fur individual elements. Such errors result from the numerical differentiation of the pseudo load vector calculated by the finite difference scheme. In the present study, an iterative method combined with mode decomposition technique is proposed to compute reliable semi-analytical design sensitivities. The improvement of design sensitivities corresponding to the rigid body mode is evaluated by exact differentiation of the rigid body modes and the error of SAM caused by numerical difference scheme is alleviated by using a Von Neumann series approximation considering the higher order terms for the sensitivity derivatives.

Extremely Low Numerical Dispersion FDTD Method Based on H(2, 4) Scheme for Lossy Material

  • Oh, Il-Young;Hong, Yongjun;Yook, Jong-Gwan
    • Journal of electromagnetic engineering and science
    • /
    • 제13권3호
    • /
    • pp.158-164
    • /
    • 2013
  • This paper expands a previously proposed optimized higher order (2, 4) finite-difference time-domain scheme (H(2, 4) scheme) for use with lossy material. A low dispersion error is obtained by introducing a weighting factor and two scaling factors. The weighting factor creates isotropic dispersion, and the two scaling factors dramatically reduce the numerical dispersion error at an operating frequency. In addition, the results confirm that the proposed scheme performs better than the H(2, 4) scheme for wideband analysis. Lastly, the validity of the proposed scheme is verified by calculating a scattering problem of a lossy circular dielectric cylinder.

Numerical study on the performance of semicircular and rectangular submerged breakwaters

  • Barzegar, Mohammad;Palaniappan, D.
    • Ocean Systems Engineering
    • /
    • 제10권2호
    • /
    • pp.201-226
    • /
    • 2020
  • A systematic numerical comparative study of the performance of semicircular and rectangular submerged breakwaters interacting with solitary waves is the basis of this paper. To accomplish this task, Nwogu's extended Boussinesq model equations are employed to simulate the interaction of the wave with breakwaters. The finite difference technique has been used to discretize the spatial terms while a fourth-order predictor-corrector method is employed for time discretization in our numerical model. The proposed computational scheme uses a staggered-grid system where the first-order spatial derivatives have been discretized with fourth-order accuracy. For validation purposes, five test cases are considered and numerical results have been successfully compared with the existing analytical and experimental results. The performances of the rectangular and semicircular breakwaters have been examined in terms of the wave reflection, transmission, and dissipation coefficients (RTD coefficients) denoted by KR, KT, KD. The latter coefficient KD emerges due to the non-energy conserving KR and KT. Our computational results and graphical illustrations show that the rectangular breakwater has higher reflection coefficients than semicircular breakwater for a fixed crest height, but as the wave height increases, the two reflection coefficients approach each other. un the other hand, the rectangular breakwater has larger dissipation coefficients compared to that of the semicircular breakwater and the difference between them increases as the height of the crest increases. However, the transmission coefficient for the semicircular breakwater is greater than that of the rectangular breakwater and the difference in their transmission coefficients increases with the crest height. Quantitatively, for rectangular breakwaters the reflection coefficients KR are 5-15% higher while the diffusion coefficients KD are 3-23% higher than that for the semicircular breakwaters, respectively. The transmission coefficients KT for rectangular breakwater shows the better performance up to 2.47% than that for the semicircular breakwaters. Based on our computational results, one may conclude that the rectangular breakwater has a better overall performance than the semicircular breakwater. Although the model equations are non-dissipative, the non-energy conserving transmission and reflection coefficients due to wave-breakwater interactions lead to dissipation type contribution.

MLS 차분법의 결정 변수에 따른 정확도 분석 및 혼합변분이론을 통한 미분근사 성능향상 (On the Improvement of the Accuracy of Higher Order Derivatives in the MLS(Moving Least Square) Difference Method via Mixed Formulation)

  • 김현영;김준식
    • 한국전산구조공학회논문집
    • /
    • 제33권5호
    • /
    • pp.279-286
    • /
    • 2020
  • 본 연구에서는 점근해석 및 논로컬 이론에서 요구하는 4차 이상의 고차 미분근사를 수행하기 위하여 계방정식에 혼합변분이론을 적용하여 MLS 차분법으로부터 구해지는 고차 미분근사의 정확도를 큰 폭으로 향상시킨다. 또한, MLS 차분법에 존재하는 세 가지 조건변수에 따른 고차미분근사의 정확도를 비교·분석한다. 혼합변분이론의 합응력을 후처리하여 변위의 미분을 근사할 경우 기존의 변위장 기반 계방정식의 차분 결과에 비해 미분 차수가 2차 낮아진다. 해석 범위내 절점 수가 과도하게 많거나 기저 차수가 클 경우 MLS 차분법의 영향영역 내에서 과적합(overfitting)이 발생한다. 또한 영향영역이 최적 범위 이상으로 넓어질 경우 근사의 정확도가 떨어진다. 위 내용을 사인 하중을 받는 단순지지보 수치예제로부터 확인하였다.

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

  • 강신정;타나하시 마모루;미야우치 토시오;이영호
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 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.

  • PDF

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

  • 강신정;타나하시 마모루;미야우치 토시오;모장오;이영호
    • 대한기계학회:학술대회논문집
    • /
    • 대한기계학회 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.

  • PDF

New implicit higher order time integration for dynamic analysis

  • Alamatian, Javad
    • Structural Engineering and Mechanics
    • /
    • 제48권5호
    • /
    • pp.711-736
    • /
    • 2013
  • In this paper new implicit time integration called N-IHOA is presented for dynamic analysis of high damping systems. Here, current displacement and velocity are assumed to be functions of the velocities and accelerations of several previous time steps, respectively. This definition causes that only one set of weighted factors is calculated from the Taylor series expansion which leads to a simple approach and reduce the computational efforts. Moreover a comprehensive study on stability of the proposed method i.e., N-IHOA compared with IHOA integration which is performed based on amplification matrices proves the ability of the N-IHOA in high damping vibrations such as control systems. Also, wide range of numerical examples which contains single/multi degrees of freedom, damped/un-damped, free/forced vibrations from finite element/finite difference demonstrate that the accuracy and efficiency of the proposed time integration is more than the common approaches such as the IHOA, the Wilson-${\theta}$ and the Newmark-${\beta}$.

On Bounds for Moments of Unimodal Distributions

  • Sharma, R.;Bhandaria, R.
    • Communications for Statistical Applications and Methods
    • /
    • 제21권3호
    • /
    • pp.201-212
    • /
    • 2014
  • We provide a simple basic method to find bounds for higher order moments of unimodal distributions in terms of lower order moments when the random variable takes value in a given finite real interval. The bounds for moments in terms of the geometric mean of the distribution are also derived. Both continuous and discrete cases are considered. The bounds for the ratio and difference of moments are obtained. The special cases provide refinements of several well-known inequalities, such as Kantorovich inequality and Krasnosel'skii and Krein inequality.