• International Journal of Fuzzy Logic and Intelligent Systems
• /
• 제16권1호
• /
• pp.64-71
• /
• 2016

This paper proposes nonlinear behavior in a love model for Romeo and Juliet with an external force of discontinuous time. We investigated the periodic motion and chaotic behavior in the love model by using time series and phase portraits with respect to some variable and fixed parameters. The computer simulation results confirmed that the proposed love model with an external force of discontinuous time shows periodic motion and chaotic behavior with respect to parameter variation.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2008년도 춘계학술대회 논문집
• /
• pp.231-234
• /
• 2008

This paper is concerned with the analysis of elasto-plastic stress waves by a time discontinuous variational integrator based on Hamiltonian in order to more accurate results in one dimensional dynamic problem. The proposed algorithm adopts both time-discontinuous variational integrator and space-continuous Hamiltonian so as to capture discontinuities of stress waves. This study enables to preserve total mechanical energy such as internal energy, kinetic energy and dissipative energy due to plastic deformation for long integration time. Finite element analysis of elasto-plastic stress waves is carried out in order to demonstrate the accuracy of the proposed algorithm.

• Structural Engineering and Mechanics
• /
• 제85권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.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2008년도 추계학술대회 논문집
• /
• pp.263-266
• /
• 2008

This paper is concerned with the analysis of elasto-plastic stress waves in a mode I semi-infinite cracked solid subjected to Heaviside pulse load. This study adopts a time-discontinuous variational integrator based on Hamiltonian in order to reduce the numerical dispersive and dissipative errors. This also utilizes an integration scheme of the constitutive model with 2nd-order accuracy which is formulated on the strain space for a rate and temperature dependent material model. Finite element analyses of elasto-plastic stress waves are carried out in order to compare the accuracy between a conventional Galerkin method and the time- discontinuous variational integrator.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제23권4호
• /
• pp.301-327
• /
• 2019

The proper orthogonal decomposition (POD) method for time-dependent problems significantly reduces the computational time as it reduces the original problem to the lower dimensional space. Even a higher degree of reduction can be reached if the solution is smooth in space and time. However, if the solution is discontinuous and the discontinuity is parameterized e.g. with time, the POD approximations are not accurate in the reduced space due to the lack of ability to represent the discontinuous solution as a finite linear combination of smooth bases. In this paper, we propose to post-process the sample solutions and re-initialize the POD approximations to deal with discontinuous solutions and provide accurate approximations while the computational time is reduced. For the post-processing, we use the Gegenbauer reconstruction method. Then we regularize the Gegenbauer reconstruction for the construction of POD bases. With the constructed POD bases, we solve the given PDE in the reduced space. For the POD approximation, we re-initialize the POD solution so that the post-processed sample solution is used as the initial condition at each sampling time. As a proof-of-concept, we solve both one-dimensional linear and nonlinear hyperbolic problems. The numerical results show that the proposed method is efficient and accurate.

• 전산구조공학
• /
• 제6권4호
• /
• pp.83-88
• /
• 1993

시간변수에 대하여 불연속성을 주는 시간불연속 Galerkin 방법을 유한요소법으로 해석하였다. 이 방법은 미분방정식 관점에서 지금까지 요소간에 연속성을 준 일반적 유한요소법과 다르게 임의의 시간요소를 선택, 매 시간단계에서 요소경계에 불연속을 허락함으로서 해의 정확성을 높이고 무조건의 안정을 주는 상미분방정식의 해법인 것이다.

• 한국전산유체공학회지
• /
• 제16권4호
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

• 대한수학회보
• /
• 제47권1호
• /
• pp.89-102
• /
• 2010

This paper concerns the problem of global robust stability of a time-delay discontinuous system with a positive-defined connection matrix under polytopic-type uncertainty. In order to give the stability condition, we firstly address the existence of solution and equilibrium point based on the properties of M-matrix, Lyapunov-like approach and the theories of differential equations with discontinuous right-hand side as introduced by Filippov. Second, we give the delay-independent and delay-dependent stability condition in terms of linear matrix inequalities (LMIs), and based on Lyapunov function and the properties of the convex sets. One numerical example demonstrate the validity of the proposed criteria.

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

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제24권2호
• /
• pp.143-159
• /
• 2020

We propose a consistent numerical method for elliptic interface problems with nonhomogeneous jumps. We modify the discontinuous bubble immersed finite element method (DB-IFEM) introduced in (Chang et al. 2011), by adding a consistency term to the bilinear form. We prove optimal error estimates in L² and energy like norm for this new scheme. One of the important technique in this proof is the Bramble-Hilbert type of interpolation error estimate for discontinuous functions. We believe this is a first time to deal with interpolation error estimate for discontinuous functions. Numerical examples with various interfaces are provided. We observe optimal convergence rates for all the examples, while the performance of early DB-IFEM deteriorates for some examples. Thus, the modification of the bilinear form is meaningful to enhance the performance.