Consider the potentially stiff system For large values, the system becomes stiff as it contains components

Question:

Consider the potentially stiff systemdyi dt dy2 = -yi(t), = -y(t), dt y (0) = 1, 1/2 (0) = 1, t = [0, 1]. The analytical solution to this problem

For large α values, the system becomes stiff as it contains components that vary with different speeds, i.e. y2(t) approaches zero much faster than y1(t) does. The solutions to this problem for α = 2 and α = 200 000 are shown in Figure 6.11. By analyzing the
number of steps required by the different methods, it becomes clear how the stiff solvers outperform the explicit methods. Calculations have been performed for different values of α with different error tolerances; the results are summarized in Table 6.4.0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.2 0.4 t (a) 0.6 0.8 - 1 0.96 0.8 0.7 0.6 0.59 0.4 0.3 0.2 0.1 0 0.2

Table 6.4= 2 Relative Tolerance 10-; 10-6 ode23 11 ode45 10 ode23s 11 ode 15s  = 200 105 73 53 14 75  = 200 000 79 620

Fantastic news! We've Found the answer you've been seeking!

Step by Step Answer:

Related Book For  book-img-for-question

Mathematical Modeling In Chemical Engineering

ISBN: 9781107049697

1st Edition

Authors: Anders Rasmuson, Bengt Andersson, Louise Olsson, Ronnie Andersson

Question Posted: