The subject of the thesis is the investigation of asymptotic properties of solutions of the Cauchy problem for structurally damped sigma-evolution operators with time dependent, monotonous, dissipation term. An appropriate energy for solutions of the sigma-evolution equations is defined and some estimates for energies of higher order are proved. In the scale invariant case the optimality of these estimates is shown. Further, the influence of properties of the time dependent dissipation on L^p-L^q estimates for the energy with p and q bigger or equal to 2 and from the conjugate line is clarified. Also smoothing properties of the operators under consideration are investigated. The connection between the regularity of the data and the regularity of the solution in terms of L^2 based Gevrey spaces is considered. Finally, L^1-L^1-estimates in the special case delta = sigma/2 and decreasing dissipative coefficient. / Thema der vorliegenden Dissertation ist die Untersuchung asymptotischer Eigenschaften von Lösungen des Cauchy Problems für strukturell gedämpfte sigma-Evolutions-Operatoren mit zeitabhängigem, monotonen Dissipationskoeffizienten. Es wird eine geeignete Energie definiert und für diese Abschätzungen, auf für entsprechende Energien höherer Ordnung gezeigt. Darüber hinaus wird der Einfluss des Dissipationskoeffizienten auf L^p-L^q Abschätzungen auf und entfernt von der konjugierten Linie untersucht. Im skaleninvarianten Fall wird die Schärfe der Abschätzungen bewiesen. Weiterhin wird der Zusammenhang zwischen der Regularität der Daten und der der Lösung in Termen von L^2-basierten Gevrey-Räumen untersucht. Schließlich werden L^1-L^1-Abschätzungen für den Spezialfall delta = sigma/2 und monoton fallenden Dissipationskoeffizienten gezeigt.

The thesis consists of two parts, a theoretical part and an applied part, and in addition an Appendix. Except for a very short chapter in the applied part and the appendix we only present previously unknown results leading to a very concise style. In the theoretical part we study rates of decay for vector-valued functions and semigroups of operators depending on a real and positive variable. Under boundedness assumptions on the function/semigroup itself and under analytic extendability assumptions of its Laplace transform/resolvent across the imaginary axis we provide (almost) sharp rates of decay. Our results improve known results in this very active area of research. In the second part of the thesis we apply our results to specific examples (from the field of PDEs): local energy decay for wave equations on exterior domains, energy decay for damped wave equations on bounded domains and decay for a viscoelastic boundary damping model for sound waves. Many more examples can be found in the vast literature.:Part 1 Quantified Tauberian theorems and decay of C0-semigroups 1 Decay of vector-valued functions 2 Optimal decay for C0-semigroups on Hilbert spaces Part 2 Applications: decay of waves 3 Local decay for waves in exterior domains 4 Waves on a square with constant damping on a strip 5 A viscoelastic boundary damping model

Lineare hyperbolische partielle Differentialgleichungen in homogenen Medien, beispielsweise die Wellengleichung, die die Ausbreitung und die Streuung akustischer Wellen beschreibt, können im Zeitbereich mit Hilfe von Randintegralgleichungen formuliert werden. Im ersten Hauptteil dieser Arbeit stellen wir eine effiziente Möglichkeit vor, numerische Approximationen solcher Gleichungen zu implementieren, wenn das Huygens-Prinzip nicht gilt. Wir nutzen die Faltungsquadraturmethode für die Zeitdiskretisierung und eine Galerkin-Randelement-Methode für die Raumdiskretisierung. Mit der Faltungsquadraturmethode geht eine diskrete Faltung der Faltungsgewichte mit der Randdichte einher. Bei Gültigkeit des Huygens-Prinzips konvergieren die Gewichte exponentiell gegen null, sofern der Index hinreichend groß ist. Im gegenteiligen Fall, das heißt bei geraden Raumdimensionen oder wenn Dämpfungseffekte auftreten, kann kein Verschwinden der Gewichte beobachtet werden. Das führt zu Schwierigkeiten bei der effizienten numerischen Behandlung. Im ersten Hauptteil dieser Arbeit zeigen wir, dass die Kerne der Faltungsgewichte in gewisser Weise die Fundamentallösung im Zeitbereich approximieren und dass dies auch zutrifft, wenn beide bezüglich der räumlichen Variablen abgeleitet werden. Da die Fundamentallösung zudem für genügend große Zeiten, etwa nachdem die Wellenfront vorbeigezogen ist, glatt ist, schließen wir Gleiches auch in Bezug auf die Faltungsgewichte, die wir folglich mit hoher Genauigkeit und wenigen Interpolationspunkten interpolieren können. Darüber hinaus weisen wir darauf hin, dass zur weiteren Einsparung von Speicherkapazitäten, insbesondere bei Langzeitexperimenten, der von Schädle et al. entwickelte schnelle Faltungsalgorithmus eingesetzt werden kann. Wir diskutieren eine effiziente Implementierung des Problems und zeigen Ergebnisse eines numerischen Langzeitexperimentes. Im zweiten Hauptteil dieser Arbeit beschäftigen wir uns mit Transmissionsproblemen der Wellengleichung im Freiraum. Solche Probleme werden gewöhnlich derart behandelt, dass der Freiraum, wenn nötig durch Einführen eines künstlichen Randes, in ein unbeschränktes Außengebiet und ein beschränktes Innengebiet geteilt wird mit dem Ziel, eventuelle Inhomogenitäten oder Nichtlinearitäten des Materials vollständig im Innengebiet zu konzentrieren. Wir werden eine Lösungsstrategie vorstellen, die es erlaubt, die aus der Teilung resultierenden Teilprobleme so weit wie möglich unabhängig voneinander zu behandeln. Die Kopplung der Teilprobleme erfolgt über Transmissionsbedingungen, die auf dem ihnen gemeinsamen Rand vorgegeben sind. Wir diskutieren ein Kopplungsverfahren, das auf verschiedene Diskretisierungsschemata für das Innen- und das Außengebiet zurückgreift. Wir werden insbesondere ein explizites Verfahren im Innengebiet einsetzen, im Gegensatz zum Außengebiet, bei dem wir ein auf ein Mehrschrittverfahren beruhendes Faltungsquadraturverfahren nutzen. Die Kopplung erfolgt nach der Strategie von Johnson und Nédélec, bei der die direkte Randintegralmethode zum Einsatz kommt. Diese Strategie führt auf ein unsymmetrische System. Wir analysieren das diskrete Problem hinsichtlich Stabilität und Konvergenz und unterstreichen die Einsatzfähigkeit des Kopplungsalgorithmus mit der Durchführung numerischer Experimente.

In this thesis we study the global existence of small data solutions to the Cauchy problem for semilinear damped wave equations with an effective dissipation term, where the data are supposed to belong to different classes of regularity. We apply these results to the Cauchy problem for weakly coupled systems of semilinear effectively damped waves with respect to the defined classes of regularity for different power nonlinearities. We also presented blow-up results for semi-linear systems with weakly coupled damped waves.

The PhD thesis deals with global in time existence results and blow-up result for a semilinear wave model with scale-invariant damping and mass. Since the time-dependent coefficients for the considered model make somehow the damping and the mass a threshold term between effective and non-effective terms, it turns out that a fundamental role in the description of qualitative properties of solutions to this semilinear model and to the corresponding linear homogeneous Cauchy problem is played by the multiplicative constants appearing in those coefficients. For coefficients that make the damping term dominant, we can use the standard approach for the classical damped wave model with L^2 − L^2 estimates and the so-called test function method. On the other hand, when the interaction among those coefficients is balanced, then, it is possible to observe how typical tools for hyperbolic models, as for example Kato’s lemma, provide sharp global in time existence results and sharp blow-up results for super- and sub-Strauss type exponents, respectively.

No description available.

In this thesis, we are interested in damped wave models with time-dependent propagation speed and time-dependent damping term both having a time-dependent oscillation term. The main goal of this thesis is to understand the influence of strong time-dependent oscillations on Sobolev solutions to the linear models and consequently, to the semilinear models. Especially, due to the deteriorating influence of oscillations on solutions, a stabilization condition and higher-order regularity of the time-dependent coefficients may compensate 'bad behaviors' arising from oscillations.:1. Introduction 2. The influence of oscillations on linear damped wave equation with time-dependent coefficients 3. Global in time existence results for damped wave models with power nonlinearity 4. Global in time existence results for damped wave models with different power nonlinearities 5. Lp-Lq estimates for wave equations with strong time-dependent oscillations 6. Further research topics A. Basic tools B. List of symbols and abbreviations Bibliography

This work addresses the question of the efficient numerical solution of time-domain boundary integral equations with retarded potentials arising in the problems of acoustic and electromagnetic scattering. The convolutional form of the time-domain boundary operators allows to discretize them with the help of Runge-Kutta convolution quadrature. This method combines Laplace-transform and time-stepping approaches and requires the explicit form of the fundamental solution only in the Laplace domain to be known. Recent numerical and analytical studies revealed excellent properties of Runge-Kutta convolution quadrature, e.g. high convergence order, stability, low dissipation and dispersion. As a model problem, we consider the wave scattering in three dimensions. The convolution quadrature discretization of the indirect formulation for the three-dimensional wave equation leads to the lower triangular Toeplitz system of equations. Each entry of this system is a boundary integral operator with a kernel defined by convolution quadrature. In this work we develop an efficient method of almost linear complexity for the solution of this system based on the existing recursive algorithm. The latter requires the construction of many discretizations of the Helmholtz boundary single layer operator for a wide range of complex wavenumbers. This leads to two main problems: the need to construct many dense matrices and to evaluate many singular and near-singular integrals. The first problem is overcome by the use of data-sparse techniques, namely, the high-frequency fast multipole method (HF FMM) and H-matrices. The applicability of both techniques for the discretization of the Helmholtz boundary single-layer operators with complex wavenumbers is analyzed. It is shown that the presence of decay can favorably affect the length of the fast multipole expansions and thus reduce the matrix-vector multiplication times. The performance of H-matrices and the HF FMM is compared for a range of complex wavenumbers, and the strategy to choose between two techniques is suggested. The second problem, namely, the assembly of many singular and nearly-singular integrals, is solved by the use of the Huygens principle. In this work we prove that kernels of the boundary integral operators $w_n^h(d)$ ($h$ is the time step and $t_n=nh$ is the time) exhibit exponential decay outside of the neighborhood of $d=nh$ (this is the consequence of the Huygens principle). The size of the support of these kernels for fixed $h$ increases with $n$ as $n^a,a<1$, where $a$ depends on the order of the Runge-Kutta method and is (typically) smaller for Runge-Kutta methods of higher order. Numerical experiments demonstrate that theoretically predicted values of $a$ are quite close to optimal. In the work it is shown how this property can be used in the recursive algorithm to construct only a few matrices with the near-field, while for the rest of the matrices the far-field only is assembled. The resulting method allows to solve the three-dimensional wave scattering problem with asymptotically almost linear complexity. The efficiency of the approach is confirmed by extensive numerical experiments.

The subject of the thesis is the investigation of asymptotic properties of solutions of the Cauchy problem for structurally damped sigma-evolution operators with time dependent, monotonous, dissipation term. An appropriate energy for solutions of the sigma-evolution equations is defined and some estimates for energies of higher order are proved. In the scale invariant case the optimality of these estimates is shown. Further, the influence of properties of the time dependent dissipation on L^p-L^q estimates for the energy with p and q bigger or equal to 2 and from the conjugate line is clarified. Also smoothing properties of the operators under consideration are investigated. The connection between the regularity of the data and the regularity of the solution in terms of L^2 based Gevrey spaces is considered. Finally, L^1-L^1-estimates in the special case delta = sigma/2 and decreasing dissipative coefficient. / Thema der vorliegenden Dissertation ist die Untersuchung asymptotischer Eigenschaften von Lösungen des Cauchy Problems für strukturell gedämpfte sigma-Evolutions-Operatoren mit zeitabhängigem, monotonen Dissipationskoeffizienten. Es wird eine geeignete Energie definiert und für diese Abschätzungen, auf für entsprechende Energien höherer Ordnung gezeigt. Darüber hinaus wird der Einfluss des Dissipationskoeffizienten auf L^p-L^q Abschätzungen auf und entfernt von der konjugierten Linie untersucht. Im skaleninvarianten Fall wird die Schärfe der Abschätzungen bewiesen. Weiterhin wird der Zusammenhang zwischen der Regularität der Daten und der der Lösung in Termen von L^2-basierten Gevrey-Räumen untersucht. Schließlich werden L^1-L^1-Abschätzungen für den Spezialfall delta = sigma/2 und monoton fallenden Dissipationskoeffizienten gezeigt.

This work addresses the question of the efficient numerical solution of time-domain boundary integral equations with retarded potentials arising in the problems of acoustic and electromagnetic scattering. The convolutional form of the time-domain boundary operators allows to discretize them with the help of Runge-Kutta convolution quadrature. This method combines Laplace-transform and time-stepping approaches and requires the explicit form of the fundamental solution only in the Laplace domain to be known. Recent numerical and analytical studies revealed excellent properties of Runge-Kutta convolution quadrature, e.g. high convergence order, stability, low dissipation and dispersion. As a model problem, we consider the wave scattering in three dimensions. The convolution quadrature discretization of the indirect formulation for the three-dimensional wave equation leads to the lower triangular Toeplitz system of equations. Each entry of this system is a boundary integral operator with a kernel defined by convolution quadrature. In this work we develop an efficient method of almost linear complexity for the solution of this system based on the existing recursive algorithm. The latter requires the construction of many discretizations of the Helmholtz boundary single layer operator for a wide range of complex wavenumbers. This leads to two main problems: the need to construct many dense matrices and to evaluate many singular and near-singular integrals. The first problem is overcome by the use of data-sparse techniques, namely, the high-frequency fast multipole method (HF FMM) and H-matrices. The applicability of both techniques for the discretization of the Helmholtz boundary single-layer operators with complex wavenumbers is analyzed. It is shown that the presence of decay can favorably affect the length of the fast multipole expansions and thus reduce the matrix-vector multiplication times. The performance of H-matrices and the HF FMM is compared for a range of complex wavenumbers, and the strategy to choose between two techniques is suggested. The second problem, namely, the assembly of many singular and nearly-singular integrals, is solved by the use of the Huygens principle. In this work we prove that kernels of the boundary integral operators $w_n^h(d)$ ($h$ is the time step and $t_n=nh$ is the time) exhibit exponential decay outside of the neighborhood of $d=nh$ (this is the consequence of the Huygens principle). The size of the support of these kernels for fixed $h$ increases with $n$ as $n^a,a<1$, where $a$ depends on the order of the Runge-Kutta method and is (typically) smaller for Runge-Kutta methods of higher order. Numerical experiments demonstrate that theoretically predicted values of $a$ are quite close to optimal. In the work it is shown how this property can be used in the recursive algorithm to construct only a few matrices with the near-field, while for the rest of the matrices the far-field only is assembled. The resulting method allows to solve the three-dimensional wave scattering problem with asymptotically almost linear complexity. The efficiency of the approach is confirmed by extensive numerical experiments.