A new convergence proof of the Adomian decomposition method for a mixed hyperbolic elliptic system of conservation laws (Q618129)

The authors consider a first order initial value problem for a system of nonlinear partial differential equations of the form \[ \frac{\partial U}{\partial t} = N(U) + H(x,t), \qquad U(x,0) = U_0(x), \tag{1} \] where \(N(U)\) is a nonlinear operator, \(H(x,t) \in \mathbb{R}^2\) is bounded on \((\alpha,\beta) \times [0,T]\), and construct the solution \(U=(u,v)^T\) as an infinite series \(U(x,t) = \sum_{k \geq 0} U_k(x,t)\) whose terms are determined recursively by means of the so-called Adomian polynomials. For practical computations, the series is conveniently truncated. A new algorithm for computing Adomian polynomials is proposed leading to a series with a faster convergence rate to the exact solution than the standard procedure. An error estimate is also obtained in this case. Since equation (1) includes in particular \(2 \times 2\) systems of conservation laws, a particular example of this class is chosen to illustrate the main features of the new approach. It is shown, in particular, that the method is stable inside the elliptic region present in the system, in contrast with other finite difference schemes.