{"id":142,"date":"2014-10-22T15:20:43","date_gmt":"2014-10-22T20:20:43","guid":{"rendered":"http:\/\/charles-oneill.com\/blog\/?p=142"},"modified":"2014-10-22T15:32:59","modified_gmt":"2014-10-22T20:32:59","slug":"solving-a-convection-pde","status":"publish","type":"post","link":"https:\/\/charles-oneill.com\/blog\/solving-a-convection-pde\/","title":{"rendered":"Solving a Convection PDE"},"content":{"rendered":"<p>The problem under consideration is a linear convection PDE.<\/p>\n<p>\\(\\) From Farlow&#8217;s &#8220;Partial Differential Equations for Scientists and Engineers&#8221;, chapter 15 problem 3 is $$u_t = -2u_x \\quad -\\infty &lt; x &lt; \\infty \\quad 0&lt;t&lt;\\infty $$$$ u(x,0)=e^{-x^2}$$<\/p>\n<p><strong>Intuition<\/strong><\/p>\n<p>We see that the governing equation is a linear convection problem. The characteristic velocity is 2. We expect the solution to be a shifted initial condition $$u(x,t)=e^{-(x-2t)^{2}}$$<\/p>\n<p><!--more--><\/p>\n<p><strong>Laplace<\/strong><\/p>\n<p>The Laplace transform of the governing equation with respect to the temporal coordinate t\u00a0 is $$L[u_{t}]=L[-2u_{x}]$$<br \/>\nExpanding with definitions gives $$U(x,s)-u(x,0)=-2\\frac{dU(x,s)}{dx}$$<br \/>\nThe transform of the IC is not found in Farlow&#8217;s tables $$L[u(x,0)]\u00a0\u00a0 \u00a0=\u00a0\u00a0 \u00a0L[e^{-x^{2}}]\u00a0\u00a0 =\u00a0\u00a0 \u00a0?$$<br \/>\nThe Laplace ODE looks like the non-homogeneous equation $$\\frac{dU(x,s)}{dx}+\\frac{1}{2}sU(x,s)=\\frac{1}{2}u(x,0)$$<br \/>\nTwo difficult steps appear before we even see the inverse Laplace transform.<br \/>\nAt this point, you should bail out and try another approach.<br \/>\n<strong>Fourier<\/strong><\/p>\n<p>The Fourier transform of the governing equation with respect to spatial coordinate x is $$F(u_{t})=F(-2u_{x})$$<br \/>\nExpanding with identities gives $$\\frac{dU(t)}{dt}=-2i\\zeta U(t)$$<br \/>\nThe solution to U(t) , being a 1st order ODE, is $$U(t)\u00a0\u00a0 \u00a0=\u00a0\u00a0 \u00a0U(0)e^{-2i\\zeta t}$$<br \/>\nConvolution $$u(x,t)\u00a0\u00a0 \u00a0=\u00a0\u00a0 \u00a0F^{-1}(U(0)e^{-2i\\zeta t})<br \/>\n=\u00a0\u00a0 \u00a0e^{-x^{2}}\\star F^{-1}(e^{-2i\\zeta t})$$<br \/>\nFinding $$e^{-2i\\zeta t}$$<br \/>\nin Farlow&#8217;s table A gives $$f(x-a)$$<br \/>\nwhen $$a=2t$$<\/p>\n<p>The solution is $$u(x,t)\u00a0\u00a0 \u00a0=\u00a0\u00a0 \u00a0e^{-(x-a)^{2}}<br \/>\n=\u00a0\u00a0 \u00a0e^{-(x-2t)^{2}}$$<\/p>\n","protected":false},"excerpt":{"rendered":"<p>The problem under consideration is a linear convection PDE. From Farlow&#8217;s &#8220;Partial Differential Equations for Scientists and Engineers&#8221;, chapter 15 problem 3 is $$u_t = -2u_x \\quad -\\infty &lt; x &lt; \\infty \\quad 0&lt;t&lt;\\infty $$$$ u(x,0)=e^{-x^2}$$ Intuition We see that the governing equation is a linear convection problem. The characteristic velocity is 2. We expect [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[150,11,151,149],"_links":{"self":[{"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/posts\/142"}],"collection":[{"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/comments?post=142"}],"version-history":[{"count":13,"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/posts\/142\/revisions"}],"predecessor-version":[{"id":165,"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/posts\/142\/revisions\/165"}],"wp:attachment":[{"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/media?parent=142"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/categories?post=142"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/charles-oneill.com\/blog\/wp-json\/wp\/v2\/tags?post=142"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}