• Journal of applied mathematics & informatics
• /
• 제16권1_2호
• /
• pp.19-35
• /
• 2004

We develop in this paper some preconditioners for sparse non-symmetric M-matrices, which combine a recursive two-level block I LU factorization with multigrid method, we compare these preconditioners on matrices arising from discretized convection-diffusion equations using up-wind finite difference schemes and multigrid orderings, some comparison theorems and experiment results are demonstrated.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1778-1782
• /
• 2004

A new algorithm, CBIR (Correlation-Based Image Registration) was proposed to improve the resolution of image registration for PSP (Pressure-Sensitive Paint). The local displacement vectors were obtained by finding the displacement which maximizes the cross-correlation between two interrogation windows of 'wind-off' and 'wind-on' images. A recursive multigrid processing was employed to increase the non-linear spatial resolutions. The variations of image were precisely measured without identifying the control points.

• 한국윤활학회:학술대회논문집
• /
• 한국윤활학회 1999년도 제30회 추계학술대회
• /
• pp.323-330
• /
• 1999

The conventional analysis for the numerical computation of fluid film thickness with elastic deformation of contact region. is performed by Newton-Rephson method for its 18st convergence characteristics. However, both high load and relatively low sliding velocity frequently make it impossible for Newton-Rahpson method to get both converged and stable solutions. In particular, this method cannot provide converged Solution under the condition of high load above 1.0 GPa which frequently occurs in line contact of EHL problem. Multigird multi-level method for the solver of non-linear partial differential equation including solid deformation is preferred to Newton-Rshpson method for better convergence and stability and is applied to line contact EHL behavior in this study.

• 대한수학회지
• /
• 제46권5호
• /
• pp.1007-1025
• /
• 2009

First-order least-squares method of a distributed optimal control problem for the incompressible Stokes equations is considered. An optimality system for the optimal solution are reformulated to the equivalent first-order system by introducing the vorticity and then the least-squares functional corresponding to the system is defined in terms of the sum of the squared $H^{-1}$ and $L^2$ norms of the residual equations of the system. Finite element approximations are studied and optimal error estimates are obtained. Resulting linear system of the optimality system is symmetric and positive definite. The V-cycle multigrid method is applied to the system to test computational efficiency.

• 한국수학교육학회지시리즈B:순수및응용수학
• /
• 제21권2호
• /
• pp.129-139
• /
• 2014

In this paper, we present a detailed comparison of the performance of the numerical solvers such as the biconjugate gradient stabilized, operator splitting, and multigrid methods for solving the two-dimensional Black-Scholes equation. The equation is discretized by the finite difference method. The computational results demonstrate that the operator splitting method is fastest among these solvers with the same level of accuracy.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 12th Asian Congress of Fluid Mechanics
• /
• pp.33.3-33.3
• /
• 2008

• Tribology and Lubricants
• /
• 제19권3호
• /
• pp.123-132
• /
• 2003

The effect of surface texture on elastohydrodynamic lubrication (EHL) point contact of a ball Joint mechanism in small reciprocating compressors is studied numerically by using multigrid method. Pressure and film thickness profiles have been calculated for surface roughness with waviness of different orientations and transverse ridge and dent at minimum and maximum Hoes M parameter conditions. The influence of the amplitude and the wavelength of the surface roughness was also studied. Results show that the oblique waviness with orientation angle of 30$^{\circ}$generates the smallest minimum film thickness as compared with those of longitudinal, transverse, and other oblique roughness. The influence of transverse waviness on the minimum film thickness is smaller than for the longitudinal waviness case.

• 대한수학회논문집
• /
• 제20권3호
• /
• pp.563-578
• /
• 2005

In [Z. Cai and B. C. Shin, SIAM J. Numer. Anal. 40 (2002), 307-318], we developed the discrete first-order system least squares method for the second-order elliptic boundary value problem by directly approximating $H(div){\cap}H(curl)-type$ space based on the Helmholtz decomposition. Under general assumptions, error estimates were established in the $L^2\;and\;H^1$ norms for the vector and scalar variables, respectively. Such error estimates are optimal with respect to the required regularity of the solution. In this paper, we study solution methods for solving the system of linear equations arising from the discretization of variational formulation which possesses discrete biharmonic term and focus on numerical results including the performances of multigrid preconditioners and the finite element accuracy.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• 제26권1호
• /
• pp.1-14
• /
• 1998

Various numerical methods for the two dimensional shallow water equations have been applied to the problems of flood routing, tidal circulation, storm surges, and atmospheric circulation. These methods are often based on the Alternating Direction Implicity(ADI) method. However, the ADI method results in inaccuracies for large time steps when dealing with a complex geometry or bathymetry. Since this method reduces the performance considerably, a fully implicit method developed by Wilders et al. (1998) is used to improve the accuracy for a large time step. Finite Difference Methods are defined on a rectangular grid. Two drawbacks of this type of grid are that grid refinement is not possibile locally and that the physical boundary is sometimes poorly represented by the numerical model boundary. Because of the second deficiency several purely numerical boundary effects can be involved. A boundary fitted curvilinear coordinate transformation is used to reduce these difficulties. It the curvilinear coordinate transformation is used to reduce these difficulties. If the coordinate transformation is orthogonal then the transformed shallow water equations are similar to the original equations. Therefore, an orthogonal coorinate transformation is used for defining coordinate system. A multigrid (MG) method is widely used to accelerate the convergence in the numerical methods. In this study, a technique using a MG method is proposed to reduce the computing time and to improve the accuracy for the orthogonal to reduce the computing time and to improve the accuracy for the orthogonal grid generation and the solutions of the shallow water equations.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제18권1호
• /
• pp.27-41
• /
• 2014

In this paper, we present an efficient numerical method for multiphase image segmentation using a multiphase-field model. The method combines the vector-valued Allen-Cahn phase-field equation with initial data fitting terms containing prescribed interface width and fidelity constants. An efficient numerical solution is achieved using the recently developed hybrid operator splitting method for the vector-valued Allen-Cahn phase-field equation. We split the modified vector-valued Allen-Cahn equation into a nonlinear equation and a linear diffusion equation with a source term. The linear diffusion equation is discretized using an implicit scheme and the resulting implicit discrete system of equations is solved by a multigrid method. The nonlinear equation is solved semi-analytically using a closed-form solution. And by treating the source term of the linear diffusion equation explicitly, we solve the modified vector-valued Allen-Cahn equation in a decoupled way. By decoupling the governing equation, we can speed up the segmentation process with multiple phases. We perform some characteristic numerical experiments for multiphase image segmentation.