Process of whiskey maturation & Current Efforts to improve Comparison Goals Modeling nonlinear Diffusion
|
|
- Anastasia Barnett
- 5 years ago
- Views:
Transcription
1 Brian England
2 Objective Process of whiskey maturation & Current Efforts to improve Comparison Goals Modeling nonlinear Diffusion Cylindrical Coordinates Initial and Boundary Conditions Methodologies and Computational Results Finite Difference Finite Volume Function Space Final Comparison and Conclusion
3 Maturation of Whiskey Driven by two processes Diffusion of Oak barrel goodness Short time scales Chemical reaction Long time scales Current Efforts Diffusion Wooden inserts, tastes more processed and is lacking Chemical Reactions Pressure vessels / burners to shift % yields and reaction rates Much better and they re improving on this front
4 Comparison Goal General Comparisons Accuracy Stability Computational Efficiency I will primarily deal with Computational Efficiency Validate a methodology for utilization in future work
5 Modeling Nonlinear Diffusion What makes the Diffusion Nonlinear? Changes in differential operators Best coordinate system Ellipsoidal Much more Complex Cylindrical Coordinates Easier & exact solution are readily available
6 Cylindrical Coordinates Differential Operator Adjustments Gradient Divergence Curl Laplacian U 1 U U U,, r r z F 1 rfr 1F F r r r z 1Fz F Fr Fz 1 rf F r F,, r z z r r r 1 ru 1 U U U r r r r z z
7 Initial and Boundary Conditions First Comparisons Boundary Conditions Chosen to ensure an equilibrium U( t,1,, z) cos sin z U( t, r,, 1) cosr Initial Condition U(0, r,, z) 0
8 Finite Difference Approach General PDE t r r r r z U 1 U U 1 U U Centered Difference Spatial Discretizations U r U U U i1, j, k i, j, k i1, j, k Forward Euler in Time r U f U r U U t f U r t U r n1 n (, ) ( ) (, ) U U i1, j, k i1, j, k r
9 Stability The discretization n n n U U n1 n Ui 1U i U i1 U U ( t) t x x Requirement for Stability t x For our PDE, the largest linearized coefficient is the most restrictive 4 1 U 1 r t r min r
10 Good news (->Stability) Can get away with ~ dr = 0.05 & dt = Wall Clock time = Seconds Tolerance set to Flux/function values of 0.01
11 Finite Volume Method More work apriori! Approach PDE with Gauss Law & Integral Form Our PDE in Integral Form V U dv K U dv t V The Integral Form allows us to represent the average value of our function U 1 U dv UdV V UdV V t t t V t V V V
12 Final Exact Formulation Utilizing Gauss Law V FdV F ds S Making appropriate substitutions V KU dv K U ds U 6 K ds We obtain the final form U S t V i 1 S i i
13 Finite Volume Approximations Difference in time K Net Flux n 1 n U U t V Differences on the boundaries U r U n i1 U r n i Calculation of Fluxes S ( r dr / )( U U ) ( r dr / )( U U ) Net Radial Sr U U U Net Azimuthal r j1 j i1 r i1 i i i1 S U U U Net Axial Z j1 j i1
14 Cell Volume & Surface Integrals Volume 1 V rout rin z rzr Axial Surface Radial Surfaces 1 SZ rout rin rr S rz r Azimuthal Surfaces S r z
15 Items of Note Method is fully conservative Only approximations lie in the derivative at the cell surfaces When proper substitutions of the volume / surfaces is performed, it s almost identical to finite difference What benefits are there? We re still approximating derivatives Leads to order of accuracies Computations is about the same
16 Results Same scenario as finite difference Stability Required dt = /6 th the dt There may have been corner issues Wall Clock 359 Seconds 10 times as long Even if the dt was matched, this method would still be slower For now
17 Function Space Methods Rely on Eigen function products Coefficients are solved via inner products Integrals are approximated by function evaluations at collocation points Our Problem - Diffusion Equation in Cylindrical Coordinates Choice of Eigenfunctions Radial - Bessel functions of the first kind Azimuthal Trigonometric Axial Trigonometric Temporal - Exponential
18 Exact Solution Fundamental Technique Principle of Superposition Handle a separate problem for each boundary and initial state for a steady state solution and transient U( t, r,,1) f ( r, ) top U( t,1,, z) f (, z) Lat U( t, r,, 1) f ( r, ) bot U(0, r,, z) f ( r,, z) Lat
19 General Solution Step by step Separation of Variables First Time Which Gives Separating and solving T() t e t U( t, r,, z) T( t) ( r,, z) Tt () t ( r,, z) T( t) ( r,, z) r r r r z Letting ( r,, z) R( r) ( ) Z( z) 1 R( r) 1 R( r) 1 ( ) 1 Z( z) rr( r) r R( r) r r ( ) Z( z) z
20 General Solution You can show that our ODE s 1 ( ) ( ) l 1 Z( z) Z() z z m Provides R( r) R( r) r r ( m) r l R( r) 0 r r General Solution depends on our BC-s
21 Solution Primary Case - Initial Conditions with homogenous BC=0 lmnt U ( t, r,, z) Almne cos mz Jl m r S in l l0 m0 n0 l0 m0 n0 l = is an integer greater than or equal to 0. k ln,, is the n th root of the Bessel function of the first kind of order L m is an integer greater than or equal to 0. Finally lmnt Blmne cos mz Jl m r Cos l k m lmn l, n
22 Constants The constants are solved via A lmn U (0, r,, z)cos mz Jl m r sin l rdrddz cos mz Jl m r sin l rdrddz B lmn U (0, r,, z)cos mz Jl m r cos l rdrddz cos mz Jl m r cos l rdrddz
23 Collocation Points Azimuthal Fourier j Evenly distributed Axial fourier Chebyshev nodal locations Radial Bessel / N j 1 cos N z Roots of m th order Bessel function Best used in tabular form to save on computations x j
24 Collocation Methods Numerous Integration techniques Quadrature at collocation points for A and B Standard Gaussian Quadrature Radial and Axial Weights Polynomial Approximations Azimuthal Weights Trigonometric Approximation We ll use nodal locations, i,j,k, to calculate weights f ( r,, z) rdrddz w f ( r,, z ) i1 j1 k1 ijk i j k
25 A Separation of Variables lmn Through some integration and appropriate initial conditions we can state B lmn l Z ( z)cos mz dz ( )sin l d R ( r) J m r rdr J l We can now approach each integral with our approximation m r rdr l Z ( z)cos mz dz ( )cos l d R ( r) J m r rdr J l m r rdr
26 Weight Calculation (Axial) A bit of matrix manipulation and brute force With fixed nodes 1 inverse and done! c1 1 1 x1 x... xn 1 x N c N N N N 1 x1 x... xn 1 xn cn 1 ( N 1) N 1 N 1 N 1 N 1 1 x1 x... xn 1 x N c N N n
27 Weight Calculation (Azimuthal) Trigonometric Gaussian Quadrature Maximal Trigonometric Degree of Exactness New Quadrature Methodologies Trigonometric Orthogonal Systems and Quadrature Formulae with Maximal Trigonometric Degree of Exactness Gradimir V. Milovanovic Easier assumption Requires proper initial conditions w1 1 1 sin x1 sin x... sin xn 1 sin x N w N N N N 1 sin x1 sin x... sin xn 1 sin xn wn 1 ( N 1) N 1 N 1 N 1 N 1 1 sin x1 sin x... sin xn 1 sin x N w N N sin( n ) d wif ( i) 1 n1 i1
28 Weight Calculation (Radial) Requires Bessel functions 1 0 J ( r) dr w f ( r ) n i n i i1 m mn x m! ( m n 1) m! ( m n 1) m n * J 0 n() r dr J N dr 0 m0 m0 J ( x ) J ( x )... J ( x ) J ( x ) w J J x J x J x J x w J J ( x ) J ( x )... J ( x ) J ( x ) w J J x J x J x J x w J * N1 0 N 1 1 * 1( 1) 1( )... 1( N1) 1( N) * N 1 1 N 1 N 1 N 1 N 1 N N 1 N 1 * N ( 1) N ( )... N ( N 1) N ( N ) N N m
29 Conclusion Function Space methods Very time consuming apriori Didn t get to finish Expectation of fast convergence Approximation lies only with initial / boundary conditions Exact solution for all time after (for the given approx)
2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationBOUNDARY PROBLEMS IN HIGHER DIMENSIONS. kt = X X = λ, and the series solutions have the form (for λ n 0):
BOUNDARY PROBLEMS IN HIGHER DIMENSIONS Time-space separation To solve the wave and diffusion equations u tt = c u or u t = k u in a bounded domain D with one of the three classical BCs (Dirichlet, Neumann,
More informationSteady and unsteady diffusion
Chapter 5 Steady and unsteady diffusion In this chapter, we solve the diffusion and forced convection equations, in which it is necessary to evaluate the temperature or concentration fields when the velocity
More informationBessel Functions. A Touch of Magic. Fayez Karoji 1 Casey Tsai 1 Rachel Weyrens 2. SMILE REU Summer Louisiana State University
Bessel s Function A Touch of Magic Fayez Karoji 1 Casey Tsai 1 Rachel Weyrens 2 1 Department of Mathematics Louisiana State University 2 Department of Mathematics University of Arkansas SMILE REU Summer
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More informationChapter 3 - Vector Calculus
Chapter 3 - Vector Calculus Gradient in Cartesian coordinate system f ( x, y, z,...) dr ( dx, dy, dz,...) Then, f f f f,,,... x y z f f f df dx dy dz... f dr x y z df 0 (constant f contour) f dr 0 or f
More informationIV. Transport Phenomena Lecture 18: Forced Convection in Fuel Cells II
IV. Transport Phenomena Lecture 18: Forced Convection in Fuel Cells II MIT Student (and MZB) As discussed in the previous lecture, we are interested in forcing fluid to flow in a fuel cell in order to
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationDUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs
DUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs MATH 4354 Fall 2005 December 5, 2005 1 Duhamel s Principle
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
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 information3: Mathematics Review
3: Mathematics Review B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 015 Sept.-Dec. 015 September 1 Review of: Table of Contents Co-ordinate systems (Cartesian,
More informationMa 530. Partial Differential Equations - Separation of Variables in Multi-Dimensions
Ma 530 Partial Differential Equations - Separation of ariables in Multi-Dimensions Temperature in an Infinite Cylinder Consider an infinitely long, solid, circular cylinder of radius c with its axis coinciding
More informationSeparation of Variables
Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us
More informationSession 5 Heat Conduction in Cylindrical and Spherical Coordinates I
Session 5 Heat Conduction in Cylindrical and Spherical Coordinates I 1 Introduction The method of separation of variables is also useful in the determination of solutions to heat conduction problems in
More informationSummary for Vector Calculus and Complex Calculus (Math 321) By Lei Li
Summary for Vector alculus and omplex alculus (Math 321) By Lei Li 1 Vector alculus 1.1 Parametrization urves, surfaces, or volumes can be parametrized. Below, I ll talk about 3D case. Suppose we use e
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
More informationMATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued).
MATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued). Bessel s differential equation of order m 0: z 2 d2 f dz 2 + z df dz + (z2 m 2 )f = 0 The equation is considered on the
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 18 Properties of Sound Sound Waves Requires medium for propagation Mainly
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 informationCHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 5-1 Road Map of the Lecture V Laplace Transform and Transfer
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationVector Calculus. A primer
Vector Calculus A primer Functions of Several Variables A single function of several variables: f: R $ R, f x (, x ),, x $ = y. Partial derivative vector, or gradient, is a vector: f = y,, y x ( x $ Multi-Valued
More informationAn Introduction to Partial Differential Equations in the Undergraduate Curriculum
An Introduction to Partial Differential Equations in the Undergraduate Curriculum John C. Polking LECTURE 12 Heat Transfer in the Ball 12.1. Outline of Lecture The problem The problem for radial temperatures
More informationDisclaimer: This Final Exam Study Guide is meant to help you start studying. It is not necessarily a complete list of everything you need to know.
Disclaimer: This is meant to help you start studying. It is not necessarily a complete list of everything you need to know. The MTH 234 final exam mainly consists of standard response questions where students
More information송석호 ( 물리학과 )
http://optics.hanyang.ac.kr/~shsong 송석호 ( 물리학과 ) Field and Wave Electromagnetics, David K. Cheng Reviews on (Week 1). Vector Analysis 3. tatic Electric Fields (Week ) 4. olution of Electrostatic Problems
More informationMathematics of Physics and Engineering II: Homework problems
Mathematics of Physics and Engineering II: Homework problems Homework. Problem. Consider four points in R 3 : P (,, ), Q(,, 2), R(,, ), S( + a,, 2a), where a is a real number. () Compute the coordinates
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationIntroduction to Vector Calculus (29) SOLVED EXAMPLES. (d) B. C A. (f) a unit vector perpendicular to both B. = ˆ 2k = = 8 = = 8
Introduction to Vector Calculus (9) SOLVED EXAMPLES Q. If vector A i ˆ ˆj k, ˆ B i ˆ ˆj, C i ˆ 3j ˆ kˆ (a) A B (e) A B C (g) Solution: (b) A B (c) A. B C (d) B. C A then find (f) a unit vector perpendicular
More informationSecond-Order Linear ODEs (Textbook, Chap 2)
Second-Order Linear ODEs (Textbook, Chap ) Motivation Recall from notes, pp. 58-59, the second example of a DE that we introduced there. d φ 1 1 φ = φ 0 dx λ λ Q w ' (a1) This equation represents conservation
More informationEigenfunctions on the surface of a sphere. In spherical coordinates, the Laplacian is. u = u rr + 2 r u r + 1 r 2. sin 2 (θ) + 1
Eigenfunctions on the surface of a sphere In spherical coordinates, the Laplacian is u = u rr + 2 r u r + 1 r 2 [ uφφ sin 2 (θ) + 1 sin θ (sin θ u θ) θ ]. Eigenfunctions on the surface of a sphere In spherical
More informationENGI Gradient, Divergence, Curl Page 5.01
ENGI 94 5. - Gradient, Divergence, Curl Page 5. 5. The Gradient Operator A brief review is provided here for the gradient operator in both Cartesian and orthogonal non-cartesian coordinate systems. Sections
More informationM445: Heat equation with sources
M5: Heat equation with sources David Gurarie I. On Fourier and Newton s cooling laws The Newton s law claims the temperature rate to be proportional to the di erence: d dt T = (T T ) () The Fourier law
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 informationMATH 241 Practice Second Midterm Exam - Fall 2012
MATH 41 Practice Second Midterm Exam - Fall 1 1. Let f(x = { 1 x for x 1 for 1 x (a Compute the Fourier sine series of f(x. The Fourier sine series is b n sin where b n = f(x sin dx = 1 = (1 x cos = 4
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationCHAPTER 4. Introduction to the. Heat Conduction Model
A SERIES OF CLASS NOTES FOR 005-006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 4 A COLLECTION OF HANDOUTS ON PARTIAL DIFFERENTIAL EQUATIONS
More informationTransient Heat Conduction in a Circular Cylinder
Transient Heat Conduction in a Circular Cylinder The purely radial 2-D heat equation will be solved in cylindrical coordinates using variation of parameters. Assuming radial symmetry the solution is represented
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationQuadrature Formula for Computed Tomography
Quadrature Formula for Computed Tomography orislav ojanov, Guergana Petrova August 13, 009 Abstract We give a bivariate analog of the Micchelli-Rivlin quadrature for computing the integral of a function
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 informationPartial Differential Equations
Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives with respect to those variables. Most (but
More information3.3 Unsteady State Heat Conduction
3.3 Unsteady State Heat Conduction For many applications, it is necessary to consider the variation of temperature with time. In this case, the energy equation for classical heat conduction, eq. (3.8),
More informationFinal Examination Solutions
Math. 42 Fulling) 4 December 25 Final Examination Solutions Calculators may be used for simple arithmetic operations only! Laplacian operator in polar coordinates: Some possibly useful information 2 u
More informationChapter 5 Time-Dependent Conduction
Chapter 5 Time-Dependent Conduction 5.1 The Lumped Capacitance Method This method assumes spatially uniform solid temperature at any instant during the transient process. It is valid if the temperature
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
More informationGauss s Law. Chapter 22. Electric Flux Gauss s Law: Definition. Applications of Gauss s Law
Electric Flux Gauss s Law: Definition Chapter 22 Gauss s Law Applications of Gauss s Law Uniform Charged Sphere Infinite Line of Charge Infinite Sheet of Charge Two infinite sheets of charge Phys 2435:
More informationBENG 221 DRUG DIFFUSION THROUGH A COATED STENT
BENG 221 DRUG DIFFUSION THROUGH A COATED STENT 12 TH NOVEMBER, 2010 1. INTRODUCTION: In this report, a bioengineering application of the one- dimensional diffusion equation, solved in both rectangular
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationElectromagnetic waves in free space
Waveguide notes 018 Electromagnetic waves in free space We start with Maxwell s equations for an LIH medum in the case that the source terms are both zero. = =0 =0 = = Take the curl of Faraday s law, then
More informationMajor Ideas in Calc 3 / Exam Review Topics
Major Ideas in Calc 3 / Exam Review Topics Here are some highlights of the things you should know to succeed in this class. I can not guarantee that this list is exhaustive!!!! Please be sure you are able
More informationLiquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses
Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses by W. A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine Collaborators:
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 informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationMath 2930 Worksheet Final Exam Review
Math 293 Worksheet Final Exam Review Week 14 November 3th, 217 Question 1. (* Solve the initial value problem y y = 2xe x, y( = 1 Question 2. (* Consider the differential equation: y = y y 3. (a Find the
More informationENGI Duffing s Equation Page 4.65
ENGI 940 4. - Duffing s Equation Page 4.65 4. Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing
More informationMultiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas
Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas Michael Barnes Plasma Science & Fusion Center Massachusetts Institute of Technology Collaborators:
More informationA First Course on Kinetics and Reaction Engineering Unit D and 3-D Tubular Reactor Models
Unit 34. 2-D and 3-D Tubular Reactor Models Overview Unit 34 describes two- and three-dimensional models for tubular reactors. One limitation of the ideal PFR model is that the temperature and composition
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 informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
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 informationMcGill University April 16, Advanced Calculus for Engineers
McGill University April 16, 2014 Faculty of cience Final examination Advanced Calculus for Engineers Math 264 April 16, 2014 Time: 6PM-9PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More information1 Curvilinear Coordinates
MATHEMATICA PHYSICS PHYS-2106/3 Course Summary Gabor Kunstatter, University of Winnipeg April 2014 1 Curvilinear Coordinates 1. General curvilinear coordinates 3-D: given or conversely u i = u i (x, y,
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More informationStaple or bind all pages together. DO NOT dog ear pages as a method to bind.
Math 3337 Homework Instructions: Staple or bind all pages together. DO NOT dog ear pages as a method to bind. Hand-drawn sketches should be neat, clear, of reasonable size, with axis and tick marks appropriately
More informationPhysics 116C Solutions to the Practice Final Exam Fall The method of Frobenius fails because it works only for Fuchsian equations, of the type
Physics 6C Solutions to the Practice Final Exam Fall 2 Consider the differential equation x y +x(x+)y y = () (a) Explain why the method of Frobenius method fails for this problem. The method of Frobenius
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationAND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)
More informationNST1A: Mathematics II (Course A) End of Course Summary, Lent 2011
General notes Proofs NT1A: Mathematics II (Course A) End of Course ummar, Lent 011 tuart Dalziel (011) s.dalziel@damtp.cam.ac.uk This course is not about proofs, but rather about using different techniques.
More informationWaves on 2 and 3 dimensional domains
Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation
More information. (70.1) r r. / r. Substituting, we have the following equation for f:
7 Spherical waves Let us consider a sound wave in which the distribution of densit velocit etc, depends only on the distance from some point, ie, is spherically symmetrical Such a wave is called a spherical
More informationBoundary value problems
1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the
More informationEELE 3331 Electromagnetic I Chapter 3. Vector Calculus. Islamic University of Gaza Electrical Engineering Department Dr.
EELE 3331 Electromagnetic I Chapter 3 Vector Calculus Islamic University of Gaza Electrical Engineering Department Dr. Talal Skaik 2012 1 Differential Length, Area, and Volume This chapter deals with integration
More information(a) Determine the general solution for φ(ρ) near ρ = 0 for arbitary values E. (b) Show that the regular solution at ρ = 0 has the series expansion
Problem 1. Curious Wave Functions The eigenfunctions of a D7 brane in a curved geometry lead to the following eigenvalue equation of the Sturm Liouville type ρ ρ 3 ρ φ n (ρ) = E n w(ρ)φ n (ρ) w(ρ) = where
More informationAppendix A Vector Analysis
Appendix A Vector Analysis A.1 Orthogonal Coordinate Systems A.1.1 Cartesian (Rectangular Coordinate System The unit vectors are denoted by x, ŷ, ẑ in the Cartesian system. By convention, ( x, ŷ, ẑ triplet
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationMath 211, Fall 2014, Carleton College
A. Let v (, 2, ) (1,, ) 1, 2, and w (,, 3) (1,, ) 1,, 3. Then n v w 6, 3, 2 is perpendicular to the plane, with length 7. Thus n/ n 6/7, 3/7, 2/7 is a unit vector perpendicular to the plane. [The negation
More informationAdiabatic expansion Isothermal compression Adiabatic compression
Tropical Cyclones: Steady State Physics 1 Energy Production 2 Carnot Theorem: Maximum efficiency results from a pa rticular energy e cycle: Isothermal expansion Adiabatic expansion Isothermal compression
More informationNotes on Special Functions
Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.
More informationPONDI CHERRY UNI VERSI TY
B.Sc. ALLIED MATHEMATICS SYLLABUS 2009-2010 onwards PONDI CHERRY UNI VERSI TY PUDUCHERRY 605 014 B.Sc. ALLIED MATHEMATICS Syllabus for Allied Mathematics for B.Sc. Physics Main/Chemistry Main/Electronics
More informationMath 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning
More informationFraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, Darmstadt, Germany
Scale Space and PDE methods in image analysis and processing Arjan Kuijper Fraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, 64283
More informationThe Divergence Theorem
The Divergence Theorem 5-3-8 The Divergence Theorem relates flux of a vector field through the boundary of a region to a triple integral over the region. In particular, let F be a vector field, and let
More informationLegendre s Equation. PHYS Southern Illinois University. October 13, 2016
PHYS 500 - Southern Illinois University October 13, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 13, 2016 1 / 10 The Laplacian in Spherical Coordinates The Laplacian is given
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 informationBASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA
1 BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA This part of the Basic Exam covers topics at the undergraduate level, most of which might be encountered in courses here such as Math 233, 235, 425, 523, 545.
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationMultiple Integrals. Chapter 4. Section 7. Department of Mathematics, Kookmin Univerisity. Numerical Methods.
4.7.1 Multiple Integrals Chapter 4 Section 7 4.7.2 Double Integral R f ( x, y) da 4.7.3 Double Integral Apply Simpson s rule twice R [ a, b] [ c, d] a x, x,..., x b, c y, y,..., y d 0 1 n 0 1 h ( b a)
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
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 informationIntroduction to First Order Equations Sections
A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Introduction to First Order Equations Sections 2.1-2.3 Dr. John Ehrke Department of Mathematics Fall 2012 Course Goals The
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationINVARIANCE OF THE LAPLACE OPERATOR.
INVARIANCE OF THE LAPLACE OPERATOR. The goal of this handout is to give a coordinate-free proof of the invariance of the Laplace operator under orthogonal transformations of R n (and to explain what this
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More information