# Applied/ACMS/absS22: Difference between revisions

Line 13: | Line 13: | ||

Abstract: I’ll discuss a particular high-dimensional system that displays subtle behaviour found in the continuum limit. The only catch is that it formally shouldn’t, which raises a few questions. When is a discrete system large enough to be called continuous? When are approximate (broken) symmetries good enough to be treated like the real thing? When and why does a fluid approximation work as well as we like to assume? What does all this say about observables and the approach to equilibria? The particular system I have in mind is a large ideal pendulum chain, and it’s cousin the continuous flexible string. I propose that the pendulum chain is a perfect model system to study notoriously difficult phenomena such as vortical turbulence, waves, cascades and thermalisation, but with many fewer degrees of freedom than a three-dimensional fluid. | Abstract: I’ll discuss a particular high-dimensional system that displays subtle behaviour found in the continuum limit. The only catch is that it formally shouldn’t, which raises a few questions. When is a discrete system large enough to be called continuous? When are approximate (broken) symmetries good enough to be treated like the real thing? When and why does a fluid approximation work as well as we like to assume? What does all this say about observables and the approach to equilibria? The particular system I have in mind is a large ideal pendulum chain, and it’s cousin the continuous flexible string. I propose that the pendulum chain is a perfect model system to study notoriously difficult phenomena such as vortical turbulence, waves, cascades and thermalisation, but with many fewer degrees of freedom than a three-dimensional fluid. | ||

=== Xiangxiong Zhang (Purdue) === | |||

Title: Recent Progress on Q^k Spectral Element Method: Accuracy, Monotonicity and Applications | |||

Abstract: In the literature, spectral element methods usually refer to finite element methods with high order polynomial basis. The Q^k spectral element method has been a popular high order method for solving second order PDEs, e.g., wave equations, for more than three decades, obtained by continuous finite element method with tenor product polynomial of degree k and with at least (k+1)-point Gauss-Lobatto quadrature. In this talk, I will present some brand new results of this classical scheme, including its accuracy, monotonicity (stability), and examples of using monotonicity to construct high order bound-preserving schemes in various applications including the Allen-Cahn equation coupled with an incompressible velocity field, Keller-Segel equation for chemotaxis, and nonlinear eigenvalue problem for Gross–Pitaevskii equation. 1) Accuracy: when the least accurate (k+1)-point Gauss-Lobatto quadrature is used, the spectral element method is also a finite difference (FD) scheme, and this FD scheme can sometimes be (k+2)-th order accurate for k>=2. This has been observed in practice but never proven before in terms of rigorous error estimates. We are able to prove it for linear elliptic, wave, parabolic and Schrödinger equations for Dirichlet boundary conditions. For Neumann boundary conditions, (k+2)-th order can be proven if there is no mixed second order derivative. Otherwise, only (k+3/2)-th order can be proven and some order loss is indeed observed in numerical tests. The accuracy result also applies to spectral element method on any curvilinear mesh that can be smoothly mapped to a rectangular mesh, e.g., solving a wave equation on an annulus region with a curvilinear mesh generated by polar coordinates. 2) Monotonicity: consider solving the Poisson equation, then a scheme is called monotone if the inverse of the stiffness matrix is entrywise non-negative. It is well known that second order centered difference or P1 finite element method can form an M-matrix thus they are monotone, and high order accurate schemes in general are not M-matrices thus not monotone. But there are exceptions. In particular, we have proven that the fourth order accurate FD scheme (Q^2 spectral element method) is a product of two M-matrices thus monotone for a variable coefficient diffusion operator: this is the first time that a high order accurate scheme is proven monotone for a variable coefficient operator. We have also proven the fifth order accurate FD scheme (Q^3 spectral element method) is a product of three M-matrices thus monotone for the Poisson equation: this is the first time that a fifth order accurate discrete Laplacian is proven monotone in two dimensions (all previously known high order monotone discrete Laplacian in 2D are fourth order accurate). |

## Revision as of 18:53, 7 February 2022

# ACMS Abstracts: Spring 2022

### Alex Townsend (Cornell)

Title: What networks of oscillators spontaneously synchronize?

Abstract: Consider a network of identical phase oscillators with sinusoidal coupling. How likely are the oscillators to spontaneously synchronize, starting from random initial phases? One expects that dense networks of oscillators have a strong tendency to pulse in unison. But, how dense is dense enough? In this talk, we use techniques from numerical linear algebra, computational algebraic geometry, and dynamical systems to derive the densest known networks that do not synchronize and the sparsest ones that do. We will find that there is a critical network density above which spontaneous synchrony is guaranteed regardless of the network's topology, and prove that synchrony is omnipresent for random networks above a lucid threshold. This is joint work with Martin Kassabov, Steven Strogatz, and Mike Stillman.

Prof. Alex Townsend is an associate professor at Cornell University in the Mathematics Department. His research is in Applied Mathematics and mainly focuses on spectral methods, low-rank techniques, fast transforms, and theoretical aspects of deep learning. Prior to Cornell, he was an Applied Math instructor at MIT (2014-2016) and a DPhil student at the University of Oxford (2010-2014). He was awarded an NSF CAREER in 2021, a SIGEST paper award in 2019, the SIAG/LA Early Career Prize in applicable linear algebra in 2018, and the Leslie Fox Prize in numerical analysis in 2015.

### Geoffrey Vasil (Sydney)

Title: The mechanics of a large pendulum chain

Abstract: I’ll discuss a particular high-dimensional system that displays subtle behaviour found in the continuum limit. The only catch is that it formally shouldn’t, which raises a few questions. When is a discrete system large enough to be called continuous? When are approximate (broken) symmetries good enough to be treated like the real thing? When and why does a fluid approximation work as well as we like to assume? What does all this say about observables and the approach to equilibria? The particular system I have in mind is a large ideal pendulum chain, and it’s cousin the continuous flexible string. I propose that the pendulum chain is a perfect model system to study notoriously difficult phenomena such as vortical turbulence, waves, cascades and thermalisation, but with many fewer degrees of freedom than a three-dimensional fluid.

### Xiangxiong Zhang (Purdue)

Title: Recent Progress on Q^k Spectral Element Method: Accuracy, Monotonicity and Applications

Abstract: In the literature, spectral element methods usually refer to finite element methods with high order polynomial basis. The Q^k spectral element method has been a popular high order method for solving second order PDEs, e.g., wave equations, for more than three decades, obtained by continuous finite element method with tenor product polynomial of degree k and with at least (k+1)-point Gauss-Lobatto quadrature. In this talk, I will present some brand new results of this classical scheme, including its accuracy, monotonicity (stability), and examples of using monotonicity to construct high order bound-preserving schemes in various applications including the Allen-Cahn equation coupled with an incompressible velocity field, Keller-Segel equation for chemotaxis, and nonlinear eigenvalue problem for Gross–Pitaevskii equation. 1) Accuracy: when the least accurate (k+1)-point Gauss-Lobatto quadrature is used, the spectral element method is also a finite difference (FD) scheme, and this FD scheme can sometimes be (k+2)-th order accurate for k>=2. This has been observed in practice but never proven before in terms of rigorous error estimates. We are able to prove it for linear elliptic, wave, parabolic and Schrödinger equations for Dirichlet boundary conditions. For Neumann boundary conditions, (k+2)-th order can be proven if there is no mixed second order derivative. Otherwise, only (k+3/2)-th order can be proven and some order loss is indeed observed in numerical tests. The accuracy result also applies to spectral element method on any curvilinear mesh that can be smoothly mapped to a rectangular mesh, e.g., solving a wave equation on an annulus region with a curvilinear mesh generated by polar coordinates. 2) Monotonicity: consider solving the Poisson equation, then a scheme is called monotone if the inverse of the stiffness matrix is entrywise non-negative. It is well known that second order centered difference or P1 finite element method can form an M-matrix thus they are monotone, and high order accurate schemes in general are not M-matrices thus not monotone. But there are exceptions. In particular, we have proven that the fourth order accurate FD scheme (Q^2 spectral element method) is a product of two M-matrices thus monotone for a variable coefficient diffusion operator: this is the first time that a high order accurate scheme is proven monotone for a variable coefficient operator. We have also proven the fifth order accurate FD scheme (Q^3 spectral element method) is a product of three M-matrices thus monotone for the Poisson equation: this is the first time that a fifth order accurate discrete Laplacian is proven monotone in two dimensions (all previously known high order monotone discrete Laplacian in 2D are fourth order accurate).