• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.2
• /
• pp.276-284
• /
• 2003

The preconditioned conjugate gradient method is an efficient iterative solution scheme for large size finite element problems. As preconditioning method, we choose an incomplete Cholesky factorization which has efficiency and easiness in implementation in this paper. The incomplete Cholesky factorization mettled sometimes leads to breakdown of the computational procedure that means pivots in the matrix become minus during factorization. So, it is inevitable that a reduction process fur stabilizing and this process will guarantee robustness of the algorithm at the cost of a little computation. Recently incomplete factorization that enhances robustness through increasing diagonal dominancy instead of reduction process has been developed. This method has better efficiency for the problem that has rotational degree of freedom but is sensitive to parameters and the breakdown can be occurred occasionally. Therefore, this paper presents new method that guarantees robustness for this method. Numerical experiment shows that the present method guarantees robustness without further efficiency loss.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.20-30
• /
• 2001

A multigrid DADI method for low Mach number preconditioning method is presented. The eigenvalues of governing equations are modified by A low Mach number preconditioner developed by Choi & Merkle, and it results in an accurate solution and fast convergence In the low Mach number region. The convergence of numerical method is further accelerated by multigrid method. The efficient and accuracy of present method is shown by comparison with conventional solution method for the compressible flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.181-182
• /
• 2003

A preconditioned numerical method for gas-liquid to-phase flow is applied to solve cavitating flow. The present method employs a density based finite-difference method of dual time-stepping integration procedure and Roe's flux difference splitting approximation with MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. The method permits simple treatment of the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristics at low Mach number. By this method, two-dimensional internal flows through a venturi tuve and decelerating cascades are computed and discussed.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1903-1908
• /
• 2004

A preconditioned numerical method for gas-liquid two-phase flows is applied to solve cavitating flow. The present method employs a finite-difference dual time-stepping integration procedure and the MUSCLTVD scheme. A homogeneous equilibrium cavitation model is used. The present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation, large density changes and incompressible flow characteristics at low Mach number. Some internal flows such as convergent-divergent nozzles are computed using this method. Comparisons of predicted and experimental results are provided and discussed.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.5 s.54
• /
• pp.11-18
• /
• 2003

In two-dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure (HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint hence, achieving an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared (CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper and efficient iterative procedure for higher order finite element formulation has not been available, thus far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient, and useful in flow analyses, in conjunction with hierarchical elements.

• Communications of the Korean Mathematical Society
• /
• v.11 no.3
• /
• pp.867-880
• /
• 1996

Numerical approximations to the linear elasticity are traditionally based on the finite element method. In this paper we propose a new formulation based on the cubic spline collocation method for linear elastic problem on the unit square. We present several numerical results for the eigenvalues of the matrix represented by cubic collocation method and preconditioner matrix which is preconditioned by FEM and FDM. Finally we present the numerical solution for some example equation.

• Bulletin of the Korean Mathematical Society
• /
• v.38 no.1
• /
• pp.17-27
• /
• 2001

We discuss a finite difference preconditioner for the$C^1$ Lagrance quadratic spline collocation method for a uniformly elliptic operator with homogeneous Dirichlet boundary conditions. Using the generalized field of values argument, we analyzed eigenvalues of the matrix preconditioned by the matrix corresponding to a finite difference operator with zero boundary condition.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.216-221
• /
• 2004

In two dimensional incompressible flaws, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint and hence achieve an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flaw analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flaw showed the present procedure to be stable, very efficient and useful in flaw analyses in conjunction with hierarchical elements.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.155-162
• /
• 2000

In this paper, preconditioning method is parallelized on fast-ethernet PC cluster. The algorithm is based on scaling the pressure terms in the momemtum equations and preconditioning the conservation equations to circumvent numerical difficulties at low Mach numbers. Parallelization is performed using a domain decomposition technique(DDT) and message passing between sub-domains are taken from the MPI library. The results are shown to have good convergence properties at all Mach number on the circular arc Bump and are capable of reasonable predicting two-dimensional turbulent flows on DCA compressor cascade.