Adaptive finite element methods for variational inequalities: Theory and applications in finance
by Zhang, Chen-Song, Ph.D., UNIVERSITY OF MARYLAND, COLLEGE PARK, 2007, 203 pages; 3283416
 Abstract: We consider variational inequalities (VIs) in a bounded open domain Ω [special characters omitted] [special characters omitted] with a piecewise smooth obstacle constraint. To solve VIs, we formulate a fully-discrete adaptive algorithm by using the backward Euler method for time discretization and the continuous piecewise linear finite element method for space discretization. The outline of this thesis is the following. Firstly, we introduce the elliptic and parabolic variational inequalities in Hilbert spaces and briefly review general existence and uniqueness results (Chapter 1). Then we focus on a simple but important example of VI, namely the obstacle problem (Chapter 2). One interesting application of the obstacle problem is the American-type option pricing problem in finance. We review the classical model as well as some recent advances in option pricing (Chapter 3). These models result in VIs with integro-differential operators. Secondly, we introduce two classical numerical methods in scientific computing: the finite element method for elliptic partial differential equations (PDEs) and the Euler method for ordinary different equations (ODEs). Then we combine these two methods to formulate a fully-discrete numerical scheme for VIs (Chapter 4). With mild regularity assumptions, we prove optimal a priori convergence rate with respect to regularity of the solution for the proposed numerical method (Chapter 5). Thirdly, we derive an a posteriori error estimator and show its reliability and efficiency. The error estimator is localized in the sense that the size of the elliptic residual is only relevant in the approximate noncontact region, and the approximability of the obstacle is only relevant in the approximate contact region (Chapter 6). Based on this new a posteriori error estimator, we design a time-space adaptive algorithm and multigrid solvers for the resulting discrete problems (Chapter 7). In the end, numerical results for d = 1, 2 show that the error estimator decays with the same rate as the actual error when the space meshsize and the time step tend to zero. Also, the error indicators capture the correct local behavior of the errors in both the contact and noncontact regions (Chapter 8).

 Adviser Richard Nochetto School UNIVERSITY OF MARYLAND, COLLEGE PARK Source DAI/B 68-11, p. , Feb 2008 Source Type Dissertation Subjects Mathematics; Finance Publication Number 3283416
Access the complete dissertation:

» Find an electronic copy at your library.
 Use the link below to access a full citation record of this graduate work: http://gateway.proquest.com/openurl%3furl_ver=Z39.88-2004%26res_dat=xri:pqdiss%26rft_val_fmt=info:ofi/fmt:kev:mtx:dissertation%26rft_dat=xri:pqdiss:3283416 If your library subscribes to the ProQuest Dissertations & Theses (PQDT) database, you may be entitled to a free electronic version of this graduate work. If not, you will have the option to purchase one, and access a 24 page preview for free (if available).