Past Forthcoming Seminars

15 March 2001
14:00
Abstract
The sphere is an important setting for applied mathematics, yet the underlying approximation theory and numerical analysis needed for serious applications (such as, for example, global weather models) is much less developed than, for example, for the cube. \\ \\ This lecture will apply recent developments in approximation theory on the sphere to two different problems in scientific computing. \\ \\ First, in geodesy there is often the need to evaluate integrals using data selected from the vast amount collected by orbiting satellites. Sometimes the need is for quadrature rules that integrate exactly all spherical polynomials up to a specified degree $n$ (or equivalently, that integrate exactly all spherical harmonies $Y_{\ell ,k}(\theta ,\phi)$ with $\ell \le n).$ We shall demonstrate (using results of M. Reimer, I. Sloan and R. Womersley in collaboration with W. Freeden) that excellent quadrature rules of this kind can be obtained from recent results on polynomial interpolation on the sphere, if the interpolation points (and thus the quadrature points) are chosen to be points of a so-called extremal fundamental system. \\ \\ The second application is to the scattering of sound by smooth three-dimensional objects, and to the inverse problem of finding the shape of a scattering object by observing the pattern of the scattered sound waves. For these problems a methodology has been developed, in joint work with I.G. Graham, M. Ganesh and R. Womersley, by applying recent results on constructive polynomial approximation on the sphere. (The scattering object is treated as a deformed sphere.)
  • Computational Mathematics and Applications Seminar
Prof Mark Stadtherr
Abstract
Continuing advances in computing technology provide the power not only to solve increasingly large and complex process modeling and optimization problems, but also to address issues concerning the reliability with which such problems can be solved. For example, in solving process optimization problems, a persistent issue concerning reliability is whether or not a global, as opposed to local, optimum has been achieved. In modeling problems, especially with the use of complex nonlinear models, the issue of whether a solution is unique is of concern, and if no solution is found numerically, of whether there actually exists a solution to the posed problem. This presentation focuses on an approach, based on interval mathematics, that is capable of dealing with these issues, and which can provide mathematical and computational guarantees of reliability. That is, the technique is guaranteed to find all solutions to nonlinear equation solving problems and to find the global optimum in nonlinear optimization problems. The methodology is demonstrated using several examples, drawn primarily from the modeling of phase behavior, the estimation of parameters in models, and the modeling, using lattice density-functional theory, of phase transitions in nanoporous materials.
  • Computational Mathematics and Applications Seminar
22 February 2001
14:00
Dr Oliver Ernst
Abstract
This talk reviews some recent joint work with Michael Eiermann and Olaf Schneider which introduced a framework for analyzing some popular techniques for accelerating restarted Krylov subspace methods for solving linear systems of equations. Such techniques attempt to compensate for the loss of information due to restarting methods like GMRES, the memory demands of which are usually too high for it to be applied to large problems in unmodified form. We summarize the basic strategies which have been proposed and present both theoretical and numerical comparisons.
  • Computational Mathematics and Applications Seminar
Dr Colin Campbell
Abstract
Support Vector Machines are a new and very promising approach to machine learning. They can be applied to a wide range of tasks such as classification, regression, novelty detection, density estimation, etc. The approach is motivated by statistical learning theory and the algorithms have performed well in practice on important applications such as handwritten character recognition (where they currently give state-of-the-art performance), bioinformatics and machine vision. The learning task typically involves optimisation theory (linear, quadratic and general nonlinear programming, depending on the algorithm used). In fact, the approach has stimulated new questions in optimisation theory, principally concerned with the issue of how to handle problems with a large numbers of variables. In the first part of the talk I will overview this subject, in the second part I will describe some of the speaker's contributions to this subject (principally, novelty detection, query learning and new algorithms) and in the third part I will outline future directions and new questions stimulated by this research.
  • Computational Mathematics and Applications Seminar
18 January 2001
14:00
Prof Francisco Marques
Abstract
The flow in a cylinder with a rotating endwall has continued to attract much attention since Vogel (1968) first observed the vortex breakdown of the central core vortex that forms. Recent experiments have observed a multiplicity of unsteady states that coexist over a range of the governing parameters. In spite of numerous numerical and experimental studies, there continues to be considerable controversy with fundamental aspects of this flow, particularly with regards to symmetry breaking. Also, it is not well understood where these oscillatory states originate from, how they are interrelated, nor how they are related to the steady, axisymmetric basic state. \\ \\ In the aspect ratio (height/radius) range 1.6 < $\Lambda$ < 2.8, the primary bifurcation is to an axisymmetric time-periodic flow (a limit cycle). We have developed a suite of numerical techniques, exploiting the biharmonic formulation of the problem in the axisymmetric case, that allows us to compute the nonlinear time evolution, the basic state, and its linear stability in a consistent and efficient manner. We show that the basic state undergoes a succession of Hopf bifurcations and the corresponding eigenvalues and eigenvectors of these excited modes describe most of the characteristics of the observed time-dependent states. \\ \\ The primary bifurcation is non-axisymmetric, to pure rotating wave, in the ranges $\Lambda$ <1.6 and $\Lambda$ > 2.8. An efficient and accurate numerical scheme is presented for the three-dimensional Navier-Stokes equations in primitive variables in a cylinder. Using these code, primary and secondary bifurcations breaking the SO(2) symmetry are analyzed. \\ \\ We have located a double Hopf bifurcation, where an axisymmetric limit cycle and a rotating wave bifurcate simultaneously. This codimension-2 bifurcation is very rich, allowing for several different scenarios. By a comprehensive two-parameter exploration about this point we have identified precisely to which scenario this case corresponds. The mode interaction generates an unstable two-torus modulate rotating wave solution and gives a wedge-shaped region in parameter space where the two periodic solutions are both stable. \\ \\ For aspect ratios around three, experimental observations suggest that the first mode of instability is a precession of the central vortex core, whereas recent linear stability analysis suggest a Hopf bifurcation to a rotating wave at lower rotation rates. This apparent discrepancy is resolved with the aid of the 3D Navier-Stokes solver. The primary bifurcation to an m=4 traveling wave, detected by the linear stability analysis, is located away from the axis, and a secondary bifurcation to a modulated rotating wave with dominant modes m=1 and 4, is seen mainly on the axis as a precessing vortex breakdown bubble. Experiments and the linear stability analysis detected different aspects of the same flow, that take place in different spatial locations.
  • Computational Mathematics and Applications Seminar

Pages