• Proceedings of the KWS Conference
• /
• 2001.10a
• /
• pp.116-118
• /
• 2001

We propose an implicit numerical implementation for Leblond's transformation plasticity constitutive equations , which are widely used in welded steel structure. We apply Euler backward scheme rule to integrate the equations and determine the consistent tangent modulus. The implementation may be used with updated Lagrangian formulation. we test a simple butt-welding process to compare with SYSWELD and discuss the accuracy.

• Communications of the Korean Mathematical Society
• /
• v.18 no.4
• /
• pp.767-779
• /
• 2003

We consider a model of population dynamics whose mortality function is unbounded. We approximate the solution of the model using a discontinuous Galerkin finite element for the age variable and a backward Euler for the time variable. We present several numerical examples. It is experimentally shown that the scheme converges at the rate of $h^{3/2}$ in the case of piecewise linear polynomial space.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.1
• /
• pp.70-82
• /
• 1995

A multi-layer simulation program is developed to estimate the ocean current considering sea bottom geometry. The so-called $\sigma$ coordinate system is introduced in vertical direction to describe sea bottom topography more accurately and effectively. Leapfrog scheme combined with Euler backward scheme is used to reduce computation error which may be possibly accumulated in time evolution by Leapfrog scheme alone. In this paper, very simple examples of rectangular basins with various bottom geometries were taken and the effect of sea bottom geometry on vertical structure of the ocean tidal current and its direction were investigated. Through comparisons between the present three dimensional calculation in which bottom topography is directly taken into consideration and the two dimensional calculation in which depth average concept is employed, it was found that magnitude of surface current and its direction could be largely affected by the sea bottom topography, particularly in shallow region with complex bottom shape.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2007.05a
• /
• pp.272-275
• /
• 2007

This paper is concerned with modified integration algorithm on the strain-space for rate and temperature dependent elasto-plastic constitutive relations in order to obtain more accurate results in numerical implementation. The proposed algorithm is integrated analytically using integration by part and chain rule and then is applied to the 2-stage Lobatto IIIA with second-order accuracy. It has advantage that is able to consider the convective stress rates on the yield surface of the strain-space. Also this paper is carried out the iteration procedure using the Newton-Raphson method to enforce consistency at the end of the step. And the performance of the proposed algorithm for rate and temperature dependent constitutive relation is illustrated by means of analysis of adiabatic shear bands.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.32 no.4
• /
• pp.319-326
• /
• 2008

In this study, we investigate the deformation behavior of F.C.C. single crystals containing micro- or submicron-sized voids by using three dimensional finite element methods. The locally homogeneous constitutive model for the rate-dependent crystal plasticity is integrated based on the backward Euler method and implemented into a finite element program (ABAQUS) by means of user-defined subroutine (UMAT). The unit cell analysis has been investigated to study the effect of stress triaxiality and crystallographic orientations on the growth and coalescence of voids in F.C.C. single crystals.

• Proceedings of the KIEE Conference
• /
• 2002.07b
• /
• pp.882-884
• /
• 2002

The transient voltage characteristics of capacitor self-exited induction generator are analyzed by the state equation which is obtained from the d-q axis equivalent circuit of stationary reference frame and torque equation. The d-q equivalent circuit is composed using the condition of stationary reference frame. The mutual inductance is only considered as a function of magnetizing current in the equivalent circuit. The characteristics are analyzed and discussed by the backward Euler method for various load conditions under specified initial conditions and input.

• Journal of the Chungcheong Mathematical Society
• /
• v.20 no.4
• /
• pp.419-432
• /
• 2007

Double-diffusive convection inside a triangular porous cavity is studied numerically. Galerkin finite element method is adopted to derive the discrete form of the governing differential equations. The first-order backward Euler scheme is used for temporal discretization with the second-order Adams-Bashforth scheme for the convection terms in the energy and species conservation equations. The Boussinesq-Oberbeck approximation is used to calculate the density dependence on the temperature and concentration fields. A parametric study is performed with the Lewis number, the Rayleigh number, the buoyancy ratio, and the shape of the triangle. The effect of gravity orientation is considered also. Results obtained include the flow, temperature, and concentration fields. The differences induced by varying physical parameters are analyzed and discussed. It is found that the heat transfer rate is sensitive to the shape of the triangles. For the given geometries, buoyancy ratio and Rayleigh numbers are the dominating parameters controlling the heat transfer.

• Proceedings of the KIEE Conference
• /
• 2008.04c
• /
• pp.106-109
• /
• 2008

Piecewise linear diode models are widely used for large-signal circuit analyses, especially power electronic circuit simulations. When using a piecewise linear diode model for simulation, a switching method to select a proper one among linear models is needed. The conventional switching method keeps the previous ON, OFF state information, and applies different switching conditions according to the state. However, this method has difficulties especially in extending to multi-piecewise linear models. This paper presents a switching method which appropriately divides the v-i plane into regions and select a linear model according to the region where the operating point(the voltage and the current of the diode) belongs. This switching method is easily extended to multi-Piecewise linear models. An example using the tableau analysis and the backward Euler integration is presented for verification.

• 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.15 no.3
• /
• pp.233-251
• /
• 2011

The Gauge-Uzawa method [GUM] in [9] which is a projection type algorithm to solve evolution Navier-Stokes equations has many advantages and superior performance. But this method has been studied for backward Euler time discrete scheme which is the first order technique, because the classical second order GUM requests rather strong stability condition. Recently, the second order time discrete GUM was modified to be unconditionally stable and estimated errors in [12]. In this paper, we contemplate several GUMs which can be derived by the same manner within [12], and we dig out properties of them for both stability and accuracy. In addition, we evaluate an stability condition for the classical GUM to construct an adaptive GUM for time to make free from strong stability condition of the classical GUM.