• Title/Summary/Keyword: Backward Euler method

Search Result 34, Processing Time 0.018 seconds

ANALYSIS OF A MESHFREE METHOD FOR THE COMPRESSIBLE EULER EQUATIONS

  • Kim, Yong-Sik;Pahk, Dae-Hyeon
    • Journal of the Korean Mathematical Society
    • /
    • v.43 no.5
    • /
    • pp.1081-1098
    • /
    • 2006
  • Mathematical analysis is made on a mesh free method for the compressible Euler equations. In particular, the Moving Least Square Reproducing Kernel (MLSRK) method is employed for space approximation. With the backward-Euler method used for time discretization, existence of discrete solution and it's $L^2-error$ estimate are obtained under a regularity assumption of the continuous solution. The result of numerical experiment made on the biconvex airfoil is presented.

CUBIC B-SPLINE FINITE ELEMENT METHOD FOR THE ROSENAU-BURGERS EQUATION

  • Xu, Ge-Xing;Li, Chun-Hua;Piao, Guang-Ri
    • East Asian mathematical journal
    • /
    • v.33 no.1
    • /
    • pp.53-65
    • /
    • 2017
  • Numerical solutions of the Rosenau-Burgers equation based on the cubic B-spline finite element method are introduced. The backward Euler method is used for discretization in time, and the obtained nonlinear algebraic system is changed to a linear system by the Newton's method. We show that those methods are unconditionally stable. Two test problems are studied to demonstrate the accuracy of the proposed method. The computational results indicate that numerical solutions are in good agreement with exact solutions.

An Assembly Sequence Planning of a Chip Mounter Using Transportation Algorithm (수송알고리즘에 의한 칩마운터의 조립순서계획)

  • Park, Tae-Hyung;Kim, Cheol-Han
    • Journal of Institute of Control, Robotics and Systems
    • /
    • v.6 no.9
    • /
    • pp.836-843
    • /
    • 2000
  • A sequence planning method is proposed to reduce the assembly time of gantey-type chip mounters with single head. The overall path of the chip mounter is divided into forward and backward path, and formulate the optimization problem is formulated as an transpoetation problem and an Euler's tour problem. The transportation alforithm is applied to find optimal backward path, and Euler's tour algorithm used to generate an assembly sequence. Simulation results are presented to verify the usefulness of the proposed method.

  • PDF

Multiaxial ratcheting assessment of Z2CND18.12N steel using modified A-V hardening rule

  • Xiaohui Chen;Yang Zhou;Wenwu Liu;Xu Zhao
    • Steel and Composite Structures
    • /
    • v.49 no.1
    • /
    • pp.1-17
    • /
    • 2023
  • Based on Ahmadzadeh-Varvani hardening rule (A-V model), multiaxial ratcheting effect of Z2CND18.12N austenitic stainless steel is simulated by ABAQUS with user subroutine UMAT. The results show that the predicted results of the origin multiaxial A-V model are lower than the experimental data, and it is difficult to control ratcheting strain rate. In order to improve the predicted capability of A-V model, the A-V model is modified. In this study. Moreover, under the assumption of the von Mises yield criterion and normal plasticity flow rule, we develop a numerical algorithm of plastic strain with the improved model to implement the finite element calculation of the model. Internal iteration in the numerical algorithm was implemented with the Euler backward method, which calculated the trial strain for each equilibrium iteration using the consistent tangent matrix. With a user subroutine, the proposed model is programmed into ABAQUS for a user - executable version. By simulating the uniaxial ratcheting of a round bar made of Z2CND18.12N austenitic stainless steel, we observe that the predicted results simulated by ABAQUS with UMAT are compared with the experimental data. The predicted results of the improved multiaxial A-V model are consistent well with the experimental data.

CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS

  • Chung, S.K.;Pani, A.K.;Park, M.G.
    • Journal of the Korean Mathematical Society
    • /
    • v.34 no.3
    • /
    • pp.515-531
    • /
    • 1997
  • In this paper, finite difference method is applied to approximate the generalized solutions of Sobolev equations. Using the Steklov mollifier and Bramble-Hilbert Lemma, a priori error estimates in discrete $L^2$ as well as in discrete $H^1$ norms are derived frist for the semidiscrete methods. For the fully discrete schemes, both backward Euler and Crank-Nicolson methods are discussed and related error analyses are also presented.

  • PDF

A DISCONTINUOUS GALERKIN METHOD FOR A MODEL OF POPULATION DYNAMICS

  • Kim, Mi-Young;Yin, Y.X.
    • 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.

DEVELOPMENT OF A HIGH-ORDER IMPLICIT DISCONTINUOUS GALERKIN METHOD FOR SOLVING COMPRESSIBLE NAVIER-STOKES EQUATIONS (압축성 Navier-Stokes 방정식 해를 위한 고차 정확도 내재적 불연속 갤러킨 기법의 개발)

  • Choi, J.H.;Lee, H.D.;Kwon, O.J.
    • Journal of computational fluids engineering
    • /
    • v.16 no.4
    • /
    • pp.72-83
    • /
    • 2011
  • A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

AN OVERVIEW OF BDF2 GAUGE-UZAWA METHODS FOR INCOMPRESSIBLE FLOWS

  • Pyo, Jae-Hong
    • 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.

Piecewise Linear Diode Models by Region Division for Circuit Simulations (회로 시뮬레이션을 위한 영역 분할식 구분적 선형 다이오드 모델)

  • Park, In-Gyu
    • 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.

  • PDF

A Configuration Design Sensitivity Analysis for Kinematically driven Mechanical Systems

  • Kim, D.W.;Yang, S.M.;Kim, H.W.;Bae, D.S.
    • Journal of the Korean Society of Manufacturing Technology Engineers
    • /
    • v.7 no.3
    • /
    • pp.110-117
    • /
    • 1998
  • A continuum-based configuration design sensitivity analysis method is developed for kinematically driven mechanical systems. The configuration design variable for mechanical systems is defined. The 3-1-3 Euler angle is employed as the orientation design variable. Kinematic admissibility conditions of configuration design change. Direct differentiation method is used to derive the governing equations of the design sensitivity. Numerical examples are presented to demonstrate the validity and effectiveness of the proposed method.

  • PDF