• Title/Summary/Keyword: Galerkin methods

Search Result 132, Processing Time 0.027 seconds

PRECONDITIONED ITERATIVE METHOD FOR PETROV-GALERKIN PROCEDURE

  • Chung, Seiyoung;Oh, Seyoung
    • Journal of the Chungcheong Mathematical Society
    • /
    • v.10 no.1
    • /
    • pp.57-70
    • /
    • 1997
  • In this paper two preconditioned GMRES and QMR methods are applied to the non-Hermitian system from the Petrov-Galerkin procedure for the Poisson equation and compared to each other. To our purpose the ILUT and the SSOR preconditioners are used.

  • PDF

Optimal Test Function Petrov-Galerkin Method (최적시행함수 Petrov-Galerkin 방법)

  • Sung-Uk Choi
    • Journal of Korea Water Resources Association
    • /
    • v.31 no.5
    • /
    • pp.599-612
    • /
    • 1998
  • Numerical analysis of convection-dominated transport problems are challenging because of dual characteristics of the governing equation. In the finite element method, a strategy is to modify the test function to weight more in the upwind direction. This is called as the Petrov-Galerkin method. In this paper, both N+1 and N+2 Petrov-Galerkin methods are applied to transport problems at high grid Peclet number. Frequency fitting algorithm is used to obtain optimal levels of N+2 upwinding, and the results are discussed. Also, a new Petrov-Galerkin method, named as "Optimal Test Function Petrov-Galerkin Method," is proposed in this paper. The test function of this numerical method changes its shape depending upon relative strength of the convection to the diffusion. A numerical experiment is carried out to demonstrate the performance of the proposed method.

  • PDF

hp-DISCONTINUOUS GALERKIN METHODS FOR THE LOTKA-MCKENDRICK EQUATION: A NUMERICAL STUDY

  • Jeong, Shin-Ja;Kim, Mi-Young;Selenge, Tsendanysh
    • Communications of the Korean Mathematical Society
    • /
    • v.22 no.4
    • /
    • pp.623-640
    • /
    • 2007
  • The Lotka-McKendrick model which describes the evolution of a single population is developed from the well known Malthus model. In this paper, we introduce the Lotka-McKendrick model. We approximate the solution to the model using hp-discontinuous Galerkin finite element method. The numerical results show that the presented hp-discontinuous Galerkin method is very efficient in case that the solution has a sharp decay.

An Analysis of Transmission Line Structure by Combining Image Mode and Galerkin Methods (영상-모드 및 Galerkin법을 이용한 전송선 문제 해석)

  • 신규현;정현교;한송엽
    • The Transactions of the Korean Institute of Electrical Engineers
    • /
    • v.40 no.12
    • /
    • pp.1296-1301
    • /
    • 1991
  • This paper presents a hybrid image mode Galerkin method for the analysis of the transmission line structures suspended between infinite parallel ground planes. A Green's function that consists of numerically accelerated image mode terms is developed, which is used in boundary integral equation. Transmission lines of arbitrary cross section are analyzed using Galerkin's method. Two kinds of configurations of transmission lines are studied in sample problems.

  • PDF

Time-discontinuous Galerkin quadrature element methods for structural dynamics

  • Minmao, Liao;Yupeng, Wang
    • 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.

HIGHER ORDER FULLY DISCRETE SCHEME COMBINED WITH $H^1$-GALERKIN MIXED FINITE ELEMENT METHOD FOR SEMILINEAR REACTION-DIFFUSION EQUATIONS

  • S. Arul Veda Manickam;Moudgalya, Nannan-K.;Pani, Amiya-K.
    • Journal of applied mathematics & informatics
    • /
    • v.15 no.1_2
    • /
    • pp.1-28
    • /
    • 2004
  • We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by an $H^1$-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index one. A priori error estimates for semidiscrete scheme are derived for both differ-ential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.

Analysis of a strip footing on a homogenous soil using element free Galerkin method

  • Ganaiea, Aashiq H.;Sawant, Vishwas A.
    • Coupled systems mechanics
    • /
    • v.4 no.4
    • /
    • pp.365-383
    • /
    • 2015
  • Strip footing is an important type of shallow foundations and is commonly used beneath the walls. Analysis of shallow foundation involves the determination of stresses and deformations. Element free Galerkin method, one of the important mesh free methods, is used for the determination of stresses and deformations. Element free Galerkin method is an efficient and accurate method as compared to finite element method. The Element Free Galerkin method uses only a set of nodes and a description of model boundary is required to generate the discrete equation. Strip footing of width 2 m subjected to a loading intensity of 200 kPa is studied. The results obtained are agreeing with the values obtained using analytical solutions available in the literature. Parametric study is done and the effect of modulus of deformation, Poisson's ratio and scaling parameter on deformation and stresses are determined.

REVIEW AND IMPLEMENTATION OF STAGGERED DG METHODS ON POLYGONAL MESHES

  • KIM, DOHYUN;ZHAO, LINA;PARK, EUN-JAE
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.25 no.3
    • /
    • pp.66-81
    • /
    • 2021
  • In this paper, we review the lowest order staggered discontinuous Galerkin methods on polygonal meshes in 2D. The proposed method offers many desirable features including easy implementation, geometrical flexibility, robustness with respect to mesh distortion and low degrees of freedom. Discrete function spaces for locally H1 and H(div) spaces are considered. We introduce special properties of a sub-mesh from a given star-shaped polygonal mesh which can be utilized in the construction of discrete spaces and implementation of the staggered discontinuous Galerkin method. For demonstration purposes, we consider the lowest case for the Poisson equation. We emphasize its efficient computational implementation using only geometrical properties of the underlying mesh.

NUMERICAL METHODS FOR A STIFF PROBLEM ARISING FROM POPULATION DYNAMICS

  • Kim, Mi-Young
    • 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.

  • PDF