• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제13권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.

• Journal of Mechanical Science and Technology
• /
• 제17권7호
• /
• pp.986-998
• /
• 2003

In this paper, an adaptive numerical integration scheme, which does not need non-overlapping and contiguous integration meshes, is proposed for the MLPG (Meshless Local Petrov-Galerkin) method. In the proposed algorithm, the integration points are located between the neighboring nodes to properly consider the irregular nodal distribution, and the nodal points are also included as integration points. For numerical integration without well-defined meshes, the Shepard shape function is adopted to approximate the integrand in the local symmetric weak form, by the values of the integrand at the integration points. This procedure makes it possible to integrate the local symmetric weak form without any integration meshes (non-overlapping and contiguous integration domains). The convergence tests are performed, to investigate the present scheme and several numerical examples are analyzed by using the proposed scheme.

• Journal of applied mathematics & informatics
• /
• 제29권1_2호
• /
• pp.311-326
• /
• 2011

In this paper, we adopt discontinuous Galerkin methods with penalty terms namely symmetric interior penalty Galerkin methods, to solve nonlinear viscoelasticity type equations. We construct finite element spaces and define an appropriate projection of u and prove its optimal convergence. We construct extrapolated fully discrete discontinuous Galerkin approximations for the viscoelasticity type equation and prove ${\ell}^{\infty}(L^2)$ optimal error estimates in both spatial direction and temporal direction.

• Korean Journal of Mathematics
• /
• 제31권4호
• /
• pp.419-431
• /
• 2023

In this paper, we propose an iterative method for a symmetric interior penalty Galerkin method for heterogeneous elliptic problems. The iterative method consists mainly of two parts based on a non-overlapping domain decomposition approach. One is an intermediate preconditioner constructed by understanding the properties of the discontinuous finite element functions and the other is a preconditioning related to the dual-primal finite element tearing and interconnecting (FETI-DP) methodology. Numerical results for the proposed method are presented, which demonstrate the performance of the iterative method in terms of various parameters associated with the elliptic model problem, the finite element discretization, and non-overlapping subdomain decomposition.

• 대한수학회지
• /
• 제51권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
• /
• 제13권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.

• Structural Engineering and Mechanics
• /
• 제35권4호
• /
• pp.387-407
• /
• 2010

In this paper, using perturbation and Galerkin method, the response of a resonant viscoelastic microbeam to an electric actuation is obtained. The microbeam is under axial load and electrical load. It is assumed that midplane is stretched, when the beam is deflected. The equation of motion is derived using the Newton's second law. The viscoelastic model is taken to be the Kelvin-Voigt model. In the first section, the static deflection is obtained using the Galerkin method. Exact linear symmetric mode shape of a straight beam and its deflection function under constant transverse load are used as admissible functions. So, an analytical expression that describes the static deflection at all points is obtained. Comparing the result with previous research show that using deflection function as admissible function decreases the computation errors and previous calculations volume. In the second section, the response of a microbeam resonator system under primary and secondary resonance excitation has been obtained by analytical multiple scale perturbation method combined with the Galerkin method. It is shown, that a small amount of viscoelastic damping has an important effect and causes to decrease the maximum amplitude of response, and to shift the resonance frequency. Also, it shown, that an increase of the DC voltage, ratio of the air gap to the microbeam thickness, tensile axial load, would increase the effect of viscoelastic damping, and an increase of the compressive axial load would decrease the effect of viscoelastic damping.

• 대한기계학회논문집A
• /
• 제34권5호
• /
• pp.629-635
• /
• 2010

SGBEM-FEM 교호법이 Nikishkov, Park 및 Atluri 에 의하여 제안되었었다. 제안된 방법을 사용하면 임의 형태의 평면 혹은 비평면 삼차원 균열에 대하여 복합 모드의 응력강도계수를 구할 수 있다. 그러나 현장에서의 적용을 위해서는 이 방법의 정확성 및 신뢰성에 대한 검토가 더욱 필요하다. 따라서 본 논문에서는 응력강도계수에 영향을 주는 주요한 몇 가지 인자를 검토하였다. 그리고 원통의 내부 및 외부에 존재하는 원주방향 표면균열에 대한 응력강도계수를 구하여 기존의 해와 비교하였다. 그 결과 SGBEM-FEM 교호법은 이들 균열에 대하여 정확한 해를 주고 있음을 확인하였다.

• 대한기계학회논문집A
• /
• 제33권2호
• /
• pp.145-152
• /
• 2009

SGBEM(Symmetric Galerkin Boundary Element Method)-FEM alternating method has been proposed by Nikishkov, Park and Atluri. In the proposed method, arbitrarily shaped three-dimensional crack problems can be solved by alternating between the crack solution in an infinite body and the finite element solution without a crack. In the previous study, the SGBEM-FEM alternating method was extended further in order to solve elastic-plastic crack problems and to obtain elastic-plastic stress fields. For the elastic-plastic analysis the algorithm developed by Nikishkov et al. is used after modification. In the algorithm, the initial stress method is used to obtain elastic-plastic stress and strain fields. In this paper, elastic-plastic J integrals for three-dimensional cracks are obtained using the method. For that purpose, accurate values of displacement gradients and stresses are necessary on an integration path. In order to improve the accuracy of stress near crack surfaces, coordinate transformation and partitioning of integration domain are used. The coordinate transformation produces a transformation Jacobian, which cancels the singularity of the integrand. Using the developed program, simple three-dimensional crack problems are solved and elastic and elastic-plastic J integrals are obtained. The obtained J integrals are compared with the values obtained using a handbook solution. It is noted that J integrals obtained from the alternating method are close to the values from the handbook.

• Steel and Composite Structures
• /
• 제32권2호
• /
• pp.253-264
• /
• 2019

This paper presents the combination resonances of FG porous (FGP) cylindrical shell under two-term excitation. The effect of structural damping on the system response is also considered. With regard to classical plate theory of shells, von-$K{\acute{a}}rm{\acute{a}}n$ equation and Hook law, the relations of stress-strain is derived for shell. According to the Galerkin method, the discretized motion equation is obtained. The combination resonances are obtained by using the method of multiple scales. Four types of FGP distributions consist of uniform porosity, non-symmetric porosity soft, non-symmetric porosity stiff and symmetric porosity distribution are considered. The influence of various porosity distributions, porosity coefficients of cylindrical shell and amplitude excitations on the combination resonances for FGP cylindrical shells is investigated.