• Proceedings of the Korea Water Resources Association Conference
• /
• 1982.07a
• /
• pp.91-98
• /
• 1982

The formulation of depth-averaged two-dimensional mathematical model for the analysis of tide induced circulation in a harbor by the Galerkin finite element techique is presented. In integration of the Galerkin approach in time both explicit and implicit method have been tested for one and two dimentional water bodies, and the two step Lax-Wendroff explicit method is found to be effective than the implicit in reducing computing time. The essential characteristics of the tide induced flow in Busan Harbor with two open boundaries has been foccud to be reproduceable in the numerical model and the simulated results encourage that the model can be used as a predictive tool.

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

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

• Journal of Aerospace System Engineering
• /
• v.14 no.3
• /
• pp.17-23
• /
• 2020

In this paper, the efficient periodic unsteady flow analysis is developed by using a Chebyshev collocation operator applied to the time differential term of the governing equations. The partial implicit time integration method was also applied in the governing equation for a fluid, which means flux terms were implicitly processed for a time integration and the time derivative terms were applied explicitly in the form of the source term by applying the Chebyshev collocation operator. To verify this method, we applied the 1D unsteady Burgers equation and the 2D oscillating airfoil. The results were compared with the existing unsteady flow frequency analysis technique, the Harmonic Balance Method, and the experimental data. The Chebyshev collocation operator can manage time derivatives for periodic and non-periodic problems, so it can be applied to non-periodic problems later.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.11a
• /
• pp.848-854
• /
• 2006

In this study, a computational analysis system has been developed in order to investigate flow-induced vibration(FIV) phenomenon for general stator-rotor cascade configurations. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. Fluid domains are modeled using the unstructured grid system with dynamic moving and local deforming methods. Unsteady, Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-\omega$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used flow computing the coupled governing equations of the fluid-structure interaction problem. Detailed FIV responses for different flow conditions are presented with respect to time and vibration characteristics are also physically investigated in the time domain.

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

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• Transactions of Materials Processing
• /
• v.3 no.1
• /
• pp.84-96
• /
• 1994

Simulation of 3 dimensional large irregularly shaped stamping process by a dynamic approach, based on an explicit time integration scheme, has been shown to be highly efficient and robust in comparison to traditional, implicit, quasi-static ones. The objective of the work is to evaluate the results from explicit code in application to deep drawing of rectangular cup and stamping of automotive front fender, in which deformation, force, thickness distribution are calculated. The method of predicting spring back and formability by and explicit code are suggested and applied to the processes.

• Journal of Mechanical Science and Technology
• /
• v.14 no.7
• /
• pp.789-795
• /
• 2000

Free-surface flows with an arbitrary deformation, induced by a submerged hydrofoil, are simulated numerically, considering two-fluid flows of both water and air. The computation is performed by a finite volume method using unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell wise local mesh refinement. The integration in space is of second order, based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels. The linear equations are solved by conjugate gradient type solvers, and the non-linearity of equations is accounted for through Picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations, the continuity equation, the conservation equation of one species, and the equations for two turbulence quantities. Finally, a comparison is quantitatively made at the same speed between the computation and experiment in which the grid sensitivity is numerically checked.