ZIB-Report
Refine
Year of publication
Document Type
- ZIB-Report (1754)
- Doctoral Thesis (10)
- Habilitation (10)
- Article (5)
- In Proceedings (1)
Keywords
- integer programming (33)
- mixed integer programming (32)
- optimal control (29)
- KOBV (28)
- Kooperativer Bibliotheksverbund Berlin-Brandenburg (24)
- Bibliotheksverbund (20)
- Mixed Integer Programming (20)
- Integer Programming (17)
- line planning (14)
- mixed-integer programming (14)
Institute
- ZIB Allgemein (909)
- Mathematical Optimization (429)
- Numerical Mathematics (196)
- Visual and Data-centric Computing (92)
- Visual Data Analysis (88)
- Mathematical Optimization Methods (71)
- Mathematics of Transportation and Logistics (66)
- Computational Medicine (62)
- Digital Data and Information for Society, Science, and Culture (55)
- Applied Algorithmic Intelligence Methods (47)
SC-90-15
The mathematical modeling of macromolecular reactions leads to countable (possibly infinite) systems of ordinary differential equations (CODE's). This paper reviews two recent developments of the so-called discrete Galerkin method, which has been developed for the numerical treatment of countable systems, which arise e.g. in polymer chemistry. The first approach can be considered as a method of lines with moving basis functions and has been implemented recently in the program package MACRON. The second type of the Galerkin method is characterized by a so-called outer time discretization of the complete problem and an appropriate and efficient solution of the arising subproblems. This method is realized in the research code CODEX.
SC-90-14
This paper presents the new program package MACRON for the simulation of macromolecular kinetics including standard chemical reactions. Such problems lead to countable (possibly) infinite systems of ordinary differential equations (CODE's), which are numerically treated by the so-called discrete Galerkin method here. By a chemical compiler the required analytical preprocessing is performed, such that the complete reaction system, standard kinetics as well as macromolecular reactions, can be entered in the chemical formalism. Typical macromolecular reaction steps are chain addition, termination, chain transfer and degradation (cracking). In order to ensure efficiency and reliability, high sophisticated numerical routines are built within the package. MACRON can be used without a detailed knowledge of the used numerical methods. As an illustration the application of MACRON to some realistic problems is presented.
SC-90-13
In continuation of part I this paper develops a variable-order time discretization in Hilbert space based on a multiplicative error correction. Matching of time and space errors as explained in part I allows to construct an adaptive multilevel discretization of the parabolic problem. In contrast to the extrapolation method in time, which has been used in part I, the new time discretization allows to separate space and time errors and further to solve fewer elliptic subproblems with less effort, which is essential in view of the application to space dimension greater than one. Numerical examples for space dimension one are included which clearly indicate the improvement.
SC-90-12
In the present paper, the improvement of an incomplete factorization of a non-symmetric matrix A is discussed. Starting from the ideas of sparsity preserving quasi-Newton methods, an algorithm is developed which improves the approximation of A by the incomplete factorization maintaining the sparsity structure of the matrices. No renumbering of the unknowns or the admittance of additional fill-in is necessary. The linear convergence of the algorithm is proved under the assumption, that $ L $ and $ U $* have the same sparsity structure and an incomplete factorization with some reasonable approximation property exits. In combination with this algorithm, the method of incomplete factorization and its several modifications are applicable to a wider class of problems with improved convergence qualities. This is shown by a numerical example. {\bf Key Words:} non-symmetric linear system, sparse secant method, incomplete factorization. AMS(MOS) {\bf Subject Classifications:} 65F10, 65N20, 65N30.
SC-90-11
Parameter-dependent systems of nonlinear equations with symmetry are treated by a combination of symbolic and numerical computations. In the symbolic part of the algorithm the complete analysis of the symmetry occurs, and it is here where symmetrical normal forms, symmetry reduced systems, and block diagonal Jacobians are computed. Given a particular problem, the symbolic algorithm can create and compute through the list of possible bifurcations thereby forming a so-called tree of decisions correlated to the different types of symmetry breaking bifurcation points. The remaining part of the algorithm deals with the numerical pathfollowing based on the implicit reparametrisation as suggested and worked out by Deuflhard/Fiedler/Kunkel. The symmetry preserving bifurcation points are computed using recently developed augmented systems incorporating the use of symmetry. {\bf Keywords:} pathfollowing, mixed symbolic-numeric algorithm, parameter-dependent, nonlinear systems, linear representations.
SC-90-10
The potential distribution of reverse biased pn-junctions can be described by a double obstacle problem for the Laplacian. This problem is solved by a self adaptive Finite Element Method involving automatic termination criteria for the iterative solver, local error estimation and local mesh refinement. Special attention is paid to the efficient resolution of the geometries typically arising in semiconductor device simulation. The algorithm is applied to a reverse biased pn- junction with multi-step field plate and stop- electrode to illustrate its efficiency and reliability.
SC-90-09
The paper presents a theoretical characterization of the often observed asymptotic mesh independence of Newton's method, which means that Newton's method applied to discretized operator equations behaves essentially the same for all sufficiently fine discretizations. The theory does not need any uniform Lipschitz assumptions that were necessary in comparable earlier treatments. The refined Newton-Mysovskii theorem, which will be of interest in a wider context, gives both existence and uniqueness of the solution and quadratic convergence for sufficiently good starting points. Attention is restricted to Galerkin approximations even though similar results should hold for finite difference methods - but corresponding proofs would certainly be more technical. As an illustrative example, adaptive 1-D collocation methods are discussed.
SC-90-08
Automatic Generation of Reaction Mechanisms for Description of Oxidation of Higher Hydrocarbons.
(1990)
Oxidation mechanisms even for rather simple hydrocarbons like heptane consist due to the occurrence of many isomeric structures of thousands of reactions of hundreds of species. The automatic generation of these reaction mechanisms using artificial intelligence means is described. Results are presented for n-heptane-air mixtures, where a hand-written reaction mechanism tested against experimental data is available.
SC-90-07
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.
SC-90-06
Symplectic difference schemes have been shown to be a right formalism for numerical computation of Hamiltonian systems. They are suitable to long time computation and of good qualitative properties. These properties are ensured by the fact that a symplectic difference scheme approximating to a time-independent Hamiltonian system can be regarded as a perturbed time-dependent Hamiltonian system of the original one. That is, a solution of a symplectic difference scheme is a solution of a certain perturbed time dependent Hamiltonian system evaluated at discrete (time) points. This is the main result of the paper. Moreover, linear symplectic difference schemes approximating to a linear time-independent Hamiltonian system can be regarded as a perturbed time-independent Hamiltonian system. So it has all properties that a linear Hamiltonian system has. Based on these results, stochastic webs and chaos in symplectic difference schemes are also discussed. They will appear in numerical simulation for Hamiltonian systems, even with one degree of freedom.