<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7232.11">
<TITLE>2d DC Conduction</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->

<P><FONT SIZE=2 FACE="Arial">I apologise in advance  for sending such a newbie question to the list.  In Payment if someone could</FONT>

<BR><FONT SIZE=2 FACE="Arial">set me on my way I will add my simple example(s) to the main example list.</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">In an effort to walk (more of a crawl) before I run I am trying to simply</FONT>

<BR><FONT SIZE=2 FACE="Arial">calculate the voltage and current density distribution in a flat bar of copper (no temperature effects assumed)</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">Fragments for this exist in the mail archives as well as documentation.  However,</FONT>

<BR><FONT SIZE=2 FACE="Arial">I must be doing something wrong as my results are always zero when using Dirichlet BC at both ends</FONT>

<BR><FONT SIZE=2 FACE="Arial">and constant (but non-zero) when I use dirichlet at one end and Neumann BC at the other.</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">I assume that the null Neumann BC along the long sides is implicit.</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">I have attached my .geo, .pro and solver.par files for the Dirichlet BC (1V at one end 0v at the other)</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">THe Formulation was :</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">Formulation {</FONT>

<BR><FONT SIZE=2 FACE="Arial">  { Name Voltage_Steady ; Type FemEquation ;</FONT>

<BR><FONT SIZE=2 FACE="Arial">    Quantity { </FONT>

<BR><FONT SIZE=2 FACE="Arial">      { Name Volt ; Type Local ; NameOfSpace Vgrad ; }</FONT>

<BR><FONT SIZE=2 FACE="Arial">    }</FONT>

<BR><FONT SIZE=2 FACE="Arial">    Equation {</FONT>

<BR><FONT SIZE=2 FACE="Arial">      Galerkin { [ invrho[] *  Dof {d Volt} ,  {d Volt} ] ;  </FONT>

<BR><FONT SIZE=2 FACE="Arial">                 In Bar ; Jacobian JVol ; Integration GradGrad ; }</FONT>

<BR><FONT SIZE=2 FACE="Arial">    }</FONT>

<BR><FONT SIZE=2 FACE="Arial">  }</FONT>

<BR><FONT SIZE=2 FACE="Arial">}</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">When I run the Dirichlet at both ends I receive a Warning: Null right Hand Side in Linear Equation</FONT>

<BR><FONT SIZE=2 FACE="Arial">which I assume is correct given the problem formulations.  I also get an indication that NaN is</FONT>

<BR><FONT SIZE=2 FACE="Arial">being generated in the SOlver (GMRES).</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">Not surprisingly the Warning goes away when I use the Neumann BC.  THe Solver does not report any</FONT>

<BR><FONT SIZE=2 FACE="Arial">NaN problems.<FONT FACE="Arial" SIZE=2 COLOR="#000000"> <<simpleBar.geo>> </FONT><FONT FACE="Arial" SIZE=2 COLOR="#000000"> <<simpleBar.pre>> </FONT><FONT FACE="Arial" SIZE=2 COLOR="#000000"> <<solver.par>> </FONT></FONT>
</P>
<BR>

<P><FONT SIZE=2 FACE="Arial">Thanks,  the mail archives have been a big help</FONT>
</P>

<P><FONT SIZE=2 FACE="Arial">Laird Bolt</FONT>
</P>

</BODY>
</HTML>