Refine
Year of publication
Language
- English (182)
Keywords
Project
Application Area
- D (182) (remove)
An optimal control problem is studied for a quasilinear Maxwell equation of nondegenerate parabolic
type. Well-posedness of the quasilinear state equation, existence of an optimal control, and weak G\^ateaux-differentiability of the control-to-state mapping are proved. Based on these results, first-order necessary optimality conditions and
an associated adjoint calculus are derived.
A mathematical model is set up that can be useful for controlled voltage excitation in time-dependent electromagnetism.
The well-posedness of the model is proved and an associated optimal control problem is investigated. Here, the control
function is a transient voltage and the aim of the control is the best approximation of desired electric and magnetic fields in
suitable $L^2$-norms.
Special emphasis is laid on an adjoint calculus for first-order necessary optimality conditions.
Moreover, a {peculiar attention is devoted to propose a formulation for which the computational complexity of the finite element solution method is substantially reduced}.
We report the cancellation of the soliton self-frequency shift in nonlinear optical fibers. A soliton which interacts with a group velocity matched low intensity dispersive pump pulse, experiences a continuous blue-shift in frequency, which counteracts the soliton self- frequency shift due to Raman scattering. The soliton self-frequency shift can be fully compensated by a suitably prepared dispersive wave. We quantify this kind of soliton-dispersive wave interaction by an adiabatic approach and demonstrate that the compensation is stable in agreement with numerical simulations.
We consider the phenomenon of an optical soliton controlled (e.g. amplified) by a much weaker second pulse which is efficiently scattered at the soliton. An important problem in this context is to quantify the small range of parameters at which the interaction takes place. This has been achieved by using adiabatic ODEs for the soliton characteristics, which is much faster than an empirical scan of the full propagation equations for all parameters in question.
We consider scattering of small-amplitude dispersive waves at an intense optical soliton which constitutes a nonlinear perturbation of the refractive index. Specifically, we consider a single-mode optical fiber and a group velocity matched pair: an optical soliton and a nearly perfectly reflected dispersive wave, a fiber-optical analogue of the event horizon. By combining (i) an adiabatic approach that is used in soliton perturbation theory and (ii) scattering theory from Quantum Mechanics, we give a quantitative account for the evolution of all soliton parameters. In particular, we quantify the increase in the soliton peak power that may result in spontaneous appearance of an extremely large, so-called champion soliton. The presented adiabatic theory agrees well with the numerical solutions of the pulse propagation equation. Moreover, for the first time we predict the full frequency band of the scattered dispersive waves and explain an emerging caustic structure in the space-time domain.
Physics of short optical pulses is an important and active research area in nonlinear optics. In what follows we theoretically consider the most extreme representatives of short pulses that contain only several oscillations of electromagnetic field. Description of such pulses is traditionally based on envelope equations and slowly varying envelope approximation, despite the fact that the envelope is not “slow” and, moreover, there is no clear definition of such a “fast” envelope. This happens due to another paradoxical feature: the standard (envelope) generalized nonlinear Schrödinger equation yields very good correspondence to numerical solutions of full Maxwell equations even for few-cycle pulses, a thing that should not be. In what follows we address ultrashort optical pulses using Hamiltonian framework for nonlinear waves. As it appears, the standard optical envelope equation is just a reformulation of general Hamiltonian equations. In a sense, no approximations are required, this is why the generalized nonlinear Schrödinger equation is so effective. Moreover, the Hamiltonian framework greatly contributes to our understanding of “fast” envelope, ultrashort solitons, stability and radiation of optical pulses. Even the inclusion of dissipative terms is possible making the Hamiltonian approach an universal theoretical tool also in extreme nonlinear optics.
We derive an optoelectronic model based on a gradient formulation for the relaxation of electron-, hole- and photon-densities to their equilibrium state. This leads to a coupled system of partial and ordinary differential equations, for which we discuss the isothermal and the non-isothermal scenario separately.
We show the existence of solutions to a system of elliptic PDEs, that was recently introduced to describe the electrothermal behavior of organic semiconductor devices. Here, two diffculties appear: (i) the elliptic term in the current-flow equation is of p(x)-Laplacian type with discontinuous exponent p, which limits the use of standard methods, and (ii) in the heat equation, we have to deal with an a priori L1 term on the right hand side describing the Joule heating in the device. We prove the existence of a weak solution under very weak assumptions on the data. Our existence proof is based on Schauder’s fixed point theorem and the concept of entropy solutions for the heat equation. Here, the crucial point is the continuous dependence of the entropy solutions on the data of the problem.
Classical gradient systems have a linear relation between rates and driving forces. In generalized gradient systems we allow for arbitrary relations derived from general non-quadratic dissipation potentials. This paper describes two natural origins for these structures.
A first microscopic origin of generalized gradient structures is given by the theory of large-deviation principles. While Markovian diffusion processes lead to classical gradient structures, Poissonian jump processes give rise to cosh-type dissipation potentials.
A second origin arises via a new form of convergence, that we call EDP-convergence. Even when starting with classical gradient systems, where the dissipation potential is a quadratic functional of the rate, we may obtain a generalized gradient system in the evolutionary -limit. As examples we treat (i) the limit of a diffusion equation having a thin layer of low diffusivity, which leads to a membrane model, and (ii) the limit of difusion over a high barrier, which gives a reaction-diffusion system.
Large-area organic light-emitting diodes are thin-film multilayer devices that show pronounced self-heating and brightness inhomogeneities at high currents. As these high currents are typical for lighting applications, a deeper understanding of the mechanisms causing these inhomogeneities is necessary. We discuss the modeling of the interplay between current flow, self-heating, and heat transfer in such devices using a system of partial differential equations of thermistor type, that is capable of explaining the development of luminance inhomogeneities. The system is based on the heat equation for the temperature coupled to a p(x)-Laplace-type equation for the electrostatic potential with mixed boundary conditions. The p(x)-Laplacian allows to take into account non-Ohmic electrical behavior of the different organic layers. Moreover, we present analytical results on the existence, boundedness, and regularity of solutions to the system. A numerical scheme based on the finite-volume method allows for efficient simulations of device structures.
We show the existence of solutions to a system of elliptic PDEs, that was recently introduced to describe the electrothermal behavior of organic semiconductor devices. Here, two difficulties appear: (i) the elliptic term in the current-flow equation is of $p(x)$-Laplacian-type with discontinuous exponent $p$, which limits the use of standard methods, and (ii) in the heat equation, we have to deal with an a priori $L^1$ term on the right hand side describing the Joule heating in the device. We prove the existence of a weak solution under very weak assumptions on the data. Our existence proof is based on Schauder’s fixed point theorem and the concept of entropy solutions for the heat equation. Here, the crucial point is the continuous dependence of the entropy solutions on the data of the problem.
Group-velocity matched cross-phase modulation between a fundamental soliton and a dispersive wave-packet has been previously suggested for optical switching applications similar to an optical transistor. Moreover, the nonlinear interaction in the resulting group-velocity horizon can be exploited for adiabatic compression of the soliton down into the few-cycle regime. Here we show that both mechanisms can be combined. In such a transient compressor, parameters of the dispersive wave may then serve to actively control the soliton compression and adjust the pulse duration in the presence of disturbances. While a certain amount of control is already enabled by the delay between soliton and dispersive wave, the means of controlling the compression process are substantially enhanced by additionally manipulating the chirp of the dispersive wave. Moreover, controlling the chirp of the dispersive wave also enables correction for limitations of the compression scheme due to a self-frequency shift of the soliton or for uncompensated dispersion in the scheme. This substantially widens the practicality of the compression scheme and other applications of the highly efficient nonlinear interaction at the group-velocity horizon.
We discuss analytical and numerical methods for the optimization of optoelectronic devices by performing optimal control of the PDE governing the carrier transport with respect to the doping profile. First, we provide a cost functional that is a sum of a regularization and a contribution, which is motivated by the modal net gain that appears in optoelectronic models of bulk or quantum- well lasers. Then, we state a numerical discretization, for which we study optimized solutions for different regularizations and for vanishing weights.
We derive thermodynamically consistent models of reaction-diffusion equations
coupled to a heat equation. While the total energy is conserved, the total
entropy serves as a driving functional such that the full coupled system is a
gradient flow. The novelty of the approach is the Onsager structure, which is
the dual form of a gradient system, and the formulation in terms of the
densities and the internal energy. In these variables it is possible to assume
that the entropy density is strictly concave such that there is a unique
maximizer (thermodynamical equilibrium) given linear constraints on the total
energy and suitable density constraints.
We consider two particular systems of this type, namely, a
diffusion-reaction bipolar energy transport system, and a
drift-diffusion-reaction energy
transport system with confining potential. We prove corresponding
entropy-entropy production inequalities with explicitly calculable constants
and establish the convergence to thermodynamical equilibrium, at first in
entropy and further in $L^1$ using Cziszar-Kullback-Pinsker type inequalities.
In these notes we discuss two approaches to evolutionary Γ- convergence of gradient systems in Hilbert spaces. The formulation of the gradient system is based on two functionals, namely the energy functional and the dissipation potential, which allows us to employ Γ- convergence methods. In the first approach we consider families of uni- formly convex energy functionals such that the limit passage of the time-dependent problems can be based on the theory of evolutionary variational inequalities as developed by Daneri and Savar ́e 2010. The second approach uses the equivalent formulation of the gradient system via the energy-dissipation principle and follows the ideas of Sandier and Serfaty 2004.
We apply both approaches to rigorously derive homogenization limits for Cahn–Hilliard-type equations. Using the method of weak and strong two-scale convergence via periodic unfolding, we show that the energy and dissipation functionals Γ-converge. In conclusion, we will give specific examples for the applicability of each of the two approaches.
We study a stationary thermistor model describing the electrothermal behavior of organic semiconductor devices featuring non-Ohmic current-voltage laws and self-heating effects. The coupled system consists of the current-flow equation for the electrostatic potential and the heat equation with Joule heating term as source. The self-heating in the device is modeled by an Arrhenius-like temperature dependency of the electrical conductivity. Moreover, the non-Ohmic electrical behavior is modeled by a power law such that the electrical conductivity depends nonlinearly on the electric field. Notably, we allow for functional substructures with different power laws, which gives rise to a $p(x)$-Laplace-type problem with piecewise constant exponent.
We prove the existence and boundedness of solutions in the two-dimensional case. The crucial point is to establish the higher integrability of the gradient of the electrostatic potential to tackle the Joule heating term. The proof of the improved regularity is based on Caccioppoli-type estimates, Poincar\'e inequalities, and a Gehring-type Lemma for the $p(x)$-Laplacian. Finally, Schauder's fixed-point theorem is used to show the existence of solutions.
We consider a coupled system of two elliptic PDEs, where the elliptic term in the first equation shares the properties of the $p(x)$-Laplacian with discontinuous exponent, while in the second equation we have to deal with an a~priori $L^1$ term on the right hand side. Such a system of equations is suitable for the description of various electrothermal effects, in particular those, where the non-Ohmic behavior can change dramatically with respect to the spatial variable. We prove the existence of a weak solution under very weak assumptions on the data and also under general structural assumptions on the constitutive equations of the model. The main difficulty consists in the fact that we have to overcome simultaneously two obstacles - the discontinuous variable exponent (which limits the use of standard methods) and the $L^1$ right hand side of the heat equation. Our existence proof based on Galerkin approximation is highly constructive and therefore seems to be suitable also for numerical purposes.
In a strong constant electric field, a dielectric particle
immersed in a weakly conducting fluid exhibits spontaneous rotations. This phenomenon is known under the name of the Quincke effect. In the original setup the particle was suspended on a silk thread and performed torsional oscillations of remarkably high amplitude.
We derive the governing equations for this experiment, and
ascertain that onset of oscillations from the quiescent state corresponds to the supercritical Hopf bifurcation.
For the case of a soft thread, we characterize the regime of large-scale torsional relaxation oscillations: explicit estimates are derived for their period and amplitude,
effects of bifurcation delay are described. In a stronger electric field, these relaxation oscillations yield to small-scale erratic rotations of the pendulum.
Uncertainty is inevitable when solving science and engineering application problems. In the face of
uncertainty, it is essential to determine robust and risk-averse solutions. In this work,
we consider a class of PDE-constrained optimization problems in which the PDE coefficients
and inputs may be uncertain. We introduce two approximations for minimizing the
conditional value-at-risk for such PDE-constrained optimization problems. These approximations are based
on the primal and dual formulations of the conditional value-at-risk. For the primal problem,
we introduce a smooth approximation of the conditional value-at-risk in order to utilize
derivative-based optimization algorithms and to take advantage of the convergence properties
of quadrature-based discretizations. For this smoothed conditional value-at-risk, we prove
differentiability as well as consistency of our approximation. For the dual problem, we
regularize the inner maximization problem, rigorously derive optimality conditions, and demonstrate
the consistency of our approximation. Furthermore, we propose a fixed-point iteration that takes
advantage of the structure of the regularized optimality conditions and provides a means of calculating
worst-case probability distributions based on the given probability level. We conclude with numerical
results.
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.
In large-area Organic Light-Emitting Diodes (OLEDs) spatially inhomogeneous luminance at high power due to inhomogeneous current flow and electrothermal feedback can be observed. To describe these self-heating effects in organic semiconductors we present a stationary thermistor model based on the heat equation for the temperature coupled to a p-Laplace-type equation for the electrostatic potential with mixed boundary conditions. The p-Laplacian describes the non-Ohmic electrical behavior of the organic material. Moreover, an Arrhenius-like temperature dependency of the electrical conductivity is considered.
We introduce a finite-volume scheme for the system and discuss its relation to recent network models for OLEDs. In two spatial dimensions we derive a priori estimates for the temperature and the electrostatic potential and prove the existence of a weak solution by Schauder's fixed point theorem.
The numerical simulation of the band structure of three-dimensional dispersive metallic photonic crystals with face-centered cubic lattices leads to large-scale nonlinear eigenvalue problems, which are very challenging due to a high dimensional subspace associated with the eigenvalue zero and the fact that the desired eigenvalues (with smallest real part) cluster near the zero eigenvalues. For
the solution of the eigenvalue problem, a Newton-type iterative method is proposed and the nullspace-free method is applied to exclude the zero eigenvalues from the associated generalized eigenvalue problem. To find the successive eigenvalue/eigenvector pairs, we propose a new non-equivalence deflation method to transform converged eigenvalues to infinity, while all other eigenvalues remain unchanged. The deflated problem is then solved by the same Newton-type method, which uses a hybrid method that combines the Jacobi-Davidson, the shift-invert residual Arnoldi and nonlinear Arnoldi methods to compute the clustered eigenvalues. Numerical results illustrate that the method is robust even for the case of computing many eigenvalues in very large problems.
Several classes of optimal control of electromagnetic fields are considered. Special emphasis is
laid on a non-standard $H$-based formulation of the equations of electromagnetism in multiply connected conductors. By this technique, the Maxwell equations can be solved with reduced computational complexity. While the magnetic field $H$ in the conductor is obtained from an elliptic equation
with the $\curl \sigma^{-1} \curl$ operator, an elliptic equation with the $\div \mu \nabla$ operator is set up for a potential $\psi$ in the isolator.
Both equations are coupled by appropriate interface conditions. In all problems, the
electrical current is controlled in the conducting domain. Several types of control functions are discussed. In particular, the problem of sparse optimal control is investigated in a package of electrical wires. For all problems, the associated sensitivity
analysis is performed.
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.
Using a standard first-order optimality condition for nonsmooth optimization problems, a general framework for a descent method is developed. This setting is applied to a class of mathematical programs with equilibrium constraints in function space from which a new algorithm is derived. Global convergence of the algorithm is demonstrated in function space and the results are then illustrated by numerical experiments.
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.
The long standing problem is discussed of how to deflate the part associated with the eigenvalue infinity in a structured matrix pencil using structure preserving unitary transformations. We derive such a deflation procedure and apply this new technique to symmetric, Hermitian or alternating pencils and in a modified form to (anti)-palindromic pencils. We present a detailed error and perturbation analysis of this and other deflation procedures and demonstrate the properties of the new algorithm with several numerical examples.
We study the approximation of Wasserstein gradient structures
by their finite-dimensional analog. We show that simple
finite-volume discretizations of the linear Fokker-Planck
equation exhibit the recently established entropic gradient-flow
structure for reversible Markov chains. Then we reprove the
convergence of the discrete scheme in the limit of vanishing
mesh size using only the involved gradient-flow structures.
In particular, we make no use of the linearity of the equations
nor of the fact that the Fokker-Planck equation is of second order.
Resolving thin conducting sheets for shielding or even skin layers inside by the mesh of numerical methods like the finite element method (FEM) can be avoided by using impedance transmission conditions (ITCs). Those ITCs shall provide an accurate approximation for small sheet thicknesses $d$, where the accuracy is best possible independent of the conductivity or the frequency being small or large -- this we will call robustness. We investigate the accuracy and robustness of popular and recently developed ITCs, and propose robust ITCs which are accurate up to $O(d^2)$.
We investigate the convergence of an implicit Voronoi finite volume
method for reaction-diffusion problems including nonlinear diffusion
in two space dimensions. The model allows to handle heterogeneous
materials and uses the chemical potentials of the involved species as
primary variables. The numerical scheme uses boundary conforming Delaunay
meshes and preserves positivity and the dissipative property of the
continuous system. Starting from a result on the global stability of
the scheme (uniform, mesh-independent global upper and lower bounds),
we prove strong convergence of the chemical activities and their gradients
to a weak solution of the continuous problem. In order to illustrate the
preservation of qualitative properties by the numerical scheme, we present
a long-term simulation of the Michaelis-Menten-Henri system. Especially,
we investigate the decay properties of the relative free energy and the
evolution of the dissipation rate over several magnitudes of time, and
obtain experimental orders of convergence for these quantities.
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.
Shielding sheets are commonly used in the protection of electronic devices. With their large aspect ratios they become a serious issue for the direct application of boundary element methods due to the occurring almost singular integrals. Impedance transmission conditions allow for boundary element formulations where only sheet mid-line has to be triangulated. We propose and analyse boundary element methods of second kind for impedance transmission conditions of two different types for the time-harmonic eddy current problem in two dimensions. In numerical experiments we study the effect on the magnetic shielding efficiency of the discretisation error due the boundary element discretisation and the modelling error due to the
impedance transmission condition.
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.
Enforcing solvability of a nonlinear matrix equation and estimation of multivariate ARMA time series
(2013)
The matrix equation $X+AX^{-1}A^T=B$, arising in parameter estimation of certain time series models,
is solvable only for certain values of the matrices $A,B$.
We present a numerical method to modify $A,B$ in order to make the matrix equation solvable.
Since solvability depends on the location of the eigenvalues of the palindromic matrix polynomial $\lambda^2 A+\lambda B+A^T$,
our method works by moving those eigenvalues to specified locations using first order spectral perturbation theory.
The method is heuristic but works in practice, as is supported by several compelling numerical examples.
These examples arise from parameter estimation of a common time series model, the multivariate ARMA(1,1).
We introduce a new operator for stabilizing error that arises from the weak enforcement of mass conservation in finite element simulations of incompressible flow problems. We show this new operator has a similar positive effect on velocity error as the well-known and very successful grad-div stabilization operator, but the new operator is more attractive from an implementation standpoint because it yields a sparser block structure matrix. That is, while grad-div produces fully coupled block matrices (i.e. block-full), the matrices arising from the new operator are block-upper triangular in two dimensions, and in three dimensions the 2,1 and 3,1 blocks are empty. Moreover, the diagonal blocks of the new operator's matrices are identical to those of grad-div. We provide error estimates and numerical examples for finite element simulations with the new operator, which reveals the significant improvement in accuracy it can provide. Solutions found using the new operator are also compared to those using usual grad-div stabilization, and in all cases, solutions are found to be very similar.
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.
Periodic Solutions to Dissipative Hyperbolic Systems. II: Hopf Bifurcation for Semilinear Problems
(2013)
We consider boundary value problems for semilinear hyperbolic systems of the type
$$
\partial_tu_j + a_j(x,\la)\partial_xu_j + b_j(x,\la,u) = 0, \; x\in(0,1), \;j=1,\dots,n
$$
with smooth coefficient functions $a_j$
and $b_j$
such that
$b_j(x,\la,0) = 0$ for all $x \in [0,1]$, $\la \in \R$, and $j=1,\ldots,n$.
We state conditions for Hopf bifurcation, i.e.,
for existence, local uniqueness (up to phase shifts), smoothness and smooth dependence
on $\la$
of time-periodic solutions bifurcating from the zero stationary solution. Furthermore,
we derive a formula which determines the bifurcation direction.
The proof is done by means of a Liapunov-Schmidt reduction procedure.
For this purpose, Fredholm properties of the linearized
system and implicit function
theorem techniques are used.
There are at least two distinguishing features of Hopf bifurcation theorems for hyperbolic PDEs in comparison with those for parabolic PDEs or for ODEs:
First, the question if a non-degenerate time-periodic solution depends smoothly on the system parameters
is much more delicate. And second,
a sufficient amount of dissipativity is needed in the system, and a priori
it is not clear how to verify this in terms of the data of the PDEs and of the boundary conditions.
Periodic Solutions to Dissipative Hyperbolic Systems. I: Fredholm Solvability of Linear Problems
(2013)
This paper concerns linear first-order hyperbolic systems in one space dimension of the type
$$
\partial_tu_j + a_j(x,t)\partial_xu_j + \sum\limits_{k=1}^nb_{jk}(x,t)u_k = f_j(x,t),\; x \in (0,1),\; j=1,\ldots,n,
$$
with periodicity conditions in time and reflection boundary conditions in space. We state a kind of dissipativity condition (depending on the coefficients $a_j$ and $b_{jj}$ and the boundary reflection coefficients), which implies Fredholm solvability of the problem, i.e., either there is a nontrivial solution to the homogeneous problem (in this case the space of such solutions has finite dimension) or the nonhomogeneous problem is uniquely solvable for any right-hand side (in this case the solution depends continuously on the right-hand side). In particular, under those conditions no small denominator effects occur.
Our results work for many non-strictly hyperbolic systems, but they are new even in the case of strict hyperbolicity.
Finally, in the case that all coefficients $a_j$ are $t$-independent, we show that the solutions are $C^\infty$-smooth if the data are $C^\infty$-smooth.
We examine robustness of exponential dichotomies of boundary value problems for general linear first-order one-dimensional hyperbolic systems. The boundary conditions are supposed to be of types ensuring smoothing solutions in finite time, which includes reflection boundary conditions. We show that the dichotomy survives in the space of continuous functions under small perturbations of all coefficients in the differential equations.
We give an exposition of recent results on regularity and Fredholm properties for first-order one-dimensional hyperbolic PDEs. We show that large classes of boundary operators cause an effect that smoothness increases with time. This property is the key in finding regularizers
(parametrices) for hyperbolic problems. We construct regularizers for periodic problems for dissipative first-order linear hyperbolic PDEs and show that these problems are modeled by Fredholm operators of index zero.
We consider systems of reaction-diffusion equations as gradient systems with respect to an entropy functional and a dissipation metric given in terms of a so-called Onsager operator, which is a sum of a diffusion part of Wasserstein type and a reaction part. We provide methods for establishing geodesic $\lambda$-convexity of the entropy functional by purely differential methods, thus circumventing arguments from mass transportation. Finally, several examples, including a drift-diffusion system, provide a survey on the applicability of the theory.
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.
In incompressible flows with vanishing
normal velocities at the boundary, irrotational forces in the momentum
equations should be balanced
completely by the pressure gradient.
Unfortunately, nearly all available discretization methods for incompressible flows violate this property.
The origin of the problem is that discrete velocity approximations
of incompressible flows are usually not
divergence-free. Hence, the use of divergence-free velocity reconstructions is
proposed wherever an $L^2$ scalar product appears in the discrete
variational formulation.
The approach is illustrated and applied to a nonconforming MAC-like discretization for unstructured Delaunay grids.
It is numerically demonstrated that a divergence-free velocity reconstruction based on the lowest-order Raviart-Thomas element
increases the robustness and accuracy of an existing convergent discretization, when irrotational forces appear in the momentum equations.
Recent research has shown that
in some practically relevant situations like multi-physics flows[11]
divergence-free mixed finite elements may have a significantly
smaller discretization error than standard non-divergence-free
mixed finite elements. In order to judge the overall performance of
divergence-free mixed finite elements, we
investigate linear solvers for the saddle point linear systems arising in $((P_k)^d,P_{k-1}^{disc})$ Scott-Vogelius finite element implementations of the incompressible Navier-Stokes equations. We investigate both direct and iterative solver methods.
Due to discontinuous pressure elements in the case of Scott-Vogelius elements, considerably more solver strategies seem to deliver promising results than in the case of standard mixed finite elements like
Taylor-Hood elements. For direct methods, we extend recent preliminary work using sparse banded solvers on the penalty method formulation to finer meshes, and discuss extensions. For iterative methods, we test augmented Lagrangian and H-LU preconditioners with GMRES, on both full and statically condensed systems.
Several numerical experiments are provided that show these classes of solvers are well suited for use with Scott-Vogelius elements, and could deliver an interesting overall performance in several applications.
We consider discretizations for reaction-diffusion systems with nonlinear
diffusion in two space dimensions. The applied model allows to handle heterogeneous
materials and uses the chemical potentials of the involved species as primary variables.
We propose an implicit Voronoi finite volume discretization on regular Delaunay
meshes that allows to prove uniform, mesh-independent global upper and lower L1
bounds for the chemical potentials. These bounds provide the main step for a convergence
analysis for the full discretized nonlinear evolution problem. The fundamental
ideas are energy estimates, a discrete Moser iteration and the use of discrete
Gagliardo-Nirenberg inequalities. For the proof of the Gagliardo-Nirenberg inequalities
we exploit that the discrete Voronoi finite volume gradient norm in 2d coincides
with the gradient norm of continuous piecewise linear finite elements.
Persistence of rogue waves in extended nonlinear Schrödinger equations: Integrable Sasa-Satsuma case
(2012)
We present the lowest order rogue wave solution of the Sasa-Satsuma equation (SSE) which is one of the integrable extensions of the nonlinear Schrödinger equation (NLSE). In contrast to the Peregrine solution of the NLSE, it is significantly more involved and contains polynomials of fourth order rather than second order in the corresponding expressions. The correct limiting case of Peregrine solution appears when the extension parameter of the SSE is reduced to zero.
When simulating isolated resonators, the application of transparent boundary conditions causes the approximated spectrum to be polluted with spurious solutions. Distinguishing these artificial solutions from solutions with a physical meaning is often difficult and requires a priori knowledge of the spectrum or the expected field distribution of resonant states. We present an implementation of the pole condition that distinguishes between incoming and outgoing waves by the location of the poles of their Laplace transform as transparent boundary condition. This implementation depends on one tuning parameter. We will use the sensitivity of the computed solutions to perturbations of this parameter as a means to identify spurious solutions. To obtain global statements, we will combine this technique with a convergence monitor for the boundary condition.