• Journal of Mechanical Science and Technology
• /
• v.17 no.4
• /
• pp.538-544
• /
• 2003

Although many methodologies exist for determining the constrained equations of motion, most of these methods depend on numerical approaches such as the Lagrange multiplier's method expressed in differential/algebraic systems. In 1992, Udwadia and Kalaba proposed explicit equations of motion for constrained systems based on Gauss's principle and elementary linear algebra without any multipliers or complicated intermediate processes. The generalized inverse method was the first work to present explicit equations of motion for constrained systems. However, numerical integration results of the equation of motion gradually veer away from the constraint equations with time. Thus, an objective of this study is to provide a numerical integration scheme, which modifies the generalized inverse method to reduce the errors. The modified equations of motion for constrained systems include the position constraints of index 3 systems and their first derivatives with respect to time in addition to their second derivatives with respect to time. The effectiveness of the proposed method is illustrated by numerical examples.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.16 no.2
• /
• pp.21-27
• /
• 2008

An inverse method is introduced to construct the problem for the error analysis of the numerical solution of initial value problem. These problems constructed through this method have a known exact solution, even though analytical solutions are generally not obtainable. The process leading to the exact solution makes use of an initially available approximate numerical solution. A smooth interpolation of the approximate solution is forced to exactly satisfy the differential equation by analytically deriving a small forcing function to absorb all of the errors in the interpolated approximate solution. Using this special case exact solution, it is possible to investigate the relationship between global errors of a candidate numerical solution process and the associated tuning parameters for a given problem. Under the assumption that the original differential equation is well-posed with respect to the small perturbations, we thereby obtain valuable information about the optimal choice of the tuning parameters and the achievable accuracy of the numerical solution.

• Journal of applied mathematics & informatics
• /
• v.39 no.5_6
• /
• pp.655-676
• /
• 2021

A parameter uniform numerical scheme is proposed for solving singularly perturbed parabolic partial differential-difference convection-diffusion equations with a small delay and advance parameters in reaction terms and spatial variable. Taylor's series expansion is applied to approximate problems with the delay and advance terms. The resulting singularly perturbed parabolic convection-diffusion equation is solved by utilizing the implicit Euler method for the temporal discretization and finite difference method for the spatial discretization on a uniform mesh. The proposed numerical scheme is shown to be an ε-uniformly convergent accurate of the first order in time and second-order in space directions. The efficiency of the scheme is proved by some numerical experiments and by comparing the results with other results. It has been found that the proposed numerical scheme gives a more accurate approximate solution than some available numerical methods in the literature.

• Journal of Korean Port Research
• /
• v.10 no.1
• /
• pp.53-61
• /
• 1996

This paper describes the effect of wave control using submerged flat plate by the numerical calculation and the hydraulic model test. The boundary element method is used to develop a numerical solution for the flow field caused by monochromatic oblique waves incident upon an infinitely long, sumerged flat plate situated in arbitrary water depth. The effect of wave blocking is examined according to the change of length, submerged depth of flat plate and incident angles. Numerical results show that longer length, shallower submergence of flat plate and larger incident angles enhance the effect of wave blocking. To validate numerical analysis method, hydraulic model test was conducted in 2-D wave flume with 60 cm metal sheet. Reflected waves are extracted from water surface elevation in front of the location of a submerged plate by least square method with 3 wave gages. From comparing experimental results with numerical results, efficiency of numerical analysis method by this study could be confirmed well within wide ranges of wave frequencies.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.15 no.1
• /
• pp.1-8
• /
• 1990

In this paper, the transient electromagneti wave scattering at dielectric cylinder is studied by new numerical analysis method. Basic formulation of boundary integral equation (BIE) for numerical method is started weighted residual technique. BIE is made to two simultaneous equation at surface inner and outside point of dielectric cylinder in extended boundary condition (EBC) and surface boundary condition (SBC). Numerical method is used Boundary element method (BEM) that is two form, one is direct method and the other is indirect method, so that this method that transformes operator inversion martics is used numerical analysis. A good agreement of this numerical solution and the other results is obtained.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.04a
• /
• pp.76-83
• /
• 2001

In this study, a new method coupling of Element-Free Galerkin(EFG) method and Infinite Elements(IE) method is presented for extending application of the EFG method to engineering problems in unbounded domain. EFG method and IE method are briefly reviewed, and then the coupling procedure of the two methods is proposed. Numerical Algorithm by way of EFG-lE coupling technique is also developed. Numerical results illustrate the performance of the proposed technique. The accuracy of numerical solutions by the developed algorithm is guaranteed in comparing those of the other methods.

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

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.04a
• /
• pp.3-10
• /
• 1997

Meshless particle method is a numerical technique which does not use the concept of element. This method can easily handle special engineering problems which cause difficulty in the use of finite element method, however it has a drawback that essential boundary condition is not satisfied. In this paper, several studies for satisfying essential boundary conditions and enhancing the accuracy of solutions are discussed. Particular emphasis is placed on a new numerical technique in which finite elements are used on the boundaries to satisfy the essential boundary conditions and meshless particle method is used in the interior domain. For coupling of the two methods interface elements are introduced into the zone between the subdomains using meshless particle method and finite element method. The shape functions and the approximated displacement functions of the interface element are derived with the ramp function based on the shape function of finite elements. The whole numerical procedures are formulated by Galerkin method. Several numerical examples for enhancing the accuracy of solution in the meshless particle method and a new coupling method are presented.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.215-224
• /
• 2001

We present the whole basis of numerical method and useful formulae for general relativistic magnetohydrodynamic simulations in Kerr space-time.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.504-512
• /
• 2003

Simulation programs have been developed and used as an attempt to improve the accuracy of Non-Destructive Evaluation(NDE) techniques. The applicability of these programs is very limited, however, because it is difficult to describe the delicacy of the propagation of stress waves. To investigate the applicability of the finite element analysis for stress wave-based NDE techniques numerical simulation for Impact-Echo method and SASW method is performed. The numerical studies are performed to determine the essential parameters such as contact time of impact load, mesh size and time step size. These studies show that the choice of parameter is very important for improving the accuracy and confidence of the numerical procedure and, thereby, the applicability of the numerical analysis for stress wave-based NDE techniques