German Title: Numerische Methoden zur Simulation und Optimierung von reagierenden Strömungen in katalytischen Monolithen

Preview

PDF, English Download (11MB) | Terms of use

Abstract

The aim of this work is to develop numerical methods and software for simulation and optimization of complex processes in catalytic monoliths to achieve better understanding of the physic-chemical processes in catalytic reactors. The fluid dynamics are modelled by the boundary layer equations (BLEs), which are a large system of parabolic partial differential equations (PDEs) with highly nonlinear boundary conditions arising from the coupling of surface processes with the flow field inside the channel. The BLEs are obtained by simplifying the comprehensive model described by the Navier-Stokes equations and applying the boundary approximation theory. The surface and gas-phase chemical reactions are described by detailed models. The PDEs are semi-discretized using the method of lines leading to a structured system of differential-algebraic equations (DAEs). The DAEs are integrated by an implicit method, based on backward differentiation formulas (BDF). We develop a new BDF code with tailored efficient and robust numerical methods by exploiting the structure, and by an appropriate scaling for ill-conditioned iteration matrices, and by computing consistent initial values. Efficient methods for computation of partial derivatives in the framework of automatic differentiation and of finite differences are introduced and compared. Our newly developed simulation tool is more stable than the existing simulation tool, and faster than by a factor of ten to more than 60, depending on the applications. To improve the performance of catalytic reactors (e.g., maximizing gas conversion or selectivity) we can control certain process conditions, such as temperature at the catalyst wall or the ratio of catalytic active surface area to the geometric surface area or the gas composition, the temperature, or the velocity at the inlet of the catalyst. It is the first time that this problem is generally formulated as an optimal control problem constrained by a system of PDEs describing the chemical fluid dynamics process and additional constraints. The direct shooting approach in combination with sequential quadratic programming (SQP) method is used for solving the resulting optimal control problem. An efficient numerical method for computation of the derivatives required by the SQP method is introduced. In addition, error analysis for the numerical Newton method is investigated in detail. We introduce a new error model. Based on our error model and analysis, the limiting accuracy of the solution of nonlinear equations by the numerical Newton method can be obtained. Our newly developed software package for simulation and optimization can be applied to different reaction mechanisms and channel settings with different initial/boundary conditions. This software is applied to two practical applications: catalytic combustion of methane and conversion of ethane to ethylene. The numerical results are presented. The simulation software provides a useful tool for the validation of reactions mechanisms. The software package allows, e.g., for a better design and operation of the conversion of natural gas to higher hydrocarbons or the improvement of exhaust treatment in cars.

Translation of abstract (German)

Ziel dieser Arbeit ist die Entwicklung numerischer Methoden und Programme zur Simulation und Optimierung komplexer Prozesse in katalytischen Monolithen, um die physikalisch-chemischen Prozesse in katalytischen Reaktoren besser verstehen zu können. Die Strömungen werden mittels der Grenzschichtgleichungen modelliert. Sie bilden ein großes System von partiellen Differentialgleichungen (PDEs) mit hochgradig nichtlinearen Randbedingungen, die sich aus der Kopplung der Oberflächenprozesse mit dem Strömungsfeld innerhalb des Kanalsergeben. Die Grenzschichtgleichungen werden abgeleitet, indem das Navier-Stokes-Modell vereinfacht und die Grenzschichtnäherung angewendet wird. Die Beschreibung der Gasphasen- und Oberflächenreaktionen erfolgt durch detaillierte Modelle. Die PDEs werden mit der Hilfe der Linienmethode semi-diskretisiert. Daraus ergibt sich ein differential-algebraisches Gleichungssystem. Das differential-algebraische Gleichungssystem wird durch eine implizite Methode integriert, die auf den "Backward-Differentiation-Formulae" (BDF) beruht. Es wird ein neuer BDF-Code mit speziell zugeschnittenen, effizienten und robusten numerischen Methoden entwickelt, der insbesondere alle Strukturen ausnutzt, schlecht-konditionierte Iterationsmatrizen geeignet skaliert und konsistente Anfangswerte berechnet. Effiziente Methoden zur Berechnung der partiellen Ableitungen im Rahmen der automatischen Differenzierung und der finiten Differenzen werden eingeführt und miteinander gekoppelt. Das neu entwickelte Simulationswerkzeug ist stabiler als das existierende und in Abhängigkeit der Anwendung 10 bis 60-mal schneller. Durch Variation bestimmter Prozessparameter lässt sich das Verhalten katalytischer Reaktoren verbessern (z.B. durch Maximierung von Umsatz oder Selektivität). Dazu zählen die Temperatur an der Katalysatorwand, das Verhältnis von katalytisch aktiver und Gesamtoberfläche sowie die Gaszusammensetzung, die Temperatur und die Geschwindigkeit am Eingang des Katalysators. Dieses Problem wird zum ersten Mal als Optimierungsproblem allgemein formuliert, das durch ein System von PDEs dargestellt wird. Dabei beschreiben die PDEs die reaktive Strömung sowie zusätzliche Bedingungen. Zur Lösung des sich ergebenden Optimierungsproblems wird ein "direktes Schieß verfahren" in Verbindung mit der Methode der sequentiellen quadratischen Programmierung (SQP) benutzt. Eine effiziente Vorgehensweise zur Berechnung der für die SQP-Methode erforderlichen Ableitungen wird dargestellt. Ein neues Modell zur Fehleranalyse der Newton-Methode wird eingeführt. Dadurch lässt sich die maximal erzielbare Genauigkeit der Lösung von nicht-linearen Gleichungen besser abschätzen. Das neu entwickelte Softwarepaket zur Simulation und Optimierung eignet sich für verschiedene Reaktionsmechanismen und Kanäle mit verschiedenen Anfangs- und/oder Randbedingungen. Die Software wird exemplarisch für zwei Anwendungen eingesetzt: katalytische Verbrennung von Methan und Umsetzung von Ethan zu Ethylen. Die numerischen Ergebnisse werden dargestellt. Diese Simulationssoftware ist geeignet zur Validierung von Reaktionsmechanismen. Sie ermöglicht die Optimierung chemischer Prozesse, wie zum Beispiel die Umsetzung von Erdgas in wertvolle Kohlenwasserstoffe oder die Abgasnachbehandlung in Kraftfahrzeugen.