Snoopy Version 6.0 (Albert)
|
#include <math.h>
#include <complex.h>
#include "common.h"
#include "gfft.h"
#include "debug.h"
#include "forcing.h"
#include "particles.h"
#include "ltide.h"
Functions | |
void | timestep (struct Field dfldo, struct Field fldi, const double t, const double tremap, const double dt) |
void | sgs_dissipation (struct Field fldi, const double t, const double dt) |
void | implicitstep (struct Field fldi, const double t, const double dt) |
void implicitstep | ( | struct Field | fldi, |
const double | t, | ||
const double | dt | ||
) |
Implicit steps of the integrator (essentially linear diffusion terms)
This is an implicit model: fldi is modified by this routine
fldi,: | (input and output) current status of the flow |
t,: | current time of the simulation |
dt,: | current timestep size |
void sgs_dissipation | ( | struct Field | fldi, |
const double | t, | ||
const double | dt | ||
) |
Subgridscale model This is the Chollet-Lesieur Model (1981) We have nu(k)=nu_i(k)*(E(kc)/kc)^(1/2) nu_i(k)=0.267+9.21*exp(-3.03 kc/k)
NB: this subgridscale model is applied only to the velocity field, even when MHD is active!
This is an implicit model: fldi is modified by this routine
fldi,: | (input and output) current status of the flow |
t,: | current time of the simulation |
dt,: | current timestep size |
void timestep | ( | struct Field | dfldo, |
struct Field | fldi, | ||
const double | t, | ||
const double | tremap, | ||
const double | dt | ||
) |
Compute the right hand side of the INCOMPRESSIBLE dynamical equation
dfldo,: | (output) right hand side of the dynamical equation |
fldi,: | (input) current status of the flow |
t,: | current time of the simulation |
tremap,: | current remap time (only when shear is on) |
dt,: | current timestep size |
Braginskii viscosity. *********** /*
HALL EFFECT