On 07/20/2017 04:23 PM, Baoyun Ge wrote:

I am using dealii to solve the electrostatic field governed by Laplace's equation. Right now I have the solution vector, which contains the electrical potential information throughout all the degrees of freedom. I am stuck at the post processing phase, where I want to evaluate the charges on the electrodes. The way I am seeking it is in two steps:

1) evaluating electric field using the electrical potential information;
2) integrating electric field over the surface enclosing the electrode under consideration (Gauss's Law).

Unfortunately, I haven't found too much useful information to guide me to accomplish these two steps. Any information and guidance from you guys would be appreciated!

Baoyun, the electrical field is just the gradient of the potential, so this is easy to evaluate via FEValues::get_function_gradients(). For the integration, what surface do you want to integrate over? Any surface? A surface that is composed of faces of cells? The former is difficult, the latter easy to achieve with the FEFaceValues class.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
For more options, visit https://groups.google.com/d/optout.

Reply via email to