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?

Thanks

Jeremie