• 한국전산유체공학회지
• /
• 제15권2호
• /
• pp.55-61
• /
• 2010

A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

• Korea-Australia Rheology Journal
• /
• 제21권1호
• /
• pp.47-58
• /
• 2009

We have investigated the transient behavior of 1D fully developed Poiseuille viscoelastic flow under finite pressure gradient described by the Oldroyd-B and Leonov constitutive equations. For analysis we employ a simple $2^{nd}$ order discretization scheme such as central difference for space and the Crank-Nicolson for time approximation. For the analysis of the Oldroyd-B model, we also apply the analytical solution, which is obtained again in this work in terms of elementary solution procedure simpler than the previous one (Waters and King, 1970). Both models demonstrate qualitatively similar solutions, but their eventual steady flowrate exhibits noticeable difference due to the absence or presence of shear thinning behavior. In the inertialess flow, the flowrate instantaneously attains a large value corresponding to the Newtonian creeping flow and then decreases to its steady value when the applied pressure gradient is low. However with finite liquid density the flow field shows severe fluctuation even accompanying reversals of flow directions. As the assigned pressure gradient increases, the flowrate achieves its steady value significantly higher than its value during oscillations after quite long period of time. We have also illustrated comparison between 1D and 2D results and possible mechanism of complex 2D flow rearrangement employing a previous solution of [mite element computation. In addition, we discuss some mathematical points regarding missing boundary conditions in 2D modeling due to the change of the type of differential equations when varying from inertialess to inertial flow.

• 물과 미래
• /
• 제25권1호
• /
• pp.101-110
• /
• 1992

유역내 공간적 및 시간적 분포특성을 가진 이동강우를 해석하기 위하여 분포형 모델을 개발하였으며 이 유역모델은 지표면 흐름과 하천망 흐름으로 나누어 유출을 해석하였다. 지표면의 흐름은 2차원의 연속방정식과 운동량 방정식을 이요하였는데 kinematic 이론을 도입하여 운동량 방정식을 간략화 하였으며 하천망 흐름은 하천망을 일체로 하여 1차원의 연속방정식과 운동량방정식들을 이용하였다. 기본방정식들의 수치해석법으로 지표면의 흐름은 유한요소법을 이용하였으며 하천망에 대한 해석은 음해법의 유한차분법을 이용하였다. 모델은 특히 이동강우에 있어서 중요한 특색인 공간적 및 시간적 특성을 효과적으로 해석할 수 있도록 개발되었다. 또한 모델은 구성된 행렬의 특징을 이용하였는데, 지표면 유출모델은 Gauss 소거법을 이용하여 그해를 구하였으며 하천망 해석은 double sweep 방법을 적용하기 위한 여러 종류의 순환계수방정식을 제안하였고 이를 이용하여 그 해를 구하였다.

• Structural Engineering and Mechanics
• /
• 제17권1호
• /
• pp.1-14
• /
• 2004

Numerical solution to linear bending analysis of circular plates is obtained by the method of harmonic differential quadrature (HDQ). In the method of differential quadrature (DQ), partial space derivatives of a function appearing in a differential equation are approximated by means of a polynomial expressed as the weighted linear sum of the function values at a preselected grid of discrete points. The method of HDQ that was used in the paper proposes a very simple algebraic formula to determine the weighting coefficients required by differential quadrature approximation without restricting the choice of mesh grids. Applying this concept to the governing differential equation of circular plate gives a set of linear simultaneous equations. Bending moments, stresses values in radial and tangential directions and vertical deflections are found for two different types of load. In the present study, the axisymmetric bending behavior is considered. Both the clamped and the simply supported edges are considered as boundary conditions. The obtained results are compared with existing solutions available from analytical and other numerical results such as finite elements and finite differences methods. A comparison between the HDQ results and the finite difference solutions for one example plate problem is also made. The method presented gives accurate results and is computationally efficient.

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

In sensitivity analysis, semi-analytical method(SAM) reveals severe inaccuracy problem when relatively large rigid body motions are identified for individual elements. Recently such errors of SAM resulted by the finite difference scheme have been improved by the separation of rigid body mode. But the eigenvalue should be obtained first before the sensitivity analysis is performed and it takes much time in the case that large system is considered. In the present study, by constructing a reduced one from the original system, iterative method combined with mode decomposition technique is proposed to compute reliable semi-analytical design sensitivities. The sensitivity analysis is performed by the eigenvector acquired from the reduced system. The error of SAM caused by difference scheme is alleviated by Von Neumann series approximation.

• Computers and Concrete
• /
• 제16권5호
• /
• pp.669-682
• /
• 2015

Degradation of reinforced concrete (RC) structures due to chloride penetration followed by reinforcement corrosion has been a serious problem in civil engineering for many years. The numerical simulation methods at present are mainly finite element method (FEM) and finite difference method (FDM), which are based on mesh. Mesh generation in engineering takes a long time. In the present article, the numerical solution of chloride transport in concrete is analyzed using radial point interpolation method (RPIM) and element-free Galerkin (EFG). They are all meshless methods. RPIM utilizes radial polynomial basis, whereas EFG uses the moving least-square approximation. A Galerkin weak form on global is used to attain the discrete equation, and four different numerical examples are presented. MQ function and appropriate parameters have been proposed in RPIM. Numerical simulation results are compared with those obtained from the finite element method (FEM) and analytical solutions. Two case of chloride transport in full saturated and unsaturated concrete are analyzed to test the practical applicability and performance of the RPIM and EFG. A good agreement is obtained among RPIM, EFG, and the experimental data. It indicates that RPIM and EFG are reliable meshless methods for prediction of chloride concentration in concrete structures.

• Korea-Australia Rheology Journal
• /
• 제16권4호
• /
• pp.227-233
• /
• 2004

Simplified model of slide-fed curtain coating flow has been developed and tested in this study. It rests on the sheet profile equations for curtain thickness in curtain flow and its trajectory derived by the integral momentum balance approach of Higgins and Scriven (1979) and Kistler (1983). It also draws on the film profile equation of film thickness variation in flow down a slide. The equations have been solved in finite difference approximation by Newton iteration with continuation. The results show that how inertia (Rey­nolds number), surface tension (capillary number), inclination angle of the slide, and air pressure difference across the curtain affect sheet trajectory and thickness profile. It has been revealed that approximate models can be useful to easily analyze coating flow dynamics without complex computations, giving qualitative agreement with full theory and with experiment.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집A
• /
• pp.686-691
• /
• 2001

Structural optimization often require the evaluation of design sensitivities. The Semi Analytic method(SAM) is popular for shape optimization because this method has several advantages. But when relatively large rigid body motions are identified for individual elements, the SA method shows severe inaccuracy. In this paper, the improvement of design sensitivities corresponding to the rigid body mode is evaluated by exact differentiation of the rigid body modes. Moreover, the error of the SA method caused by numerical difference scheme is alleviated by using a series approximation for the sensitivity derivatives and considering the higher order terms.

• 한국해안해양공학회지
• /
• 제6권3호
• /
• pp.216-225
• /
• 1994

Berkhoff등(1982)의 수리실험 결과를 사용하여 PCGM파 포물형근사식(PA) 수치모형의 정밀도와 계산효율을 비교하였다. 두 모형 결과는 수리실험의 실측치와 잘 부합되었으나 비교적 도고가 큰 천퇴 뒷편의 일부 단면에서는 PCGM모형 결과에 비해 PA모형 결과가 실험치에 가까운 것으로 나타났다. PA모형은 천해역의 넓은 지역에 대한 파랑변형을 예측하는 데 유용한 모형이나 반사파를 무시할 수 있는 경우에만 적용 가능하다. 한편 완경사 파랑식으로부터 직접 차분화하여 반사파를 고려할 수 있는 PCGM모형의 계산효율을 증대시킬 필요가 있다. PA모형 결과를 PCGM모형의 초기 입력조건으로 사용함으로써 약 40%의 계산시간 감소를 가져왔다.

• 대한기계학회논문집
• /
• 제14권6호
• /
• pp.1509-1514
• /
• 1990

본 연구에서는 비압축성 유체가 온도의 영향을 받아 변화하는 경우에 대한 시 일링 간극내에서의 압력분포를 유한차분법을 이용하여 해석하였다. 여러기서 얻어진 결과를 이용하여 시일의 축력과 모멘트를 해석함으로써 시일의 동적 불안정성에 대하 여 논하였다. 이 때 기계평면시일의 형상은 코닝이 있고, 시일의 중심축이 경사진 경우를 고려하였다.