electrostaticForce

Syntax

Defined in couplingProperties dictionary.

forceModels
(
    electrostaticForce
);
electrostaticForceProps
{
    relativePermittivity   scalar;
    vacuumPermittivity     scalar;
    useParticleChargeField Switch;

    // (optional) forceSubModel switches
    verbose                Switch;
};

relativePermittivity = relative permittivity of the fluid $\epsilon:$ (dimensionless) field (default: “1”, relative permittivity of vacuum)
vacuumPermittivity = the physical constant $\epsilon_0$ (default: $8.854 \times 10^{-12} \mathrm{F/m}$ )
useParticleChargeField = consider charged particles when calculating the electric potential field (default: true)

This forceModel reads the following forceSubModel switches and overwrites the defaults as indicated in parentheses:

verbose
interpolation

Examples

forceModels
(
    electrostaticForce
);
electrostaticForceProps
{
}

Description

Calculates the electrostatic potential and field in the domain from the particle charge density $\rho_p$ and boundary conditions by solving Poisson’s equation

$\Delta \phi & = - \frac{\rho_q}{\epsilon \epsilon_0}$

$\vec{E} & = - \nabla \phi$

and communicates the field values to Aspherix, where the actual force is computed. This model requires the solver command enable_electric_field command to work, force computation is detailed in the according documentation. Particle charges are only taken into account if useParticleChargeField is on.

$\vec{F} = q \vec{E}$

Note

The user must provide the electric potential field “electricPotential” (units [1, 2, -3, 0, 0, -1, 0]) and can, in addition, provide the electric field “E” (units: [1, 1, -3, 0, 0, -1, 0]).

Restrictions

The model is under development and therefore in experimental state.