pnp.h File Reference

Go to the source code of this file.

Functions
void	ohmic_flux (scalar *c, int *z, double dt, vector *K=NULL)

Variables
scalar	phi
	The electric potential and the volume charge density are scalars while the permittivity and conductivity are face vectors.

Function Documentation

ohmic_flux()

void ohmic_flux	(	scalar *	c,
		int *	z,
		double	dt,
		vector *	K = NULL
	)

If the volume conductivity is not provided it is set to the value of the valence.

The fluxes of each specie through each face due to ohmic transport are

The specie concentration is updated using the net amount of that specie leaving/entering each cell through the face in the interval \(dt\)

Definition at line 19 of file pnp.h.

References _i, c, dimension, dt, f, i, K, k, phi, s, vectors_append(), x, and z.

Here is the call graph for this function:

Variable Documentation

phi

scalar phi

extern

The electric potential and the volume charge density are scalars while the permittivity and conductivity are face vectors.

Ohmic conduction flux of charged species

This function computes the fluxes due to ohmic conduction appearing in the Nernst–Planck equation. The species charge concentrations are then updated using the explicit scheme

\[ c^{n+1}_i = c^n_i +\Delta t \, \nabla \cdot( K_i c^n_i \nabla \phi^n) \]

where \(c_i\) is the volume density of the \(i\)-specie, \(K_i\) its volume electric conductivity and \(\phi\) the electric potential.

In mgphi we will store the statistics for the multigrid resolution of the electric poisson equation.

Definition at line 34 of file implicit.h.

Referenced by box_matrix(), correct_qz(), dphidt(), event_acceleration(), event_cleanup(), event_defaults(), event_init(), event_metric(), event_pressure(), event_tracer_diffusion(), fraction(), if(), if(), levelset_to_vof(), matrix(), matrix1(), ohmic_flux(), refine_cm_lonlat(), refine_face_y_lonlat(), relax_nh(), relax_nh(), relax_nh1(), residual_nh(), residual_nh(), residual_nh2(), residual_nh3(), and solid().