• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.10a
• /
• pp.11-18
• /
• 2001

In this paper, a new algorithm of coupling Element-Free Galerkin Method(EFGM) and Boundary Element Method(BEM) using the variational formulation is presented. A global variational coupling formulation of EFGM-BEM is achieved by combining the variational form on each subregion. In the formulation, Lagrange multiplier method is introduced to satisfy the compatibility conditions between EFGM subregion and BEM subregion. Some numerical examples are studied to verify accuracy and efficiency of the proposed method, in which numerical performance of the method is compared with that of conventional method such as EFGM-BEM direct coupling method, EFGM and BEM. The proposed method incorporating the merits of EFGM and BEM is expected to be applied to special engineering problems such as the crack propogation problems in very large domain, and underground structures with joints.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.507-512
• /
• 2007

A new type of extended element-free Galerkin method (XFEM) is proposed on this paper. The blending region which was inevitable in the extended finite element method and the extended meshfree method is removed in this method. For this end, two different techniques are developed. The first one is the modification of the domain of influence so that the crack tip is always placed on the edge of a domain of influence. The second method is the use of the Lagrange multiplier. The crack is virtually extended beyond the actual crack tip. The virtual extension was forced close by the Lagrange multiplier. The first method can be applied to two dimensional problems only Lagrange multiplier method can be used in both two and three dimensions.

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

• International Journal of Control, Automation, and Systems
• /
• v.4 no.6
• /
• pp.689-696
• /
• 2006

This paper presents a new algorithm for the closed-loop $H_{\infty}$ composite control of weakly coupled bilinear systems with time-varying parameter uncertainties and exogenous disturbance using the successive Galerkin approximation(SGA). By using weak coupling theory, the robust $H_{\infty}$ control can be obtained from two reduced-order robust $H_{\infty}$ control problems in parallel. The $H_{\infty}$ control theory guarantees robust closed-loop performance but the resulting problem is difficult to solve for uncertain bilinear systems. In order to overcome the difficulties inherent in the $H_{\infty}$ control problem, two $H_{\infty}$ control laws are constructed in terms of the approximated solution to two independent Hamilton-Jacobi-Isaac equations using the SGA method. One of the purposes of this paper is to design a closed-loop parallel robust $H_{\infty}$ control law for the weakly coupled bilinear systems with parameter uncertainties using the SGA method. The other is to reduce the computational complexity when the SGA method is applied to the high order systems.

• 제어로봇시스템학회:학술대회논문집
• /
• 2004.08a
• /
• pp.407-412
• /
• 2004

This paper presents a new algorithm for the closed-loop $H_{\infty}$ composite control of singularly perturbed nonlinear systems with a exogenous disturbance, using the successive Galerkin approximation(SGA). The singularly perturbed nonlinear system is decomposed into two subsystems of a slow-time scale and a fast-time scale via singular perturbation theory, and two $H_{\infty}$ control laws are obtained to each subsystem by using the SGA method. The composite control law that consists of two $H_{\infty}$ control laws of each subsystem is designed. One of the purposes of this paper is to design the closed-loop $H_{\infty}$ composite control law for the singularly perturbed nonlinear systems via the SGA method. The other is to reduce the computational complexity when the SGA method is applied to the high order systems.

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

• The Journal of the Korea Contents Association
• /
• v.13 no.3
• /
• pp.428-434
• /
• 2013

With the increase of the frequency in large-scale floods and natural disasters, the demands for highly accurate numerical river models are also rapidly growing. Generally, flows in rivers are modeled with previously developed and well-established numerical models based on shallow water equations. However, the so-far-developed models reveal a lot of limitations in the analysis of discontinuous flow or flow which needs accurate modeling. In this study, the numerical shallow water model based on the discontinuous Galerkin method was applied to the simulation of one-dimensional transcritical flow, including dam break flows and a flow over a hump. The favorable agreement was observed between numerical solutions and analytical solutions.

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

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.12 no.3
• /
• pp.133-141
• /
• 2007

A finite element model is developed for the extended Boussinesq equations that is capable of simulating the dynamics of long and short waves. Galerkin weighted residual method and the introduction of auxiliary variables for 3rd spatial derivative terms in the governing equations are used for the model development. The Adams-Bashforth-Moulton Predictor Corrector scheme is used as a time integration scheme for the extended Boussinesq finite element model so that the truncation error would not produce any non-physical dispersion or dissipation. This developed model is applied to the problems of solitary wave propagation. Predicted results is compared to available analytical solutions and laboratory measurements. A good agreement is observed.

• Structural Engineering and Mechanics
• /
• v.61 no.6
• /
• pp.765-773
• /
• 2017

The quadratic B-spline finite element method yields mass and stiffness matrices which are half the size of matrices obtained by the conventional finite element method. We solve the free vibration problem of a rotating Rayleigh beam using the quadratic B-spline finite element method. Rayleigh beam theory includes the rotary inertia effects in addition to the Euler-Bernoulli theory assumptions and presents a good mathematical model for rotating beams. Galerkin's approach is used to obtain the weak form which yields a system of symmetric matrices. Results obtained for the natural frequencies at different rotating speeds show an accurate match with the published results. A comparison with Euler-Bernoulli beam is done to decipher the variations in higher modes of the Rayleigh beam due to the slenderness ratio. The results are obtained for different values of non-uniform parameter ($\bar{n}$).