Past Computational Mathematics and Applications Seminar

14 June 2018
14:00
Prof. Joel Tropp
Abstract

Random matrices now play a role in many areas of theoretical, applied, and computational mathematics. Therefore, it is desirable to have tools for studying random matrices that are flexible, easy to use, and powerful. Over the last fifteen years, researchers have developed a remarkable family of results, called matrix concentration inequalities, that balance these criteria. This talk offers an invitation to the field of matrix concentration inequalities and their applications.

  • Computational Mathematics and Applications Seminar
7 June 2018
14:00
Prof. Max Gunzburger
Abstract

We first consider multilevel Monte Carlo and stochastic collocation methods for determining statistical information about an output of interest that depends on the solution of a PDE with inputs that depend on random parameters. In our context, these methods connect a hierarchy of spatial grids to the amount of sampling done for a given grid, resulting in dramatic acceleration in the convergence of approximations. We then consider multifidelity methods for the same purpose which feature a variety of models that have different fidelities. For example, we could have coarser grid discretizations, reduced-order models, simplified physics, surrogates such as interpolants, and, in principle, even experimental data. No assumptions are made about the fidelity of the models relative to the “truth” model of interest so that unlike multilevel methods, there is no a priori model hierarchy available. However, our approach can still greatly accelerate the convergence of approximations.

  • Computational Mathematics and Applications Seminar
24 May 2018
14:00
Prof. Michael Ferris
Abstract


In the past few decades, power grids across the world have become dependent on markets that aim to efficiently match supply with demand at all times via a variety of pricing and auction mechanisms. These markets are based on models that capture interactions between producers, transmission and consumers. Energy producers typically maximize profits by optimally allocating and scheduling resources over time. A dynamic equilibrium aims to determine prices and dispatches that can be transmitted over the electricity grid to satisfy evolving consumer requirements for energy at different locations and times. Computation allows large scale practical implementations of socially optimal models to be solved as part of the market operation, and regulations can be imposed that aim to ensure competitive behaviour of market participants.

Questions remain that will be outlined in this presentation.

Firstly, the recent explosion in the use of renewable supply such as wind, solar and hydro has led to increased volatility in this system. We demonstrate how risk can impose significant costs on the system that are not modeled in the context of socially optimal power system markets and highlight the use of contracts to reduce or recover these costs. We also outline how battery storage can be used as an effective hedging instrument.

Secondly, how do we guarantee continued operation in rarely occuring situations and when failures occur and how do we price this robustness?

Thirdly, how do we guarantee appropriate participant behaviour? Specifically, is it possible for participants to develop strategies that move the system to operating points that are not socially optimal?

Fourthly, how do we ensure enough transmission (and generator) capacity in the long term, and how do we recover the costs of this enhanced infrastructure?
 

  • Computational Mathematics and Applications Seminar
17 May 2018
14:00
Dr. Fehmi Cirak
Abstract

Advances in manufacturing technologies, most prominently in additive manufacturing or 3d printing, are making it possible to fabricate highly optimised products with increasing geometric and hierarchical complexity. This talk will introduce our ongoing work on design optimisation that combines CAD-compatible geometry representations, multiresolution geometry processing techniques and immersed finite elements with classical shape and topology calculus. As example applications,the shape optimisation of mechanical structures and electromechanical components, and the topology optimisation of lattice-skin structures will be discussed.

  • Computational Mathematics and Applications Seminar
10 May 2018
14:00
Abstract

The development of reduced order models for complex applications, offering the promise for rapid and accurate evaluation of the output of complex models under parameterized variation, remains a very active research area. Applications are found in problems which require many evaluations, sampled over a potentially large parameter space, such as in optimization, control, uncertainty quantification and applications where near real-time response is needed.

However, many challenges remain to secure the flexibility, robustness, and efficiency needed for general large-scale applications, in particular for nonlinear and/or time-dependent problems.

After giving a brief general introduction to reduced order models, we discuss developments in two different directions. In the first part, we discuss recent developments of reduced methods that conserve chosen invariants for nonlinear time-dependent problems. We pay particular attention to the development of reduced models for Hamiltonian problems and propose a greedy approach to build the basis. As we shall demonstrate, attention to the construction of the basis must be paid not only to ensure accuracy but also to ensure stability of the reduced model. Time permitting, we shall also briefly discuss how to extend the approach to include more general dissipative problems through the notion of port-Hamiltonians, resulting in reduced models that remain stable even in the limit of vanishing viscosity and also touch on extensions to Euler and Navier-Stokes equations.

The second part of the talk discusses the combination of reduced order modeling for nonlinear problems with the use of neural networks to overcome known problems of on-line efficiency for general nonlinear problems. We discuss the general idea in which training of the neural network becomes part of the offline part and demonstrate its potential through a number of examples, including for the incompressible Navier-Stokes equations with geometric variations.

This work has been done with in collaboration with B.F. Afkram (EPFL, CH), N. Ripamonti EPFL, CH) and S. Ubbiali (USI, CH).

  • Computational Mathematics and Applications Seminar
3 May 2018
14:00
Prof. Espen Jakobsen
Abstract


In this talk we will introduce and analyse a class of robust numerical methods for nonlocal possibly nonlinear diffusion and convection-diffusion equations. Diffusion and convection-diffusion models are popular in Physics, Chemistry, Engineering, and Economics, and in many models the diffusion is anomalous or nonlocal. This means that the underlying “particle" distributions are not Gaussian, but rather follow more general Levy distributions, distributions that need not have second moments and can satisfy (generalised) central limit theorems. We will focus on models with nonlinear possibly degenerate diffusions like fractional Porous Medium Equations, Fast Diffusion Equations, and Stefan (phase transition) Problems, with or without convection. The solutions of these problems can be very irregular and even possess shock discontinuities. The combination of nonlinear problems and irregular solutions makes these problems challenging to solve numerically.
The methods we will discuss are monotone finite difference quadrature methods that are robust in the sense that they “always” converge. By that we mean that under very weak assumptions, they converge to the correct generalised possibly discontinuous generalised solution. In some cases we can also obtain error estimates. The plan of the talk is: 1. to give a short introduction to the models, 2. explain the numerical methods, 3. give results and elements of the analysis for pure diffusion equations, and 4. give results and ideas of the analysis for convection-diffusion equations. 
 

  • Computational Mathematics and Applications Seminar
26 April 2018
14:00
Prof. Ron DeVore
Abstract


A very common problem in Science is that we have some Data Observations and we are interested in either approximating the function underlying the data or computing some quantity of interest about this function.  This talk will discuss what are best algorithms for such tasks and how we can evaluate the performance of any such algorithm.
 

  • Computational Mathematics and Applications Seminar
8 March 2018
14:00
Gabriel Barrenechea
Abstract

In this talk I will review recent results on the analysis of shock-capturing-type methods applied to convection-dominated problems. The method of choice is a variant of the Algebraic Flux-Correction (AFC) scheme. This scheme has received some attention over the last two decades due to its very satisfactory numerical performance. Despite this attention, until very recently there was no stability and convergence analysis for it. Thus, the purpose of the works reviewed in this talk was to bridge that gap. The first step towards the full analysis of the method is a rewriting of it as a nonlinear edge-based diffusion method. This writing makes it possible to present a unified analysis of the different variants of it. So, minimal assumptions on the components of the method are stated in such a way that the resulting scheme satisfies the Discrete Maximum Principle (DMP) and is convergence. One property that will be discussed in detail is the linearity preservation. This property has been linked to the good performance of methods of this kind. We will discuss in detail its role and the impact of it in the overall convergence of the method. Time permitting, some results on a posteriori error estimation will also be presented. 
This talk will gather contributions with A. Allendes (UTFSM, Chile), E. Burman (UCL, UK), V. John (WIAS, Berlin), F. Karakatsani (Chester, UK), P. Knobloch (Prague, Czech Republic), and 
R. Rankin (U. of Nottingham, China).

  • Computational Mathematics and Applications Seminar
1 March 2018
14:00
Prof Jan Hesthaven
Abstract

The development of reduced order models for complex applications, offering the promise for rapid and accurate evaluation of the output of complex models under parameterized variation, remains a very active research area. Applications are found in problems which require many evaluations, sampled over a potentially large parameter space, such as in optimization, control, uncertainty quantification and applications where near real-time response is needed.

However, many challenges remain to secure the flexibility, robustness, and efficiency needed for general large-scale applications, in particular for nonlinear and/or time-dependent problems.

After giving a brief general introduction to reduced order models, we discuss developments in two different directions. In the first part, we discuss recent developments of reduced methods that conserve chosen invariants for nonlinear time-dependent problems. We pay particular attention to the development of reduced models for Hamiltonian problems and propose a greedy approach to build the basis. As we shall demonstrate, attention to the construction of the basis must be paid not only to ensure accuracy but also to ensure stability of the reduced model. Time permitting, we shall also briefly discuss how to extend the approach to include more general dissipative problems through the notion of port-Hamiltonians, resulting in reduced models that remain stable even in the limit of vanishing viscosity and also touch on extensions to Euler and Navier-Stokes equations.

The second part of the talk discusses the combination of reduced order modeling for nonlinear problems with the use of neural networks to overcome known problems of on-line efficiency for general nonlinear problems. We discuss the general idea in which training of the neural network becomes part of the offline part and demonstrate its potential through a number of examples, including for the incompressible Navier-Stokes equations with geometric variations.

This work has been done with in collaboration with B.F. Afkram (EPFL, CH), N. Ripamonti EPFL, CH) and S. Ubbiali (USI, CH).

  • Computational Mathematics and Applications Seminar

Pages