Refine
Year of publication
Language
- English (182)
Keywords
Project
Application Area
- D (182) (remove)
We study the class of linear differential-algebraic m-input m-output systems
which have a transfer function with proper inverse.
A sufficient condition for the transfer function to have proper inverse
it that the system has 'strict and non-positive relative degree'.
We present two main results:
First, a so called 'zero dynamics form' is derived: this form is - within the class of system equivalence - a simple "almost normal" form of the DAE; it is a counterpart to the well-known Byrnes-Isidori form for
ODE systems with strictly proper transfer function.
The 'zero dynamics form' is exploited to characterize structural properties such as
asymptotically stable zero dynamics,
minimum phase, and high-gain stabilizability.
The zero dynamics are characterized by (A,E,B)-invariant subspaces.
Secondly, it is shown that the 'funnel controller' (that is a static nonlinear output error feedback) achieves, for all DAE systems with asymptotically stable zero dynamics and transfer function with proper inverse, tracking of a reference signal by the output signal within a pre-specified funnel. This funnel determines the transient behaviour.
We prove the existence, uniqueness, regularity and smooth dependence
of the weak solution on the initial data for a certain class of semilinear first order
dissipative hyperbolic systems with spacially discontinuous coefficients. Such
kind of hyperbolic problems have succesfully been used to describe the dynamics
of distributed feedback multisection semiconductor lasers in recent years. We
show that in a suitable function space of continuous functions the weak solutions
generate a smooth semiflow.
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.
For a system of globally pulse-coupled phase-oscillators,
we derive conditions for stability of the completely synchronous
state and all possible two-cluster states and explain
how the different states are naturally connected via
bifurcations. The coupling is modeled using the phaseresponse-
curve (PRC), which measures the sensitivity of
each oscillator’s phase to perturbations. For large systems
with a PRC, which turns to zero at the spiking threshold,
we are able to find the parameter regions where multiple
stable two-cluster states coexist and illustrate this by an
example. In addition, we explain how a locally unstable
one-cluster state may form an attractor together will its
homoclinic connections. This leads to the phenomenon
of intermittent, asymptotic synchronization with abating
beats away from the perfect synchrony.
We use the traveling wave model for simulating and analyzing
nonlinear dynamics of complex semiconductor ring laser devices.
This modeling allows to consider temporal-spatial distributions
of the coun\-ter-pro\-pa\-ga\-ting slowly varying optical fields
and the carriers, what can be important when studying
non-homogeneous ring cavities, propagation of short pulses or fast switching.
By performing numerical integration of the model
equations we observe several dynamic regimes as well as transitions
between them. The computation of ring cavity modes explains some
peculiarities of these regimes.
A new approach to derive transparent boundary conditions (TBCs) for wave, Schrödinger, heat and drift-diffusion equations is presented. It relies on the pole condition and distinguishes between physical reasonable and unreasonable solutions by the location of the singularities of the spatial Laplace transform of the exterior solution. To obtain a numerical algorithm, a Möbius transform is applied to map the Laplace transform onto the unit disc. In the transformed coordinate the solution is expanded into a power series. Finally, equations for the coefficients of the power series are derived. These are coupled to the equation in the interior, and yield transparent boundary conditions.
Numerical results are presented in the last section, showing that the error introduced by the new approximate TBCs decays exponentially in the number of coefficients.
The pole condition approach for deriving transparent boundary conditions is extended to the time-dependent, two-dimensional case. Non-physical modes of the solution are identified by the position of poles of the solution's spatial Laplace transform in the complex plane. By requiring the Laplace transform to be analytic on some
problem dependent complex half-plane, these modes can be
suppressed. The resulting algorithm computes a finite number of coefficients of a series expansion of the Laplace transform, thereby providing an approximation to the exact boundary condition. The resulting error decays super-algebraically with the number of coefficients, so relatively few additional degrees of freedom are
sufficient to reduce the error to the level of the discretization error in the interior of the computational domain. The approach shows good results for the Schroedinger and the drift-diffusion equation
but, in contrast to the one-dimensional case, exhibits instabilities for the wave and Klein-Gordon equation. Numerical examples are shown that demonstrate the good performance in the former and the instabilities in the latter case.
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 show that many couplings between parabolic systems for processes in solids can be formulated as a gradient system with respect to the total free energy or the total entropy. This includes Allen-Cahn, Cahn-Hilliard, and reaction-diffusion systems and the heat equation. For this, we write the coupled system as an Onsager system $(X,\Phi,K)$ defining the evolution $\dot{U} = -K(U)D\Phi(U)$. Here $\Phi$ is the driving functional, while the Onsager operator $K(U)$ is symmetric and positive semidefinite. If the inverse $G = K^{-1}$ exists, the triple $(X,\Phi,G)$ defines a gradient system. Onsager systems are well suited to model bulk-interface interactions by using the dual dissipation potential $\Psi^*(U,\Xi) = 1/2 <\Xi,K(U)\Xi>$. Then, the two functionals $\Phi$ and $\Psi^*$ can be written as a sum of a volume integral and a surface integral, respectively. The latter may contain interactions of the driving forces in the interface as well as the traces of the driving forces from the bulk. Thus, capture and escape mechanisms like thermionic emission appear naturally in Onsager systems, namely simply through integration by parts.
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.
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.
We show that the spectrum of linear delay differential equations with
large delay splits into two different parts. One part, called the
strong spectrum, converges to isolated points when the delay parameter
tends to infinity. The other part, called the pseudocontinuous spectrum,
accumulates near criticality and converges after rescaling to a set
of spectral curves, called the asymptotic continuous spectrum. We
show that the spectral curves and strong spectral points provide a
complete description of the spectrum for sufficiently large delay
and can be comparatively easily calculated by approximating expressions.
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.
The perturbation and ADAE index of a degenerated hyperbolic system modelling a heat exchanger
(2007)
The heat exchanger in a heat pump can be modelled by the zero Mach-number limit of the Euler equations of compressible fluid flow. This system turns out to be a coupled hyperbolic/parabolic equation with coupled, time-dependent boundary conditions. Using the theory of abstract differential-algebraic equations it is shown that the frozen coefficient system has ADAE index 1. Moreover, the much stronger result is proven that the system has time-perturbation index one and space-perturbation index two even in the case of time-dependent boundary conditions. The results are stated in terms of the original physical variables. The estimates agree well with numerical experiments.
The well-known Kalman-Yakubovich-Popov Lemma establishes an equivalence between dissipativity and the solvability of a linear matrix inequality. In this paper we strengthen this result by showing the equivalence of dissipativity to the solvability of a so-called Lur'e equation, which mainly is a linear matrix inequality with a rank minimizing condition. Finally, we apply the result to standard systems to obtain the well-known result about the solvability of the algebraic Riccati equation.
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.
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.
Classical stability properties of solutions
that are well-known for ordinary differential
equations (ODEs) are generalized to differential-algebraic equations (DAEs).
A new test equation is derived for the analysis of numerical methods applied
to DAEs with respect to the stability of the numerical approximations.
Morevover, a stabilization technique is developed to improve the stability of classical DAE integration methods. The stability regions for these stabilized discretization methods are determined and it is shown that they much better reproduce the stability properties known for the ODE case
than in the unstabilized form.
Movies that depict the stability regions for several methods are included for interactive use.