Temperature evolution of an oceanic fracture zone

Size: px
Start display at page:

Download "Temperature evolution of an oceanic fracture zone"

Transcription

1 Temperature evolution of an oceanic fracture one Xiaohua Xu & Zhao Chen October 29, 212

2 Outline Background Introduction to the Problem Green's Function and the Analytical Solution for Temperature Introduction to Green s Function Solution of Temperature by method of Green s Function Model Topography V.S. Real Topography Solution the Model Topography Topography map with the evolution of Temperature Comparison with the real Topography Conclusions and Discussions

3 1 Background

4 1 Background

5 1 Background The Equation of Heat Conduction: Initial Condition: Boundary Condition: x t T T T T( x,, t T ( x m T( x,, t Tmerf ( ( x 2 t T( x,, t = T( x,, t Tm x Automatically Satisfied

6 2.1 Green's Function The Green s Function can be treated as the temperature at (x,y, at the time t, due to an instantaneous point source of strength unity generated at the point P(x,y, at time t. We solve G and make a Convolution with the Initial Condition T, we will get the solution. 2 2 G ( G G ( x x 2 2 ( ( t t t x The Green s Function also satisfy the same B.C. as T does. G( x,, t = To satisfy this boundary condition, we can make an analogy with what we ve learned in Electromagnetics, which means that we need a Imaginary source that absorbs heat at the mirror point of the original heat source.

7 2.1 Green's Function -q -Q U= x T= x E Heat Flow q Q The equilibrium function of U and T are all Laplace s Equation. To satisfy T= at =, what we need is to make the final solution of the Green s Function to be G(- -G(+ (since the upper source is negative.

8 2.2 Solution to the Green s Function 2 2 G ( G G ( x x ( ( t t 2 2 t x Take a Fourier Transform to Green s Function. G i ( kxxk 4 ( k x k G+e ( t t t Apply the formula for the solution of First Order ODE. For y' P( x y Q( x, the solution will be: So, G x x y e Q( x e dx t t x x t 2 i( kxx k t 4 ( k k dt 4 ( k k dt e e ( t t e dt e x P( t dt P( t dt kx k t t i kxx k e 4 ( ( 2 ( x

9 2.2 Solution to the Green s Function Take a Inverse Fourier Transform to the solution ( kx k ( tt 2 i( kxx k 2 i( kxxk G e e e dkxdk ( t t i( xx i( [2 ( tt k ] [2 ( tt k ] e 2 2 x 4 ( tt 4 ( tt 2 2 ( xx ( 4 ( tt e d(2 ( t t k d(2 ( t t x k 1 4 ( t t e 2 2 ( xx ( 4 ( tt Finally, ( xx ( ( xx ( 4 ( tt 4 ( tt ( 1 G e e 4 ( t t

10 2.3 Solution to Temperature Evolution The Initial Condition of T is T '( x,, t (1 H( x Tm TmH ( x erf ( 2 t The Temperature will be the Convolution of T and G T ( x,, t T '( I. C.* G T '( x,, t G( x x,, t t dx d ( TmGdx Tmerf ( Gdx d 2 t x x 1, 2, 3 4 ( t t 4 ( t t 4 ( t t

11 2.3 Solution to Temperature Evolution 2 2 ( xx ( 2 2 ( xx ( ( 4 ( tt 4 ( tt 1 T Gdx d T e e dx d m m 4 ( t t T m 4 T 4 T 4 m m 2 ( 2 ( erfc( x 2 4 ( tt 4 ( tt x 2 4 ( tt 4 ( tt e e d d e e d d x erfc( 4 ( t t 4 ( t t Tm x erfc( erfc( 4 4 ( t t 4 ( t t Tm x erfc( erf ( 2 4 ( t t 4 ( t t 1 3 Solution to the first part

12 2.3 Solution to Temperature Evolution T erf ( Gdx d m 2 t 2 2 ( ( x ( tt 4 ( tt ( ( ( Tm erfc 4 e e erf d 4 ( t t 4 ( t t 2 t Consider the first part: erf ( * G 2 t = F F erf G 2 ( 4 ( tt e erf ( d 4 ( t t 2 t [ ( ( ] 2 t Fourier Transform just come to my mind

13 2.3 Solution to Temperature Evolution F F F 2 t 2ik 2 t 2ik 2 ( erf ( ( erf '( ( ( 2 t e d / F ( 2ik t 2ik t 2 2 ( ( 2 t 2 t 2ik e e e d 1 2 t 2 ik t 2 t tk e i k 1 erf G e e 2 t i k ( 2 i tk 4 tk 2 t e d i t k tk t t k ( ( * F [ ] F 1 i k 2 2 tk 1 4 [ e ] erf ( 2 t ( 2

14 2.3 Solution to Temperature Evolution Source Imginary Source 2 tt ( 4 ( 2 1 e erf ( d erf ( 4 ( t t 2 t 2 t 2 tt ( 4 ( 2 1 e erf ( d erf ( 4 ( t t 2 t 2 t T erf ( Gdx d m 2 t 2 2 ( ( x ( tt 4 ( tt ( ( ( Tm erfc e e erf d 4 4 ( t t 4 ( t t 2 t Tm x erfc( erf ( 2 4 ( t t 2 t Solution to the second part

15 2.3 Solution to Temperature Evolution T ( x,, t ( T Gdx T erf ( Gdx d m m 2 t Tm x [ erfc( erf ( 2 4 ( t t 4 ( t t x erfc( erf ( ] 4 ( t t 2 t

16 3.1 Bathymetry: Solution of Model Topography 2D model for Ocean Floor Topography by Isostasy. (P175 w 2 ( T T t m v m m w 1/2 3D model for Ocean Floor Topography by Isostasy: Plug Temperature to the model and do the Integral on P175 1/2 m vtm t x t x ( t t w erfc( erfc( m w 2 t 2 t Arbitrarily defination: v m K, 1 m s, Tm 13 K, kg m kg m t Myr , w 1, 2

17 3.2 Model Temperature and Topography t, t 2 Myr, Temperature

18 3.2 Model Temperature and Topography t 1 Myr, t 2 Myr, Temperature

19 3.2 Model Temperature and Topography t 5 Myr, t 2 Myr, Temperature

20 3.2 the Model Temperature and Topography t 1 Myr, t 2 Myr, Temperature

21 3.2 Model Temperature and Topography t 2 Myr, t 2 Myr, Temperature

22 3.2 Model Temperature and Topography t 3 Myr, t 2 Myr, Temperature

23 3.2 Model Temperature and Topography t, t 2 Myr, Topography

24 3.2 Model Temperature and Topography t 1 Myr, t 2 Myr, Topography

25 3.2 Model Temperature and Topography t 5 Myr, t 2 Myr, Topography

26 3.2 Model Temperature and Topography t 1 Myr, t 2 Myr, Topography

27 3.2 Model Temperature and Topography t 2 Myr, t 2 Myr, Topography

28 3.2 Model Temperature and Topography t 3 Myr, t 2 Myr, Topography

29 3.3 Comparison with real Topography

30 3.3 Comparison with real Topography t 5 Myr, t 28 Myr, Topography Model Reality Long:129W Lat:36-42N Data all set as in the paper

31 3.3 Comparison with real Topography t 15 Myr, t 28 Myr, Topography Model Reality Long:133W Lat:36-42N

32 3.3 Comparison with real Topography t 25 Myr, t 28 Myr, Topography Model Reality Long:137W Lat:36-42N

33 3.3 Comparison with real Topography t 35 Myr, t 28 Myr, Topography Model Reality Long:141W Lat:36-42N

34 3.3 Comparison with real Topography t 45 Myr, t 28 Myr, Topography Model Reality Long:144W Lat:36-42N

35 4 Conclusions and Discussion (1.It works well for this model on Difference, but not quite well on actually depth of the topography. Top of the ridge? Sediments? (2.The model is smoothed but real topography is uneven and have many spikes on it. The model of Isostasy is too ideal. Neglected factors? Solution in the paper based on elastic model can explain the spikes to some extent. (3.About the initial condition. ERF function?

36 Thank You for Your Attention!

An Introduction to Partial Differential Equations

An Introduction to Partial Differential Equations An Introduction to Partial Differential Equations Ryan C. Trinity University Partial Differential Equations Lecture 1 Ordinary differential equations (ODEs) These are equations of the form where: F(x,y,y,y,y,...)

More information

Fracture Zone Flexure

Fracture Zone Flexure November 16, 2011 Outline 1 Ocean Floor Topography Reality and model 2 3 Assumptions Results Conclusion Ocean Floor Topography Ocean Floor Topography Reality and model Figure: From Geodynamics Turcotte

More information

TEMPERATURE EVOLUTION ACROSS AN OCEAN FRACTURE ZONE

TEMPERATURE EVOLUTION ACROSS AN OCEAN FRACTURE ZONE TEMPERATURE EVOLUTION ACROSS AN OCEAN FRACTURE ZONE Christine Chesley 31 Oct 2016 https://www.britannica.com/science/submarine-fracture-zone OUTLINE Posing the Problem Equation and BCs Solution Visual

More information

Differential equations, comprehensive exam topics and sample questions

Differential 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 information

Math Ordinary Differential Equations Sample Test 3 Solutions

Math Ordinary Differential Equations Sample Test 3 Solutions Solve the following Math - Ordinary Differential Equations Sample Test Solutions (i x 2 y xy + 8y y(2 2 y (2 (ii x 2 y + xy + 4y y( 2 y ( (iii x 2 y xy + y y( 2 y ( (i The characteristic equation is m(m

More information

Math 4263 Homework Set 1

Math 4263 Homework Set 1 Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that

More information

PH2130. Week 1. Week 2. Week 3. Questions for contemplation. Why differential equations? Why usually linear diff eq n s? Why usually 2 nd order?

PH2130. Week 1. Week 2. Week 3. Questions for contemplation. Why differential equations? Why usually linear diff eq n s? Why usually 2 nd order? PH130 week 3, page 1 PH130 Questions for contemplation Week 1 Why differential equations? Week Why usually linear diff eq n s? Week 3 Why usually nd order? PH130 week 3, page Aims of Wk 3 Lect 1 Recognise

More information

Jim Lambers MAT 285 Spring Semester Practice Exam 2 Solution. y(t) = 5 2 e t 1 2 e 3t.

Jim Lambers MAT 285 Spring Semester Practice Exam 2 Solution. y(t) = 5 2 e t 1 2 e 3t. . Solve the initial value problem which factors into Jim Lambers MAT 85 Spring Semester 06-7 Practice Exam Solution y + 4y + 3y = 0, y(0) =, y (0) =. λ + 4λ + 3 = 0, (λ + )(λ + 3) = 0. Therefore, the roots

More information

Partial Differential Equations

Partial 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 information

(b) What is the amplitude at the altitude of a satellite of 400 km?

(b) What is the amplitude at the altitude of a satellite of 400 km? Practice final quiz 2015 Geodynamics 2015 1 a) Complete the following table. parameter symbol units temperature T C or K thermal conductivity heat capacity density coefficient of thermal expansion volumetric)

More information

Solutions of differential equations using transforms

Solutions of differential equations using transforms Solutions of differential equations using transforms Process: Take transform of equation and boundary/initial conditions in one variable. Derivatives are turned into multiplication operators. Solve (hopefully

More information

KINEMATICS & DYNAMICS

KINEMATICS & DYNAMICS KINEMATICS & DYNAMICS BY ADVANCED DIFFERENTIAL EQUATIONS Question (**+) In this question take g = 0 ms. A particle of mass M kg is released from rest from a height H m, and allowed to fall down through

More information

Partial Differential Equations

Partial 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 information

e (x y)2 /4kt φ(y) dy, for t > 0. (4)

e (x y)2 /4kt φ(y) dy, for t > 0. (4) Math 24A October 26, 2 Viktor Grigoryan Heat equation: interpretation of the solution Last time we considered the IVP for the heat equation on the whole line { ut ku xx = ( < x

More information

PDEs, part 1: Introduction and elliptic PDEs

PDEs, part 1: Introduction and elliptic PDEs PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,

More information

6. Continuous Release - Point Source

6. Continuous Release - Point Source 1 6. Continuous Release - Point Source A scalar released continuously into a moving fluid and from a discrete point will form a plume that grows in the lateral dimension through diffusion and extends downstream

More information

Solving First Order PDEs

Solving First Order PDEs Solving Ryan C. Trinity University Partial Differential Equations Lecture 2 Solving the transport equation Goal: Determine every function u(x, t) that solves u t +v u x = 0, where v is a fixed constant.

More information

EXAM. Exam #1. Math 3350 Summer II, July 21, 2000 ANSWERS

EXAM. Exam #1. Math 3350 Summer II, July 21, 2000 ANSWERS EXAM Exam #1 Math 3350 Summer II, 2000 July 21, 2000 ANSWERS i 100 pts. Problem 1. 1. In each part, find the general solution of the differential equation. dx = x2 e y We use the following sequence of

More information

MATH 251 Final Examination August 14, 2015 FORM A. Name: Student Number: Section:

MATH 251 Final Examination August 14, 2015 FORM A. Name: Student Number: Section: MATH 251 Final Examination August 14, 2015 FORM A Name: Student Number: Section: This exam has 11 questions for a total of 150 points. Show all your work! In order to obtain full credit for partial credit

More information

Ishik University / Sulaimani Architecture Department. Structure. ARCH 214 Chapter -5- Equilibrium of a Rigid Body

Ishik University / Sulaimani Architecture Department. Structure. ARCH 214 Chapter -5- Equilibrium of a Rigid Body Ishik University / Sulaimani Architecture Department 1 Structure ARCH 214 Chapter -5- Equilibrium of a Rigid Body CHAPTER OBJECTIVES To develop the equations of equilibrium for a rigid body. To introduce

More information

Math 240: Spring-mass Systems

Math 240: Spring-mass Systems Math 240: Spring-mass Systems Ryan Blair University of Pennsylvania Tuesday March 1, 2011 Ryan Blair (U Penn) Math 240: Spring-mass Systems Tuesday March 1, 2011 1 / 15 Outline 1 Review 2 Today s Goals

More information

GEOL5690 Class notes #2: Passive Margins & Thermal Subsidence

GEOL5690 Class notes #2: Passive Margins & Thermal Subsidence GEOL569 Class notes #2: Passive Margins & Thermal Subsidence Reference: Turcotte and Schubert, sec. 4-15, 4-16, 4-23, 4-25, last is most relevant, earlier readings provide background. [Alternative: Stüwe,

More information

Numerical Methods for PDEs

Numerical Methods for PDEs Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh

More information

Differential Equations

Differential Equations Electricity and Magnetism I (P331) M. R. Shepherd October 14, 2008 Differential Equations The purpose of this note is to provide some supplementary background on differential equations. The problems discussed

More information

Partial Differential Equations (PDEs)

Partial Differential Equations (PDEs) C H A P T E R Partial Differential Equations (PDEs) 5 A PDE is an equation that contains one or more partial derivatives of an unknown function that depends on at least two variables. Usually one of these

More information

Sturm-Liouville Theory

Sturm-Liouville Theory More on Ryan C. Trinity University Partial Differential Equations April 19, 2012 Recall: A Sturm-Liouville (S-L) problem consists of A Sturm-Liouville equation on an interval: (p(x)y ) + (q(x) + λr(x))y

More information

MATH 2250 Final Exam Solutions

MATH 2250 Final Exam Solutions MATH 225 Final Exam Solutions Tuesday, April 29, 28, 6: 8:PM Write your name and ID number at the top of this page. Show all your work. You may refer to one double-sided sheet of notes during the exam

More information

Functions of Several Variables

Functions of Several Variables Functions of Several Variables Partial Derivatives Philippe B Laval KSU March 21, 2012 Philippe B Laval (KSU) Functions of Several Variables March 21, 2012 1 / 19 Introduction In this section we extend

More information

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018 Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent

More information

Power Series and Analytic Function

Power Series and Analytic Function Dr Mansoor Alshehri King Saud University MATH204-Differential Equations Center of Excellence in Learning and Teaching 1 / 21 Some Reviews of Power Series Differentiation and Integration of a Power Series

More information

Pithy P o i n t s Picked I ' p and Patljr Put By Our P e r i p a tetic Pencil Pusher VOLUME X X X X. Lee Hi^h School Here Friday Ni^ht

Pithy P o i n t s Picked I ' p and Patljr Put By Our P e r i p a tetic Pencil Pusher VOLUME X X X X. Lee Hi^h School Here Friday Ni^ht G G QQ K K Z z U K z q Z 22 x z - z 97 Z x z j K K 33 G - 72 92 33 3% 98 K 924 4 G G K 2 G x G K 2 z K j x x 2 G Z 22 j K K x q j - K 72 G 43-2 2 G G z G - -G G U q - z q - G x) z q 3 26 7 x Zz - G U-

More information

1.061 / 1.61 Transport Processes in the Environment

1.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. 6.

More information

Geotherms. Reading: Fowler Ch 7. Equilibrium geotherms One layer model

Geotherms. Reading: Fowler Ch 7. Equilibrium geotherms One layer model Geotherms Reading: Fowler Ch 7 Equilibrium geotherms One layer model (a) Standard model: k = 2.5 W m -1 C -1 A = 1.25 x 10-6 W m -3 Q moho = 21 x 10-3 W m -2 shallow T-gradient: 30 C km -1 deep T-gradient:

More information

Generating and characteristic functions. Generating and Characteristic Functions. Probability generating function. Probability generating function

Generating and characteristic functions. Generating and Characteristic Functions. Probability generating function. Probability generating function Generating and characteristic functions Generating and Characteristic Functions September 3, 03 Probability generating function Moment generating function Power series expansion Characteristic function

More information

Math 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

Math 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 information

Math 240: Spring/Mass Systems II

Math 240: Spring/Mass Systems II Math 240: Spring/Mass Systems II Ryan Blair University of Pennsylvania Monday, March 26, 2012 Ryan Blair (U Penn) Math 240: Spring/Mass Systems II Monday, March 26, 2012 1 / 12 Outline 1 Today s Goals

More information

A proof for the full Fourier series on [ π, π] is given here.

A proof for the full Fourier series on [ π, π] is given here. niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at

More information

Ex. 1. Find the general solution for each of the following differential equations:

Ex. 1. Find the general solution for each of the following differential equations: MATH 261.007 Instr. K. Ciesielski Spring 2010 NAME (print): SAMPLE TEST # 2 Solve the following exercises. Show your work. (No credit will be given for an answer with no supporting work shown.) Ex. 1.

More information

Solving First Order PDEs

Solving First Order PDEs Solving Ryan C. Trinity University Partial Differential Equations January 21, 2014 Solving the transport equation Goal: Determine every function u(x, t) that solves u t +v u x = 0, where v is a fixed constant.

More information

Partial Differential Equations Summary

Partial Differential Equations Summary Partial Differential Equations Summary 1. The heat equation Many physical processes are governed by partial differential equations. temperature of a rod. In this chapter, we will examine exactly that.

More information

1 f. result from periodic disturbance same period (frequency) as source Longitudinal or Transverse Waves Characterized by

1 f. result from periodic disturbance same period (frequency) as source Longitudinal or Transverse Waves Characterized by result from periodic disturbance same period (frequency) as source Longitudinal or Transverse Waves Characterized by amplitude (how far do the bits move from their equilibrium positions? Amplitude of MEDIUM)

More information

LECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS

LECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS LECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS 1. Regular Singular Points During the past few lectures, we have been focusing on second order linear ODEs of the form y + p(x)y + q(x)y = g(x). Particularly,

More information

Fundamental Solution

Fundamental Solution Fundamental Solution onsider the following generic equation: Lu(X) = f(x). (1) Here X = (r, t) is the space-time coordinate (if either space or time coordinate is absent, then X t, or X r, respectively);

More information

UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH

UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:

More information

Mathematics 2 for Business Schools Topic 7: Application of Integration to Economics. Building Competence. Crossing Borders.

Mathematics 2 for Business Schools Topic 7: Application of Integration to Economics. Building Competence. Crossing Borders. Mathematics 2 for Business Schools Topic 7: Application of Integration to Economics Building Competence. Crossing Borders. Spring Semester 2017 Learning objectives After finishing this section you should

More information

Topics in Fourier analysis - Lecture 2.

Topics in Fourier analysis - Lecture 2. Topics in Fourier analysis - Lecture 2. Akos Magyar 1 Infinite Fourier series. In this section we develop the basic theory of Fourier series of periodic functions of one variable, but only to the extent

More information

Final Exam May 4, 2016

Final 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 information

Chap 4. State-Space Solutions and

Chap 4. State-Space Solutions and Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations

More information

The Pseudospectral Method

The Pseudospectral Method The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 43 Outline 1 Introduction Motivation History

More information

Computer Vision & Digital Image Processing. Periodicity of the Fourier transform

Computer Vision & Digital Image Processing. Periodicity of the Fourier transform Computer Vision & Digital Image Processing Fourier Transform Properties, the Laplacian, Convolution and Correlation Dr. D. J. Jackson Lecture 9- Periodicity of the Fourier transform The discrete Fourier

More information

Kinematics of fluid motion

Kinematics of fluid motion Chapter 4 Kinematics of fluid motion 4.1 Elementary flow patterns Recall the discussion of flow patterns in Chapter 1. The equations for particle paths in a three-dimensional, steady fluid flow are dx

More information

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered

More information

Elementary ODE Review

Elementary ODE Review Elementary ODE Review First Order ODEs First Order Equations Ordinary differential equations of the fm y F(x, y) () are called first der dinary differential equations. There are a variety of techniques

More information

Final: Solutions Math 118A, Fall 2013

Final: Solutions Math 118A, Fall 2013 Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or

More information

Introduction of Partial Differential Equations and Boundary Value Problems

Introduction 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 information

Math 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv

Math 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium

More information

Mathematical Methods - Lecture 9

Mathematical Methods - Lecture 9 Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr

More information

FOURIER TRANSFORM METHODS David Sandwell, January, 2013

FOURIER TRANSFORM METHODS David Sandwell, January, 2013 1 FOURIER TRANSFORM METHODS David Sandwell, January, 2013 1. Fourier Transforms Fourier analysis is a fundamental tool used in all areas of science and engineering. The fast fourier transform (FFT) algorithm

More information

Thermal and compositional structure of the Mantle and Lithosphere

Thermal and compositional structure of the Mantle and Lithosphere Chapter 1 Thermal and compositional structure of the Mantle and Lithosphere 1.1 Primordial heat of the Earth The most widely accepted planetary formation theory says that the solar system accreted from

More information

3. Continental Heat Flow. Ge 163 4/3/15

3. Continental Heat Flow. Ge 163 4/3/15 3. Continental Heat Flow Ge 163 4/3/15 Outline 1. Measurement of heat flow and Fourier s law 2. Overview of major variations 3. Sclater Histograms 4. Heat-flow - heat-production relationship 5. Continental

More information

Sample Solutions of Assignment 3 for MAT3270B: 2.8,2.3,2.5,2.7

Sample Solutions of Assignment 3 for MAT3270B: 2.8,2.3,2.5,2.7 Sample Solutions of Assignment 3 for MAT327B: 2.8,2.3,2.5,2.7 1. Transform the given initial problem into an equivalent problem with the initial point at the origin (a). dt = t2 + y 2, y(1) = 2, (b). dt

More information

Elastodynamic single-sided homogeneous Green's function representation: Theory and examples

Elastodynamic single-sided homogeneous Green's function representation: Theory and examples Delft University of Technology Elastodynamic single-sided homogeneous Green's function representation: Theory and examples Wapenaar, Kees; Reinicke Urruticoechea, Christian DOI 10.3997/2214-4609.201701220

More information

Darcy s law in 3-D. K * xx K * yy K * zz

Darcy s law in 3-D. K * xx K * yy K * zz PART 7 Equations of flow Darcy s law in 3-D Specific discarge (vector) is calculated by multiplying te ydraulic conductivity (second-order tensor) by te ydraulic gradient (vector). We obtain a general

More information

Wave Equation Modelling Solutions

Wave Equation Modelling Solutions Wave Equation Modelling Solutions SEECS-NUST December 19, 2017 Wave Phenomenon Waves propagate in a pond when we gently touch water in it. Wave Phenomenon Our ear drums are very sensitive to small vibrations

More information

The Heat Equation John K. Hunter February 15, The heat equation on a circle

The Heat Equation John K. Hunter February 15, The heat equation on a circle The Heat Equation John K. Hunter February 15, 007 The heat equation on a circle We consider the diffusion of heat in an insulated circular ring. We let t [0, ) denote time and x T a spatial coordinate

More information

Synthesizing Geometries for 21st Century Electromagnetics

Synthesizing Geometries for 21st Century Electromagnetics ECE 5322 21 st Century Electromagnetics Instructor: Office: Phone: E Mail: Dr. Raymond C. Rumpf A 337 (915) 747 6958 rcrumpf@utep.edu Lecture #18b Synthesizing Geometries for 21st Century Electromagnetics

More information

Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine

Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine Lecture 2 The wave equation Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.0 28/09/2018 1 Learning objectives of this lecture Understand the fundamental properties of the wave equation

More information

The Hopf equation. The Hopf equation A toy model of fluid mechanics

The Hopf equation. The Hopf equation A toy model of fluid mechanics The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van

More information

Anna University May/June 2013 Exams ME2151 Engineering Mechanics Important Questions.

Anna University May/June 2013 Exams ME2151 Engineering Mechanics Important Questions. Anna University May/June 2013 Exams ME2151 Engineering Mechanics Important Questions 1. Find the resultant force and its direction for the given figure 2. Two forces are acting at a point O as shown in

More information

dy dt = 1 y t 1 +t 2 y dy = 1 +t 2 dt 1 2 y2 = 1 2 ln(1 +t2 ) +C, y = ln(1 +t 2 ) + 9.

dy dt = 1 y t 1 +t 2 y dy = 1 +t 2 dt 1 2 y2 = 1 2 ln(1 +t2 ) +C, y = ln(1 +t 2 ) + 9. Math 307A, Winter 2014 Midterm 1 Solutions Page 1 of 8 1. (10 points Solve the following initial value problem explicitly. Your answer should be a function in the form y = g(t, where there is no undetermined

More information

OPAC102. The Acoustic Wave Equation

OPAC102. The Acoustic Wave Equation OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating

More information

Homework for Math , Fall 2016

Homework for Math , Fall 2016 Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant

More information

M445: Heat equation with sources

M445: 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 information

Error functions. Nikolai G. Lehtinen. April 23, erf x = 2 x. e t2 dt (1) π. erf ( ) = 1, erf (+ ) = 1 erf (x ) = [erf (x)]

Error functions. Nikolai G. Lehtinen. April 23, erf x = 2 x. e t2 dt (1) π. erf ( ) = 1, erf (+ ) = 1 erf (x ) = [erf (x)] Error functions Nikolai G. Lehtinen April 3, 1 1 Error function erf x and complementary error function erfc x Gauss) error function is and has properties erf x = x e t dt 1) π erf ) = 1, erf + ) = 1 erf

More information

Solutions of Spring 2008 Final Exam

Solutions of Spring 2008 Final Exam Solutions of Spring 008 Final Exam 1. (a) The isocline for slope 0 is the pair of straight lines y = ±x. The direction field along these lines is flat. The isocline for slope is the hyperbola on the left

More information

kx m x B N 1 C L, M Mg θ

kx m x B N 1 C L, M Mg θ .004 MODELING DYNAMICS AND CONTROL II Spring 00 Solutions to Problem Set No. 7 Problem 1. Pendulum mounted on elastic support. This problem is an execise in the application of momentum principles. Two

More information

Computational Neuroscience. Session 1-2

Computational Neuroscience. Session 1-2 Computational Neuroscience. Session 1-2 Dr. Marco A Roque Sol 05/29/2018 Definitions Differential Equations A differential equation is any equation which contains derivatives, either ordinary or partial

More information

Lecture 9: Kinetic Energy and Work 1

Lecture 9: Kinetic Energy and Work 1 Lecture 9: Kinetic Energy and Work 1 CHAPTER 6: Work and Kinetic Energy The concept of WORK has a very precise definition in physics. Work is a physical quantity produced when a Force moves an object through

More information

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012 MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to

More information

Series Solution of Linear Ordinary Differential Equations

Series Solution of Linear Ordinary Differential Equations Series Solution of Linear Ordinary Differential Equations Department of Mathematics IIT Guwahati Aim: To study methods for determining series expansions for solutions to linear ODE with variable coefficients.

More information

1 Model for gas-liquid reactions Based on film model

1 Model for gas-liquid reactions Based on film model 1 Model for gas-liquid reactions Based on film model A(g >l) +νb(l) Products Governing equations are as follows: d 2 C A dx 2 = k 2 C A C B (1) D B d 2 C B dx 2 = νk 2 C A C B (2) Here x is the actual

More information

A SPECIAL SOLUTION OF CONSTANT COEFFICIENTS PARTIAL DERIVATIVE EQUATIONS WITH FOURIER TRANSFORM METHOD

A SPECIAL SOLUTION OF CONSTANT COEFFICIENTS PARTIAL DERIVATIVE EQUATIONS WITH FOURIER TRANSFORM METHOD Advanced Mathematical Models & Applications Vol.3, No., 208, pp.85-93 A SPECIAL SOLUTION OF CONSTANT COEFFICIENTS PARTIAL DERIVATIVE EQUATIONS WITH FOURIER TRANSFORM METHOD Murat Düz * Department of Mathematics

More information

G: Uniform Convergence of Fourier Series

G: Uniform Convergence of Fourier Series G: Uniform Convergence of Fourier Series From previous work on the prototypical problem (and other problems) u t = Du xx 0 < x < l, t > 0 u(0, t) = 0 = u(l, t) t > 0 u(x, 0) = f(x) 0 < x < l () we developed

More information

Do not turn over until you are told to do so by the Invigilator.

Do not turn over until you are told to do so by the Invigilator. UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2013 14 CALCULUS AND MULTIVARIABLE CALCULUS MTHA4005Y Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.

More information

Chapter 15. Mechanical Waves

Chapter 15. Mechanical Waves Chapter 15 Mechanical Waves A wave is any disturbance from an equilibrium condition, which travels or propagates with time from one region of space to another. A harmonic wave is a periodic wave in which

More information

20D - Homework Assignment 4

20D - Homework Assignment 4 Brian Bowers (TA for Hui Sun) MATH 0D Homework Assignment November, 03 0D - Homework Assignment First, I will give a brief overview of how to use variation of parameters. () Ensure that the differential

More information

Section 5.6 Integration by Parts

Section 5.6 Integration by Parts .. 98 Section.6 Integration by Parts Integration by parts is another technique that we can use to integrate problems. Typically, we save integration by parts as a last resort when substitution will not

More information

Equations of motion for the Pendulum and Augmented-Reality Pendulum sketches

Equations of motion for the Pendulum and Augmented-Reality Pendulum sketches Equations of motion for the Pendulum and Augmented-Reality Pendulum sketches Ludovico Carbone Issue: Date: February 3, 20 School of Physics and Astronomy University of Birmingham Birmingham, B5 2TT Introduction

More information

Scientific Computing: An Introductory Survey

Scientific 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 information

to the potential V to get V + V 0 0Ψ. Let Ψ ( x,t ) =ψ x dx 2

to the potential V to get V + V 0 0Ψ. Let Ψ ( x,t ) =ψ x dx 2 Physics 0 Homework # Spring 017 Due Wednesday, 4/1/17 1. Griffith s 1.8 We start with by adding V 0 to the potential V to get V + V 0. The Schrödinger equation reads: i! dψ dt =! d Ψ m dx + VΨ + V 0Ψ.

More information

Find the value of λ. (Total 9 marks)

Find the value of λ. (Total 9 marks) 1. A particle of mass 0.5 kg is attached to one end of a light elastic spring of natural length 0.9 m and modulus of elasticity λ newtons. The other end of the spring is attached to a fixed point O 3 on

More information

Math 201 Assignment #11

Math 201 Assignment #11 Math 21 Assignment #11 Problem 1 (1.5 2) Find a formal solution to the given initial-boundary value problem. = 2 u x, < x < π, t > 2 u(, t) = u(π, t) =, t > u(x, ) = x 2, < x < π Problem 2 (1.5 5) Find

More information

6 Non-homogeneous Heat Problems

6 Non-homogeneous Heat Problems 6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems

More information

10550 PRACTICE FINAL EXAM SOLUTIONS. x 2 4. x 2 x 2 5x +6 = lim x +2. x 2 x 3 = 4 1 = 4.

10550 PRACTICE FINAL EXAM SOLUTIONS. x 2 4. x 2 x 2 5x +6 = lim x +2. x 2 x 3 = 4 1 = 4. 55 PRACTICE FINAL EXAM SOLUTIONS. First notice that x 2 4 x 2x + 2 x 2 5x +6 x 2x. This function is undefined at x 2. Since, in the it as x 2, we only care about what happens near x 2 an for x less than

More information

Part A: From (14.19) in your text (see notes Stability1 ), we have

Part A: From (14.19) in your text (see notes Stability1 ), we have Prob14.4 14.4 Consider a round-rotor generator delivering a steady-state power P G =.5 to an infinite bus through a transmission line with reactance X L =.4. Assume that E a =1.8, V =1/_, H=5sec, and X

More information

IV. Transport Phenomena Lecture 18: Forced Convection in Fuel Cells II

IV. 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 information

Partial Differential Equations

Partial Differential Equations M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1

More information

DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt

More information

Harmonic Oscillator. Mass-Spring Oscillator Resonance The Pendulum. Physics 109 Experiment Number 12

Harmonic Oscillator. Mass-Spring Oscillator Resonance The Pendulum. Physics 109 Experiment Number 12 Harmonic Oscillator Mass-Spring Oscillator Resonance The Pendulum Physics 109 Experiment Number 12 Outline Simple harmonic motion The vertical mass-spring system Driven oscillations and resonance The pendulum

More information

Slope Fields, Differential Equations... and Maybe Euler s Method, Too

Slope Fields, Differential Equations... and Maybe Euler s Method, Too Slope Fields, Differential Equations... and Maybe Euler s Method, Too Chuck Garner, Ph.D. Department of Mathematics Rockdale Magnet School for Science and Technology October 10, 2014 / Georgia DOE AP Content

More information