This operator implements an electric field flux of the form
where 
is calculated using a quasi-neutral assumption
Parameters: C D sNA
C is a field, function, or variable for the mobile concentration
D is a field, function, or variable for the diffusivity
sNA is the scaled net active dopant concentration
( 
).