The current config["Boundaries"]["Postprocessing"]["Capacitance"] assumes the boundary for integration is an infinitely thin rectangular prism and thus calculates the 2-sided integral of D \cdot n to get the charge. If the given boundary actually encloses a finite volume, the integral should be 1-sided as it is for a true external boundary.
This might be able to automatically detected by checking if the boundary elements for the given attribute list form a watertight surface.
The current
config["Boundaries"]["Postprocessing"]["Capacitance"]
assumes the boundary for integration is an infinitely thin rectangular prism and thus calculates the 2-sided integral ofD \cdot n
to get the charge. If the given boundary actually encloses a finite volume, the integral should be 1-sided as it is for a true external boundary.This might be able to automatically detected by checking if the boundary elements for the given attribute list form a watertight surface.