Marcus Waurick

A Hilbert space approach to fractional differential equations

We study fractional differential equations of Riemann-Liouville and Caputo type in Hilbert spaces. Using exponentially weighted spaces of functions defined on $\mathbb{R}$ , we define fractional operators by means of a functional calculus using the Fourier transform. Main tools are extrapolation- and interpolation spaces. Main results are the existence and uniqueness of solutions and the causality of solution operators for non-linear fractional differential equations.

Evolutionary Equations

This is the final version of the lecture notes of the 23rd Internet Seminar on Evolutionary Equations, see also https://www.mat.tuhh.de/isem23/.

Maximal Regularity for Non-Autonomous Evolutionary Equations

We discuss the issue of maximal regularity for evolutionary equations with non-autonomous coefficients. Here evolutionary equations are abstract partial-differential algebraic equations considered in Hilbert spaces. The catch is to consider time-dependent partial differential equations in an exponentially weighted Hilbert space. In passing, one establishes the time derivative as a continuously invertible, normal operator admitting a functional calculus with the Fourier--Laplace transformation providing the spectral representation. Here, the main result is then a regularity result for well-posed evolutionary equations solely based on an assumed parabolic-type structure of the equation and estimates of the commutator of the coefficients with the square root of the time derivative. We thus simultaneously generalise available results in the literature for non-smooth domains. Examples for equations in divergence form, integro-differential equations, perturbations with non-autonomous and rough coefficients as well as non-autonomous equations of eddy current type are considered.

Nonlocal homogenisation theory for curl-div-systems

We study the curl-div-system with variable coefficients and a nonlocal homogenisation problem associated with it. Using, in part refining, techniques from nonlocal $H$ -convergence for closed Hilbert complexes, we define the appropriate topology for possibly nonlocal and non-periodic coefficients in curl-div systems to model highly oscillatory behaviour of the coefficients on small scales. We address curl-div systems under various boundary conditions and analyse the limit of the ratio of small scale over large scale tending to zero. Already for standard Dirichlet boundary conditions and local coefficients the limit system is nontrivial and unexpected. Furthermore, we provide an analysis of highly oscillatory local coefficients for a curl-div system with impedance type boundary conditions relevant in scattering theory for Maxwell's equations and relate the abstract findings to local $H$ -convergence and weak $*$ -convergence of the coefficients.

On a Class of Degenerate Abstract Parabolic Problems and Applications to Some Eddy Current Models

We present an abstract framework for parabolic type equations which possibly degenerate on certain spatial regions. The degeneracies are such that the equations under investigation may admit a type change ranging from parabolic to elliptic type problems. The approach is an adaptation of the concept of so-called evolutionary equations in Hilbert spaces and is eventually applied to a degenerate eddy current type model. The functional analytic setting requires quite minimal assumptions on the boundary and interface regularity. The degenerate eddy current model is justified as a limit model of non-degenerate hyperbolic models of Maxwell's equations.

The Index of Some Mixed Order Dirac-Type Operators and Generalised Dirichlet-Neumann Tensor Fields

We revisit a construction principle of Fredholm operators using Hilbert complexes of densely defined, closed linear operators and apply this to particular choices of differential operators. The resulting index is then computed with the help of explicitly describing the dimension of the cohomology groups of generalised (`harmonic') Dirichlet and Neumann tensor fields. The main results of this contribution are the computation of the indices of Dirac-type operators associated to the elasticity complex and the newly found biharmonic complex, relevant for the biharmonic equation, elasticity, and for the theory of general relativity. The differential operators are of mixed order and cannot be seen as leading order type with relatively compact perturbation. As a by-product we present a comprehensive description of the underlying generalised Dirichlet-Neumann vector and tensor fields defining the respective cohomology groups, including an explicit construction of bases in terms of topological invariants, which are of both analytical and numerical interest. Though being defined by certain projection mechanisms, we shall present a way of computing these basis functions by solving certain PDEs given in variational form. For all of this we rephrase core arguments in the work of Rainer Picard [1982] applied to the de Rham complex and use them as a blueprint for the more involved cases presented here. In passing, we also provide new vector-analytical estimates of generalised Poincare-Friedrichs type useful for elasticity or the theory of general relativity.

preprint2018arXiv

Some remarks on the notions of boundary systems and boundary triple(t)s

In this note we show that if a boundary system in the sense of (Schubert et al. 2015) gives rise to any skew-self-adjoint extension, then it induces a boundary triplet and the classification of all extensions given by (Schubert et al. 2015) coincides with the skew-symmetric version of the classical characterization due to (Gorbachuk et al. 1991). On the other hand we show that for every skew-symmetric operator there is a natural boundary system which leads to an explicit description of at least one maximal dissipative extension. This is in particular also valid in the case that no boundary triplet exists for this operator.

A Solution Theory for a General Class of SPDEs

In this article we present a way of treating stochastic partial differential equations with multiplicative noise by rewriting them as stochastically perturbed evolutionary equations in the sense of \cite{picardbook}, where a general solution theory for deterministic evolutionary equations has been developed. This allows us to present a unified solution theory for a general class of SPDEs which we believe has great potential for further generalizations. We will show that many standard stochastic PDEs fit into this class as well as many other SPDEs such as the stochastic Maxwell equation and time-fractional stochastic PDEs with multiplicative noise on sub-domains of $\Rd$. The approach is in spirit similar to the approach in \cite{dapratozabczyk}, but complementing it in the sense that it does not involve semi-group theory and allows for an effective treatment of coupled systems of SPDEs. In particular, the existence of a (regular) fundamental solution or Green's function is not required.

Continuous dependence on the coefficients for a class of non-autonomous evolutionary equations

The continuous dependence of solutions to certain (non-autonomous, partial, integro-differential-algebraic, evolutionary) equations on the coefficients is addressed. We give criteria that guarantee that convergence of the coefficients in the weak operator topology implies weak convergence of the respective solutions. We treat three examples: A homogenization problem for a Kelvin-Voigt model for elasticity, the discussion of continuous dependence of the coefficients for acoustic waves with impedance type boundary conditions and a singular perturbation problem for a mixed type equation. By means of counter examples we show optimality of the results obtained.

Numerical methods for changing type systems

In this note we develop a numerical method for partial differential equations with changing type. Our method is based on a unified solution theory found by Rainer Picard for several linear equations from mathematical physics. Parallel to the solution theory already developed, we frame our numerical method in a discontinuous Galerkin approach in space-time with certain exponentially weighted spaces.

On Maximal Regularity for a Class of Evolutionary Equations

The issue of so-called maximal regularity is discussed within a Hilbert space framework for a class of evolutionary equations. Viewing evolutionary equations as a sums of two unbounded operators, showing maximal regularity amounts to establishing that the operator sum considered with its natural domain is already closed. For this we use structural constraints of the coefficients rather than semi-group strategies or sesqui-linear form methods, which would be difficult to come by for our general problem class. Our approach, although limited to the Hilbert space case, complements known strategies for approaching maximal regularity and extends them in a different direction. The abstract findings are illustrated by re-considering some known maximal regularity results within the framework presented.

On the continuous dependence on the coefficients of evolutionary equations

In an abstract Hilbert space setting, we discuss many linear phenomena of mathematical physics. The functional analytic framework presented is used to address continuous dependence of the solution operators $\mathcal{S}(\mathcal{M})$ of certain (linear partial differential) equations on the coefficients $\mathcal{M}$ . For this, we introduce a particular class of coefficients $\mathcal{M}$ and study the (nonlinear) mapping $\mathcal{M}\mapsto \mathcal{S}(\mathcal{M})$ . We provide criteria that guarantee the continuity of $\mathcal{S}(\cdot)$ under the norm, the strong, and the weak operator topology. We exemplify our findings in non-autonomous electro-magnetic theory, thermodynamics and acoustics.

On the homogenization of partial integro-differential-algebraic equations

We present a Hilbert space perspective to homogenization of standard linear evolutionary boundary value problems in mathematical physics and provide a unified treatment for (non-)periodic homogenization problems in thermodynamics, elasticity, electro-magnetism and coupled systems thereof. The approach permits the consideration of memory problems as well as differential-algebraic equations. We show that the limit equation is well-posed and causal. We rely on techniques from functional analysis and operator theory only.

On the weighted Gauß--Radau Quadrature

In this short note, we collect some facts on the weighted Gauß--Radau quadrature. In particular, we focus on the location of the Gauß--Radau points being a continuous function of the $L^1$ -weighting function.

On Wellposedness for Some Thermo-Piezo-Electric Coupling Models

There is an increasing reliance on mathematical modelling to assist in the design of piezoelectric ultrasonic transducers since this provides a cost-effective and quick way to arrive at a first prototype. Given a desired operating envelope for the sensor the inverse problem of obtaining the associated design parameters within the model can be considered. It is therefore of practical interest to examine the well-posedness of such models. There is a need to extend the use of such sensors into high temperature environments and so this paper shows, for a broad class of models, the well-posedness of the magneto-electro-thermo-elastic problem. Due to its widespread use in the literature, we also show the well-posedness of the quasi-electrostatic case

Stabilization via Homogenization

In this short note we treat a 1+1-dimensional system of changing type. On different spatial domains the system is of hyperbolic and elliptic type, that is, formally, $\partial_t^2 u_n-\partial_x^2 u_n = \partial_t f$ and $u_n-\partial_x^2 u_n= f$ on the respective spatial domains $\bigcup_{j\in \{1,\ldots,n\}} \big(\frac{j-1}{n},\frac{2j-1}{2n}\big)$ and $\bigcup_{j\in \{1,\ldots,n\}} \big(\frac{2j-1}{2n},\frac{j}{n}\big)$ . We show that $(u_n)_n$ converges weakly to $u$ , which solves the exponentially stable limit equation $\partial_t^2 u +2\partial_t u + u -4\partial_x^2 u = 2(f+\partial_t f)$ on $[0,1]$ . If the elliptic equation is replaced by a parabolic one, the limit equation is \emph{not} exponentially stable.

The Callias Index Formula Revisited

We revisit the Callias index formula for Dirac-type operators $L$ in odd space dimension $n$ , and prove that \begin{align} \text{ind} \, (L) =\bigg(\frac{i}{8π}\bigg)^{\frac{n-1}{2}}\frac{1}{2(\frac{n-1}{2})!} \lim_{Λ\to\infty}\frac{1}{Λ}\sum_{i_{1},\dots,i_{n} = 1}^n \varepsilon_{i_{1}\dots i_{n}} \int_{ΛS_{n-1}}\text{tr}_{\mathbb{C}^d}\, (U(x)(\partial_{i_{1}}U)(x)\dots (\partial_{i_{n-1}}U)(x)) x_{i_{n}}\, d^{n-1} σ(x), \, (*) \end{align} where $U(x) = \text{sgn} \,(Φ(x))$ and $L$ in $L^{2}(\mathbb{R}^{n})^{2^{\widehat n}d}$ is of the form $L= \mathcal{Q} + Φ,$ where $\mathcal{Q} = \bigg(\sum_{j=1}^{n}γ_{j,n}\partial_{j}\bigg) I_d,$ with $γ_{j,n}$ elements of the Euclidean Dirac algebra, and $n=2{\widehat n}$ or $n=2{\widehat n}+1$ . Here $Φ$ is assumed to satisfy the following conditions: \begin{align} & Φ\in C_{b}^{2}\big(\mathbb{R}^{n};\mathbb{C}^{d\times d}\big), \quad d \in \mathbb{N}, \\ & Φ(x)=Φ(x)^{*}, \end{align} there exists $c>0$ , $R\geq0$ such that $|Φ(x)|\geq c I_d, \quad x\in\mathbb{R}^{n}\backslash B(0,R),$ and there exists $\varepsilon> 1/2$ such that for all $α\in\mathbb{N}_{0}^{n}$ , $|α|<3$ , there is $κ>0$ such that $\|(\partial^αΦ)(x)\|\leq \begin{cases} κ(1+|x|)^{-1}, & |α|=1,\\ κ(1+ |x|)^{-1-\varepsilon}, & |α|=2, \end{cases}\quad x\in\mathbb{R}^{n}.$ These conditions on $Φ$ render $L$ a Fredholm operator, and appear to be the most general conditions known to date for which Callias' index formula has been derived. Generalizations of the index formula $(*)$ to certain classes of non-Fredholm operators $L$ invoking the (generalized) Witten index are also discussed.

A Note on a Two-Temperature Model in Linear Thermoelasticity

We discuss the so-called two-temperature model in linear thermoelasticity and provide a Hilbert space framework for proving well-posedness of the equations under consideration. With the abstract perspective of evolutionary equations, the two-temperature model turns out to be a coupled system of the elastic equations and an abstract ode. Following this line of reasoning, we propose another model being entirely an abstract ode. We highlight also an alternative way for a two-temperature model, which might be of independent interest.

A solution decomposition for a singularly perturbed fourth-order problem

We consider a singularly perturbed fourth-order problem with third-order terms on the unit square. With a formal power series approach, we decompose the solution into solutions of reduced (third-order) problems and various layer parts. The existence of unique solutions for the problem itself and for the reduced third-order problems is also addressed.

On Abstract $\mathrm{grad}-\mathrm{div}$ Systems

For a large class of dynamical problems from mathematical physics the skew-selfadjointness of a spatial operator of the form $A=\left(\begin{array}{cc} 0 & -C^{*}\\ C & 0 \end{array}\right)$ , where $C:D\left(C\right)\subseteq H_{0}\to H_{1}$ is a closed densely defined linear operator, is a typical property. Guided by the standard example, where $C=\mathrm{grad}=\left(\begin{array}{c} \partial_{1}\\ \vdots\\ \partial_{n} \end{array}\right)$ (and $-C^{*}=\mathrm{div}$ , subject to suitable boundary constraints), an abstract class of operators $C=\left(\begin{array}{c} C_{1}\\ \vdots\\ C_{n} \end{array}\right)$ is introduced (hence the title). As a particular application we consider a non-standard coupling mechanism and the incorporation of diffusive boundary conditions both modeled by setting associated with a skew-selfadjoint spatial operator $A$ .

Perturbations of positive semigroups on $L_p$ -spaces

We give a characterization of a variation of constants type estimate relating two positive semigroups on (possibly different) $L_p$ -spaces to one another in terms of corresponding estimates for the respective generators and of estimates for the respective resolvents. The results have applications to kernel estimates for semigroups induced by accretive and non-local forms on $σ$ -finite measure spaces.

A differential equation with state-dependent delay from cell population biology

We analyze a differential equation with a state-dependent delay that is implicitly defined via the solution of an ODE. The equation describes an established though little analyzed cell population model. Based on theoretical results of Hartung, Krisztin, Walther and Wu we elaborate conditions for the model ingredients, in particular vital rates, that guarantee the existence of a local semiflow. Here proofs are based on implicit function arguments. To show global existence, we adapt a theorem from a classical book on functional differential equations by Hale and Lunel, which gives conditions under which - if there is no global existence - closed and bounded sets are left for good, to the $C^1$ -topology, which is the natural setting when dealing with state-dependent delays. The proof is based on an older result for semiflows on metric spaces.

A note on causality in Banach spaces

In this note we provide examples that show that a common notion of causality for linear operators on Banach spaces does not carry over to the closure of the respective operators. We provide an alternative definition for causality, which is equivalent to the usual definition for closed linear operators but does carry over to the closure.

Mother Operators and their Descendants

A mechanism deriving new well-posed evolutionary equations from given ones is inspected. It turns out that there is one particular spatial operator from which many of the standard evolutionary problems of mathematical physics can be generated by this abstract mechanism using suitable projections. The complexity of the dynamics of the phenomena considered can be described in terms of suitable material laws. The idea is illustrated with a number of concrete examples.

On a Connection between the Maxwell System, the Extended Maxwell System, the Dirac Operator and Gravito-Electromagnetism

Maxwell's equation, Dirac's equation and the equation of gravito-electromagnetism are shown to be particular instances of the extended Maxwell system. The equations are discussed in the framework of the theory of evolutionary equations. Their formal relationship are systematically analyzed. Applications to coupled systems such as the Maxwell-Dirac system are also discussed.

On Some Models in Linear Thermo-Elasticity with Rational Material Laws

We shall consider some common models in linear thermo-elasticity within a common structural framework. Due to the flexibility of the structural perspective we will obtain well-posedness results for a large class of generalized models allowing for more general material properties such as anisotropies, inhomogeneities, etc.

Well-posedness via Monotonicity. An Overview

The idea of monotonicity (or positive-definiteness in the linear case) is shown to be the central theme of the solution theories associated with problems of mathematical physics. A "grand unified" setting is surveyed covering a comprehensive class of such problems. We elaborate the applicability of our scheme with a number examples. A brief discussion of stability and homogenization issues is also provided.

Boundary systems and (skew-)self-adjoint operators on infinite metric graphs

We generalize the notion of Lagrangian subspaces to self-orthogonal subspaces with respect to a (skew-)symmetric form, thus characterizing (skew-)self-adjoint and unitary operators by means of self-ortho-gonal subspaces. By orthogonality preserving mappings, these characterizations can be transferred to abstract boundary value spaces of (skew-)symmetric operators. Introducing the notion of boundary systems we then present a unified treatment of different versions of boundary triples and related concepts treated in the literature. The application of the abstract results yields a description of all (skew-)self-adjoint realizations of Laplace and first derivative operators on graphs.

G-convergence of linear differential equations

We discuss $G$ -convergence of linear integro-differential-algebaric equations in Hilbert spaces. We show under which assumptions it is generic for the limit equation to exhibit memory effects. Moreover, we investigate which classes of equations are closed under the process of $G$ -convergence. The results have applications to the theory of homogenization. As an example we treat Maxwell's equation with the Drude-Born-Fedorov constitutive relation.

Homogenization in fractional elasticity

In this note we treat the equations of fractional elasticity. After establishing well-posedness, we show a compactness result related to the theory of homogenization. For this, a previous result in (abstract) homogenization theory of evolutionary equations has to be improved. The approach also permits the consideration of non-local operators (in time and space).

On a Comprehensive Class of Linear Control Problems

We discuss a class of linear control problems in a Hilbert space setting. This class encompasses such diverse systems as port-Hamiltonian systems, Maxwell's equations with boundary control or the acoustic equations with boundary control and boundary observation. The boundary control and observation acts on abstract boundary data spaces such that the only geometric constraint on the underlying domain stems from requiring a closed range constraint for the spatial operator part, a requirement which for the wave equation amounts to the validity of a Poincare-Wirtinger-type inequality. We also address the issue of conservativity of the control problems under consideration.

On Evolutionary Equations with Material Laws Containing Fractional Integrals

A well-posedness result for a time-shift invariant class of evolutionary operator equations involving material laws with fractional time-integrals of order $α\in]0,1[$ is considered and exemplified by an application to a Kelvin-Voigt type model.

On Non-Autonomous Evolutionary Problems

The paper extends well-posedness results of a previously explored class of time-shift invariant evolutionary problems to the case of non-autonomous media. The Hilbert space setting developed for the time-shift invariant case can be utilized to obtain an elementary approach to non-autonomous equations. The results cover a large class of evolutionary equations, where well-known strategies like evolution families may be difficult to use or fail to work. We exemplify the approach with an application to a Kelvin-Voigt-type model for visco-elastic solids.

On Non-Autonomous Integro-Differential-Algebraic Evolutionary Problems

In this article, we show that a technique for showing well-posedness results for evolutionary equations in the sense of [13] established in [16] applies to a broader class of non-autonomous integro-differential-algebraic equations. Using the concept of evolutionary mappings we prove that the respective solution operators do not depend on certain parameters describing the underlying spaces in which the well-posedness results are established.

On Some Models for Elastic Solids with Micro-Structure

We review the concept of well-posedness in the context of evolutionary problems from mathematical physics for a particular subclass of problems from elasticity theory. The complexity of physical phenomena appears as encoded in so called material laws. The usefulness of the structural perspective developed is illustrated by showing that many initial boundary value problems in the theory of elastic solids share the same type of solution theory. Moreover, interconnections of the respective models are discussed via a previously introduced mother/descendant mechanism.

A Functional Analytic Perspective to Delay Differential Equations

We generalize the solution theory for a class of delay type differential equations developed in a previous paper, dealing with the Hilbert space case, to a Banach space setting. The key idea is to consider differentiation as an operator with the whole real line as the underlying domain as a means to incorporate pre-history data. We focus our attention on the issue of causality of the differential equations as a characterizing feature of evolutionary problems and discuss various examples. The arguments mainly rely on a variant of the contraction mapping theorem and a few well-known facts from functional analysis.

A Hilbert Space Perspective on Ordinary Differential Equations with Memory Term

We discuss ordinary differential equations with delay and memory terms in Hilbert spaces. By introducing a time derivative as a normal operator in an appropriate Hilbert space, we develop a new approach to a solution theory covering integro-differential equations, neutral differential equations and general delay differential equations within a unified framework. We show that reasonable differential equations lead to causal solution operators.

A note on elliptic type boundary value problems with maximal monotone relations

In this note we discuss an abstract framework for standard boundary value problems in divergence form with maximal monotone relations as "coefficients". A reformulation of the respective problems is constructed such that they turn out to be unitary equivalent to inverting a maximal monotone relation in a Hilbert space. The method is based on the idea of "tailor-made" distributions as provided by the construction of extrapolation spaces, see e.g. [Picard, McGhee: Partial Differential Equations: A unified Hilbert Space Approach. DeGruyter, 2011]. The abstract framework is illustrated by various examples.