velocity_moments

Pre-calculate spatial derivatives of the electron moments that will be needed for the time advance

Calculate spatial derivatives of the electron moments.

This version, for use in implicit solvers for electrons, works with a single point in r, given by ir.

Pre-calculate spatial derivatives of the moments that will be needed for the time advance

Pre-calculate spatial derivatives of the neutral moments that will be needed for the time advance

create a moment struct containing information about the electron moments

compute the isotropic pressure from the already computed ppar and pperp

computes the integral over vpa <= 0 of the integrand, using the input vpa_wgts this could be made more efficient for the case that dz/dt = vpa is time-independent, but it has been left general for the cases where, e.g., dz/dt = wpa*vth + upar varies in time

computes the integral over vpa >= 0 of the integrand, using the input vpa_wgts this could be made more efficient for the case that dz/dt = vpa is time-independent, but it has been left general for the cases where, e.g., dz/dt = wpa*vth + upar varies in time

computes the integral over vpa of the integrand, using the input vpa_wgts

compute the integral needed for the generalised Chodura condition

IChodura = (Z^2 vBohm^2 / cref^2) * int ( f bz^2 / vz^2 + dfdr*rhostar/vz ) vBohm = sqrt(Z Te/mi) with Z = 1 and mref = mi cref = sqrt(2Ti/mi) and normalise to the local ion density, appropriate to assessing the Chodura condition

IChodura <= (Te/e)d ne / dphi |(sheath entrance) = ni

to a single species plasma with Z = 1

NB: if this function is called and if dens_updated is false, then the incoming pdf is the un-normalized pdf that satisfies int dv pdf = density

calculate the updated density (dens) for a given species; should only be called when evolvedensity = false, in which case the vpa coordinate is vpa/cs

update velocity moments that are calculable from the evolved ion pdf

update velocity moments that are calculable from the evolved neutral pdf

calculate the updated density (dens) and parallel pressure (ppar) for all species this function is only used once after initialisation the function used to update moments at run time is updatederivedmoments! in time_advance.jl

update velocity moments of the evolved neutral pdf

calculate the neutral density from the neutral pdf

calculate the updated density (dens) for a given species

calculate the updated pressure in the rr direction (pr) for a given species

calculate the updated pressure in zz direction (pz) for a given species

calculate the updated pressure in the zeta zeta direction (pzeta) for a given species

calculate the updated heat flux zzz direction (qz) for a given species

calculate the updated ur (mean velocity in r) for a given species

calculate the updated uz (mean velocity in z) for a given species

calculate the updated uzeta (mean velocity in zeta) for a given species

NB: if this function is called and if ppar_updated is false, then the incoming pdf is the un-normalized pdf that satisfies int dv pdf = density

calculate the updated energy density (or parallel pressure, ppar) for a given species; which of these is calculated depends on the definition of the vpa coordinate

calculate the updated perpendicular pressure (pperp) for a given species

NB: the incoming pdf is the normalized pdf

calculate the updated parallel heat flux (qpar) for a given species

NB: if this function is called and if upar_updated is false, then the incoming pdf is the un-normalized pdf that satisfies int dv pdf = density

calculate the updated parallel flow (upar) for a given species