• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.666-670
• /
• 2008

The low-speed preconditioning technique is applied to solve the compressible Reynolds averaged Navier-Stokes equations for low-speed flows. The space discretization is based on Roe's flux-difference splitting with third-order-accurate MUSCL extrapolation. Time integration is performed employing a diagonal approximate factorization algorithm. The dual-time stepping has been incorporated to solve the unsteady flows. Full multigrid method is implemented to accelerate the convergence rate. To verify the algorithms several cases have been tested. Demonstrated the improvement on convergence and quality of the solution.

• Journal of applied mathematics & informatics
• /
• v.30 no.1_2
• /
• pp.235-244
• /
• 2012

In this paper, we solve the three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind using the full multigrid (FMG) algorithm. We derive a finite difference approximations for the biharmonic equation on a 18 point compact stencil. The unknown solution and its second derivatives are carried as unknowns at grid points. In the multigrid methods, we use a fourth order interpolation to producing a new intermediate unknown functions values on a finer grid, and the full weighting restriction operators to calculating the residuals at coarse grid points. A set of test problems gives excellent results.

• 한국가시화정보학회:학술대회논문집
• /
• 2003.11a
• /
• pp.63-66
• /
• 2003

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.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• 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.

• 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).

• Journal of Korea Multimedia Society
• /
• v.7 no.12
• /
• pp.1680-1691
• /
• 2004

One of the representative methods of optical flow is a gradient method which estimates the movement of an object based on the differential of image brightness. However, the method is ineffective for large displacement of the object and many improved methods have been proposed to copy with such limitations. One of these improved techniques is the multigrid processing, which is used in many optical flow algorithms. As an alternative novel technique we have been proposing an orthogonal functional expansion method, where whole displacements are expanded from low frequency terms. This method is expected to be applicable to flow estimation with large displacement and deformation including expansion and contraction, which are difficult to cope with by conventional optical flow methods. In the orthogonal functional expansion method, the apparent displacement field is calculated iteratively by a projection method which utilizes derivatives of the invariant constraint equations of brightness constancy. One feature of this method is that differentiation of the input image is not necessary, thereby reducing sensitivity to noise. In this paper, we apply our method to several real images in which the objects undergo large displacement and/or deformation including expansion. We demonstrate the effectiveness of the orthogonal functional expansion method by comparing with conventional methods including our optimally scaled multigrid optical flow algorithm.

• Wind and Structures
• /
• v.12 no.1
• /
• pp.1-19
• /
• 2009

In this paper, the airflow around an ideal thin plate (hereafter referred to as ITP) with various ratios of central slot is simulated by using the finite-difference-method (FDM)-based Arbitrary-Lagrangian-Eulerian descriptions for the rigid oscillating body. The numerical procedure employs the second-order projection scheme to decouple the governing equations, and the multigrid algorithm with three levels to improve the computational efficiency in evaluating of the pressure equation. The present CFD method is validated through comparing the computed flutter derivatives of the ITP without slot to Theodorsen analytical solutions. Then, the unsteady aerodynamics of the ITP with and without central slot is investigated. It is found that even a smaller ratio of central slot of the ITP has notable effects on pressure distributions of the downstream section, and the pressure distributions on the downstream section will further be significantly affected by the slot ratio and the reduced wind speeds. Continuous increase of $A_2^*$ with the increase of central slot may be the key feature of the slotted ITP. Finally, flutter analyses based on the flutter derivatives of the slotted ITP are performed, and moreover, flutter instabilities of a scaled sectional model of a twin-deck bridge with various ratios of deck slot are investigated. The results confirm that the central slot is effective to improve bridge flutter stabilities, and that the flutter critical wind speeds increase with the increase of slot ratio.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1114-1124
• /
• 2017

Acceleration/preconditioning strategies available in the SCEPTRE radiation transport code are described. A flexible transport synthetic acceleration (TSA) algorithm that uses a low-order discrete-ordinates ($S_N$) or spherical-harmonics ($P_N$) solve to accelerate convergence of a high-order $S_N$ source-iteration (SI) solve is described. Convergence of the low-order solves can be further accelerated by applying off-the-shelf incomplete-factorization or algebraic-multigrid methods. Also available is an algorithm that uses a generalized minimum residual (GMRES) iterative method rather than SI for convergence, using a parallel sweep-based solver to build up a Krylov subspace. TSA has been applied as a preconditioner to accelerate the convergence of the GMRES iterations. The methods are applied to several problems involving electron transport and problems with artificial cross sections with large scattering ratios. These methods were compared and evaluated by considering material discontinuities and scattering anisotropy. Observed accelerations obtained are highly problem dependent, but speedup factors around 10 have been observed in typical applications.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.11 no.7
• /
• pp.1341-1347
• /
• 2007

The displacement field (Optical flow) has been calculated by bottom-up approaches based on local processing. In contrast with them, in this paper, a top-down approach based on expanding in turn from the lowest order mode the whole motion in an image pair of sequential images is proposed. The intensity of medical images usually represents a quantity which is conserved during the motion. Hence sequential images are ideally related by a coordinate transformation. The displacement field can be determined from the generalized moments of the two images. The equations which transform arbitrary generalized moments from a source image to a target image are expressed as a function of the displacement field. The appareent displacement field is then computed iteratively by a projection method which utilizes the functional derivatives of the linearized moment equations. This method is demonstrated using a pair of sequential heart images. For comparative evaluation, we applied Horn and Schunck's method, a standard multigrid method, and our proposed algorithm to sequential image.

• Journal of computational fluids engineering
• /
• v.9 no.1
• /
• pp.34-40
• /
• 2004

Recently with growth in the capability of super computers and Parallel computers, shape design optimization is becoming easible for real problems. Also, Computational Fluid Dynamics(CFD) techniques have been improved for higher reliability and higher accuracy. In the shape design optimization, analysis solvers and optimization schemes are essential. In this work, the Roe's 2nd-order Upwind TVD scheme and DADI time march with multigrid were used for the flow solution with the Euler equation and FDM(Finite Differenciation Method), GA(Genetic Algorithm) and Kriging were used for the design optimization. Kriging were applied to 2-D airfoil design optimization and compared with FDM and GA's results. When Kriging is applied to the nonlinear problems, satisfactory results were obtained. From the result design optimization by Kriging method appeared as good as other methods.