1D Heat equation and a finite-difference solver
|
|
- Robyn Jessica Pierce
- 5 years ago
- Views:
Transcription
1 1D Heat equation and a finite-difference solver Guillaume Riflet MARETEC IST 1 The advection-diffusion equation The original concept applied to a property within a control volume from which is derived the integral advection-diffusion equation states as {Rate of change in time} = {Ingoing Outgoing fluxes} + {Created Destroyed}. 1) Annotated in a correct mathematical encapsulation, equation 1 yields d P d = P v n) ds K P n) ds + Sc Sk) d, 2) where P is the transported property concentration, v is the advecting field, K is the diffusivity coefficient, is the surface of the control volume, n is the outward normal to the control surface and Sc and Sk are the source and sink terms, respectively. The first term on the right hand side RHS) of equation 2 states ingoing and outgoing fluxes due to advection, the second term in the RHS states ingoing and outgoing fluxes due to diffusion, and the last term on the RHS accounts for source and sink terms. the Note that Fick s law, K P, is applied to mathematically describe diffusion [1]. By resorting to the divergence theorem, stating that the divergence of a vector field inside any finite gaussian volume is equal to its flux through the boundary of the volume, i.e..e d = E n ds, 3) 1
2 1 The advection-diffusion equation 2 being the gaussian volume, its boundary, E the vector field, and E n its flux through the boundary. n is defined as the external normal unit vector to the boundary. Hence, by applying equation3) to equation2), the finite volume formulation of equation 4) is obtained: d P d = v P ) d K P ) d + Sc Sk) d. 4) By resorting to the Leibniz integration rule, d d fx, t) fx, t)d = d, which is true as long as the integral volume and consequently x) is held fixed in time, in one hand. By joining all the integrals in into a single integral on the other, equation 4 yields { } P + v P ) + K P ) Sc + Sk d = 0. Since the above integral holds zero for all gaussian volume, then its integrand must be zero, thus yielding the differential equation of advectiondiffusion: P + v P ) + K P ) Sc + Sk = 0. 5) A few comments: first, the control volume is held fixed in time. Second, the time derivative applied to a function varying in time only in its explicit component such as P x, t)), may be annotated as a partial derivative,, without loss of generality. Third, the total derivative and the partial derivative are related by the Leibniz chain rule, d P xt), t) = P xt 0), t) + d xt) P xt 0 ), t 0 ) x + d yt) P xt 0 ), t 0 ) y + d zt) P xt 0 ), t 0 ). z
3 1 The advection-diffusion equation 3 If dx = 0, then the partial derivative equals the total derivative, which is the case of P within the fixed volume. By noting that equation 5 returns v P ) = v P + P v, P + v P + P v + K P ) Sc + Sk = 0. By noting that the material derivative is defined by equation 5 finally yields D Dt + v, DP + P v + K P ) Sc + Sk = 0. Dt In incompressible fluids we have v = 0, hence in such case DP + K P ) Sc + Sk = 0. 6) Dt Thus, to sum up, equation 2) is the integral equation of advection and diffusion with source and sink terms) in flux form, equation 4) is the integral equation of advection-diffusion in volume form, equation 5) is the differential equation of advection-diffusion where the advective field is the velocity of the fluid particles relative to a fixed reference frame) and equation 6) is the differential equation of advection-diffusion for an incompressible fluid. 1.1 The 1D Heat-equation The 1D heat equation consists of property P as being temperature T or heat applying for the unidimensional case of differential equation 5). The tridimensional case with a decay or sink term writes T + µ T ) + k T = 0, 7) where µ is the diffusivity coefficient and k is the decay coefficient. In the problem, only the 1D equation is relevant: T T µ 2 + k T = 0. 8) x2 This equation admits solutions under certain initial values and certain boundary conditions google wikipedia heat equation).
4 2 Finite-difference solvers 4 2 Finite-difference solvers In order to solve numerically equations like equation 8), it is necessary to discretize the derivative operators. After discretizing the derivative operators, it is also useful to have an estimate of the error of the solution. When one knows its analytical solution, then it suffices to compute the RMS root mean square error) in order to know exactly what is the error. However, often, numerical solvers are used because there are no known analytical solution to the PDE partial differential equation), thus there is no alternative but to estimate the error. A traditional approach consists in using constant discrete time steps and constant discrete spatial steps dx combined with the Taylor serie expansion definition of analytical functions [2]. The Taylor serie expansion of an analytical real function forward in its variable coordinate is f t + ) = f n) t) )n 9) n! n=0 and backward in its variable coordinate is f t ) = f n) t) )n. 10) n! 2.1 Explicit method n=0 The Taylor series in equation 9) expanded up to the second order with third order error), f t + ) = ft) + f t) + f t) )2 2 + o ) 3). 11) A first derivative approximation may be obtained from equation 11, f t) = f t + ) f t) + o ). 12) The first derivative approximation in equation 12) yields first order error. This numerical method is often referred to as the explicit method or the forward in time method. Note that o) 3 ) = o) 2 ), o) = o ) = o), are general properties of truncature error operator.
5 2 Finite-difference solvers Implicit method Likewise, the Taylor series in equation 10) expanded up to the second order with third order error) yields, f t ) = ft) + f t) ) + f t) )2 2 + o ) 3). 13) A first derivative approximation may be obtained from equation 13, f t) = f t) f t ) + o ). 14) This numerical method is referred to as the implicit method or the backward in time method and is also first order error. 2.3 Mid-point method By subtracting equation 13) from equation 11), the mid-point method is obtained: ft + ) ft ) = f t) 2 + o ) 3), 15) f t) = ft + ) ft ) 2 + o ) 2). 16) Equation 16) is referred to as the mid-point method and has a second order error. 2.4 Second derivative discretization By summing equation 13) from equation 11), a discretization of the second derivative is obtained: ft + ) + ft ) = 2 ft) + f t) ) 2 + o ) 4). 17) Note that the third order terms in the forward and backward Taylor series expansion cancel out. In fact, every odd order term are cancelled when summing the forward and backward Taylor series expansion, whereas its the even terms that are cancelled when subtracting the backward Taylor series from the forward Taylor series. Thus, equation 17) yields f t) = ft + ) 2 ft) + ft ) ) 2 + o ) 2). 18)
6 3 Discretizing the heat equation 6 3 Discretizing the heat equation The idea is to discretize the heat equation 8) with a numerical scheme forward in time and centred in space FTCS). Thus by using the scheme in equation 12) to the time derivative and by using the numerical scheme in equation 18) to second-derivative in space, equation 8) becomes, in its algebraic form: T x, t + t) T x, t) = t T x + x, t) 2 T x, t) + T x x, t) µ x 2 k T x, t), 19) which in turn yields or even yet, T x, t + t) = T x, t) + µ t T x + x, t) 2 T x, t) + T x x, t)) x2 k t T x, t), 20) T x, t + t) = µ t T x + x, t) x2 + 1 k t 2 µ t ) x 2 T x, t) + µ t T x x, t). 21) x2 The numerical scheme is defined positive and physically realist) if and only if every coefficient affected to T x + x, t), T x, t) and T x x, t) is positive. Which is the case in equation 21) as long as 3.1 the grid 1 k t 2 µ t 1. 22) x2 Consider an even spaced bi-dimensional grid evolving over space by increments of x and over time by increments of t. Each grid point is indexed by integers i in space and n in time. The spatial axis of the grid represents the heated bar. The time axis of the grid represents the temporal evolution
7 3 Discretizing the heat equation 7 of the heat in the bar. Thus, given the values of temperature in the bar for all values of index i at a given instant n, Ti n, the temperature in the next time instant is Ti n+1 is deduced from equation 21), T n+1 i = Dif Ti+1 n + 1 k t 2 Dif) Ti n + Dif Ti 1, n 23) where, by definition, Dif µ t x Boundary conditions For boundary conditions, a null diffusive flux is considered as the heated bar is adiabatically insulated from the environment), thus a special scheme must be considered at the boundary grid points for i = 1 and i = Ni where Ni is the total number of grid points in the spatial axis. i = 1 T1 n+1 = Dif T2 n + 1 k t Dif) T1 n, i = Ni T n+1 Ni = 1 k t Dif) TNi n + Dif TNi 1. n Notice how the central diffusive flux coefficient was cut by half at the boundaries when compared with equation 23). Indeed the net diffusive flux is represented by an ingoing diffusive flux represented by the terms T i+1 and T i 1, and by an outgoing diffusive flux represented by the term T i. Thus, the outgoing diffusive represents the amount of heat from T i that flows to both the neighbour grid points i + 1 and i 1. At the boundaries, however, the outgoing diffusive flux flows to only one neighbouring grid point, thus cutting in half the diffusion that occurs. 3.3 Initial condition Equation 23) calculates T n+1 as a function of T n. Thus, at some point, T 1 needs to be defined. It s called the initial condition. The initial condition will describe an initial state of the modelled bar. 3.4 Continuous emission Suppose that there s an heat source at some point near the bar. This heat source will input heat continuously in the bar. A possible way to model a continuous emission is as follows: first, the regular diffusivity-decay algorithm is employed, T n+1 i = Dif Ti+1 n + 1 k t 2 Dif) Ti n + Dif Ti 1, n
8 3 Discretizing the heat equation 8 Then, the continous source term C is added at the given grid point, say i = p, Ti n+1 = Ti n+1 i p Ti n+1 = Ti n+1 + C t i = p. References [1] Chapra, S. C. Surface water-quality modeling. Mcgraw-hill Series In Water Resources And Environmental Engineering 1997). [2] Fletcher, C. A. J., and Srinivas, K. Computational Techniques for Fluid Dynamics 1. Springer, 1991.
Numerical 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 information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More information2. Conservation of Mass
2 Conservation of Mass The equation of mass conservation expresses a budget for the addition and removal of mass from a defined region of fluid Consider a fixed, non-deforming volume of fluid, V, called
More information25. Chain Rule. Now, f is a function of t only. Expand by multiplication:
25. Chain Rule The Chain Rule is present in all differentiation. If z = f(x, y) represents a two-variable function, then it is plausible to consider the cases when x and y may be functions of other variable(s).
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
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 informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationCharacteristic finite-difference solution Stability of C C (CDS in time/space, explicit): Example: Effective numerical wave numbers and dispersion
Spring 015 Lecture 14 REVIEW Lecture 13: Stability: Von Neumann Ex.: 1st order linear convection/wave eqn., F-B scheme Hyperbolic PDEs and Stability nd order wave equation and waves on a string Characteristic
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationMath 212-Lecture 8. The chain rule with one independent variable
Math 212-Lecture 8 137: The multivariable chain rule The chain rule with one independent variable w = f(x, y) If the particle is moving along a curve x = x(t), y = y(t), then the values that the particle
More informationDiffusion of a density in a static fluid
Diffusion of a density in a static fluid u(x, y, z, t), density (M/L 3 ) of a substance (dye). Diffusion: motion of particles from places where the density is higher to places where it is lower, due to
More informationProblem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0139.9 NUMERICAL FLUID MECHANICS SPRING 015 Problem Set 4 Issued: Wednesday, March 18, 015 Due: Wednesday,
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationLecture II: Vector and Multivariate Calculus
Lecture II: Vector and Multivariate Calculus Dot Product a, b R ' ', a ( b = +,- a + ( b + R. a ( b = a b cos θ. θ convex angle between the vectors. Squared norm of vector: a 3 = a ( a. Alternative notation:
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationComputational Fluid Dynamics Prof. Dr. SumanChakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. SumanChakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. #11 Fundamentals of Discretization: Finite Difference
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationHow to Use Calculus Like a Physicist
How to Use Calculus Like a Physicist Physics A300 Fall 2004 The purpose of these notes is to make contact between the abstract descriptions you may have seen in your calculus classes and the applications
More informationDivergence Theorem December 2013
Divergence Theorem 17.3 11 December 2013 Fundamental Theorem, Four Ways. b F (x) dx = F (b) F (a) a [a, b] F (x) on boundary of If C path from P to Q, ( φ) ds = φ(q) φ(p) C φ on boundary of C Green s Theorem:
More information1 Introduction to MATLAB
L3 - December 015 Solving PDEs numerically (Reports due Thursday Dec 3rd, carolinemuller13@gmail.com) In this project, we will see various methods for solving Partial Differential Equations (PDEs) using
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationDivergence Theorem Fundamental Theorem, Four Ways. 3D Fundamental Theorem. Divergence Theorem
Divergence Theorem 17.3 11 December 213 Fundamental Theorem, Four Ways. b F (x) dx = F (b) F (a) a [a, b] F (x) on boundary of If C path from P to Q, ( φ) ds = φ(q) φ(p) C φ on boundary of C Green s Theorem:
More informationSummary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization:
The LSZ reduction formula based on S-5 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free theory:
More informationIdeas from Vector Calculus Kurt Bryan
Ideas from Vector Calculus Kurt Bryan Most of the facts I state below are for functions of two or three variables, but with noted exceptions all are true for functions of n variables..1 Tangent Line Approximation
More informationIntroduction to Heat and Mass Transfer. Week 8
Introduction to Heat and Mass Transfer Week 8 Next Topic Transient Conduction» Analytical Method Plane Wall Radial Systems Semi-infinite Solid Multidimensional Effects Analytical Method Lumped system analysis
More informationModule 3: BASICS OF CFD. Part A: Finite Difference Methods
Module 3: BASICS OF CFD Part A: Finite Difference Methods THE CFD APPROACH Assembling the governing equations Identifying flow domain and boundary conditions Geometrical discretization of flow domain Discretization
More information1. Consider the initial value problem: find y(t) such that. y = y 2 t, y(0) = 1.
Engineering Mathematics CHEN30101 solutions to sheet 3 1. Consider the initial value problem: find y(t) such that y = y 2 t, y(0) = 1. Take a step size h = 0.1 and verify that the forward Euler approximation
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 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 informationSKMM 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SKMM 3023 Applied Numerical Methods Numerical Differentiation ibn Abdullah Faculty of Mechanical Engineering Òº ÙÐÐ ÚºÒÙÐÐ ¾¼½ SKMM 3023 Applied Numerical Methods Numerical
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
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 information1.061 / 1.61 Transport Processes in the Environment
MIT OpenCourseWare http://ocw.mit.edu 1.061 / 1.61 Transport Processes in the Environment Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Solution
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationENGI Partial Differentiation Page y f x
ENGI 344 4 Partial Differentiation Page 4-0 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
More informationCURVILINEAR MOTION: RECTANGULAR COMPONENTS (Sections )
CURVILINEAR MOTION: RECTANGULAR COMPONENTS (Sections 12.4-12.5) Today s Objectives: Students will be able to: a) Describe the motion of a particle traveling along a curved path. b) Relate kinematic quantities
More information16.4. Power Series. Introduction. Prerequisites. Learning Outcomes
Power Series 6.4 Introduction In this Section we consider power series. These are examples of infinite series where each term contains a variable, x, raised to a positive integer power. We use the ratio
More informationAdvection, Conservation, Conserved Physical Quantities, Wave Equations
EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University Contents EP711 Supplementary
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 7
Numerical Fluid Mechanics Fall 2011 Lecture 7 REVIEW of Lecture 6 Material covered in class: Differential forms of conservation laws Material Derivative (substantial/total derivative) Conservation of Mass
More informationNumerical Analysis and Computer Science DN2255 Spring 2009 p. 1(8)
Numerical Analysis and Computer Science DN2255 Spring 2009 p. 1(8) Contents Important concepts, definitions, etc...2 Exact solutions of some differential equations...3 Estimates of solutions to differential
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
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 informationCURVILINEAR MOTION: GENERAL & RECTANGULAR COMPONENTS
CURVILINEAR MOTION: GENERAL & RECTANGULAR COMPONENTS Today s Objectives: Students will be able to: 1. Describe the motion of a particle traveling along a curved path. 2. Relate kinematic quantities in
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 informationScattering amplitudes and the Feynman rules
Scattering amplitudes and the Feynman rules based on S-10 We have found Z( J ) for the phi-cubed theory and now we can calculate vacuum expectation values of the time ordered products of any number of
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20
2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20 REVIEW Lecture 19: Finite Volume Methods Review: Basic elements of a FV scheme and steps to step-up a FV scheme One Dimensional examples d x j x j 1/2
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
More informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
More informationIntegration in the Complex Plane (Zill & Wright Chapter 18)
Integration in the omplex Plane Zill & Wright hapter 18) 116-4-: omplex Variables Fall 11 ontents 1 ontour Integrals 1.1 Definition and Properties............................. 1. Evaluation.....................................
More informationDefinition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.
5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that
More informationCURVILINEAR MOTION: GENERAL & RECTANGULAR COMPONENTS
CURVILINEAR MOTION: GENERAL & RECTANGULAR COMPONENTS Today s Objectives: Students will be able to: 1. Describe the motion of a particle traveling along a curved path. 2. Relate kinematic quantities in
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationPV Generation in the Boundary Layer
1 PV Generation in the Boundary Layer Robert Plant 18th February 2003 (With thanks to S. Belcher) 2 Introduction How does the boundary layer modify the behaviour of weather systems? Often regarded as a
More informationNotation Nodes are data points at which functional values are available or at which you wish to compute functional values At the nodes fx i
LECTURE 6 NUMERICAL DIFFERENTIATION To find discrete approximations to differentiation (since computers can only deal with functional values at discrete points) Uses of numerical differentiation To represent
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationJUST THE MATHS UNIT NUMBER INTEGRATION APPLICATIONS 10 (Second moments of an arc) A.J.Hobson
JUST THE MATHS UNIT NUMBER 13.1 INTEGRATION APPLICATIONS 1 (Second moments of an arc) by A.J.Hobson 13.1.1 Introduction 13.1. The second moment of an arc about the y-axis 13.1.3 The second moment of an
More informationMATH 3150: PDE FOR ENGINEERS FINAL EXAM (VERSION D) 1. Consider the heat equation in a wire whose diffusivity varies over time: u k(t) 2 x 2
MATH 35: PDE FOR ENGINEERS FINAL EXAM (VERSION D). Consider the heat equation in a wire whose diffusivity varies over time: u t = u k(t) x where k(t) is some positive function of time. Assume the wire
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationIntroduction to Environment System Modeling
Introduction to Environment System Modeling (3 rd week:modeling with differential equation) Department of Environment Systems, Graduate School of Frontier Sciences, the University of Tokyo Masaatsu AICHI
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
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 informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationConvection and buoyancy oscillation
Convection and buoyancy oscillation Recap: We analyzed the static stability of a vertical profile by the "parcel method"; For a given environmental profile (of T 0, p 0, θ 0, etc.), if the density of an
More informationIntroduction to Fluid Dynamics
Introduction to Fluid Dynamics Roger K. Smith Skript - auf englisch! Umsonst im Internet http://www.meteo.physik.uni-muenchen.de Wählen: Lehre Manuskripte Download User Name: meteo Password: download Aim
More information14.7 The Divergence Theorem
14.7 The Divergence Theorem The divergence of a vector field is a derivative of a sort that measures the rate of flow per unit of volume at a point. A field where such flow doesn't occur is called 'divergence
More informationx 2 y = 1 2. Problem 2. Compute the Taylor series (at the base point 0) for the function 1 (1 x) 3.
MATH 8.0 - FINAL EXAM - SOME REVIEW PROBLEMS WITH SOLUTIONS 8.0 Calculus, Fall 207 Professor: Jared Speck Problem. Consider the following curve in the plane: x 2 y = 2. Let a be a number. The portion of
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 5
.9 Numerical Fluid Mechanics Fall 011 Lecture 5 REVIEW Lecture 4 Roots of nonlinear equations: Open Methods Fixed-point Iteration (General method or Picard Iteration), with examples Iteration rule: x g(
More information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
More informationMulti-physics Modeling Using Cellular Automata
Multi-physics Modeling sing Cellular Automata Brian Vic Mechanical Engineering Department, Virginia Tech, Blacsburg, VA 246-238 This paper proposes a new modeling and solution method that is relatively
More informationThe Finite Volume Mesh
FMIA F Moukalled L Mangani M Darwish An Advanced Introduction with OpenFOAM and Matlab This textbook explores both the theoretical foundation of the Finite Volume Method (FVM) and its applications in Computational
More informationChain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics
3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)
More informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
More informationLecture Notes 6: Dynamic Equations Part A: First-Order Difference Equations in One Variable
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 54 Lecture Notes 6: Dynamic Equations Part A: First-Order Difference Equations in One Variable Peter J. Hammond latest revision 2017
More informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More informationFirst-Order Differential Equations
CHAPTER 1 First-Order Differential Equations 1. Diff Eqns and Math Models Know what it means for a function to be a solution to a differential equation. In order to figure out if y = y(x) is a solution
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
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 informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 6. Solve the completely inhomogeneous diffusion problem on the half-line
Strauss PDEs 2e: Section 3.3 - Exercise 2 Page of 6 Exercise 2 Solve the completely inhomogeneous diffusion problem on the half-line v t kv xx = f(x, t) for < x
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
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 information