Trimester Seminar

Venue: HIM lecture hall, Poppelsdorfer Allee 45
Organizers: Susanne C. Brenner, Björn Engquist, Max Gunzburger, Daniel Peterseim, Marc Alexander Schweitzer

Tuesday, February 21

15:00 - 16:00 Sören Bartels: Stable discretization of singular flows and application to a problem in optimal insulation

Thursday, February 2

15:00 - 16:00 Ivan Oseledets: Topological optimization and aposteriori estimates: is it possible to compute the constant?

Abstract: Joint work with George Ovchinnikov, Vlad Pimanov, Denis Zorin.

A major problem in topology optimization in thermal or elasticity problems is the emergence of so-called checkerboard structures, where the minimum of the discretized FEM functional does not correspond to a nice value of the original functional. Or it maybe, and the true solution is a certain microstructure. A typical approach is either to do filtering of the topology to remove bad configurations, or introduce penalty terms for them (for example, by penalizing the length of the boundary of the obtaned shape). We propose a fresh view on this problem, based on a classical connection between the exact solution and the solution of the discretized version, which leads to a computable upper bound of the functional. By minimizing this upper bound we get a guarantee for the effectiveness of the computed topology. Another challenge is that to get a final number, we do not only need the aposteriori error estimate, but we need the exact value of the constant. We provide a series of different solutions obtained by the optimization algorithm for different empirical values of the constant, but we are not yet able to numerically compute the true desired values of the functionals for such structures, thus leaving the main question about the necessity of the microstructures in the optimal topology still open.

Tuesday, January 31

15:00 - 16:00 Ivan Graham: Fast UQ for diffusion and neutron transport

Abstract: This is joint work with F. Kuo and I.H. Sloan (UNSW), D. Nuyens (Leuven) and M. Parkinson and R. Scheichl (Bath).

We shall describe some practical uncertainty quantification problems involving the diffusion equation and the Boltzmann transport equation where the uncertainty is described by stationary Gaussian random fields. Small length scale and high variance presents PDEs problems with a challenging multiscale structure. We describe a UQ algorithm for these applications which combines circulant embedding techniques for the sampling with Quasi Monte-Carlo methods for computing the required high-dimensional integrals. The method is capable of handling very high stochastic dimension and is consistently faster than Monte-Carlo methods. Multilevel variants can be used to obtain further acceleration. We describe some recent theory and computations.

Thursday, January 26

14:30 - 16:00 Wolfgang Hackbusch: Numerical Tensor Calculus

Tuesday, January 24

15:00 - 16:00 Donald L. Brown: Multiscale Methods and Stability in Some High-Frequency Helmholtz-Type Problems

Abstract: In this talk, we will discuss issues in the computation of high-frequency Helmholtz-type problems. In particular, we discuss the issue of pollution effects and how certain multiscale LOD correction methods can eliminate the effect in certain resolution regimes. This will help to motivate the issue of frequency dependent stability in continuous problems. We then present a few of the current results on frequency explicit bounds for heterogenous acoustic and elastic media.

Thursday, January 19

15:00 - 16:00 Robert Scheichl: Multilevel Quasi Monte Carlo Methods for Uncertainty Quantification

Abstract: Large-scale PDE-constrained uncertainty propagation and Bayesian inference are inherently difficult and computationally intensive problems. Sampling methods are among the most accurate and promising methods for these problems. Their cost is dimension independent, thus rendering them very suitable for the infinite-dimensional PDE setting. However, classical sampling approaches, such as Monte Carlo methods or Metropolis-Hastings MCMC, are very slow to converge and often infeasible for realistic applications. The problem becomes even more involved when the differential operator depends in a non-affine way on the random parameters and when the quantities of interest depend nonlinearly on the PDE solution. A popular model problem for this situation is the lognormal diffusion problem which is of actual practical interest in hydrology. Quasi-Monte Carlo methods and their multilevel extensions provide alternatives with potentially vastly improved computational efficiencies, in terms of cost to accuracy. We provide a dimension-independent convergence analysis for the lognormal problem for the forward uncertainty propagation case, as well as for the inverse Bayesian inference setting supported by numerical experiments.

Wednesday, January 18

15:00 - 16:00 Alexandre Ern: Finite element quasi-interpolation and best approximation

Abstract: We introduce quasi-interpolation operators for scalar- and vector-valued finite element spaces with some continuity across mesh interfaces. These operators are stable in L1, are projections, and deliver optimal local approximation estimates in Sobolev spaces. The theory is illustrated on H1-, H(curl)- and H(div)-conforming finite element spaces.

Tuesday, January 17

15:00 - 16:00 Chi-Wang Shu: Discontinuous Galerkin Finite Element Method for Multiscale Problems

Abstract: In this talk, we first give a brief introduction to the discontinuous Galerkin method, which is a finite element method using completely discontinuous basis functions, for solving hyperbolic conservation laws and parabolic and elliptic equations. We will then survey the progress in developing discontinuous Galerkin methods for multiscale problems, in three different approaches, namely using the heterogeneous multiscale method (HMM) framework, using domain decompositions, and using multiscale basis in the discontinuous Galerkin method. Numerical results will be shown to demonstrate the effectiveness of the multiscale discontinuous Galerkin methods.

Friday, January 6

14:15 - 15:00 Christian Stohrer: Finite Element Heterogeneous Multiscale Methods for Maxwell's Equations

Abstract: To approximate the effective behavior an electromagnetic wave propagating through a multiscale medium, we adapt the Finite Element Heterogeneous Multiscale Method to Maxwell's equations. The proposed method can be applied to time-harmonic and time dependent Maxwell's equations. In the talk we show how a priori-error estimates can be proven using the notion of T-coercivity for time-harmonic Maxwell's equations and a Strang-type lemma for the time-dependent case.

15:00 - 16:00 Michael Feischl: Fast random field generation with H-matrices

Abstract: We use the H-matrix technology to compute the approximate square root of a covariance matrix in linear complexity. This allows us to generate normal and log-normal random fields on general finite point sets with optimal complexity. We derive rigorous error estimates, which show that particularly for short correlation length, this new method outperforms the standard method of truncating the Karhunen-Loève expansion.

Besides getting rid of the log-factor in the complexity estimate, the method requires only mild assumptions on the covariance function and on the point set. Therefore, it might also be an alternative to circulant embedding, which is well-defined only for regular grids and stationary covariance functions.