Green's Function for the Heat Equation

Total Page:16

File Type:pdf, Size:1020Kb

Green's Function for the Heat Equation cs: O ani pe ch n e A c M c Hassan, Fluid Mech Open Acc 2017, 4:2 d e i s u s l F Fluid Mechanics: Open Access DOI: 10.4172/2476-2296.1000152 ISSN: 2476-2296 Perspective Open Access Green’s Function for the Heat Equation Abdelgabar Adam Hassan* Department of Mathematics, College of Science and Arts at Tabrjal, AlJouf University, Kingdom of Saudi Arabia Abstract The solution of problem of non-homogeneous partial differential equations was discussed using the joined Fourier- Laplace transform methods in finding the Green’s function of heat equation in different situations. Keywords: Heat equation; Green’s function; Sturm-Liouville So problem; Electrical engineering; Quantum mechanics dy d22 y dp() x dy d y dy d () =+=+() () ()() px px 22px pxbx dx dx dx dx dx dx dx Introduction Thus eqn (3) can be written as: The Green's function is a powerful tool of mathematics method dy is used in solving some linear non-homogenous PDEs, ODEs. So d px() ++ qxy()()λ rxy =0 (4) Green’s functions are derived by the specially development method of dx dx separation of variables, which uses the properties of Dirac’s function. Where q(x)=p(x)c(x) and r(x)=p(x)d(x). This method was considerable more efficient than the others well known classical methods. Equation in form (4) is known as Sturm-Liouville equation. Satisfy the boundary conditions [1-5]. The series solution of differential equation yields an infinite series which often converges slowly. Thus it is difficult to obtain an insight Regular Sturm-Liouville problem: In case p(a)≠0 and p(b)≠0, into over-all behavior of the solution. The Green’s function approach p(x),q(x),r(x) are continuous, the Sturm-Liouville eqn (4) can be would allow us to have an integral representation of the solution expressed as: instead of an infinite series. L[y]=λr(x)y (5) To obtain the filed u, caused by distributed source we calculate the d dy Where L ≡ px() + qx() (6) effect of each elementary portion of source and add (integral) them all. dx dx If G(r, r ) is the field at the observers point r caused by a unit source 0 If the above equation is associated with the following boundary at the source point r , then the field at r caused by distribution f(r ) is 0 0 condition: the integral of f(r0) G(r, r0) over the whole range of r0 occupied by the course. The function G is called Green’s function. α1y(a)+α2y(a)=0 Preliminaries β1y(b)+β2y(b)=0 (7) Where +α ≠0 and β +β ≠0 Sturm-Liouville problem 1 2 1 2 Consider a linear second order differential equation: The eqn (4) and the boundary condition (7) are called regular Sturm-Liouville problem (RSLP). 2 Ax() d y+ Bx() dy ++ Cxy()()λ Dxy =0 (1) dx2 ∂x Singular Sturm-Liouville problem: Consider the equation: Where λ is a parameter to be determined by the boundary L[y]+λr(x)y=0 a<x<b (8) conditions. A(x) is positive continuous function, then by dividing every Where L is defined by eqn (6), p(x) is smooth and r(x)is positive, term by A(x), eqn (1) can be written as: then the Sturm-Liouville problem is called singular if one of the 2 following situations is occurred. d y+() dy ++()()λ = (2) 2 bx cxy dxy 0 dx ∂x (i) If p(a)=0 or p(b)=0 or both Bx() Cx() Dx() Where bx() = , cx() = and dx() = (ii)The interval (a,b) is infinite. Ax() Ax() Ax() Let us define integrating factor p(x) by: x () = ζζ px exp∫ bd() *Corresponding author: Hassan AA, Department of Mathematics, College of a Science and Arts at Tabrjal, AlJouf University, Kingdom of Saudi Arabia, Tel: Multiplying eqn (2) by p(x), we have: +966538043646; E-mail: [email protected] 2 Received March 01, 2016; Accepted March 22, 2017; Published April 02, 2017 () d y+()() dy ++()()()()λ = (3) px2 pxbx pxcxy pxdxy 0 dx ∂x Citation: Hassan AA (2017) Green’s Function for the Heat Equation. Fluid Mech Open Acc 4: 152. doi: 10.4172/2476-2296.1000152 Since: Copyright: © 2017 Hassan AA. This is an open-access article distributed under xx bdζζ bd ζζ the terms of the Creative Commons Attribution License, which permits unrestricted ∫∫ x dp() x ddaa use, distribution, and reproduction in any medium, provided the original author and = e= e∫ b()ζζ d= pxbx()() dx dx dx a source are credited. Fluid Mech Open Acc, an open access journal Volume 4 • Issue 2 • 1000152 ISSN: 2476-2296 Citation: Hassan AA (2017) Green’s Function for the Heat Equation. Fluid Mech Open Acc 4: 152. doi: 10.4172/2476-2296.1000152 Page 2 of 6 Eigenvalue and Eigenfunction: The Eigenvalue from eqn (4) Some important properties of Dirac delta function defining by a Sturm-Liouville operator can be expressed as: Property (1): Symmetry −1 d dy λ = px( ) + qx( ) (9) δ (−x)=δ (x) (18) rx( ) dx dx Proof: Let ζ=−x, then dx=−dζ The non-trival solutions that satisfy the equation and boundary conditions are called eigenfunctions. Therefore the eigenfunction of We can write: the Sturm-Liouville problem from complete sets of orthogonal bases +∞ −∞ +∞ ∫fxdxfdfd( xf)δ(−=−) ∫∫(ζδζ) ( ) ζ= ( ζδζ) ( ) ζ= (0) (19) for the function space is which the weight function is r(x). −∞ +∞ −∞ The Dirac Delta function +∞ But, ∫ f( xf)δ ( x) dx = (0) (20) The delta function is defined as: −∞ Therefore, from eqns (19) and (20), we conclude that δ (−x)=δ (x). 0 x≠ζ δζ( x−=) (10) Property (2): Scaling ∞=x ζ δδ(ax) = 1 ( x) (21) But such that the integral of δ (x−ξ) is normalized to unit; a +∞ Proof: Let ζ=ax, then dx= 1 dζ ∫ δζ( x−=) dx 1 (11) a −∞ If a>0, then, In fact the first operator where Dirac used the delta function is the +∞ +∞ ζζ+∞ ( x)δ( ) = δζ1 ζ= 1 δζ ζ=101ff = (0) integration: ∫∫f ax dx f( ) d ∫ f( ) d −∞ −∞ a a a−∞ a aaa +∞ ∫ f( x)δζ( x− ) dx (12) +∞ 11+∞ 1 −∞ Since, ∫∫f(ζδ) ( xxf) dx= f( ) δ( x) dx = (0) −∞ aa−∞ a Where f(x) is a continuous function, we have to find the value of Therefore: δδ= 1 the integration eqn (12). (ax) a ( x) Since δ (x−ξ) is zero for x≠ξ, the limit of integration may be change Similarly for a<0, δδ(ax) = 1 ( x) , then, to ζ− and ζ+ε, where ε is a small positive number, f(x) is continuous at −a x−ξ, it’s values within the interval (ζ−ε,ζ−) will not different much from We can write: δδ(ax) = 1 ( x) f(ζ), approximately that: a +∞ ζζ++εε Property (3) ∫∫f( xx)δζ( x−=) dx f( ) δζ( x −≈) dx f( ζδζ) ∫( x −) dx (13) −∞ ζζ−εε− δ22−=1 δδ ++ ( −) ( x a) ( xa) xa (22) With the approximation improving as ε approaches zero. 2a From eqn (11), we have: Proof: The argument of this function goes to zero when x=a and x=−a, wherefore: +∞ ζ +ε I = ∫∫δζ( x−=) dx δζ( x −=) dx 1 (14) +∞ +∞ −∞ ζ −ε 2 2 ∫∫f(ζδ) ( xa−=) dx f( x)δ( x+− a)( x a) dx From all values ofε, then by letting ε→0, we can exactly have: −∞ −∞ +∞ ∫ f( xf)δζ( x−=) dx ( ζ) (15) Only at the zero of the argument of the delta function that is: −∞ +∞ −aa +εε+ ∫∫∫f( x)δ( xa2−=2) dx f( x) δδ(xa 22 −+ 22) dx f( x) (xa −) dx (23) Despite the delta function considered as fundamental role in −∞ −aa −εε− electrical engineering and quantum mechanics, but no conventional 2 2 could be found that satisfies eqns (10) and (11), then the delta function Near the two zeros x −a can be approximated as: sought to be view as the limit of the sequence of strongly peaked −2axa( +) ,x →− a 22 = x−=−+ a) ( x ax)( a) function δ (x) such that: ( +2axa( −) ,x →+ a n δδ( xx) = n ( ) limn→∞ (16) In the limit as ε0 the integral eqn (23) becomes: As: +∞ −aa +εε+ 2 2 ∫∫f( x)δδ(x− a) dx = f( x) ( −+22axa)( ) dx + ∫ f( x) δ ( axa)( −) dx δ ( ) = n −∞ −aa −εε− n x 22 −+aaεε+ π + 11 (1 nx) =∫∫f( x)δδ(xa++) dx f( x) (xa−) dx 22aa−−aaεε− 22 +∞ n −nx 1 δ ( ) = e (17) = −+ − n x ∫ f( x) δδ(xa) ( xa) dx π −∞ 2a 2 sin (nx) δ ( x) = ,etc 221 n 2 Therefore: δxa− = δδ( xa −+) ( xa −) nxπ ( ) 2a Fluid Mech Open Acc, an open access journal Volume 4 • Issue 2 • 1000152 ISSN: 2476-2296 Citation: Hassan AA (2017) Green’s Function for the Heat Equation. Fluid Mech Open Acc 4: 152. doi: 10.4172/2476-2296.1000152 Page 3 of 6 Green’s Function for the Heat −−ζτs − expxsa (29) gxs,,ζτ , = Heat equation over infinite or semi-infinite domains ( ) 2as Consider one dimensional heat equation: Taking Laplace transform of eqn (29) we obtain: 2 ∂∂uu2 2 −=−a f xt, (24) 2 ( ) H (t −τ ) −−( x ζ ) ∂∂tx Gxs,,,ζτ = exp ( ) 2 2 (30) aatπτ( − ) aa( t −τ ) Subject to boundary conditions |u(x,t)|<∞ as |x|<∞ and initial condition. Let G(x, t, ζ, τ) be the Green’s function for the one- Example (5.1.1) Consider one dimensional heat equation: dimensional heat equation, then: ∂∂2 uu−−a2 = f( xt, ) ∂∂2 ∂t ∂x2 (31) GG−=2 δ ( xt−ζτ) ∂( −), −∞< x , ζ <∞ ,0 < t , τ (25) a 2 ∂∂tx Subject to boundary condition u(0,t)=0, limuxt( , ) <∞ and initial x→0 Subject to the boundary condition |G(x, t, ζ, τ)|<∞ as |x|<∞, and the condition U(x,0)=0 initial condition G(x, t, ζ, τ)=0.
Recommended publications
  • Laplace's Equation, the Wave Equation and More
    Lecture Notes on PDEs, part II: Laplace's equation, the wave equation and more Fall 2018 Contents 1 The wave equation (introduction)2 1.1 Solution (separation of variables)........................3 1.2 Standing waves..................................4 1.3 Solution (eigenfunction expansion).......................6 2 Laplace's equation8 2.1 Solution in a rectangle..............................9 2.2 Rectangle, with more boundary conditions................... 11 2.2.1 Remarks on separation of variables for Laplace's equation....... 13 3 More solution techniques 14 3.1 Steady states................................... 14 3.2 Moving inhomogeneous BCs into a source................... 17 4 Inhomogeneous boundary conditions 20 4.1 A useful lemma.................................. 20 4.2 PDE with Inhomogeneous BCs: example.................... 21 4.2.1 The new part: equations for the coefficients.............. 21 4.2.2 The rest of the solution.......................... 22 4.3 Theoretical note: What does it mean to be a solution?............ 23 5 Robin boundary conditions 24 5.1 Eigenvalues in nasty cases, graphically..................... 24 5.2 Solving the PDE................................. 26 6 Equations in other geometries 29 6.1 Laplace's equation in a disk........................... 29 7 Appendix 31 7.1 PDEs with source terms; example........................ 31 7.2 Good basis functions for Laplace's equation.................. 33 7.3 Compatibility conditions............................. 34 1 1 The wave equation (introduction) The wave equation is the third of the essential linear PDEs in applied mathematics. In one dimension, it has the form 2 utt = c uxx for u(x; t): As the name suggests, the wave equation describes the propagation of waves, so it is of fundamental importance to many fields.
    [Show full text]
  • Modelling for Aquifer Parameters Determination
    Journal of Engineering Science, Vol. 14, 71–99, 2018 Mathematical Implications of Hydrogeological Conditions: Modelling for Aquifer Parameters Determination Sheriffat Taiwo Okusaga1, Mustapha Adewale Usman1* and Niyi Olaonipekun Adebisi2 1Department of Mathematical Sciences, Olabisi Onabanjo University, Ago-Iwoye, Nigeria 2Department of Earth Sciences, Olabisi Onabanjo University, Ago-Iwoye, Nigeria *Corresponding author: [email protected] Published online: 15 April 2018 To cite this article: Sheriffat Taiwo Okusaga, Mustapha Adewale Usman and Niyi Olaonipekun Adebisi. (2018). Mathematical implications of hydrogeological conditions: Modelling for aquifer parameters determination. Journal of Engineering Science, 14: 71–99, https://doi.org/10.21315/jes2018.14.6. To link to this article: https://doi.org/10.21315/jes2018.14.6 Abstract: Formulations of natural hydrogeological phenomena have been derived, from experimentation and observation of Darcy fundamental law. Mathematical methods are also been applied to expand on these formulations, but the management of complexities which relate to subsurface heterogeneity is yet to be developed into better models. This study employs a thorough approach to modelling flow in an aquifer (a porous medium) with phenomenological parameters such as Transmissibility and Storage Coefficient on the basis of mathematical arguments. An overview of the essential components of mathematical background and a basic working knowledge of groundwater flow is presented. Equations obtained through the modelling process
    [Show full text]
  • THE ONE-DIMENSIONAL HEAT EQUATION. 1. Derivation. Imagine a Dilute Material Species Free to Diffuse Along One Dimension
    THE ONE-DIMENSIONAL HEAT EQUATION. 1. Derivation. Imagine a dilute material species free to diffuse along one dimension; a gas in a cylindrical cavity, for example. To model this mathematically, we consider the concentration of the given species as a function of the linear dimension and time: u(x; t), defined so that the total amount Q of diffusing material contained in a given interval [a; b] is equal to (or at least proportional to) the integral of u: Z b Q[a; b](t) = u(x; t)dx: a Then the conservation law of the species in question says the rate of change of Q[a; b] in time equals the net amount of the species flowing into [a; b] through its boundary points; that is, there exists a function d(x; t), the ‘diffusion rate from right to left', so that: dQ[a; b] = d(b; t) − d(a; t): dt It is an experimental fact about diffusion phenomena (which can be under- stood in terms of collisions of large numbers of particles) that the diffusion rate at any point is proportional to the gradient of concentration, with the right-to-left flow rate being positive if the concentration is increasing from left to right at (x; t), that is: d(x; t) > 0 when ux(x; t) > 0. So as a first approximation, it is natural to set: d(x; t) = kux(x; t); k > 0 (`Fick's law of diffusion’ ): Combining these two assumptions we have, for any bounded interval [a; b]: d Z b u(x; t)dx = k(ux(b; t) − ux(a; t)): dt a Since b is arbitrary, differentiating both sides as a function of b we find: ut(x; t) = kuxx(x; t); the diffusion (or heat) equation in one dimension.
    [Show full text]
  • The Cauchy-Kovalevskaya Theorem
    CHAPTER 2 The Cauchy-Kovalevskaya Theorem This chapter deals with the only “general theorem” which can be extended from the theory of ODEs, the Cauchy-Kovalevskaya Theorem. This also will allow us to introduce the notion of non-characteristic data, principal symbol and the basic clas- sification of PDEs. It will also allow us to explore why analyticity is not the proper regularity for studying PDEs most of the time. Acknowledgements. Some parts of this chapter – in particular the four proofs of the Cauchy-Kovalevskaya theorem for ODEs – are strongly inspired from the nice lecture notes of Bruce Diver available on his webpage, and some other parts – in particular the organisation of the proof in the PDE case and the discussion of the counter-examples – are strongly inspired from the nice lectures notes of Tsogtgerel Gantumur available on his webpage. 1. Recalls on analyticity Let us first do some recalls on the notion of real analyticity. Definition 1.1. A function f defined on some open subset U of the real line is said to be real analytic at a point x if it is an infinitely di↵erentiable function 0 2U such that the Taylor series at the point x0 1 f (n)(x ) (1.1) T (x)= 0 (x x )n n! − 0 n=0 X converges to f(x) for x in a neighborhood of x0. A function f is real analytic on an open set of the real line if it is real analytic U at any point x0 . The set of all real analytic functions on a given open set is often denoted by2C!U( ).
    [Show full text]
  • LECTURE 7: HEAT EQUATION and ENERGY METHODS Readings
    LECTURE 7: HEAT EQUATION AND ENERGY METHODS Readings: • Section 2.3.4: Energy Methods • Convexity (see notes) • Section 2.3.3a: Strong Maximum Principle (pages 57-59) This week we'll discuss more properties of the heat equation, in partic- ular how to apply energy methods to the heat equation. In fact, let's start with energy methods, since they are more fun! Recall the definition of parabolic boundary and parabolic cylinder from last time: Notation: UT = U × (0;T ] ΓT = UT nUT Note: Think of ΓT like a cup: It contains the sides and bottoms, but not the top. And think of UT as water: It contains the inside and the top. Date: Monday, May 11, 2020. 1 2 LECTURE 7: HEAT EQUATION AND ENERGY METHODS Weak Maximum Principle: max u = max u UT ΓT 1. Uniqueness Reading: Section 2.3.4: Energy Methods Fact: There is at most one solution of: ( ut − ∆u = f in UT u = g on ΓT LECTURE 7: HEAT EQUATION AND ENERGY METHODS 3 We will present two proofs of this fact: One using energy methods, and the other one using the maximum principle above Proof: Suppose u and v are solutions, and let w = u − v, then w solves ( wt − ∆w = 0 in UT u = 0 on ΓT Proof using Energy Methods: Now for fixed t, consider the following energy: Z E(t) = (w(x; t))2 dx U Then: Z Z Z 0 IBP 2 E (t) = 2w (wt) = 2w∆w = −2 jDwj ≤ 0 U U U 4 LECTURE 7: HEAT EQUATION AND ENERGY METHODS Therefore E0(t) ≤ 0, so the energy is decreasing, and hence: Z Z (0 ≤)E(t) ≤ E(0) = (w(x; 0))2 dx = 0 = 0 U And hence E(t) = R w2 ≡ 0, which implies w ≡ 0, so u − v ≡ 0, and hence u = v .
    [Show full text]
  • Lecture 24: Laplace's Equation
    Introductory lecture notes on Partial Differential Equations - ⃝c Anthony Peirce. 1 Not to be copied, used, or revised without explicit written permission from the copyright owner. Lecture 24: Laplace's Equation (Compiled 26 April 2019) In this lecture we start our study of Laplace's equation, which represents the steady state of a field that depends on two or more independent variables, which are typically spatial. We demonstrate the decomposition of the inhomogeneous Dirichlet Boundary value problem for the Laplacian on a rectangular domain into a sequence of four boundary value problems each having only one boundary segment that has inhomogeneous boundary conditions and the remainder of the boundary is subject to homogeneous boundary conditions. These latter problems can then be solved by separation of variables. Key Concepts: Laplace's equation; Steady State boundary value problems in two or more dimensions; Linearity; Decomposition of a complex boundary value problem into subproblems Reference Section: Boyce and Di Prima Section 10.8 24 Laplace's Equation 24.1 Summary of the equations we have studied thus far In this course we have studied the solution of the second order linear PDE. @u = α2∆u Heat equation: Parabolic T = α2X2 Dispersion Relation σ = −α2k2 @t @2u (24.1) = c2∆u Wave equation: Hyperbolic T 2 − c2X2 = A Dispersion Relation σ = ick @t2 ∆u = 0 Laplace's equation: Elliptic X2 + Y 2 = A Dispersion Relation σ = k Important: (1) These equations are second order because they have at most 2nd partial derivatives. (2) These equations are all linear so that a linear combination of solutions is again a solution.
    [Show full text]
  • Greens' Function and Dual Integral Equations Method to Solve Heat Equation for Unbounded Plate
    Gen. Math. Notes, Vol. 23, No. 2, August 2014, pp. 108-115 ISSN 2219-7184; Copyright © ICSRS Publication, 2014 www.i-csrs.org Available free online at http://www.geman.in Greens' Function and Dual Integral Equations Method to Solve Heat Equation for Unbounded Plate N.A. Hoshan Tafila Technical University, Tafila, Jordan E-mail: [email protected] (Received: 16-3-13 / Accepted: 14-5-14) Abstract In this paper we present Greens function for solving non-homogeneous two dimensional non-stationary heat equation in axially symmetrical cylindrical coordinates for unbounded plate high with discontinuous boundary conditions. The solution of the given mixed boundary value problem is obtained with the aid of the dual integral equations (DIE), Greens function, Laplace transform (L-transform) , separation of variables and given in form of functional series. Keywords: Greens function, dual integral equations, mixed boundary conditions. 1 Introduction Several problems involving homogeneous mixed boundary value problems of a non- stationary heat equation and Helmholtz equation were discussed with details in monographs [1, 4-7]. In this paper we present the use of Greens function and DIE for solving non- homogeneous and non-stationary heat equation in axially symmetrical cylindrical coordinates for unbounded plate of a high h. The goal of this paper is to extend the applications DIE which is widely used for solving mathematical physics equations of the 109 N.A. Hoshan elliptic and parabolic types in different coordinate systems subject to mixed boundary conditions, to solve non-homogeneous heat conduction equation, related to both first and second mixed boundary conditions acted on the surface of unbounded cylindrical plate.
    [Show full text]
  • Chapter H1: 1.1–1.2 – Introduction, Heat Equation
    Name Math 3150 Problems Haberman Chapter H1 Due Date: Problems are collected on Wednesday. Submitted work. Please submit one stapled package per problem set. Label each problem with its corresponding problem number, e.g., Problem H3.1-4 . Labels like Problem XC-H1.2-4 are extra credit problems. Attach this printed sheet to simplify your work. Labels Explained. The label Problem Hx.y-z means the problem is for Haberman 4E or 5E, chapter x , section y , problem z . Label Problem H2.0-3 is background problem 3 in Chapter 2 (take y = 0 in label Problem Hx.y-z ). If not a background problem, then the problem number should match a corresponding problem in the textbook. The chapters corresponding to 1,2,3,4,10 in Haberman 4E or 5E appear in the hybrid textbook Edwards-Penney-Haberman as chapters 12,13,14,15,16. Required Background. Only calculus I and II plus engineering differential equations are assumed. If you cannot understand a problem, then read the references and return to the problem afterwards. Your job is to understand the Heat Equation. Chapter H1: 1.1{1.2 { Introduction, Heat Equation Problem H1.0-1. (Partial Derivatives) Compute the partial derivatives. @ @ @ (1 + t + x2) (1 + xt + (tx)2) (sin xt + t2 + tx2) @x @x @x @ @ @ (sinh t sin 2x) (e2t+x sin(t + x)) (e3t+2x(et sin x + ex cos t)) @t @t @t Problem H1.0-2. (Jacobian) f1 Find the Jacobian matrix J and the Jacobian determinant jJj for the transformation F(x; y) = , where f1(x; y) = f2 y cos(2y) sin(3x), f2(x; y) = e sin(y + x).
    [Show full text]
  • 10 Heat Equation: Interpretation of the Solution
    Math 124A { October 26, 2011 «Viktor Grigoryan 10 Heat equation: interpretation of the solution Last time we considered the IVP for the heat equation on the whole line u − ku = 0 (−∞ < x < 1; 0 < t < 1); t xx (1) u(x; 0) = φ(x); and derived the solution formula Z 1 u(x; t) = S(x − y; t)φ(y) dy; for t > 0; (2) −∞ where S(x; t) is the heat kernel, 1 2 S(x; t) = p e−x =4kt: (3) 4πkt Substituting this expression into (2), we can rewrite the solution as 1 1 Z 2 u(x; t) = p e−(x−y) =4ktφ(y) dy; for t > 0: (4) 4πkt −∞ Recall that to derive the solution formula we first considered the heat IVP with the following particular initial data 1; x > 0; Q(x; 0) = H(x) = (5) 0; x < 0: Then using dilation invariance of the Heaviside step function H(x), and the uniquenessp of solutions to the heat IVP on the whole line, we deduced that Q depends only on the ratio x= t, which lead to a reduction of the heat equation to an ODE. Solving the ODE and checking the initial condition (5), we arrived at the following explicit solution p x= 4kt 1 1 Z 2 Q(x; t) = + p e−p dp; for t > 0: (6) 2 π 0 The heat kernel S(x; t) was then defined as the spatial derivative of this particular solution Q(x; t), i.e. @Q S(x; t) = (x; t); (7) @x and hence it also solves the heat equation by the differentiation property.
    [Show full text]
  • 9 the Heat Or Diffusion Equation
    9 The heat or diffusion equation In this lecture I will show how the heat equation 2 2 ut = α ∆u; α 2 R; (9.1) where ∆ is the Laplace operator, naturally appears macroscopically, as the consequence of the con- servation of energy and Fourier's law. Fourier's law also explains the physical meaning of various boundary conditions. I will also give a microscopic derivation of the heat equation, as the limit of a simple random walk, thus explaining its second title | the diffusion equation. 9.1 Conservation of energy plus Fourier's law imply the heat equation In one of the first lectures I deduced the fundamental conservation law in the form ut + qx = 0 which connects the quantity u and its flux q. Here I first generalize this equality for more than one spatial dimension. Let e(t; x) denote the thermal energy at time t at the point x 2 Rk, where k = 1; 2; or 3 (straight line, plane, or usual three dimensional space). Note that I use bold fond to emphasize that x is a vector. The law of the conservation of energy tells me that the rate of change of the thermal energy in some domain D is equal to the flux of the energy inside D minus the flux of the energy outside of D and plus the amount of energy generated in D. So in the following I will use D to denote my domain in R1; R2; or R3. Can it be an arbitrary domain? Not really, and for the following to hold I assume that D is a domain without holes (this is called simply connected) and with a piecewise smooth boundary, which I will denote @D.
    [Show full text]
  • Review of the Methods of Transition from Partial to Ordinary Differential
    Archives of Computational Methods in Engineering https://doi.org/10.1007/s11831-021-09550-5 REVIEW ARTICLE Review of the Methods of Transition from Partial to Ordinary Diferential Equations: From Macro‑ to Nano‑structural Dynamics J. Awrejcewicz1 · V. A. Krysko‑Jr.1 · L. A. Kalutsky2 · M. V. Zhigalov2 · V. A. Krysko2 Received: 23 October 2020 / Accepted: 10 January 2021 © The Author(s) 2021 Abstract This review/research paper deals with the reduction of nonlinear partial diferential equations governing the dynamic behavior of structural mechanical members with emphasis put on theoretical aspects of the applied methods and signal processing. Owing to the rapid development of technology, materials science and in particular micro/nano mechanical systems, there is a need not only to revise approaches to mathematical modeling of structural nonlinear vibrations, but also to choose/pro- pose novel (extended) theoretically based methods and hence, motivating development of numerical algorithms, to get the authentic, reliable, validated and accurate solutions to complex mathematical models derived (nonlinear PDEs). The review introduces the reader to traditional approaches with a broad spectrum of the Fourier-type methods, Galerkin-type methods, Kantorovich–Vlasov methods, variational methods, variational iteration methods, as well as the methods of Vaindiner and Agranovskii–Baglai–Smirnov. While some of them are well known and applied by computational and engineering-oriented community, attention is paid to important (from our point of view) but not widely known and used classical approaches. In addition, the considerations are supported by the most popular and frequently employed algorithms and direct numerical schemes based on the fnite element method (FEM) and fnite diference method (FDM) to validate results obtained.
    [Show full text]
  • Method of Separation of Variables
    2.2. Linearity 33 2.2 Linearity As in the study of ordinary differential equations, the concept of linearity will be very important for liS. A linear operator L lJy definition satisfies Chapter 2 (2.2.1) for any two functions '/1,1 and 112, where Cl and ('2 are arhitrary constants. a/at and a 2 /ax' are examples of linear operators since they satisfy (2.2.1): Method of Separation iJ at (el'/11 + C2'/1,) of Variables a' 2 (C1'/11 +C2U 2) ax It can be shown (see Exercise 2.2.1) that any linear combination of linear operators is a linear operator. Thus, the heat operator 2.1 Introduction a iJl In Chapter I we developed from physical principles an lmderstanding of the heat at - k a2·2 eqnation and its corresponding initial and boundary conditions. We are ready is also a linear operator. to pursue the mathematical solution of some typical problems involving partial A linear equation fora it is of the form differential equations. \Ve \-vilt use a technique called the method of separation of variables. You will have to become an expert in this method, and so we will discuss L(u) = J, (2.2.2) quite a fev.; examples. v~,fe will emphasize problem solving techniques, but \ve must also understand how not to misuse the technique. where L is a linear operator and f is known. Examples of linear partial dijjerentinl A relatively simple but typical, problem for the equation of heat conduction equations are occurs for a one - dimensional rod (0 ::; x .::; L) when all the thermal coefficients are constant.
    [Show full text]