• Journal of information and communication convergence engineering
• /
• v.12 no.1
• /
• pp.39-45
• /
• 2014

A numerical method for solving Hamiltonian equations is said to be symplectic if it preserves the symplectic structure associated with the equations. Various symplectic methods are widely used in many fields of science and technology. A symplectic method preserves an approximate Hamiltonian perturbed from the original Hamiltonian. It theoretically supports the effectiveness of symplectic methods for long-term integration. Although it is also related to long-term integration, numerical stability of symplectic methods have received little attention. In this paper, we consider explicit symplectic methods defined for Hamiltonian equations with Hamiltonians of the special form, and study their numerical stability using the harmonic oscillator as a test equation. We propose a new stability criterion and clarify the stability of some existing methods that are visually based on the criterion. We also derive a new method that is better than the existing methods with respect to a Courant-Friedrichs-Lewy condition for hyperbolic equations; this new method is tested through a numerical experiment with a nonlinear wave equation.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.11b
• /
• pp.1249-1253
• /
• 2001

Numerical analysis is sometimes used to solve the problems in the engineering and natural science fields. On this reason, the faster, more practical system in computing the numerical solution is required. This paper deals with the numerical evaluation of various numerical integration methods which is frequently used in the engineering fields. This paper choices four integration methods such as Euler method, Heun's method, Runge-Kutta method and Gill's method for evaluating the each integration method. In numerical examples, the free vibration problem on an elastic foundation is chosen. As the numerical results, the natural frequencies and the running time are obtained, and these results are compared to examine the practicality of integration methods.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.63-66
• /
• 2002

The present study introduces some useful methods of implementing the Internet numerical analysis program with existing numerical codes for utilizing the Internet environment. The Internet gives developers good environment for development and release. Several methods were suggested, and some of them were implemented with an existing numerical code named SOLA-VOF, a computational fluid dynamics program to solve two-dimensional transient flow problems with free surface. This was reconstructed with Java technologies and compared with the original one. Java technologies have been applied to development of Internet applications for a long time. The objective of this work is to contrive methods of implementing Internet numerical analysis program with existing numerical codes and confirm the possibility of them. Methods using the applet-servlet communication were suggested and implemented. In addition, the Java web services with XML was introduced, which makes possible the cooperation of components. Although the concept has been suggested and developed for business applications, it can also be used for engineering softwares. Therefore, this study will be a preparation for numerical analysis to participate in engineering web services.

• KIEAE Journal
• /
• v.16 no.1
• /
• pp.81-87
• /
• 2016

Purpose: Heat transfers phenomena are described by the second order partial differential equation and its boundary conditions. In a three-dimensional structure of a building, the heat transfer phenomena generally include more than one material, and thus, become complicate. The analytic solutions are useful to understand heat transfer phenomena, but they can hardly be applied in engineering or design problems. Engineers and designers have generally been forced to use numerical methods providing reliable results. Finite volume methods with the unstructured grid system is only the suitable means of the analysis for the complex and arbitrary domains. Method: To obtain an numerical solution, a discretization method, which approximates the differential equations, and the interpolation methods for temperature and heat flux between two or more materials are required. The discretization methods are applied to small domains in space and time, and these numerical solutions form the descretized equations provide approximated solutions in both space and time. The accuracy of numerical solutions is dependent on the quality of discretizations and size of cells used. The higher accuracy, the higher numerical resources are required. The balance between the accuracy and difficulty of the numerical methods is critical for the success of the numerical analysis. A simple and easy interpolation methods among multiple materials are developed. The linear equations are solved with the BiCGSTAB being a effective matrix solver. Result: This study provides an overview of discretization methods, boundary interface, and matrix solver for the 3-dimensional numerical heat transfer including two materials.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.91-98
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• Bulletin of the Korean Space Science Society
• /
• 2008.10a
• /
• pp.26.4-27
• /
• 2008

The completion ('initiation' de facto) of the KASI Orbit Propagator and Estimator (KASIOPEA) has been delayed for several reasons unfortunately. Due to the lack of working staffs and the Division priority rearrangement, the initial plan was dismantled and ignored for many years. However, fundamental researches regarding the essential parts of KASIOPEA has been done by author. The numerical integration module of the KASIOPEA is the most sensitive part in the precision of the final output in general. There is no silver bullet in the numerical integration in an orbit propagation as a non-stiff ODE case. Many numerical integration method like single-step methods, multi-step method, and extrapolation methods have been used in overly populated orbit propagator or estimator. In this study, several popular methods from single-step, multi-step, and extrapolation methods have been tested in numerical accuracy and stability.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.110-117
• /
• 2004

The aim of a benchmark study is carried out nine methods are employed for ULS analysis which implicitly predict the ultimate strength of aluminium stiffened panels under axial compression. For this purpose, DNV PULS, experimental and numerical data on the ultimate strength of panels were collected. Comparison of these experimental / numerical, DNV PULS / numerical, results with theoretical solutions by the candidate methods is performed. Also it's compared that ALPS/ULSAP program is based on closed-form formula for the ULS of plates and grillages under axial compression. It is considered that ALPS/ULSAP methodology provides quite accurate and reasonable ULS calculations by a comparison with more refined FEA. Comparison of these experimental data, numerical, computational software results with the simplified solutions obtained by the candidate methods is then performed. The model uncertainties associated with the candidate methods are studied in terms of mean bias and COV (i.e., coefficient of variation) against experiments and numerical solutions, and the relative performance of the various methods is discussed.

• Fire Protection Technology
• /
• s.16
• /
• pp.20-23
• /
• 1994

This article investigates the different numerical methods, which are widely used for purpose of simulating a fire compartment the particular numerical methods such as finite difference, finite element, control Volume, and finite analysis are discribed in order to understand basic concepts and their applications. The fire simulations using fferent methods for the different physical geometrics have been reported in many recent literatures The convergence rate, the accuracy, and the stability are no simply dependent upon the specific method, The study of popular nu-merical methods by being compared among those is therefore significant to understand the nu-merical simulation of fire compartment.

• Journal of applied mathematics & informatics
• /
• v.33 no.5_6
• /
• pp.657-670
• /
• 2015

In this paper, we provide numerical analysis of so-called Legendre Gauss-Radau and Legendre-Gauss collocation methods for ordinary differential equations. After recasting these collocation methods as Runge-Kutta methods, we prove that the Legendre-Gauss collocation method is equivalent to the well-known Gauss method, while the Legendre-Gauss-Radau collocation method does not belong to the classes of Radau IA or Radau IIA methods in the Runge-Kutta literature. Making use of the well-established theory of Runge-Kutta methods, we study stability and accuracy of the Legendre-Gauss-Radau collocation method. Numerical experiments are conducted to confirm our theoretical results on the accuracy and numerical stability of the Legendre-Gauss-Radau collocation method, and compare Legendre-Gauss collocation method with the Gauss method.

• Structural Engineering and Mechanics
• /
• v.62 no.2
• /
• pp.151-162
• /
• 2017

Two $Chang-{\alpha}$ dissipative family methods and two $KR-{\alpha}$ family methods were developed for time integration recently. Although the four family methods are in the category of the dissipative structure-dependent integration methods, their performances may be drastically different due to the detrimental property of weak instability or overshoot for the two $KR-{\alpha}$ family methods. This weak instability or overshoot will result in an adverse overshooting behavior or even numerical instability. In general, the four family methods can possess very similar numerical properties, such as unconditional stability, second-order accuracy, explicit formulation and controllable numerical damping. However, the two $KR-{\alpha}$ family methods are found to possess a weak instability property or overshoot in the high frequency responses to any nonzero initial conditions and thus this property will hinder them from practical applications. Whereas, the two $Chang-{\alpha}$ dissipative family methods have no such an adverse property. As a result, the performances of the two $Chang-{\alpha}$ dissipative family methods are much better than for the two $KR-{\alpha}$ family methods. Analytical assessments of all the four family methods are conducted in this work and numerical examples are used to confirm the analytical predictions.