fene-p.h File Reference

#include "log-conform.h"

Include dependency graph for fene-p.h:

Go to the source code of this file.

Functions
static void	fenep (double trA, double *nu, double *eta)
void	event_defaults (void)
	Event: defaults (i = 0)
void	event_init (void)
	Event: init (i = 0)

Variables
double	L2 = 1.

Function Documentation

event_defaults()

void event_defaults

(

void

)

Event: defaults (i = 0)

Boundary conditions for VOF-advected tracers usually depend on boundary conditions for the VOF field.

Event: defaults (i = 0)

Event: defaults (i = 0)

Initialisation

We set the default values.

Electrohydrodynamic stresses

The EHD force density, \(\mathbf{f}_e\), can be computed as the divergence of the Maxwell stress tensor \(\mathbf{M}\),

\[ M_{ij} = \varepsilon (E_i E_j - \frac{E^2}{2}\delta_{ij}) \]

where \(E_i\) is the \(i\)-component of the electric field, \(\mathbf{E}=-\nabla \phi\) and \(\delta_{ij}\) is the Kronecker delta.

We need to add the corresponding acceleration to the Navier–Stokes solver.

If the acceleration vector a (defined by the Navier–Stokes solver) is constant, we make it variable.

Event: defaults (i = 0)

Event: defaults (i = 0 )

The default density field is set to unity (times the metric).

We reset the multigrid parameters to their default values.

If the viscosity is non-zero, we need to allocate the face-centered viscosity field.

We also initialize the list of tracers to be advected with the VOF function \(f\) (or its complementary function).

We set limiting.

On trees, we ensure that limiting is also applied to prolongation and refinement.

We add the interface and the density to the default display.

We switch to a pure minmod limiter by default for increased robustness.

With the MUSCL scheme we use the CFL depends on the dimension of the problem.

On trees we need to replace the default bilinear refinement/prolongation with linear so that reconstructed values also use slope limiting.

The restriction/refine attributes of the charge density are those of a tracer otherwise the conservation is not guaranteed.

By default the permittivity is unity and other quantities are zero.

Definition at line 18 of file fene-p.h.

References f_r, f_s, and fenep().

Here is the call graph for this function:

event_init()

void event_init

(

void

)

Event: init (i = 0)

Definition at line 24 of file fene-p.h.

References _i, dim, dimension, L2, trA, and x.

fenep()

static void fenep ( double  trA, double *  nu, double *  eta  )

static

Definition at line 12 of file fene-p.h.

References eta, L2, nu, and trA.

Referenced by event_defaults().

Here is the caller graph for this function:

Variable Documentation

L2

double L2 = 1.

Functions \f$f_s\f$ and \f$f_r\f$ for the FENE-P model

See log-conform.h.

Definition at line 10 of file fene-p.h.

Referenced by event_init(), and fenep().