Volume Terms Differentiated

Differentiated version of evalVolumeIntegrals. Throws an error for unsupported schemes.

Inputs

• mesh

• sbp

• eqn

• opts

• assembler

Computes the jacobian of calcEulerFlux with respect to q. Methods are available for 2D and 3D

Inputs

• params: ParamType, conservative variables only

• q: vector of conservative variables at node

• aux_vars: auxiliary variables at the node

• dir: direction vector (possibly scaled) to compute the flux jacobian in

Inputs/Outputs

• Fjac: flux jacobian, numDofPerNode x numDofPerNode

Aliasing restrictions: params.p_dot is overwritten

Computes the gradient of pressure with respect to q at a node. Methods are available in 2D and 3D

Inputs

• params: ParamType, conservative variables only

• q: vector of conservative variables at the node

Inputs/Outputs

• pdot: vector of length numDofPerNode, overwritten with derivative of p wrt q

Computes the derivative of the strong form volume terms with respect to q, ie.

d/dq (-Q * f)

but only the mesh.numDofPerNode x mesh.numDofPerNode diagonal block

Inputs

• mesh

• sbp

• eqn

• opts

• assembler: used to assemble the contribution of each element into the Jacobian

Computes the derivataive of the volume term of the Roe scheme with respect to q, ie

d/dq (Q^T * f)

Inputs

• mesh

• sbp

• eqn

• opts

• assembler: used to assemble the contribution of each element into the Jacobian