Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: measured boundary conditions with pde toolbox
Replies: 17   Last Post: Apr 8, 2014 6:00 AM

 Messages: [ Previous | Next ]
 Bruno Luong Posts: 9,822 Registered: 7/26/08
Re: measured boundary conditions with pde toolbox
Posted: Aug 10, 2009 3:54 PM

"Doug " <dhk@umd.edu> wrote in message <h5psrj\$8g9\$1@fred.mathworks.com>...
> OK, I've answered my own question, and the answer is a boundary M-file.
>
> Though wbound makes a boundary M-file that is based on a simple formula, in exactly the same format as a boundary condition matrix, it's possible to write your own boundary M-file *any way you like*. As long as it takes (p,e,u,time) as inputs and returns [q,g,h,r] as outputs, it works fine with assempde. You just pass the filename of the boundary M-file as the first argument to assempde. So I've written a boundary M-file that takes the necessary inputs (p,e,u,time) but ignores them, instead constructing [q,g,h,r] from my measurements. And it works. Reading the help file for pdebound over and over was the key!
>
> Bruno, I'm still not following your argument that Neumann boundary conditions can't work. I understand that they leave room for an arbitrary constant of integration, but that's OK for me, because in my application, the solution to Poisson's equation is a stream function. My real interest is its gradient, for which the constant of integration is irrelevant. You seem to be concerned about the fact that the nodes are discrete points, but so are my measurements, and so is every numerical grid; again I don't understand. In case we're still talking about two different things, here's what I mean by Neumann boundary conditions:
> http://en.wikipedia.org/wiki/Neumann_boundary_condition
> http://mathworld.wolfram.com/BoundaryConditions.html

You never get the right convergence property with *point-wise* Neumann condition, i.e., your PDE is ill posed; from that you might get all nasty stuff, and the numerical solution might be very far from the true solution, whatever it is.

But nevermind, nobody never get the sense why we bother to define all these Sobolev space anyway, it is not just for fun.

Bruno

Date Subject Author
8/5/09 Doug
8/5/09 Bruno Luong
8/5/09 Doug
8/6/09 Bruno Luong
8/10/09 Doug
8/10/09 Bruno Luong
8/11/09 Bruno Luong
11/17/12 Tony Kittler
11/17/12 Tony Kittler
11/17/12 Tony Kittler
11/18/12 Tony Kittler
11/18/12 Tony Kittler
11/18/12 Tony Kittler
4/2/14 Dmitry
4/2/14 Bill Greene
4/3/14 Dmitry
4/7/14 Dmitry
4/8/14 Dmitry