• Tribology and Lubricants
• /
• 제28권4호
• /
• pp.160-166
• /
• 2012

This paper presents a Reynolds equation solver for hydrostatic gas bearings, implemented to run on graphics processing units (GPUs). The original analysis code for the central processing unit (CPU) was modified for the GPU by using the compute unified device architecture (CUDA). The red-black Gauss-Seidel (RBGS) algorithm was employed instead of the original Gauss-Seidel algorithm for the iterative pressure solver, because the latter has data dependency between neighboring nodes. The implemented GPU program was tested on the nVidia GTX580 system and compared to the original CPU program on the AMD Llano system. In the iterative pressure calculation, the implemented GPU program showed 20-100 times faster performance than the original CPU codes. Comparison of the wall-clock times including all of pre/post processing codes showed that the GPU codes still delivered 4-12 times faster performance than the CPU code for our target problem.

• International Journal of Control, Automation, and Systems
• /
• 제6권5호
• /
• pp.776-784
• /
• 2008

In this paper, new upper matrix bounds for the solution of the continuous algebraic Riccati equation (CARE) are derived. Following the derivation of each bound, iterative algorithms are developed for obtaining sharper solution estimates. These bounds improve the restriction of the results proposed in a previous paper, and are more general. The proposed bounds are always calculated if the stabilizing solution of the CARE exists. Finally, numerical examples are given to demonstrate the effectiveness of the present schemes.

• 대한기계학회논문집
• /
• 제6권3호
• /
• pp.211-221
• /
• 1982

Numerical analysis has been made on the heat transfer of two dimensional turbulent channel with a slat type blockage. Especially the effects of the height of slat and Reynolds number on the heat transfer characteristics of channel wall have been investigated. The methods of accelerating the convergence of the numerical solution of governing differential equation have been also examined. Line-by-line iterative method shows higher convergence rate than point-by-point iterative method for solution of both momentum equation and energy equation. The results show that the ratio of heat transfer coefficient of the wall near the blockage to that of the fully developed flow increase with increasing the ratio of blockage to channel height and decreasing the Reynolds number. These trends of variation of heat transfer coefficient with respect to the height of slat and Reynolds number agree with those of Sparrow's experiment on the pipe flow with slat type blockage.

• 호남수학학술지
• /
• 제42권1호
• /
• pp.17-35
• /
• 2020

In this paper we introduce and consider a new class of variational inequalities with four operators. This class is called the extended general mixed quasi variational inequality. We show that the extended general mixed quasi variational inequality is equivalent to the fixed point problem. We use this alternative equivalent formulation to discuss the existence of a solution of extended general mixed quasi variational inequality and also develop several iterative methods for solving extended general mixed quasi variational inequality and its variant forms. We consider the convergence analysis of the proposed iterative methods under appropriate conditions. We also introduce a new class of resolvent equation, which is called the extended general implicit resolvent equation and establish an equivalent relation between the extended general implicit resolvent equation and the extended general mixed quasi variational inequality. Some special cases are also discussed.

• East Asian mathematical journal
• /
• 제24권3호
• /
• pp.305-315
• /
• 2008

The purpose of this paper is to study the existence and uniqueness of the fixed point of uniformly pseudo-contractive operator and the solution of equation with uniformly accretive operator, and to approximate the fixed point and the solution by the Mann iterative sequence in an arbitrary Banach space or an uniformly smooth Banach space respectively. The results presented in this paper show that if X is a real Banach space and A : X $\rightarrow$ X is an uniformly accretive operator and (I-A)X is bounded then A is a mapping onto X when A is continuous or $X^*$ is uniformly convex and A is demicontinuous. Consequently, the corresponding results which depend on the assumptions that the fixed point of operator and solution of the equation are in existence are improved.

• 대한수학회보
• /
• 제43권4호
• /
• pp.791-804
• /
• 2006

This paper is concerned with a second-order iterated differential equation of the form $c_0x'(Z)+c_1x'(z)+c_2x(z)=x(az+bx(z))+h(z)$ with the distinctive feature that the argument of the unknown function depends on the state. By constructing a convergent power series solution of an auxiliary equation, analytic solutions of the original equation are obtained.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제3권2호
• /
• pp.1-4
• /
• 1999

The Helmholtz equation is very important in physics and engineering. However, solution of the Helmholtz equation is in general known as a very difficult phenomenon. For if the ${\omega}$ is negative, the FDM discretized linear system becomes indefinite, whose solution by iterative method requires a very clever preconditioner. In this paper we assume that ${\omega}$ is nonnegative, and determine the optimal ${\rho}$ parameter for the three dimensional ADI iteration for the Helmholtz equation. The ADI(Alternating Direction Implicit) method is also getting new attentions due to the fact that it is very suitable to the vector/parallel computers, for example, as a preconditioner to the Krylov subspace methods. However, classical ADI was developed for two dimensions, and for three dimensions it is known that its convergence behaviour is quite different from that in two dimensions. So far, in three dimensions the so-called Douglas-Rachford form of ADI was developed. It is known to converge for a relatively wide range of ${\rho}$ values but its convergence is very slow. In this paper we determine the necessary conditions of the ${\rho}$ parameter for the convergence and optimal ${\rho}$ for the three dimensional ADI iteration of the Peaceman-Rachford form for the real Helmholtz equation with nonnegative ${\omega}$. Also, we conducted some experiments which is in close agreement with our theory. This straightforward extension of Peaceman-rachford ADI into three dimensions will be useful as an iterative solver itself or as a preconditioner to the the Krylov subspace methods, such as CG(Conjugate Gradient) method or GMRES(m).

• 대한기계학회논문집B
• /
• 제27권6호
• /
• pp.753-765
• /
• 2003

A finite element code for the numerical solution of the Navier-Stokes equation is parallelized by vertex-oriented domain decomposition. To accelerate the convergence of iterative solvers like conjugate gradient method, parallel block ILU, iterative block ILU, and distributed ILU methods are tested as parallel preconditioners. The effectiveness of the algorithms has been investigated when P1P1 finite element discretization is used for the parallel solution of the Navier-Stokes equation. Two-dimensional and three-dimensional Laplace equations are calculated to estimate the speedup of the preconditioners. Calculation domain is partitioned by one- and multi-dimensional partitioning methods in structured grid and by METIS library in unstructured grid. For the domain-decomposed parallel computation of the Navier-Stokes equation, we have solved three-dimensional lid-driven cavity and natural convection problems in a cube as benchmark problems using a parallelized fractional 4-step finite element method. The speedup for each parallel preconditioning method is to be compared using upto 64 processors.

• 한국전산응용수학회:학술대회논문집
• /
• 한국전산응용수학회 2002년도 학술발표대회
• /
• pp.11-11
• /
• 2002

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 2001년도 ICCAS
• /
• pp.169.3-169
• /
• 2001

This paper introduces a new approach to analysis of error convergence for a class of iterative learning control systems. First, a nonlinear plant is represented using a Takagi-Sugeno(T-S) fuzzy model. Then each iterative learning controller is designed for each linear plant in the T-S fuzzy model. From the view point of two-dimensional(2-D) system theory, we transform the proposed learning systems to a 2-D error equation, which is also established in the form of T-S fuzzy model. We analysis the error convergence in the sense of induced 2 L -norm, where the effects of disturbances and initial conditions on 2-D error are considered. The iterative learning controller design problem to guarantee the error convergence can be reduced to linear matrix inequality problems. In comparison with others, our learning algorithm ...