• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.31-35
• /
• 2006

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.10-18
• /
• 2016

Numerical boundary conditions are as important as the governing equations when analyzing the fluid flows numerically. An explicit boundary condition method updates the solutions at the boundaries with extrapolation from the interior of the computational domain, while the implicit boundary condition method in conjunction with an implicit time integration method solves the solutions of the entire computational domain including the boundaries simultaneously. The implicit boundary condition method, therefore, is more robust than the explicit boundary condition method. In this paper, steady compressible 2-Dimensional Navier-Stokes solver is developed. We present the implicit boundary condition method coupled with LU-SGS(Lower Upper Symmetric Gauss Seidel) method. Also, the explicit boundary condition method is implemented for comparison. The preconditioning Navier-Stokes equations are solved on unstructured meshes. The numerical computations for a number of flows show that the implicit boundary condition method can give accurate solutions.

• Journal of computational fluids engineering
• /
• v.12 no.3
• /
• pp.29-40
• /
• 2007

An implicit discontinuous Galerkin method for the two-dimensional Euler equations was developed on unstructured triangular meshes. The method can achieve high-order spatial accuracy by using hierachical basis functions based on Legendre polynomials. Numerical tests were conducted to estimate the convergence order of numerical solutions to the Ringleb flow and the supersonic vortex flow for which analytic solutions are available. Also, the flows around a 2-D circular cylinder and an NACA0012 airfoil were numerically simulated. The numerical results showed that the implicit discontinuous Galerkin methods couples with a high-order representation of curved solid boundaries can be an efficient method to obtain very accurate numerical solutions on unstructured meshes.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.30-40
• /
• 2007

The implicit discontinuous Galerkin method for the two-dimensional Euler equations was developed on unstructured triangular meshes, which can achieve higher-order accuracy by wing hierachical basis functions based on Legendre polynomials. Numerical tests were conducted to estimate the convergence order of numerical solutions to the Ringleb flow and the supersonic vortex flow for which analytic solutions are available. And, the flows around a circle and a NACA0012 airfoil was also numerically simulated. Numerical results show that the implicit discontinuous Galerkin methods with higher-order representation of curved solid boundaries can be an efficient higher-order method to obtain very accurate numerical solutions on unstructured meshes.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.23 no.4
• /
• pp.337-344
• /
• 2014

While implicit integration methods such as Newton's method have excellent stability for the analysis of stiff and constrained mechanical systems, they have the drawback that the evaluation and LU-factorization of the system Jacobian matrix required at every time step are time-consuming. This paper proposes a Jacobian update-free Newton's method in order to overcome these defects. Because the motions of all bodies in a vehicle model are limited with respect to the chassis body, the equations are formulated with respect to the moving chassis-body reference frame instead of the fixed inertial reference frame. This makes the system Jacobian remain nearly constant, and thus allows the Newton's method to be free from the Jacobian update. Consequently, the proposed method significantly decreases the computational cost of the vehicle dynamic simulation. This paper provides detailed generalized formulation procedures for the equations of motion, constraint equations, and generalized forces of the proposed method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.10a
• /
• pp.92-98
• /
• 1998

Nonlinear finite element analyses of one dimensional consolidation problem were performed using an anisotropic hardening constitutive model. For the analyses, the anisotropic hardening elasto-plastic constitutive model based on the generalized isotropic hardening(GIH) rule was implemented into a nonlinear finite element analysis program, PLASTIC. In order to preserve the accuracy of the finite element solution for nonlinear problems, an implicit stress integration algorithm was employed. A consistent tangent moduli could also ensure the quadratic convergence of Newton's method. As a result, the nonlinear solution was accurately calculated and was converged to be asymptotically quadratic. In a consolidation problem, the relationship between load and settlement and between settlement and time vertical was analyzed comparing with results using the Cam-clay type model and the final consolidation settlement and the duration of primary consolidation could be evaluated rigorously using the GIH constitutive model.

• Land and Housing Review
• /
• v.1 no.1
• /
• pp.59-65
• /
• 2010

The paper presents a comparison between a semi-implicit time integration linear finite element implementation and fully-implicit nonlinear Newton-Raphson finite element implementation of a triphasic small strain mixture formulation of an elastic partially saturated porous medium. The pore air phase pressure pa is assumed atmospheric, i.e., $p_a$ = 0, although the formulation and implementation are general to handle increase in pore air pressure as a result of loading, if needed. The solid skeleton phase is assumed linear isotropic elastic and partially saturated 'consolidation' in the presence of surface infiltration and traction is simulated. The verification of the implementation against an analytical solution for partially saturated pore water flow (no deformation) and comparison between the two implementations is presented and the important of the porosity-dependent nature of the partially saturated permeability is assessed on comparison with a commercial code for the partially saturated flow with deformation. As a result, the response of partially saturated permeability subjected to the porosity influences on the saturation of a soil, and the different behaviors of the partially saturated soil between staggered and monolithic coupled programs is worth of attention because the negative pore water pressure in the partially saturated soil depends on the difference.

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

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

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.9 no.5
• /
• pp.995-1001
• /
• 2005

Implicit Integration has become a standard approach to efficient cloth animation, and it guarantees the stability of the system so that large steps can be used. Therefore, it is regarded as the best method for the real-time or interactive animation of cloth. Since the implicit method was introduced for stable cloth animation, various cloth animation techniques based on the method have been proposed. It is now possible to generate the real-time animation of cloth model with thousands of mass-point in general PC environments. Although the implicit method guarantees the stability, the implementation of the implicit method is generally more difficult than that of the explicit method. Even worse, it is very difficult to parallelize the computation process of the implicit method. The cloth animation with implicit method can be formalized as a linear system solving. In this paper we propose an stable and efficient cloth animation techniques based on the implicit method. The proposed method can be easily parallelized. Self-collision is another important issue in cloth animation, we also propose an efficient self-collision avoidance techniques.