NUMERICAL STABILITY OF THE SAULYEV FINITE-DIFFERENCE ALGORITHMS FOR ELECTROCHEMICAL KINETIC SIMULATIONS - MATRIX STABILITY ANALYSIS FOR AN EXAMPLE PROBLEM INVOLVING MIXED BOUNDARY-CONDITIONS

被引:11
作者
BIENIASZ, LK [1 ]
OSTERBY, O [1 ]
BRITZ, D [1 ]
机构
[1] AARHUS UNIV,DATALOG AFDELING,DK-8000 AARHUS C,DENMARK
来源
COMPUTERS & CHEMISTRY | 1995年 / 19卷 / 04期
关键词
D O I
10.1016/0097-8485(95)00047-V
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
The stepwise numerical stability of the Saul'yev finite difference discretization of an example diffusional initial boundary value problem from electrochemical kinetics has been investigated using the matrix method of stability analysis. Special attention has been paid to the effect of the discretization of the mixed, linear boundary condition on stability, assuming the two-point, forward-difference approximation for the gradient at the left boundary (electrode). Criteria regulating the error growth in time have been identified. In particular it has been shown that, in contrast to the claims of unconditional stability of the Saul'yev algorithms, reported in the literature, the left-right variant of the Saul'yev algorithm becomes unstable for large values of the dimensionless diffusion parameter lambda = delta t/h(2), under mixed boundary conditions. This limitation is not, however, severe for most practical applications.
引用
收藏
页码:357 / 370
页数:14
相关论文
共 32 条
[1]   ON SOLUTION OF DIFFUSION EQUATIONS BY NUMERICAL METHODS [J].
BARAKAT, HZ ;
CLARK, JA .
JOURNAL OF HEAT TRANSFER, 1966, 88 (04) :421-&
[2]  
Bard A. J., 2001, ELECTROCHEMICAL METH, V2nd, P50
[3]   ELSIM - A PC PROGRAM FOR ELECTROCHEMICAL KINETIC SIMULATIONS - VERSION-2.0 - SOLUTION OF THE SETS OF KINETIC PARTIAL-DIFFERENTIAL EQUATIONS IN ONE-DIMENSIONAL GEOMETRY, USING FINITE-DIFFERENCE AND ORTHOGONAL COLLOCATION METHODS [J].
BIENIASZ, LK .
COMPUTERS & CHEMISTRY, 1993, 17 (04) :355-368
[4]   NUMERICAL STABILITY OF FINITE-DIFFERENCE ALGORITHMS FOR ELECTROCHEMICAL KINETIC SIMULATIONS - MATRIX STABILITY ANALYSIS OF THE CLASSIC EXPLICIT, FULLY IMPLICIT AND CRANK-NICOLSON METHODS, EXTENDED TO THE 3-POINT AND 4-POINT GRADIENT APPROXIMATION AT THE ELECTRODES [J].
BIENIASZ, LK ;
OSTERBY, O ;
BRITZ, D .
COMPUTERS & CHEMISTRY, 1995, 19 (04) :351-355
[5]   THE VONNEUMANN STABILITY OF FINITE-DIFFERENCE ALGORITHMS FOR THE ELECTROCHEMICAL KINETIC SIMULATION OF DIFFUSION COUPLED WITH HOMOGENEOUS REACTIONS [J].
BIENIASZ, LK .
JOURNAL OF ELECTROANALYTICAL CHEMISTRY, 1993, 345 (1-2) :13-25
[6]   ELECTROCHEMICAL KINETIC SIMULATIONS OF MIXED DIFFUSION HOMOGENEOUS REACTION PROBLEMS BY THE SAULYEV FINITE-DIFFERENCE ALGORITHMS [J].
BIENIASZ, LK ;
BRITZ, D .
ANALYTICA CHIMICA ACTA, 1993, 278 (01) :59-70
[7]   A METHOD-ORIENTED APPROACH TO THE FORMULATION OF ALGORITHMS FOR ELECTROCHEMICAL KINETIC SIMULATIONS [J].
BIENIASZ, LK .
JOURNAL OF ELECTROANALYTICAL CHEMISTRY, 1992, 340 (1-2) :19-34
[8]   EXPLICIT SOLUTION OF EQUATION OF HEAT CONDUCTION [J].
BIRTWISTLE, GM .
COMPUTER JOURNAL, 1968, 11 (03) :317-+
[9]   THE SAULYEV METHOD OF DIGITAL-SIMULATION UNDER DERIVATIVE BOUNDARY-CONDITIONS [J].
BRITZ, D ;
DASILVA, BM ;
AVACA, LA ;
GONZALEZ, ER .
ANALYTICA CHIMICA ACTA, 1990, 239 (01) :87-93
[10]   SOME NUMERICAL INVESTIGATIONS OF THE STABILITY OF ELECTROCHEMICAL DIGITAL-SIMULATION, PARTICULARLY AS AFFECTED BY 1ST-ORDER HOMOGENEOUS REACTIONS [J].
BRITZ, D ;
OSTERBY, O .
JOURNAL OF ELECTROANALYTICAL CHEMISTRY, 1994, 368 (1-2) :143-147