Approximation and computation of the solution to the magnetosphere-ionosphere coupling equation via a mixed formulation

dc.contributor.advisorKloucek, Petr
dc.contributor.advisorToffoletto, Frank R.
dc.creatorWightman, Jennifer Lee
dc.date.accessioned2009-06-04T08:38:00Z
dc.date.available2009-06-04T08:38:00Z
dc.date.issued2004
dc.description.abstractThis study develops a numerical technique for the approximation of the magnetosphere-ionosphere (MI) coupling equation, which is a crucial step in the Rice Convection Model (RCM), a physical model that treats plasma in Earth's inner and middle magnetosphere via a multi-fluid approximation. The MI coupling equation is a second-order elliptic boundary value problem that describes conservation of current between the magnetosphere and the ionosphere. The current RCM solver is based on a finite difference scheme and produces unphysical results when the ionospheric conductance has large spatial gradients. We develop an alternative finite element approximation of the MI coupling equation, applying the method of fictitious domains to treat the high-latitude boundary condition along the immersed boundary Gamma, a boundary that varies in time and does not align with the computational grid. The result of using fictitious domains is a domain decomposition problem that we solve via a mixed finite element formulation. We compare both a nonconforming and conforming finite element approach within the framework of the mixed formulation. We are able to demonstrate that both the conforming and nonconforming methods generate solutions that are compatible with the current RCM solver when actual RCM data is used. Furthermore, we demonstrate on several analytic test examples that the finite element approximation is more accurate than the finite difference approximation. Therefore, we conclude that the finite element solver is more robust than the finite difference solver. In addition, we provide convergence results for the nonconforming approximation when the conductance coefficients are bounded and measurable, and we use spectral theory from the harmonic Steklov eigenproblem to derive a precise definition of the trace space on the interface Gamma. Our overall approximation technique is generalizable to a class of elliptic boundary value problems in which the boundary varies in time or does not align with a fixed grid. Finally, our numerical solver can be modified for use in the RCM-Jupiter that is currently being developed.
dc.format.extent217 p.en_US
dc.format.mimetypeapplication/pdf
dc.identifier.callnoTHESIS MATH.SCI. 2005 WIGHTMAN
dc.identifier.citationWightman, Jennifer Lee. "Approximation and computation of the solution to the magnetosphere-ionosphere coupling equation via a mixed formulation." (2004) Diss., Rice University. <a href="https://hdl.handle.net/1911/18837">https://hdl.handle.net/1911/18837</a>.
dc.identifier.urihttps://hdl.handle.net/1911/18837
dc.language.isoeng
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.
dc.subjectMathematics
dc.titleApproximation and computation of the solution to the magnetosphere-ionosphere coupling equation via a mixed formulation
dc.typeThesis
dc.type.materialText
thesis.degree.departmentMathematical Sciences
thesis.degree.disciplineEngineering
thesis.degree.grantorRice University
thesis.degree.levelDoctoral
thesis.degree.nameDoctor of Philosophy
Files
Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
3168161.PDF
Size:
16.7 MB
Format:
Adobe Portable Document Format