• International Journal of Naval Architecture and Ocean Engineering
• /
• 제11권1호
• /
• pp.572-583
• /
• 2019

This paper aims to assess the applicability of the Runge Kutta Discontinuous Galerkin-Direct Ghost Fluid Method to the internal explosion inside a water-filled tube, which previously was studied by many researchers in separate works. Once the explosive charge located at the inner center of the water-filled tube explodes, the tube wall is subjected to an extremely high intensity fluid loading and deformed. The deformation causes a modification of the field of fluid flow in the region near the water-structure interface so that has substantial influence on the response of the structure. To connect the structure and the fluid, valid data exchanges along the interface are essential. Classical fluid structure interaction simulations usually employ a matched meshing scheme which discretizes the fluid and structure domains using a single mesh density. The computational cost of fluid structure interaction simulations is usually governed by the structure because the size of time step may be determined by the density of structure mesh. The finer mesh density, the better solution, but more expensive computational cost. To reduce such computational cost, a non-matched meshing scheme which allows for different mesh densities is employed. The coupled numerical approach of this paper has fewer difficulties in the implementation and computation, compared to gas dynamics based approach which requires complicated analytical manipulations. It can also be applied to wider compressible, inviscid fluid flow analyses often found in underwater explosion events.

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

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

• 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.

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

• Journal of applied mathematics & informatics
• /
• 제22권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.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2005년도 춘계학술대회 논문집
• /
• pp.43-46
• /
• 2005

Micro molding technology is a promising mass production technology for polymer based microstructures. Mass production technologies such as the micro injection/compression molding, hot embossing, and micro reaction molding are already in use. In the present study, we have developed a numerical analysis system to simulate three-dimensional non-isothermal cavity filling for hot embossing, with a special emphasis on the free surface capturing. Precise free surface capturing has been successfully accomplished with the level set method, which is solved by means of the Runge-Kutta discontinuous Galerkin (RKDG) method. The RKDG method turns out to be excellent from the viewpoint of both numerical stability and accuracy of volume conservation. The Stokes equations are solved by the stabilized finite element method using the equal order tri-linear interpolation function. To prevent possible numerical oscillation in temperature Held we employ the streamline upwind Petrov-Galerkin (SUPG) method. With the developed code we investigated the detailed change of free surface shape in time during the mold filling. In the filling simulation of a simple rectangular cavity with repeating protruded parts, we find out that filling patterns are significantly influenced by the geometric characteristics such as the thickness of base plate and the aspect ratio and pitch of repeating microstructures. The numerical analysis system enables us to understand the basic flow and material deformation taking place during the cavity filling stage in microstructure fabrications.

• Advances in aircraft and spacecraft science
• /
• 제7권5호
• /
• pp.387-404
• /
• 2020

The numerical simulation of shock waves in supersonic flows is challenging because of several instabilities which can affect the solution. Among them, the carbuncle phenomenon can introduce nonphysical perturbations in captured shock waves. In the present work, a hybrid numerical flux is proposed for the evaluation of the convective fluxes that avoids carbuncle and keeps high-accuracy on shocks and boundary layers. In particular, the proposed flux is a combination between an upwind approximate Riemann problem solver and the Local Lax-Friedrichs scheme. A simple strategy to mix the two fluxes is proposed and tested in the framework of a discontinuous Galerkin discretisation. The approach is investigated on the subsonic flow in a channel, on the supersonic flow around a cylinder, on the supersonic flow on a flat plate and on the flow in a overexpanded rocket nozzle.

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

• 한국시뮬레이션학회:학술대회논문집
• /
• 한국시뮬레이션학회 2001년도 The Seoul International Simulation Conference
• /
• pp.373-379
• /
• 2001

There is a big issue to generate 3D hexahedral finite element (FE) model, since a process to divide the whole domain into several simple-shaped sub-domains is required before generating a continuous mesh with mapped mesh generators. In general, it is nearly impossible to set up proper division numbers interactively to keep mesh connectivity between sub-domains on a complicated arbitrary-shaped domain. If mesh continuity between sub-domains is not required in an analysis, this complicated process can be omitted. Element-free Galerkin method (EFGM) can accept discontinuous meshes, which only requires nodal information. However it is difficult to choose a reasonable influenced domain in moving least squares scheme with non-uniformly distributed nodes in discontinuous FE models. A new FE scheme fur discontinuous mesh is proposed in this paper by applying improved EFGM with some modification to derive FE approximated function in discontinuous parts. Its validity is evaluated on linear elastic problems.

• 대한토목학회논문집
• /
• 제34권5호
• /
• pp.1383-1393
• /
• 2014

최근, 급속한 컴퓨터 하드웨어의 성능 향상과 전산유체역학 분야의 이론적 발전으로, 고차 정확도의 수치기법들이 계산수리학 분야에 적용되어 왔다. 본 연구에서는 1차원 천수방정식에 대한 수치 해법으로 TVD Runge-Kutta 불연속 갤러킨(RKDG) 유한요소법을 적용하였다. 대표적인 천이류(transcritical flow)의 예로 순간적인 댐 붕괴에 의한 댐 붕괴류(dam-break flow) 흐름과 지형변화에 의한 천이류를 모의하였다. 리만(Riemann) 근사해법으로 로컬 Lax-Friedrichs (LLF), Roe, HLL 흐름률(flux) 기법을 사용하였고, 불필요한 진동을 제거하기 위하여, 기울기 제한자로서 MUSCL 제한자를 사용하였다. 개발된 모델은 1차원 댐 붕괴류와 천이류에 적용하였다. 수치해석 결과는 해석해, 수리실험 결과와 비교하였다.