Discontinuous Galerkin Methods for Parabolic Partial Differential Equations with Random Input Data

dc.contributor.advisorRiviere, Beatrice M.en_US
dc.contributor.committeeMemberHeinkenschloss, Matthiasen_US
dc.contributor.committeeMemberSymes, William W.en_US
dc.contributor.committeeMemberVannucci, Marinaen_US
dc.creatorLiu, Kunen_US
dc.date.accessioned2013-09-16T15:51:28Zen_US
dc.date.accessioned2013-09-16T15:51:33Zen_US
dc.date.available2013-09-16T15:51:28Zen_US
dc.date.available2013-09-16T15:51:33Zen_US
dc.date.created2013-05en_US
dc.date.issued2013-09-16en_US
dc.date.submittedMay 2013en_US
dc.date.updated2013-09-16T15:51:34Zen_US
dc.description.abstractThis thesis discusses and develops one approach to solve parabolic partial differential equations with random input data. The stochastic problem is firstly transformed into a parametrized one by using finite dimensional noise assumption and the truncated Karhunen-Loeve expansion. The approach, Monte Carlo discontinuous Galerkin (MCDG) method, randomly generates $M$ realizations of uncertain coefficients and approximates the expected value of the solution by averaging M numerical solutions. This approach is applied to two numerical examples. The first example is a two-dimensional parabolic partial differential equation with random convection term and the second example is a benchmark problem coupling flow and transport equations. I first apply polynomial kernel principal component analysis of second order to generate M realizations of random permeability fields. They are used to obtain M realizations of random convection term computed from solving the flow equation. Using this approach, I solve the transport equation M times corresponding to M velocity realizations. The MCDG solution spreads toward the whole domain from the initial location and the contaminant does not leave the initial location completely as time elapses. The results show that MCDG solution is realistic, because it takes the uncertainty in velocity fields into consideration. Besides, in order to correct overshoot and undershoot solutions caused by the high level of oscillation in random velocity realizations, I solve the transport equation on meshes of finer resolution than of the permeability, and use a slope limiter as well as lower and upper bound constraints to address this difficulty. Finally, future work is proposed.en_US
dc.format.mimetypeapplication/pdfen_US
dc.identifier.citationLiu, Kun. "Discontinuous Galerkin Methods for Parabolic Partial Differential Equations with Random Input Data." (2013) Diss., Rice University. <a href="https://hdl.handle.net/1911/71989">https://hdl.handle.net/1911/71989</a>.en_US
dc.identifier.slug123456789/ETD-2013-05-431en_US
dc.identifier.urihttps://hdl.handle.net/1911/71989en_US
dc.language.isoengen_US
dc.rightsCopyright is held by the author, unless otherwise indicated. Permission to reuse, publish, or reproduce the work beyond the bounds of fair use or other exemptions to copyright law must be obtained from the copyright holder.en_US
dc.subjectParabolic PDEsen_US
dc.subjectMonte Carlo Discontinuous Galerkinen_US
dc.subjectLocally mass conservationen_US
dc.subjectRandom input dataen_US
dc.subjectKernel PCAen_US
dc.subjectRandom permeabilityen_US
dc.subjectDarcy's Lawen_US
dc.subjectCoupled flow and transporten_US
dc.titleDiscontinuous Galerkin Methods for Parabolic Partial Differential Equations with Random Input Dataen_US
dc.typeThesisen_US
dc.type.materialTexten_US
thesis.degree.departmentComputational and Applied Mathematicsen_US
thesis.degree.disciplineEngineeringen_US
thesis.degree.grantorRice Universityen_US
thesis.degree.levelDoctoralen_US
thesis.degree.nameDoctor of Philosophyen_US
Files
Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
LIU-THESIS.pdf
Size:
5.84 MB
Format:
Adobe Portable Document Format
License bundle
Now showing 1 - 1 of 1
No Thumbnail Available
Name:
license.txt
Size:
1.61 KB
Format:
Item-specific license agreed upon to submission
Description: