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

• Proceedings of the KIEE Conference
• /
• 2005.11b
• /
• pp.217-219
• /
• 2005

This paper describes implicit restated arnoldi method algorithm and Its application to small size power systems In order to observe the salient features of IRAM algorithm. Two area system with 36 state variables and England 39-bus system with 150 state variables have been tested using IRAM, and the eigenvalue results of IRAM are compared with those of the results obtained from QR method.

• Korean Journal of Computational Design and Engineering
• /
• v.12 no.3
• /
• pp.220-232
• /
• 2007

Current geometric modeling and analysis are commonly based on B-Rep modeling and a finite elements method respectively. Furthermore, it is difficult to represent an object whose material property is heterogeneous using the B-Rep method because the B-Rep is basically used for homogeneous models. In addition, meshes are required to analyze a property of a model when the finite elements method is applied. However, the process of generating meshes from B-Rep is cumbersome and sometimes difficult especially when the model is deformed as time goes by because the topology of deforming meshes are changed. To overcome those problems in modeling and analysis including homogeneous and heterogeneous materials, we suggest a unified modeling and analysis method based on implicit representation of the model using R-function which is suggested by Rvachev. For implicit modeling of an object a distance field is approximated and blended for a complex object. Using the implicit function mesh-free analysis is possible where meshes are not necessary. Generally mesh-free analysis requires heavy computational cost compared to a finite elements method. To improve the computing time of function evaluation, we utilize GPU programming. Finally, we give an example of a simple pipe design problem and show modeling and analysis process using our unified modeling and analysis method.

• Bulletin of the Korean Mathematical Society
• /
• v.29 no.1
• /
• pp.31-40
• /
• 1992

We introduce the generalized Runge-Kutta methods with the exponentially dominant order .omega. in [3], and the convergence theorems of the generalized explicit Euler method are derived in [4]. In this paper we will study the convergence of the generalized implicit Euler method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.18 no.1
• /
• pp.15-22
• /
• 1985

The study on computation time, accuracy, and convergency characteristic of the implicit finite difference method is presented with the variation of the space increment and time step in a two-dimensional transient heat conduction problem with a dirichlet boundary condition. Numerical analysis were conducted by the model having the conditions of the solution domain from 0 to 3m, thermal diffusivity of 1.26 $m^2/h$, initial condition of 272 K, and boundary condition of 255.4 K. The results obtained are summarized as follows : 1) The degree of influence with respect to the accuracy of the time step and space increment in the alternating-direction implicit method and Crank-Nicholson implicit method were relatively small, but in case of the fully implicit method showed opposite tendency. 2) To prescribe near the zero for the space increment and tine step in a two dimensional transient problem were good in a accuracy aspect but unreasonable in a computational time aspect. 3) The reasonable condition of the space increment and the time step considering accuracy and computation time could be generalized with the Fourier modulus increment, F, ana dimensionless space increment, X, irrespective of the solution domain.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.73-77
• /
• 2004

In the present study, an efficient and robust implicit operator for the LU-SGS method is proposed. Numerical experiments for supersonic flow are performed to demonstrate the performance of the proposed method.

• Bulletin of the Korean Mathematical Society
• /
• v.55 no.3
• /
• pp.881-898
• /
• 2018

This article deals with implementation of a high-order finite difference scheme for numerical solution of the incompressible Navier-Stokes equations on curvilinear grids. The numerical scheme is based on pseudo-compressibility approach. A fifth-order upwind compact scheme is used to approximate the inviscid fluxes while the discretization of metric and viscous terms is accomplished using sixth-order central compact scheme. An implicit Euler method is used for discretization of the pseudo-time derivative to obtain the steady-state solution. The resulting block tridiagonal matrix system is solved by approximate factorization based alternating direction implicit scheme (AF-ADI) which consists of an alternate sweep in each direction for every pseudo-time step. The convergence and efficiency of the method are evaluated by solving some 2D benchmark problems. Finally, computed results are compared with numerical results in the literature and a good agreement is observed.

• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.