• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.105-105
• /
• 2021

본 연구에서는 2차원 비선형 천수모형과 수심평균된 스칼라 이송모형을 해석하는 수치모형에 대해 기술하였다. 수치모형의 정확성을 보장함과 동시에 안정성을 높이기 위해 유한체적법, 플럭스 재구성 및 minmod 제한자를 사용하였다. 비선형 천수방정식의 이송항과 바닥 경사항은 계산된 수심의 양수 보존과 흐름의 정상 상태를 보장하기 위한 second order well-balanced positivity preserving central-upwind method를 이용하여 수치적으로 이산화되었다. 마찬가지로, 이송-확산 방정식 내 이송항은 동일한 2차 풍상차분법을 통해 수치적으로 풀이하였다. 격자점 경계면에서의 불연속으로 인한 수치진동을 방지하기 위해 이송항의 계산에 포함된 보존항의 차이로 인해 발생하는 스칼라의 수치확산을 최소화하기 위해 무차원의 비소산함수를 도입하였다. 또한, 확산항은 유한차분법을 이용하여 이산화하였다. 제안된 수치모형은 시간미분항의 계산을 위해 오일러 기법을 적용하여 계산된 수심 및 스칼라의 양수 보존여부와 함께 정지된 흐름의 정상 상태의 보존여부를 확인하였다. 제안된 수치모형의 해석 정확성을 평가하기 위해 1, 2차원 공간 내 다양한 흐름 조건에서의 해석해를 이용한 3개의 벤치마크 테스트를 수행하였다. 평균 제곱근 오차(Root Mean Squared Error, RMSE)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.4A
• /
• pp.693-700
• /
• 2006

To detect and evaluate the damage present in bridge, static identification method is known to be simple and effective, compared to dynamic method. In this study, the damage detection method in steel box girder bridge using static responses including displacement, slope and curvature is examined. The static displacement is calculated using finite element analysis and the slope and curvature are determined from the displacement using central difference method. The location of damage is detected using the absolute differences of these responses in intact and damaged bridge. Steel box girder bridge with corner crack is modeled using singular element in finite element method. The results show that these responses were significantly useful in detecting and predicting the location of damage present in bridge.

• Journal of Korean Ophthalmic Optics Society
• /
• v.15 no.2
• /
• pp.151-154
• /
• 2010

Purpose: We investigated how the movement of iris and visual axis affects the finite schematic eye Methods: Using the schematic eye with the crystalline lens in the existing forms of the radial GRIN and the spherical GRIN, the iris centre was moved 0.5 mm in nasal direction and visual axis was tilted $5^{\circ}$ in same direction, with the additional degree of 2.5 down to locate the focal point in fovea. This study analyzed performance change of the optical system, designing it same as the real eye. Results: The whole aberration distribution showed a considerable difference in performance in comparison with the real eye; the biggest difference shown at the central field of optical system. The spherical aberration showed the biggest difference, and a peripheral power error and field curvature leaned toward (+) direction in aberration distribution. Conclusions: When designing the schematic eye with the performance similar with that of the real eye by taking into consideration the iris centre and visual axis, the aberration at the center field of optical system in particular should be corrected. Spherical aberration which showed the biggest difference should be corrected in the first place. In addition, a peripheral power error and field curvature that leaned toward (+) direction should be moved toward (-) direction.

• Bulletin of the Korean Mathematical Society
• /
• v.55 no.3
• /
• pp.881-898
• /
• 2018

This article deals with implementation of a high-order finite difference scheme for numerical solution of the incompressible Navier-Stokes equations on curvilinear grids. The numerical scheme is based on pseudo-compressibility approach. A fifth-order upwind compact scheme is used to approximate the inviscid fluxes while the discretization of metric and viscous terms is accomplished using sixth-order central compact scheme. An implicit Euler method is used for discretization of the pseudo-time derivative to obtain the steady-state solution. The resulting block tridiagonal matrix system is solved by approximate factorization based alternating direction implicit scheme (AF-ADI) which consists of an alternate sweep in each direction for every pseudo-time step. The convergence and efficiency of the method are evaluated by solving some 2D benchmark problems. Finally, computed results are compared with numerical results in the literature and a good agreement is observed.

• Tribology and Lubricants
• /
• v.14 no.4
• /
• pp.58-63
• /
• 1998

This paper presents a study of the elastohydrodynamic lubrication of elliptical contacts where lubricant entrainment coincides with the major axis of the Hertzian contact ellipse. A finite difference method and the Newton-Raphson method are applied to analyze the problem. Film contours and pressure distributions are compared with the results for lubricant entrainment coincides with the minor axis. Variations of the minimum and central film thicknesses with the radius ratio are also examined. The results showed that the present numerical scheme can be used generally in the analysis of the EHL of elliptical contacts where the lubricant entraining vector did not coincide with either of the principal axis of the conjunction.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1998.04a
• /
• pp.81-87
• /
• 1998

In this paper, a global acoustic design sensitivity analysis (DSA) of field point pressure with respect to structural sizing design variables is developed. Firstly acoustic sensitivity is formulated and implemented numerically. And it is combined with continuum structural sensitivity to obtain the global acoustic, design sensitivity. For this procedure, GASA (global acoustic design sensitivity analyzer) has been developed. A half scale of automobile cavity model is considered in this paper. In order to confirm accuracy of the results of global acoustic DSA obtained by GASA, it is compared with the result of central finite difference method. In order to reduce computation time, Rayleigh approximated solution is evaluated and compared with the solution which used every nodal velocities. Also the acoustic optimization procedure is performed using design sensitivities. From these numerical studies, it can be shown that global acoustic DSA is a useful tool to improve acoustic problems.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.746-752
• /
• 2001

Design sensitivity analysis for nonlinear structural problems has been emerged in the last decade as a glowing area of engineering research. As a result, theoretical formulations and computational algorithms have already developed for design sensitivity of nonlinear structural problems. There is not enough research for practical nonlinear problems using multi-element, due to difficulties of implementation into FEA. Therefore, nonlinear response analysis for stiffened shell which consists of Mindlin plate and Timoshenko beam, was considered. Specially, it presents the backward-Euler method which is adopted to describe an exact yield state in the stress computation procedure. Then, design sensitivity analysis of nonlinear structures, particularly elasto-perfectly-plastic structure, is developed using direct differentiation method. The accuracy of the developed sensitivity analysis was compared with the central finite difference method. Finally, on the basis of above results, design improvement for stiffened shell is suggested.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1998.10a
• /
• pp.212-218
• /
• 1998

This paper presents a study of the elastohydrodynamic lubrication of elliptical contacts where lubricant entrainment coincides with the major axis of the Hertzian contact ellipse. A finite difference method and the Newton-Raphson method are applied to analyze the problem. Film contours and pressure distributions are compared with the results for lubricant entrainment coincides with the minor axis. Variations of the minimum and central film thicknesses with the radius ratio are also examined. Therefore, the present numerical scheme can be used generally in the analysis of the EHL of elliptical contacts where the lubricant entraining vector did not coincide with either of the principal axis of the conjunction.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.28 no.1
• /
• pp.79-92
• /
• 1992

This is analyzed using the finite element method which is appling excellent isoparametric curve element in the aspect of large usages of dynamic responses in which is regarding geometric and material nonlinear of a large scale shell structure of an airplane, a submarine, a ship, and an ocean structure. The solution of dynamic equations is got by direct integration method using time-stepping procedure and regarding Central Difference Method of the both solutions. But because formal matrix factorization is not necessary in each time step and it does not take less time to compute relatively, this method must be regarded very few time steps on the condition. Axisymmatric shell problems are inspected using 8 node Isoparametric element in this paper. Partial axisymmatric spherical shell is used as a model to analyze axisymmatric nonlinear dynamic behavior regarding. Total Lagrangian formulation in geometric nonlinear behavior and elastio-viscoplastic in material nonlinear behavior.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.8 s.251
• /
• pp.796-805
• /
• 2006

In this study, the problem of using the mixing index as a measure of the mixing performance for a certain flow field has been discussed. The flow model subjected to this study is the two-dimensional unsteady lid-driven cavity flow. The transport equation for the concentration within the cavity was solved by using the finite volume method where the convective terms are discretized with the central difference scheme. It was shown that both the concentration dispersion and the mixing index depend highly on the initial distribution of the concentration, and therefore the mixing index obtained from the concentration dispersion equation loses its universal applicability.