Past Computational Mathematics and Applications Seminar

22 October 2020
14:00
Victoria Howle
Abstract

In this talk, we introduce a new preconditioner for the large, structured systems appearing in implicit Runge–Kutta time integration of parabolic partial differential equations. This preconditioner is based on a block LDU factorization with algebraic multigrid subsolves for scalability.

We compare our preconditioner in condition number and eigenvalue distribution, and through numerical experiments, with others in the literature. In experiments run with implicit Runge–Kutta stages up to s = 7, we find that the new preconditioner outperforms the others, with the improvement becoming more pronounced as the spatial discretization is refined and as temporal order is increased.

 

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
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
11 June 2020
14:00
Alex Townsend
Abstract

Consider a network of identical phase oscillators with sinusoidal coupling. How likely are the oscillators to globally synchronize, starting from random initial phases? One expects that dense networks have a strong tendency to synchronize and the basin of attraction for the synchronous state to be the whole phase space. But, how dense is dense enough? In this (hopefully) entertaining Zoom talk, we use techniques from numerical linear algebra and computational Algebraic geometry to derive the densest known networks that do not synchronize and the sparsest networks that do. This is joint work with Steven Strogatz and Mike Stillman.


[To be added to our seminars mailing list, or to receive a Zoom invitation for a particular seminar, please contact trefethen@maths.ox.ac.uk.]

  • Computational Mathematics and Applications Seminar
4 June 2020
14:00
Simon Chandler-Wilde
Abstract

The boundary integral equation method is a popular method for solving elliptic PDEs with constant coefficients, and systems of such PDEs, in bounded and unbounded domains. An attraction of the method is that it reduces solution of the PDE in the domain to solution of a boundary integral equation on the boundary of the domain, reducing the dimensionality of the problem. Second kind integral equations, featuring the double-layer potential operator, have a long history in analysis and numerical analysis. They provided, through C. Neumann, the first existence proof to the Laplace Dirichlet problem in 3D, have been an important analysis tool for PDEs through the 20th century, and are popular computationally because of their excellent conditioning and convergence properties for large classes of domains. A standard numerical method, in particular for boundary integral equations, is the Galerkin method, and the standard convergence analysis starts with a proof that the relevant operator is coercive, or a compact perturbation of a coercive operator, in the relevant function space. A long-standing open problem is whether this property holds for classical second kind boundary integral equations on general non-smooth domains. In this talk we give an overview of the various concepts and methods involved, reformulating the problem as a question about numerical ranges. We solve this open problem through counterexamples, presenting examples of 2D Lipschitz domains and 3D Lipschitz polyhedra for which coercivity does not hold. This is joint work with Prof Euan Spence, Bath.

 

[To be added to our seminars mailing list, or to receive a Zoom invitation for a particular seminar, please contact trefethen@maths.ox.ac.uk.]

  • Computational Mathematics and Applications Seminar
28 May 2020
14:00
Patrick Farrell
Abstract

We discuss two recent extensions of work on Reynolds-robust preconditioners for the Navier-Stokes equations, to non-Newtonian fluids and to the equations of magnetohydrodynamics.  We model non-Newtonian fluids by means of an implicit constitutive relation between stress and strain. This framework is broadly applicable and allows for proofs of convergence under quite general assumptions. Since the stress cannot in general be solved for in terms of the strain, a three-field stress-velocity-pressure formulation is adopted. By combining the augmented Lagrangian approach with a kernel-capturing space decomposition, we derive a preconditioner that is observed to be robust to variations in rheological parameters in both two and three dimensions.  In the case of magnetohydrodynamics, we consider the stationary incompressible resistive Newtonian equations, and solve a four-field formulation for the velocity, pressure, magnetic field and electric field. A structure-preserving discretisation is employed that enforces both div(u) = 0 and div(B) = 0 pointwise. The basic idea of the solver is to split the fluid and electromagnetic parts and to employ our existing Navier-Stokes solver in the Schur complement. We present results in two dimensions that exhibit robustness with respect to both the fluids and magnetic Reynolds numbers, and describe ongoing work to extend the solver to three dimensions.

[To be added to our seminars mailing list, or to receive a Zoom invitation for a particular seminar, please contact trefethen@maths.ox.ac.uk.]

  • Computational Mathematics and Applications Seminar
21 May 2020
14:00
Mark Embree
Abstract

In 2007, Andrew Mayo and Thanos Antoulas proposed a rational interpolation algorithm to solve a basic problem in control theory: given samples of the transfer function of a dynamical system, construct a linear time-invariant system that realizes these samples.  The resulting theory enables a wide range of data-driven modeling, and has seen diverse applications and extensions.  We will introduce these ideas from a numerical analyst's perspective, show how the selection of interpolation points can be guided by a Sylvester equation and pseudospectra of matrix pencils, and mention an application of these ideas to a contour algorithm for the nonlinear eigenvalue problem. (This talk involves collaborations with Michael Brennan (MIT), Serkan Gugercin (Virginia Tech), and Cosmin Ionita (MathWorks).)

[To be added to our seminars mailing list, or to receive a Zoom invitation for a particular seminar, please contact helen.mcgregor@maths.ox.ac.uk.]

  • Computational Mathematics and Applications Seminar

Pages