Changes
/* Thermal Boundary Conditions */
<math> T(\mathbf{r}) = T_0 </math>
The Neuman Neumann boundary condition constrains the rate of heat flow through the domain boundary walls:
<math>-k \frac{\partial T}{\partial n} = -k \mathbf{\hat{n}} . \nabla T(\mathbf{r}) = q_s </math>
where q<sub>s</sub> is the heat flux passing through the boundary walls.
If q<sub>s</sub> = 0, the Neumann boundary condition is also known as the adiabatic boundary conditions:
<math>\frac{\partial T}{\partial n} = 0 </math>
== Electrostatics Analysis==