Schneebeli, Anna. Interior penalty discontinuous galerkin methods for electromagnetic and acoustic wave equations. 2006, Doctoral Thesis, University of Basel, Faculty of Science.

PDF
5Mb 
Official URL: http://edoc.unibas.ch/diss/DissB_7760
Downloads: Statistics Overview
Abstract
Introduction: In this thesis we present and analyze the numerical approximation of the second
order electromagnetic and acoustic wave equation by the interior penalty (IP)
discontinuous Galerkin (DG) finite element method (FEM). In Part I we focus
on timeharmonic Maxwell source problems in the highfrequency regime. Part
II is devoted to the study of the IP DG FEM for timedependent acoustic and
electromagnetic wave equations.
We begin by stating Maxwell's equations in time and frequency domain. We
proceed by a variational formulation of Maxwell's equations, and describe the
key challenges that are faced in the analysis of the Maxwell operator. Then,
we review conforming finite element methods to discretize the second order
Maxwell operator. We end this general introduction with some numerical results
to highlight the performance and feasibility of conforming FEM for Maxwell's
equations.
Chapter 2: In this chapter, we introduce and analyze the interior penalty discontinuous
Galerkin method for the numerical discretization of the indefinite timeharmonic
Maxwell equations in highfrequency regime. Based on suitable duality arguments,
we derive apriori error bounds in the energy norm and the L2norm. In
particular, the error in the energy norm is shown to converge with the optimal
order O(hminfs;`g) with respect to the mesh size h, the polynomial degree `, and
the regularity exponent s of the analytical solution. Under additional regularity
assumptions, the L2error is shown to converge with the optimal order O(h`+1).
The theoretical results are confirmed in a series of numerical experiments on
triangular meshes.
The thesis' author's principal contributions are the proof of the L2error
bound in Section 2.6, and the proof of Lemma 2.4.1.
Chapter 3: We present and analyze an interior penalty method for the numerical discretization
of the indefinite timeharmonic Maxwell equations in mixed form. The
method is based on the mixed discretization of the curlcurl operator developed
in [44] and can be understood as a nonstabilized variant of the approach
proposed in [63]. We show the wellposedness of this approach and derive optimal
apriori error estimates in the energynorm as well as the L2norm. The
theoretical results are confirmed in a series of numerical experiments.
The thesis' author's principal contribution is the proof of the L2error bound
in Section 3.6.
Chapter 4: The symmetric interior penalty discontinuous Galerkin finite element method
is presented for the numerical discretization of the secondorder scalar wave
equation. The resulting stiffness matrix is symmetric positive definite and the
mass matrix is essentially diagonal; hence, the method is inherently parallel
and, leads to fully explicit time integration when coupled with an explicit timestepping
scheme. Optimal a priori error bounds are derived in the energy norm
and the L2norm for the semidiscrete formulation. In particular, the error
in the energy norm is shown to converge with the optimal order O(hminfs;`g)
with respect to the mesh size h, the polynomial degree `, and the regularity
exponent s of the continuous solution. Under additional regularity assumptions,
the L2error is shown to converge with the optimal order O(h`+1). Numerical
results confirm the expected convergence rates and illustrate the versatility of
the method.
Chapter 5: We develop the symmetric interior penalty discontinuous Galerkin (DG) method
for the spatial discretization in the method of lines approach of the timedependent
Maxwell equations in secondorder form. We derive optimal apriori
estimates for the semidiscrete error in the energy norm. For smooth solutions,
these estimates hold for DG discretizations on general finite element meshes.
For lowregularity solutions that have singularities in space, the theoretical estimates
hold on conforming, affine meshes. Moreover, on conforming triangular
meshes, we derive optimal error estimates in the L2norm. Finally, we valuate
our theoretical results by a series of numerical experiments.
order electromagnetic and acoustic wave equation by the interior penalty (IP)
discontinuous Galerkin (DG) finite element method (FEM). In Part I we focus
on timeharmonic Maxwell source problems in the highfrequency regime. Part
II is devoted to the study of the IP DG FEM for timedependent acoustic and
electromagnetic wave equations.
We begin by stating Maxwell's equations in time and frequency domain. We
proceed by a variational formulation of Maxwell's equations, and describe the
key challenges that are faced in the analysis of the Maxwell operator. Then,
we review conforming finite element methods to discretize the second order
Maxwell operator. We end this general introduction with some numerical results
to highlight the performance and feasibility of conforming FEM for Maxwell's
equations.
Chapter 2: In this chapter, we introduce and analyze the interior penalty discontinuous
Galerkin method for the numerical discretization of the indefinite timeharmonic
Maxwell equations in highfrequency regime. Based on suitable duality arguments,
we derive apriori error bounds in the energy norm and the L2norm. In
particular, the error in the energy norm is shown to converge with the optimal
order O(hminfs;`g) with respect to the mesh size h, the polynomial degree `, and
the regularity exponent s of the analytical solution. Under additional regularity
assumptions, the L2error is shown to converge with the optimal order O(h`+1).
The theoretical results are confirmed in a series of numerical experiments on
triangular meshes.
The thesis' author's principal contributions are the proof of the L2error
bound in Section 2.6, and the proof of Lemma 2.4.1.
Chapter 3: We present and analyze an interior penalty method for the numerical discretization
of the indefinite timeharmonic Maxwell equations in mixed form. The
method is based on the mixed discretization of the curlcurl operator developed
in [44] and can be understood as a nonstabilized variant of the approach
proposed in [63]. We show the wellposedness of this approach and derive optimal
apriori error estimates in the energynorm as well as the L2norm. The
theoretical results are confirmed in a series of numerical experiments.
The thesis' author's principal contribution is the proof of the L2error bound
in Section 3.6.
Chapter 4: The symmetric interior penalty discontinuous Galerkin finite element method
is presented for the numerical discretization of the secondorder scalar wave
equation. The resulting stiffness matrix is symmetric positive definite and the
mass matrix is essentially diagonal; hence, the method is inherently parallel
and, leads to fully explicit time integration when coupled with an explicit timestepping
scheme. Optimal a priori error bounds are derived in the energy norm
and the L2norm for the semidiscrete formulation. In particular, the error
in the energy norm is shown to converge with the optimal order O(hminfs;`g)
with respect to the mesh size h, the polynomial degree `, and the regularity
exponent s of the continuous solution. Under additional regularity assumptions,
the L2error is shown to converge with the optimal order O(h`+1). Numerical
results confirm the expected convergence rates and illustrate the versatility of
the method.
Chapter 5: We develop the symmetric interior penalty discontinuous Galerkin (DG) method
for the spatial discretization in the method of lines approach of the timedependent
Maxwell equations in secondorder form. We derive optimal apriori
estimates for the semidiscrete error in the energy norm. For smooth solutions,
these estimates hold for DG discretizations on general finite element meshes.
For lowregularity solutions that have singularities in space, the theoretical estimates
hold on conforming, affine meshes. Moreover, on conforming triangular
meshes, we derive optimal error estimates in the L2norm. Finally, we valuate
our theoretical results by a series of numerical experiments.
Advisors:  Grote, Marcus J. 

Committee Members:  Perugia, Ilaria 
Faculties and Departments:  05 Faculty of Science > Departement Mathematik und Informatik > Mathematik > Numerik (Grote) 
UniBasel Contributors:  Grote, Marcus J. 
Item Type:  Thesis 
Thesis Subtype:  Doctoral Thesis 
Thesis no:  7760 
Thesis status:  Complete 
Number of Pages:  132 
Language:  English 
Identification Number: 

edoc DOI:  
Last Modified:  22 Jan 2018 15:50 
Deposited On:  13 Feb 2009 15:51 
Repository Staff Only: item control page