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

I Computing boundary values

  1. Aug 10, 2016 #1
    Hi, I was recently following an example shown in this link and just had a couple questions:
    http://www.scientificpython.net/pyb...e-equation-and-making-a-video-of-the-solution

    I believe I understand the steps, but was just not quite understanding the justification. In the link above, to compute values at i = 0 and/or j = 0, the terms corresponding to u at i-1 and j-1 disappear and the terms at i+1 and j+1, respectively, are doubled. The same thing is done at the end of the boundary where i = n and/or j = m; values for i+1 and j+1, respectively, are removed and the existing terms doubled. I was just wondering why exactly the author did this and the error involved in doing so? These values simply don't exist on the grid as they are outside the domain, but are there any alternative approaches possible?

    Also, would you happen to have any suggestions for better methods than the central difference scheme shown in the link for solving a 2-dimensional (or higher) wave equation?
     
  2. jcsd
  3. Aug 10, 2016 #2
    It is in order to fulfill the boundary conditions. Using this trick the slope "out" of the lattice becomes zero.

    Other methods of solution could be using Chebyshev polynomials.
     
  4. Aug 10, 2016 #3
    I'm not terribly familiar with Chebyshev polynomials, but it seems like there's quite a bit of literature on this very problem (and associated tricks). Thank you!
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted



Similar Discussions: Computing boundary values
  1. Boundary value problem (Replies: 1)

  2. Boundary Value Problem (Replies: 9)

  3. Boundary Value Problem (Replies: 7)

Loading...