Computational Mathematics and Applications Seminar

Please note that the list below only shows forthcoming events, which may not include regular events that have not yet been entered for the forthcoming term. Please see the past events page for a list of all seminar series that the department has on offer.

Past events in this series
15 October 2020
14:00
Jim Bremer
Abstract

 

One of the standard methods for the solution of elliptic boundary value problems calls for reformulating them as systems of integral equations.  The integral operators that arise in this fashion typically have singular kernels, and, in many cases of interest, the solutions of these equations are themselves singular.  This makes the accurate discretization of the systems of integral equations arising from elliptic boundary value problems challenging.

Over the last decade, Generalized Gaussian quadrature rules, which are n-point quadrature rules that are exact for a collection of 2n functions, have emerged as one of the most effective tools for discretizing singular integral equations. Among other things, they have been used to accelerate the discretization of singular integral operators on curves, to enable the accurate discretization of singular integral operators on complex surfaces and to greatly reduce the cost of representing the (singular) solutions of integral equations given on planar domains with corners.

We will first briefly outline a standard method for the discretization of integral operators given on curves which is highly amenable to acceleration through generalized Gaussian quadratures. We will then describe a numerical procedure for the construction of Generalized Gaussian quadrature rules.

Much of this is joint work with Zydrunas Gimbutas (NIST Boulder) and Vladimir Rokhlin (Yale University).

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please send email to trefethen@maths.ox.ac.uk.

  • Computational Mathematics and Applications Seminar
29 October 2020
14:00
Kirill Serkh
Abstract

It has long been known that many elliptic partial differential equations can be reformulated as Fredholm integral equations of the second kind on the boundaries of their domains. The kernels of the resulting integral equations are weakly singular, which has historically made their numerical solution somewhat onerous, requiring the construction of detailed and typically sub-optimal quadrature formulas. Recently, a numerical algorithm for constructing generalized Gaussian quadratures was discovered which, given 2n essentially arbitrary functions, constructs a unique n-point quadrature that integrates them to machine precision, solving the longstanding problem posed by singular kernels.

When the domains have corners, the solutions themselves are also singular. In fact, they are known to be representable, to order n, by a linear combination (expansion) of n known singular functions. In order to solve the integral equation accurately, it is necessary to construct a discretization such that the mapping (in the L^2-sense) from the values at the discretization points to the corresponding n expansion coefficients is well-conditioned. In this talk, we present exactly such an algorithm, which is optimal in the sense that, given n essentially arbitrary functions, it produces n discretization points, and for which the resulting interpolation formulas have condition numbers extremely close to one.

---

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please send email to trefethen@maths.ox.ac.uk.

  • Computational Mathematics and Applications Seminar
12 November 2020
14:00
Abstract

Whitney elements on simplices are perhaps the most widely used finite elements in computational electromagnetics. They offer the simplest construction of polynomial discrete differential forms on simplicial complexes. Their associated degrees of freedom (dofs) have a very clear physical meaning and give a recipe for discretizing physical balance laws, e.g., Maxwell’s equations. As interest grew for the use of high order schemes, such as hp-finite element or spectral element methods, higher-order extensions of Whitney forms have become an important computational tool, appreciated for their better convergence and accuracy properties. However, it has remained unclear what kind of cochains such elements should be associated with: Can the corresponding dofs be assigned to precise geometrical elements of the mesh, just as, for instance, a degree of freedom for the space of Whitney 1-forms belongs to a specific edge? We address this localization issue. Why is this an issue? The existing constructions of high order extensions of Whitney elements follow the traditional FEM path of using higher and higher “moments” to define the needed dofs. As a result, such high order finite k-elements in d dimensions include dofs associated to q-simplices, with k < q ≤ d, whose physical interpretation is obscure. The present paper offers an approach based on the so-called “small simplices”, a set of subsimplices obtained by homothetic contractions of the original mesh simplices, centered at mesh nodes (or more generally, when going up in degree, at points of the principal lattice of each original simplex). Degrees of freedom of the high-order Whitney k-forms are then associated with small simplices of dimension k only.  We provide an explicit  basis for these elements on simplices and we justify this approach from a geometric point of view (in the spirit of Hassler Whitney's approach, still successful 30 years after his death).   

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please send email to trefethen@maths.ox.ac.uk.

  • Computational Mathematics and Applications Seminar
Add to My Calendar