• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 추계 학술대회논문집
• /
• pp.47-50
• /
• 2006

In the aeroacoustic application of computational fluid dynamics, the physical phenomena like the crackle in the unsteady compressible jets should be based on very time-accurate numerical solution. The accuracy of the present numerical scheme is extended to the fifth order, using the WENO filter to the sixth-order central difference computation. However, the computational capacity is very restricted by the environment of computational power, so therefore the quadrilateral adaptive grids technique is introduced for this high-order accuracy scheme. The first problem is the multi-dimensional interpolation between fine and coarse grids. Some general benchmark problems are solved to show the effectiveness of this method.

• 한국해안해양공학회지
• /
• 제14권2호
• /
• pp.116-127
• /
• 2002

대부분의 지진해일에 대한 수치모의는 천수방정식을 지배방정식으로 하고 leap-frog 유한차분기법을 주로 적용한다. 이 기법은 격자간격을 적절히 선정하면 분산효과를 옳게 고려할 수 있으나, 수심이 변하는 경우 전 계산영역에서 분산효과를 모두 만족시킬 수는 없다. 본 연구에서는 균일한 격자와 계산시간간격을 사용하면서도 수심이 완만하게 변하는 지형상을 전파하는 지진해일의 국부적인 분산관계를 만족시키기 위해 기존의 2차원 완 변수심상 분산보정 유한차분기법을 개선하고, 이 수치모형의 현장 적용성을 검증하기 위해 1983년 동해중부 지진해일을 검조기록이 있는 동해안의 속초, 묵호, 포항 그리고 울산항에 대하여 수치모의하였다. 또한 1983년 지진해일에 대해 각 항만에서 측정된 검조기록과 수치모의에 의해 계산된 값을 비교 분석하여 수치모형의 정밀도를 평가하였다.

• Structural Engineering and Mechanics
• /
• 제1권1호
• /
• pp.87-102
• /
• 1993

A $C^{\circ}$ continuous finite element formulation of a higher order displacement theory is presented for predicting linear and geometrically non-linear in the sense of von Karman transient responses of composite and sandwich plates. The displacement model accounts for non-linear cubic variation of tangential displacement components through the thickness of the laminate and the theory requires no shear correction coefficients. In the time domain, the explicit central difference integrator is used in conjunction with the special mass matrix diagonalization scheme which conserves the total mass of the element and included effects due to rotary inertia terms. The parametric effects of the time step, finite element mesh, lamination scheme and orthotropy on the linear and geometrically non-linear responses are investigated. Numerical results for central transverse deflection, stresses and stress resultants are presented for square/rectangular composite and sandwich plates under various boundary conditions and loadings and these are compared with the results from other sources. Some new results are also tabulated for future reference.

• 융합보안논문지
• /
• 제5권3호
• /
• pp.93-97
• /
• 2005

MacCormack 방법은 hyperbolic 편미분 방정식의 근을 구하는데 많이 쓰이는 방법으로 그 정확도가 2차 오더가 된다. 하지만 이 방법으로 편미분방정식을 풀 경우 불연속인 점에서는 엔트로피를 만족하지 않는 경우가 있어 우리는 임의의 항을 첨가하여 근을 구해야한다. 이 임의의 항을 첨가하지 않고 직접 방정식으로부터 구하는 방법을 생각하는데 있어서 기존의 MacCormack 방법에 새 central scheme의 개념을 이용하면 전형적인 MacCormack 방법의 정확도와 장점을 보존할 수 있다. 이 새로운 방법을 이용하여 1D Burgers' 방정식과 1D Euler gas dynamic 방정식에 활용하여 그 결과를 살펴본다.

• 대한원격탐사학회:학술대회논문집
• /
• 대한원격탐사학회 2003년도 Proceedings of ACRS 2003 ISRS
• /
• pp.351-353
• /
• 2003

We present a retrieval scheme for the remote sensing of evapotranspiration (ET) over rice paddy. To perform the retrieval, high-resolution airborne imagery of multi-spectral visible and thermal infrared data, and ground-based meteorological measurements are utilized. Our ET retrieval scheme is based on the basic principal of surface energy budget, which is a result of balance in longwave and shortwave radiation, latent heat, sensible heat, and energy flux into the ground. To partition the latent and sensible heat fluxes of interest from the energy balance equation, three basic parameters are of most concern, including albedo, surface temperature, and normalized difference vegetation index (NDVI). The NDVI and albedo can be easily derived from the visible and near infrared spectral data, while the surface tem-perature can be determined through the analysis of the infrared data with the Stefan Boltzmann law. From the airborne imagery taken on 28 April 2003, we observe very good dry and wet pixels that can be easily corre-sponded to the radiation and evaporation controlled crite-ria, respectively, and, hence, for the further use in defin-ing the evaporative fraction needed to partition sensible and latent heat fluxes from the net energy flux. The de-rived ET is compared with the in situ measurements.

• 한국윤활학회:학술대회논문집
• /
• 한국윤활학회 1999년도 제29회 춘계학술대회
• /
• pp.182-188
• /
• 1999

A theoretical study of elastohydrodynamic lubrication of elliptical contacts with both rolling and spin has been carried out. A finite difference method and the Newton-Raphson method are applied to solve the problem. The velocity vectors resulting from combined spin and rolling/sliding motion lead to asymmetric pressure distributions and film shapes. Film contours and variations of the minimum and central film thicknesses are compared with various spin-roll ratios. At high spin-roll ratios the minimum film thickness is considerably reduced, whereas the central film thickness decreases less dramatically, The present numerical scheme can be used in the analysis of general elliptical contact problems.

• Tribology and Lubricants
• /
• 제23권2호
• /
• pp.49-55
• /
• 2007

A numerical analysis of elastohydrodynamic lubrication of elliptical contacts with both rolling and spinning has been carried out. A finite difference method with non-uniform grid systems and the Newton-Raphson method are applied to solve the problems. The velocity vectors resulting from combined spinning and rolling/sliding motion lead to asymmetric pressure distributions and film shapes. Pressure distributions, film contours and variations of the minimum and central film thicknesses are compared with various spin-roll ratios. Reduction of the minimum film thickness under spinning is remarkable whereas the central film thickness is relatively less. The spin motion have large effect on variations of the minimum film thickness with load parameter which are small in pure rolling/sliding cases. Therefore present numerical scheme can be used in the analysis of general elliptical contact EHL problems and further studies are required.

• Tribology and Lubricants
• /
• 제24권6호
• /
• pp.355-361
• /
• 2008

In highly stressed machine elements such as angular contact ball bearings and toroidal type traction drives, the elastohydrodynamic lubrication of elliptical contacts with both rolling and spinning motion are occur. In this paper, a finite difference method with non-uniform grid systems and the Newton-Raphson method are applied to solve the problems. Pressure distributions, film contours and variations of the minimum and central film thicknesses are compared with various ellipticity parameter, dimensionless speed and load parameter. The results showed that the spinning motion has significant influence especially on the film shapes. Reduction of the minimum film thickness under spinning is remarkable whereas the central film thickness is relatively less. Especially variations of the minimum film thicknesses with rolling velocity, load and ellipticity ratio are a great different from those of pure rolling. Therefore present numerical scheme can be used in the analysis of general elliptical contact EHL problems and further studies are required.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1996년도 춘계 학술대회논문집
• /
• pp.99-104
• /
• 1996

The jet pump is being used in many fields for several purposes because of its simple construction and easy operation. The characteristics of the geometrical variables, pressure gradient and velocity distribution of the jet pump are studied using the CFD technique. The flow calculations through a bended nozzle. a mixing chamber and a venturi are presented and phenomenological aspects are discussed. This study solve 3-D steady incompressible Navier-Stokes equations using the Iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives and 3rd-order accurate QUICK scheme for the convective terms. The Mark-and-cell concept was applied efficiently to solve continuity equation, which is differenced 2nd-order accurate central differenced scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A O-type of grid system is generated inside a nozzle and venturi of the jet pump. It has concluded that the results of present study properly agree with physical flow phenomena.