• Title/Summary/Keyword: Stiff ODEs

Search Result 2, Processing Time 0.018 seconds

HIGH ORDER IMPLICIT METHOD FOR ODES STIFF SYSTEMS

  • Vasilyeva, Tatiana;Vasilev, Eugeny
    • Journal of applied mathematics & informatics
    • /
    • v.8 no.1
    • /
    • pp.165-180
    • /
    • 2001
  • This paper presents a new difference scheme for numerical solution of stiff system of ODE’s. The present study is mainly motivated to develop an absolutely stable numerical method with a high order of approximation. In this work a double implicit A-stable difference scheme with the sixth order of approximation is suggested. Another purpose of this study is to introduce automatic choice of the integration step size of the difference scheme which is derived from the proposed scheme and the one step scheme of the fourth order of approximation. The algorithm was tested by means of solving the Kreiss problem and a chemical kinetics problem. The behavior of the gas explosive mixture (H₂+ O₂) in a closed space with a mobile piston is considered in test problem 2. It is our conclusion that a hydrogen-operated engine will permit to decrease the emitted levels of hazardous atmospheric pollutants.

A Preconditioned Time Method for Efficient Calculation of Reactive Flow (예조건화 시간차분을 통한 화학반응유동의 효율적 계산)

  • Kim, Seong-Lyong;Jeung, In-Seuck;Choi, Jeong-Yeol
    • 한국연소학회:학술대회논문집
    • /
    • 1999.10a
    • /
    • pp.219-230
    • /
    • 1999
  • The Equations of Chemical kinetics are very stiff, which forces the use of an implicit scheme. The problem of implicit scheme, however, is that the jacobian must be solved at each time step. In this paper, we examined the methodology that can be stable without full chemical jacobian, This method is derived by applying the different time steps to the chemical source term. And the lower triangular chemical jacobian is derived. This is called the preconditioned time differencing method and represents partial implicit method. We show that this method is more stable in chemical kinetics than the full implicit method and that this is more efficient in supersonic combustion problem than the full jacobian method with same accuracy.

  • PDF