Weak form of Boundary Value Problems. Simulation Methods in Acoustics
|
|
- Betty Curtis
- 5 years ago
- Views:
Transcription
1 Weak form of Boundary Value Problems Simulation Methods in Acoustics
2 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) ˆf (x) Collocation and Galerkin with test functions ψ i (x) ψ i, r = for all i = 1, 2... N Collocation: ψ i (x) = δ(x i ) Galerkin: ψi (x) = φ i (x) Other ψ φ choices are possible, but not used often in practice
3 The boundary value problem Given a BVP (PDE + BCs) A{u(x)} = f (x) x Ω (1) B i {u(x)} = ḡ i (x) x Γ i Γ i = Γ (2) A and B i are differential operators Example (Poisson equation) 2 u(x) x 2 = 1 x Ω = [, 1] (3) u(x) = x = (Dirichlet) (4) u(x) x In our example: A = 2 / x 2, f (x) 1 B1 = 1, ḡ 1 (x) B 2 = / x, ḡ 2 1 = 1 x = 1 (Neumann) (5) i
4 Weak form construction Construction steps: 1. Multiply PDE by a sufficiently well behaved 1 function ψ(x) 2. Integrate over the whole domain Ω 3. Integration by parts as much as possible (shift derivatives) 4. Apply BC 1. Multiply the PDE by the test function ψ(x) ψ(x)a{u(x)} = ψ(x)f (x) x Ω 2. If the equality holds, the integrals must be also equal ψ(x)a{u(x)} dx = ψ(x)f (x) dx ψ(x) Ω Statement: The solution of the weak form must also be the solution of the strong (original) form of the PDE Proof: The weak form must hold for all ψ(x). We can choose test functions that are nonzero only in B δ (x ) for any x Ω. Thus, the original equation must hold for all x Ω. Then, why it is a weak form? 1 The definition well behaved depends on the type of PDE and BC Ω
5 Weak form example I. Look at the Poisson equation example: 2 u(x) x 2 = f (x) x Ω = [, 1], f (x) 1 1. Multiply by the test function ψ(x) ψ(x) 2 u(x) x 2 2. The weak form is written as ψ(x) 2 u(x) x 2 dx = Ω = ψ(x)f (x) Ω ψ(x)f (x) dx 3. Apply integration by parts on the l.h.s. Recall that x 2 x 1 f g dx = [fg] x2 x 1 x 2 x 1 fg dx With f = u/ x and g = ψ(x) [ ψ(x) u(x) ] 1 1 ψ(x) u(x) + dx = x x x 1 ψ(x)f (x) dx
6 Weak form example II. 4. Apply boundary conditions u() = u(x) x = 1 x=1 ( ψ x=1 + ψ u ) + x x= 1 ψ u x x dx = 1 ψ(x)1 dx Let s compare the strong and weak forms Strong: 2nd derivative of u, Weak: 1st Strong: f continuous on Ω, Weak: integrable (f ψ integrable) The weak form imposes less strict criteria on u and f! But it requires some conditions for ψ(x). We can express these criteria using function space definitions: Strong: u C 2 (Ω), f C (Ω) Weak: u C 1 (Ω), f C 1 (Ω), u, u/ x L 2 (Ω)
7 Sobolev space definition Square-integrable functions (in Lebesgue sense) { } L 2 (Ω) = f : Ω R f 2 = f 2 dx < Sobolev space the function and its derivatives are in L p H 1 (Ω) = { f L 2 (Ω), f } x L2 (Ω) Note: derivatives are meant in a weak sense g(x) is the weak derivative of u, i.e. g = u/ x if and only if ψ(x) ψ(x)g(x) dx = u(x) dx ψ(x) C (Ω) Ω Ω x The weak derivative is unique Note: C (Ω) means infinite times differentiable functions having a compact support on Ω Ω
8 Example: solve the weak form Let s search the solution in the form: u(x) = Ax + Bx 2 ψ(x) = Âx + ˆBx 2 Find A, B such that the weak form is satisfied for all Â, ˆB Substitute: Âx + ˆBx 2 x=1 + (Âx + 2 ˆBx ) ( A + 2Bx + 1 ) ( ) (Â + 2 ˆBx A + 2Bx dx = 1 Rearrange (A and B on l.h.s.) [Â ˆB ] 1 [ ] [ ] 1 [1 ] A 2x dx = [ Â ˆB ] ( 1 2x B Evaluating the integrals [Â [ ] [ ] 1 1 A [Â ˆB] = 1 4/3 B [ ] 1/2 ˆB] 2/3 ) x= (Âx + ˆBx 2 ) dx We obtain the solution: A =, B = 1/2, thus u(x) = x 2 /2 [ ] x x 2 dx [ Â ˆB ] [ ]) 1 1
9 How did we solve? By choosing the trial functions (u(x) = Ax + Bx 2 ) the trial function space was discretized (x and x 2 are its two elements). Thus, we arrived at a finite dimensional problem (number of dimensions is d = 2 in this example) By letting ψ(x) = Âx + ˆBx 2 we used the Galerkin method. Note that we cannot apply the collocation method (ψ i = δ(x i )) as ψ/ x must be in L 2 By means of the discretization the BVP was reduced to an algebraic set of linear equations! Thus, in the end a system of linear equations was solved. In this case, we luckily got the analytical solution, however, in general we only get an approximate solution
10 Another solution Let s solve with a different right hand side: 2 u(x) x 2 = x 2 u() =, u (1) = All the same until we get [Â [ ] [ 1 1 A ˆB] 1 4/3 B By dropping [ Â ˆB ] [ ] [ ] 1 1 A = 1 4/3 B ] = [ Â ˆB ] 1 [ ] log(2)/2 π/4 1 [ ] x 1+x 2 x 2 dx [ Â ˆB ] 1+x 2 [ ] A = B [ ] The discretized weak form gives an approximate solution
11 Visualization of the example The approximation gets better if we increase the dimensionality of our function space, e.g. by adding higher order polynomials u(x) Order = 2 Weak solution u(x) -d 2 u(x)/dx 2 f = -1/(x 2 + 1) x u(x) Order = 5 Weak solution u(x) -d 2 u(x)/dx 2 f = -1/(x 2 + 1) x
12 Beam example I. Solve the following BVP (static deflection of a clamped free Euler Bernoulli beam) EIu (x) = f (x) Construct the weak form EI u() = u () = (Dirichlet) u (L) = u (L) = (Neumann) L ψ(x)u (x) dx = L Integration by parts once... ( ) L EI [ψu ] L ψ u dx = ψ(x)f (x) dx L... and twice ( ) L EI [ψu ] L [ψ u ] L + ψ u dx = f ψ dx L f ψ dx
13 Beam example II. Choose the test and trial function spaces as ( ) iπ φ i = ψ i = sin i = 1... n 2L ( ) (i n)π φ i = ψ i = cos i = n n = N 2L Note: by letting ψ i = φ i we use the Galerkin method The approximations using the above spaces reads as u(x) û(x) = ψ(x) ˆψ(x) = N q j φ j (x) = φ(x)q j=1 N w j φ j (x) = w T φ T (x) j=1 Note: φ is a row vector, q and w are column vectors
14 Beam example III. Notice that the expressions [ ψu ] L [ ψ u ] L Can all be written in a similar form φ (α) 1 w T φ(β) 1 φ (α) φ (α) φ (α) N φ(β) L 1 φ(β) N N φ(β) N ψ u dx q 1. Integration L dx or evaluation of φ(α) i j (product of shape functions and their derivatives) at a fixed point is needed. Each term gives an N N matrix in the square brackets. On the r.h.s. we have a column vector in the brackets L φ 1 (x)f (x) w T. dx = w T b φ N (x)f (x) φ (β) q N
15 Beam example IV. Write the system in matrix notation: EI [ w T (A 1 + A 2 + A 3 ) ] q = w T b Note that this is a scalar equation This must hold for all vectors w T, which is only possible if EI Aq = b Thus, we arrive at a system of linear algebraic equations
16 Beam example V. Satisfying the BCs Contrary to the previous example, our trial (shape) functions do not satisfy the Dirichlet BCs. (φ i (), φ i () ) Constraints must be imposed: [ ] φ1 () φ 2 () φ N () q 2 { φ 1 () φ 2 () φ N () =. } The same in matrix form: A c q = We already know some methods for constraints let s use the Lagrange method [ ] { } { } EI A A T c q b = A c λ This can finally be solved to get q and λ q 1 q N
17 Beam Solution The weak solution with n = 5 (N = 1) and two discrete loads as shown by the arrows in the figure.4.2 Weak solution Forces u(x) x Just to remember: the weak solution is a function ( πx û (x) =.6926 sin 2 ( πx.8171 cos 2 ) sin (πx).7489 sin ) cos (πx).5373 cos ( 3πx 2 ( 3πx 2 ) sin (2πx).217 sin ).839 cos (2πx) +.53 cos ( ) 5πx 2 ( ) 5πx 2
Finite-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 informationWeighted Residual Methods
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
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 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 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 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 informationMethod of Finite Elements I
Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 and Boundary Conditions, 26 March, 2018 Institute of Structural Engineering Method of Finite Elements I 1 Outline 1 2 Penalty method
More informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
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 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 informationWeighted 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 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 informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
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 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 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 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 informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
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 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 informationSolving the Heat Equation (Sect. 10.5).
Solving the Heat Equation Sect. 1.5. Review: The Stationary Heat Equation. The Heat Equation. The Initial-Boundary Value Problem. The separation of variables method. An example of separation of variables.
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 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 informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationEigenvalue Problem. 1 The First (Dirichlet) Eigenvalue Problem
Eigenvalue Problem A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram July 06 The method of separation variables for solving the heat equation
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 informationUndamped Vibration of a Beam
Undamped Vibration of a Beam Louie L. Yaw Walla Walla University Engineering Department PDE Class Presentation June 5, 2009 Problem - Undamped Transverse Beam Vibration +u p(x, t) 0 m(x), EI(x) +x V L
More informationDynamics of beams. Modes and waves. D. Clouteau. September 16, Department of Mechanical and Civil Engineering Ecole Centrale Paris, France
Dynamics of and waves Department of Mechanical and Civil Engineering Ecole Centrale Paris, France September 16, 2008 Outline 1 2 Pre-stressed Bernoulli-Euler in moving frames Kinematical assumption Resultant
More informationSome Aspects of Solutions of Partial Differential Equations
Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
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 informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
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 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 information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
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 informationSystems of Algebraic Equations and Systems of Differential Equations
Systems of Algebraic Equations and Systems of Differential Equations Topics: 2 by 2 systems of linear equations Matrix expression; Ax = b Solving 2 by 2 homogeneous systems Functions defined on matrices
More informationPARTIAL FRACTIONS. Introduction
Introduction PARTIAL FRACTIONS Writing any given proper rational expression of one variable as a sum (or difference) of rational expressions whose denominators are in the simplest forms is called the partial
More informationTo factor an expression means to write it as a product of factors instead of a sum of terms. The expression 3x
Factoring trinomials In general, we are factoring ax + bx + c where a, b, and c are real numbers. To factor an expression means to write it as a product of factors instead of a sum of terms. The expression
More informationAP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period:
WORKSHEET: Series, Taylor Series AP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period: 1 Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationSOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4
SOLUTIONS TO SELECTED POBLEMS FOM ASSIGNMENTS 3, 4 Problem 5 from Assignment 3 Statement. Let be an n-dimensional bounded domain with smooth boundary. Show that the eigenvalues of the Laplacian on with
More informationswapneel/207
Partial differential equations Swapneel Mahajan www.math.iitb.ac.in/ swapneel/207 1 1 Power series For a real number x 0 and a sequence (a n ) of real numbers, consider the expression a n (x x 0 ) n =
More information2.2 The Limit of a Function
2.2 The Limit of a Function Introductory Example: Consider the function f(x) = x is near 0. x f(x) x f(x) 1 3.7320508 1 4.236068 0.5 3.8708287 0.5 4.1213203 0.1 3.9748418 0.1 4.0248457 0.05 3.9874607 0.05
More informationThe geometry of least squares
The geometry of least squares We can think of a vector as a point in space, where the elements of the vector are the coordinates of the point. Consider for example, the following vector s: t = ( 4, 0),
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 informationChapter 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 informationThis Week. Professor Christopher Hoffman Math 124
This Week Sections 2.1-2.3,2.5,2.6 First homework due Tuesday night at 11:30 p.m. Average and instantaneous velocity worksheet Tuesday available at http://www.math.washington.edu/ m124/ (under week 2)
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 informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationSeparation of variables
Separation of variables Idea: Transform a PDE of 2 variables into a pair of ODEs Example : Find the general solution of u x u y = 0 Step. Assume that u(x,y) = G(x)H(y), i.e., u can be written as the product
More informationSection 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series
Section 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series Power Series for Functions We can create a Power Series (or polynomial series) that can approximate a function around
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 information2x 1 7. A linear congruence in modular arithmetic is an equation of the form. Why is the solution a set of integers rather than a unique integer?
Chapter 3: Theory of Modular Arithmetic 25 SECTION C Solving Linear Congruences By the end of this section you will be able to solve congruence equations determine the number of solutions find the multiplicative
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 informationHomework 6 Math 309 Spring 2016
Homework 6 Math 309 Spring 2016 Due May 18th Name: Solution: KEY: Do not distribute! Directions: No late homework will be accepted. The homework can be turned in during class or in the math lounge in Pedelford
More informationCHALLENGE! (0) = 5. Construct a polynomial with the following behavior at x = 0:
TAYLOR SERIES Construct a polynomial with the following behavior at x = 0: CHALLENGE! P( x) = a + ax+ ax + ax + ax 2 3 4 0 1 2 3 4 P(0) = 1 P (0) = 2 P (0) = 3 P (0) = 4 P (4) (0) = 5 Sounds hard right?
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
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 informationIntroduction to Finite Element Method. Dr. Aamer Haque
Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline Euler-Bernoulli Beams Assumptions
More informationEXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS
EXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS Course code: 8D Date: Thursday April 8 th, Time: 4h 7h Place: AUD 3 Read this first! Write your name and student identification number on each paper
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More information1. Partial differential equations. Chapter 12: Partial Differential Equations. Examples. 2. The one-dimensional wave equation
1. Partial differential equations Definitions Examples A partial differential equation PDE is an equation giving a relation between a function of two or more variables u and its partial derivatives. The
More informationQuintic 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 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 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 information1 Wave Equation on Finite Interval
1 Wave Equation on Finite Interval 1.1 Wave Equation Dirichlet Boundary Conditions u tt (x, t) = c u xx (x, t), < x < l, t > (1.1) u(, t) =, u(l, t) = u(x, ) = f(x) u t (x, ) = g(x) First we present the
More informationAnalysis III for D-BAUG, Fall 2017 Lecture 11
Analysis III for D-BAUG, Fall 2017 Lecture 11 Lecturer: Alex Sisto (sisto@math.ethz.ch) 1 Beams (Balken) Beams are basic mechanical systems that appear over and over in civil engineering. We consider a
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 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 informationFinite Fields. Saravanan Vijayakumaran Department of Electrical Engineering Indian Institute of Technology Bombay
1 / 25 Finite Fields Saravanan Vijayakumaran sarva@ee.iitb.ac.in Department of Electrical Engineering Indian Institute of Technology Bombay September 25, 2014 2 / 25 Fields Definition A set F together
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationDeflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering
Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering 008/9 Dr. Colin Caprani 1 Contents 1. Introduction... 3 1.1 General... 3 1. Background... 4 1.3 Discontinuity Functions...
More information(1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3
Math 127 Introduction and Review (1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3 MATH 127 Introduction to Calculus III
More informationTHE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t)))
THE WAVE EQUATION The aim is to derive a mathematical model that describes small vibrations of a tightly stretched flexible string for the one-dimensional case, or of a tightly stretched membrane for the
More informationDeflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering
Deflection of Flexural Members - Macaulay s Method 3rd Year Structural Engineering 009/10 Dr. Colin Caprani 1 Contents 1. Introduction... 4 1.1 General... 4 1. Background... 5 1.3 Discontinuity Functions...
More informationEquations in Quadratic Form
Equations in Quadratic Form MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: make substitutions that allow equations to be written
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
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 informationWaves and the Schroedinger Equation
Waves and the Schroedinger Equation 5 april 010 1 The Wave Equation We have seen from previous discussions that the wave-particle duality of matter requires we describe entities through some wave-form
More informationMCE503: Modeling and Simulation of Mechatronic Systems. Modeling of Continuous Beams: Finite-Mode Approach
MCE503: Modeling and Simulation of Mechatronic Systems MCE503 p./23 Modeling of Continuous Beams: Finite-Mode Approach Reading: KMR Chapter 0 Cleveland State University Mechanical Engineering Hanz Richter,
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 informationTwitter: @Owen134866 www.mathsfreeresourcelibrary.com Prior Knowledge Check 1) Simplify: a) 3x 2 5x 5 b) 5x3 y 2 15x 7 2) Factorise: a) x 2 2x 24 b) 3x 2 17x + 20 15x 2 y 3 3) Use long division to calculate:
More informationMA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation
MA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation Lecture 12 Lecture 12 MA 201, PDE (2016) 1 / 24 Formal Solution of
More informationbe the set of complex valued 2π-periodic functions f on R such that
. Fourier series. Definition.. Given a real number P, we say a complex valued function f on R is P -periodic if f(x + P ) f(x) for all x R. We let be the set of complex valued -periodic functions f on
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
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 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 information6x 3 12x 2 7x 2 +16x 7x 2 +14x 2x 4
2.3 Real Zeros of Polynomial Functions Name: Pre-calculus. Date: Block: 1. Long Division of Polynomials. We have factored polynomials of degree 2 and some specific types of polynomials of degree 3 using
More informationGeneral 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 information12.7 Heat Equation: Modeling Very Long Bars.
568 CHAP. Partial Differential Equations (PDEs).7 Heat Equation: Modeling Very Long Bars. Solution by Fourier Integrals and Transforms Our discussion of the heat equation () u t c u x in the last section
More informationn 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 informationHomework 2 - Math 603 Fall 05 Solutions
Homework 2 - Math 603 Fall 05 Solutions 1. (a): In the notation of Atiyah-Macdonald, Prop. 5.17, we have B n j=1 Av j. Since A is Noetherian, this implies that B is f.g. as an A-module. (b): By Noether
More informationCalculus 2502A - Advanced Calculus I Fall : Local minima and maxima
Calculus 50A - Advanced Calculus I Fall 014 14.7: Local minima and maxima Martin Frankland November 17, 014 In these notes, we discuss the problem of finding the local minima and maxima of a function.
More informationLecture 3: Function Spaces I Finite Elements Modeling. Bruno Lévy
Lecture 3: Function Spaces I Finite Elements Modeling Bruno Lévy Overview 1. Motivations 2. Function Spaces 3. Discretizing a PDE 4. Example: Discretizing the Laplacian 5. Eigenfunctions Spectral Mesh
More informationStructural Mechanics Column Behaviour
Structural Mechanics Column Behaviour 008/9 Dr. Colin Caprani, 1 Contents 1. Introduction... 3 1.1 Background... 3 1. Stability of Equilibrium... 4. Buckling Solutions... 6.1 Introduction... 6. Pinned-Pinned
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More information