• Journal of Korea Game Society
• /
• v.18 no.5
• /
• pp.123-132
• /
• 2018

Various methods have been proposed to efficiently animate the motion of soft objects in realtime. In order to maintain the topology between the elements of the objects, it is required to employ constraint forces, which limit the size of the time steps for the numerical integration and reduce the efficiency. To tackle this, an implicit method with larger steps was proposed. However, the method is, in essence, a linear system with a large matrix, of which solution requires heavy computations. Several approximate methods have been proposed, but the approximation is obtained with an increased damping and the loss of accuracy. In this paper, new integration method based on harmonic oscillation with better stability was proposed, and it was further stabilized with the hybridization with approximate implicit method. GPU parallelism can be easily implemented for the method, and large-scale soft objects can be simulated in realtime.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.48-54
• /
• 1998

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

• Journal of The Korean Astronomical Society
• /
• v.40 no.4
• /
• pp.91-97
• /
• 2007

The Fokker-Planck (FP) model is one of the commonly used methods for studies of the dynamical evolution of dense spherical stellar systems such as globular clusters and galactic nuclei. The FP model is numerically stable in most cases, but we find that it encounters numerical difficulties rather often when the effects of tidal shocks are included in two-dimensional (energy and angular momentum space) version of the FP model or when the initial condition is extreme (e.g., a very large cluster mass and a small cluster radius). To avoid such a problem, we have developed a new integration scheme for a two-dimensional FP equation by adopting an Alternating Direction Implicit (ADI) method given in the Douglas-Rachford split form. We find that our ADI method reduces the computing time by a factor of ${\sim}2$ compared to the fully implicit method, and resolves problems of numerical instability.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.113-117
• /
• 1995

Viscous analysis on incompressible flows is performed using unstructured triangular meshes. A two-dimensional and axisymmetric incompressible Navier-Stokes equations are solved in time-marching form by artificial compressibility method. The governing equations are discretized by a cell-centered based finite-volume method. and a centered scheme is used for inviscid and viscous fluxes with fourth order artificial dissipation. An explicit multi-stage Runge-Kutta method is used for the time integration with local time stepping and implicit residual smoothing. Convergence properties are examined and solution accuracies are also validated with benchmark solution and experiment.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.12
• /
• pp.1433-1441
• /
• 2009

In multibody dynamic analysis, one of the most important problems is to reduce computation times for real-time simulation. This paper presents the derivation procedure of equations of motion of a 6${\times}$6 autonomous vehicle in terms of chassis local coordinates which do not require coordinates transformation matrix to enhance efficiency for real-time dynamic analysis. Also, equations of motion are derived using the VT(velocity transformation) technique and symbolic computation method coded by MATLAB. The Jacobian matrix of the equations of motion of a system is derived from symbolic operations to apply the implicit integration method. The analysis results were compared with ADAMS results to verify the accuracy and approve the feasibility of real time analysis.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1996.06a
• /
• pp.97-109
• /
• 1996

In the present work computations are carried out for analysis of complicated sheet metal forming process such as forming of a rear hinge. Finite element formulation using dynamic explicit time integration scheme and step-wise combined Implicit/Explicit scheme are introduced for numerical analysis of sheet metal forming process. The rigid-plastic finite element method based on membrane elements has long been employed as a useful numerical technique for the analysis of sheet metal forming because of its time effectiveness. The explicit scheme in general use is based on the elastic-plastic modelling of material requiring large computation time. In finite element simulation of sheet metal forming processes, the robustness and stability of computation are important requirements since the computation time and convergency become major points of consideration besides the solution accuracy due to the complexity of geometry and boundary conditions. The implicit scheme employs a more reliable and rigorous scheme in considering the equilibrium at each step of deformation, while in the explicit scheme the problem of convergency is eliminated at the cost of solution accuracy. The explicit approach and the implicit approach have merits and demerits, respectively. In order to combine the merits of these two methods a step-wise combined implicit/explicit scheme has been developed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.

• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• Proceedings of the IEEK Conference
• /
• 2003.11b
• /
• pp.125-128
• /
• 2003

This paper proposes efficient cloth animation and collision handling methods. There have been various techniques for the generation of cloth behavior. However, the cloth animation is still a challenging subject in real-time environments. This paper presents an efficient animation method based on implicit integration. The proposed method can efficiently animate virtual cloth object with complex geometry. In addition, this paper also introduces an efficient collision handling method. The collision resolution is another important issue in cloth animation since deformable objects has special collision problem called self-collision. In this paper, the self-collision was successfully avoided in real - time environments.