Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

[Mathematica] Solving Heat Equation in Spherical Coordinates

  1. Jul 28, 2012 #1
    Hello Folks,

    I have this equation to solve (expressed in LaTeX):

    \frac{\partial{h}}{\partial t} = \frac{1}{n} \left[ \frac{1}{r^2 \sin^2{\phi}} \frac{\partial}{\partial \theta} \left( K \frac{\partial h}{\partial \theta} \right) + \frac{1}{r^2 \sin \phi} \frac{\partial}{\partial \phi} \left( K \sin \phi \frac{\partial h}{\partial \phi}\right) + s(\theta,\phi,t) \right]

    This is similar to heat equation expressed in spherical coordinates, using mathematical convention for \phi and \theta and where s is a source term (but comes from data and do not need to be computed), and n is constant (does not depend on time) and again this is something we know (or assume), and finally, as you can read, there is no gradient in the radial direction.

    I'd like to use Mathematica as I need to resolve this equation numerically for different dataset (for K,n and s, that I have) on a sphere, and it seems to be THE tool for such task. But I don't have any experience (a very a little actually) with Mathematica.

    Could someone help me with this?

    thanks a lot!

    Last edited: Jul 28, 2012
  2. jcsd
  3. Jul 28, 2012 #2
    I might have something close to what I need, but in Maple:

    PDE := n*(diff(h(theta, phi, t), t)) = (diff(K*(diff(h(theta, phi, t),
    theta)), theta))/(r^2*sin(phi)^2)+(diff(K*sin(phi)*(diff(h(theta, phi,
    t), phi)), phi))/(r^2*sin(phi)) ;

    sol := pdsolve(PDE, h(theta, phi, t));

    but I really want to do it with Mathematica, as my Maple license will expire soon and only Mathematica will be available in my place.

    Any ideas how I can convert this into Mathematica?

Share this great discussion with others via Reddit, Google+, Twitter, or Facebook