Charge Relaxation Theory
COMSOL Multiphysics includes physics interfaces for the modeling of static electric fields and currents. Deciding what specific physics interface and study type to select for a particular modeling situation requires a basic understanding of the charge dynamics in conductors.
The different physics interfaces involving only the scalar electric potential can be interpreted in terms of the charge relaxation process. The fundamental equations involved are Ohm’s law
the equation of continuity
and Gauss’ law
By combining these, one can deduce the following differential equation for the space charge density in a homogeneous medium
This equation has the solution
where
is called the charge relaxation time. For a good conductor like copper, τ is of the order of 1019 s, whereas for a good insulator like silica glass, it is of the order of 103 s. For a pure insulator, it becomes infinite.
When modeling real-world devices, there is not only the intrinsic time scale of the charge relaxation time but also an external time scale t at which a device is energized or the observation time. It is the relation between the external time scale and the charge relaxation time that determines what physics interface and study type to use. The results are summarized in Table 4-3 below,
τ>>t
τ<<t
τ~t
First Case: τ >> t
If the external time scale is short compared to the charge relaxation time, the charges do not have time to redistribute to any significant degree. Thus the charge distribution can be considered as a given model input. The best approach is to solve the Electrostatics formulation using the electric potential V.
By combining the definition of the potential with Gauss’ law, you can derive the classical Poisson’s equation. Under static conditions, the electric potential V is defined by the equivalence E = −∇V. Using this together with the constitutive relation D = ε0E + P between D and E, you can rewrite Gauss’ law as a variant of Poisson’s equation
This equation is used in the Electrostatics interface. It is worth noting that Gauss’ law does not require the charge distribution to be static. Thus, provided dynamics are slow enough that induced electric fields can be neglected and hence a scalar electric potential is justified, the formulation can be used also in the Time Dependent study type. That typically involves either prescribing the charge dynamics or coupling a separate formulation for this.
Second Case: τ <<t
If the external time scale is long compared to the charge relaxation time, the stationary solution to the equation of continuity has been reached. In a stationary coordinate system, a slightly more general form of Ohm’s law than above states that
where Je is an externally generated current density. The static form of the equation of continuity then reads
To handle current sources, the equation can be generalized to
This equation is used in the static study type for the Electric Currents interface.
General Case: Charge Dynamics
If the charge relaxation time is comparable to the external time scale, the Time Dependent or Frequency Domain study types for the Electric Currents interface must be used.
Combining the time-harmonic equation of continuity
with the equation (∇⋅ D = ρ) yields the following equation for the frequency domain study type:
For the Time Dependent study type, use the transient equation of continuity
and the resulting equation becomes
These dynamic formulations are valid as long as induced electric fields can be ignored and hence the electric field is essentially curl free. This condition is fulfilled provided that skin effect and wave propagation effects can be ignored. The skin depth must be much larger than the geometrical dimensions of the modeled device and so must the wavelength. Note also that these formulations can be used to model dielectric regions of capacitive/resistive devices even though the interior of electrodes sometimes does not meet the large skin depth condition. In that case, the electrodes must only be represented as boundary conditions (fixed or floating potential). The interior, metallic domains are not included in the analysis. Obviously, this is only a valid approach for devices where metallic electrodes do not entirely bypass (short circuit) the capacitive/resistive layers. If metallic electrodes short circuit the capacitive/resistive layers, the time evolution of the current is determined by inductive and resistive effects with very little influence from the capacitive layers. Then the Magnetic Fields interface is the appropriate modeling tool.