• Structural Engineering and Mechanics
• /
• v.3 no.4
• /
• pp.391-410
• /
• 1995

The application of adaptive finite element method to dynamic problems is investigated. Both the kinetic and strain energy errors induced by space and time discretization were estimated in a consistent manner and controlled by the simultaneous use of the adaptive mesh generation and the automatic time stepping. Also an optimal ratio of spatial discretization error to temporal discretization error was discussed. In this study it was found that the best performance can be obtained when the specified spatial and temporal discretization errors have the same value. Numerical examples are carried out to verify the performance of the procedure.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.22 no.1
• /
• pp.95-102
• /
• 2014

This paper presents zero crossing point(ZCP) estimation methods considering discretization error for a high speed brushless DC(BLDC) motor drive. The ZCP is estimated by detecting the change of back-EMF polarity for the BLDC sensorless drive, and the discretization error exist on the estimated ZCP. The discretization error of the ZCP is a cause of the delay of a commutation timing of current and increment of a current ripple factor. Besides a delay of a ZCP estimation brings on the limitation of a speed range for the BLDC sensorless drive. The compensation method based on the error analysis with probability theory for reducing the effects of the discretization error of the ZCP is proposed. Also a ZCP estimation method according to the Back-EMF patterns is proposed to widen the speed range for the BLDC sensorless drive. The proposed methods are verified by the experiment.

• Journal of applied mathematics & informatics
• /
• v.31 no.3_4
• /
• pp.479-490
• /
• 2013

In this paper, we investigate a priori error estimates and superconvergence of varitional discretization for nonlinear parabolic optimal control problems with control constraints. The time discretization is based on the backward Euler method. The state and the adjoint state are approximated by piecewise linear functions and the control is not directly discretized. We derive a priori error estimates for the control and superconvergence between the numerical solution and elliptic projection for the state and the adjoint state and present a numerical example for illustrating our theoretical results.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.4
• /
• pp.203-214
• /
• 2017

In this paper we propose and analyze two a posteriori error estimators for the stabilized $P_1/P_1$ finite element discretization of the Stokes equations. These error estimators are computed by solving local Poisson or Stokes problems on elements of the underlying triangulation. We establish their asymptotic exactness with respect to the velocity error under certain conditions on the triangulation and the regularity of the exact solution.

• Structural Engineering and Mechanics
• /
• v.6 no.3
• /
• pp.339-345
• /
• 1998

The computation size and accuracy in the boundary element method are mutually coupled and strongly influenced by the formulations in boundary discretization and integration. This aspect is studied numerically for two-dimensional elastodynamic problems in the frequency-domain. The localized nature of error is observed in the computed results. A boundary discretization criterion is examined. The number of integration points in the boundary integration is studied to find the optimum number for accuracy. Useful information is obtained concerning the optimization in boundary discretization and integration.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.34D no.8
• /
• pp.15-25
• /
• 1997

The integration method of carier concentrations to redcue the discretization error of th box integratio method used in the discretization of the three-dimensional poisson's equation is presented. The carrier concentration is approximated in the closed form as an exponential function of the linearly varying potential in the element. The presented method is implemented in the three-dimensional poisson's equation solver running under the windows 95. The accuracy and the convergence chaacteristics of the three-dimensional poisson's equation solver are compared with those of DAVINCI for the PN junction diode and the n-MOSFET under the thermal equilibrium and the DC reverse bias. The potential distributions of the simulatied devices from the three-dimensional poisson's equation solver, compared with those of DAVINCI, has a relative error within 2.8%. The average number of iterations needed to obtain the solution of the PN junction diode and the n-MOSFET using the presented method are 11.47 and 11.16 while the those of DAVINCI are 21.73 and 23.0 respectively.

• Journal of the Korean Mathematical Society
• /
• v.60 no.2
• /
• pp.273-302
• /
• 2023

In this paper, a fully discrete numerical scheme for the viscoelastic Oldroyd flow is considered with an introduced auxiliary variable. Our scheme is based on the finite element approximation for the spatial discretization and the backward Euler scheme for the time discretization. The integral term is discretized by the right trapezoidal rule. Firstly, we present the corresponding equivalent form of the considered model, and show the relationship between the origin problem and its equivalent system in finite element discretization. Secondly, unconditional stability and optimal error estimates of fully discrete numerical solutions in various norms are established. Finally, some numerical results are provided to confirm the established theoretical analysis and show the performances of the considered numerical scheme.

• Structural Engineering and Mechanics
• /
• v.75 no.6
• /
• pp.685-699
• /
• 2020

Polygonal finite element provides a great flexibility in mesh generation of crack propagation problems where the topology of the domain changes significantly. However, the control of the discretization error in such problems is a main concern. In this paper, a polygonal-FEM is presented in modeling of crack propagation problems via an automatic adaptive mesh refinement procedure. The adaptive mesh refinement is accomplished based on the Zienkiewicz-Zhu error estimator in conjunction with a weighted SPR technique. Adaptive mesh refinement is employed in some steps for reduction of the discretization error and not for tracking the crack. In the steps that no adaptive mesh refinement is required, local modifications are applied on the mesh to prevent poor polygonal element shapes. Finally, several numerical examples are analyzed to demonstrate the efficiency, accuracy and robustness of the proposed computational algorithm in crack propagation problems.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2067-2077
• /
• 1996

In adaptive finite element analysis, r- and h-methods are generally used on the basis of a discretization error estimator. In this paper, an rh-method is proposed as a new adaptive method which can improve the adaptivity performance by using both of them. This suggested rh-method moves nodal coordinates of initially given model to adjust element discretization errors and thereafter performes the h-method tdo obtain the specified accuracy of finite element solutions. Numerical experiments for various plane problems were performed using 4-noded isoparametric quadrilateral elements. As a result, the rh-method has been shown to be an accurate and efficient adaptive analysis method to obtain as improved solution.

• Journal of applied mathematics & informatics
• /
• v.32 no.5_6
• /
• pp.707-719
• /
• 2014

In this paper, we investigate a variational discretization approximation of elliptic optimal control problems with control constraints. The state and the co-state are approximated by piecewise linear functions, while the control is not directly discretized. By using some proper intermediate variables, we derive a second-order convergence in $L^2$-norm and superconvergence between the numerical solution and elliptic projection of the exact solution in $H^1$-norm or the gradient of the exact solution and recovery gradient in $L^2$-norm. Then we construct a posteriori error estimates by using the superconvergence results and do some numerical experiments to confirm our theoretical results.