Where: Kirwan Hall 3206

Speaker: Francisco-Javier Sayas (Department of Mathematical Sciences, University of Delaware) - http://www.math.udel.edu/~fjsayas/

Abstract: I will explain ongoing work with my team (Tom Brown, Shukai Du, and Hasan Eruslu) on Finite Element simulation of elastic wave propagation in media that are modeled using a strain-to-stress relation that keeps track of the past evolution of the solid. The family of models I will discuss includes the classical (Zener) differential viscoelastic model, a fractional derivative version thereof, and combinations of the above with pure elastic behavior. The analysis will be presented using transfer function techniques, but I will explain how in some cases refined results can be found using techniques of semigroup theory.

Where: Kirwan Hall 3206

Speaker: Abner J. Salgado (University of Tennessee, Knoxville) - http://www.math.utk.edu/~abnersg/

Abstract: We propose and analyze a two-scale finite element method for the Isaacs equation. By showing the consistency of the approximation and that the method satisfies the discrete maximum principle we establish convergence to the viscosity solution. By properly choosing each of the scales, and using the recently derived discrete Alexandrov Bakelman Pucci estimate we can deduce rates of convergence.

Where: Kirwan Hall 3206

Speaker: Franziska Weber (Department of Mathematics, University of Maryland, College Park) - https://terpconnect.umd.edu/~frweber/

Abstract: We present a convergent finite difference method for approximating wave maps into the sphere. The method is based on a reformulation of the second order wave map equation as a first order system by using the angular momentum as an auxiliary variable. This enables us to preserve the length constraint as well as the energy inherit in the system of equations at the discrete level. The method is shown to converge to a weak solution of the wave map equation as the discretization parameters go to zero. Moreover, it is fast in the sense that O(N log N) operations are required in each time step (where N is the number of grid cells) and a linear CFL-condition is sufficient for stability and convergence. The performance of the method is illustrated by numerical experiments.

The method can be extended to a convergent scheme for the damped wave map equation and the heat map flow. If time permits, I will also discuss possible extensions of the method to applications for liquid crystal dynamics.

Where: Kirwan Hall 3206

Speaker: Rob Stevenson (University of Amsterdam) - https://staff.fnwi.uva.nl/r.p.stevenson/

Abstract: We consider a Fictitious Domain formulation of an elliptic PDE, and solve the arising saddle-point problem by an inexact preconditioned Uzawa iteration.

Solving the arising `inner' elliptic problems with an adaptive finite element method, we prove that the overall method converges with the best possible rate.

So far our results apply to two-dimensional domains and lowest order finite elements (continuous piecewise linears on the fictitious domain, and piecewise constants on the boundary of the physical domain).

Joint work with S. Berrone (Torino), A. Bonito (Texas A&M), and M. Verani (Milano).

Where: Kirwan Hall 3206

Speaker: Matthias Maier (University of Minnesota) - http://www-users.math.umn.edu/~msmaier/

Abstract: In the terahertz frequency range, the effective (complex-valued) surface conductivity of atomically thick 2D materials such as graphene has a positive imaginary part that is considerably larger than the real part. This feature allows for the propagation of slowly decaying electromagnetic waves, called surface plasmon-polaritons (SPPs), that are confined near the material interface with wavelengths much shorter than the wavelength of the free-space radiation. SPPs are promising ingredients in the design of novel optical applications promising "subwavelength optics" beyond the diffraction limit. There is a compelling need for controllable numerical schemes which, placed on firm mathematical grounds, can reliably describe SPPs in a variety of geometries.

In this talk we present a higher-order finite element approach for the simulation of SPP structures on a conducting sheet, excited by a plane-wave or electric Hertzian dipole sources. Aspects of the numerical treatment such as absorbing, perfectly matched layers, local refinement and a-posteriori error control are discussed. Corresponding analytical results are briefly presented as well.

Where: Kirwan Hall 3206

Speaker: Shawn Walker (Louisiana State University) - https://www.math.lsu.edu/~walker/

Where: Kirwan Hall 3206

Speaker: Jie Shen (Purdue University) - https://www.math.purdue.edu/~shen/

Where: Kirwan Hall 3206

Speaker: Christian Glusa (Sandia National Laboratories) -

Abstract: We explore the connection between fractional order partial differential

equations in two or more spatial dimensions with boundary integral

operators to develop techniques that enable one to efficiently tackle

the integral fractional Laplacian. We develop all of the components

needed to construct an adaptive finite element code that can be used to

approximate fractional partial differential equations, on non-trivial

domains in \(d\geq 1\) dimensions. Our main approach consists of taking

tools that have been shown to be effective for adaptive boundary element

methods and, where necessary, modifying them so that they can be applied

to the fractional PDE case. Improved a priori error estimates are

derived for the case of quasi-uniform meshes which are seen to deliver

sub-optimal rates of convergence owing to the presence of singularities.

Attention is then turned to the development of an a posteriori error

estimate and error indicators which are suitable for driving an adaptive

refinement procedure. We assume that the resulting refined meshes are

locally quasi-uniform and develop efficient methods for the assembly of

the resulting linear algebraic systems and their solution using

iterative methods, including the multigrid method. The storage of the

dense matrices along with efficient techniques for computing the dense

matrix vector products needed for the iterative solution is also

considered. Importantly, the approximation does not make any strong

assumptions on the shape of the underlying domain and does not rely on

any special structure of the matrix that could be exploited by fast

transforms. The performance and efficiency of the resulting algorithm is

illustrated for a variety of examples.

This is joint work with Mark Ainsworth, Brown University.

Where: Kirwan Hall 3206

Speaker: Simon Foucart (Texas A&M University) - http://www.math.tamu.edu/~foucart/

Where: Kirwan Hall 3206

Speaker: Alex Townsend (Cornell University) - http://www.math.cornell.edu/~ajt/

Where: Kirwan Hall 3206

Speaker: Mert Gurbuzbalaban (Rutgers University) - https://mert.lids.mit.edu