Get Help

Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Community
- :
- PTC Mathcad
- :
- PTC Mathcad
- :
- Neumann Boundary Conditions

Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Highlighted

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

05-15-2014
05:11 PM

05-15-2014
05:11 PM

Hello folks,

I'm trying to learn how to specify Neumann boundary conditions using **numol** to solve a PDE.

As you may be aware, Prime 3.0 does not feature **Pdesolve** so if we want to solve a PDE, we need to use either **numol**, **relax** or **multigrid**.

I'm using a rather simple example of Pdesolve in which we have a flat plate perfectly insulated at x=0 and with convection at x=L.

Now, the problem I'm having is that I cannot decipher how to express these boundaries conditions when using **numol**.

If someone knows how to do this, please solve the problem in the attached file using numol and post the solution in the reply.

Thanks a lot for you help!

**P.S.: I saved the file in Mathcad 11 format so that a bigger audience may be able to open it and contribute. Thanks again.**

Solved! Go to Solution.

Labels:

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

05-15-2014
09:02 PM

6 REPLIES 6

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

05-15-2014
07:28 PM

05-15-2014
07:28 PM

Re: Neumann Boundary Conditions

I looked at your pdesolve example and your second boundary condition is called "mixed", not Neumann, because it involves both the dependent variable and its derivative. Numol does not accept mixed boundary conditions. That is why we need pdesolve back in Prime, but hopefully in a vector/matrix form.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

05-15-2014
07:51 PM

05-15-2014
07:51 PM

Re: Neumann Boundary Conditions

Thanks Harvey for your reply.

Now, if we change the second boundary condition to a true Neumann BC, say, a heat flux of 5. How do we set up the boundary conditions?

The example that is included in the Mathcad 15 help is, in my opinion, confusing, at least in terms of the boundary conditions.

Thank you again.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

05-15-2014
09:02 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

05-15-2014
09:20 PM

05-15-2014
09:20 PM

Re: Neumann Boundary Conditions

I think I had the init argument wrong, so an updated file is attached. It doesn't change the results because the init values are not a function of x.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

05-15-2014
09:33 PM

05-15-2014
09:33 PM

Re: Neumann Boundary Conditions

Oops. The x and t arguments were still incorrect in places. But because of your problem, nothing changed.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

05-15-2014
09:38 PM

05-15-2014
09:38 PM

Re: Neumann Boundary Conditions

Harvey, thanks a lot!

Yes, I noticed the minor issue with the argument of bc_func, but that much I knew how to figure out. As you mentioned, the result is unaltered as the slab is initially at a constant temperature.

Pdesolve is DEFINITELY needed in Prime as numol is not quite as versatile as pdesolve (not that pdesolve is perfect, but considering...).

Anyway, thanks again.

Best,

Josué

Top Tags