Up-Downwind 1D element for perfect gaz cases
Instead of the standard linear 1D element, a parameter pe is introduced to control the slope. The higher pe the steepest curve => it can even lead to a quasi-discontinuous function. If pe is zero, then the shape functions are the standard ones.
se is linked to the flow direction: se=1 for upwind, se=-1 for downwind.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer(kind=I4), | intent(in) | :: | numn | local node number |
||
real(kind=R8), | intent(in) | :: | ksi | |||
real(kind=R8), | intent(in) | :: | pe | parameter controlling the shape function slope |
||
real(kind=R8), | intent(in) | :: | se | +1 -> upwind, -1 -> downwind |