• Structural Engineering and Mechanics
• /
• 제17권6호
• /
• pp.735-749
• /
• 2004

This work presents a direct integration scheme, based on a fourth order finite difference approach, for elastodynamics. The proposed scheme was chosen as an alternative for attenuating the errors due to the use of the central difference method, mainly when the time-step length approaches the critical time-step. In addition to eliminating the spurious numerical oscillations, the fourth order finite difference scheme keeps the advantages of the central difference method: reduced computer storage and no requirement of factorisation of the effective stiffness matrix in the step-by-step solution. A study concerning the stability of the fourth order finite difference scheme is presented. The Finite Element Method and the Boundary Element Method are employed to solve elastodynamic problems. In order to verify the accuracy of the proposed scheme, two examples are presented and discussed at the end of this work.

• International Journal of Aeronautical and Space Sciences
• /
• 제18권4호
• /
• pp.816-826
• /
• 2017

In this paper, a generalized discretization scheme is proposed that can derive general-order finite difference equations representing the joint probability density function of dynamic response of stochastic systems. The various order of finite difference equations are applied to solutions of the Fokker-Planck-Kolmogorov (FPK) equation. The finite difference equations derived by the proposed method can greatly increase accuracy even at the tail parts of the probability density function, giving accurate reliability estimations. Compared with exact solutions and finite element solutions, the generalized finite difference method showed increasing accuracy as the order increases. With the proposed method, it is allowed to use different orders and types (i.e. forward, central or backward) of discretization in the finite difference method to solve FPK and other partial differential equations in various engineering fields having requirements of accuracy or specific boundary conditions.

• 전기전자학회논문지
• /
• 제23권3호
• /
• pp.817-825
• /
• 2019

마이크로프로세서를 이용한 제어알고리즘 구현에서 차분방정식이 매우 유용하게 사용된다. 샘플링 데이터로부터 미분 값을 추정하기 위해 전향, 후향 및 중심 차분 방식이 사용되어왔다. 차분 값을 계산하기 위해서는 차분계수가 매우 중요하다. 본 논문에서는 유한 차분 계수를 계산하기 위한 새로운 방식을 제시하고자 한다. 제안된 방식의 유효성을 입증하기 위해 RLS 알고리즘을 적용한 파라미터 추정에 대하여 적용하였다.

• 한국전산유체공학회지
• /
• 제11권1호
• /
• pp.36-42
• /
• 2006

The existing numerical approximations of convection flux, especially the spatial higher-order difference schemes, in unstructured cell-centered finite volume methods are examined in detail with each other and evaluated with respect to the accuracy through their application to a 2-D benchmark problem. Six higher-order schemes are examined, which include two second-order upwind schemes, two central difference schemes and two hybrid schemes. It is found that the 2nd-order upwind scheme by Mathur and Murthy(1997) and the central difference scheme by Demirdzic and Muzaferija(1995) have more accurate prediction performance than the other higher-order schemes used in unstructured cell-centered finite volume methods.

• Environmental Engineering Research
• /
• 제23권4호
• /
• pp.442-448
• /
• 2018

This paper presents the numerical simulation of advection-diffusion mechanism of BOD concentration which was used as an indicator of waste only in one flow-direction of waste stabilization ponds (1-dimension (1-D)). This model was represented in partial differential equation order 2. The purpose of this paper was to determine the simulation of the model 1-D of wastewater transport phenomena based advection-diffusion mechanism and did validate the model. Numerical methods which was used for the solution of this model is finite difference method with Forward Time Central Space scheme. The simulation results which was obtained would be compared with field observation data as a validation model. Collection of field data was carried out in the Wastewater Treatment Plant Sewon, Bantul, D.I. Yogyakarta. The results of numerical simulations were indicate that the advection-diffusion mechanism takes place continuously over time. Then validation of the model was state that there was a difference between the calculation results with the field data, with a correlation value of 0.998.

• Geomechanics and Engineering
• /
• 제18권3호
• /
• pp.225-234
• /
• 2019

A numerical stepwise approach for cavity expansion problem in strain-softening rock or soil mass is investigated, which is compatible with Mohr-Coulomb and generalized Hoek-Brown failure criteria. Based on finite difference method, plastic region is divided into a finite number of concentric rings whose thicknesses are determined internally to satisfy the equilibrium and compatibility equations, the material parameters of the rock or soil mass are assumed to be the same in each ring. For the strain-softening behavior, the strength parameters are assumed to be a linear function of deviatoric plastic strain (${\gamma}p^*$) for each ring. Increments of stress and strain for each ring are calculated with the finite difference method. Assumptions of large-strain for soil mass and small-strain for rock mass are adopted, respectively. A new numerical stepwise approach for limited pressure and plastic radius are obtained. Comparisons are conducted to validate the correctness of the proposed approach with Vesic's solution (1972). The results show that the perfectly elasto-plastic model may underestimate the displacement and stresses in cavity expansion than strain-softening coefficient considered. The results of limit expansion pressure based on the generalised H-B failure criterion are less than those obtained based on the M-C failure criterion.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2007년도 정기총회 및 학술발표대회
• /
• pp.395-398
• /
• 2007

In this study, new dispersion-correction terms are added to leap-frog finite difference scheme for the linear shallow-water equations with the purpose of considering the dispersion effects of the linear Boussinesq equations for the propagation of tsunamis. The new model is applied to near Gadeok island in Pusan about The Central East Sea Tsunami in 1983 and The Hokkaldo Nansei Oki Earthquake Tsunami in 1993 one simulated in the study.

• 한국전산구조공학회논문집
• /
• 제31권6호
• /
• pp.331-337
• /
• 2018

본 논문은 MLS(moving least squares) 차분법의 1차 미분 근사함수를 바탕으로 시간에 따른 수치해석이 가능한 해석기법을 제시한다. 오직 1차 미분 근사함수로만 지배방정식을 이산화했으며, 근사함수를 조립하는 형태로 전체 시스템 방정식을 구성하여 차분법으로 이산화된 운동방정식이 유한요소법(finite element method)과 유사한 모습을 갖게 되었다. 운동방정식을 시간적분하기 위해서 중앙차분법(central difference method)을 사용하였다. 유한요소 알고리즘을 통해서 MLS 차분법과 유한요소법의 고유진동 해석을 수행하였으며, 두 해석결과를 비교하였다. 또한, 동적해석 결과를 기존의 2차 미분 근사함수를 활용한 해석결과와 함께 도시함으로써 제안된 수치기법의 정확성을 검증하였다. 1차 미분 근사함수를 조립하는 과정에서 해석결과의 떨림현상이 억제되었으며 상대적으로 균일한 응력분포를 구할 수 있었다.

• 전산구조공학
• /
• 제9권3호
• /
• pp.187-197
• /
• 1996

Explicit 직접적분법 알고리듬을 사용하여 Euler기둥의 동적 좌굴거동을 해석할 수 있는 수치해석법을 제시하였다. 평면뼈대 유한요소를 기하학적 비선형 거동과 전체좌굴의 영향을 고려할 수 있도록 보의 대변위 이론으로부터 유도하였고, central difference method를 바탕으로 해석 알고리듬을 개발하였다. 다양한 형상, 크기, 재하시간을 갖는 충격하중에 대하여 Euler기둥의 동적좌굴거동과 고유치 문제를 해석하였다. 수치해석 예제를 통하여 본 연구의 결과를 검증하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• pp.1827-1832
• /
• 2003

In this research, the simulation method for acoustic sounds by a uniform flow around a two-dimensional circular cylinder by using the finite difference lattice Boltzmann model is explained. To begin with, we examine the boundary condition which determined with the distribution function $f_i^{(0)}$ concerning with density, velocity and internal energy at boundary node. Very small acoustic pressure fluctuation, with same frequency as that of Karman vortex street, is compared with the pressure fluctuation around a circular cylinder. The acoustic sound' propagation velocity shows that acoustic approa ching the upstream, due to the Doppler effect in the uniform flow, slowly propagated. For the do wnstream, on the other hand, it quickly propagates. It is also apparently the size of sound pressure was proportional to the central distance $r^{-1/2}$ of the circular cylinder. The lattice BGK model for compressible fluids is shown to be one of powerful tool for simulation of gas flows.