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

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

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.680-687
• /
• 2007

The paper presents a numerical simulation of flow of a turbulent boundary layer, representing a typical wind environment and matching a series of wind tunnel observations. The simulations are carried out at a Reynolds number of 20,000, based on the velocity U at a pseudo-height h, and large enough that the flow be effectively Reynolds number independent. Some wall models are proposed for the LES(Large Eddy Simulation) of the turbulent boundary layer over a rough surface. The Jenson number, $J=h/z_0$, based on the roughness length $z_0$, is 600 to match the wind tunnel data. The computational mesh is uniform with a spacing of h/32, as this aids rapid convergence of the multigrid solver, and the governing equations are discretised using second order finite differences within a parallel multiblock environment. The results presented include the comparison between wind tunnel measurements and LES computations of the turbulent boundary layer over rough surface.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

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

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

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.7-14
• /
• 2004

It is well known that preconditioning methods are efficient for convergence acceleration at compressible low Mach number flows. In this study, the original Euler equations and three preconditioners nondimensionalized differently are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as flux discretization and time integration respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning one produces Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

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

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