• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.2 no.1
• /
• pp.51-62
• /
• 2008

This paper proposes a second-order implicit constraint enforcement method which yields enhanced controllability compared to a first-order implicit constraints enforcement method. Although the proposed method requires solving a linear system twice, it yields superior accuracy from the constraints error perspective and guarantees the precise and natural movement of objects, in contrast to the first-order method. Thus, the proposed method is the most suitable for exact prediction simulations. This paper describes the numerical formulation of second-order implicit constraints enforcement. To prove its superiority, the proposed method is compared with the firstorder method using a simple two-link simulation. In this paper, there is a reasonable discussion about the comparison of constraints error and the analysis of dynamic behavior using kinetic energy and potential energy.

• Journal of Institute of Control, Robotics and Systems
• /
• v.20 no.4
• /
• pp.389-394
• /
• 2014

This paper presents the effect of a reflective force computed from a first-order-hold method on the stability of a haptic system. A haptic system is composed of a haptic device with a mass and a damper, a virtual spring, a sampler and a sample-and-hold. The boundary condition of the maximum virtual stiffness is analytically derived by using the Routh-Hurwitz criterion and the condition shows that the maximum virtual stiffness is proportional to the square root of the mass and the damper of a haptic device and also is inversely proportional to the sampling time to the power of three over two. The effectiveness of the derived condition is evaluated by the simulation. When the reflective forces are computed by using the first-order-hold method, the maximum available stiffness to guarantee the stability is increased several hundred times as large as when the zero-order-hold method is applied.

• 한국기계연구소 소보
• /
• s.14
• /
• pp.135-144
• /
• 1985

A method is described to obtain the first order force and second order steady force on the fixed two dimensional submerged or semisubmerged cylinders at infinite depth of water due to regular waves. The first order diffraction wave velocity potential which describes the flow diffracted by a body is obtained numerically using source distribution method on the mean wetted surface. And a technique to remove the irregular frequency phenomena of the source distribution method is also applied. The second order steady force is calculates by means of direct integration of the pressures on the body as derived from the first order velocity potential and is also computed by means of reflection wave height derives from momentum conservation theory. The results are compared with those of published works, and show good agreement.

• Journal of the Korean Mathematical Society
• /
• v.55 no.2
• /
• pp.471-506
• /
• 2018

In this article, some projection methods (or fractional-step methods) are proposed and analyzed for the micropolar Navier-Stokes equations (MNSE). These methods allow us to decouple the MNSE system into two sub-problems at each timestep, one is the linear and angular velocities system, the other is the pressure system. Both first-order and second-order projection methods are considered. For the classical first-order projection scheme, the stability and error estimates for the linear and angular velocities and the pressure are established rigorously. In addition, a modified first-order projection scheme which leads to some improved error estimates is also proposed and analyzed. We also present the second-order projection method which is unconditionally stable. Ample numerical experiments are performed to confirm the theoretical predictions and demonstrate the efficiency of the methods.

• Journal of Institute of Convergence Technology
• /
• v.3 no.2
• /
• pp.23-29
• /
• 2013

When a human operator interacts with a virtual wall that is modeled as a virtual spring model, the lager the stiffness of the virtual spring is, the more realistic the operator feels that the virtual wall is. In the previous studies, it is shown that the maximum available stiffness of a virtual spring to guarantee the stability can be increased when the first-order-hold method is applied, however the effects of a human impedance on the stability are not considered. This paper presents the effects of a human impedance on stability of haptic system with a virtual spring and a first-order-hold (FOH) method. The human impedance model is modeled as a linear second-order system model. The relations between the maximum available stiffness of a virtual spring and the human impedance such as a mass, a damping and a stiffness are analyzed through the MATLAB simulation. It is shown that the maximum available stiffness is proportional to the square root of the human mass or damping respectively.

• Journal of the Korean Mathematical Society
• /
• v.59 no.6
• /
• pp.1171-1184
• /
• 2022

We study the real-analytic continuation of local real-analytic solutions to the Cauchy problems of quasi-linear partial differential equations of first order for a scalar function. By making use of the first integrals of the characteristic vector field and the implicit function theorem we determine the maximal domain of the analytic extension of a local solution as a single-valued function. We present some examples including the scalar conservation laws that admit global first integrals so that our method is applicable.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.391-399
• /
• 2008

In this paper, we further generalize the work of Lin and Abel to the case of the first and the second order derivatives of energy release rates for two-dimensional, multiply cracked systems. The direct integral expressions are presented for the energy release rates and their first and second order derivatives. The salient feature of this numerical method is that the energy release rates and their first and second order derivatives can be computed in a single analysis. It is demonstrated through a set of examples that the proposed method gives expectedly decreasing, but acceptably accurate results for the energy release rates and their first and second order derivatives. The computed errors were approximately 0.5% for the energy release rates, $3\sim5%$ for their first order derivatives and $10\sim20%$ for their second order derivatives for the mesh densities used in the examples. Potential applications of the present method include a universal size effect model and a probabilistic fracture analysis of cracked structures.

• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• Communications for Statistical Applications and Methods
• /
• v.19 no.3
• /
• pp.333-344
• /
• 2012

We consider some statistical properties of the first order difference-based error variance estimator in nonparametric regression models with a single outlier. So far under an outlier(s) such difference-based estimators has been rarely discussed. We propose the first order difference-based estimator using the leave-one-out method to detect a single outlier and simulate the outlier detection in a nonparametric regression model with the single outlier. Moreover, the outlier detection works well. The results are promising even in nonparametric regression models with many outliers using some difference based estimators.