PDE Namespace Reference

Variables
str	compute_primitive_variables
str	max_eigenvalue
str	flux
str	source_term

Variable Documentation

compute_primitive_variables

str PDE.compute_primitive_variables

Initial value:

=  r"""
  const auto irho = 1.0 / Q[Shortcuts::rho];
  const auto u0 = Q[Shortcuts::rhoU + 0] * irho;
  const auto u1 = Q[Shortcuts::rhoU + 1] * irho;
#if DIMENSIONS == 3
  const auto u2 = Q[Shortcuts::rhoU + 2] * irho;
#endif
 
#if DIMENSIONS == 3
  const auto uSq = u0 * u0 + u1 * u1 + u2 * u2;
  const auto u_n = (normal == 0) ? u0 : (normal == 1) ? u1 : u2;
#else
  const auto uSq = u0 * u0 + u1 * u1;
  const auto u_n = (normal == 0) ? u0 : u1;
#endif
 
  const auto internalE = Q[Shortcuts::rhoE] - 0.5 * Q[Shortcuts::rho] * uSq;
  const auto p = (GAMMA - 1.0) * internalE;
"""

Definition at line 4 of file PDE.py.

flux

str PDE.flux

Initial value:

=  r"""
  {compute_primitive_variables}
 
  F[Shortcuts::rho] = Q[Shortcuts::rhoU + normal];
 
  F[Shortcuts::rhoU + 0] = Q[Shortcuts::rhoU + 0] * u_n;
  F[Shortcuts::rhoU + 1] = Q[Shortcuts::rhoU + 1] * u_n;
#if DIMENSIONS == 3
  F[Shortcuts::rhoU + 2] = Q[Shortcuts::rhoU + 2] * u_n;
#endif
 
  F[Shortcuts::rhoU + normal] += p;
 
  F[Shortcuts::rhoE] = (Q[Shortcuts::rhoE] + p) * u_n;
""".format(
    compute_primitive_variables=compute_primitive_variables
)

Definition at line 34 of file PDE.py.

max_eigenvalue

str PDE.max_eigenvalue

Initial value:

=  r"""
  {compute_primitive_variables}
  const auto speedOfSound = sqrt(GAMMA * p * irho);
  auto result = fmax(0.0, fabs(u_n - speedOfSound));
  result = fmax(result, fabs(u_n + speedOfSound));
  return result;
""".format(
    compute_primitive_variables=compute_primitive_variables
)

Definition at line 24 of file PDE.py.

source_term

str PDE.source_term

Initial value:

=  r"""
  S[Shortcuts::rho]       = 0.0;
  S[Shortcuts::rhoU + 0]  = 0.0;
  S[Shortcuts::rhoU + 1]  = -Q[Shortcuts::rho] * GRAVITY;
#if DIMENSIONS == 3
  S[Shortcuts::rhoU + 2]  = 0.0;
#endif
  S[Shortcuts::rhoE]      = -Q[Shortcuts::rhoU + 1] * GRAVITY;
"""

Definition at line 52 of file PDE.py.