• Korean Journal of Mathematics
• /
• v.13 no.2
• /
• pp.161-176
• /
• 2005

We consider a model of population dynamics whose mortality function is unbounded. We note that the regularity of the solution depends on the growth rate of the mortality near the maximum age. We propose Gauss-Legendre methods along the characteristics to approximate the solution when the solution is smooth enough. It is proven that the scheme is convergent at fourth-order rate in the maximum norm. We also propose discontinuous Galerkin finite element methods to approximate the solution which is not smooth enough. The stability of the method is discussed. Several numerical examples are presented.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.1
• /
• pp.39-78
• /
• 2020

In this paper, we propose a decoupled local discontinuous Galerkin method for solving the Klein-Gordon-Schrödinger (KGS) equations. The KGS equations is a model of the Yukawa interaction of complex scalar nucleons and real scalar mesons. The advantage of our scheme is that the computation of the nucleon and meson field is fully decoupled, so that it is especially suitable for parallel computing. We present the conservation property of our fully discrete scheme, including the energy and Hamiltonian conservation, and establish the optimal error estimate.

• Journal of applied mathematics & informatics
• /
• v.18 no.1_2
• /
• pp.113-126
• /
• 2005

The Cahn-Hilliard equation is modeled to describe the dynamics of phase separation in glass and polymer systems. A priori error estimates for the Cahn-Hilliard equation have been studied by the authors. In order to control accuracy of approximate solutions, a posteriori error estimation of the Cahn-Hilliard equation is obtained by discontinuous Galerkin method.

• The Journal of the Korea Contents Association
• /
• v.13 no.3
• /
• pp.428-434
• /
• 2013

With the increase of the frequency in large-scale floods and natural disasters, the demands for highly accurate numerical river models are also rapidly growing. Generally, flows in rivers are modeled with previously developed and well-established numerical models based on shallow water equations. However, the so-far-developed models reveal a lot of limitations in the analysis of discontinuous flow or flow which needs accurate modeling. In this study, the numerical shallow water model based on the discontinuous Galerkin method was applied to the simulation of one-dimensional transcritical flow, including dam break flows and a flow over a hump. The favorable agreement was observed between numerical solutions and analytical solutions.

• East Asian mathematical journal
• /
• v.31 no.5
• /
• pp.685-693
• /
• 2015

In this paper, we introduce fully discrete mixed discontinuous Galerkin approximations for parabolic problems. And we analyze the error estimates in $l^{\infty}(L^2)$ norm for the primary variable and the error estimates in the energy norm for the primary variable and the flux variable.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.4
• /
• pp.295-308
• /
• 2016

We propose a projection-based analysis of a new hybridizable discontinuous Gale-rkin method for second order elliptic equations. The method is more advantageous than the standard HDG method in a sense that the new method has higher-order accuracy and lower computational cost, and is more flexible. Notable distinctions of our new method, when compared to the standard HDG emthod, are that our method uses $L^2$-projection and suitable stabilization parameter depending on a mesh size for superconvergence. We show that the error for the solution of the equation converges with order p + 2 when we only use polynomials of degree p + 1 as a finite element space without postprocessing. After establishing the theory, we carry out numerical tests to demonstrate and ensure that the proposed method is effective and accurate in practice.

• Advances in aircraft and spacecraft science
• /
• v.4 no.5
• /
• pp.555-571
• /
• 2017

The description of transitional flows by means of RANS equations is sometimes based on non-local approaches which require the computation of some boundary layer properties. In this work a non-local Laminar Kinetic Energy model is used to predict transitional and separated flows. Usually the non-local term of this model is evaluated along the grid lines of a structured mesh. An alternative approach, which does not rely on grid lines, is introduced in the present work. This new approach allows the use of fully unstructured meshes. Furthermore, it reduces the grid-dependence of the predicted results. The approach is employed to study the transitional flows in the T106c turbine cascade and around a NACA0021 airfoil by means of a discontinuous Galerkin method. The local nature of the discontinuous Galerkin reconstruction is exploited to implement an adaptive algorithm which automatically refines the mesh in the most significant regions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.4
• /
• pp.267-280
• /
• 2009

In this paper we investigate a simple two-level additive Schwarz preconditioner for the P1 symmetric interior penalty Galerkin method of the Poisson equation on rectangular meshes. The construction is based on the decomposition of the global space of piecewise linear polynomials into the sum of local subspaces, each of which corresponds to an element of the underlying mesh, and the global coarse subspace consisting of piecewise constants. This preconditioner is a direct combination of the block Jacobi iteration and the cell-centered finite difference method, and thus very easy to implement. Explicit upper and lower bounds for the maximum and minimum eigenvalues of the preconditioned matrix system are derived and confirmed by some numerical experiments.

• Journal of applied mathematics & informatics
• /
• v.22 no.1_2
• /
• pp.317-329
• /
• 2006

In this paper, a singularly perturbed ordinary differential equation with non-smooth data is considered. The numerical method is generated by means of a Petrov-Galerkin finite element method with the piecewise-exponential test function and the piecewise-linear trial function. At the discontinuous point of the coefficient, a special technique is used. The method is shown to be first-order accurate and singular perturbation parameter uniform convergence. Finally, numerical results are presented, which are in agreement with theoretical results.

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.1-11
• /
• 2016

Considerable efforts are required to develop a monotone, robust and stable high-order numerical scheme for solving the hyperbolic system. The discontinuous Galerkin(DG) method is a natural choice, but elimination of the spurious oscillations from the high-order solutions demands a new development of proper limiters for the DG method. There are several available limiters for controlling or removing unphysical oscillations from the high-order approximate solution; however, very few studies were directed to analyze the exact role of the limiters in the hyperbolic systems. In this study, the performance of the several well-known limiters is examined by comparing the high-order($p^1$, $p^2$, and $p^3$) approximate solutions with the exact solutions. It is shown that the accuracy of the limiter is in general problem-dependent, although the Hermite WENO limiter and maximum principle limiter perform better than the TVD and generalized moment limiters for most of the test cases. It is also shown that application of the troubled cell indicators may improve the accuracy of the limiters under some specific conditions.