src/fene-p.h

Functions fs and fr for the FENE-P model

See log-conform.h.

double L2 = 1.;

static void fenep (double trA, double * nu, double * eta) {
  *η = 1;
  *ν = 1./(1. - trA/L2);
  return;
} 

event defaults (i = 0) {
  f_s = fenep;
  f_r = fenep;
}

event init (i = 0) {
#if AXI
  double dim = 3;
#else
  double dim = dimension;
#endif  
  scalar trac = trA;
  foreach()
    trac[] = dim*L2/(dim + L2);
}

Usage

Tests