<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Markeringsbobletekst Tegn";
margin:0cm;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
span.MarkeringsbobletekstTegn
{mso-style-name:"Markeringsbobletekst Tegn";
mso-style-priority:99;
mso-style-link:Markeringsbobletekst;
font-family:"Tahoma","sans-serif";
mso-fareast-language:DA;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:3.0cm 2.0cm 3.0cm 2.0cm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="DA" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hi John<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">The short answer is that from you brief description I cannot see what is wrong.
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">According to Gauss’ theorem the flux through a closed surface equals the sources inside the surface. So I guess you cannot fix both charge and
flux in your problem.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">KR Lars<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">Fra:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> John_V [mailto:jvillar.john@gmail.com]
<br>
<b>Sendt:</b> 6. august 2011 20:12<br>
<b>Til:</b> Lars Rindorf<br>
<b>Cc:</b> getdp@geuz.org<br>
<b>Emne:</b> Re: [Getdp] Von Neumann misconception?<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal" style="margin-bottom:12.0pt">Thanks Lars. If I understand you, you're telling me that my answer A was the right one. (Eimposed was a very simple function in my case. It was constant.) But in my solution, the field at UpperEndCap is much
different from Eimposed. Why is that?<br>
<br>
I am suspecting the real answer is B. This suspicion is, however, based only on observation. I don't know whether this makes sense in terms of how the program is supposed to operate. Domain_Ele does not exclude UpperEndCap. That is, the nodes on UpperEndCap
(a surface) are also nodes of Domain_Ele, which is the entire mesh. Does that offer a way to explain why the field is not what I expected?<br>
<span style="color:#888888"><br>
John</span><o:p></o:p></p>
<div>
<p class="MsoNormal">On Thu, Aug 4, 2011 at 2:23 AM, Lars Rindorf <<a href="mailto:lrf@teknologisk.dk">lrf@teknologisk.dk</a>> wrote:<o:p></o:p></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:11.0pt;color:#1F497D">Hi John</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:11.0pt;color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:11.0pt;color:#1F497D">The Neumann (von Neumann is a different mathematician) boundary condition fixes the ‘influx’ grad(V) normal to the boundary,
and it arises when grad(V) = -E is known and can be replaced by a function, such as Eimposed. A dirichlet boundary condition, e.g. V=0, fixes the magnitude.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:11.0pt;color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:11.0pt;color:#1F497D">Be aware that in the formulation the ‘Eimposed’ is the normal component of the incoming field/flux and it is thus a scalar
not a vector.</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:11.0pt;color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:11.0pt;color:#1F497D">KR Lars</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:11.0pt;color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span style="font-size:10.0pt">Fra:</span></b><span style="font-size:10.0pt">
<a href="mailto:getdp-bounces@ace20.montefiore.ulg.ac.be" target="_blank">getdp-bounces@ace20.montefiore.ulg.ac.be</a> [mailto:<a href="mailto:getdp-bounces@ace20.montefiore.ulg.ac.be" target="_blank">getdp-bounces@ace20.montefiore.ulg.ac.be</a>]
<b>På vegne af </b>John_V<br>
<b>Sendt:</b> 3. august 2011 14:21<br>
<b>Til:</b> <a href="mailto:getdp@geuz.org" target="_blank">getdp@geuz.org</a><br>
<b>Emne:</b> [Getdp] Von Neumann misconception?</span><o:p></o:p></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt">Consider an electrostatic problem in a cylindrical volume, bottom end cap (surface) constrained to V=0, source terms in some volume elements, and top end cap with a von Neumann constraint
implemented as shown in the formulation below. Does the 3rd Galerkin term (the von Neumann term)<o:p></o:p></p>
<div style="margin-left:30.0pt">
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">a) Require that the electric field through the upper end cap be equal to Eimposed?<br>
b) Require that the electric field through the upper end cap be equal to the sum: field produced by charges in SourceDomain + Eimposed?<br>
c) Something else?<o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt"><br>
Formulation { <br>
{ Name Electrostatics_v; Type FemEquation;<br>
Quantity {<br>
{ Name v; Type Local; NameOfSpace Hgrad_v_Ele; }<br>
}<br>
Equation {<br>
Galerkin { [ epsr[] * Dof{d v} , {d v} ]; In Domain_Ele; <br>
Jacobian Vol; Integration GradGrad; }<br>
Galerkin { [ -q[]*chargeUnit/eps0/ElementVol[] , {v} ]; In SourceDomain; <br>
Jacobian Vol; Integration GradGrad; }<br>
Galerkin { [ Eimposed , {v} ]; In UpperEndCap; <br>
Jacobian Sur; Integration GradGrad; }<br>
}<br>
}<br>
}<br>
<br>
The mesh is a cylindrical volume oriented along z. Domain_Ele refers to all volume elements of my mesh. SourceDomain is a set of volume elements in which there are charges. (These are near the lower end cap of the cylindrical mesh.) The 3rd term was meant to
impose my von Neumann condition. When I wrote it I thought I was doing (a) above, but the result suggests otherwise.
<br>
<br>
John<o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>