Alternative Boundary Condition Implementations for Crank Nicolson Solution to the Heat Equation
|
|
- Augusta Horn
- 5 years ago
- Views:
Transcription
1 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 ME 448/548: Alternative BC Implementation for the Heat Equation
2 Overview. Goal is to allow Dirichlet, Neumann and mixed boundary conditions 2. Use ghost node formulation Preserve spatial accuracy of O( x 2 ) Preserve tridiagonal structure to the coefficient matrix 3. Implement in a code that uses the Crank-Nicolson scheme. 4. Demonstrate the technique on sample problems ME 448/548: Alternative BC Implementation for the Heat Equation page
3 Heat Transfer Boundary Conditions T w i = x = 2 3 Constant wall temperature, T w, is prescribed. The solution gives the internal temperatures and the wall heat flux, q w. q w i = x = 2 3 Constant wall heat flux, q w, is prescribed. The solution gives the internal temperatures temperatures and wall temperature, T w h, T i = x = 2 3 Convective boundary condition where h and T are prescribed. The solution gives the internal temperatures, the wall temperature, T w, and the wall heat flux, q w. ME 448/548: Alternative BC Implementation for the Heat Equation page 2
4 Convective Boundary Condition The general form of a convective boundary condition is u x = g + h u () x= This is also known as a Robin boundary condition or a boundary condition of the third kind. The simplistic implementation is to replace the derivative in Equation () with a one-sided difference u k+ 2 u k+ = g + h u k+ (2) x
5 Introduce a Ghost Node Imagine that there is a node û that is outside of the domain u u u 2 u 3... x x x x x 2 this node is used to enforce the boundary condition from Equation (). The value û does not explicitly appear in the numerical scheme. ME 448/548: Alternative BC Implementation for the Heat Equation page 4
6 Use the BC to compute û by extrapolation Use a central difference approximation at x = (x = x ) to impose the boundary condition. The value of û condition is u 2 û 2 x = g + h u. (3) consistent with the boundary Edge of physical domain u u 2... û = u 2 2 x(g + h u ). (4) x x x x x 2 ME 448/548: Alternative BC Implementation for the Heat Equation page 5
7 Equation for u Evaluate the finite difference form of the heat equation at x = x. u k+ u k t = θα 2u k+ + u k+ ] [ûk+ x 2 2 [ûk + ( θ)α 2uk + uk 2 x 2 ] Choose θ = /2 and use the formulas for û at time step k and time step k + u k+ u k t = θα x 2 [ u k+ 2 2 x(g k+ + h k+ u k+ ) 2u k+ + u k+ 2 ( θ)α + x 2 The terms in boxes are from the boundary condition [ u k 2 2 x(gk + hk uk ) 2uk + uk 2 ] ] ME 448/548: Alternative BC Implementation for the Heat Equation page 6
8 Rearrange the Equation for u Algebraically rearranging the preceding equation gives where a u k+ + b u k+ 2 = d (5) a = t + 2θα x2( + xhk+ ) (6) b = 2θα (7) x [ 2 ] 2( θ)α d = ( + xh k t x 2 ) u k (8) 2( θ)α + u k x 2 2 2α [ ] θg k+ + ( θ)g k x These equations define the terms for the first row in the system of equations ME 448/548: Alternative BC Implementation for the Heat Equation page 7
9 Data structure for implementing alternative BC in the Matlab code Store the boundary condition specification in a 2 3 matrix. The first row has data for x = and the second row has data for x = L. The first column is a flag with the boundary condition type. ubc = Type Value Value 2 Type Value Value 2 b = at the x = boundary and b = 2 at the x = L boundary ubc(b, ) = : u(x b ) = value ubc(b, 2) = value of u at boundary ubc(b, 3) = not used ubc(b, ) = 2: u/ x xb = g + hu(x b ) ubc(b, 2) = g ubc(b, 3) = h ME 448/548: Alternative BC Implementation for the Heat Equation page 8
10 Verification: Solve the toy problem on half of the domain The toy problem used to test the codes u t = α 2 u x 2 t >, x L u(, t) = u(l, t) = ; u(x, ) = sin(πx/l) only needs to be solved on one half of the domain u t = u α 2 t >, x L/2 x 2 u u(, t) = ; x = L/2 u(x, ) = sin(πx/l) u t = u α 2 t >, L/2 x L x 2 u x = u(l, t) = ; L/2 u(x, ) = sin(πx/l) ME 448/548: Alternative BC Implementation for the Heat Equation page 9
11 Verification: Solve the toy problem on half of the domain Output of democnbc.5 IC CN Exact.5 IC CN Exact u x x ME 448/548: Alternative BC Implementation for the Heat Equation page
BTCS Solution to the Heat Equation
BTCS Solution to the Heat Equation ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@mepdxedu ME 448/548: BTCS Solution to the Heat Equation Overview
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Models for CFD ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: 1D Convection-Diffusion
More informationFinite Difference Methods (FDMs) 2
Finite Difference Methods (FDMs) 2 Time- dependent PDEs A partial differential equation of the form (15.1) where A, B, and C are constants, is called quasilinear. There are three types of quasilinear equations:
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationIntroduction to Finite Di erence Methods
Introduction to Finite Di erence Methods ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: Introduction to Finite Di erence Approximation
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 informationMTF071 Computational Fluid Dynamics of Turbulent
CHALMERS TEKNISKA HÖGSKOLA Termo- och Fluiddynamik 42 96 Göteborg MTF07 Computational Fluid Dynamics of Turbulent Flow http://wwwtfdchalmersse/gr-kurs/mtf07 Task K2 Lars Davidson 2003-02-03 In Task K you
More informationOverview of Convection Heat Transfer
Overview of Convection Heat Transfer ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu February 23, 2012 ME 448/548: Convection Heat Transfer
More informationFinite Difference Methods (FDMs) 1
Finite Difference Methods (FDMs) 1 1 st - order Approxima9on Recall Taylor series expansion: Forward difference: Backward difference: Central difference: 2 nd - order Approxima9on Forward difference: Backward
More informationOverview of Convection Heat Transfer
Overview of Convection Heat Transfer ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu February 23, 2012 ME 448/548: Convection Heat Transfer
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 information1 Finite difference example: 1D implicit heat equation
1 Finite difference example: 1D implicit heat equation 1.1 Boundary conditions Neumann and Dirichlet We solve the transient heat equation ρc p t = ( k ) (1) on the domain L/2 x L/2 subject to the following
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
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 informationThe Control-Volume Finite-Difference Approximation to the Diffusion Equation
The Control-Volume Finite-Difference Approimation to the Diffusion Equation ME 448/548 Notes Gerald Recktenwald Portland State Universit Department of Mechanical Engineering gerr@mepdedu ME 448/548: D
More informationInitial and Boundary Conditions
Initial and Boundary Conditions Initial- and boundary conditions are needed For a steady problems correct initial conditions is important to reduce computational time and reach convergence Boundary conditions
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
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 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 informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationBoundary conditions. Diffusion 2: Boundary conditions, long time behavior
Boundary conditions In a domain Ω one has to add boundary conditions to the heat (or diffusion) equation: 1. u(x, t) = φ for x Ω. Temperature given at the boundary. Also density given at the boundary.
More informationSalinity Calibration fit with Matlab
Salinity Calibration fit with Matlab ME 121 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 121: Salinity calibration fit Overview These slides
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
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 informationModule 7: The Laplace Equation
Module 7: The Laplace Equation In this module, we shall study one of the most important partial differential equations in physics known as the Laplace equation 2 u = 0 in Ω R n, (1) where 2 u := n i=1
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationA Review of Linear Algebra
A Review of Linear Algebra Gerald Recktenwald Portland State University Mechanical Engineering Department gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab: Implementations
More informationMATH 333: Partial Differential Equations
MATH 333: Partial Differential Equations Problem Set 9, Final version Due Date: Tues., Nov. 29, 2011 Relevant sources: Farlow s book: Lessons 9, 37 39 MacCluer s book: Chapter 3 44 Show that the Poisson
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationarxiv: v2 [math.na] 11 Jun 2018
DISCRETIZATIONS OF THE SPECTRAL FRACTIONAL LAPLACIAN ON GENERAL DOMAINS WITH DIRICHLET, NEUMANN, AND ROIN OUNDARY CONDITIONS NICOLE CUSIMANO, FÉLIX DEL TESO, LUCA GERARDO-GIORDA, AND GIANNI PAGNINI arxiv:78.362v2
More informationHigh-order finite volume method for curved boundaries and non-matching domain-mesh problems
High-order finite volume method for curved boundaries and non-matching domain-mesh problems May 3-7, 0, São Félix, Portugal Ricardo Costa,, Stéphane Clain,, Gaspar J. Machado, Raphaël Loubère Institut
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
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 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 informationAn Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations.
An Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations by Tong Chen A thesis submitted in conformity with the requirements
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
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 informationDirect Numerical Simulation of Two-Phase Flows with phase change
Direct Numerical Simulation of Two-Phase Flows with phase change Sébastien Tanguy, Lucia Rueda Villegas, Catherine Colin, Grégory Huber Institut de Mécanique des Fluides de Toulouse Maître de conférence
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationChapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
More informationApproximations of diffusions. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 3b Approximations of diffusions Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.1 04/10/2018 1 Learning objectives Become aware of the existence of stability conditions for the
More informationMethod of Lines. Received April 20, 2009; accepted July 9, 2009
Method of Lines Samir Hamdi, William E. Schiesser and Graham W. Griffiths * Ecole Polytechnique, France; Lehigh University, USA; City University,UK. Received April 20, 2009; accepted July 9, 2009 The method
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
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 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 informationPartial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction
Page of 8 Partial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction FILE:Chap 3 Partial Differential Equations-V6. Original: May 7, 05 Revised: Dec 9, 06, Feb 0, 07,
More informationA CCD-ADI method for unsteady convection-diffusion equations
A CCD-ADI method for unsteady convection-diffusion equations Hai-Wei Sun, Leonard Z. Li Department of Mathematics, University of Macau, Macao Abstract With a combined compact difference scheme for the
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 informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More information1 Goal: Use Mass Balances to Predict Control Corrections
ME 121 February 21, 2018 Mass Additions to Control Fish Tank Salinity Gerald Recktenwald gerry@pdx.edu 1 Goal: Use Mass Balances to Predict Control Corrections A mass balance is at the core of a model
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 informationNumerical Solution of f(x) = 0
Numerical Solution of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@pdx.edu ME 350: Finding roots of f(x) = 0 Overview Topics covered in these slides
More informationOn the Discretization Time-Step in the Finite Element Theta-Method of the Discrete Heat Equation
On the Discretization Time-Step in the Finite Element Theta-Method of the Discrete Heat Equation Tamás Szabó Eötvös Loránd University, Institute of Mathematics 1117 Budapest, Pázmány P. S. 1/c, Hungary
More informationCode: 101MAT4 101MT4B. Today s topics Finite-difference method in 2D. Poisson equation Wave equation
Code: MAT MTB Today s topics Finite-difference method in D Poisson equation Wave equation Finite-difference method for elliptic PDEs in D Recall that u is the short version of u x + u y Dirichlet BVP:
More informationQuaternion Cubic Spline
Quaternion Cubic Spline James McEnnan jmcennan@mailaps.org May 28, 23 1. INTRODUCTION A quaternion spline is an interpolation which matches quaternion values at specified times such that the quaternion
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 informationMat1062: Introductory Numerical Methods for PDE
Mat1062: Introductory Numerical Methods for PDE Mary Pugh January 13, 2009 1 Ownership These notes are the joint property of Rob Almgren and Mary Pugh 2 Boundary Conditions We now want to discuss in detail
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
More informationNumerical Solutions to PDE s
Introduction Numerical Solutions to PDE s Mathematical Modelling Week 5 Kurt Bryan Let s start by recalling a simple numerical scheme for solving ODE s. Suppose we have an ODE u (t) = f(t, u(t)) for some
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationGREEN S IDENTITIES AND GREEN S FUNCTIONS
GREEN S IENTITIES AN GREEN S FUNCTIONS Green s first identity First, recall the following theorem. Theorem: (ivergence Theorem) Let be a bounded solid region with a piecewise C 1 boundary surface. Let
More informationLECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES
LECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES Figure 3.47 a. Two-mass, linear vibration system with spring connections. b. Free-body diagrams. c. Alternative free-body
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II Supplementary Material on Finite Different Methods for Two-Point Boundary-Value Problems Xiangmin Jiao Stony Brook University Xiangmin Jiao Stony Brook University AMS527:
More informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
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 informationPartial Differential Equations
Next: Using Matlab Up: Numerical Analysis for Chemical Previous: Ordinary Differential Equations Subsections Finite Difference: Elliptic Equations The Laplace Equations Solution Techniques Boundary Conditions
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationOverview of Convection Heat Transfer
Overview of Convection Heat Transfer Lecture Notes for ME 448/548 Gerald Recktenwald March 16, 2006 1 Introduction CFD is a powerful tool for engineering design involving fluid flow and heat transfer.
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter Boundary-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
More informationApplications of the Maximum Principle
Jim Lambers MAT 606 Spring Semester 2015-16 Lecture 26 Notes These notes correspond to Sections 7.4-7.6 in the text. Applications of the Maximum Principle The maximum principle for Laplace s equation is
More informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
More informationChapter 5. Numerical Methods: Finite Differences
Chapter 5 Numerical Methods: Finite Differences As you know, the differential equations that can be solved by an explicit analytic formula are few and far between. Consequently, the development of accurate
More information1D heat conduction problems
Chapter 1D heat conduction problems.1 1D heat conduction equation When we consider one-dimensional heat conduction problems of a homogeneous isotropic solid, the Fourier equation simplifies to the form:
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationOne dimensional steady state diffusion, with and without source. Effective transfer coefficients
One dimensional steady state diffusion, with and without source. Effective transfer coefficients 2 mars 207 For steady state situations t = 0) and if convection is not present or negligible the transport
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 16
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 16 In the last lectures, we have seen one-dimensional boundary value
More informationDual Reciprocity Boundary Element Method for Magma Ocean Simulations
Dual Reciprocity Boundary Element Method for Magma Ocean Simulations Tyler W. Drombosky drombosk@math.umd.edu Saswata Hier-Majumder saswata@umd.edu 28 April 2010 Physical Motivation Earth s early history
More informationVerification of Sound Absorption Characteristics Constituted Porous Structure
Verification of Sound Absorption Characteristics Constituted Porous Structure Toru Yoshimachi 1, Ryo Ishii 1, Kuniharu Ushijima 2, Naoki Masuda 2, Takao Yamaguchi 3, Yun Huang 4, Zhe Cui 4 1 JSOL Corporation,
More informationExplicit and implicit FEM-FCT algorithms with flux linearization
Explicit and implicit FEM-FCT algorithms with flux linearization Dmitri Kuzmin Institute of Applied Mathematics (LS III), Dortmund University of Technology Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationHigher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives
International Mathematical Forum, 2, 2007, no. 67, 3339-3350 Higher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives Jürgen Geiser Department of Mathematics
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationDiffusion Processes. Lectures INF2320 p. 1/72
Diffusion Processes Lectures INF2320 p. 1/72 Lectures INF2320 p. 2/72 Diffusion processes Examples of diffusion processes Heat conduction Heat moves from hot to cold places Diffusive (molecular) transport
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 informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationLinear Regression Linear Least Squares
Linear Regression Linear Least Squares ME 120 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu ME120: Linear Regression Introduction Introduction
More information