• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.207-216
• /
• 2023

Three time-discontinuous Galerkin quadrature element methods (TDGQEMs) are developed for structural dynamic problems. The weak-form time-discontinuous Galerkin (TDG) statements, which are capable of capturing possible displacement and/or velocity discontinuities, are employed to formulate the three types of quadrature elements, i.e., single-field, single-field/least-squares and two-field. Gauss-Lobatto quadrature rule and the differential quadrature analog are used to turn the weak-form TDG statements into a system of algebraic equations. The stability, accuracy and numerical dissipation and dispersion properties of the formulated elements are examined. It is found that all the elements are unconditionally stable, the order of accuracy is equal to two times the element order minus one or two times the element order, and the high-order elements possess desired high numerical dissipation in the high-frequency domain and low numerical dissipation and dispersion in the low-frequency domain. Three fundamental numerical examples are investigated to demonstrate the effectiveness and high accuracy of the elements, as compared with the commonly used time integration schemes.

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

• Communications of the Korean Mathematical Society
• /
• v.18 no.4
• /
• pp.767-779
• /
• 2003

We consider a model of population dynamics whose mortality function is unbounded. We approximate the solution of the model using a discontinuous Galerkin finite element for the age variable and a backward Euler for the time variable. We present several numerical examples. It is experimentally shown that the scheme converges at the rate of $h^{3/2}$ in the case of piecewise linear polynomial space.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.57-70
• /
• 2008

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.57-70
• /
• 2008

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.

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

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.83-88
• /
• 1993

A time-discontinuous Galerkin method based upon using a finite element formulation in time has evolved. This method, working from the differential equation viewpoint, is different from those which have been generally used. They admit discontinuities with respect to the time variable at each time step. In particular, the elements can be chosen arbitrarily at each time step with no connection with the elements corresponding to the previous step. Interpolation functions and weighting functions are taken to be discontinuous across inter-element boundaries. These methods lead to a unconditional stable higher-order accurate ordinary differential equation solver.

• Journal of the Korean Mathematical Society
• /
• v.51 no.4
• /
• pp.665-678
• /
• 2014

We present a new space-time discontinuous Galerkin (DG) method for solving the time dependent, positive symmetric hyperbolic systems. The main feature of this DG method is that the discrete equations can be solved semi-explicitly, layer by layer, in time direction. For the partition made of triangle or rectangular meshes, we give the stability analysis of this DG method and derive the optimal error estimates in the DG-norm which is stronger than the $L_2$-norm. As application, the wave equation is considered and some numerical experiments are provided to illustrate the validity of this DG method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.3
• /
• pp.169-180
• /
• 2009

A discontinuous Galerkin method with interior penalty terms is presented for linear Sobolev equation. On appropriate finite element spaces, we apply a symmetric interior penalty Galerkin method to formulate semidiscrete approximate solutions. To deal with a damping term $\nabla{\cdot}({\nabla}u_t)$ included in Sobolev equations, which is the distinct character compared to parabolic differential equations, we choose special test functions. A priori error estimate for the semidiscrete time scheme is analyzed and an optimal $L^\infty(L^2)$ error estimation is derived.