65N30 Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods
Refine
Year of publication
Language
- English (45)
Keywords
- finite element method (6)
- eigenvalue (4)
- AFEM (3)
- finite elements (3)
- geodesic finite elements (3)
- mixed finite elements (3)
- optimal control (3)
- Asymptotic Expansions (2)
- Currents (2)
- Matching (2)
With symmetric local absorbing boundary conditions for the Helmholtz equation scattering problems can be solved on a truncated domain, where the outgoing radiation condition is approximated by a Dirichlet-to-Neumann map with higher tangential derivatives on its outer boundary. Feng's conditions are symmetric local absorbing boundary conditions, which are based on an asymptotic expansion of the coefficients of the exact Dirichlet-to-Neumann map for large radia of the circular outer boundary. In this article we analyse the well-posedness of variational formulations with symmetric local absorbing boundary conditions in general and show how the modelling error introduced by Feng's conditions depends on the radius of the truncated domain.
Adaptive Numerical Solution of Eigenvalue Problems arising from Finite Element Models. AMLS vs. AFEM
(2015)
We discuss adaptive numerical methods for the solution of eigenvalue problems arising either
from the finite element discretization of a partial differential equation (PDE) or from discrete finite element modeling.
When a model is described by a partial differential equation, the adaptive finite element method
starts from a coarse finite element mesh which, based on a posteriori error estimators, is adaptively refined
to obtain eigenvalue/eigenfunction approximations of prescribed accuracy. This method is well established for classes of elliptic PDEs,
but is still in its infancy for more complicated PDE models.
For complex technical systems, the typical approach is to directly derive finite element models
that are discrete in space and are combined with macroscopic models to describe certain phenomena like damping or friction.
In this case one typically starts with a fine uniform mesh and computes eigenvalues and eigenfunctions using
projection methods from numerical linear algebra that are often combined with the algebraic multilevel
substructuring to achieve an adequate performance.
These methods work well in practice but their convergence and error analysis is rather difficult.
We analyze the relationship between these two extreme approaches. Both approaches have their pros and cons which are discussed in detail.
Our observations are demonstrated with several numerical examples.
The efficient and reliable computation of guided modes in photonic crystal wave-guides is of great importance for designing optical devices. Transparent boundary conditions based on Dirichlet-to-Neumann operators allow for an exact computation of well-confined modes and modes close to the band edge in the sense that no modelling error is introduced. The well-known super-cell method, on the other hand, introduces a modelling error which may become prohibitively large for guided modes that are not well-confined. The Dirichlet-to-Neumann transparent boundary conditions are, however, not applicable for all frequencies as they are not uniquely defined and their computation is unstable for a countable set of frequencies that correspond to so called Dirichlet eigenvalues. In this work we describe how to overcome this theoretical difficulty introducing Robin-to-Robin transparent boundary conditions whose construction do not exhibit those forbidden frequencies. They seem, hence, well suited for an exact and reliable computation of guided modes in photonic crystal wave-guides.
In this paper we present an efficient algorithm for the calculation of photonic crystal band structures and band structures of photonic crystal waveguides. Our method relies on the fact that the dispersion curves of the band structure are smooth functions of the quasi-momentum in the one-dimensional Brillouin zone. We show the derivation and computation of the group velocity, the group velocity dispersion, and any higher derivative of the dispersion curves. These derivatives are then employed in a Taylor expansion of the dispersion curves. We control the error of the Taylor expansion with the help of a residual estimate and introduce an adaptive scheme for the selection of nodes in the one-dimensional Brillouin zone at which we solve the underlying eigenvalue problem and compute the derivatives of the dispersion curves. The proposed algorithm is not only advantageous as it decreases the computational effort to compute the band structure but also because it allows
for the identification of crossings and anti-crossings of dispersion curves, respectively. This identification is not possible with the standard approach of solving the underlying eigenvalue problem at a discrete set of values of the quasi-momentum without taking the mode parity into account.
In this work we present a complete algorithm for the exact computation of the guided mode band structure in photonic crystal (PhC) wave-guides. In contrast to the supercell method, the used approach does not introduce any modelling error and is hence independent of the confinement of the modes. The approach is based on Dirichlet-to-Neumann (DtN) transparent boundary conditions that yield a nonlinear eigenvalue problem. For the solution of this nonlinear eigenvalue problem we present a direct technique using Chebyshev interpolation that requires a band gap calculation of the PhC in advance. For this band gap calculation –- we introduce as a very efficient tool –- a Taylor expansion of the PhC band structure. We show that our algorithm –- like the supercell method –- converges exponentially, however, its computational costs –- in comparison to the supercell method –- only increase moderately since the size of the matrix to be inverted remains constant.
In this paper, we introduce a high order finite element (FEM) implementation using perfectly matched layer (PML) for the scattering by plasmonic structures inside layered media. The PML is proven to be very accurate and efficient by a comparative analysis with a commercial FEM software and the Multiple Multipole Program (MMP). A convergence analysis using hp-adaptive refinement inside the PML layer
shows that adaptive mesh refinement inside the PML layer is most efficient. Based on this convergence analysis an hp-strategy is proposed, which shows a remarkable error reduction for small additional computational costs.
According to the Helmholtz decomposition,
the irrotational parts of the momentum balance equations of
the incompressible Navier-Stokes equations are balanced by the
pressure gradient.
Unfortunately, nearly all mixed methods for incompressible flows
violate this fundamental property, resulting in the well-known
numerical instability of poor mass conservation.
The origin of this problem is the
lack of $L^2$-orthogonality between discretely divergence-free velocities
and irrotational vector fields.
In order to cure this,
a new variational crime using divergence-free
velocity reconstructions is proposed.
Applying lowest order Raviart-Thomas velocity reconstructions
to the nonconforming Crouzeix-Raviart element
allows to construct
a cheap flow discretization for general 2d and 3d
simplex meshes that possesses the same
advantageous robustness properties like divergence-free flow
solvers.
In the Stokes case,
optimal a-priori error estimates for the velocity gradients
and the pressure are derived. Moreover, the discrete velocity
is independent of the continuous pressure.
Several detailed linear and nonlinear
numerical examples illustrate the theoretical findings.
We present new residual estimates based on Kato's square root theorem for spectral approximations of diagonalizable non-self-adjoint differential operators of convection-diffusion-reaction type. These estimates are incorporated as part of an hp-adaptive finite element algorithm for practical spectral computations, where it is shown that the
resulting a posteriori error estimates are reliable. Provided experiments demonstrate the efficiency and reliability of our approach.
The computation of guided modes in photonic crystal wave-guides is a key issue in the process of designing devices in photonic communications. Existing methods, such as the super-cell method, provide an efficient computation of well-confined modes. However, if the modes are not well-confined, the modelling error of the super-cell method becomes prohibitive and advanced methods applying transparent boundary conditions for periodic media are needed. In this work we demonstrate the numerical realization of a recently proposed Dirichlet-to-Neumann approach and compare the results with those of the super-cell method. For the resulting non-linear eigenvalue problem we propose an iterative solution based on Newton's method and a direct solution using Chebyshev interpolation of the non-linear operator. Based on the Dirichlet-to-Neumann approach, we present a formula for the group velocity of guided modes that can serve as an objective function in the optimization of photonic crystal wave-guides.
Grad-div stabilization has been proved to be a very useful tool in discretizations
of incompressible flow problems. Standard error analysis for inf-sup stable conforming pairs of
finite element spaces predicts that the stabilization parameter should be optimally chosen
to be $\mathcal O(1)$. This paper revisits this choice for the Stokes equations on the basis
of minimizing the $H^1(\Omega)$ error of the velocity and the $L^2(\Omega)$ error of the pressure.
It turns out, by applying a refined error analysis, that the optimal parameter choice is more subtle
than known so far in the literature. It depends on the used norm,
the solution, the family of finite
element spaces, and the type of mesh. Depending on the situation, the
optimal
stabilization parameter might range from being very small to very large.
The analytic results
are supported by numerical examples.