Positivity preserving solutions of partial integro-differential equations
Author(s)Lewis, Alexander M. (Alexander McDowell)
Massachusetts Institute of Technology. Dept. of Chemical Engineering.
Gregory J. McRae.
MetadataShow full item record
Differential equations are one of the primary tools for modeling phenomena in chemical engineering. While solution methods for many of these types of problems are well-established, there is growing class of problems that lack standard solution methods: partial integro-differential equations. The primary challenges in solving these problems are due to several factors, such as large range of variables, non-local phenomena, multi-dimensionality, and physical constraints. All of these issues ultimately determine the accuracy and solution time for a given problem. Typical solution techniques are designed to handle every system using the same methods. And often the physical constraints of the problem are not addressed until after the solution is completed if at all. In the worst case this can lead to some problems being over-simplified and results that provide little physical insight. The general concept of exploiting solution domain knowledge can address these issues. Positivity and mass-conservation of certain quantities are two conditions that are difficult to achieve in standard numerical solution methods. However, careful design of the discretizations can achieve these properties with a negligible performance penalty. Another important consideration is the stability domain. The eigenvalues of the discretized problem put restrictions on the size of the time step. For "stiff' systems implicit methods are generally used but the necessary matrix inversions are costly, especially for equations with integral components. By better characterizing the system it is possible to use more efficient explicit methods.(cont.) This work improves upon and combines several methods to develop more efficient methods. There are a vast number of systems that be solved using the methods developed in this work. The examples considered include population balances, neural models, radiative heat transfer models, among others. For the capstone portion, financial option pricing models using "jump-diffusion" motion are considered. Overall, gains in accuracy and efficiency were demonstrated across many conditions.
Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, 2009."May 15th, 2009."Includes bibliographical references (leaves 246-249).
DepartmentMassachusetts Institute of Technology. Department of Chemical Engineering
Massachusetts Institute of Technology