Fractional Spectral and Spectral Element Methods

Size: px
Start display at page:

Download "Fractional Spectral and Spectral Element Methods"

Transcription

1 Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis of Mohsen Zayernouri) George Em Karniadakis Division of Applied Mathematics, Brown University & Department of Mechanical Engineering, MIT & Pacific Northwest National Laboratory, CM4/DoE The CRUNCH group:

2

3 Outline Fractional Sturm-Liouville Problem (FSLP): Singular & Regular, each of kind-i and kind-ii Fractional Ordinary Differential Equations (FODEs): PG Spectral Methods & Discontinuous Spectral/hp Element Method Discontinuous Petrov-Galerkin (DPG) for Fractional PDEs: Notion of exponentially-accurate time-integration Fractional Collocation Method (Linear/Nonlinear/Multi-Term FPDEs): A new class of fractional Lagrange interpolants A Unified Spectral Method for FPDEs: Hyperbolic, Parabolic, and Elliptic FPDEs in (1+d)-dimensions Conclusion

4 Singular Sturm-Liouville Problem (Integer-Order): Jacques François Sturm ( ) Joseph Liouville ( ) Local Operator Local Boundary Conditions, (Jacobi polynomials), (Quadratic)

5 Singular Fractional Sturm-Liouville problem: Global Operator Non-local Boundary Conditions i =1: SFSLP of Kind-I i =2: SFSLP of Kind-II

6 Singular Fractional Sturm-Liouville problem: Jacobi Polyfractonomials : Theorem: The exact eigenfunctions of SFSLP-I (i=1) and SFSLP-II (i=2) are: and correspondingly, the exact eigenvalues are

7 Orthogonality Properties: Recurrence relations Fractional derivatives

8 Recall: SFSLP-I: SFSLP-II: What if? Then, FSLP is no longer singular and it becomes Regular

9 Regular Fractional Sturm-Liouville problem (RFSLP): Regular Fractional Operator: for RFSLP-I for RFSLP-II

10 Theorem: The exact eigenfunctions of RFSLP-I (i=1) and RFSLP-II (i=2) are: Regular Fractional Sturm-Liouville problem (RFSLP): Jacobi Polyfractonomials : and correspondingly the exact eigenvalues are:

11 Regular Eigenvalues

12 Eigen-solutions of RFSLP-I Eigen-solutions of RFSLP-II The same number of zeros Sharp gradient near Dirichlet end

13 Approximation Properties of the Polyfractonomials

14 Fractional Ordinary Differential Equations (FODEs) Model Problem 1: Fractional Initial-Value Problem (RL derivative) Model Problem 2: Fractional Delay Equation (RL derivative)

15 Petrov-Galerkin Spectral Method (FODEs) BASIS functions: (eigenfunctions of FSLP -I) Taking ( ) Property: TEST functions: (eigenfunctions of FSLP -II) Property:

16 Petrov-Galerkin Spectral Method (FODEs) Model Problem 1: Fractional Integration-by-Parts Taking Diagonal Stiffness Matrix!

17 Exponential Decay of the L2-Error

18 Fractional Ordinary Differential Equations (FODEs) What if the solution is Not Smooth, or, it is Piecewise smooth? How to treat Long-Time Integration?

19 Discontinuous Spectral/hp Element Method (FODEs) Partitioning the time-domain into non-overlapping elements: Time BASIS: TEST: How to construct them?

20 Discontinuous Spectral/hp Element Method (FODEs) BASIS functions: TEST functions:

21 Discontinuous Spectral/hp Element Method (FODEs) Consider the general case: Taking

22 Discontinuous Spectral/hp Element Method (FODEs) Fractional Integration-by-Parts : History Load Diagonal stiffness Mass and the delay Mass matrices are exact

23 A=B=0, Non-Smooth Solution: h-ref More Effective than p-ref

24 CPU Time: PG Methods vs. FDM A=B=0

25 To fade or not to fade?

26 Fractional Spectral Collocation Method How to do even better?! Thinking of Employing Nodal rather than previously used Modal expansions Avoiding Quadrature costs (in the construction of the matrices) Efficient treatment of Multi-Term fractional FPDEs Efficient treatment of Nonlinear FPDEs

27 Fractional Spectral Collocation Method Nodal Expansion: : Fractional Lagrange Interpolants at some arbitrary interpolation/collocation points: We note the difference between and Polynomial/Fourier interpolants!

28 Fractional Spectral Collocation Method Differentiation Matrices, : Linear Systems

29 II) Space-Fractional Nonlinear Burgers Equation:

30 A Unified Spectral Method for FPDEs A Unified Spectral Method, which treats any (1+d)-dimensional FPDE (i.e., hyperbolic, parabolic, and elliptic) with the same ease. Time-Space Hypercube and To also develop a Unified Fast Solver for this general scheme!

31 A Unified Spectral Method for FPDEs Basis Space: Based on Polyfractonomials of FSLP-I

32 A Unified Spectral Method for FPDEs Test Space: Based on Polyfractonomials of FSLP-II

33 A Unified Spectral Method for FPDEs General Lyapunov System

34 A Unified Spectral Method for FPDEs This choice of basis/test functions leads to Diagonal Temporal Stiffness Diagonal Spatial Stiffness Tridiagonal Spatial Stiffness when when Symmetric Temporal/Spatial Mass Matrices and

35 A Unified Spectral Method for FPDEs Theorem. A closed-form solution for the multi-dimensional matrix of unknown coefficient matrix is explicitly given as Remark. The fractional order is usually the same in all dimensions, so the eigenvalue problem is usually solved only once for all!

36 Hyperbolic FPDEs: A Unified Spectral Method for FPDEs

37 Parabolic FPDEs: A Unified Spectral Method for FPDEs

38 Concluding Remarks Exact non-polynomial eigenfunctions of the singular/regular fractional Sturm-Liouville problems were obtained and introduced as new BASIS and TEST functions. A Petrov-Galerkin spectral method and a discontinuous spectral element method were developed for FODEs (FIVPs and Delay Problems) An exponentially accurate DPG method for FPDES, also for time-integration method was proposed for first-order in time operators was developed. A fractional spectral collocation method was developed for efficient treatment of multi-term and nonlinear FPDEs A unified spectral method along with a general fast linear solver was developed for FPDEs including hyperbolic, parabolic and elliptic problems Less memory storage and CPU time in the proposed schemes Exponential decay of error was obtained in several test-cases for FODEs and FPDEs.

39 Special Issue of Journal of Computational Physics on FPDEs 20 papers Submission: April 30 th, 2014 Production: August 30 th, 2014

40 Publications 1) M. Zayernouri and G.E. Karniadakis, Fractional Sturm-Liouville Eigen-Problems: Theory and Numerical Approximation, J. Comput. Phys. vol. 252, (2013), Pages ) M. Zayernouri and G.E. Karniadakis, Exponentially Accurate Spectral and Spectral Element Methods for Fractional ODEs, J. Comput. Phys. vol. 257-Part A (2014), Pages ) M. Zayernouri and G.E. Karniadakis, Fractional Spectral Collocation Method, SIAM J. Scientific Computing, (2013), Accepted. 4) M. Zayernouri, W. Cao, and G.E. Karniadakis, Spectral and Discontinuous Spectral Element Methods for Fractional Delay Equations, Submitted to SIAM J. Scientific Computing, (2013). 5) M. Zayernouri, I. Alvey, and G.E. Karniadakis, Discontinuous Petrov-Galerkin (DPG) Methods for Time- & Space-Fractional Advection Equation, Submitted to SIAM J. Scientific Computing, (2013). 6) M. Zayernouri, M. Ainsworth, and G.E. Karniadakis, A Unified Petrov-Galerkin Spectral Method for Fractional PDEs, to be Submitted to CMAME, (2013).

41 I) Time- and Space-Fractional Multi-Term FPDEs : Lyapunov equation:

Spectral and Spectral Element Methods for Fractional PDEs

Spectral and Spectral Element Methods for Fractional PDEs Spectral and Spectral Element Methods for Fractional PDEs by Mohsen Zayernouri B.Sc., Mechanical Engineering, Azad University, Iran, 2004 M.Sc., Mechanical Engineering, Tehran Polytechnic, Iran, 2006 Ph.D.,

More information

Downloaded 01/21/14 to Redistribution subject to SIAM license or copyright; see

Downloaded 01/21/14 to Redistribution subject to SIAM license or copyright; see SIAM J. SCI. COMPUT. Vol. 6, o. 1, pp. A40 A62 c 2014 Society for Industrial and Applied Mathematics FRACTIOAL SPECTRAL COLLOCATIO METHOD MOHSE ZAYEROURI AD GEORGE EM KARIADAKIS Abstract. We develop an

More information

Fast Structured Spectral Methods

Fast Structured Spectral Methods Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia

More information

Contents as of 12/8/2017. Preface. 1. Overview...1

Contents as of 12/8/2017. Preface. 1. Overview...1 Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite

More information

Lecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.

Lecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form. Lecture IX Abstract When solving PDEs it is often necessary to represent the solution in terms of a series of orthogonal functions. One way to obtain an orthogonal family of functions is by solving a particular

More information

arxiv: v1 [math.na] 27 Oct 2016

arxiv: v1 [math.na] 27 Oct 2016 A Petrov-Galerin Spectral Element Method for Fractional Elliptic Problems arxiv:1610.08608v1 [math.na] 7 Oct 016 Ehsan Kharazmi a,b, Mohsen Zayernouri a,b, and George Em Karniadais c a Department of Computational

More information

Index. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2

Index. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2 Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604

More information

Numerical Mathematics

Numerical Mathematics Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector

More information

Final Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications

Final Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications Final Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications PI: George Em Karniadakis Division of Applied Mathematics, Brown University April 25, 2005 1 Objectives

More information

256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.

256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim. 56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =

More information

Introduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA

Introduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC

More information

Applied Numerical Analysis

Applied Numerical Analysis Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory

More information

Overlapping Schwarz preconditioners for Fekete spectral elements

Overlapping Schwarz preconditioners for Fekete spectral elements Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,

More information

PARTIAL DIFFERENTIAL EQUATIONS

PARTIAL DIFFERENTIAL EQUATIONS MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL

More information

FRACTIONAL-ORDER UNIAXIAL VISCO-ELASTO-PLASTIC MODELS FOR STRUCTURAL ANALYSIS

FRACTIONAL-ORDER UNIAXIAL VISCO-ELASTO-PLASTIC MODELS FOR STRUCTURAL ANALYSIS A fractional-order uniaxial visco-elasto-plastic model for structural analysis XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J. Owen,

More information

Downloaded 08/06/14 to Redistribution subject to SIAM license or copyright; see

Downloaded 08/06/14 to Redistribution subject to SIAM license or copyright; see SIAM J. SCI. COMPUT. Vol. 36, No. 4, pp. B684 B707 c 2014 Society for Industrial and Applied Mathematics DISCONTINUOUS SPECTRAL ELEMENT METHODS FOR TIME- AND SPACE-FRACTIONAL ADVECTION EQUATIONS MOHSEN

More information

Preface. Figures Figures appearing in the text were prepared using MATLAB R. For product information, please contact:

Preface. Figures Figures appearing in the text were prepared using MATLAB R. For product information, please contact: Linear algebra forms the basis for much of modern mathematics theoretical, applied, and computational. The purpose of this book is to provide a broad and solid foundation for the study of advanced mathematics.

More information

MATHEMATICAL METHODS INTERPOLATION

MATHEMATICAL METHODS INTERPOLATION MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU

More information

PART IV Spectral Methods

PART IV Spectral Methods PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,

More information

Beyond Wiener Askey Expansions: Handling Arbitrary PDFs

Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1

More information

Partial Differential Equations with Numerical Methods

Partial Differential Equations with Numerical Methods Stig Larsson Vidar Thomée Partial Differential Equations with Numerical Methods May 2, 2003 Springer Berlin Heidelberg New York Barcelona Hong Kong London Milan Paris Tokyo Preface Our purpose in this

More information

Index. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems

Index. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,

More information

Spectral Processing. Misha Kazhdan

Spectral Processing. Misha Kazhdan Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry

More information

Upon successful completion of MATH 220, the student will be able to:

Upon successful completion of MATH 220, the student will be able to: MATH 220 Matrices Upon successful completion of MATH 220, the student will be able to: 1. Identify a system of linear equations (or linear system) and describe its solution set 2. Write down the coefficient

More information

NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING

NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical

More information

Performance Evaluation of Generalized Polynomial Chaos

Performance Evaluation of Generalized Polynomial Chaos Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu

More information

Chapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs

Chapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u

More information

Lecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.

Lecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C. Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference

More information

Numerical Methods for Engineers and Scientists

Numerical Methods for Engineers and Scientists Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL

More information

NUMERICAL MATHEMATICS AND COMPUTING

NUMERICAL MATHEMATICS AND COMPUTING NUMERICAL MATHEMATICS AND COMPUTING Fourth Edition Ward Cheney David Kincaid The University of Texas at Austin 9 Brooks/Cole Publishing Company I(T)P An International Thomson Publishing Company Pacific

More information

Partial Differential Equations

Partial Differential Equations Next: Using Matlab Up: Numerical Analysis for Chemical Previous: Ordinary Differential Equations Subsections Finite Difference: Elliptic Equations The Laplace Equations Solution Techniques Boundary Conditions

More information

Special Functions of Mathematical Physics

Special Functions of Mathematical Physics Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table

More information

Applied Linear Algebra

Applied Linear Algebra Applied Linear Algebra Peter J. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@math.umn.edu http://www.math.umn.edu/ olver Chehrzad Shakiban Department of Mathematics University

More information

Preface. 2 Linear Equations and Eigenvalue Problem 22

Preface. 2 Linear Equations and Eigenvalue Problem 22 Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits

More information

1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel

1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel 1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,

More information

A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere

A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background

More information

Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications

Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535

More information

Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs

Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation

More information

Rational Chebyshev pseudospectral method for long-short wave equations

Rational Chebyshev pseudospectral method for long-short wave equations Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84

More information

FINITE-DIMENSIONAL LINEAR ALGEBRA

FINITE-DIMENSIONAL LINEAR ALGEBRA DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup

More information

NUMERICAL METHODS FOR ENGINEERING APPLICATION

NUMERICAL METHODS FOR ENGINEERING APPLICATION NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto

More information

ADVANCED ENGINEERING MATHEMATICS

ADVANCED ENGINEERING MATHEMATICS ADVANCED ENGINEERING MATHEMATICS DENNIS G. ZILL Loyola Marymount University MICHAEL R. CULLEN Loyola Marymount University PWS-KENT O I^7 3 PUBLISHING COMPANY E 9 U Boston CONTENTS Preface xiii Parti ORDINARY

More information

Journal of Computational Physics

Journal of Computational Physics Journal of Computational Physics 23 (211) 8573 862 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp An eigen-based high-order

More information

Computing Spectra of Linear Operators Using Finite Differences

Computing Spectra of Linear Operators Using Finite Differences Computing Spectra of Linear Operators Using Finite Differences J. Nathan Kutz Department of Applied Mathematics University of Washington Seattle, WA 98195-2420 Email: (kutz@amath.washington.edu) Stability

More information

HI CAMBRIDGE n S P UNIVERSITY PRESS

HI CAMBRIDGE n S P UNIVERSITY PRESS Infinite-Dimensional Dynamical Systems An Introduction to Dissipative Parabolic PDEs and the Theory of Global Attractors JAMES C. ROBINSON University of Warwick HI CAMBRIDGE n S P UNIVERSITY PRESS Preface

More information

The Conjugate Gradient Method

The Conjugate Gradient Method The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large

More information

Numerical Computation of Sturm-Liouville Problem with Robin Boundary Condition

Numerical Computation of Sturm-Liouville Problem with Robin Boundary Condition Numerical Computation of Sturm-Liouville Problem with Robin Boundary Condition Theddeus T. Akano, Omotayo A. Fakinlede Abstract The modelling of physical phenomena, such as the earth s free oscillations,

More information

Introduction to PDEs and Numerical Methods: Exam 1

Introduction to PDEs and Numerical Methods: Exam 1 Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently

More information

D. Shepard, Shepard functions, late 1960s (application, surface modelling)

D. Shepard, Shepard functions, late 1960s (application, surface modelling) Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in

More information

Generalised Summation-by-Parts Operators and Variable Coefficients

Generalised Summation-by-Parts Operators and Variable Coefficients Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation

More information

T u

T u WANG LI-LIAN Assistant Professor Division of Mathematical Sciences, SPMS Nanyang Technological University Singapore, 637616 T 65-6513-7669 u 65-6316-6984 k lilian@ntu.edu.sg http://www.ntu.edu.sg/home/lilian?

More information

Partial Differential Equations

Partial Differential Equations Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This

More information

BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA

BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA 1 BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA This part of the Basic Exam covers topics at the undergraduate level, most of which might be encountered in courses here such as Math 233, 235, 425, 523, 545.

More information

Numerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by

Numerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester

More information

High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation

High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with

More information

Tyn Myint-U Lokenath Debnath. Linear Partial Differential Equations for Scientists and Engineers. Fourth Edition. Birkhauser Boston Basel Berlin

Tyn Myint-U Lokenath Debnath. Linear Partial Differential Equations for Scientists and Engineers. Fourth Edition. Birkhauser Boston Basel Berlin Tyn Myint-U Lokenath Debnath Linear Partial Differential Equations for Scientists and Engineers Fourth Edition Birkhauser Boston Basel Berlin Preface to the Fourth Edition Preface to the Third Edition

More information

THE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS

THE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS THE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS 2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. BY DR. LOTHAR COLLATZ

More information

Preliminary Examination in Numerical Analysis

Preliminary Examination in Numerical Analysis Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify

More information

arxiv: v2 [math.na] 8 Sep 2017

arxiv: v2 [math.na] 8 Sep 2017 arxiv:1704.06339v [math.na] 8 Sep 017 A Monte Carlo approach to computing stiffness matrices arising in polynomial chaos approximations Juan Galvis O. Andrés Cuervo September 3, 018 Abstract We use a Monte

More information

A MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY

A MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,

More information

Physics 6303 Lecture 8 September 25, 2017

Physics 6303 Lecture 8 September 25, 2017 Physics 6303 Lecture 8 September 25, 2017 LAST TIME: Finished tensors, vectors, and matrices At the beginning of the course, I wrote several partial differential equations (PDEs) that are used in many

More information

SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS

SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration

More information

Syllabus (Session )

Syllabus (Session ) Syllabus (Session 2016-17) Department of Mathematics nstitute of Applied Sciences & Humanities AHM-1101: ENGNEERNG MATHEMATCS Course Objective: To make the students understand the concepts of Calculus,

More information

Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS

Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general

More information

Finite Element Solver for Flux-Source Equations

Finite Element Solver for Flux-Source Equations Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University

More information

General elastic beam with an elastic foundation

General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

Numerical Methods for PDEs

Numerical Methods for PDEs Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh

More information

ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS

ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis

More information

Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions

Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New

More information

n 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes

n 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and

More information

Introduction to Partial Differential Equations

Introduction to Partial Differential Equations Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define

More information

CAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS

CAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS CAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS Preliminaries Round-off errors and computer arithmetic, algorithms and convergence Solutions of Equations in One Variable Bisection method, fixed-point

More information

Runge Kutta Chebyshev methods for parabolic problems

Runge Kutta Chebyshev methods for parabolic problems Runge Kutta Chebyshev methods for parabolic problems Xueyu Zhu Division of Appied Mathematics, Brown University December 2, 2009 Xueyu Zhu 1/18 Outline Introdution Consistency condition Stability Properties

More information

Pseudospectral and High-Order Time-Domain Forward Solvers

Pseudospectral and High-Order Time-Domain Forward Solvers Pseudospectral and High-Order Time-Domain Forward Solvers Qing H. Liu G. Zhao, T. Xiao, Y. Zeng Department of Electrical and Computer Engineering Duke University DARPA/ARO MURI Review, August 15, 2003

More information

182 George Street Providence, RI 02912, U.S.A.

182 George Street   Providence, RI 02912, U.S.A. Curriculum Vitae Guang Lin CONTACT INFORMATION Division of Applied Mathematics Phone: 401-863-3694 Box F, Brown University glin@dam.brown.edu 182 George Street http://www.dam.brown.edu/people/glin/, U.S.A.

More information

Course Description - Master in of Mathematics Comprehensive exam& Thesis Tracks

Course Description - Master in of Mathematics Comprehensive exam& Thesis Tracks Course Description - Master in of Mathematics Comprehensive exam& Thesis Tracks 1309701 Theory of ordinary differential equations Review of ODEs, existence and uniqueness of solutions for ODEs, existence

More information

Numerical Methods for Partial Differential Equations CAAM 452. Spring 2005

Numerical Methods for Partial Differential Equations CAAM 452. Spring 2005 Numerical Methods for Partial Differential Equations Instructor: Tim Warburton Class Location: Duncan Hall 1046 Class Time: 9:5am to 10:40am Office Hours: 10:45am to noon in DH 301 CAAM 45 Spring 005 Homeworks

More information

Research Article An Extension of the Legendre-Galerkin Method for Solving Sixth-Order Differential Equations with Variable Polynomial Coefficients

Research Article An Extension of the Legendre-Galerkin Method for Solving Sixth-Order Differential Equations with Variable Polynomial Coefficients Mathematical Problems in Engineering Volume 2012, Article ID 896575, 13 pages doi:10.1155/2012/896575 Research Article An Extension of the Legendre-Galerkin Method for Solving Sixth-Order Differential

More information

MA3025 Course Prerequisites

MA3025 Course Prerequisites MA3025 Course Prerequisites MA 3025 (4-1) MA3025 (4-1) Logic and Discrete Mathematics: Provides a rigorous foundation in logic and elementary discrete mathematics. Topics from logic include modeling English

More information

Numerical Analysis Comprehensive Exam Questions

Numerical Analysis Comprehensive Exam Questions Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order

More information

Scientific Computing I

Scientific Computing I Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to

More information

A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS

A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor

More information

Preliminary Examination, Numerical Analysis, August 2016

Preliminary Examination, Numerical Analysis, August 2016 Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any

More information

Kernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.

Kernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt. SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University

More information

On the Non-linear Stability of Flux Reconstruction Schemes

On the Non-linear Stability of Flux Reconstruction Schemes DOI 10.1007/s10915-011-9490-6 TECHNICA NOTE On the Non-linear Stability of Flux econstruction Schemes A. Jameson P.E. Vincent P. Castonguay eceived: 9 December 010 / evised: 17 March 011 / Accepted: 14

More information

Integration of Vlasov-type equations

Integration of Vlasov-type equations Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma

More information

AND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup

AND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)

More information

An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations

An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander

More information

Numerical Analysis and Methods for PDE I

Numerical Analysis and Methods for PDE I Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling

More information

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

Solution of Non Linear Singular Perturbation Equation. Using Hermite Collocation Method

Solution of Non Linear Singular Perturbation Equation. Using Hermite Collocation Method Applied Mathematical Sciences, Vol. 7, 03, no. 09, 5397-5408 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ams.03.37409 Solution of Non Linear Singular Perturbation Equation Using Hermite Collocation

More information

Introduction to Numerical Analysis

Introduction to Numerical Analysis J. Stoer R. Bulirsch Introduction to Numerical Analysis Second Edition Translated by R. Bartels, W. Gautschi, and C. Witzgall With 35 Illustrations Springer Contents Preface to the Second Edition Preface

More information

METHODS OF ENGINEERING MATHEMATICS

METHODS OF ENGINEERING MATHEMATICS METHODS OF ENGINEERING MATHEMATICS Edward J. Hang Kyung K. Choi Department of Mechanical Engineering College of Engineering The University of Iowa Iowa City, Iowa 52242 METHODS OF ENGINEERING MATHEMATICS

More information

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,

More information

Multi-Element Probabilistic Collocation Method in High Dimensions

Multi-Element Probabilistic Collocation Method in High Dimensions Multi-Element Probabilistic Collocation Method in High Dimensions Jasmine Foo and George Em Karniadakis Division of Applied Mathematics, Brown University, Providence, RI 02912 USA Abstract We combine multi-element

More information

Partial Differential Equations and the Finite Element Method

Partial Differential Equations and the Finite Element Method Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION

More information

Differential Equations with Boundary Value Problems

Differential Equations with Boundary Value Problems Differential Equations with Boundary Value Problems John Polking Rice University Albert Boggess Texas A&M University David Arnold College of the Redwoods Pearson Education, Inc. Upper Saddle River, New

More information

Contents. Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information

Contents. Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information Contents Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information xi xiv xvii xix 1 Preliminaries 1 1.0 Introduction.............................

More information

A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations

A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for

More information

Recurrence Relations and Fast Algorithms

Recurrence Relations and Fast Algorithms Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations

More information