• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.553-556
• /
• 2002

A numerical study on the use of the momentum interpolation mettled for flows with a large body force is presented. The inherent problems of the momentum interpolation method are discussed first. Numerical experiments are performed for a typical flow involving a large body force. The tact that the momentum interpolation method may result in physically unrealistic solutions is demonstrated. Numerical experiments changing the numerical grid have shown that a simple way of removing the physically unrealistic solution is a proper grid refinement where there is a large pressure gradient. An effective way of specifying the pressure and pressure correction at the boundary by a local mass conservation near the boundary is proposed, and it is shown that this method can effectively remove the inherent problem of the specification of pressure and pressure correction at the boundary when one uses the momentum interpolation method.

• 한국전산유체공학회지
• /
• 제7권2호
• /
• pp.8-16
• /
• 2002

A numerical study on the use of the momentum interpolation method for flows with a large body force is presented. The inherent problems of the momentum interpolation method are discussed first. The origins of problems of the momentum interpolation methods are the validity of linear assumptions employed for the evaluation of the cell-face velocities, the enforcement of mass conservation for the cell-centered velocities and the specification of pressure and pressure correction at the boundary. Numerical experiments are performed for a typical flow involving a large body force. The numerical results are compared with those by the staggered grid method. The fact that the momentum interpolation method may result in physically unrealistic solutions is demonstrated. Numerical experiments changing the numerical grid have shown that a simple way of removing the physically unrealistic solution is a proper grid refinement where there is a large pressure gradient. An effective way of specifying the pressure and pressure correction at the boundary by a local mass conservation near the boundary is proposed, and it is shown that this method can effectively remove the inherent problem of the specification of pressure and pressure correction at the boundary when one uses the momentum interpolation method.

• Nuclear Engineering and Technology
• /
• 제32권2호
• /
• pp.169-179
• /
• 2000

A detailed numerical analysis of the evolution of thermal stratification in a curved piping system in a nuclear power plant is performed. A finite volume based thermal-hydraulic computer code has been developed employing a body-fitted, non-orthogonal curvilinear coordinate for this purpose. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of momentum interpolation method. The SIMPLE algorithm is employed for the pressure and velocity coupling, and the convection terms are approximated by a higher-order bounded scheme. The thermal-hydraulic computer code developed in the present study has been applied to the analysis of thermal stratification in a curved duct and some of the predicted results are compared with the available experimental data. It is shown that the predicted results agree fairly well with the experimental measurements and the transient formation of thermal stratification in a curved duct is also well predicted.

• 대한기계학회논문집B
• /
• 제28권6호
• /
• pp.619-627
• /
• 2004

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method. The results agree well with previous numerical and experimental results. This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications. The usefulness of this method will be confirmed when we solve the complex geometries and moving bodies.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.916-921
• /
• 2003

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method(FVM). This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications.

• 한국항만학회지
• /
• 제13권1호
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• Journal of Mechanical Science and Technology
• /
• 제14권7호
• /
• pp.789-795
• /
• 2000

Free-surface flows with an arbitrary deformation, induced by a submerged hydrofoil, are simulated numerically, considering two-fluid flows of both water and air. The computation is performed by a finite volume method using unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell wise local mesh refinement. The integration in space is of second order, based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels. The linear equations are solved by conjugate gradient type solvers, and the non-linearity of equations is accounted for through Picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations, the continuity equation, the conservation equation of one species, and the equations for two turbulence quantities. Finally, a comparison is quantitatively made at the same speed between the computation and experiment in which the grid sensitivity is numerically checked.

• Journal of Advanced Marine Engineering and Technology
• /
• 제32권8호
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• 대한기계학회논문집B
• /
• 제23권2호
• /
• pp.243-253
• /
• 1999

The dual reciprocity boundary element method (DRBEM) is used to solve the Graetz problem of laminar flow inside circular duct. In this method the domain integral tenn of boundary integral equation resulting from source term of governing equation is transformed into equivalent boundary-only integrals by using the radial basis interpolation function, and therefore complicate domain discretization procedure Is completely removed. Velocity profile is obtained by solving the momentum equation first and then, using this velocities as Input data, energy equation Is solved to get the temperature profile by advancing from duct entrance through the axial direction marching scheme. DRBEM solution is tested for the uniform temperature and heat flux boundary condition cases. Local Nusselt number, mixed mean temperature and temperature profile inside duct at each dimensionless axial location are obtained and compared with exact solutions for the accuracy test Solutions arc in good agreement at the entry region as well as fully developed region of circular duct, and their accuracy are verified from error analysis.

• 대한조선학회논문집
• /
• 제29권3호
• /
• pp.59-70
• /
• 1992

2차원 비압축성, 정상, 층류 Navier-Stokes 방정식을 일반 곡선 좌표계에서 계산하기 위한 유한체적법을 개발하였다. 수치해석은 정규 격자구조를 채택하였으며, 이때 압력의 불안정한 거동은 모멘텀 보간법에 의하여 제거하였다. NACA0012 날개 단면 주위의 유동을 개발된 컴퓨터 프로그램으로 계산한 결과들은 실험 및 다른 계산결과와 잘 일치하였다.