# Past Computational Mathematics and Applications Seminar

9 March 2000
14:00
Abstract
This lecture is about the extension of our CAD-Free platform to the simulation and sensitivity analysis for design and control of multi-model configurations. We present the different ingredients of the platform using simple models for the physic of the problem. This presentation enables for an easy evaluation of coupling, design and control strategies. \\ \\ Sensitivity analysis has been performed using automatic differentiation in direct or reverse modes and finite difference or complex variable methods. This former approach is interesting for cases where only one control parameter is involved as we can evaluate the state and sensitivity in real time with only one evaluation. \\ \\ We show that for some class of applications, incomplete sensitivities can be evaluated dropping the state dependency which leads to a drastic cost reduction. \\ \\ The design concerns the structural characteristic of the model and our control approach is based in perturbating the inflow incidence by a time dependent flap position described by a dynamic system encapsulating a gradient based minimization algorithm expressed as dynamic system. This approach enables for the definition of various control laws for different minimization algorithm. \\ \\ We present dynamic minimization algorithms based on the coupling of several heavy ball method. This approach enables for global minimization at a cost proportional to the number of balls times the cost of one steepest descent minimization. \\ \\ Two and three dimension flutter problem simulation and control are presented and the sensitivity of the different parameters in the model is discussed.
• Computational Mathematics and Applications Seminar
24 February 2000
14:00
--
Abstract
This seminar has been cancelled.
• Computational Mathematics and Applications Seminar
Dr Kurt Lust
Abstract
There is a growing interest in the study of periodic phenomena in large-scale nonlinear dynamical systems. Often the high-dimensional system has only low-dimensional dynamics, e.g., many reaction-diffusion systems or Navier-Stokes flows at low Reynolds number. We present an approach that exploits this property in order to compute branches of periodic solutions of the large system of ordinary differential equations (ODEs) obtained after a space discretisation of the PDE. We call our approach the Newton-Picard method. Our method is based on the recursive projection method (RPM) of Shroff and Keller but extends this method in many different ways. Our technique tries to combine the performance of straightforward time integration with the advantages of solving a nonlinear boundary value problem using Newton's method and a direct solver. Time integration works well for very stable limit cycles. Solving a boundary value problem is expensive, but works also for unstable limit cycles. \\ \\ We will present some background material on RPM. Next we will explain the basic features of the Newton-Picard method for single shooting. The linearised system is solved by a combination of direct and iterative techniques. First, we isolate the low-dimensional subspace of unstable and weakly stable modes (using orthogonal subspace iteration) and project the linearised system on this subspace and on its (high-dimensional) orthogonal complement. In the high-dimensional subspace we use iterative techniques such as Picard iteration or GMRES. In the low-dimensional (but "hard") subspace, direct methods such as Gaussian elimination or a least-squares are used. While computing the projectors, we also obtain good estimates for the dominant, stability-determining Floquet multipliers. We will present a framework that allows us to monitor and steer the convergence behaviour of the method. \\ \\ RPM and the Newton-Picard technique have been developed for PDEs that reduce to large systems of ODEs after space discretisation. In fact, both methods can be applied to any large system of ODEs. We will indicate how these methods can be applied to the discretisation of the Navier-Stokes equations for incompressible flow (which reduce to an index-2 system of differential-algebraic equations after space discretisation when written in terms of velocity and pressure.) \\ \\ The Newton-Picard method has already been extended to the computation of bifurcation points on paths of periodic solutions and to multiple shooting. Extension to certain collocation and finite difference techniques is also possible.
• Computational Mathematics and Applications Seminar
10 February 2000
14:00
Prof Tom Hughes
Abstract
• Computational Mathematics and Applications Seminar
3 February 2000
14:00
Prof Henk van der Vorst
Abstract
Krylov subspace methods offer good possibilities for the solution of large sparse linear systems of equations.For general systems, some of the popular methods often show an irregular type of convergence behavior and one may wonder whether that could be improved or not. Many suggestions have been made for improvement and the question arises whether these corrections are cosmetic or not. There is also the question whether the irregularity shows inherent numerical instability. In such cases one should take extra care in the application of smoothing techniques. We will discuss strategies that work well and strategies that might have been expected to work well.
• Computational Mathematics and Applications Seminar
27 January 2000
14:00
Prof Neil Sandham
Abstract
This work forms part of a larger research project to develop efficient low-dissipative high-order numerical techniques for high-speed turbulent flow simulation, including shock wave interactions with turbulence. The requirements on a numerical method are stringent.For the turbulence the method must be capable of resolving accurately a wide range of length scales, whilst for shock waves the method must be stable and not generate excessive local oscillations. Conventional methods are either too dissipative, or incapable of shock capturing. Higher-order ENO, WENO or hybrid schemes are too expensive for practical computations. Previous work of Yee, Sandham & Djomehri (1999) developed high-order shock-capturing schemes which minimize the use of numerical dissipation away from shock waves. The objective of the present study is to further minimize the use of numerical dissipation for shock-free compressible turbulence simulations.
• Computational Mathematics and Applications Seminar
20 January 2000
14:00
Prof Bruce Christianson
Abstract
In 1983 Pantoja described a stagewise construction of the exact Newton direction for a discrete time optimal control problem. His algorithm requires the solution of linear equations with coefficients given by recurrences involving second derivatives, for which accurate values are therefore required. \\ \\ Automatic differentiation is a set of techniques for obtaining derivatives of functions which are calculated by a program, including loops and subroutine calls, by transforming the text of the program. \\ \\ In this talk we show how automatic differentiation can be used to evaluate exactly the quantities required by Pantoja's algorithm, thus avoiding the labour of forming and differentiating adjoint equations by hand. \\ \\ The cost of calculating the newton direction amounts to the cost of solving one set of linear equations, of the order of the number of control variables, for each time step. The working storage cost can be made smaller than that required to hold the solution.
• Computational Mathematics and Applications Seminar
2 December 1999
15:00
--
Abstract
No seminar today
• Computational Mathematics and Applications Seminar