• East Asian mathematical journal
• /
• v.15 no.1
• /
• pp.121-133
• /
• 1999

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.474-479
• /
• 2004

According to our previous study, it is confirmed that the Petrov-Galerkin natural element method (PGNEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method (BG-NEM). This paper is an extension of PG-NEM to two-dimensional nonlinear dynamic problem. For the analysis, a constant average acceleration method and a linearized total Lagrangian formulation is introduced with the PG-NEM. At every time step, the grid points are updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates the nonlinear dynamic problem.

• Journal of applied mathematics & informatics
• /
• v.28 no.3_4
• /
• pp.953-966
• /
• 2010

In this paper, we develop discontinuous Galerkin methods with penalty terms, namaly symmetric interior penalty Galerkin methods to solve nonlinear parabolic equations. By introducing an appropriate projection of u onto finite element spaces, we prove the optimal convergence of the fully discrete discontinuous Galerkin approximations in ${\ell}^2(L^2)$ normed 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.

• Nuclear Engineering and Technology
• /
• v.49 no.1
• /
• pp.29-42
• /
• 2017

The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1274-1279
• /
• 2003

In this paper, a new meshfree technique which improves the numerical integration accuracy is introduced. This new method called the Petrov-Galerkin natural element(PG-NE) is based on the Voronoi diagram and the Delaunay triangulation which is based on the same concept used for conventional natural element method called the Bubnov-Galerkin natural element(BG-NE). But, unlike BG-NE method, the test shape function is differently chosen from the trial shape function. The proposed technique ensures that the numerical integration error is remarkably reduced.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.6 s.237
• /
• pp.703-710
• /
• 2005

A finite element discretization of the advection and redistancing equations of level set method has been studied. It has been shown that Galerkin spatial discretization combined with Crank-Nicolson temporal discretization of the advection equation of level set yields a good result and that consistent streamline upwind Petrov-Galerkin(CSUPG) discretization of the redistancing equation gives satisfactory solutions for two test problems while the solutions of streamline upwind Petrov-Galerkin(SUPG) discretization are dissipated by the numerical diffusion added for the stability of a hyperbolic system. Furthermore, it has been found that the solutions obtained by CSUPG method are comparable to those by second order ENO method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.8 s.179
• /
• pp.2091-2099
• /
• 2000

The object of the present study is to present an adaptive wavelet-Galerkin method for the analysis of thin-walled box beam. Due to good localization properties of wavelets, wavelet methods emerge as alternative efficient solution methods to finite element methods. Most structural applications of wavelets thus far are limited in fixed-scale, non-adaptive frameworks, but this is not an appropriate use of wavelets. On the other hand, the present work appears the first attempt of an adaptive wavelet-based Galerkin method in structural problems. To handle boundary conditions, a fictitous domain method with penalty terms is employed. The limitation of the fictitious domain method is also addressed.

• 한국전산유체공학회:학술대회논문집
• /
• 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.50 no.3
• /
• pp.963-982
• /
• 2013

The Galerkin method has been developed mainly for 2nd order differential equations. To get numerical solutions, there are some choices of Riesz bases for the approximation subspace $V_j{\subset}L^2$. In this paper we shall propose a uniform approach to find suitable Riesz bases for higher order differential equations. Especially for the beam equation (4-th order equation), we also report numerical results.