Refine
Year of publication
Document Type
- ZIB-Report (25)
- In Proceedings (8)
- Article (3)
- Book chapter (2)
- Doctoral Thesis (1)
Is part of the Bibliography
- no (39)
Keywords
Institute
- ZIB Allgemein (24)
- Numerical Mathematics (6)
- Computational Medicine (1)
New, highly efficient techniques for the numerical simulation and parameter identification of large chemical reaction systems are surveyed. The survey particularly addresses to chemists, which are understood to be potential users of the distributed software packages. These packages are written in the form of interactive dialogue systems - thus enabling chemists to concentrate on the chemistry of their problem. Large scale examples from chemical research environments are included.
This report presents the final realization and implementation of a global inexact Newton method proposed by Deuflhard. In order to create a complete piece of software, a recently developed iterative solver (program GBIT) due to Deuflhard, Freund, Walter is adapted and serves as the standard iterative linear solver. Alternative linear iterative solvers may be adapted as well, e.g. the widely distributed code GMRES. The new software package GIANT (Global Inexact Affine Invariant Newton Techniques) allows an efficient and robust numerical solution of very large scale highly nonlinear systems. Due to the user friendly interface and its modular design, the software package is open for an easy adaptation to specific problems. Numerical experiments for some selected problems illustrate performance and usage of the package.
The paper surveys three aspects of chemical computing, which seem to play a role in recent developments. First, extrapolation methods for the numerical treatment of differential- algebraic equations are introduced. The associated extrapolation code LIMEX has reached a certain level of sophistication, which makes it a real competitor to the elsewhere widely used multi-step code DASSL of Petzold. Second, adaptive methods of lines for partial differential equations such as those arising in combustion problems are treated. Both static and dynamic regridding techniques are discussed in some detail. Finally, some new ideas about the treatment of the kinetic equations arising from polymer reactions are presented. The new feature of the suggested approach is the application of a Galerkin procedure using sets of orthogonal polynomials over a discrete variable (which, of course, in the case of polymer reactions is the polymer degree). The new approach may open the door to a new reliable low dimensional treatment of complex polymer reactions.
The rapidly increasing number of AIDS cases requires a realistic estimation of the future development of the HIV/AIDS disease. For that purpose we develop a large system of coupled nonlinear differential equations describing simultaneously the dynamics of the development of the disease, the population size, the gender and age structure. A set of 1650 coupled equations are linked by balanced parameters. The balancing procedure is described by a set of (formally) 2,178,000 additional algebraic conditions. As the considered system is stiff, it requires new special extrapolation methods combined with techniquees of dynamical sparsing for the solution of sparsely filled systems. According to our simulations we expect 19,0,000 deaths caused by AIDS in the Federal Republic of Germany (former territories) in the year 2000. Such an epidemical spread would tie up about 4-7 percent of the actual health care workers.