I am trying to simulate the electric field produced by a charge localized deposition into a dielectric.
The charge deposition is constant in time and is a function of space. The charge conduction must be calculated to allow charges to flow down to a ground element. I used the following simple geometry: a dielectric cube with a grounded boundary.
Typically, the Gauss’s law, the Ohm’s law and the continuity equation are combined to give the following equation to solve (with a source term):
- ∇.E = ρ/ε
- J = σE
- ∂ρ/∂t + ∇.J = ρ_s
I used the Static Current Conduction model (StatCurrentSolveVec) to calculate the potential as a function of time. I used a user function to inject charges in only one element of the volume. This user function is called into the body force section of the .sif. Then I used the FluxSolver to calculate the electric field from the potential.
The potential and the electric field are well calculated for each step time but they are constant. I think the problem comes from the body force: the charges do not build up and do not flow.
Does the StatCurrentSolveVec allows the transport of charges?
How to set the body force to inject charge into an element?