Advection Volume Integrals

AdvectionEquationMod.evalVolumeIntegrals

Evaluates the volume integrals of the weak form. eqn.res is updated with the result. Both the precompute_volume_flux and non precompute_volume_flux versions are contained within this function

Inputs:

mesh : mesh type sbp : Summation-by-parts operator eqn : Advection equation object opts : options dictionary

Outputs

None

Populates eqn.flux_parametric. Repeatedly calls the other method of this function.

Inputs:

mesh
sbp
eqn
opts

Calculates the advection flux in the parametric directions at a node.

Inputs:

params: a ParamType object
q: the solution value at the node
alphas_xy: the advection velocities in the x-y directions, vector of length
           Tdim
dxidx: scaled mapping jacobian at the node, Tdim x Tdim matrix

Inputs/Outputs:

flux: vector of length Tdim to populate with the flux in the parametric
      directions