Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
|
|
- Suzanna Shana Dennis
- 5 years ago
- Views:
Transcription
1 Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with J. M. Melenk (TU Wien) and L. Oberbroeckling (Loyola Univ. MD)
2 The Model Problem Find U( x) u( x), v( x) such that T ( ) ( ) ( ) ( ) in (0,1) U(0) U(1) 0, E U x A x U x F x I where E, 2 0 2,0 1, 0 and A a ( x) a ( x) f( x) a21( x) a22( x) gx ( ) ( x), F( x) are given. 2
3 Assumptions: (1) The functions a ij (x), f (x),g(x) are analytic, and C, C,,,, C, f g f g a a ( n) n f C, I f fn n ( n) n g C, I g gn n ( n) n ij, I a a!, 0 0 0!, a C n! n, i, j 1,2 (2) The matrix A is pointwise positive definite, i.e. fixed > 0, such that T T 2 ξ Aξ ξ ξ ξ x I 3
4 Remark: Both components of the solution have boundary layers at x = 0 and x = 1, of width O(ln), but the second component has an additional sublayer of width O(ln). 4
5 Remark: Both components of the solution have boundary layers at x = 0 and x = 1, of width O(ln), but the second component has an additional sublayer of width O(ln). This is illustrated in the figure below, which shows the solution corresponding to A, f ( x), 10, / 2 2 4
6 5
7 Variational Formulation Find I 2 1 U H0 such that 1 U, V V V : (, ) T 2 B F u v H I 0 where, with, the usual L 2 (Ι) inner product, 2 2 UV B, u, u v, v + a u a v, u a u a v, v V,, F f u g v
8 It follows that the bilinear form is coercive, i.e. 2 1 U U U U 2 B H I, E 0 where U u v u v E 1, I 1, I 0, I 0, I denotes the energy norm. We also have the a-priori estimate U max 1, A, I 2 2 E 0, I 0, I f g 7
9 Scale separation: The relationship between and determines the number and nature of the layers. Correspondingly, there are four cases: (I) (II) (III) (IV) The no scale separation case which occurs when neither μ/1 nor ε/μ is small. The 3-scale case in which all scales are separated and occurs when μ/1 is small and ε/μ is small. The first 2-scale case which occurs when μ/1 is not small but ε/μ is small. The second 2-scale case which occurs when μ/1 is small but ε/μ is not small. 8
10 Theorem 1: There exist constants C, b, δ, q, γ > 0 independent of ε and μ, such that the following assertions are true for the solution U: 9
11 Theorem 1: There exist constants C, b, δ, q, γ > 0 independent of ε and μ, such that the following assertions are true for the solution U: ( I) U C max n, ( n ) 1/2 n 1, I n 9
12 Theorem 1: There exist constants C, b, δ, q, γ > 0 independent of ε and μ, such that the following assertions are true for the solution U: ( I) U C max n, ( n ) 1/2 n 1, I n (II) U W U Uˆ R where W C n, R R C e e, I E, I U ( n) n n b/ b / ( n) BL, I BL n ( x) C e BL ˆ ( n) n n dist( x, I )/ BL ( x) C e U n dist( x, I )/ 9
13 Theorem 1: There exist constants C, b, δ, q, γ > 0 independent of ε and μ, such that the following assertions are true for the solution U: ( I) U C max n, ( n ) 1/2 n 1, I (II) U W U Uˆ R where W C n, R R C e e, I E, I U ( n) n n b/ b / ( n) BL, I BL BL n ( x) C e ˆ ( n) n n dist( x, I )/ BL ( x) C e U n n dist( x, I )/ Additionally, the second component vˆ of sharper estimate Uˆ BL satisfies the 9
14 2 ( n) n n dist( x, I )/ vˆ ( x) C e 10
15 (III) 2 ( n) n n dist( x, I )/ vˆ ( x) C e, I E, I If / q, U W Uˆ R where R W b/ ( n) n 1, I Ce BL max n, ˆ ( n) n n dist( x, I )/ BL ( x) C e U R C n 10
16 (III) 2 ( n) n n dist( x, I )/ vˆ ( x) C e, I E, I If / q, U W Uˆ R where R W b/ ( n) n 1, I R C Ce BL max n, ˆ ( n) n n dist( x, I )/ BL ( x) C e U n Additionally, the second component vˆ of satisfies the sharper estimate Uˆ 2 ( n) n n dist( x, I )/ vˆ ( x) C e BL 10
17 (IV) U W U R R W U, I E, I ( n) n n ( n) BL, I R BL C n where C( / ) 2 b/ ( x) C / / e e n n dist ( x, I )/ 11
18 (IV) U W U R R W U, I E, I ( n) n n ( n) BL, I R BL C n where C( / ) 2 b/ ( x) C / / e e n n dist ( x, I )/ Remark: The approximation of U will be constructed according to the above regularity results, i.e. the mesh and polynomial degree distribution will be chosen appropriately, based on the relationship between ε and μ. 11
19 Discretization As usual, we seek U V 2 H 1 I 2 and we have N N 0 s. t. U, V V V 2 B F V N U U V U V N N 2 V E E N 12
20 Discretization As usual, we seek U V 2 H 1 I 2 and we have N N 0 s. t. U, V V V 2 B F V N U U V U V N N 2 V E E N The space V N is defined as follows: 12
21 Discretization As usual, we seek U V 2 H 1 I 2 and we have N N 0 s. t. U, V V V 2 B F V N N U U V U V N 2 V E E N The space V N is defined as follows: Let P n (t) denote the n th Legendre polynomial and define 1 1 1( ) 1, 2( ) i 3 1 ( ) P ( t) dt P ( ) P ( ), i 3,..., p 1 i i2 i1 i (2i 3) 12
22 Then with Π p the space of polynomials of degree p over [ 1, 1], we have,,,, p span p 1 13
23 Then with Π p the space of polynomials of degree p over [ 1, 1], we have,,,, p span p 1 Now, partition the domain I = (0, 1) by and set 0 1 x0 x1 x M I x, x, h x x, j 1,..., M j j1 j j j j1 13
24 Then with Π p the space of polynomials of degree p over [ 1, 1], we have,,,, p span p 1 Now, partition the domain I = (0, 1) by and set 0 1 x0 x1 x M I x, x, h x x, j 1,..., M j j1 j j j j1 Also define the standard (or master) element I ST = ( 1, 1), and note that it can be mapped onto the j th element by the linear mapping 13
25 1 1 x Q ( ) 1 x 1 x 2 2 j j1 j 14
26 1 1 x Q ( ) 1 x 1 x 2 2 j j1 j The finite element space V N is then defined as 1 p 0 V, u H I : u Q ( ), j 1,, M N j p j where p p p p,,, M 1 2 degrees assigned to the elements. We have dim VN, p pi 1 is the vector of polynomial M i1 14
27 An hp finite element method error estimates Definition: Spectral Boundary Layer Mesh For 0, p and 0 1 define 1 2 S, p : VN, p H0 with 0,1, p 1/ 2 0, p, p,1 p,1 p,1, p 1/ 2 0, p,1 p,1, p 1/ 2 p The polynomial degree is taken to be uniformly p over all elements. 15
28 In practice, the mesh is constructed as follows: If both ε and μ are large (and no layers are present), then [0,1], p If both ε and μ are small (i.e. 0 < ε < μ << 1and both layers are present), then 0, p, p,1 p,1 p,1, p In all other cases, 0, p,1 p,1, p 16
29 Theorem 2: There exist constants C,, > 0, depending only on the input data, such that U U N Ce p E where U is the exact solution and U N is the finite element solution computed using the Spectral Boundary Layer Mesh. Sketch of Proof: We utilize the previously stated regularity results, along with the following approximation results, the first one being used for the approximation of the smooth part(s) and the boundary layers (within the layer). 17
30 Lemma 1: Let I j be an interval of length h j and let V C (I ST ) satisfy for some C u, γ u > 0, K 1, V ( n), I j n C max n, K n 1,2,3... u u Then there exist η, β, C > 0, depending only on γ u such hk j that, under the condition, we have p 1/2 hj K h V V V V CC e p 1 j p j p j u 0, I j 0, I where p j j : H 1 (I j ) Π p is a linear operator that satisfies V( I ) V( I ). p j j j n j j p j, 18
31 The next result is used for the approximation of the boundary layers outside the layer. Lemma 2: Let ν > 0 and let u satisfy u x u x C e x I dist( x, I )/ ( ) ( ) u. Let Δ be an arbitrary mesh on I with mesh points ξ and 1 ξ where ξ (0, ½). Then the piecewise linear interpolant π 1 u Satisfies on (ξ, 1 ξ ): / u 1u u 1u CC 0,(,1 ) ue 0,(,1 ) for some C > 0 independent of ν., 19
32 The proof is separated in four cases corresponding to the four cases stated in the regularity results. In Case I (asymptotic case), Theorem 1 and Lemma 1 give the desired result. In Case II (3 scale separation), Theorem 1 and Lemma 1 allow us to handle the smooth part and the remainder in the expansion. For the layers, we use Lemma 2 for their approximation outside the layer region and Lemma 1 within. Cases III and IV (2 scale separation) follow with similar combinations of Theorem 1 and Lemmas 1 and 2. 20
33 Numerical Results We consider the problem with A 2 1 1, 1 2 F 1 We are computing Error 100 U U EXACT FEM EI, U EXACT EI, (An exact solution is available.) 21
34 We will be comparing the following methods: The h version on a uniform mesh, with p = 1, 2, 3 The p version on a single element, with p = 1, 2, The hp version on the 5 element (variable) mesh, with p = 1, 2, The h version on a Shishkin mesh, with p = 1, 2, 3 The h version on an exponentially graded mesh, with p = 1, 2, 3, where the mesh points are given by 22
35 with M x j x j x j j0 1 sj x j (2 p 1)ln 1, j 0,, M 2 M 2 s 1 exp (2 p 1)
36 Percentage Relative Error in the Energy Norm = 0.4, = 1 h-ver, unif mesh, p = 1 h-ver, unif mesh, p = 2 h-ver, unif mesh, p = 3 p-ver, 1 elem slope -1 slope slope DOF 24
37 Percentage Relative Error in the Energy Norm = 0.01, = 0.1 h-ver, unif mesh, p = 1 h-ver, unif mesh, p = 2 h-ver, unif mesh, p = 3 p-ver, 1 elem slope -0.9 slope slope DOF 25
38 Percentage Relative Error in the Energy Norm = 10-7/ , = 0.01 slope -1.5 h-ver, unif mesh, p = 1 h-ver, unif mesh, p = 2 h-ver, unif mesh, p = 3 p-ver, 1 elem slope DOF 26
39 10 2 = 10-7/ , = 0.01 Percentage Relative Error in the Energy Norm h-ver, unif mesh, p = 1 p-ver, 1 elem hp-ver, 5 elem h-ver, exp mesh, p = 1 h-ver, Shishkin mesh, p = DOF slope -1 slope -0.7 slope
40 10 1 = 10-7/ , = 0.01 Percentage Relative Error in the Energy Norm h-ver, unif mesh, p = 1 h-ver, unif mesh, p = 2 h-ver, unif mesh, p = 3 h-ver, Shishkin mesh, p = 1 h-ver, Shishkin mesh, p = 2 h-ver, Shishkin mesh, p = DOF slope -1 slope -3 slope -2 slope slope slope
41 10 1 = 10-7/ , = 0.01 Percentage Relative Error in the Energy Norm h-ver, unif mesh, p = 1 h-ver, unif mesh, p = 2 h-ver, unif mesh, p = 3 h-ver, Shishkin mesh, p = 1 h-ver, Shishkin mesh, p = 2 h-ver, Shishkin mesh, p = 3 hp-ver, 5 elem DOF 29
42 10 2 = 10-5, = 10-3 Percentage Relative Error in the Energy Norm h-ver, unif mesh, p = 1 p-ver, 1 elem hp-ver, 5 elem h-ver, exp mesh, p = 1 h-ver, Shishkin mesh, p = DOF slope -1 slope slope
43 hp-version, 5 elements Percentage Relative Error in the Energy Norm , = 0.01 = 10-5, = 10-3 = 10-6, = DOF 31
44 Example 2: A cos x/ 4 2.2e 2 2 x x 1 x, x 2e 0 f, u(0) u(1) 10x 1 0 (An exact solution is NOT available, hence we use a reference solution.) 32
45 Percentage Relative Error in the Energy Norm = 10-7/ , = 0.01 hp ver, 5 elem h ver, Shishkin mesh, p=1 h ver, exp mesh, p=1 slope -1 slope DOF 33
46 Percentage Relative Error in the Energy Norm 10 2 = 10-5, = hp ver, 5 elem h ver, Shishkin mesh, p=1 h ver, exp mesh, p=1 slope -1 slope DOF 34
47 Percentage Relative Error in the Energy Norm hp-version, 5 elements , = 0.01 = 10-5, = 10-3 = 10-6, = DOF
48 Closing Remarks The hp FEM on the Spectral Boundary Layer Μesh yields robust exponential convergence in the energy norm, for the entire range 0 ε μ 1. The regularity theory developed was crucial in the proof of the main approximation result. Extending the current approach to systems with more equations, while conceptually straight forward, appears to be cumbersome. Systems of two equations of convection-diffusion type are currently being investigated. 36
An exponentially graded mesh for singularly perturbed problems
An exponentially graded mesh for singularly perturbed problems Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with P. Constantinou (UCY), S. Franz (TU Dresden)
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationParameter robust methods for second-order complex-valued reaction-diffusion equations
Postgraduate Modelling Research Group, 10 November 2017 Parameter robust methods for second-order complex-valued reaction-diffusion equations Faiza Alssaedi Supervisor: Niall Madden School of Mathematics,
More informationA PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM
A PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN Abstract. In this paper a singularly perturbed reaction diffusion partial
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationNumerical Solutions of Laplacian Problems over L-Shaped Domains and Calculations of the Generalized Stress Intensity Factors
WCCM V Fifth World Congress on Computational Mechanics July 7-2, 2002, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Numerical Solutions of Laplacian Problems over L-Shaped Domains
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationA first-order system Petrov-Galerkin discretisation for a reaction-diffusion problem on a fitted mesh
A first-order system Petrov-Galerkin discretisation for a reaction-diffusion problem on a fitted mesh James Adler, Department of Mathematics Tufts University Medford, MA 02155 Scott MacLachlan Department
More informationFinite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems
Finite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems Zhimin Zhang Department of Mathematics Wayne State University Detroit, MI 48202 Received 19 July 2000; accepted
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
4 Interpolation 4.1 Polynomial interpolation Problem: LetP n (I), n ln, I := [a,b] lr, be the linear space of polynomials of degree n on I, P n (I) := { p n : I lr p n (x) = n i=0 a i x i, a i lr, 0 i
More informationConvergence on Layer-adapted Meshes and Anisotropic Interpolation Error Estimates of Non-Standard Higher Order Finite Elements
Convergence on Layer-adapted Meshes and Anisotropic Interpolation Error Estimates of Non-Standard Higher Order Finite Elements Sebastian Franz a,1, Gunar Matthies b, a Department of Mathematics and Statistics,
More informationStochastic 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 informationSuperconvergence analysis of the SDFEM for elliptic problems with characteristic layers
Superconvergence analysis of the SDFEM for elliptic problems with characteristic layers S. Franz, T. Linß, H.-G. Roos Institut für Numerische Mathematik, Technische Universität Dresden, D-01062, Germany
More informationA uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction diffusion problems
IMA Journal of Numerical Analysis (2003 23, 627 644 A uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction diffusion problems NIALL MADDEN Department of
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationApproximation of Geometric Data
Supervised by: Philipp Grohs, ETH Zürich August 19, 2013 Outline 1 Motivation Outline 1 Motivation 2 Outline 1 Motivation 2 3 Goal: Solving PDE s an optimization problems where we seek a function with
More informationCholesky factorisations of linear systems coming from a finite difference method applied to singularly perturbed problems
Cholesky factorisations of linear systems coming from a finite difference method applied to singularly perturbed problems Thái Anh Nhan and Niall Madden The Boundary and Interior Layers - Computational
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationEfficient hp-finite elements
Efficient hp-finite elements Ammon Washburn July 31, 2015 Abstract Ways to make an hp-finite element method efficient are presented. Standard FEMs and hp-fems with their various strengths and weaknesses
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationEfficient Reduced Order Modeling of Low- to Mid-Frequency Vibration and Power Flow in Complex Structures
Efficient Reduced Order Modeling of Low- to Mid-Frequency Vibration and Power Flow in Complex Structures Yung-Chang Tan Graduate Student Research Assistant Matthew P. Castanier Assistant Research Scientist
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationIterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations
Iterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations P.A. Farrell 1, P.W. Hemker 2, G.I. Shishkin 3 and L.P. Shishkina 3 1 Department of Computer Science,
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationhp-dg timestepping for time-singular parabolic PDEs and VIs
hp-dg timestepping for time-singular parabolic PDEs and VIs Ch. Schwab ( SAM, ETH Zürich ) (joint with O. Reichmann) ACMAC / Heraklion, Crete September 26-28, 2011 September 27, 2011 Motivation Well-posedness
More informationQuintic B-spline method for numerical solution of fourth order singular perturbation boundary value problems
Stud. Univ. Babeş-Bolyai Math. 63(2018, No. 1, 141 151 DOI: 10.24193/subbmath.2018.1.09 Quintic B-spline method for numerical solution of fourth order singular perturbation boundary value problems Ram
More informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationMultiscale methods for time-harmonic acoustic and elastic wave propagation
Multiscale methods for time-harmonic acoustic and elastic wave propagation Dietmar Gallistl (joint work with D. Brown and D. Peterseim) Institut für Angewandte und Numerische Mathematik Karlsruher Institut
More informationEfficiency-based h- and hp-refinement strategies for finite element methods
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 7; : 5 [Version: /9/8 v.] Efficiency-based h- and hp-refinement strategies for finite element methods H. De Sterck, Thomas A. Manteuffel,
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationImplementation 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 informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationAN OPTIMAL UNIFORM A PRIORI ERROR ESTIMATE FOR AN UNSTEADY SINGULARLY PERTURBED PROBLEM
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 11, Number 1, Pages 24 33 c 2014 Institute for Scientific Computing and Information AN OPTIMAL UNIFORM A PRIORI ERROR ESTIMATE FOR AN UNSTEADY
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationInterior Layers in Singularly Perturbed Problems
Interior Layers in Singularly Perturbed Problems Eugene O Riordan Abstract To construct layer adapted meshes for a class of singularly perturbed problems, whose solutions contain boundary layers, it is
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationCondition number estimates for matrices arising in the isogeometric discretizations
www.oeaw.ac.at Condition number estimates for matrices arising in the isogeometric discretizations K. Gahalaut, S. Tomar RICAM-Report -3 www.ricam.oeaw.ac.at Condition number estimates for matrices arising
More informationAnalysis of a Streamline-Diffusion Finite Element Method on Bakhvalov-Shishkin Mesh for Singularly Perturbed Problem
Numer. Math. Theor. Meth. Appl. Vol. 10, No. 1, pp. 44-64 doi: 10.408/nmtma.017.y1306 February 017 Analysis of a Streamline-Diffusion Finite Element Method on Bakhvalov-Shishkin Mesh for Singularly Perturbed
More informationNEW A PRIORI FEM ERROR ESTIMATES FOR EIGENVALUES
NEW A PRIORI FEM ERROR ESTIMATES FOR EIGENVALUES ANDREW V. KNYAZEV AND JOHN E. OSBORN Abstract. We analyze the Ritz method for symmetric eigenvalue problems and prove a priori eigenvalue error estimates.
More informationApplied Mathematics and Computation 169 (2005)
Applied Mathematics and Computation 169 (2005) 485 499 www.elsevier.com/locate/amc Solving Laplacian problems with boundary singularities: a comparison of a singular function boundary integral method with
More informationHIGH ORDER GALERKIN METHODS WITH GRADED MESHES FOR TWO-DIMENSIONAL REACTION-DIFFUSION PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 3, Pages 319 343 c 2016 Institute for Scientific Computing and Information HIGH ORDER GALERKIN METHODS WITH GRADED MESHES FOR
More informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationON THE ENERGY DECAY OF TWO COUPLED STRINGS THROUGH A JOINT DAMPER
Journal of Sound and Vibration (997) 203(3), 447 455 ON THE ENERGY DECAY OF TWO COUPLED STRINGS THROUGH A JOINT DAMPER Department of Mechanical and Automation Engineering, The Chinese University of Hong
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationQuarkonial frames of wavelet type - Stability, approximation and compression properties
Quarkonial frames of wavelet type - Stability, approximation and compression properties Stephan Dahlke 1 Peter Oswald 2 Thorsten Raasch 3 ESI Workshop Wavelet methods in scientific computing Vienna, November
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
More informationINVESTIGATION OF STABILITY AND ACCURACY OF HIGH ORDER GENERALIZED FINITE ELEMENT METHODS HAOYANG LI THESIS
c 2014 Haoyang Li INVESTIGATION OF STABILITY AND ACCURACY OF HIGH ORDER GENERALIZED FINITE ELEMENT METHODS BY HAOYANG LI THESIS Submitted in partial fulfillment of the requirements for the degree of Master
More informationThe hp Finite Element Method For Singularly Perturbed Problems in Nonsmooth Domains
The hp Finite Element Method For Singularly Perturbed Problems in Nonsmooth Domains Christos Xenophontos Department of Mathematics and Computer Science Clarkson University Potsdam, New York 13699-5815
More informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More informationThe Fourier series for a 2π-periodic function
The Fourier series for a 2π-periodic function Let f : ( π, π] R be a bounded piecewise continuous function which we continue to be a 2π-periodic function defined on R, i.e. f (x + 2π) = f (x), x R. The
More informationInterpolation. Chapter Interpolation. 7.2 Existence, Uniqueness and conditioning
76 Chapter 7 Interpolation 7.1 Interpolation Definition 7.1.1. Interpolation of a given function f defined on an interval [a,b] by a polynomial p: Given a set of specified points {(t i,y i } n with {t
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationLie Groups for 2D and 3D Transformations
Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and
More informationScientific 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 informationLinear maps. Matthew Macauley. Department of Mathematical Sciences Clemson University Math 8530, Spring 2017
Linear maps Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 8530, Spring 2017 M. Macauley (Clemson) Linear maps Math 8530, Spring 2017
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationFree probabilities and the large N limit, IV. Loop equations and all-order asymptotic expansions... Gaëtan Borot
Free probabilities and the large N limit, IV March 27th 2014 Loop equations and all-order asymptotic expansions... Gaëtan Borot MPIM Bonn & MIT based on joint works with Alice Guionnet, MIT Karol Kozlowski,
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationAn inverse elliptic problem of medical optics with experimental data
An inverse elliptic problem of medical optics with experimental data Jianzhong Su, Michael V. Klibanov, Yueming Liu, Zhijin Lin Natee Pantong and Hanli Liu Abstract A numerical method for an inverse problem
More informationACCURATE NUMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRANSFER
FDS-2000 Conference, pp. 1 10 M. Sapagovas et al. (Eds) c 2000 MII ACCURATE UMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRASFER B. GAHA 1, J. J. H. MILLER 2, and G. I. SHISHKI
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationis the desired collar neighbourhood. Corollary Suppose M1 n, M 2 and f : N1 n 1 N2 n 1 is a diffeomorphism between some connected components
1. Collar neighbourhood theorem Definition 1.0.1. Let M n be a manifold with boundary. Let p M. A vector v T p M is called inward if for some local chart x: U V where U subsetm is open, V H n is open and
More informationOptimization and Optimal Control in Banach Spaces
Optimization and Optimal Control in Banach Spaces Bernhard Schmitzer October 19, 2017 1 Convex non-smooth optimization with proximal operators Remark 1.1 (Motivation). Convex optimization: easier to solve,
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationParameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition
Parameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition P. A. Farrell 1,A.F.Hegarty 2, J. J. H. Miller 3, E. O Riordan 4,andG.I. Shishkin 5 1
More informationGALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS
GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More information, applyingl Hospital s Rule again x 0 2 cos(x) xsinx
Lecture 3 We give a couple examples of using L Hospital s Rule: Example 3.. [ (a) Compute x 0 sin(x) x. To put this into a form for L Hospital s Rule we first put it over a common denominator [ x 0 sin(x)
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationQuestion 9: PDEs Given the function f(x, y), consider the problem: = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1. x 2 u. u(x, 0) = u(x, 1) = 0 for 0 x 1
Question 9: PDEs Given the function f(x, y), consider the problem: 2 u x 2 u = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1 u(x, 0) = u(x, 1) = 0 for 0 x 1 u(0, y) = u(1, y) = 0 for 0 y 1. a. Discuss how you
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
More informationKey words. Incompressible magnetohydrodynamics, mixed finite element methods, discontinuous Galerkin methods
A MIXED DG METHOD FOR LINEARIZED INCOMPRESSIBLE MAGNETOHYDRODYNAMICS PAUL HOUSTON, DOMINIK SCHÖTZAU, AND XIAOXI WEI Journal of Scientific Computing, vol. 40, pp. 8 34, 009 Abstract. We introduce and analyze
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationSpectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem
Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v
More informationThe Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge
The Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge Vladimir Kozlov (Linköping University, Sweden) 2010 joint work with A.Nazarov Lu t u a ij
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationDepartment of Aerospace Engineering AE602 Mathematics for Aerospace Engineers Assignment No. 6
Department of Aerospace Engineering AE Mathematics for Aerospace Engineers Assignment No.. Find the best least squares solution x to x, x 5. What error E is minimized? heck that the error vector ( x, 5
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
Bath, October 28th 2005 1 High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde Steve Arden Funded by: Leverhulme Trust University
More informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More information1 Orthonormal sets in Hilbert space
Math 857 Fall 15 1 Orthonormal sets in Hilbert space Let S H. We denote by [S] the span of S, i.e., the set of all linear combinations of elements from S. A set {u α : α A} is called orthonormal, if u
More information