Inviscid solid-wall (no-penetration) face flux through unit normal nrm.
At a solid wall v.n = 0, so the Euler flux collapses to the pressure term:
F = [ 0, pnx, pny, 0 ],
with p the interior-side pressure and (nx,ny) the boundary-face normal.
Mass and energy flux are EXACTLY zero by construction (a "weak"/pressure
wall), so no spurious mass crosses the wall and no spurious entropy is
generated by an acoustic over-pressure -- unlike a ghost-reflection wall,
whose Riemann star pressure carries an O(rhoc|v_n|) error wherever the
reconstructed interior normal velocity is non-zero (always, on a curved
wall). q is the interior-side (reconstructed) conserved state; gam the
ratio of specific heats. nrm need not be normalised in sign: the residual's
+xi/+eta divergence orientation makes the resulting wall pressure force act
into the domain at both min and max edges.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | q(neq2d) | |||
| real(kind=wp), | intent(in) | :: | nrm(2) | |||
| real(kind=wp), | intent(in) | :: | gam | |||
| real(kind=wp), | intent(out) | :: | flux(neq2d) |