Commit 39e9a0e4c71dfb8f35928feb38198173e62751e0

Authored by Chris Smartt
1 parent 036ba562

Small update the Theory Manual related to Neumann boundary condition

DOCUMENTATION/SACAMOS_TheoryManual.pdf
No preview for this file type
DOCUMENTATION/THEORY_MANUAL/Tex/Laplace_Solver_Documentation.tex
@@ -123,7 +123,7 @@ The Neumann boundary condition is expressed as @@ -123,7 +123,7 @@ The Neumann boundary condition is expressed as
123 \begin{equation} 123 \begin{equation}
124 \frac{\partial \phi(r) }{\partial r}=0 124 \frac{\partial \phi(r) }{\partial r}=0
125 \end{equation} 125 \end{equation}
126 -This boundary condition enforces a condition that the normal derivative of the potential is zero and hence the charge on the boundary is zero everywhere. The zero charge propertiy is required by the capacitance matrix calculation and using the Neumann boundary condition leads to convergence of the capacitance matrix elements to the correct values as the boundary distance and the mesh density are increased. 126 +This boundary condition enforces a condition that the normal derivative of the potential is zero and hence the charge inside the boundary is zero (by application of the integral form of Gauss' law). The zero charge propertiy is required by the capacitance matrix calculation and using the Neumann boundary condition leads to convergence of the capacitance matrix elements to the correct values as the boundary distance and the mesh density are increased.
127 127
128 \subsubsection{Asymptotic boundary conditions} 128 \subsubsection{Asymptotic boundary conditions}
129 129
@@ -148,7 +148,7 @@ Expression immediately implies that the potential for $r \gg r'$, with $r'$ the @@ -148,7 +148,7 @@ Expression immediately implies that the potential for $r \gg r'$, with $r'$ the
148 \begin{equation}\frac{\partial }{{\partial r}}\phi \left( {r,\theta } \right) = \frac{1}{{r\ln \left( r \right)}}\phi \left( {r,\theta } \right)\end{equation} 148 \begin{equation}\frac{\partial }{{\partial r}}\phi \left( {r,\theta } \right) = \frac{1}{{r\ln \left( r \right)}}\phi \left( {r,\theta } \right)\end{equation}
149 \\ 149 \\
150 150
151 -It has been observed that the use of the asymptotic boundary condition described here can have an effect on the convergence of the capacitance matrix calculation due to a net charge on the outer boundary which is not taken into account in the capacitance matrix calculation described in section \ref{LC_energy}. Currently it is not known whether a 'generalised capacitance matrix' type solution can be used along with this boundary condition as the outer boundary is not an equipotential (as it is in the Dirichelet case). 151 +It has been observed that the use of the asymptotic boundary condition described here can have an effect on the convergence of the capacitance matrix calculation due to the total charge within the boundary not being null, which impacts the capacitance matrix calculation described in section \ref{LC_energy}. Currently it is not known whether a 'generalised capacitance matrix' type solution can be used along with this boundary condition as the outer boundary is not an equipotential (as it is in the Dirichelet case).
152 As a consequence of this the Neumann boundary condition is recommended and used in the SACAMOS software as the default open boundary condition. 152 As a consequence of this the Neumann boundary condition is recommended and used in the SACAMOS software as the default open boundary condition.
153 153
154 154