A numerical method for a free boundary value problem arising from chemical kinetics

✍ Scribed by Nicholas Batens; Roger Van Keer

• Publisher: Elsevier Science
• Year: 1999
• Tongue: English
• Weight: 217 KB
• Volume: 111
• Category: Article
• ISSN: 0377-0427
• DOI: 10.1016/s0377-0427(99)00142-9

✦ Synopsis

In this paper we deal with a numerical method for a complex free boundary value problem (FBP) in 1D, arising from chemical kinetics and involving two moving internal boundary points and six unknown concentration proÿles. Basically, the method consists of three steps: (1) a suitable ÿxed domain transformation for each of the three time-varying intervals, which results in a strongly nonlinear, nonlocal boundary value problem (BVP); (2) a nonstandard central di erence method with respect to the space variable, that takes properly into account the various transition conditions, in particular by using quadratic Lagrange interpolation to approximate the involved concentration proÿles near the relevant grid points; (3) a time discretisation method for the resulting initial value problem (IVP) for the nonlinear system of 1st order ODEs, that takes fully proÿt of the special structure of the mass matrix, implying that only a full 8 × 8 submatrix must be inverted at each discrete time point. A numerical example is presented.