Discharges in Liquids
Governing Equations
The Electric Discharge interface provides built-in charge transport models in liquids. The default model consists of a group of transport equations for each modeled charge carrier (Ref. 3):
where
e, p, and n denote electrons, positive ions, and negative ions
ni is the number density of the charge carrier (SI unit: 1/m3)
E is the electric field (SI unit: V/m)
zi denotes the carrier charge (SI unit: 1)
μi denotes the carrier mobility (SI unit: m2/(V·s))
wi is the drift velocity in the electric field (SI unit: m/s)
Di is the diffusion coefficient (SI unit: m2/s)
Ri is the reaction rate (SI unit: 1/(m3·s))
τa is the attachment time constant (SI unit: s)
βep is the electron–ion recombination coefficient (SI unit: m3/s)
βpn is the ion–ion recombination coefficient (SI unit: m3/s)
SF is the field ionization (SI unit: 1/(m3·s))
e is the electric charge (SI unit: C)
nioni is the number density of ionizable species (SI unit: 1/m3)
a denotes the molecular separation distance (SI unit: m)
m* denotes the effective electron mass (SI unit: kg)
φΔ and φγ are ionization potential parameters (SI unit: V)
The above transport equations are fully coupled with Poisson’s equation through the electric field and the space charge:
Dissociation
Under a high electric field, dielectric liquids undergo dissociation, commonly referred to as the Wien effect, as described by the following equations (Ref. 9):
where
Sd denotes dissociation rate (generating free ions, SI unit: 1/(m3·s))
βpn is the ion–ion recombination coefficient (SI unit: m3/s)
n0 is the zero-field number density (SI unit: 1/m3)
σ0 is the zero-field electric conductivity (SI unit: S/m)
μp is the positive ion mobility (SI unit: m2/(V·s))
μn is the positive ion mobility (SI unit: m2/(V·s))
F is the Onsager function
I1 is the modified Bessel function of the first kind and order 1
lB is the Bjerrum length (SI unit: m)
lO is the Onsager length (SI unit: m)