• Journal of the Korean Mathematical Society
• /
• v.61 no.4
• /
• pp.659-681
• /
• 2024

We propose V-cycle multigrid methods for vector field problems arising from the lowest order hexahedral Nédélec finite element. Since the conventional scalar smoothing techniques do not work well for the problems, a new type of smoothing method is necessary. We introduce new smoothers based on substructuring with nonoverlapping domain decomposition methods. We provide the convergence analysis and numerical experiments that support our theory.

• Communications of the Korean Mathematical Society
• /
• v.12 no.3
• /
• pp.813-827
• /
• 1997

Multigrid methods for two first-order system least squares (FOSLS) using bilinear finite elements are developed for the pure traction problem of planar linear elasticity. They are two-stage algorithms that first solve for the gradients of displacement, then for the displacement itself. In this paper, concentration is given on solving for the gradients of displacement only. Numerical results show that the convergences are uniform even as the material becomes nearly incompressible. Computations for convergence rates are included.

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

In this paper, the preconditioned multistage time stepping methods which are popular multigrid smoothers is studied for the compressible flow calculations. Fourier analysis on the local time stepping and block-Jacobi preconditioned residual operators is performed using the linearized 2-D Navier-Stokes equations. It fumed out that block-Jacobi preconditioner has better performance in eigenvalue clustering. They are implemented in the 2-D compressible Euler and Wavier-Stokes calculations with multigrid methods to verify that the block-Jacobi preconditioned multistage time stepping shows better performance in convergence acceleration.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.4 no.2
• /
• pp.135-142
• /
• 2000

We consider a multigrid algorithm for higher order finite difference scheme for the Poisson problem on rectangular domain. Several smoothers including Jacobi, Red-black Gauss-Seidel are tested and compared. Since higher order scheme gives much more accurate result then 5-point scheme, one may use small number of levels with higher order scheme and thus the overall cost is reduced quite a lot. The numerical experiment compares the two cases.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.59-64
• /
• 2002

In this paper, the preconditioned multistage time stepping methods which are popular multigrid smoothers is implemented for the compressible Navier-Stokes calculation with full-coarsening multigrid method. The convergence characteristic of the point-Jacobi and Alternating direction line Jacobi(DDADI) preconditioners are studied. The performance of 2nd order upwind numerical fluxes such as 2nd order upwind TVD scheme and MUSCL-type linear reconstruction scheme are compared in the inviscid and viscous turbulent flow caculations.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.26 no.4
• /
• pp.323-332
• /
• 2022

This paper presents a numerical study on multigrid algorithms of V-cycle type for problems posed in the Hilbert space H(curl) in three dimensions. The multigrid methods are designed for discrete problems originated from the discretization using the hexahedral Nédélec edge element of the lowest-order. Our suggested methods are associated with smoothers constructed by substructuring based on domain decomposition methods of nonoverlapping type. Numerical experiments to demonstrate the robustness and the effectiveness of the suggested algorithms are also provided.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1392-1396
• /
• 2007

The governing equations in generalized curvilinear coordinates for 3D laminar flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms. and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. Multigrid methods are also applied because solving the equations on the coarse grids requires much less computational effort per iteration than on the fine grid. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Humphrey et al., 1977).

• Tribology and Lubricants
• /
• v.25 no.6
• /
• pp.367-373
• /
• 2009

Many research of elastohydrodynamic lubrication (EHL) has been performed under the condition of steady state loading. However, mechanical elements undergo severe high loads that are in the fluctuating modes of frequency and amplitude. Conventional numerical method for the circular contact of EHL study has the difficulty in making the film thickness and pressure of EHL converged in high loads of steady state as well as fluctuating loading conditions. In this work, multigrid multilevel method are used for the stable convergence of film pressure and thickness under the conditions of high as well as varying loads, and very stable solutions of film behaviors with elastic deformation are obtained. Several results of dynamic loading condition are shown and compared with those of steady state condition in the aspects of circular EHL film thickness and pressure.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.6 no.2
• /
• pp.9-24
• /
• 2002

Total Variation(TV) regularization method is effective for reconstructing "blocky", discontinuous images from contaminated image with noise. But TV is represented by highly nonlinear integro-differential equation that is hard to solve. There have been much effort to obtain stable and fast methods. C. Vogel introduced "the Fixed Point Lagged Diffusivity Iteration", which solves the nonlinear equation by linearizing. In this paper, we apply multigrid(MG) method for cell centered finite difference (CCFD) to solve system arise at each step of this fixed point iteration. In numerical simulation, we test various images varying noises and regularization parameter $\alpha$ and smoothness $\beta$ which appear in TV method. Numerical tests show that the parameter ${\beta}$ does not affect the solution if it is sufficiently small. We compute optimal $\alpha$ that minimizes the error with respect to $L^2$ norm and $H^1$ norm and compare reconstructed images.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2001.06a
• /
• pp.218-223
• /
• 2001

Many research of elastohydrodynamic lubrication (EHL) has been performed under the condition of steady state loding. However, mechanical elements undergo severe high loads that are fluctuating in general. Conventional numerical method for the circular contact of EHL study has a difficulty in converging the film pressure and thickness especially in high load of steady state. In this work, multigrid multilevel method expels the convergence problem under the condition of high load and very stable convergence is obtained under the dynamic loading condition over 1.0GPa. Several results of dynamic loading condition are shown and compared with those of steady state condition.