• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.60-67
• /
• 2008

In this paper we consider an application of the method of electromagnetic (EM) migration to the interpretation of a typical marine controlled-source (MCSEM) survey consisting of a set of sea-bottom receivers and a moving electrical bipole transmitter. Three-dimensional interpretation of MCSEM data is a very challenging problem because of the enormous number of computations required in the case of the multi-transmitter and multi-receiver data acquisition systems used in these surveys. At the same time, we demonstrate that the MCSEM surveys with their dense system of transmitters and receivers are extremely well suited for application of the migration method. In order to speed up the computation of the migration field, we apply a fast form of integral equation (IE) solution based on the multigrid quasi-linear (MGQL) approximation which we have developed. The principles of migration imaging formulated in this paper are tested on a typical model of a sea-bottom petroleum reservoir.

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

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.690-694
• /
• 2008

In this presentation, we study numerically an optimization problem of a fire fighting robot cooling system. The governing equation for the system is the unsteady heat equation with source term. We use a multigrid method for numerical solutions in three-dimensional space. We investigated the effects of various parameters and the results will be presented.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.103-121
• /
• 2015

In this paper, we briefly review and describe a projection algorithm for numerically computing the two-dimensional time-dependent incompressible Navier-Stokes equation. The projection method, which was originally introduced by Alexandre Chorin [A.J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comput., 22 (1968), pp. 745-762], is an effective numerical method for solving time-dependent incompressible fluid flow problems. The key advantage of the projection method is that we do not compute the momentum and the continuity equations at the same time, which is computationally difficult and costly. In the projection method, we compute an intermediate velocity vector field that is then projected onto divergence-free fields to recover the divergence-free velocity. Numerical solutions for flows inside a driven cavity are presented. We also provide the source code for the programs so that interested readers can modify the programs and adapt them for their own purposes.

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

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

A dual-time stepping algorithm combined with a parallelized multigrid DADI method is presented to predict the dynamic damping coefficients. The Basic Finner model is chosen to validate the prediction capability of the present unsteady Euler method. The linearity of the pitch- and roll-damping coefficients is shown in the low angular rates and the interesting large drop and stiff increment in transonic region for roll-damping coefficients are explained in detail. Through the analysis for the pressure distributions at Mach number 1.0 to 1.2, the sudden drop results from the normal shock and the stiff increment of roll-damping reflects the transition of the normal shock to the oblique shock. The results also show that the Euler equations can give the damping coefficients with a comparable accuracy.

• Communications of the Korean Mathematical Society
• /
• v.20 no.1
• /
• pp.179-193
• /
• 2005

We present results of fully nonlinear time-dependent simulations of a thin liquid film flowing up an inclined plane. Equations of the type $h_t+f_y(h) = -{\in}^3{\nabla}{\cdot}(M(h){\nabla}{\triangle}h)$ arise in the context of thin liquid films driven by a thermal gradient with a counteracting gravitational force, where h = h(x, t) is the fluid film height. A hybrid scheme is constructed for the solution of two-dimensional higher-order nonlinear diffusion equations. Problems in the fluid dynamics of thin films are solved to demonstrate the accuracy and effectiveness of the hybrid scheme.

• Tribology and Lubricants
• /
• v.25 no.5
• /
• pp.335-341
• /
• 2009

Elastohydrodynamic lubrication (EHL) analysis shows that film thickness is very flat in the contact area and pressure distribution is somehow similar to that of Hertzian contact pressure except the outlet region with pressure spike. These typical patterns of EHL film thickness and pressure are the cases under the steady contact conditions of applied loads and speeds. However, many engineering contacts are rather under the conditions of varying loads and contact speeds, and therefore the predictions for endurance life and performance of machine elements with steady EHL analysis are not suitable in many occasions. This study shows the differences in film thickness formation and pressure distribution between steady and transient contact conditions in several contact cases.

• Proceedings of the Korean Fiber Society Conference
• /
• 2003.10a
• /
• pp.47-48
• /
• 2003

Dilute bubble suspensions are prepared by introducing carbon dioxide bubbles into polyol resin. The apparent shear viscosity is measured with a wide gap parallel plate rheometer. A numerical simulation for deformation of a single bubble suspended in a Newtonian fluid is conducted by using a finite volume method (FVM) where multigrid algorithms are incorporated. Transient and steady results of bubble deformation were obtained and were in good agreement with experimental results. At high capillary number, viscosity of the suspension increases as the volume fraction increases, while at low capillary number, the viscosity decreases as the volume fraction increases.

• Tribology and Lubricants
• /
• v.21 no.6
• /
• pp.272-277
• /
• 2005

Continuously variable transmission (CVT) of toroidal type has a elliptical shape of contact zone under the elastohydrodynamic lubrication (EHL) condition, where the power is transmitted only by shearing the lubricant. Due to the small contact area of elliptical shape, the traction of the shear behaviors of lubricant over the contact zone is under extremely high contact pressure over 1.0GPa. During the power transmission by shearing the fluid, many kinds of mechanical movements occur such as squeezing, sliding, rolling and spin. Among the movements, the spin effect that is the most undesirable contact behavior in transmitting the power frequently makes significant abnormal wear damage. In this work, the analysis of elliptical contact of EHL with spin effect is performed, which will give very useful information to understand the traction behaviors in toroidal type of CVT system.