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

• 한국전산유체공학회지
• /
• 제21권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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.311-317
• /
• 2011

The present paper deals with the continuous works of extending the multi-dimensional limiting process (MLP) for compressible flows, which has been quite successful in finite volume methods, into discontinuous Galerkin (DG) methods. From the series of the previous, it was observed that the MLP shows several superior characteristics, such as an efficient controlling of multi-dimensional oscillations and accurate capturing of both discontinuous and continuous flow features. Mathematically, fundamental mechanism of oscillation-control in multiple dimensions has been established by satisfaction of the maximum principle. The MLP limiting strategy is extended into DG framework, which takes advantage of higher-order reconstruction within compact stencil, to capture detailed flow structures very accurately. At the present, it is observed that the proposed approach yields outstanding performances in resolving non-compressive as well as compressive flaw features. In the presentation, further numerical analyses and results are going to be presented to validate that the newly developed DG-MLP methods provide quite desirable performances in controlling numerical oscillations as well as capturing key flow features.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제25권4호
• /
• pp.173-195
• /
• 2021

A new implicit discontinuous Galerkin spectral element method (DGSEM) based on the first order hyperbolic system(FOHS) is presented for solving elliptic type partial different equations, such as the Poisson problems. By utilizing the idea of hyperbolic formulation of Nishikawa[1], the original Poisson equation was reformulated in the first-order hyperbolic system. Such hyperbolic system is solved implicitly by the collocation type DGSEM. The steady state solution in pseudo-time, which is the solution of the original Poisson problem, was obtained by the implicit solution of the global linear system. The optimal polynomial orders of 𝒪(𝒽^𝑝+1)) are obtained for both the solution and gradient variables from the test cases in 1D and 2D regular grids. Spectral accuracy of the solution and gradient variables are confirmed from all test cases of using the uniform grids in 2D.

• 한국해안·해양공학회논문집
• /
• 제32권6호
• /
• pp.569-574
• /
• 2020

천수방정식에 대한 불연속 갤러킨 기법 (DG)은 주로 양해법 기반으로 개발되어 적용되어 왔으나, 바닥마찰항의 처리, 과도한 CFL 조건 등의 불리한 점이 지적되어 왔다. 이에 대한 대안으로써, 본 연구에서는 음해법 기반의 모형을 개발하고 이를 적용하여 향후 가능성을 입증하였다. 본 논문에서 연구한 사례에서는 선형 삼각형 요소를 사용하였고, 수치흐름률로서 Roe 흐름률을 이용하였으며, TVD 특성 보존을 위한 기울기 제한자를 적용하였다. 적용 사례로서 실린더 주변의 흐름과 댐 붕괴류 문제 등에 대하여 적용하고, 기존의 실험치, 수치해와 비교하여 잘 일치함을 확인하였다.

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

• 한국전산유체공학회지
• /
• 제12권4호
• /
• pp.14-19
• /
• 2007

We present a direct numerical simulation technique and some preliminary results of the capillary spreading of a droplet containing particles on the solid substrate. We used the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension and employed the discontinuous Galerkin method for the stabilization of the interface advection equation. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles. We investigated the droplet spreading by the capillary force and discussed effects of the presence of particles on the spreading behavior. It has been observed that a particulate drop spreads less than the pure liquid drop. The amount of spread of a particulate drop has been found smaller than that of the liquid with effectively the same viscosity as the particulate drop.

• 한국지반환경공학회 논문집
• /
• 제15권3호
• /
• pp.41-48
• /
• 2014

준설토에 대한 연구는 주로 준설토의 1차원 침강 및 자중압밀 특성을 파악하는 실험적 연구가 진행되었다. 하지만 양질의 준설지반 확보를 위한 효과적인 투기장의 설계와 배출수에 의한 환경오염을 최소화하기 위해서는 준설토의 투기에 의한 유동특성의 체계적인 연구가 필요하다. 본 연구에서는 준설토 투기장의 펌핑에 의한 토사의 유동 형상을 모사하기 위하여 준설토사를 단일상으로 가정하고 연속 방정식을 유도하여 좌표축에 따른 힘 평형 방정식을 유도하였다. 준설토장의 3차원 거동 해석을 위한 컴퓨터 연산 부하와 모델링 소요시간을 최적화하기 위하여, 토체의 깊이 방향으로 적분을 수행하는 깊이 적분 방법을 지배 방정식에 적용하여, 3차원적 지형조건을 고려할 수 있도록 하였다. 지배 방정식의 보간함수를 이용한 공간분할에서 Petrov-Galerkin 수식화 기법을 적용하였다. 일반화된 사다리꼴 법칙으로 시간적분을 수행하고 Newton의 반복과정을 이용할 수 있도록 근사화시켰다. 가중행렬은 DG과 CDG 기법을 적용하였으며, 준설토 유동해석에서 가중행렬에 따른 수치적인 안정성을 평가하기 위하여 사각형 기둥 슬럼프 시뮬레이션을 수행하였다. 수치기법에 대한 비교 분석 결과는 DG 기법을 적용한 SU/PG 수식화가 유사진동을 최소화시키는 가장 안정적인 수치해석결과를 보여주는 것으로 나타났다.

• 한국전산유체공학회지
• /
• 제19권3호
• /
• pp.52-61
• /
• 2014

The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

• Advances in aircraft and spacecraft science
• /
• 제9권5호
• /
• pp.463-477
• /
• 2022

This study compares various ways of calculating flows for the problems with the presence of shock waves by first-order schemes and higher-order DG method on the tests from the Quirk list, namely: Quirk's problem and its modifications, shock wave diffraction at a 90 degree corner, the problem of double Mach reflection. It is shown that the use of HLLC and Godunov's numerical schemes flows in calculations can lead to instability, the Rusanov-Lax-Friedrichs scheme flow can lead to high dissipation of the solution. The most universal in heavy production calculations are hybrid schemes flows, which allow the suppression of the development of instability and conserve the accuracy of the method.