Numerical solution of non-linear algebraic equations with discontinuities

The solution of systems of non-linear algebraic equations with discontinuities in the solution search domain is considered. It is demonstrated that such problems are often very difficult to solve, even with the state of the art numerical solvers, and even when initial guesses close to the solutions are used. The application of constrained solution methods that do not require evaluation of function and derivative values outside of a predefined (feasible) subspace of the variables for solving such systems is considered. An algorithm is presented for identifying and handling of sub-expressions that introduce discontinuities. These are either removed by algebraic manipulations, or defined as boundaries of a feasible subspace. Using the proposed approach, it is demonstrated that a feasible solution for originally unsolvable problems can be found.

[1]  Vijay M. Naik,et al.  A continuum model for polymer adsorption at the solid-liquid interface , 1999 .

[2]  William H. Press,et al.  Numerical recipes in C (2nd ed.): the art of scientific computing , 1992 .

[3]  William R. Paterson.  A new method for solving a class of nonlinear equations , 1986 .

[4]  T. Yonemoto,et al.  Investigating the differences in acid separation behaviour on weak base ion exchange resins , 2000 .

[5]  E. Rogel,et al.  A lattice fluid approach to complex mixtures: natural gas and crude oil , 1997 .

[6]  Jorge R. Paloschi,et al.  Bounded homotopies to solve systems of sparse algebraic nonlinear equations , 1997 .

[7]  Ross E. Swaney,et al.  Robust solution of algebraic process modelling equations , 1994 .

[8]  Modéalisation de la copolyméarisation éathylégne α‐oléafine par catalyse ziegler—natta , 1991 .

[9]  Yash P. Gupta.  Bracketing Method for Online Solution of Low-Dimensional Nonlinear Algebraic Equations , 1995 .

[10]  S. Gupta,et al.  Steady state simulation of continuous‐flow stirred‐tank slurry propylene polymerization reactors , 1992 .

[11]  Lorenz T. Biegler,et al.  Iterative linear programming strategies for constrained simulation , 1991 .

[12]  Orit Shacham,et al.  Finding boundaries of the domain of definition for functions along a one-dimensional ray , 1990, TOMS.

[13]  J. M. Smith.  Chemical Engineering Kinetics , 1980 .

[14]  Abul Hasnat† and,et al.  Microphase-enhanced Reactions: Simultaneous Effects of Ion Coupling and Counterion Binding , 1999 .

[15]  Neima Brauner,et al.  A web-based library for testing performance of numerical software for solving nonlinear algebraic equations , 2001 .

[16]  Mordechai Shacham,et al.  Numerical solution of constrained non‐linear algebraic equations , 1986 .

[17]  Alexander P. Morgan,et al.  Chemical equilibrium systems as numerical test problems , 1990, TOMS.

[18]  William H. Press,et al.  Numerical recipes , 1990 .

[19]  David J. Goodman,et al.  Personal Communications , 1994, Mobile Communications.