Group Abstract Group Abstract

Message Boards Message Boards

Solving the Laplace equation with only Neumann boundary conditions

Posted 11 years ago
POSTED BY: Jens Rix
3 Replies
Posted 11 years ago
POSTED BY: Jens Rix

I want to solve the Laplace equation with Neumann boundary conditions on all boundaries, The solution seems random when I do not include any DirichletConditions though.

You can't have Neumann boundary conditions on all boundaries. This is ill-posed problem as there is no unique solution in this case. One of the boundaries must be Dirichlet.

POSTED BY: Nasser M. Abbasi

First, what does it mean to say

NeumannValue[r, z == 0]

taken into account that r is an independent variable?

Then there are not enough conditions to avoid randomness in a way. To see that take manual's example

NDSolveValue[{\!\(
\*SubsuperscriptBox[\(\[Del]\), \({x, y}\), \(2\)]\(u[x, y]\)\) == 
   NeumannValue[1., x >= 0.35], 
  DirichletCondition[u[x, y] == 0., x <= -0.3]}, u, {x, y} \[Element] 
  Disk[]]
Plot3D[%[x, y], {x, y} \[Element] Disk[]]

and kick the Dirichlet condition out. The result is (in 2D)

In[9]:= NDSolveValue[{\!\(
\*SubsuperscriptBox[\(\[Del]\), \({x, y}\), \(2\)]\(u[x, y]\)\) == 
   NeumannValue[1., x >= -1.]}, u, {x, y} \[Element] Disk[]]

During evaluation of In[9]:= NDSolveValue::femibcnd: No DirichletCondition or Robin-type NeumannValue was specified; the result may be off by a constant value. >>

Out[9]= InterpolatingFunction[{{-1., 1.}, {-1., 1.}}, <>]

and this NDSolveValue::femibcnd is pretty reasonable. In 3D the question is: Where are the conditions for the shank?

POSTED BY: Udo Krause
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard