• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

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

• Journal of the Korea Computer Graphics Society
• /
• v.5 no.1
• /
• pp.27-33
• /
• 1999

In this paper, we propose an efficient technique for the animation of flexible thin objects. Mass-spring model was employed to represent the flexible objects. Till now, many techniques have used the mass-spring model to generate plausible animation of soft objects. A straight-forward approach to the animation with mass-spring model is explicit Euler method, but the explicit Euler method has serious disadvantage that it suffers from 'instability problem'. The implicit integration method is a possible solution to overcome the instability problem. However, the most critical flaw of the implicit method is that it involves a large linear system. This paper presents a fast animation technique for mass-spring model with approximated implicit method. The proposed technique stably updates the state of n mass-points in O(n) time when the number of total springs are O(n). We also consider the interaction of the flexible object and air in order to generate plausible results.

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

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

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.186-191
• /
• 1999

An Euler solver is developed to predict accurate aerodynamic data such as lift coefficient, drag coefficient, and moment coefficient. The conservation law form of the compressible Euler equations are used in the generalized curvilinear coordinates system. The Euler solver uses a finite volume method and the second order Roe's flux difference splitting scheme with min-mod flux limiter to calculate the fluxes accurately. An implicit scheme which includes the boundary conditions is implemented to accelerate the convergence rate. The multi-block grid is integrated into the flow solver for complex geometry. The flowfields are analyzed around NACA 0012 airfoil in the cases of $M_{\infty}=0.75,\;\alpha=2.0\;and\;M_{\infty}=0.80,\;\alpha=1.25$. The numerical results are compared with other numerical results from the literature. The final goal of this research is to prepare a robust and an efficient Navier-Stokes solver eventually.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.93-98
• /
• 2001

An adaptive Cartesian grid method having the best elements of structured, unstructured, and Cartesian grids is developed to solve the steady two-dimensional Euler equations. The solver is based on a cell-centered finite-volume method with Roe's flux-difference splitting and implicit point Gauss-seidel time integration method. Calculations of several compressible flows are carried out to show the efficiency of the developed computer code. The results were generally in good agreements with existing data in the literature and the developed code has the good ability to capture important feature of the flows.

• 한국전산유체공학회:학술대회논문집
• /
• 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 computational fluids engineering
• /
• v.2 no.1
• /
• pp.13-20
• /
• 1997

A fundamental study for the numerical scheme to simulate unsteady nonlinear waves by solving Euler equations is presented. First a conservation form and a non-conservation form of the Euler equations with a free surface fitted coordinate system are compared. Next, a time splitting fractional step method and an alternating direction implicit(ADI) method for the time integration are compared. For the comparative study, flow calculations around a bottom bump in a channel and a NACA 0012 hydrofoil in a flume are performed. The results show that the ADI method with a third order upwind differencing scheme is very efficient in reducing the computing time with keeping the accuracy, And, there is no distinct difference between two expression forms except that the non-conservative form shows faster wave propagating velocity than the conservation form. Some results are compared with experiments and show good agreement.

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

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.