• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.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.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.6
• /
• pp.569-574
• /
• 2020

Though the discontinuous Galerkin (DG) method has been developed and applied to shallow water equations mainly in explicit schemes, they have been criticized for the limitation in treatment of bottom friction terms and severe CFL conditions. In this study, an implicit scheme is devised and applied to some representative benchmark problems. The linear triangular elements were employed and the Roe numerical fluxes were adopted for convective fluxes. To preserve TVD property, the slope limiter was employed. As the case studies, the model is applied to the flow around the cylinders and the dam-break flow. Then, the results are compared with the experimental and numerical data of previous studies and good agreements were observed.

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

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

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.05a
• /
• 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
• /
• v.7 no.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.

• Journal of computational fluids engineering
• /
• v.12 no.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.

• Proceedings of the Korea Society for Simulation Conference
• /
• 2001.10a
• /
• 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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1383-1393
• /
• 2014

Recently, with rapid improvement in computer hardware and theoretical development in the field of computational fluid dynamics, high-order accurate schemes also have been applied in the realm of computational hydraulics. In this study, numerical solutions of 1D shallow water equations are presented with TVD Runge-Kutta discontinuous Galerkin (RKDG) finite element method. The transcritical flows such as dam-break flows due to instant dam failure and transcritical flow with bottom elevation change were studied. As a formulation of approximate Riemann solver, the local Lax-Friedrichs (LLF), Roe, HLL flux schemes were employed and MUSCL slope limiter was used to eliminate unnecessary numerical oscillations. The developed model was applied to 1D dam break and transcritical flow. The results were compared to the exact solutions and experimental data.