Weighted Residual Methods
|
|
- Lawrence Todd
- 5 years ago
- Views:
Transcription
1 Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences Warsaw Poland
2 Outline 1 Problem definition Boundary-Value Problem Boundary conditions
3 Outline 1 Problem definition Boundary-Value Problem Boundary conditions 2 Weighted Residual Method General idea Approximation Error functions Minimization of errors System of algebraic equations Categories of WRM
4 Outline 1 Problem definition Boundary-Value Problem Boundary conditions 2 Weighted Residual Method General idea Approximation Error functions Minimization of errors System of algebraic equations Categories of WRM 3 ODE example A simple BVP approached by WRM Numerical solution Another numerical solution
5 Outline 1 Problem definition Boundary-Value Problem Boundary conditions 2 Weighted Residual Method General idea Approximation Error functions Minimization of errors System of algebraic equations Categories of WRM 3 ODE example A simple BVP approached by WRM Numerical solution Another numerical solution
6 Boundary-Value Problem Let B be a domain with the boundary B, and: L(.) be a (second order) differential operator, f = f (x) be a known source term in B, n = n(x) be the unit vector normal to the boundary B.
7 Boundary-Value Problem Let B be a domain with the boundary B, and: L(.) be a (second order) differential operator, f = f (x) be a known source term in B, n = n(x) be the unit vector normal to the boundary B. Boundary-Value Problem Find u = u(x) =? satisfying PDE L(u) = f in B and subject to (at least one of) the following boundary conditions u = û on B 1, u u n = ˆγ on B2, x x n + ˆα u = ˆβ on B 3, where û = û(x), ˆγ = ˆγ(x), ˆα = ˆα(x), and ˆβ = ˆβ(x) are known fields prescribed on adequate parts of the boundary B = B 1 B 2 B 3 Remarks: the boundary parts are mutually disjoint, for f 0 the PDE is called homogeneous.
8 Types of boundary conditions There are three kinds of boundary conditions: 1 the first kind or Dirichlet b.c.: u = û on B 1,
9 Types of boundary conditions There are three kinds of boundary conditions: 1 the first kind or Dirichlet b.c.: u = û on B 1, 2 the second kind or Neumann b.c.: u x n = ˆγ on B 2,
10 Types of boundary conditions There are three kinds of boundary conditions: 1 the first kind or Dirichlet b.c.: u = û on B 1, 2 the second kind or Neumann b.c.: 3 the third kind or Robin b.c.: u x n = ˆγ on B 2, u x n + ˆα u = ˆβ on B 3, also known as the generalized Neumann b.c., it can be presented as u x n = ˆγ + ˆα( û u ) on B 3. Indeed, this form is obtained for ˆβ = ˆγ + ˆα û.
11 Outline 1 Problem definition Boundary-Value Problem Boundary conditions 2 Weighted Residual Method General idea Approximation Error functions Minimization of errors System of algebraic equations Categories of WRM 3 ODE example A simple BVP approached by WRM Numerical solution Another numerical solution
12 General idea of the method Weighted Residual Method (WRM) assumes that a solution can be approximated analytically or piecewise analytically. In general, a solution to a PDE can be expressed as a linear combination of a base set of functions where the coefficients are determined by a chosen method, and the method attempts to minimize the approximation error.
13 General idea of the method Weighted Residual Method (WRM) assumes that a solution can be approximated analytically or piecewise analytically. In general, a solution to a PDE can be expressed as a linear combination of a base set of functions where the coefficients are determined by a chosen method, and the method attempts to minimize the approximation error. In fact, WRM represents a particular group of methods where an integral error is minimized in a certain way. Depending on this way the WRM can generate: the finite volume method, finite element methods, spectral methods, finite difference methods.
14 Approximation Assumption: the exact solution, u, can be approximated by a linear combination of N (linearly-independent) analytical functions, that is, N u(x) ũ(x) = U s φ s (x) s=1 Here: ũ is an approximated solution, and U s are unknown coefficients, the so-called degrees of freedom, φ s = φ s (x) form a base set of selected functions (often called as trial functions or shape functions). This set of functions generates the space of approximated solutions. s = 1,... N where N is the number of degrees of freedom.
15 Residua or error functions In general, an approximated solution, ũ, does not satisfy exactly the PDE and/or some (or all) boundary conditions. The generated errors can be described by the following error functions: 0 the PDE residuum R 0 (ũ) = L(ũ) f,
16 Residua or error functions In general, an approximated solution, ũ, does not satisfy exactly the PDE and/or some (or all) boundary conditions. The generated errors can be described by the following error functions: 0 the PDE residuum R 0 (ũ) = L(ũ) f, 1 the Dirichlet condition residuum R 1 (ũ) = ũ û,
17 Residua or error functions In general, an approximated solution, ũ, does not satisfy exactly the PDE and/or some (or all) boundary conditions. The generated errors can be described by the following error functions: 0 the PDE residuum R 0 (ũ) = L(ũ) f, 1 the Dirichlet condition residuum R 1 (ũ) = ũ û, 2 the Neumann condition residuum R 2 (ũ) = ũ x n ˆγ,
18 Residua or error functions In general, an approximated solution, ũ, does not satisfy exactly the PDE and/or some (or all) boundary conditions. The generated errors can be described by the following error functions: 0 the PDE residuum R 0 (ũ) = L(ũ) f, 1 the Dirichlet condition residuum R 1 (ũ) = ũ û, 2 the Neumann condition residuum 3 the Robin condition residuum R 2 (ũ) = ũ x n ˆγ, R 3 (ũ) = ũ x n + ˆαũ ˆβ.
19 Minimization of errors Requirement: Minimize the errors in a weighted integral sense R 0 (ũ) ψ r + R 1 (ũ) ψ r + R 2 (ũ) ψ r + R 3 (ũ) ψ r = 0. B B 1 B 2 Here, { ψ 0 } { 1 } { 2 } { 3 } r, ψ r, ψ r, and ψ r (r = 1,... M) are sets of weight functions. Note that M weight functions yield M conditions (or equations) from which to determine the N coefficients U s. To determine these N coefficients uniquely we need N independent conditions (equations). B 3
20 Minimization of errors Requirement: Minimize the errors in a weighted integral sense R 0 (ũ) ψ r + R 1 (ũ) ψ r + R 2 (ũ) ψ r + R 3 (ũ) ψ r = 0. B B 1 B 2 Here, { ψ 0 } { 1 } { 2 } { 3 } r, ψ r, ψ r, and ψ r (r = 1,... M) are sets of weight functions. Note that M weight functions yield M conditions (or equations) from which to determine the N coefficients U s. To determine these N coefficients uniquely we need N independent conditions (equations). Now, using the formulae for residua results in B L(ũ) 0 ψ r + B 1 1 ũ ψ r + B B 2 ũ x n 2 ψ r + 0 = f ψ r + B 1 B 3 1 û ψ r + B 2 B 3 ( ũ x n + ˆα ũ ) 3 ψ r 2 ˆγ ψ r + B 3 3 ˆβ ψ r.
21 System of algebraic equations By applying the approximation ũ = N U s φ s, and using the properties of linear operators, L(ũ) = N U s L(φ s ), s=1 s=1 ũ N x n = φ s U s x n, the following system of algebraic equations is obtained: N A rs U s = B r s=1 s=1 A rs = B B r = B L(φ s ) ψ 0 r + f 0 ψ r + B 1 B 1 φ s 1 û ψ r + 1 ψ r + B 2 B 2 2 ˆγ ψ r + φ s x n ψ 2 r + B 3 3 ˆβ ψ r. B 3 ( ) φs x n + ˆα φ 3 s ψ r,
22 Categories of WRM generated by weight functions There are four main categories of weight functions which generate the following categories of WRM: Subdomain method. Collocation method. Least squares method. Galerkin method.
23 Categories of WRM generated by weight functions Main categories: Subdomain method. Here the domain is divided in M subdomains B r where { 0 1 x B r, ψ r (x) = 0 outside, such that this method minimizes the residual error in each of the chosen subdomains. Note that the choice of the subdomains is free. In many cases an equal division of the total domain is likely the best choice. However, if higher resolution (and a corresponding smaller error) in a particular area is desired, a non-uniform choice may be more appropriate. Collocation method. Least squares method. Galerkin method.
24 Categories of WRM generated by weight functions Main categories: Subdomain method. Collocation method. In this method the weight functions are chosen to be Dirac delta functions 0 ψ r (x) = δ(x x r ). such that the error is zero at the chosen nodes x r. Least squares method. Galerkin method.
25 Categories of WRM generated by weight functions Main categories: Subdomain method. Collocation method. Least squares method. This method uses derivatives of the residual itself as weight functions in the form 0 ψ r (x) = R 0(ũ(x)). U r The motivation for this choice is to minimize B R2 0 of the computational domain. Note that (if the boundary conditions are satisfied) this choice of the weight function implies ( ) R 2 0 = 0 U r for all values of U r. Galerkin method. B
26 Categories of WRM generated by weight functions Main categories: Subdomain method. Collocation method. Least squares method. Galerkin method. In this method the weight functions are chosen to be identical to the base functions. 0 ψ r (x) = φ r (x). In particular, if the base function set is orthogonal (i.e., B φ r φ s = 0 if r s), this choice of weight functions implies that the residual R 0 is rendered orthogonal with the minimization condition for all base functions. B R 0 0 ψ r = 0
27 Outline 1 Problem definition Boundary-Value Problem Boundary conditions 2 Weighted Residual Method General idea Approximation Error functions Minimization of errors System of algebraic equations Categories of WRM 3 ODE example A simple BVP approached by WRM Numerical solution Another numerical solution
28 A simple BVP approached by WRM Boundary Value Problem (for an ODE) Find u = u(x) =? satisfying d 2 u dx du = 0 in B = [a, b], 2 dx subject to boundary conditions on B = B 1 B 2 = {a} {b}: u du x=a = û (Dirichlet), dx = ˆγ (Neumann). x=b
29 A simple BVP approached by WRM Boundary Value Problem (for an ODE) Find u = u(x) =? satisfying d 2 u dx du = 0 in B = [a, b], 2 dx subject to boundary conditions on B = B 1 B 2 = {a} {b}: u du x=a = û (Dirichlet), dx = ˆγ (Neumann). x=b WRM approach: Residua for an approximated solution ũ R 0(ũ) = d2 ũ dx dũ 2 dx, R1(ũ) = ũ x=a û, R 2(ũ) = dũ dx ˆγ. x=b
30 A simple BVP approached by WRM Boundary Value Problem (for an ODE) Find u = u(x) =? satisfying d 2 u dx du = 0 in B = [a, b], 2 dx subject to boundary conditions on B = B 1 B 2 = {a} {b}: u du x=a = û (Dirichlet), dx = ˆγ (Neumann). x=b WRM approach: Residua for an approximated solution ũ R 0(ũ) = d2 ũ dx dũ 2 dx, R1(ũ) = ũ x=a û, R 2(ũ) = dũ dx Minimization of weighted residual error b a ˆγ. x=b ( d2ũ dx dũ ) 0 [ (ũ ) ] [( ) ] 1 dũ ψ 2 r + û ψ r + dx x=a dx ˆγ 2 ψ r = 0. x=b
31 A simple BVP approached by WRM Boundary Value Problem (for an ODE) Find u = u(x) =? satisfying d 2 u dx du = 0 in B = [a, b], 2 dx subject to boundary conditions on B = B 1 B 2 = {a} {b}: u du x=a = û (Dirichlet), dx = ˆγ (Neumann). x=b WRM approach: System of algebraic equations: A rs = b a ( d 2 φ s dx 2 ] B r = [û ψ 1 r x=a dφs dx ) N A rs U s = B r, s=1 0 ψ r + ] + [ˆγ ψ 2 r. x=b [φ s 1 ψ r ] x=a [ dφs + dx ] 2 ψ r x=b,
32 Numerical (and exact) solution Boundary limits and values: a = 0, û = 1, b = 1, ˆγ = 2. u(x) a = b = 1 x
33 Numerical (and exact) solution Boundary limits and values: a = 0, û = 1, b = 1, ˆγ = 2. Shape functions (s = 1, 2): { φs } = { 1, e x }. 2.5 u(x) φ 2 (x) = e x φ 1 (x) = a = b = 1 x
34 Numerical (and exact) solution Boundary limits and values: a = 0, û = 1, b = 1, ˆγ = 2. Shape functions (s = 1, 2): { φs } = { 1, e x }. Weight functions (r = 1, 2): { 0 ψ r } = { 1 ψ r } = { 2 ψ r } = { 1, x }. u(x) φ 2 (x) = e x ψ 1 (x) = ψ 1 1 (x) = ψ 2 1 (x) = 1 φ 1 (x) = ψ 2 (x) = ψ 1 2 (x) = ψ 2 2 (x) = x x a = b = 1
35 Numerical (and exact) solution Boundary limits and values: a = 0, û = 1, b = 1, ˆγ = 2. Shape functions (s = 1, 2): { φs } = { 1, e x }. Weight functions (r = 1, 2): { 0 ψ r } = { 1 ψ r } = { 2 ψ r } = { 1, x }. System of equations: [ ] [ ] [ ] 1 (1 + e) U1 3 =. 0 e 2 Coefficients: Approximated solution: U 2 U 1 = 1 2 e, U2 = 2 e u(x) φ 2 (x) = e x ψ 1 (x) = ψ 1 1 (x) = ψ 2 1 (x) = 1 φ 1 (x) = 1 0 ũ(x) = 2ex +e 2 e (exact solution) ψ 2 (x) = ψ 1 2 (x) = ψ 2 2 (x) = x x a = b = 1 ũ = U 1 + U 2 e x = 2ex + e 2 e.
36 Another numerical solution Boundary limits and values: a = 0, û = 1, b = 1, ˆγ = 2. u(x) 2 u(x) = 2ex +e 2 e (exact solution) a = b = 1 x
37 Another numerical solution Boundary limits and values: a = 0, û = 1, b = 1, ˆγ = 2. Shape and weight functions (s = 1, 2, 3): { φs } = { 0 ψ s } = { 1 ψ s } = { 2 ψ s } = { 1, x, x 2 }. u(x) 2 u(x) = 2ex +e 2 e (exact solution) φ 2 (x) = x φ 1 (x) = 1 φ 3 (x) = x 2 a = b = 1 x
38 Another numerical solution Boundary limits and values: a = 0, û = 1, b = 1, ˆγ = 2. Shape and weight functions (s = 1, 2, 3): { φs } = { 0 ψ s } = { 1 ψ s } = { 2 ψ s } = { 1, x, x 2 }. System of equations: Coefficients: U 1 U 2 U 3 3 = 3. 2 U 1 = , U2 =, U3 = Approximated solution: ũ = U 1 + U 2 x + U 3 x 2 = 3 ( 5 + 4x + 4x 2 ). 17 u(x) φ 2 (x) = x ũ(x) = 3 ( x + 4x 2 ) u(x) = 2ex +e 2 e (exact solution) φ 1 (x) = 1 φ 3 (x) = x 2 a = b = 1 x
Weighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Problem definition. oundary-value Problem..................
More informationRitz Method. Introductory Course on Multiphysics Modelling
Ritz Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences Warsaw Poland
More informationAn Introductory Course on Modelling of Multiphysics Problems
An Introductory Course on Modelling of Multiphysics Problems Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationGalerkin Finite Element Model for Heat Transfer
Galerkin Finite Element Model for Heat Transfer Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents 1 Notation remarks 1 2 Local differential
More informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationWeak form of Boundary Value Problems. Simulation Methods in Acoustics
Weak form of Boundary Value Problems Simulation Methods in Acoustics Note on finite dimensional description of functions Approximation: N f (x) ˆf (x) = q j φ j (x) j=1 Residual function: r(x) = f (x)
More informationFINITE ELEMENT METHOD: APPROXIMATE SOLUTIONS
FINITE ELEMENT METHOD: APPROXIMATE SOLUTIONS I Introduction: Most engineering problems which are expressed in a differential form can only be solved in an approximate manner due to their complexity. The
More informationPART 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 informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationD(u(x)) = p(x). a i ϕ i (2.1) i=1
Chapter 2 Method of Weighted Residuals Prior to development of the Finite Element Method, there existed an approximation technique for solving differential equations called the Method of Weighted Residuals
More informationLecture 24: Starting to put it all together #3... More 2-Point Boundary value problems
Lecture 24: Starting to put it all together #3... More 2-Point Boundary value problems Outline 1) Our basic example again: -u'' + u = f(x); u(0)=α, u(l)=β 2) Solution of 2-point Boundary value problems
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
More informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationAlternative Boundary Condition Implementations for Crank Nicolson Solution to the Heat Equation
Alternative Boundary Condition Implementations for Crank Nicolson Solution to the Heat Equation ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu
More informationNumerical 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 informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationConsequences of Orthogonality
Consequences of Orthogonality Philippe B. Laval KSU Today Philippe B. Laval (KSU) Consequences of Orthogonality Today 1 / 23 Introduction The three kind of examples we did above involved Dirichlet, Neumann
More informationElements of linear algebra
Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign Read section 8. to see where equations of type (au x ) x = f show up and their (exact) solution
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in 2-D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf 2 Matthias K. Gobbert 2 1 Department of Mechanical Engineering 2 Department of Mathematics and Statistics
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationElectrodynamics PHY712. Lecture 3 Electrostatic potentials and fields. Reference: Chap. 1 in J. D. Jackson s textbook.
Electrodynamics PHY712 Lecture 3 Electrostatic potentials and fields Reference: Chap. 1 in J. D. Jackson s textbook. 1. Poisson and Laplace Equations 2. Green s Theorem 3. One-dimensional examples 1 Poisson
More informationBoundary-value Problems in Rectangular Coordinates
Boundary-value Problems in Rectangular Coordinates 2009 Outline Separation of Variables: Heat Equation on a Slab Separation of Variables: Vibrating String Separation of Variables: Laplace Equation Review
More informationA brief introduction to finite element methods
CHAPTER A brief introduction to finite element methods 1. Two-point boundary value problem and the variational formulation 1.1. The model problem. Consider the two-point boundary value problem: Given a
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationSection 12.6: Non-homogeneous Problems
Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More informationSpectral Representation of Random Processes
Spectral Representation of Random Processes Example: Represent u(t,x,q) by! u K (t, x, Q) = u k (t, x) k(q) where k(q) are orthogonal polynomials. Single Random Variable:! Let k (Q) be orthogonal with
More informationCourse and Wavelets and Filter Banks
Course 18.327 and 1.130 Wavelets and Filter Bans Numerical solution of PDEs: Galerin approximation; wavelet integrals (projection coefficients, moments and connection coefficients); convergence Numerical
More informationMEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy
MEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy e-mail: jnreddy@tamu.edu Texas A&M University College Station, TX 77843-3123 USA
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationFinite Element Method for Ordinary Differential Equations
52 Chapter 4 Finite Element Method for Ordinary Differential Equations In this chapter we consider some simple examples of the finite element method for the approximate solution of ordinary differential
More informationv(x, 0) = g(x) where g(x) = f(x) U(x). The solution is where b n = 2 g(x) sin(nπx) dx. (c) As t, we have v(x, t) 0 and u(x, t) U(x).
Problem set 4: Solutions Math 27B, Winter216 1. The following nonhomogeneous IBVP describes heat flow in a rod whose ends are held at temperatures u, u 1 : u t = u xx < x < 1, t > u(, t) = u, u(1, t) =
More informationConcepts. 3.1 Numerical Analysis. Chapter Numerical Analysis Scheme
Chapter 3 Concepts The objective of this work is to create a framework to implement multi-disciplinary finite element applications. Before starting, it is necessary to explain some basic concepts of the
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationSPECTRAL 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 informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017 Institute of Structural Engineering Method of Finite Elements II 1 Course
More informationLECTURE # 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 informationLecture 26: Putting it all together #5... Galerkin Finite Element Methods for ODE BVPs
Lecture 26: Putting it all together #5... Galerkin Finite Element Methods for ODE BVPs Outline 1) : Review 2) A specific Example: Piecewise linear elements 3) Computational Issues and Algorithms A) Calculation
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 informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
More informationMA 201: Method of Separation of Variables Finite Vibrating String Problem Lecture - 11 MA201(2016): PDE
MA 201: Method of Separation of Variables Finite Vibrating String Problem ecture - 11 IBVP for Vibrating string with no external forces We consider the problem in a computational domain (x,t) [0,] [0,
More informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods The Connection to Green s Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 1 Outline 1 Introduction
More informationTopic 12 Overview of Estimation
Topic 12 Overview of Estimation Classical Statistics 1 / 9 Outline Introduction Parameter Estimation Classical Statistics Densities and Likelihoods 2 / 9 Introduction In the simplest possible terms, the
More informationIndex. 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 informationFINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like
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 informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationTrefftz-DG solution to the Helmholtz equation involving integral equations
Trefftz-DG solution to the Helmholtz equation involving integral equations H. Barucq, A. Bendali, M. Fares, V. Mattesi, S. Tordeux Magique 3D Inria Bordeaux Sud Ouest LMA UMR CNRS 5142 INSA of Toulouse
More informationarxiv: v5 [math.na] 1 Sep 2018
High-Order Adaptive Extended Stencil FEM (AES-FEM) Part I: Convergence and Superconvergence Xiangmin Jiao Rebecca Conley Tristan J. Delaney arxiv:1603.09325v5 [math.na] 1 Sep 2018 Abstract Finite elements
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationIn this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation.
Eigen Function Expansion and Applications. In this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation. a/ The theory. b/ Example: Solving the Euler equation in two ways.
More informationModel Reduction, Centering, and the Karhunen-Loeve Expansion
Model Reduction, Centering, and the Karhunen-Loeve Expansion Sonja Glavaški, Jerrold E. Marsden, and Richard M. Murray 1 Control and Dynamical Systems, 17-81 California Institute of Technology Pasadena,
More informationSOLVING ELLIPTIC PDES
university-logo SOLVING ELLIPTIC PDES School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 POISSON S EQUATION Equation and Boundary Conditions Solving the Model Problem 3 THE LINEAR ALGEBRA PROBLEM
More informationK. BLACK To avoid these diculties, Boyd has proposed a method that proceeds by mapping a semi-innite interval to a nite interval [2]. The method is co
Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{2 c 1998 Birkhauser-Boston Spectral Element Approximations and Innite Domains Kelly Black Abstract A spectral-element
More informationVariational Principles for Equilibrium Physical Systems
Variational Principles for Equilibrium Physical Systems 1. Variational Principles One way of deriving the governing equations for a physical system is the express the relevant conservation statements and
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationMidterm 2: Sample solutions Math 118A, Fall 2013
Midterm 2: Sample solutions Math 118A, Fall 213 1. Find all separated solutions u(r,t = F(rG(t of the radially symmetric heat equation u t = k ( r u. r r r Solve for G(t explicitly. Write down an ODE for
More informationUsing Green s functions with inhomogeneous BCs
Using Green s functions with inhomogeneous BCs Using Green s functions with inhomogeneous BCs Surprise: Although Green s functions satisfy homogeneous boundary conditions, they can be used for problems
More informationBoundary regularity of solutions of degenerate elliptic equations without boundary conditions
Boundary regularity of solutions of elliptic without boundary Iowa State University November 15, 2011 1 2 3 4 If the linear second order elliptic operator L is non with smooth coefficients, then, for any
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationTwo new enriched multiscale coarse spaces for the Additive Average Schwarz method
346 Two new enriched multiscale coarse spaces for the Additive Average Schwarz method Leszek Marcinkowski 1 and Talal Rahman 2 1 Introduction We propose additive Schwarz methods with spectrally enriched
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
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 informationLecture 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 informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationLecture 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 information5. The Indirect Approach To Domain Decomposition
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c2003 DDM.org 5. The Indirect Approach To Domain Decomposition
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
More informationTopic 16 Interval Estimation
Topic 16 Interval Estimation Additional Topics 1 / 9 Outline Linear Regression Interpretation of the Confidence Interval 2 / 9 Linear Regression For ordinary linear regression, we have given least squares
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationChapter 4. Two-Dimensional Finite Element Analysis
Chapter 4. Two-Dimensional Finite Element Analysis general two-dimensional boundary-value problems 4.1 The Boundary-Value Problem 2nd-order differential equation to consider α α β φ Ω (4.1) Laplace, Poisson
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More information256 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 informationPROJECTION METHODS FOR DYNAMIC MODELS
PROJECTION METHODS FOR DYNAMIC MODELS Kenneth L. Judd Hoover Institution and NBER June 28, 2006 Functional Problems Many problems involve solving for some unknown function Dynamic programming Consumption
More informationNUMERICAL SOLUTIONS OF SYSTEM OF SECOND ORDER BOUNDARY VALUE PROBLEMS USING GALERKIN METHOD
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 37 (2017) 161-174 NUMERICAL SOLUTIONS OF SYSTEM OF SECOND ORDER BOUNDARY VALUE PROBLEMS USING GALERKIN METHOD Mahua Jahan Rupa 1 and Md. Shafiqul Islam 2
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationVirtual Work and Variational Principles
Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More information1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM)
1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM) 1.1 Introduction In this work, performances of two most widely
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More information