• 한국전산유체공학회지
• /
• 제18권3호
• /
• pp.34-41
• /
• 2013

We develop a novel immersed boundary (IB) method based on implicit direct forcing scheme for incompressible flows. The proposed IB method is based on an iterative procedure for calculating the direct forcing coupled with the momentum equations in order to satisfy no-slip boundary conditions on IB surfaces. We perform simulations of two-dimensional flows over a circular cylinder for low and moderate Reynolds numbers. The present method shows that the errors for estimated velocities on IB surfaces are significantly reduced even for low Reynolds number with a fairly large time step while the previous methods based on direct forcing failed to provide no-slip boundary conditions on IB surfaces.

• 대한수학회논문집
• /
• 제11권2호
• /
• pp.495-501
• /
• 1996

There is currently a regain of interest in ADI (Alternating Direction Implicit) method as a preconditioner for iterative Method for solving large sparse linear systems, because of its suitability for parallel computation. However the classical ADI is not applicable to FE(Finite Element) matrices. In this paper wer propose a Block-ADI method, which is applicable to Finite Element metrices. The new approach is a combination of classical ADI method and domain decompositi on. Also, we provide a partial proof of the convergence based on the results from the regular splittings, in case the discretization metrix is symmetric positive definite.

• East Asian mathematical journal
• /
• 제40권1호
• /
• pp.109-118
• /
• 2024

This paper proposes a numerical method based on domain decomposition to find approximate solutions for one-dimensional convection-diffusion equations with Neumann boundary conditions. First, the equations are transformed into convection-diffusion equations with Dirichlet conditions. Second, the author introduces the Prediction/Correction Domain Decomposition (PCDD) method and estimates errors for the interface prediction scheme, interior scheme, and correction scheme using known error estimations. Finally, the author compares the PCDD algorithm with the fully explicit scheme (FES) and the fully implicit scheme (FIS) using three examples. In comparison to FES and FIS, the proposed PCDD algorithm demonstrates good results.

• Earthquakes and Structures
• /
• 제7권2호
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2007년도 춘계학술대회A
• /
• pp.512-517
• /
• 2007

In the previous development of the recursive thermostat chained fully flexible cell molecular dynamics simulation, implicit time integration method such as generalized leapfrog integration is used. The implicit algorithm is very much complicated and not easy to show time reversibility because it is solved by the nonlinear iterative procedure. Thus we develop simple, explicit symplectic time integration formula for the recursive thermostat chained fully flexible unit cell simulation. Uniaxial tension test is performed to verify the present explicit algorithm. We check that the present simulation satisfies the ergodic hypothesis for various values of fictitious mass and coefficient of multiple thermostat system. The proposed method should be helpful to predict mechanical and thermal behavior of nano-scale structure.

• 대한기계학회논문집
• /
• 제18권6호
• /
• pp.1421-1429
• /
• 1994

An algorithm is developed for solving the inverse kinematic problem of a 6-degree-of-freedom robot with a wrist offset for which the closed form inverse solutions are not obtainable, but knowledge of one joint variable allows closed form solutions of the remaining joint variables. The algorithm does not require Forward Kinematics nor Jacobian but uses the implicit kinematic relationships between joint variables and the given hand position. An iterative back substitution method is used to solve the inversion and the optimal conditions of the convergence are incoporated. An example is given to illustrate the concepts, the solution procedure and its convergency.

• Structural Engineering and Mechanics
• /
• 제56권4호
• /
• pp.549-567
• /
• 2015

A first-order moment method (FORM) reliability analysis is commonly used for structural stability analysis. It requires the values and partial derivatives of the performance to function with respect to the random variables for the design. These calculations can be cumbersome when the performance functions are implicit. A Gaussian process (GP)-based response surface is adopted in this study to approximate the limit state function. By using a trained GP model, a large number of values and partial derivatives of the performance functions can be obtained for conventional reliability analysis with a FORM, thereby reducing the number of stability analysis calculations. This dynamic renewed knowledge source can provide great assistance in improving the predictive capacity of GP during the iterative process, particularly from the view of machine learning. An iterative algorithm is therefore proposed to improve the precision of GP approximation around the design point by constantly adding new design points to the initial training set. Examples are provided to illustrate the GP-based response surface for both structural and non-structural reliability analyses. The results show that the proposed approach is applicable to structural reliability analyses that involve implicit performance functions and structural response evaluations that entail time-consuming finite element analyses.

• Journal of Mechanical Science and Technology
• /
• 제15권1호
• /
• pp.125-138
• /
• 2001

An efficient correction storage scheme on a structured grid is applied to a sequence of approximate Jacobian systems arising at each time step from a linearization of the discrete nonlenear system of equations, obtained by the implicit time discretization of the conservation laws for unsteady fluid flows. The contribution of freezing the Jacobian matrix to computing costs is investigated within the correction storage scheme. The performance of the procedure is exhibited by measuring CPU time required to obtain a fully developed laminar vortex shedding flow past a circular cylinder, and is compared with that of a collective iterative method on a single grid. In addition, some computed results of the flow are presented in terms of some functionals along with measured data. The computational test shows that the computing costs may be saved in favor of the correction storage scheme with the frozen Jacobian matrix, to a great extent.

• Journal of Mechanical Science and Technology
• /
• 제16권10호
• /
• pp.1276-1286
• /
• 2002

The numerical techniques are presented to solve the static and dynamic contact problems of deformable bodies having large rotations of the contact surfaces. The contact conditions on the possible contact surfaces are enforced by using the contact error vector, and an iterative scheme similar to augmented Lagrange multiplier method is employed to reduce the contact error vector monotonically. For dynamic contact problems using implicit time integration, a contact error vector is also defined by combining the displacement, velocity, and acceleration on the contact surface. The suggested iterative technique is implemented to ABAQUS by using the UEL subroutine UEL. In this work, after the computing procedures to solve the frictional contact problems are explained, the numerical examples are presented to compare the present solutions with those obtained by ABAQUS.