Videos and presentation materials from other INI events are also available.
Event  When  Speaker  Title  Presentation Material 

HOP 
17th January 2007 14:30 to 15:30 
An overview of oscillatory integrals and integral operators in high frequency scattering In this talk we review recent research on integral equation methods for high frequency timeharmonic scattering. A number of interesting questions arise in this context, including:
We will indicate the (exciting but limited) progress that has been made internationally in addressing these questions in the last few years, and give some idea of the methods that have been used, focussing particularly on the simplest problem of scattering of timeharmonic acoustic waves by a bounded sound soft obstacle in two and three dimensions. We will also list some of the many open, difficult problems in this area. 

HOP 
24th January 2007 14:30 to 15:30 
Circlets, vectorands and real linear operators  
HOP 
30th January 2007 14:30 to 15:30 
What can happen to the spectrum in lowrank perturbations?  
HOP 
30th January 2007 16:00 to 17:00 
K Sankaran  Finitevolume timedomain modelling of the Maxwell system and domaintruncation techniques  
HOP 
31st January 2007 14:30 to 15:30 
Numerical simulation of BoseEinstein condensation  
HOP 
1st February 2007 16:00 to 17:00 
S Kamvissis  On nonlinear steepest descent  
HOP 
7th February 2007 14:30 to 15:30 
A Abdulle  Macrotomicro methods for highly oscillatory partial differential equations  
HOP 
13th February 2007 14:30 to 15:30 
Filontype methods for highly oscillatory integrals  
HOP 
13th February 2007 16:00 to 17:00 
A Iserles  Computing modified Fourier approximations in cubes  
HOP 
14th February 2007 11:30 to 12:30 
D Levin  Infinite oscillatory integrals  
HOP 
14th February 2007 14:30 to 15:30 
Computational high frequency wave propagation The numerical approximation of high frequency wave propagation is important in many applications, including geophysics, electromagnetics and acoustics. When the wavelength is short compared to the overall size of the computational domain direct simulation using standard wave equations is very expensive. Fortunately, there are computationally much less costly models, that are good approximations of many wave equations precisely for very high frequencies. In this talk I will focus on the geometrical optics approximation, which is the infinite frequency limit of wave equations. Geometrical optics has traditionally been simulated using ray tracing. More recently numerical procedures based on various partial differential equations have been introduced. In this talk I will survey such numerical methods for geometrical optics and also briefly discuss finite frequency corrections and some other models 

HOP 
15th February 2007 11:30 to 12:30 
Momentfree numerical integration of highly oscillatory functions  
HOP 
15th February 2007 14:30 to 15:30 
A numerical approach to the method of steepest descent  
HOP 
21st February 2007 11:30 to 12:30 
Some criteria for convergence of Fourier series in different function spaces  
HOP 
21st February 2007 14:30 to 15:30 
Regularized combined field integral equations  
HOP 
26th February 2007 14:30 to 15:30 
O Runborg  A fast phase space method for computing creeping waves  
HOP 
27th February 2007 16:00 to 17:00 
S Fulton  Highly accurate computation of the DirichletNeumann map  
HOP 
28th February 2007 11:30 to 12:30 
Exponential operator splitting for timedependent Schroedinger equation  
HOP 
28th February 2007 14:30 to 15:30 
P Monk  The use of plane waves to approximate wave problems  
HOP 
1st March 2007 14:30 to 15:30 
Algorithms for sparse approximation and applications  
HOP 
2nd March 2007 14:30 to 15:30 
V Smyshlyaev  High frequency asymptotics in scattering by smooth and nonsmooth obstacles: analytical and numerical aspects  
HOP 
6th March 2007 10:00 to 11:00 
Robust boundary integral methods in high frequency acoustic scattering  
HOP 
7th March 2007 11:30 to 12:30 
Analysis of plane wave discontinuous Galerkin methods  
HOP 
7th March 2007 14:30 to 15:30 
Integrating exponentials over a simplex: a surprising link with exponential Brownian motion  
HOP 
12th March 2007 14:30 to 15:30 
High frequency scattering by a collection of convex bodies  
HOP 
13th March 2007 11:30 to 12:30 
Nonlinear diffusions as limits of BGKtype kinetic equations  
HOP 
14th March 2007 11:30 to 12:30 
R Melnik  Cayley Transform Techniques in PDEs and their Numerical Approximations (on the example of convectiondiffusion equations)  
HOP 
14th March 2007 14:30 to 15:30 
A fast multipole method for oscillatory kernels  
HOP 
15th March 2007 16:00 to 17:00 
On asymptotic regimes for the MaxwellDirac system  
HOP 
20th March 2007 14:30 to 15:30 
Backward error analysis and its limitations  
HOP 
20th March 2007 16:00 to 17:00 
Molecular dynamics, thermostats and the accuracy of averages  
HOP 
21st March 2007 14:30 to 15:30 
Longtimestep methods for highly oscillatory Hamiltonian problems  
HOP 
21st March 2007 16:00 to 17:00 
S Reich  Fluids, parcel dynamics and oscillations  
HOP 
22nd March 2007 16:00 to 17:00 
Analysis of splitting methods for reactiondiffusion equations in the light of stochastic calculus  
HOPW02 
26th March 2007 10:00 to 11:00 
Three recent results on asymptotics of oscillations The results are separate, and apparently paradoxical, and have implications for physics. First, when two exponentials compete, their interference can be dominated by the contribution with smaller exponent. Second, repeated differentiation of almost all functions in a wide class generates trigonometric oscillations (almost all functions tend to cosx). Third, it is possible to find bandlimited functions that oscillate arbitrarily faster than their fastest Fourier component (superoscillations). 

HOPW02 
26th March 2007 11:30 to 12:15 
Heterogeneous multiscale methods Continuum simulations of solids or fluids for which some atomistic information is needed are typical example of multiscale problems with very large ranges of scales. For such problems it is necessary to restrict the simulations on the microscale to a smaller subset of the full computational domain. The heterogeneous multiscale method is a framework for developing and analyzing numerical methods that couple computations from very different scales. Local microscale simulations on small domains supply missing data to a macroscale simulation on the full domain. Examples are local molecular dynamics computations that produce data to a continuum macroscale model, or a highly oscillatory dynamical system for which a local estimate of resonances is enough to supply data for a smoother evolution of averages. We will focus on dynamical systems. 

HOPW02 
26th March 2007 14:15 to 15:00 
Hamiltonian slow manifolds with internal oscillation Many Hamiltonian systems possess special families of solutions which can be described approximately as slowly drifting periodic orbits. Examples include the gravitational threebody problem, the interaction of two identical charged particles in a magnetic field, and the propagation and interaction of discrete breathers (timeperiodic spatially localised motions in Hamiltonian networks of units). Theory will be presented for how to improve a zeroth order manifold of approximate solutions to rth order for any r > 0, meaning one that contains all nearby periodic orbits of nearby period and has error of order the rth power of the drift field (even with a small constant for r=1). In the normally hyperbolic case an exactly invariant nearby submanifold can be constructed. If there are normally elliptic directions, however, this is impossible in general but the above rth order approximations can be achieved provided that the normal frequencies avoid all multiples of that of the approximately periodic motion and are fast compared with the drift. An effective Hamitlonian is derived to describe the drift of the orbits. Applications to the above fields will be given. An introduction has been published as R.S.MacKay, Slow Manifolds, in Energy localisation and transfer, eds T.Dauxois, A.LitvakHinenzon, R.S.MacKay, A.Spanoudaki (World Sci Publ Co, 2004), 149192. 

HOPW02 
26th March 2007 15:30 to 16:15 
Dispersive and dissipative behaviour of Galerkin approximation using high order polynomials We consider the dispersive properties of Galerkin finite element methods for wave propagation. The dispersive properties of conforming finite element scheme are analysed in the setting of the Helmholtz equation and an explicit form the discrete dispersion relation is obtained for elements of arbitrary order. It is shown that the numerical dispersion displays three different types of behaviour depending on the order of the polynomials used relative to the meshsize and the wave number. Quantitative estimates are obtained for the behaviour and rates of decay of the dispersion error in the differing regimes. 

HOPW02 
26th March 2007 16:15 to 17:00 
S Dobrokhotov 
Generalization of the Maslov Theory for Localized Asymptotics and Tsunami waves We suggest a new asymptotic representation for the solutions to the multidimensional wave equations with variable velocity with localized initial data. This representation is the generalization of the Maslov canonical operator based also on a simple relationship between fast decaying and fast oscillating solutions, and on boundary layer ideas. It establishes the connection between initial localized perturbations and wave profiles near the wave fronts including the neighborhood of backtracking (focal or turning) and self intersection points. We show that wave profiles are related with a form of initial sources and also with the Lagrangian manifolds organized by the rays and wavefronts. In particular we discuss the influence of such topological characteristics like the Maslov and Morse indices to metamorphosis of the profiles after crossing the focal points. We apply these formulas to the problem of a propagation of tsunami waves in the frame of socalled ``piston model''. Finally we suggest a fast asymptoticallynumerical algorithm for simulation of tsunami wave over nonuniform bottom. Different scenarios of the distribution of the waves are considered, the wave profiles of the front are obtained in connection with the different shapes of the source and with the diverse rays generating the fronts. It is possible to use suggested algorithm to predict in real time the zones of the beaches where the amplitude of the tsunami wave has dangerous high values. In this connection we also discuss the following questions: the problems of the regularization of the wave field near focal points; illpossed problems appearing in the geometry of the wavefronts; the inverse problem connected with the possibility of reconstruction of the source via the measurement of the tsunami wave profile on the shelf etc. This work was done together with S.SekerzhZenkovich, B.Tirozzi, B.Volkov and was partially supported by RFBR grant N 050100968 and Agreement Between University "La Sapienza", Rome and Institute for Problems in Mechanics RAS, Moscow. Bibliography [1] S.Yu. Dobrokhotov, S.Ya SekerzhZenkovich, B. Tirozzi, T.Ya. Tudorovskiy, The description of tsunami waves propagation based on the Maslov canonical operator, Doklady Mathematics, 2006, v.74, N 1, pp. 592596 [2] S.Yu. Dobrokhotov, S.Ya SekerzhZenkovich, B. Tirozzi, T.Ya. Tudorovskiy, Asymptotic theory of tsunami waves: geometrical aspects and the generalized Maslov representation, Publications of Kyoto Research Mathematical Institute, Vol.4, page 118153, 2006, ISSN 18802818. [3] S.Dobrokhotov, S.SekerzhZenkovich, B.Tirozzi, B.Volkov Explicit asymptotics for tsunami waves in framework of the piston model, Russ. Journ. Earth Sciences, 2006, v.8, ES403, pp.112 [4] S.Dobrokhotov, S.Sinitsyn, B.Tirozzi, Asymptotics of Localized Solutions of the OneDimensional Wave Equation with Variable Velocity. I. The Cauchy Problem, Russ. Jour.Math.Phys., v.14, N1, 2007, pp.2856 

HOPW02 
26th March 2007 17:00 to 17:30 
Numerical study of oscillatory regimes in the KadomtsevPetviashvili equation The aim of this talk is the accurate numerical study of the KP equation. In particular we are concerned with the small dispersion limit of this model, where no comprehensive analytical description exists so far. To this end we first study a similar highly oscillatory regime for asymptotically small solutions, which can be described via the DaveyStewartson system. In a second step we investigate numerically the small dispersion limit of the KP model in the case of large amplitudes. Similarities and differences to the much better studied Kortewegde Vries situation are discussed as well as the dependence of the limit on the additional transverse coordinate. Related Links 

HOPW02 
27th March 2007 09:00 to 09:45 
Numerical oscillations over long times  
HOPW02 
27th March 2007 09:45 to 10:30 
Conservation of energy and actions in numerical discretizations of nonlinear wave equations For numerical discretizations of nonlinearly perturbed wave equations the longtime nearconservation of energy, momentum, and harmonic actions is studied. The proofs are based on the technique of modulated Fourier expansions in time. Rigorous statements on the longtime conservation properties are shown under suitable numerical nonresonance conditions and under a CFL condition. The time step need not be small compared to the inverse of the largest frequency in the spacediscretized system. This is joint work with Christian Lubich and David Cohen. 

HOPW02 
27th March 2007 10:30 to 11:00 
Highly oscillatory Hamiltonian systems with nonconstant mass matrix We will present a class of numerical methods (based on the trigonometric methods) for such Hamiltonian problems. We will then present a frequency expansion of the numerical solution: the modulated Fourier expansion. The system that determines the coefficients of this expansion has two formal invariants which are related to the total energy and the oscillatory energy of the original system. This allows us to prove the nearconservation of the total and the oscillatory energy for the numerical schemes over very long time intervals. Related Links


HOPW02 
27th March 2007 11:30 to 12:15 
S Reich 
Highly oscillatory PDEs, slow manifolds and regularized PDE formations The main motivation of my talk is provided by geophysical fluid dynamics. The underlying Euler or NavierStokes equations display oscillatory wave dynamics on a wide range of temporal and spatial scales. Simplified models are often based on the idea of balance and the concept of a slow manifold. Examples are provided by hydrostatic and geostrophic balance. One would also like to exploit these concepts on a computational level. However, slow manifolds are idealized objects that do not fully characterize the complex fluid behavior. I will describe a novel regularization technique that makes use of balance and slow manifolds in an adaptive manner. The regularization approach is based on a reinterpretation of the (linearly) implicit midpoint rule as an explicit timestepping method applied to a regularized set of Euler equations. Adaptivity can be achieved by means of a predictorcorrector interpretation of the regularization. 

HOPW02 
27th March 2007 14:15 to 15:00 
Tunnel effect for KramersFokkerPlanck type operators For a large class of KramersFokkerPlanck type operators, we determine in the semiclassical (here the low temperature) limit the full asymptotic expansion of the splitting between the lowest eigenvalue (0) and the next one. In a previous work with Herau and C. Stolk we did so in the case when a certain potential (or exponent of a Maxwellian) has preciely one local minimum and then the splitting is "large". In this new work we treat the case when the potential has two local minima. Then the splitting is exponentially small and related to a tunnel effect between the minima. Our most direct source of inspiration has been works by HerauF.Nier, B.HelfferNier, but our methods are quite different. Joint work with F Herau and M Hitrik. 

HOPW02 
27th March 2007 15:30 to 16:15 
Fast singular oscillating limits of the 3D NavierStokes and Euler equations: Global regularity and threedimensional Euler dynamics In the first part of the talk, we study highly oscillatory problems for the incompressible 3D NavierStokes and prove existence on infinite time intervals of regular strong solutions; smoothness assumptions for initial data are the same as in local existence theorems. There are no conditional assumptions on the properties of solutions at later times, nor are the global solutions close to any 2D manifold. The global existence is proven using techniques of fast singular oscillating limits, Lemmas on restricted convolutions and the LittlewoodPaley dyadic decomposition. In the second part of the talk, we analyze regularity and dynamics of highly oscillatory problems for the 3D Euler equations. Detailed proofs can be found in the following references: Asymptotic Analysis, vol. 15, No. 2, p. 103150, 1997 (with A. Babin and B. Nicolaenko). Indiana University Mathematics Journal, vol. 48, No. 3, p. 11331176, 1999 (with A. Babin and B. Nicolaenko). Indiana University Mathematics Journal, vol. 50, p. 135, 2001 (with A. Babin and B. Nicolaenko). Russ. Mathematical Surveys, vol. 58, No. 2 (350), p. 287318, 2003 (with B. Nicolaenko). Methods and Applications of Analysis, vol. 11, No. 4, p. 605633, 2004 (with B. Nicolaenko, C. Bardos and F. Golse). Hokkaido Mathematical Journal, vol. 35, No. 2, p. 321364, 2006 (with Y. Giga, K. Inui and S. Matsui). Annals of Math. Studies, Princeton University Press (Editors: J. Bourgain and C. Kenig), In Press (with Y. Giga and B. Nicolaenko). 

HOPW02 
27th March 2007 16:15 to 17:00 
Some remarks concerning exponential integrators TBA 

HOPW02 
27th March 2007 17:00 to 17:30 
Spectral results for membrane with perturbed stiffness and mass density We study the spectrum of a nonhomogeneous membrane that consists of two parts with strongly different stiffness and mass density. The small parameter describes the quotient of stiffness coefficients. The Mth power of the parameter is comparable to the ratio of densities. We show that the asymptotic behaviour of eigenvalues and eigenfunctions depends on rate M. We distinguish three cases M<1, M=1 and M>1. The strong resolvent convergence of perturbed operators leads to loss of the completeness of limit eigenfunction system in both cases when M is different from 1. Therefore the limit operators describe only a part of the prelimit membrane vibrations. With an eye to close this gap we use the WKBasymptotic expansions with a quantized small parameter to prove the existence of other kind of eigenvibrations, namely high frequency vibrations. In the critical case M=1 the limit operator is a nonselfadjoint one, nevertheless the perturbed operators are selfadjoint in a certain topology. The multiplicity of spectrum and structure of root spaces are investigated. Complete asymptotic expansions for eigenelements are constructed and justified in each case of the perturbations. 

HOPW02 
28th March 2007 09:00 to 09:45 
Nanomechanics of biomolecules: whether oscillatory details are important or not  
HOPW02 
28th March 2007 09:45 to 10:30 
M Hochbruck 
Exponential integrators for oscillatory secondorder differential equations In this talk, we analyse a family of exponential integrators for secondorder differential equations in which highfrequency oscillations in the solution are generated by a linear part. We characterise methods which allow secondorder error bounds by presenting a unified error analysis for the whole family of methods. A major advantage of our analysis is that it does not require bounds for pointwise products of matrices and therefore, generalises to abstract differential equations, where the linear part is an unbounded operator with infinitely many large eigenvalues directly. This is joint work with Volker Grimm 

HOPW02 
28th March 2007 10:30 to 11:00 
Small dispersion limit of the KdV and semiclassical limit of focusing NLS We compare numerically the small dispersion limit of the Korteweg de Vries equation with the asymptotic formulas obtained in the works of LaxLevermore and DeifVenakides and Zhou. We compare numerically the semiclassical limit of the focusing nonlinear Schrodinger equation with the asymptotic formula obtained in the works of KamvissisMcLaughlinMiller and TovbisVenakidesZhou. As a results we outline the regions of the (x,t) where the above asymptotic formulas give a satisfactory description of the respective equations. 

HOPW02 
28th March 2007 11:30 to 12:15 
Beyond the adiabatic approximation: exponentially small coupling terms For multilevel timedependent quantum systems one can construct superadiabatic representations in which the coupling between separated levels is exponentially small in the adiabatic limit. We explicitly determine the asymptotic behavior of the exponentially small coupling term for generic twostate systems with realsymmetric Hamiltonian. The superadiabatic coupling term takes a universal form and depends only on the location and the strength of the complex singularities of the adiabatic coupling function. First order perturbation theory for the Hamiltonian in the superadiabatic representation then allows to describe the timedevelopment of exponentially small adiabatic transitions and thus to rigorously confirm Michael Berrys predictions on the universal form of adiabatic transition histories (Joint work with V. Betz from Warwick). Related Links


HOPW02 
28th March 2007 14:00 to 14:30 
Approximate energy conservation for symplectic time semidiscretizations of semilinear wave equations We prove that Astable symplectic RungeKutta time semidiscretizations applied to semilinear wave equations with periodic boundary conditions, analytic nonlinearity and analytic initial data conserve a modified energy p to an exponentially small error. This modified energy is O(h^p)close to the original energy where p is the order of the method and h the timestepsize. Standard backward error analysis can not be applied because of the occurrence of unbounded operators in the construction of the modified vectorfield. This loss of regularity in the construction can be taken care of by projecting the PDE to a finitedimensional space and by coupling the number of excited modes as well as the number of terms in the expansion of the vectorfield with the stepsize. This way we obtain exponential estimates of the form O(\exp(C/ h^(1/2) )). A similar technique has been used for averaging of rapidly forced Hamiltonian PDEs by [Matthies and Scheel, 2003]. As a sideproduct, we also provide a convergence analysis of RungeKutta methods in Hilbert spaces. Related Links


HOPW02 
28th March 2007 14:30 to 15:00 
Numerical approximation of multivariate highly oscillatory integrals The aim of this talk is to describe several methods for numerically approximating the integral of a multivariate highly oscillatory function. We begin with a review of the asymptotic and Filontype methods developed by Iserles and Nørsett. Using a method developed by Levin as a point of departure we will construct a new method that uses the same information as a Filontype method, and obtains the same asymptotic order, while not requiring moments. This allows us to integrate over nonsimplicial domains, and with complicated oscillators. We also present a method for approximating oscillatory integrals with stationary points. 

HOPW02 
28th March 2007 15:30 to 16:15 
Asymptotics for the CamassaHolm equation I will present recent results on asymptotic behaviors for the Camassa–Holm (CH) equation ut  utxx + 2?ux + 3uux = 2uxuxx + uuxxx on the line, ? being a nonnegative parameter. Firstly, I will describe the longtime asymptotic behavior of the solution u? (x, t), ? > 0 of the initialvalue problem with fast decaying initial data u0(x). It appears that u? (x, t) behaves differently in different sectors of the (x, t)halfplane. Then I will analyse the behavior of u? (x, t) as ? 0. The methods are inverse scattering in a matrix RiemannHilbert approach and Deift and Zhou’s nonlinear steepest descent method. Work in collaboration with Dmitry Shepelsky. 

HOPW02 
28th March 2007 16:15 to 17:00 
Norms and condition numbers of oscillatory integral operators in acoustic scattering In this talk we discuss domain and boundary integral operators arising in the theory and numerical treatment by integral equation methods of the Helmholtz equation or time harmonic Maxwell equations. These integral operators are increasingly oscillatory as the wave number k increases (k proportional to the frequency of the time harmonic incident field). An interesting theoretical question, also of practical significance, is the dependence of the norms of these integral operators and their inverses on k. We investigate this question, in particular for classical single and doublelayer potential operators for the Helmholtz equation on the boundary of bounded Lipschitz domains. The results and techniques used depend on the domain geometry. In certain 2D cases (for example where the boundary is a starlike polygon) bounds which are sharp in their dependence on k can be obtained, but there are many open problems for more general geometries and higher dimension. Related Links


HOPW02 
28th March 2007 17:00 to 17:30 
Exponential asymptotics for the primitive equations Following the work of Matthies (2001), it is shown how Gevrey (exponential) regularity of the solution and a classical method used together to prove an exponentiallyaccurate approximation result for a singular perturbation problem with a small parameter. The model considered is the viscous primitive equations of the ocean, although the method is applicable more generally. 

HOPW02 
29th March 2007 09:00 to 09:45 
Finite difference approximation of homogenization problems for elliptic equations In this talk, the problem of the approximation by finite differences of solutions to elliptic problems with rapidly oscillating coefficients and periodic boundary conditions will be discussed. The meshsize is denoted by $h$ while $\ee$ denotes the period of the rapidly oscillating coefficient. Using Bloch wave decompositions, we analyze the case where the ratio $h/\ee$ is rational. We show that if $h/\ee$ is kept fixed, being a rational number, even when $h,\ee\to 0$, the limit of the numerical solution does not coincide with the homogenized one obtained when passing to the limit as $\ee\to 0$ in the continuous problem. Explicit error estimates are given showing that, as the ratio $h/\ee$ approximates an irrational number, solutions of the finite difference approximation converge to the solutions of the homogenized elliptic equation. We consider both the 1d and the multidimensional case. Our analysis yields a quantitative version of previous results on numerical homogenization by Avellaneda, Hou and Papanicolaou, among others. This is a joint work with Rafael Orive 

HOPW02 
29th March 2007 09:45 to 10:30 
S Venakides 
Rigorous leading order asymptotic solutions of the semiclassical focusing NLS equation Solving an integrable nonlinear differential equation reduces to solving the matrix RiemannHilbert problem (RHP) that arises from the invese scattering procedure. Typically, this means finding a matrix m(z) that is analytic in the closed complex plane except on a given oriented contour on which m(z) has a jump discontinuity: the value of m(z) left of the contour equals its value on the right multiplied by a given "jump matrix'' V(z). The spacetime variables appear as parameters in the jump matrix; the complex variable z is the spectral variable of the linear differential operator that effects the linearization of the nonlinear system through the Lax pair associated with the system. We derive an explicit formula for the leading asymptotic (oscillatory) solution of the semiclassical focusing nonlinear Schroedinger equation for a class of initial data and prove an error estimate. We extend the analysis when time is large, the asymptotic formula is then expressed in terms of elementary functions. We first outline the notion of a Lax pair and how the RHP arises naturally.We then describe the steepest descent procedure that leads to an explicitly solvable asymptotic RHP through jump matrix factorizations and contour deformations. We also outline the gfunciton mechanism and describe how it determines the successful deformations. 

HOPW02 
29th March 2007 10:30 to 11:00 
SJA Malham 
Numerically evaluating the Evans function by Magnus integration We use Magnus methods to compute the Evans function for spectral problems as arise when determining the linear stability of travelling wave solutions to reactiondiffusion and related partial differential equations. In a typical application scenario, we need to repeatedly sample the solution to a system of linear nonautonomous ordinary differential equations for different values of one or more parameters as we detect and locate the zeros of the Evans function in the right half of the complex plane. In this situation, a substantial portion of the computational effortthe numerical evaluation of the iterated integrals which appear in the Magnus seriescan be performed independent of the parameters and hence needs to be done only once. More importantly, for any given tolerance Magnus integrators possess lower bounds on the step size which are uniform across large regions of parameter space and which can be estimated a priori. We demonstrate, analytically as well as through numerical experiment, that these features render Magnus integrators extremely robust and, depending on the regime of interest, efficient in comparison with standard ODE solvers. 

HOPW02 
29th March 2007 11:30 to 12:15 
Y Efendiev 
Multiscale numerical methods using limited global information In this talk, I will describe multiscale numerical methods for flows in heterogeneous porous media. The main idea of these methods is to construct local basis functions that can capture the small scale information when they are coupled via some global formulation. I will discuss the use of local boundary conditions, oversampling methods and the use of global information in constructing basis functions. Applications of these methods to stochastic equations will be also discussed. This is joint work with J. Aarnes, T. Hou, L. Jiang, V. Ginting. 

HOPW02 
29th March 2007 14:00 to 14:30 
Exponential integrators and functions of the matrix exponential Exponential integrators are the most efficient class of methods for the timestepping of stiff, semilinear, oscillatory PDEs such as the KdV equation. They solve the stiff, linear part of the PDE exactly. In the case of periodic boundary conditions, a Fourier spectral method can be used, so the linear part is diagonal and the methods can be applied straightforwardly. For other spatial discretizations, functions of the matrix exponential are required, which are susceptible to rounding errors. Several methods for evaluating these functions will be discussed. 

HOPW02 
29th March 2007 14:30 to 15:00 
VH Hoang 
High dimensional finite element methods for elliptic problems with highly oscillating coefficients Elliptic homogenization problems in a domain $\Omega \subset \R^d$ with $n+1$ separated scales are reduced to elliptic onescale problems in dimension $(n+1)d$. These onescale problems are discretized by a sparse tensor product finite element method (FEM). We prove that this sparse FEM has accuracy, work and memory requirement comparable to standard FEM for single scale problems in $\Omega$ while it gives numerical approximations of the correct homogenized limit as well as of all first order correctors, throughout the physical domain with performance independent of the physical problem's scale parameters. Numerical examples for model diffusion problems with two and three scales confirm our results. 

HOPW02 
29th March 2007 15:30 to 16:15 
A multiscale method for stiff ordinary differential equations with resonance We introduce a multiscale method to compute the effective behavior of a class of stiff and highly oscillatory ODEs. The oscillations may be in resonance with one another and thereby generate some hidden slow dynamics. Our method relies on correctly tracking a set of slow variables whose dynamics is effectively closed, and is sufficient to approximate the effective behavior of the ODEs. This set of variables is found by our numerical methods. We demonstrate our algorithms by a few examples that include a commonly studied problem of Fermi, Pasta, and Ulam (FPU). 

HOPW02 
29th March 2007 16:15 to 17:00 
VP Smyshlyaev 
Localization and propagation in high contrast highly oscillatory media via 'nonclassical' homogenization We discuss fundamental mathematical tools for analysis of localization and propagation effects in highly oscillatory media with high contrasts. With the underlying two small parameters of the oscillations and the contrast, there is a "critical" scaling when the phenomena at the micro and macro scales are coupled in a nontrivial way, with "unusual" effects observed in an asymptotically explicit way. The related mathematical tool is that of a "nonclassical" (highcontrast) homogenization, accounting for "highfrequency" oscillations, as opposes to the "classical" homogenization whose scope is limited by dealing in effect with low frequencies only. Those tools include "nonclassical" twoscale asymptotic expansions, twoscale operator and spectral convergence, and twoscale compactness (with the latter building on, among others, recent deep ideas of V.V. Zhikov [1]). We illustrate this on the problem of wave localization in high contrast periodic media with a defect [2] (a problem relevant to photonic crystal fibres). We also discuss the use of these techniques to the problem of "slowing down" of wave packets in high contrast highlyoscillatory media (the socalled "slow light" effect, with relevance to coupled resonances and metastability), and other prospects. Joint work with Ilia V. Kamotski. [1] V.V. Zhikov, On an extension of the method of twoscale convergence and its applications, (Russian) Mat. Sbornik 191 (2000), 3172; English translation in Sbornik Math. 191 (2000), 9731014; V.V. Zhikov, Gaps in the spectrum of some elliptic operators in divergent form with periodic coefficients. (Russian) Algebra i Analiz 16 (2004), 3458; [2] I.V. Kamotski and V.P. Smyshlyaev, Localised eigenstates due to defects in high contrast periodic media via homogenisation. BICS preprint 3/06. http://www.bath.ac.uk/mathsci/bics/preprints/BICS06_3.pdf (2006) Related Links


HOPW02 
29th March 2007 17:00 to 17:30 
Exponential estimates in averaging and homogenisation Many partial differential equations with rapid spatial or temporal scales have effective descriptions which can be derived by homogenisation or averaging. In this talk we deal with examples, where quantitative estimates of the error is possible for higher order homogenisation and averaging. In particular, we provide theorems, which allow homogenisation and averaging beyond all orders by giving exponential estimates of appropriately averaged and homogenised descriptions. Methods include iterated averaging transformations, optimal truncation of asymptotic expansions and highly regular solutions (Gevrey regularity). Prototypical examples are reactiondiffusion equations with heterogeneous reaction terms or rapid external forcing, nonlinear Schr\"odinger equations describing dispersion management, and secondorder linear elliptic equations. Related Links


HOPW02 
30th March 2007 09:00 to 09:45 
A Iserles 
From high oscillation to rapid approximation In this talk we address our recent work on the theory of modified Fourier expansions and their generalisation for the approximation of smooth functions in one and several variables. In particular, we debate rapid evaluation of expansion coefficients by means of Filonlike and GaussTuranlike quadrature techniques and various acceleration techniques. 

HOPW02 
30th March 2007 09:45 to 10:30 
On the zeromass constraint for KP type equations For a rather general class of equations of KadomtsevPetviashvili (KP) type, we prove that the zeromass constraint (in x) is satisfied at any non zero time even if it is not satisfied at initial time zero. Our results are based on a precise analysis of the fundamental solution of the linear part and of its anti xderivative. 

HOPW02 
30th March 2007 10:30 to 11:00 
High frequency scattering by convex polygons Standard boundary integral methods for high frequency scattering problems, with piecewise polynomial approximation spaces, suffer from the debilitating restriction that the number of degrees of freedom required to achieve a prescribed level of accuracy grows at least linearly with respect to the frequency. For problems of acoustic scattering by sound soft convex polygons, it has recently been demonstrated by ChandlerWilde and Langdon that by including plane wave basis functions supported on a graded mesh in the approximation space, with smaller elements closer to the corners of the polygon, convergence rates that depend only logarithmically on the frequency can be achieved. In this talk we review this approach, and consider some improvements and extensions to more complicated scattering problems. Related Links 

HOPW02 
30th March 2007 11:30 to 12:00 
JM SanzSerna 
Mollified impulse methods revisited We introduce a family of impulselike methods for the integration of highly oscillatory secondorder differential equations whose forces can be split into a fast and a slow part. Methods of this family are specified by two weight functions: one is used to average positions and the other to mollify the force. In cases where the fast forces are conservative, the new family inlcudes as particular cases the mollified impulse methods introduced by GarcíaArchilla, Skeel and the present author. On the other hand the methods here extend to nonlinear fast forces a wellknown class of exponential integrators introduced by Hairer and Lubich. A convergence analysis will be presented that provides insight into the role played by the processes of averaging and mollification. A simple condition on the weight functions is shown to be both necessary and sufficient to avoid order reduction. 

HOPW02 
30th March 2007 14:00 to 14:30 
On the localisation principle of high frequency scattering problems in integral equation discretisations The numerical simulation of scattering problems at moderate to high frequencies is a challenging problem. The discretization of any suitable mathematical model usually requires resolving the oscillations, which naturally leads to large and, in the context of integral equations, dense discretization matrices. Yet, high frequency scattering problems have a very local nature. The localization principle states that the reflection, refraction or diffraction of a wave that hits an obstacle is governed mainly by the geometry of the scattering object which is local to the point of contact. This principle is exploited by asymptotic methods, such as geometrical optics and the geometrical theory of diffraction. Suddenly, higher frequencies are desirable and lower frequencies become problematic. In this talk we examine how the localization principle can be exploited numerically in a more classical finite element setting. In particular, we may arrive at a sparse discretization matrix for integral equations by the use of Filontype quadrature rules for oscillatory integrals. We discuss the advantages and limitations of this approach and we examine the asymptotic nature of this sparse representation. Related Links


HOP 
10th April 2007 14:30 to 15:30 
G Dujardin  Normal form and long time analysis of splitting schemes for the linear Schroedinger equation  
HOP 
11th April 2007 14:30 to 15:30 
Symplectic integrators for highly oscillatory Hamiltonian systems  
HOP 
12th April 2007 14:30 to 15:30 
Stochastic Liegroup integrators  
HOP 
17th April 2007 14:30 to 15:30 
Symplectic Runge  Kutta Methods  
HOP 
18th April 2007 16:00 to 17:00 
E Zuazua  Wave propagation on planar graphs  
HOP 
25th April 2007 14:30 to 15:30 
A construction of mixed finite element spaces on polyhedral meshes  
HOP 
26th April 2007 11:30 to 12:30 
An efficient numerical solution of SturnLiouville problems  
HOP 
1st May 2007 14:30 to 15:30 
R Melnik  Timedependent models of phase transformations and their numerical approximations  
HOP 
2nd May 2007 14:30 to 15:30 
Improving Poincar\'e: exponential asymptotics for ODEs and integrals with saddles  
HOP 
3rd May 2007 14:30 to 15:30 
The higher order Stokes phenomenon: when is a Stokes line not a Stokes line?  
HOP 
8th May 2007 11:30 to 12:30 
Oscillations and rough paths  
HOP 
8th May 2007 14:30 to 15:30 
Dispersive numerical schemes for linear and nonlinear Schr\"odinger equations  
HOP 
9th May 2007 14:30 to 15:30 
Numerical computation of spacetime patterns of PDEs with symmetry  
HOP 
16th May 2007 14:30 to 15:30 
Liegroup variational integrators and their applications to geometric optimal control theory  
HOP 
17th May 2007 16:00 to 17:00 
A Its  Generalised FisherHartwig and Szeg\"oWidom asymptotics  
HOP 
18th May 2007 14:30 to 15:30 
J Vanneste  Inertiagravitywave generation by slow motion in the atmosphere and oceans  
HOP 
22nd May 2007 16:30 to 17:30 
Geometric integrators for the CamassaHolm equation  
HOP 
23rd May 2007 14:30 to 15:30 
Lifting the curse of dimensionality  numerical integrations in very high dimensions  
HOP 
24th May 2007 11:30 to 12:30 
A Ostermann  Exponential integrators for advection and diffusion problems  
HOP 
24th May 2007 14:30 to 15:30 
Approximate energy conservation for symplectic time semidiscretizations of semilinear wave equations  
HOP 
29th May 2007 14:30 to 15:30 
Fast multiscale denoising of surface meshes in 3D with intrinsic texture  
HOP 
30th May 2007 14:30 to 15:30 
Weak convergence, realisation of holonomic constraints, and the quantum adiabatic theorem  
HOP 
31st May 2007 14:30 to 15:30 
E Celledoni  Simulation of rigid body dynamics using Jacobi elliptic functions and integrals  
HOP 
7th June 2007 16:00 to 17:00 
A mathematical study of the CarParrinello method  
HOP 
11th June 2007 14:30 to 15:30 
P Olver  Invariantisation of numerical schemes  
HOP 
13th June 2007 11:30 to 12:30 
On the approximation of spectra and pseudospectra of nonselfadjoint Hilbert space operators  
HOP 
13th June 2007 14:30 to 15:30 
P Deift  Riemann Hilbert problems: theory and applications  
HOP 
14th June 2007 14:30 to 15:30 
M Ganesh  A spectrally accurate algorithm for electromagnetic scattering in three dimensions  
HOPW07 
19th June 2007 10:30 to 11:00 
Efficient evaluation of highfrequency propagation and scattering with applications to propagation in nonsphericallysymmetric atmospheres  
HOPW07 
19th June 2007 11:00 to 11:30 
High frequency scattering by simple ice crystal shapes  
HOPW07 
19th June 2007 11:30 to 12:00 
M Ganesh  An integral equation algorithm for surface scattering  
HOPW07 
19th June 2007 12:00 to 12:30 
P Childs  Aspects of seismic inversion  
HOPW07 
19th June 2007 14:00 to 14:30 
Asymptotics for highfrequency multiple scattering  
HOPW07 
19th June 2007 14:30 to 15:00 
Strategies for the evaluation of boundary integrals arising in 3D PUBEM for Helmholtz problems  
HOPW07 
19th June 2007 15:00 to 15:30 
Modelling of high intensity focused ultrasound propagation for cancer therapy  
HOPW07 
19th June 2007 16:00 to 16:30 
High frequency scattering problems from industrial applications  
HOP 
20th June 2007 14:30 to 15:00 
Reconsidering trigonometric integrators  
HOP 
20th June 2007 15:00 to 15:30 
M Khanamirian  Numerical solution to highly oscillatory systems of ODEs  
HOP 
20th June 2007 15:30 to 16:00 
From the wave packet estimate to exact observability  
HOP 
21st June 2007 16:00 to 17:00 
The fast evaluation of matrix functions for exponential integrators  
HOP 
22nd June 2007 15:30 to 16:30 
Homogenisation: Analysis and Computation  
HOP 
27th June 2007 14:30 to 15:00 
Interpolation in special orthogonal groups  
HOP 
27th June 2007 15:00 to 15:30 
Fourier analysis of the 1d discontinuous Galerkin methods  
HOP 
27th June 2007 15:30 to 16:00 
Spectral methods and modified Fourier series  
HOP 
28th June 2007 14:30 to 15:30 
VH Hoang  Sparse finite element method for nonlinear monotone problems with multiple scales  
HOPW05 
2nd July 2007 10:00 to 11:00 
Fast highorder highfrequency solvers in computational acoustics and electromagnetics The numerical solution of highly oscillatory wavepropagation and scattering problems presents a variety of significant challenges: these problems require high discretization densities and often give rise to poorly conditioned numerics; realistic engineering configurations, further, usually require consideration of geometries of great complexity and large extent. In this talk we will consider a number of methodologies that were introduced recently to address these difficulties. We will thus discuss algorithms that can solve, with highorder accuracy, problems of scattering for complex threedimensional geometriesincluding, possibly, singular elements such as wires, corners, edges and open screens. In particular, we will describe solutions achieved for two realistic threedimensional problems of very high frequencysurface scattering and atmospheric GPS propagationwhich previous threedimensional solvers could not address adequately. For added efficiency, these solvers, which are based on integral equations, highorder integration and fast Fourier transforms, can be used in conjunction with new regularized combined field equationswhich require much smaller numbers of iterations in a iterative linear algebra solver than combined field equations available previously. We will also describe a new class of highorder surface representation methods which, starting from point clouds or CAD data, can produce highorderaccurate surface parametrizations of complex engineering surfaces, as required by highorder solvers. Time permitting, applications of these methodologies to solution of timedomain problems and fast evaluation of fullynonlocal and convergent computational boundary conditions for timedomain problems will be mentioned. In all cases these algorithms exhibit highorder convergence, they run on low memories and reduced operation counts, and they can produce solutions with a high degree of accuracy. 

HOPW05 
2nd July 2007 11:30 to 12:00 
Domain decomposition for multiscale elliptic PDEs In the MonteCarlo simulation of fluid flow in stochastic media one typically has to solve elliptic PDEs with highly oscillatory coefficients. In practical applications in hydrogeology, these (random) oscillatory coefficients can have wavelength of order 10^{3} or smaller and amplitude of the order of 10^{8}. Accurate finite element approximation in 2D requires of the order of 10^{8} degrees of freedom and the resulting linear systems can have condition number close to 10^{16}. In this talk we will discuss domain decomposition preconditioning for such linear systems. Our overall aim is to solve such problems in a time close to the time required for solving a discretisation of a standard Poisson problem with constant coefficient with the same number of degrees of freedom. The essential step in constructing two or multilevel preconditioners is to replace the discretisation on the finest grid with a suitable discretisation on a coarser grid or grids. In the present application, standard (piecewise polynomial based) coarsening fails because, even if the fine mesh resolves the oscillations in the coefficient, coarser meshes typically fail to do so. By extending the classical domain decomposition theory to this case, we show that a suitable coarsening strategy for heterogeneous media involves the construction of low energy coarse space basis functions. This naturally suggests that multiscale finite element methods can provide good coarse spaces, and leads to a new class of domain decomposition preconditioners. Recent results of Scheichl and Vainikko used the same theoretical technique to explain the robustness of certain algebraically defined preconditioners. The theoretical results are illustrated by numerical examples on deterministic and random problems. Reference: I.G. Graham, P. Lechner and R. Scheichl, Domain Decomposition for Multiscale PDEs, Numer. Math. DOI 10.1007/s0021100700741 (2007) . 

HOPW05 
2nd July 2007 14:00 to 15:00 
E Zuazua 
An alternate direction descent method for the control of flows in the presence of shocks We present a new optimization strategy to compute numerical approximations of minimizers for optimal control problems governed by scalar conservation laws in the presence of shocks. We focus on the 1d inviscid Burgers equation. We first prove the existence of minimizers and, by a Gammaconvergence argument, the convergence of discrete minima obtained by means of numerical approximation schemes satisfying the so called onesided Lipschitz condition (OSLC). Then we address the problem of developing efficient descent algorithms. We first consider and compare the existing two possible approaches. The first one, the socalled discrete approach, based on a direct computation of gradients in the discrete problem and the socalled continuous one, where the discrete descent direction is obtained as a discrete copy of the continuous one. When optimal solutions have shock discontinuities, both approaches produce highly oscillating minimizing sequences and the effective descent rate is very weak. As a remedy we propose a new method that uses the recent developments of generalized tangent vectors and the linearization around discontinuous solutions. We develop a new descent stratagey, that we shall call "alternating descent method", distinguishing descent directions that move the shock and those that perturb the profile of the solution away of it. As we shall see, a suitable alternating combination of these two classes of descent directions allows building very efficient and fast descent algorithms. 

HOPW05 
2nd July 2007 15:30 to 16:00 
Condition number estimates for oscillatory integral operators We present results obtained during the programme on the conditioning of the standard combined potential boundary integral operators in timeharmonic acoustic scattering (cf. the talk by Oscar Bruno), in particular addressing behaviour as the wavenumber k tends to infinity, when the integral operator becomes increasingly oscillatory. While study of this topic goes back to Kress and Spassov (Numer. Math. 1983), the focus previously has been on the canonical case of a circle/sphere for which spherical harmonics are the eigenfunctions and the singular values are known explicitly. However, even for this case, it is only recently (Dominguez, Graham, Smyshlyaev, Preprint NI07004HOP and Numer. Math. 2007) that rigorous upper bounds have been obtained on the operator and its inverse as a function of k. For nonspherical scatterers the only result is a recent upper bound on the inverse operator for piecewise smooth starlike domains (ChandlerWilde & Monk, to appear SIAM J. Math. Anal.). In this talk we derive a range of lower bounds on the operator and its inverse, which show that the behaviour for large k depends subtly on the geometry. The sharpness of these lower bounds is demonstrated by numerical simulation and, in some instances, by provable upper bounds. The main computational message is that, while the condition number grows only mildly (like k^{1/3}) for a circle/sphere, behaviour can be much worse (like k^{5/4}) for nonstarlike domains. 

HOPW05 
2nd July 2007 16:00 to 16:30 
New Galerkin methods for highfrequency scattering simulations Recently developed integral equation methods for surface scattering simulations that combine the advantages of rigorous solvers (error controllability) with those of asymptotic methods (frequencyindependent discretizations) have displayed the potential of delivering scattering returns within a prescribed errortolerance in times that do not depend on the wavenumber k. The nearoptimal characteristics of these novel schemes, pioneered by Bruno, Geuzaine, Monro and Reitich, have rapidly generated great interest and significant number of work in recent years. In singlescattering configurations, an actual proof that provides a rigorous upper bound for the operation count of O(k^{1/9}) in the case of circular/spherical boundaries was recently established by Dominguez, Graham and Smyshlyaev for a pversion boundary element implementation of a similar approach. As in the original algorithm, they take profound advantage of the exponential decay (with increasing wavenumber k) of the corresponding physical density in the deep shadow region, and approximate this quantity by zero there. In this talk, we present two improved Galerkin schemes for the solution of singlescattering problems, and we show that, within the framework of our first scheme, the error in best approximation of the surface current grows at most at O(k^a) (for any a>0) over the "entire" boundary. Moreover, as we show, our second approach based on a novel change of variables around the "transition regions" reduces this dependency to O(log(k)). 

HOPW05 
2nd July 2007 16:30 to 17:00 
Plane wave discontinuous Galerkin methods Standard low order Lagrangian finite element discretization of boundary value problems for the Helmholtz equation $\Delta u  \omega^{2}u=f$ are afflicted with the socalled pollution phenomenon: though for sufficiently small $h\omega$ an accurate approximation of $u$ is possible, the Galerkin procedure fails to provide it. Attempts to remedy this have focused on incorporating extra information in the form of plane wave functions $\boldsymbol{x}\mapsto \exp(i\omega\boldsymbol{d}\cdot\boldsymbol{x})$, $\boldsymbol{d}=1$, into the trial spaces. Prominent examples of such methods are the plane wave partition of unity finite element method of Babuska and Melenk, and the ultraweak Galerkin discretization due to Cessenat and Despres. Both perform well in computations. It turns out that the latter method can be recast as a special socalled discontinuous Galerkin (DG) method employing local trial spaces spanned by a few plane waves. This perspective paves the way for marrying plane wave approximation with many of the various DG methods developed for 2ndorder elliptic boundary value problems. We have pursued this for a generic mixed DG method and a primal DG method which generalized the ultraweak scheme. For these methods we have developed a convergence analysis for the $h$version, which achieves convergence through mesh refinement. Key elements are approximation estimates for plane waves and sophisticated duality techniques. The latter entail estimating how well local plane waves can approximate the solution of a dual problem. Unfortunately, we could not help invoking general polynomial estimates in Sobolev spaces for this purpose. This incurs unsatisfactory pollutionaffected final error bounds $O(\omega^{2}h)$. On the other hand, a more detailed 1D analysis confirms that the plane wave DG method does not suffer from pollution. To bridge the gap is an open theoretical challenge. 

HOPW05 
3rd July 2007 09:00 to 10:00 
Integrators for Highly oscillatory Hamiltonian systems: an homogenisation approach We introduce a class of symplectic (and in fact also non symplectic) schemes for the numerical integration of highly oscillatory Hamiltonian systems. The bottom line for the approach is to exploit the HamiltonJacobi form of the equations of motion. Because we perform a twoscale expansion of the solution of the HamiltonJacobi equations itself, we readily obtain, after an appropriate discretization, symplectic integration schemes. Adequate modifications also provide non symplectic schemes. The efficiency of the approach is demonstrated using several variants. This is joint work with F. Legoll (LAMIENPC, France) 

HOPW05 
3rd July 2007 10:00 to 10:30 
Generalised polynomial chaos for analysing the stability If a dynamical system exhibits a periodic response, analysing the stability of this state yields crucial information. We consider oscillators modelled by systems of ordinary differential equations or differential algebraic equations. Hence we focus on stability properties of periodic solutions with respect to perturbations in corresponding initial values. Floquet theory represents a local concept for analysing the stability. Alternatively, we consider a stochastic perturbation following some probability distribution to obtain global information on stability. This strategy yields a system with stochastic input parameters. Thus results concerning the expected value and the variance of the corresponding solution are of interest. Monte Carlo methods can be used to compute these key figures approximately, where often a huge number of realisations is required. We apply the alternative approach of generalised polynomial chaos to obtain according approximations. Numerical simulations of oscillators using this strategy are presented. 

HOPW05 
3rd July 2007 11:00 to 11:30 
Adaptive parameterization in molecular dynamics I will describe work on formulation and numerical methods for molecular dynamics. Increasingly, in tackling complex models, it is necessary to manipulate the model "onthefly" during simulation, examples being quantummechanical/molecular mechanical models (QM/MM) and coarsegraining in molecular dynamics. I will discuss a new technique for adaptive parameterization of molecular landscapes. 

HOPW05 
3rd July 2007 11:30 to 12:00 
Splitting methods for oscillatory nonautonomous linear systems Several relevant physical phenomena are modelled (totally or partially) by linear PDEs that, once spatially discretized, give rise to systems of coupled harmonic oscillators. To diagonalise these systems is usually prohibitively expensive. Then, to use splitting methods (involving matrixvector products and possibly FFTs) is a valid alternative. A theoretical analysis about the stability and accuracy of splitting methods on the harmonic oscillator allows us to build new methods which outperform the existing methods from the literature. Nonautonomous problems are also of great interest but, the most efficient methods for the autonomous case are not valid in this setting. From the Magnus series expansion (as a formal solution to the nonautonomous problem) we show how to adapt these methods by treating the "time" separately from the coordinates. This technique allows us to build new methods whose performance is tested on the Schrödinger equation with timedependent potentials. Related Links 

HOPW05 
3rd July 2007 14:00 to 15:00 
Treatment of oscillations in classical molecular dynamics Atomic oscillations present in classical molecular dynamics (MD) restrict the step size that can be used. These oscillations are decidedly nonlinear, and there is only a modest separation of temporal scales. For reasons to be explained, amplitudes of the highest frequency oscillations are small, and, hence, MD is an example of a (mildly) stiff oscillatory problem. Actually, getting accurate trajectories is not usually the aim of MD; instead, time averages or time correlation functions are sought, and the implications of this are examined. Two techniques for lengthening the step size have found general use: multiple time stepping and bondlength constraints. Accuracy and stability of these techniques is discussed. Related Links


HOPW05 
3rd July 2007 15:30 to 16:00 
The Magnus expansion in the adiabatic picture In this talk we review the main aspects of the Magnus expansion for treating timedependent systems evolving in a nearadiabatic regime. We discuss its applicability and the convergence domain. We will show, in particular, how this procedure alllows us to obtain very accurate approximations to the transition probability for timedependent twostate quantum systems even far from the adiabatic limit. 

HOPW05 
3rd July 2007 16:00 to 16:30 
The separable shadow hybrid Monte Carlo (S2HMC) method for improved performence over hybrid Monte Carlo Hybrid Monte Carlo (HMC) is a rigorous sampling method that uses molecular dynamics (MD) as a global Monte Carlo move. The acceptance rate of HMC decays exponentially with system size. The Shadow Hybrid Monte Carlo (SHMC) was previously introduced to overcome this performance degradation by sampling instead from the shadow Hamiltonian defined for MD when using a symplectic integrator. However SHMC's performance is limited by the need to generate momenta for the MD step from a nonseparable shadow Hamiltonian. The Separable Shadow Hybrid Monte Carlo (S2HMC) method, based on a separable formulation of the shadow Hamiltonian, allows allows efficient generation of momenta and retains the advantage of SHMC. 

HOPW05 
3rd July 2007 16:45 to 17:15 
Level set methods for capturing semiclassical dynamics of Schr\"odinger equations with different potentials In this talk we present newly developed level methods for capturing semiclassical dynamics of Schr\"{o}dinger equations with different potentials. We discuss the essential ideas behind the techniques, the coupling of these techniques to handle several canonical potentials, including the phase space based level set method for given smooth potentials; the field space based level set method for selfconsistent potentials governed by the Poisson equation; as well as the Blochband based level set method for periodic potentials. The relations between computed multivalued solutions and desirable physical observables are established. Numerical examples are presented to validate the numerical methods. 

HOPW05 
3rd July 2007 17:15 to 17:45 
Adiabatic invariance and geometric phase in slowly deforming domains We consider the evolution of a 2d perfect fluid as its domain is deformed slowly in a prescribed fashion. Subject to some assumptions, the leadingorder Eulerian flow is found to be steady and depend only on the instantaneous form of the boundary; it is thus an adiabatic invariant of the system. Also to leading order, the Lagrangian particle trajectory is found to consist of a dynamical and a geometric component, in the fashion of the HannayBerry phase. Related Links


HOPW05 
4th July 2007 09:00 to 10:00 
MCMC methods for multiscale measure In many application domains there is growing interest in the fitting of stochastic differential equations (SDEs, diffusions) to data. The data may come from experimental observations or from large scale computer simulations. In many case the data has a multiscale character which is incompatible with a (or the desired) diffusion process at small scales. However it may be compatible at intermediate scales. In order to understand this situation I will study the fitting of SDEs to data generated by multiscale diffusion processes, in situations where averaging and homogenization apply. The parametric model will be the averaged or homogenized equation; the data, however, will be chosen from the multiscale model. Understanding this mismatch between data and model will shed light on the original problem of incompatibility between model and data at small scales. 

HOPW05 
4th July 2007 10:00 to 10:30 
Multiscale numerical methods using limited global information and applications In this talk, I will describe multiscale numerical methods for flows in heterogeneous porous media. The main idea of these methods is to construct local basis functions that can capture the small scale information when they are coupled via some global formulation. I will discuss the use of local boundary conditions, oversampling methods and the use of global information in constructing basis functions. Applications of these methods to stochastic equations will be also discussed. This is joint work with J. Aarnes, T. Hou, L. Jiang, V. Ginting. 

HOPW05 
4th July 2007 11:00 to 11:30 
A Boutet de Monvel 
Longtime asymptotics for the focusing NLS equation with timeperiodic boundary condition We consider the focusing NLS equation on the quarter plane. Initial data are vanishing at infinity while boundary date are timeperiodic (ae^[2i\omega t]). The main tool is the asymptotic analysis of the associated matrix RiemannHilbert problem. We will show that we obtain 4 different asymptotics in different regions: region 1: a ZakharovManakov vanishing asymptotics region 2: a train of asymptotics solitons region 3: a modulated elliptic wave region 4: a plane wave. Related Links 

HOPW05 
4th July 2007 11:30 to 12:00 
S Reich 
A Machuniform algorithm based on regularised fluid equations A Machuniform algorithm is an algorithm suitable for simulations at any level of the Mach number. Semiimplicit time discretization methods and pressure correction algorithms are among the popular choices to obtain Machuniform algorithms with applications to meteorology and combustion simulations. Here we suggest an alternative (although related) approach based on a regularized pressure formulation of the Euler equations. The equations are integrated in time by a staggered timestepping method which avoids artificial numerical damping of large scale accoustic waves. The regularization is motivated by a simple spatial multiscale analysis. Related Links


HOPW05 
4th July 2007 14:00 to 15:00 
A Abdulle 
SROCK: explicit methods for stiff stochastic problems Abstract: In this talk, we discuss a new class of methods for the solution of stiff stochastic differential equations. (This is a joint work with S. Cirilli). 

HOPW05 
4th July 2007 15:30 to 16:00 
S Smitheman 
The Dirichlet to Neumann map for the modified Helmholtz and Helmholtz equations with complex boundary data We present a spectral collocation type method for computing the Dirichlet to Neumann map for the modified Helmholtz equation. For regular and irregular polygons, we demonstrate quadratic convergence for sine basis functions and exponential convergence for Chebyshev basis functions. We go on to outline how our method can be extended to the Helmholtz equation, for which we also present numerical results. Our work is an extension of previous results of Prof. Fokas and collaborators for the Laplace equation (J. of Comput. and Appl. Maths. 167, 465483 (2004)). 

HOPW05 
4th July 2007 16:00 to 16:30 
Mathematical analysis and numerical simulation of BoseEinstein condensation In this talk, I review the mathematical results of the dynamcis of BoseEinstein condensate (BEC) and present some efficient and stable numerical methods to compute ground states and dynamics of BEC. As preparatory steps, we take the 3D GrossPitaevskii equation (GPE) with an angular momentum rotation, scale it to obtain a fourparameter model and show how to reduce it to 2D GPE in certain limiting regimes. Then we study numerically and asymptotically the ground states, excited states and quantized vortex states as well as their energy and chemical potential diagram in rotating BEC. Some very interesting numerical results are observed. Finally, we study numerically stability and interaction of quantized vortices in rotating BEC. Some interesting interaction patterns will be reported. Related Links


HOPW05 
4th July 2007 16:45 to 17:15 
Radial basis function collocation for Schr\"odinger's equation Gaussian radial basis functions may provide an efficient alternative to standard function spaces used in computational quantum mechanics. Combining them with a collocation spatial discretization may give simple and efficient integrators for the time dependent Schrodinger equation. My talk will describe an exploration of these possibilities. 

HOPW05 
4th July 2007 17:15 to 17:45 
A Zanna 
Exact and approximate methods for the free rigid body We consider both exact and approximate methods for the free rigid body equations. In particular, we are interested in understanding whether and when exact algorithms can be competitive with approximate ones. Both the exact and the approximate methods are extensively tested for several problems, including satellite and molecular dynamics simulations. 

HOPW05 
5th July 2007 09:00 to 10:00 
S Vandewalle 
Numerical integration of highly oscillatory functions based on analytic continuation We consider the integration of highly oscillatory functions. Based on analytic continuation, rapidly converging quadrature rules are derived for a fairly general class of oscillatory integrals with an analytic integrand. The accuracy of the quadrature increases both for the case of a fixed number of points and increasing frequency, and for the case of an increasing number of points and fixed frequency. These results are then used to obtain quadrature rules for more general oscillatory integrals, i.e., for functions that exhibit some smoothness but that are not analytic. The approach described in this paper is related to the steepest descent or saddle point method from complex analysis. However, it does not employ asymptotic expansions. It can be used for small or moderate frequencies as well as for very high frequencies. We consider both the onedimensional case, and the multidimensional case. Finally, we briefly elaborate on the use of the new integration rules in the context of solving highly oscillatory integral equations. 

HOPW05 
5th July 2007 10:00 to 10:30 
H Brunner 
Open problems in the computational solution of volterra functional equations with highly oscillatory solutions The approximation of solutions to Volterra integral equations by collocation or discontinuous Galerkin methods leads to a set of `semidiscretised' equations that in general are not amenable to numerical computation: an additional discretisation process that is able efficiently and accurately to cope with the highly oscillatory nature of the kernel of the given Volterra integral operator is needed. Here, the use of Filontype quadrature is an obvious possibility; however, it is not yet clear how best to do this when the kernel is weakly singular. In this talk I will describe current work related to the above problem, and especially to collocation methods for various types of Volterra functional equations, including equations with variable (and possibly vanishing) delay arguments. It will also be shown that in the case of Volterra integral equations of the first kind, the choice of the quadrature scheme in discontinuous Galerkin methods can have a major effect on the convergence properties of the approximate solution. 

HOPW05 
5th July 2007 11:00 to 11:30 
Asymptotic least squares approximation for highly oscillatory differential equations This talk presents a new approach for approximating highly oscillatory ordinary differential equations. By using the asymptotic expansion in a least squares system, we are able to obtain a result that preserves the asymptotic accuracy of the expansion, while converging rapidly to the exact solution. We are thus able to accurately approximate such differential equations by solving a very small linear system. We apply this method to the computation of highly oscillatory integrals, as well as second order oscillatory differential equations. 

HOPW05 
5th July 2007 11:30 to 12:00 
A Iserles 
Modified fourier expansions and spectral problems for highly oscillatory Fredholm operators Although highly oscillatory Fredholm operators are compact and have point spectrum, their calculation by standard means, e.g. the finite section method, is notoriously difficult. As an alternative, we propose expanding the underlying eigenfunctions in modified Fourier series. This leads to infinitedimensional algebraic eigenvalue problems that exhibit intriguing structure and rapid decay of coefficients. This is exploited in an effective numerical algorithm. 

HOPW05 
5th July 2007 14:00 to 15:00 
Approximation of nonadiabatic quantum dynamics by surface hopping The timedependent Schroedinger equation provides the fundamental description of quantum mechanical molecular dynamics. Its multiscale character suggests a splitting in two coupled subproblems, the so called BornOppenheimer approximation: One solves a family of stationary Schroedinger equations in the electronic degrees of freedom (one equation for each nucelonic configuration) and subsequently a timedependent Schroedinger equation in the nucleonic degrees of freedom, whose potential has been determined by the electronic problem. This splitting fails to provide an approximation, if different electronic eigenvalues are not uniformly separated for all nucleonic configurations. The talk explains a microlocal point of view on this nonadiabatic coupling between electronic and nucleonic degrees of freedom and derives an associated deterministic surface hopping algorithm. Its numerical realization crucially relies on the sampling of highly oscillatory initial data on highdimensional configuration spaces, which can be tackled by a Monte Carlo approach. The presented results are joint work with C. Fermanian, S. Kube, and M. Weber. 

HOPW05 
5th July 2007 15:30 to 16:00 
C Sparber 
Asymptotics for linear and nonlinear Schr\"odinger equations with periodic potentials We consider semiclassically scaled linear and nonlinear Schrödinger equation with highly oscillatory periodic potentials. We will discuss rigorous asymptotic descriptions for such problems as well as a recently developed numerical approach based on Bloch decomposition. 

HOPW05 
5th July 2007 16:00 to 16:30 
Computational quantum dynamics using Hagedorn wavepackets This talk, which is based on joint work with Erwan Faou and Vasile Gradinaru, reports on work in progress on a newly developed numerical approach to manybody quantum dynamics in the semiclassical regime. We present a symmetric splitting integrator for the propagation of multidimensional extensions of GaussHermite wavepackets appearing in analytical work by Hagedorn. The integrator evolves positions and momenta of the wavepackets according to the StoermerVerlet integrator of classical mechanics, and gains its computational feasibility and efficiency for many particles by the possibility of thinning out the moving basis sets according to a hyperbolic cross approximation or a Hartreetype approximation in a moving frame. The algorithm reduces to the Strang splitting of the Schroedinger equation in the limit of the full set of orthonormal basis functions, and it is robust in the semiclassical limit. 

HOPW05 
5th July 2007 16:45 to 17:15 
On macroscale variables for multiple time scale problems  
HOPW05 
5th July 2007 17:15 to 17:45 
A multiscale method for computing highly oscillatory ODEs We prsent a multiscale method for computing the effective behavior of a class of stiff and highly oscillatory ordinary differential equations. The oscillations may be in resonance with one another and thereby generate hidden slow dynamics. The proposed method relies on correctly tracking a set of slow variables whose dynamics is closed up to $\epsilon$ perturbation, and is sufficient to approximate any variable and functional that are slow under the dynamics of the ODE. This set of variables is detected numerically as a preprocessing step in the numerical methods. Error and complexity estimates are obtained. The advantages of the method is demonstrated with a few examples, including a commonly studied problem of Fermi, Pasta, and Ulam. 

HOPW05 
6th July 2007 09:00 to 10:00 
Partitioned RungeKutta methods in space and time In a partitioned RK (PRK) method, different coefficients are used for different variables. This allows, for example, explicit symplectic integrators for separable Hamiltonians; in fact, these are a special case of splitting methods. There also exist implicit symplectic PRK methods, but these seem to have no clear application, since non partitioned (Gauss) methods dominate them. Recently, a possible application has emerged for implicit symplectic PRK methods in the spatial discretization of Hamiltonian PDEs. Here they can yield stable multisymplectic methods, while by contrast, explicit PRK methods are unconditionally unstable. 

HOPW05 
6th July 2007 10:00 to 10:30 
High oscillations versus parasitic solutions The theory of modulated Fourier expansions is a powerful tool for the study of the longtime behaviour of differential equations with highly oscillatory solutions (conservation of energy, momentum, and harmonic actions). There is a discrete analogue that permits to study the longtime behaviour of linear multistep methods applied to (nonoscillatory) Hamiltonian systems. The parasitic solutions of the difference equations play the role of harmonic oscillations. In this talk we explain the common ideas of both theories. This is jointwork with Christian Lubich. Related Links 

HOPW05 
6th July 2007 11:00 to 11:30 
Towards symplectic Lie group integrators We will discuss recent progress towards analyzing and constructing symplectic Lie group integrators. Potential applications for Lie Poisson systems acted upon by the coadjoint action, and other systems with a symplectic group action. 

HOPW05 
6th July 2007 11:30 to 12:00 
Eulerian and semiLagrangian exponential integrators for convection dominated problems We consider a new class of integration methods particularly suited for convection diffusion problems with dominating convection. These methods are exponential integrators and their peculiarity is that they allow for the computation of exponentials of the linearized convection term. The main reason for developing this type of methods is that as it turns out they can be applied to the numerical integration of the considered PDEs in a semiLagrangian fashion. The main challenge in the numerical approximation of convection dominated phenomena is to avoid the occurrence of spurious oscillations in the numerical solution, (numerical dispersion), without adding diffusion. This task is achieved nicely by semiLagrangian methods. In these methods linear convective terms are integrated exactly by computing first the characteristics corresponding to the grid ponts of the adopted discretization, and then producing the numerical approximation via a simple although expensive interpolation procedure. 

HOPW05 
6th July 2007 14:00 to 14:30 
High order coefficient approximations methods for the SturmLiouville boundary value problem The determination of the eigenvalues of SturmLiouville problems is of great interest in mathematics and its applications. However most eigenvalue problems cannot be solved analytically, and computationally efficient approximation techniques are of great applicability. An important class of methods obtain eigenvalue approximations by applying an integrator based on coefficient approximation in a shooting process. These coefficient approximation methods replace the coefficient functions of the SturmLiouville equation by simpler approximations and then solve the approximating problem. The standard reference in the piecewise constant approximation case is due to S. Pruess [1], and therefore the methods are often referred to as Pruess methods. The Pruess method has some significant advantages. While a naive integrator is forced to make increasingly smaller steps in the search for large eigenvalues (due to the increasingly oscillatory nature of the solution), the stepsize is not restricted by the oscillations in the solution for a Pruess method. A drawback of the Pruess methods is the difficulty in obtaining higher order methods; unless Richardson extrapolation is used the method is only second order. Higher order methods based on coefficient approximation can be realized using a perturbation technique. This approach leads to the socalled Piecewise Perturbation Methods (PPM) [2]. The PPM add some perturbation corrections to the solution of the approximating problem in order to obtain a more accurate approximation to the solution of the original problem. High order PPM were found to be well suited to be used in a shooting procedure to compute the eigenvalues efficiently and accurately. This resulted in a Matlab software package which can be used to compute the eigenvalues of a SturmLiouville or Schr\"odinger problem up to high accuracy ({\sc MATSLISE} [3]). Recently it was shown that the piecewise perturbation approach may be viewed as the application of a modified Neumann expansion [4]. The excellent performance of piecewise perturbation methods for the SturmLiouville problem can thus be seen as a convincing illustration of the power and potential of the Neumann series integrators. Another integral series which has been recognized as a very effective computational tool for problems with highly oscillatory solution, is the Magnus expansion. Also integrators based on this Magnus expansion can be combined with coefficient approximation and form another extension of the Pruess ideas to high order approximations. \\ [1] Pruess, S. Estimating the eigenvalues of SturmLiouville problems by approximating the differential equation. SIAM J. Numer. Anal. 10 (1973). [2] Ixaru, L. Gr., De Meyer, H. and Vanden Berghe, G. CP methods for the Schr\"odinger equation revisited, J. Comput. Appl. Math. 88 (1997). [3] Ledoux, V., Van Daele, M., and Vanden Berghe, G. Matslise: A matlab package for the numerical solution of SturmLiouville and Schrodinger equations. ACM Trans. Math. Software 31 (2005). [4] Degani, I., AND Schiff, J. RCMS: Right Correction Magnus Series approach for oscillatory ODEs. J. Comput. Appl. Math. 193 (2006). 

HOPW05 
6th July 2007 14:30 to 15:00 
Numerical investigation of the conjugate locus for the Euler top Conjugate and cut loci of geodesic flows have a significant interest in Riemannian geometry, but few examples are known especially on manifolds of dimension greater than 2. An interesting example from mechanics is given by the flow of the Euler top, namely, a geodesic flow on SO(3) with a left invariant metric. The (first) conjugate locus can be determined analytically if two of the three moments of inertia of the body are equal. In that case, the conjugate locus is either a segment or circle (if the body is oblate) or a noninjective mapping of an astroid of revolution (if the body is prolate) [Bates and Fasso` 2006, see the link below]. The analytic study of the generic case of distinct moments of inertia is much more difficult, if not even prohibitive. We thus resort to the numerical construction of the conjugate locus, based on accurate numerical integrations of the flow and of its tangent map. The dependency of the conjugate locus on the moments of inertia is studied in a deformation scenario from the symmetric case. (Joint collaboration with L. Bates). Related Links


HOPW05 
6th July 2007 15:30 to 16:00 
Lie group and homogeneous variational integrators and their applications to geometric optimal control theory The geometric approach to mechanics serves as the theoretical underpinning of innovative control methodologies in geometric control theory. These techniques allow the attitude of satellites to be controlled using changes in its shape, as opposed to chemical propulsion, and are the basis for understanding the ability of a falling cat to always land on its feet, even when released in an inverted orientation. We will discuss the application of geometric structurepreserving numerical schemes to the optimal control of mechanical systems. In particular, we consider Lie group variational integrators, which are based on a discretization of Hamilton's principle that preserves the Lie group structure of the configuration space. In contrast to traditional Lie group integrators, issues of equivariance and orderofaccuracy are independent of the choice of retraction in the variational formulation. The importance of simultaneously preserving the symplectic and Lie group properties is also demonstrated. In addition, we will introduce a numerically robust shooting based optimization algorithm that relies on the conservation properties of geometric integrators to accurately compute sensitivity derivatives, thereby yielding an optimization algorithm for the control of mechanical systems that is exceptionally efficient. The role of geometric phases in these control algorithms will also be addressed. Recent extensions to homogeneous spaces yield intrinsic methods for Hamiltonian flows on the sphere, and have potential applications to the simulation of geometric exact rods, structures and mechanisms. The research has been supported in part by NSF grant DMS0726263 and DMS0504747. Related Links


HOPW05 
6th July 2007 16:00 to 16:30 
R Melnik 
Spurious solutions in applications of effective envelope theory: roots of the problem and methodologies for its remedy The appearance of spurious solutions is a long standing problem in the application of the envelope function approximation. The subject of this talk is the exact envelope function theory developed by M.G. Burt and further extended by B.A. Foreman. Using this theory, multiband models can be derived for nanoscale heterostructures, similar to the effective mass models developed by J.M. Luttinger and W. Kohn that opened the door to modern treatments of important classes of bulk semiconductor materials and nanostructures. In applying this approach, we have to undergo several levels of approximations of the Hamiltonian. In this talk we analyze such approximations critically and show how spurious nonphysical solutions in the resulting coupled systems of partial differential equations can be avoided. The discussion is motivated by our research into lowdimensional semiconductor nanostructures, in particular quantum dots. Related Links 

HOPW05 
6th July 2007 16:30 to 17:00 
Backscatter diffraction coefficients for surfacebreaking cracks: Comparison of two semianalytical approaches Evaluation of diffraction coefficients for surfacebreaking cracks has been an outstanding problem in mathematical modelling of ultrasonic inspection for over twenty years. Several groups of applied mathematicians have been involved in developing suitable semianalytical solutions over an even longer period. Recently two different approaches, one based on the Budaev and Bogy's technique and one based on the Gautesen's method have been developed and coded. We report on their successful crossvalidation and validation against experimental data. Related Links


HOPW08 
13th September 2010 08:55 to 09:00 
Welcome by ESF Rapporteur  
HOPW08 
13th September 2010 09:00 to 10:00 
Non polynomial approximations of wave problems  
HOPW08 
13th September 2010 10:00 to 11:00 
Computing Slowly Advancing Features in FastSlow Systems without Scale Separation  A Young Measure Approach  
HOPW08 
13th September 2010 11:30 to 12:30 
Fourier series on triangles and tetrahedra  
HOPW08 
13th September 2010 14:00 to 14:30 
L Gauckler  Modulated Fourier expansions for the longtime analysis of Hamiltonian PDEs  
HOPW08 
13th September 2010 14:30 to 15:00 
R Quispel  IntegralPreserving Integrators  
HOPW08 
13th September 2010 15:00 to 15:30 
On the numerical solution of firstkind Volterra integral equations with highly oscillatory kernels  
HOPW08 
13th September 2010 16:00 to 17:00 
Multiple integrals under differential constraints: twoscale convergence and homogenization  
HOPW08 
13th September 2010 17:00 to 18:00 
On Wigner and Bohmian Measures  
HOPW08 
13th September 2010 18:00 to 18:30 
E Karatsuba  On approximation of special oscillating series of quantum theory  
HOPW08 
14th September 2010 09:00 to 10:00 
A Addulle  Numerical integration for numerical homogenizations methods  
HOPW08 
14th September 2010 10:00 to 11:00 
Multiscale spaces for multiscale highcontrast problems and their applications  
HOPW08 
14th September 2010 11:30 to 12:30 
Computing expectation values for molecular quantum dynamics  
HOPW08 
14th September 2010 14:00 to 14:30 
Variance reduction in stochastic homogenization using antithetic variables  
HOPW08 
14th September 2010 14:30 to 15:00 
Oscillatory systems with three separated timescales  
HOPW08 
14th September 2010 15:00 to 15:30 
Developing Integrators for Highly Oscillatory Hamiltonian Systems Using Homogenization  
HOPW08 
14th September 2010 16:30 to 17:00 
Preconditioning the Helmholtz equation using rapidly oscillating functions  
HOPW08 
14th September 2010 17:00 to 17:30 
Coercivity, Nonnormality and Numerical Range of boundary integral operators in highfrequency scattering  
HOPW08 
14th September 2010 17:30 to 18:00 
Convergent highfrequency algorithms for single and multiple scattering  
HOPW08 
14th September 2010 18:00 to 18:30 
RungaKutta convolution quadrature: convergence theory and applications to timedomain BIE of acoustic scattering  
HOPW08 
15th September 2010 09:00 to 10:00 
Modulated Fourier expansions and the FermiPastaUlam problem  
HOPW08 
15th September 2010 10:00 to 11:00 
T Lelievre  Metastability, rare events and sampling problems in molecular dynamics  
HOPW08 
15th September 2010 11:30 to 12:00 
Approximation by plane waves  
HOPW08 
15th September 2010 12:00 to 12:30 
The generalised eigenfunction method and timedependent linear waterwave impact on a vertical elastic plate  
HOPW08 
15th September 2010 14:00 to 14:30 
Numerical solution of Riemann–Hilbert problems: Painlevé II  
HOPW08 
15th September 2010 14:30 to 15:00 
J Geier  Efficient finite difference schemes for highly oscillatory linear ODE  
HOPW08 
15th September 2010 15:00 to 15:30 
Highly Oscillatory Integrals and their Applications  
HOPW08 
15th September 2010 16:00 to 17:00 
C SanzSerna 
Optimal Tuning of Hybrid Monte Carlo
(This is joint work with A. Beskos, N.S. Pillai, G.O. Roberts and A.M. Stuart).


HOPW08 
15th September 2010 17:00 to 18:00 
The Hybrid Monte Carlo Algorithm in High Dimensions  
HOPW08 
15th September 2010 18:00 to 18:30 
Error Estimates for Gaussian Beam Superposition  
HOPW08 
16th September 2010 09:00 to 10:00 
H Owhadi  Homogenization with nonseparated scales and high contrast  
HOPW08 
16th September 2010 10:00 to 11:00 
A new transform method and some of its numerical implementations  
HOPW08 
16th September 2010 11:30 to 12:00 
CB Schoenlieb  Higherorder total variation for oscillating patterns  
HOPW08 
16th September 2010 12:00 to 12:30 
Oscillation Analysis of Some Hybrid Dynamical Systems of Transmission Pipelines  
HOPW08 
16th September 2010 14:00 to 14:30 
Tensor product method for radiative transfer  
HOPW08 
16th September 2010 14:30 to 15:00 
W Hao  Bloch DecompositionBased Gaussian Beam Method for the Schrodinger equation with Periodic Potentials  
HOPW08 
16th September 2010 15:00 to 15:30 
Strichartz estimates for the Schroedinger equation on a tree and applications  
HOPW08 
16th September 2010 16:00 to 16:30 
Upscaling of network models for porous media  
HOPW08 
16th September 2010 16:30 to 17:00 
Thin domains with highly oscillating boundaries  
HOPW08 
16th September 2010 17:00 to 17:30 
Oscillatory diffusion equation  
HOPW08 
16th September 2010 17:30 to 18:00 
An Effective Finite Difference Approach for Optical Beam Propagations  
HOPW08 
16th September 2010 18:00 to 18:30 
Global and Blowup patterns of the Cauchy problem of a fourthorder thin film equation  
HOPW08 
17th September 2010 09:00 to 10:00 
Applications of High Frequency Methods in Aeronautics: status and needs  
HOPW08 
17th September 2010 10:10 to 10:40 
Nonlinear hyperbolic conservation laws: diffusivedispersive limits  
HOPW08 
17th September 2010 10:40 to 11:10 
Coercivity of boundary integral operators in high frequency scattering  
HOPW08 
17th September 2010 11:30 to 12:30 
Existence of approximate solitary waves in simplectic algorithms of integration
(coauthors: Erwan Faou and Benoit Grebert)
